3. 矿山点云预处理与误差敏感环节排查

在矿山环境下，激光雷达原始点云并不是“拿来就能配准”的干净数据。粉尘、水雾、飞石、车辆扬尘、皮带机架、管线、支护网以及大面积裸露岩壁，都会让点云中混入大量非结构点或伪结构点。预处理环节承担着把原始点云整理成“可配准、可约束、可判断”的几何观测的任务。问题在于，预处理既是抑制误差的第一道关口，也是最容易把误差“放大后传给后端”的敏感环节：滤波参数不当会削弱真实结构，地面分割失败会让位姿约束建立在错误平面上，特征提取退化则会让 ICP 或特征匹配在错误方向上收敛。本节围绕这三类敏感环节展开，并给出一次可复现的预处理参数对比实验。

一、预处理链路与误差放大路径

典型的矿山点云预处理链路为：原始点云 → 体素降采样/统计滤波 → 地面分割 → 特征提取（面、线、角点或曲率显著点）→ 配准与位姿估计。误差放大并非发生在单点，而是沿链路逐级传递：降采样体素过大，会把巷道壁面细节抹平，导致面特征法向估计不稳；统计滤波的邻域与标准差阈值设置过激，会把远处稀疏但真实的巷道轮廓当作离群点剔除，使配准在远端失去约束；地面分割阈值不当，会把斜坡道、台阶、堆料坡误判为地面或把真实地面漏掉，导致 z 向与俯仰角出现系统性偏差；特征提取在长直巷道、空旷采场等退化场景中提取到的点几乎共线或共面，配准的雅可比矩阵在某个方向上接近奇异，位姿在该方向上的不确定性被迅速放大，表现为沿巷道方向的漂移。

二、滤波参数不当：从“去噪”到“去结构”

统计滤波的核心参数是邻域点数 k 与标准差倍数 α，判定条件为：若某点到其 k 个近邻的平均距离超过全局均值加 α 倍标准差，则判为离群点。矿山巷道远端点云稀疏，平均近邻距离天然偏大，若 α 取 0.5 这类偏小值，会把远端真实轮廓大量剔除，配准在远端失去约束，长直巷道方向漂移加剧。体素降采样则要兼顾细节与效率：体素边长 0.1 m 可保留支护网与管线细节，但点数大、耗时长；0.5 m 会抹平壁面起伏，使面特征法向抖动。经验上，巷道场景体素取 0.15～0.25 m，k 取 20～40，α 取 1.0～1.5，并优先保留远端稀疏结构。

三、地面分割失败：错误平面带来的系统性偏差

常用做法是拟合地面平面或采用坡度/高度栅格判定。矿山地面并非理想水平面：斜坡道坡度可达 8°～15°，台阶与堆料坡局部更陡。若坡度阈值设为 5°，斜坡会被判为非地面而漏分，地面约束缺失；若阈值放宽到 25°，堆料坡与车辆顶面又会被误判为地面，把错误平面引入位姿估计，造成 z 向与俯仰角的系统性偏差。合理做法是分区域设定阈值，并对分割结果做平面度与内点比例检查：内点比例低于 60% 或拟合残差标准差大于 0.1 m 时，应判定分割可疑并回退到局部地面块拟合。

四、特征退化：配准在“看不见的方向”上漂移

长直巷道中，壁面点近似共面、顶底板点近似平行，沿巷道轴向的几何约束极弱，配准信息矩阵在该方向上的特征值很小，微小噪声即可引起较大位姿偏移。排查时应统计特征点在三个主方向上的分布：若某方向方差占比超过 85%，说明该方向约束退化，需要引入里程计、UWB 或已知控制点作为补充约束，或改用面特征加权的配准策略，而不是单纯调大迭代次数。

五、预处理参数对比实验（示例）

以下为一次教学示例实验设计，数据为合理假设的巷道段点云，用于说明对比方法而非实测结论。固定配准算法与初值，仅改变预处理参数，记录配准后的相对位姿误差与特征退化指标。

实验操作上，可用如下命令批量跑四组参数并输出误差统计（示例命令，参数按实际工具调整）：

for cfg in A B C D; do
  run_preprocess --voxel $V --k $K --alpha $A --slope $S \
    --input tunnel_seg.pcd --output out_$cfg.pcd
  run_registration --src out_$cfg.pcd --ref ref.pcd \
    --report err_$cfg.txt
done
grep -E "axial_err|z_bias|inlier_ratio" err_*.txt

判读要点：轴向误差主要受降采样与滤波影响，z 向偏差主要受地面分割影响，内点比例下降往往先于误差增大出现，可作为早期预警指标。若某组内点比例低于 60%，即使当前误差尚可，也不应作为长期参数使用。