3. 点云数据预处理与断面分析准备

原始激光点云从矿山现场采集回来后,并不能直接用于断面偏差和超欠挖计算。原因很直接:一个测站的点云只覆盖局部区域,多个测站之间存在坐标基准差异;扫描过程中会混入粉尘、水雾、飞鸟、行人、车辆、支护设备等非目标点;密集点云数据量巨大,直接参与断面拟合会显著拖慢计算并放大噪声影响。因此,断面分析前的预处理本质上是一次“把原始观测数据整理成可比较、可计算、可复现的几何基准”的过程。本节围绕去噪、抽稀、坐标统一、多测站拼接、质量检查与断面提取前的数据整理展开,重点说明每一步的实际操作对象、参数选择依据和结果判断方法。

一、去噪:先分清“噪声类型”,再选滤波策略

矿山点云噪声大致分为三类。第一类是孤立飞点,通常由粉尘、水雾反射或设备边缘衍射产生,表现为远离主体表面的稀疏点;第二类是非目标物体点,如行人、车辆、管线、临时支护材料,它们本身是真实回波,但不属于巷道或采场表面;第三类是表面毛刺,表现为附着在岩壁表面附近的随机抖动点,幅度一般在几厘米以内。三类噪声的处理方式不同,不能一律用同一种滤波参数。

对孤立飞点,常用统计离群点滤波(SOR)。其思路是统计每个点到最近 k 个邻域点的平均距离,若该距离超过全局均值加若干倍标准差,则判为离群点。教学示例中可取邻域点数 k=20、标准差倍数 1.5~2.0;倍数越小去噪越激进,容易误删真实细节(如锚杆头、电缆挂钩),倍数越大保留细节但残留飞点多。对表面毛刺,可用半径滤波或体素栅格内的中值/均值平滑,但要注意平滑会削弱断面上的真实棱角和超挖尖点,因此断面分析前一般只做轻度平滑,不做强平滑。

非目标物体点最可靠的处理方式是“先分割、再人工剔除”。可先按高程、距巷道中线的横向距离或连通域尺寸做粗分割,把明显不属于巷道轮廓的大块点云(如停放的铲运机)分离出来,再在可视化软件中框选删除。需要强调的是,去噪不是越干净越好:支护锚杆、挂网、管线本身会影响断面净空判断,是否保留应依据计算目的决定——计算“超挖”时通常关注岩面,喷浆面和支护构件应单独标记,避免与岩面混算。

二、抽稀:在保留几何特征的前提下降低冗余

地面激光扫描在近距离处点间距可能只有几毫米,而断面拟合并不需要如此高的密度。抽稀的常用方法是体素栅格下采样:设定一个体素边长(如 2 cm 或 5 cm),每个体素内保留一个代表点(重心点或距体素中心最近的点)。体素边长应与后续断面提取的点间距需求匹配——若断面按 10 cm 间隔采样,体素取 2~5 cm 通常足够;若需要识别小尺度超挖尖点,体素不宜大于 5 cm。

抽稀的另一个作用是均匀化点密度。多测站拼接后,近站点云密、远站点云稀,直接拟合会使近站区域权重过大。体素下采样能在一定程度上均衡密度分布,使断面拟合更稳定。抽稀后应记录体素边长和保留点数,作为后续复现的参数依据。

三、坐标统一与多测站拼接

多测站拼接的核心是把各测站独立坐标系下的点云统一到同一个工程坐标系。常用做法分两步:先以标靶球或标靶纸做粗配准,再用 ICP(迭代最近点)做精配准。粗配准提供初始变换,精配准在重叠区域上迭代最小化对应点距离。教学示例中,若两站重叠区域点云充足,ICP 的对应点距离阈值可设为 5~10 cm,迭代 30~50 次,最终配准残差(RMS)宜控制在 1~2 cm 以内;若残差明显偏大,应检查标靶识别是否正确、重叠区域是否被移动物体污染。

坐标统一还包括把拼接结果转换到矿山工程坐标系(如以巷道中线为 X 轴、竖直向上为 Z 轴)。这一步必须明确参照约定:断面偏差的正负号、超欠挖的判定方向都依赖坐标轴定义。建议在数据整理阶段就固定“以设计断面为基准,实测点位于设计轮廓外侧为正偏差(超挖),内侧为负偏差(欠挖)”,并在后续所有图和表中保持一致,避免中途切换参照系导致符号混乱。

原始多测站点云 → 可计算断面基准的数据整理链 原始点云 多测站 / 含噪声 密度不均 ① 去噪 SOR 去飞点 分割剔除车辆/行人 ② 抽稀 体素下采样 边长 2~5 cm ③ 多测站拼接 标靶粗配准 + ICP 残差 RMS ≤ 2 cm ④ 坐标统一 转工程坐标系 固定正负号约定 ⑤ 质量检查 重叠区残差 空洞 / 密度 / 覆盖 ⑥ 断面提取准备 设定断面位置/间距 裁切计算范围 可计算数据集 统一坐标 + 干净点云 参数记录完整 残差超限 → 回查标靶与重叠区 预处理参数记录(示例,需按实际工程填写) SOR:k=20,σ 倍数=1.5~2.0 | 体素边长:2~5 cm | ICP 距离阈值:5~10 cm,迭代 30~50 次 配准残差 RMS:≤ 2 cm | 断面间距:按设计桩号 5~10 m | 正负号约定:外侧为正(超挖),内侧为负(欠挖)

四、质量检查:用可量化的指标判断“能不能用”

预处理完成后必须做质量检查,而不是凭肉眼“看着差不多”。建议至少检查四项:其一,重叠区配准残差,取相邻测站重叠带内对应点距离的 RMS,若超过 2 cm 应回查标靶和重叠区;其二,点云空洞,统计断面范围内沿巷道轴向的连续无点区长度,若某段超过 0.5 m 且无合理解释(如设备遮挡),该段断面不宜参与偏差统计;其三,密度均匀性,可在近站和远站各取 1 m² 区域统计点数,若差异超过 3 倍,说明抽稀或拼接存在问题;其四,覆盖完整性,检查设计断面轮廓是否被实测点云完整包络,特别是拱顶和墙脚,缺失部位会导致超欠挖计算出现虚假极值。

五、断面提取前的数据整理

最后一步是把点云整理成“按断面可切、按桩号可查”的状态。具体做法是:沿巷道中线按设计桩号设定断面位置(教学示例取 5~10 m 间距),在每个断面位置建立垂直于中线的切平面,裁切一定厚度(如 10~20 cm)的点云带作为该断面的计算数据;同时把设计断面轮廓线、中线、高程基准一并导入,确保实测点与设计轮廓在同一坐标系下比较。整理完成后应形成一份可追溯的记录:测站编号、拼接残差、去噪参数、体素边长、断面桩号列表、正负号约定。这份记录既是质量凭证,也是后续复核和复现的依据。

判断预处理是否合格的实用标准:把整理后的点云按断面切开,若断面轮廓连续、无明显飞点、拱顶与墙脚均有实测点、且相邻断面之间轮廓变化平滑,则基本满足超欠挖计算要求;若出现断面缺口或轮廓跳变,应先回到质量检查环节定位原因,不要直接进入偏差计算。

需要提醒的是,预处理各步骤之间存在耦合:去噪过强会削弱真实超挖尖点,抽稀过粗会丢失小尺度偏差,拼接残差偏大会使整段断面的偏差整体偏移。因此参数选择应围绕“断面偏差计算需要多高的几何保真度”来定,而不是追求点云越干净、越稀疏越好。把每一步的参数和判断依据记录下来,才能让后续的超欠挖计算结果可解释、可复核。