地面点滤波实现原理一、功能该功能机载/无人机点云的地面点提取模块。实现目标不是“训练模型”而是在已知三维点集上依据几何与形态学规则把点划分为地面inlier与非地面outlier为稀疏 DEM、DSM 等栅格产品提供可靠输入。工程上采用两阶段策略先用规则格网 形态学开运算得到较稳的初始地面种子粗分类再用渐进三角网密化PTD在二维 Delaunay 三角网上迭代“生长”地面点精分类。二者通过统一的平面格网索引格网单元尺寸、角点对齐、邻域查询耦合超内存场景由分块 重叠带 角点对齐保证块间索引与标签一致。二、实现架构类职责与内存模型2.1 分层设计层次职责关键状态点云容器抽象统一三维点存储、地面布尔标记、PLY/矢量读写点数组、地面标记数组、内外点索引列表格网索引将连续 XY 离散为单元每单元存点索引列表最小角点、单元边长、行列数、平均密度、格网到点索引映射形态学滤波在格网 DEM 上做开运算与多尺度/自适应扩展格网最低 Z、空洞掩膜、形态学 DEM 层PTD 主滤波种子 虚拟边界点 OpenCV Subdiv2D 迭代密化种子索引表、三角网细分器、坡度正弦/余弦预计算大规模编排按内存预算切 XY 矩形、扩展带、块内掩膜、合并标签核心矩形列表、扩展矩形列表、块可用性掩膜工具链粗差剔除、空洞 TPS 插值、精度评估、可选图割变体与主路径解耦按需挂接主路径实现类为PTD 滤波器形态学为独立类但在种子阶段被主类内嵌调用非继承多态而是组合。2.2 地面标签的生命周期形态学阶段写入布尔地面标记数组可多尺度取并。PTD 阶段在标记数组上继续置真迭代结束后同步维护内点/外点索引列表供导出。分块阶段块内产生完整标记数组再与块核心区域可用性掩膜按位与写回全局结果仅核心区内点参与最终地面判定扩展带内点只服务于边界三角网稳定。实现上地面判定以布尔标记数组为准索引列表在滤波结束时由标记数组扫描生成。三、类型一平面格网索引3.1 为何需要格网PTD 每轮迭代要在“新接受地面点”影响区域内找未分类候选点形态学要在每个单元取代表高程通常最低 Z。若对全点云做 KD 树邻域搜索大规模点云 IO 与 cache 都不友好。实现选择固定边长的规则格网空间换时间且与形态学、分块对齐天然一致。3.2 离散化规则给定单元边长dBlkSize原点为(dMinX, dMinY)由点云四至向下取整或角点对齐分块时传入统一角点保证各块格网 I/J 一致。点(x,y)的格网索引idxI floor((x - dMinX) / dBlkSize)实现中 I/J 与 X/Y 对应关系以代码为准读写时勿混淆行列。容器为mapIdxIJ, vector点索引稀疏格网只存有数据的单元。角点对齐设置角点坐标与倍数关系后强制dMinX/dMinY对齐到统一网格避免分块边界处同一物理位置落入不同 I/J造成候选点漏检或重复处理。3.3 初始化与密度估计自动参数自动格网尺寸逻辑实现者在调参时需理解用试探性较细格网插入全点统计有效格网数有至少一个点的单元。平均密度 ≈ 总点数 / (有效格网数 × 单元面积)。若密度低于阈值认为点过于稀疏自动参数失败。目标“每格约默认点数如 6 个”反解边长dBlkSize ≈ sqrt(目标点数/密度)并钳制在[0.2m, 0.1×最大地物尺度]。若空洞格网占比 1%用形态学外圈填充估计空洞率按比率放大边长减少“碎格网”。山地判定在较粗格网边长取最大地物尺度上比较邻格高差超过约tan(5°)×地物尺度则标记为坡地供后续图割变体或坡度标记使用。3.4 栅格化接口形态学与 PTD 共用最低 Z 栅格每格有点则取 Z 最小者无点则标记 void形态学种子、PTD 边界扩展用。最高 Z 栅格用于特定分析或剔除异常低点。邻域查询给定中心格 (i,j) 与窗口半径 NN收集邻域内所有点索引或要求邻域内至少MinNum个点才返回虚拟种子赋 Z 时用。3.5 实现注意点边长过小 → 格网数量爆炸代码中有约 1e8 格网上限检查。插入时检查边长合法性若单格点数过多自动参数回退。OpenMP格网填充、邻域收集等循环可并行格网映射插入需注意线程安全实现中部分路径关闭 OMP 或使用 critical。四、类型二形态学初始地面种子4.1 算法思想形态学路径把点云投影为二维 DEM 格网每格一个高程在栅格上做开运算先膨胀后腐蚀结构元素尺寸与最大地物尺度相关理论上移除尺度小于结构元素的凸起地物保留近似“地面包络”。实现不是单次开运算结束而是填洞 → 去异常低点 → 迭代开运算 → 再插值 → 可选格网域“密化”与 PTD 同名但发生在栅格上。4.2 栅格构建流程与 PTD 相同格网插入点云。每格取最低 Z生成浮点高程矩阵无点格标记 void。Inpainting小窗口插值补全 void避免开运算在空洞处发散。TPS 多尺度插值可选路径按格网尺寸与比例参数对 void 区做薄板样条插值适合大面积无点格。RemoveLowerOutliers剔除相对邻域过低的高程异常实现用固定阈值量级防止坑洞点拉低开运算面。4.3 开运算迭代OpeningIteration结构元素边长像素≈最大地物尺度 / 格网边长取奇数。对高程格网膨胀邻域内取最大填平小坑。腐蚀邻域内取最小削掉小凸起。输出形态学 DEM uchar 标记有效地面格。新 void 掩膜被判定为非地面的格。迭代可重复自适应路径使用MultiScaleGrid在多个放大倍率的格网分辨率上构建层叠格网在粗层上判断地形特征平地/坡地再指导细层开运算参数——实现上把“多尺度”放在格网金字塔而非点云金字塔。4.4 从栅格回到点云标签开运算结束后对每个有效格网若该格为地面格则格内所有点标记为地面候选布尔数组置真。多尺度外层循环PTD 种子阶段对不同尺度因子如 1.0、2.5、5.0 倍边长重复调用形态学类结果逻辑或合并提高复杂地形召回率。扫描线分析大最大地物尺度时另一种子路径用多方向扫描线分析格网高程适合特定生产配置。4.5 与 PTD 的衔接PTD 主滤波器在滤波开头小最大地物尺度 → 自适应形态学同时得到坡地标记数组。否则 → 多尺度形态学或扫描线。形态学不维护三角网只输出初始布尔标记数组。PTD 随后重建格网边长可能与形态学阶段不同以 PTD 参数为准。五、类型三渐进三角网密化 PTD核心5.1 理论对应与实现选择PTDProgressive TIN Densification经典流程由稀疏地面种子构建 TIN在三角面片上对未分类点做坡度 高差检验通过则加入地面并更新 TIN。本实现三角网为OpenCV Subdiv2D平面 Delaunay顶点存于种子索引表前若干为虚拟边界点索引 0 表示虚拟点本身其后为真实点云索引或特殊编码。插入点在 XY 平面定位Z 参与生长判据。5.2 边界虚拟种子ExtendingSeeds问题Delaunay 在凸包外三角面片无数据边缘大量候选点无法被合理三角形覆盖易误判或停滞。实现边界虚拟种子扩展在点云平面四至矩形上按格网边长生成边界格网序列边距与最大地物、格网边长相关。在边界格网位置布置平面虚拟点对每点在邻域格内找最近真实种子点复制其Z避免虚拟点 Z 任意。虚拟点写入虚拟种子数组同时扩展 Delaunay 的 rect 边界整数像素矩形略大于点云范围。真实种子再插入三角网先虚拟点再所有布尔标记为真的点。虚拟点不参与最终地面导出索引表映射约定仅用于撑开TIN。5.3 三角网初始化与插入语义初始化 Delaunay 矩形与预留顶点数约为点数一半经验。首次插入遍历虚拟种子 → 遍历当前地面标记为真的点。后续迭代插入仅插入本轮新接受的地面点记录每条插入边在细分中的编号用于计算受影响矩形在格网 I/J 四至范围内传播。保持边引用一致失败插入则回滚种子索引表。5.4 单次迭代的三步DensifyingIteration实现为固定循环最大迭代次数可配置为 0 时仅输出种子步骤 1 — 插入新地面点到 TIN将上一轮接受的地面点依次插入 Subdiv2D对每点记录插入是否成功与退化、重复、边界有关。影响矩形沿三角网邻接边 walk更新 I/J 最小最大范围实现中存为连续四个整数 minI,minJ,maxI,maxJ。若本轮无新点且非首轮影响矩形为空 →迭代终止收敛。步骤 2 — 收集候选点FindCandidates对每个影响矩形格网四至范围枚举矩形内格网 → 汇总格内尚未分类的点索引布尔标记为假。去重后得到候选列表。这是 PTD局部性的关键不用全图扫描只扫受影响格网。步骤 3 — 密化判定Densifying对每个候选点在 XY 平面定位所在三角形Submotion locate取三角形三顶点坐标通过种子索引表映射到三维点虚拟点用虚拟数组。计算三角形平面方程参数axbyczd0归一化法向。对每个相关三角形实现中常处理与候选点关联的边相邻的三角面执行PointGrowing。PointGrowing 判据标准模式平面法向坡度预筛若|c| sin(坡度阈值)法向近水平该三角形不参与生长过陡坡面。点在三角形内XY 重心坐标或等价 inside 测试。竖直投影在三角形内VerticalPointInTriangle将点沿 Z 投影到平面投影点须在三角形 XY 投影内避免跨越垂直遮挡。IsGroundPoint点到平面距离|axbyczd|。点到三顶点最短欧氏距离的平方minDist2。判据距离/√minDist2 ≤ sin(坡度阈值)的等价形式实现用预计算 sin/cos 坡度阈值。另高差dfThr_DZ在平面距离中隐含使用。通过则标记该候选点为地面。记录其关联三角网边编号供下轮插入。每个三角形还带最低 Z 候选与当前最低点两个索引实现用于在三角面内优先生长最低点。边缘保持模式bPreserveEdge若标准判据失败查当前边邻接三个三角形的法向与连通性IsConnectedToFacet要求与至少一个邻接面片在法向夹角上“足够水平”cos 阈值。高程差不超过垂直容差。用于陡坎处避免把坎下点强行拉入缓坡三角面。并行按三角形分块 OpenMP每线程局部 map 合并边—新点映射。若本轮新地面点数为 0 → 迭代结束。5.5 可选后处理栅格插值bFillOuliers若开启填充且提供范围矢量在格网标记“空洞”单元无地面点、在 shp 范围内。用 TPS 等方式在空洞格网插值生成合成地面点存入插值点共享容器供 DSM 产品补洞。与形态学里的 TPS 插值类似但发生在分类完成之后。5.6 可选变体图割迭代FilterByGraphCut实现中存在实验/备用路径先构建全点 Delaunay在三角网边上算平滑项权重法向夹角相关顶点上算数据项权重与地面可能性相关用 OpenCV GCGraph 做 s-t 最小割迭代更新地面标记。生产主路径仍以坡度高差生长为准图割体现“全局一致性”思路与 PTD 局部生长形成对照。实现者扩展算法时可参考权重计算与坡地标记数组的联动。六、类型四大规模分块6.1 动机单块加载全点云时内存约正比于点数PTD 还需 Submotion、候选缓冲、临时 map。实现通过XY 平面分块控制峰值内存。6.2 块尺寸与扩展带由用户内存上限MB与平均密度估算核心块边长允许点数 ≈ 内存系数 × 1024² / 每点字节估算。块面积 ≈ 点数/密度边长为面积开方。扩展带宽度≈最大地物尺度 格网边长实现中两者相加。核心块为互不重叠除最后一行/列边界处理的矩形扩展块 核心块每边外扩扩展带宽度。扩展带内点参与该块 PTD三角网、候选搜索但写回全局结果时仅核心矩形内点保留地面标签可用性掩膜与结果按位与。6.3 块边界与索引一致性全图统一角点坐标传入各块 PTD 实例保证格网 I/J 对齐。最后一行/列块在点归属上用左闭右开或闭区间区分避免角点重复计数或遗漏实现中对 lastX/lastY 有四分支。6.4 块级执行流程对每块扩展矩形从 PLY/SHP按包围盒裁剪读点只加载扩展矩形内点。构造块内 PTD 对象若参数未完整则块内自动参数估计。块内关闭内层 OMP若外层块间并行已开启避免过度订阅 CPU。执行完整 PTD 滤波流程种子→格网→虚拟边界→迭代。将块内地面标记与核心掩膜合并追加写地面点/矢量文件分块输出模式。块间并行OpenMP parallel for 块编号写文件用 critical 保护。6.5 密集点云标志若估算块边长小于约4×扩展带置密集点云标志后续可调整策略警告日志——实现者对超密航高场景应预期分块更碎、内存更高。七、类型五辅助模块模块实现要点粗差剔除基于离散粗糙度等可在 PTD 前挂接主路径中默认注释掉空洞填补独立类TPS/格网插值与 PTD 后处理类似精度评估分块加载、对比参考地面服务算法调参形态学工具结构元素生成、栅格形态学、扫描线等供 Morpho 调用几何运算三角形内判定、平面参数、连通 facet 判据供 PTD 调用八、关键参数在实现中的含义参数实现绑定格网间距格网索引边长、形态学格网分辨率、候选搜索粒度自动估计时与密度挂钩最大地物尺度形态学结构元素尺寸、虚拟边界扩展宽度、分块扩展带主分量、自动格网上限钳制坡度阈值预计算 sin/cos三角形法向筛选与距离/水平距离比高差阈值平面距离绝对值边缘保持邻接面片高程差最大迭代次数密化循环上限为 0 则仅形态学种子保持边缘切换 PointGrowing 与 PointGrowing_PresEdge山地标志自动参数中坡地检测形态学自适应关联坡地标记分块开关入口选择 Massive 或 Simple 路径块内存上限反推核心块几何尺寸块间并行块间 OpenMP块内关闭 OMP填充开关是否 TPS 插值合成点需范围矢量九、并行与内存策略PTD 迭代三角形级 OpenMP合并 map 需注意。候选收集影响矩形数量大时矩形串行循环格网数量可控。分块块间并行 写盘 critical单块内倾向串行 OMP 避免嵌套。DEM 行优化若将来在类似 GDAL 变换中复用同一行点共享 DEM 缓冲可降 IO——PTD 本体是内存点云。Subdiv2D 插入大量退化边时插入失败会回滚种子索引日志需关注。十、与稀疏 DEM 产品线的衔接上游典型消费地面点点云PLY 或 GPKG 等。可选插值合成点填充空洞。布尔标记或分类字段写入矢量。它不生成 DEM 本身但地面点的 XY 分布与 Z 质量直接决定后续栅格化插值是否稳定。分块接缝处若角点未对齐会出现条带状地面遗漏扩展带宽度小于地物尺度时坎边三角网不稳定。十一、五类任务与主流程对照实现者速查点云输入 → [可选] 粗差剔除 → 形态学格网最低Z → 开运算 → 点标记种子 → 格网索引 Insert 全点 → 边界虚拟种子 Delaunay 初始化 → 循环插入新地面 → 影响矩形 → 候选点 → 三角面生长判据 → [可选] 空洞 TPS 插值点 → 输出地面标记 / 导出类型核心数据结构核心数值操作格网索引格网→点索引映射离散化、邻域枚举、最低 Z 栅格形态学种子高程格网 uchar/float开运算、插值、格网→点 ORPTD 密化Submotion 种子索引表平面定位、坡度/高差、迭代插入大规模分块核心/扩展矩形列表裁剪读点、掩膜写回、角点对齐辅助工具—几何判定、评估、填洞十二、实现扩展与排错建议种子过少查形态学 void 比例、开运算结构元素是否过大、格网边长是否过大。迭代不收敛/极慢影响矩形是否过大单轮接受点过多最大迭代次数三角网退化。块接缝角点、扩展带、lastX/lastY 归属逻辑写回是否只用核心掩膜。陡坡漏分坡度阈值、边缘保持是否关闭、山地自动参数是否生效。城区建筑残留减小最大地物尺度或格网边长加强形态学多尺度。十三、总结它的实现本质在规则格网上用形态学开运算得到初始地面种子在 OpenCV 平面 Delaunay 三角网上用坡度与高差约束迭代生长地面点并用格网索引把每轮搜索限制在局部区域超大数据通过分块与扩展带在内存可控下重复同一套逻辑。扩展算法应优先改种子质量、生长判据、虚拟边界与分块几何四者而非仅替换对外接口。
地面点滤波实现原理
发布时间：2026/5/20 7:48:07
地面点滤波实现原理一、功能该功能机载/无人机点云的地面点提取模块。实现目标不是“训练模型”而是在已知三维点集上依据几何与形态学规则把点划分为地面inlier与非地面outlier为稀疏 DEM、DSM 等栅格产品提供可靠输入。工程上采用两阶段策略先用规则格网 形态学开运算得到较稳的初始地面种子粗分类再用渐进三角网密化PTD在二维 Delaunay 三角网上迭代“生长”地面点精分类。二者通过统一的平面格网索引格网单元尺寸、角点对齐、邻域查询耦合超内存场景由分块 重叠带 角点对齐保证块间索引与标签一致。二、实现架构类职责与内存模型2.1 分层设计层次职责关键状态点云容器抽象统一三维点存储、地面布尔标记、PLY/矢量读写点数组、地面标记数组、内外点索引列表格网索引将连续 XY 离散为单元每单元存点索引列表最小角点、单元边长、行列数、平均密度、格网到点索引映射形态学滤波在格网 DEM 上做开运算与多尺度/自适应扩展格网最低 Z、空洞掩膜、形态学 DEM 层PTD 主滤波种子 虚拟边界点 OpenCV Subdiv2D 迭代密化种子索引表、三角网细分器、坡度正弦/余弦预计算大规模编排按内存预算切 XY 矩形、扩展带、块内掩膜、合并标签核心矩形列表、扩展矩形列表、块可用性掩膜工具链粗差剔除、空洞 TPS 插值、精度评估、可选图割变体与主路径解耦按需挂接主路径实现类为PTD 滤波器形态学为独立类但在种子阶段被主类内嵌调用非继承多态而是组合。2.2 地面标签的生命周期形态学阶段写入布尔地面标记数组可多尺度取并。PTD 阶段在标记数组上继续置真迭代结束后同步维护内点/外点索引列表供导出。分块阶段块内产生完整标记数组再与块核心区域可用性掩膜按位与写回全局结果仅核心区内点参与最终地面判定扩展带内点只服务于边界三角网稳定。实现上地面判定以布尔标记数组为准索引列表在滤波结束时由标记数组扫描生成。三、类型一平面格网索引3.1 为何需要格网PTD 每轮迭代要在“新接受地面点”影响区域内找未分类候选点形态学要在每个单元取代表高程通常最低 Z。若对全点云做 KD 树邻域搜索大规模点云 IO 与 cache 都不友好。实现选择固定边长的规则格网空间换时间且与形态学、分块对齐天然一致。3.2 离散化规则给定单元边长dBlkSize原点为(dMinX, dMinY)由点云四至向下取整或角点对齐分块时传入统一角点保证各块格网 I/J 一致。点(x,y)的格网索引idxI floor((x - dMinX) / dBlkSize)实现中 I/J 与 X/Y 对应关系以代码为准读写时勿混淆行列。容器为mapIdxIJ, vector点索引稀疏格网只存有数据的单元。角点对齐设置角点坐标与倍数关系后强制dMinX/dMinY对齐到统一网格避免分块边界处同一物理位置落入不同 I/J造成候选点漏检或重复处理。3.3 初始化与密度估计自动参数自动格网尺寸逻辑实现者在调参时需理解用试探性较细格网插入全点统计有效格网数有至少一个点的单元。平均密度 ≈ 总点数 / (有效格网数 × 单元面积)。若密度低于阈值认为点过于稀疏自动参数失败。目标“每格约默认点数如 6 个”反解边长dBlkSize ≈ sqrt(目标点数/密度)并钳制在[0.2m, 0.1×最大地物尺度]。若空洞格网占比 1%用形态学外圈填充估计空洞率按比率放大边长减少“碎格网”。山地判定在较粗格网边长取最大地物尺度上比较邻格高差超过约tan(5°)×地物尺度则标记为坡地供后续图割变体或坡度标记使用。3.4 栅格化接口形态学与 PTD 共用最低 Z 栅格每格有点则取 Z 最小者无点则标记 void形态学种子、PTD 边界扩展用。最高 Z 栅格用于特定分析或剔除异常低点。邻域查询给定中心格 (i,j) 与窗口半径 NN收集邻域内所有点索引或要求邻域内至少MinNum个点才返回虚拟种子赋 Z 时用。3.5 实现注意点边长过小 → 格网数量爆炸代码中有约 1e8 格网上限检查。插入时检查边长合法性若单格点数过多自动参数回退。OpenMP格网填充、邻域收集等循环可并行格网映射插入需注意线程安全实现中部分路径关闭 OMP 或使用 critical。四、类型二形态学初始地面种子4.1 算法思想形态学路径把点云投影为二维 DEM 格网每格一个高程在栅格上做开运算先膨胀后腐蚀结构元素尺寸与最大地物尺度相关理论上移除尺度小于结构元素的凸起地物保留近似“地面包络”。实现不是单次开运算结束而是填洞 → 去异常低点 → 迭代开运算 → 再插值 → 可选格网域“密化”与 PTD 同名但发生在栅格上。4.2 栅格构建流程与 PTD 相同格网插入点云。每格取最低 Z生成浮点高程矩阵无点格标记 void。Inpainting小窗口插值补全 void避免开运算在空洞处发散。TPS 多尺度插值可选路径按格网尺寸与比例参数对 void 区做薄板样条插值适合大面积无点格。RemoveLowerOutliers剔除相对邻域过低的高程异常实现用固定阈值量级防止坑洞点拉低开运算面。4.3 开运算迭代OpeningIteration结构元素边长像素≈最大地物尺度 / 格网边长取奇数。对高程格网膨胀邻域内取最大填平小坑。腐蚀邻域内取最小削掉小凸起。输出形态学 DEM uchar 标记有效地面格。新 void 掩膜被判定为非地面的格。迭代可重复自适应路径使用MultiScaleGrid在多个放大倍率的格网分辨率上构建层叠格网在粗层上判断地形特征平地/坡地再指导细层开运算参数——实现上把“多尺度”放在格网金字塔而非点云金字塔。4.4 从栅格回到点云标签开运算结束后对每个有效格网若该格为地面格则格内所有点标记为地面候选布尔数组置真。多尺度外层循环PTD 种子阶段对不同尺度因子如 1.0、2.5、5.0 倍边长重复调用形态学类结果逻辑或合并提高复杂地形召回率。扫描线分析大最大地物尺度时另一种子路径用多方向扫描线分析格网高程适合特定生产配置。4.5 与 PTD 的衔接PTD 主滤波器在滤波开头小最大地物尺度 → 自适应形态学同时得到坡地标记数组。否则 → 多尺度形态学或扫描线。形态学不维护三角网只输出初始布尔标记数组。PTD 随后重建格网边长可能与形态学阶段不同以 PTD 参数为准。五、类型三渐进三角网密化 PTD核心5.1 理论对应与实现选择PTDProgressive TIN Densification经典流程由稀疏地面种子构建 TIN在三角面片上对未分类点做坡度 高差检验通过则加入地面并更新 TIN。本实现三角网为OpenCV Subdiv2D平面 Delaunay顶点存于种子索引表前若干为虚拟边界点索引 0 表示虚拟点本身其后为真实点云索引或特殊编码。插入点在 XY 平面定位Z 参与生长判据。5.2 边界虚拟种子ExtendingSeeds问题Delaunay 在凸包外三角面片无数据边缘大量候选点无法被合理三角形覆盖易误判或停滞。实现边界虚拟种子扩展在点云平面四至矩形上按格网边长生成边界格网序列边距与最大地物、格网边长相关。在边界格网位置布置平面虚拟点对每点在邻域格内找最近真实种子点复制其Z避免虚拟点 Z 任意。虚拟点写入虚拟种子数组同时扩展 Delaunay 的 rect 边界整数像素矩形略大于点云范围。真实种子再插入三角网先虚拟点再所有布尔标记为真的点。虚拟点不参与最终地面导出索引表映射约定仅用于撑开TIN。5.3 三角网初始化与插入语义初始化 Delaunay 矩形与预留顶点数约为点数一半经验。首次插入遍历虚拟种子 → 遍历当前地面标记为真的点。后续迭代插入仅插入本轮新接受的地面点记录每条插入边在细分中的编号用于计算受影响矩形在格网 I/J 四至范围内传播。保持边引用一致失败插入则回滚种子索引表。5.4 单次迭代的三步DensifyingIteration实现为固定循环最大迭代次数可配置为 0 时仅输出种子步骤 1 — 插入新地面点到 TIN将上一轮接受的地面点依次插入 Subdiv2D对每点记录插入是否成功与退化、重复、边界有关。影响矩形沿三角网邻接边 walk更新 I/J 最小最大范围实现中存为连续四个整数 minI,minJ,maxI,maxJ。若本轮无新点且非首轮影响矩形为空 →迭代终止收敛。步骤 2 — 收集候选点FindCandidates对每个影响矩形格网四至范围枚举矩形内格网 → 汇总格内尚未分类的点索引布尔标记为假。去重后得到候选列表。这是 PTD局部性的关键不用全图扫描只扫受影响格网。步骤 3 — 密化判定Densifying对每个候选点在 XY 平面定位所在三角形Submotion locate取三角形三顶点坐标通过种子索引表映射到三维点虚拟点用虚拟数组。计算三角形平面方程参数axbyczd0归一化法向。对每个相关三角形实现中常处理与候选点关联的边相邻的三角面执行PointGrowing。PointGrowing 判据标准模式平面法向坡度预筛若|c| sin(坡度阈值)法向近水平该三角形不参与生长过陡坡面。点在三角形内XY 重心坐标或等价 inside 测试。竖直投影在三角形内VerticalPointInTriangle将点沿 Z 投影到平面投影点须在三角形 XY 投影内避免跨越垂直遮挡。IsGroundPoint点到平面距离|axbyczd|。点到三顶点最短欧氏距离的平方minDist2。判据距离/√minDist2 ≤ sin(坡度阈值)的等价形式实现用预计算 sin/cos 坡度阈值。另高差dfThr_DZ在平面距离中隐含使用。通过则标记该候选点为地面。记录其关联三角网边编号供下轮插入。每个三角形还带最低 Z 候选与当前最低点两个索引实现用于在三角面内优先生长最低点。边缘保持模式bPreserveEdge若标准判据失败查当前边邻接三个三角形的法向与连通性IsConnectedToFacet要求与至少一个邻接面片在法向夹角上“足够水平”cos 阈值。高程差不超过垂直容差。用于陡坎处避免把坎下点强行拉入缓坡三角面。并行按三角形分块 OpenMP每线程局部 map 合并边—新点映射。若本轮新地面点数为 0 → 迭代结束。5.5 可选后处理栅格插值bFillOuliers若开启填充且提供范围矢量在格网标记“空洞”单元无地面点、在 shp 范围内。用 TPS 等方式在空洞格网插值生成合成地面点存入插值点共享容器供 DSM 产品补洞。与形态学里的 TPS 插值类似但发生在分类完成之后。5.6 可选变体图割迭代FilterByGraphCut实现中存在实验/备用路径先构建全点 Delaunay在三角网边上算平滑项权重法向夹角相关顶点上算数据项权重与地面可能性相关用 OpenCV GCGraph 做 s-t 最小割迭代更新地面标记。生产主路径仍以坡度高差生长为准图割体现“全局一致性”思路与 PTD 局部生长形成对照。实现者扩展算法时可参考权重计算与坡地标记数组的联动。六、类型四大规模分块6.1 动机单块加载全点云时内存约正比于点数PTD 还需 Submotion、候选缓冲、临时 map。实现通过XY 平面分块控制峰值内存。6.2 块尺寸与扩展带由用户内存上限MB与平均密度估算核心块边长允许点数 ≈ 内存系数 × 1024² / 每点字节估算。块面积 ≈ 点数/密度边长为面积开方。扩展带宽度≈最大地物尺度 格网边长实现中两者相加。核心块为互不重叠除最后一行/列边界处理的矩形扩展块 核心块每边外扩扩展带宽度。扩展带内点参与该块 PTD三角网、候选搜索但写回全局结果时仅核心矩形内点保留地面标签可用性掩膜与结果按位与。6.3 块边界与索引一致性全图统一角点坐标传入各块 PTD 实例保证格网 I/J 对齐。最后一行/列块在点归属上用左闭右开或闭区间区分避免角点重复计数或遗漏实现中对 lastX/lastY 有四分支。6.4 块级执行流程对每块扩展矩形从 PLY/SHP按包围盒裁剪读点只加载扩展矩形内点。构造块内 PTD 对象若参数未完整则块内自动参数估计。块内关闭内层 OMP若外层块间并行已开启避免过度订阅 CPU。执行完整 PTD 滤波流程种子→格网→虚拟边界→迭代。将块内地面标记与核心掩膜合并追加写地面点/矢量文件分块输出模式。块间并行OpenMP parallel for 块编号写文件用 critical 保护。6.5 密集点云标志若估算块边长小于约4×扩展带置密集点云标志后续可调整策略警告日志——实现者对超密航高场景应预期分块更碎、内存更高。七、类型五辅助模块模块实现要点粗差剔除基于离散粗糙度等可在 PTD 前挂接主路径中默认注释掉空洞填补独立类TPS/格网插值与 PTD 后处理类似精度评估分块加载、对比参考地面服务算法调参形态学工具结构元素生成、栅格形态学、扫描线等供 Morpho 调用几何运算三角形内判定、平面参数、连通 facet 判据供 PTD 调用八、关键参数在实现中的含义参数实现绑定格网间距格网索引边长、形态学格网分辨率、候选搜索粒度自动估计时与密度挂钩最大地物尺度形态学结构元素尺寸、虚拟边界扩展宽度、分块扩展带主分量、自动格网上限钳制坡度阈值预计算 sin/cos三角形法向筛选与距离/水平距离比高差阈值平面距离绝对值边缘保持邻接面片高程差最大迭代次数密化循环上限为 0 则仅形态学种子保持边缘切换 PointGrowing 与 PointGrowing_PresEdge山地标志自动参数中坡地检测形态学自适应关联坡地标记分块开关入口选择 Massive 或 Simple 路径块内存上限反推核心块几何尺寸块间并行块间 OpenMP块内关闭 OMP填充开关是否 TPS 插值合成点需范围矢量九、并行与内存策略PTD 迭代三角形级 OpenMP合并 map 需注意。候选收集影响矩形数量大时矩形串行循环格网数量可控。分块块间并行 写盘 critical单块内倾向串行 OMP 避免嵌套。DEM 行优化若将来在类似 GDAL 变换中复用同一行点共享 DEM 缓冲可降 IO——PTD 本体是内存点云。Subdiv2D 插入大量退化边时插入失败会回滚种子索引日志需关注。十、与稀疏 DEM 产品线的衔接上游典型消费地面点点云PLY 或 GPKG 等。可选插值合成点填充空洞。布尔标记或分类字段写入矢量。它不生成 DEM 本身但地面点的 XY 分布与 Z 质量直接决定后续栅格化插值是否稳定。分块接缝处若角点未对齐会出现条带状地面遗漏扩展带宽度小于地物尺度时坎边三角网不稳定。十一、五类任务与主流程对照实现者速查点云输入 → [可选] 粗差剔除 → 形态学格网最低Z → 开运算 → 点标记种子 → 格网索引 Insert 全点 → 边界虚拟种子 Delaunay 初始化 → 循环插入新地面 → 影响矩形 → 候选点 → 三角面生长判据 → [可选] 空洞 TPS 插值点 → 输出地面标记 / 导出类型核心数据结构核心数值操作格网索引格网→点索引映射离散化、邻域枚举、最低 Z 栅格形态学种子高程格网 uchar/float开运算、插值、格网→点 ORPTD 密化Submotion 种子索引表平面定位、坡度/高差、迭代插入大规模分块核心/扩展矩形列表裁剪读点、掩膜写回、角点对齐辅助工具—几何判定、评估、填洞十二、实现扩展与排错建议种子过少查形态学 void 比例、开运算结构元素是否过大、格网边长是否过大。迭代不收敛/极慢影响矩形是否过大单轮接受点过多最大迭代次数三角网退化。块接缝角点、扩展带、lastX/lastY 归属逻辑写回是否只用核心掩膜。陡坡漏分坡度阈值、边缘保持是否关闭、山地自动参数是否生效。城区建筑残留减小最大地物尺度或格网边长加强形态学多尺度。十三、总结它的实现本质在规则格网上用形态学开运算得到初始地面种子在 OpenCV 平面 Delaunay 三角网上用坡度与高差约束迭代生长地面点并用格网索引把每轮搜索限制在局部区域超大数据通过分块与扩展带在内存可控下重复同一套逻辑。扩展算法应优先改种子质量、生长判据、虚拟边界与分块几何四者而非仅替换对外接口。
