基于车载LiDAR点云数据的街景面片提取及优化方法

基于车载LiDAR点云数据的街景面片提取及优化方法
【专利摘要】本发明公开了一种基于车载LiDAR点云数据的街景面片提取及优化方法，包括以下步骤：LiDAR点云数据的面片提取，将LiDAR点云数据划分成若干个小立方体块，并将各个小立方体块拟合，得到符合预设要求的面片；对面片进行合并，合并后的面片进行面片增长，对增长后的面片再进行合并，最后对得到的面片进行几何纠正，并消除噪声面片。本发明实用性广，即使点云数据质量较差也能得到较好的结果。
【专利说明】基于车载L i DAR点云数据的街景面片提取及优化方法

【技术领域】
[0001] 本发明涉及车载LiDAR点云数据处理【技术领域】，特别是涉及一种基于车载LiDAR 点云数据的街景面片提取及优化方法。

【背景技术】
[0002] 随着街景地图的发展，实景对象的面片信息提取相关问题是当前亟待解决的问 题。车载激光扫描系统的发展使得实景对象的立面提取越来越多地使用车载LiDAR数据， 但利用车载LiDAR数据时仍然存在几个主要问题：点云数据处理时间长，计算量大；场景复 杂，目标丰富；建筑物立面层次结构复杂，立面细节特征丰富。针对这些问题的存在，基于点 云数据的面片优化是关键步骤之一。
[0003] 目前三维点云平面检测常用的方法有三种，基于随机抽样一致性算法（RANdom SAmple Consensus, RANSAC)、基于霍夫变换的算法、区域增长法，基于这三种方法也出现了 很多改进算法。
[0004] 基于随机抽样一致性算法：是一种模型参数估计方法，将其应用到点云的平面检 测，并在随机采样过程中引入最小采样半径来限制采样范围，解决了外点密度和分布变化 带来的模型估计不稳定问题，并且采用随机子集对估计到的模型进行评价，减少了运算量， 在保持形状模型检测精度的同时提高了计算速度。
[0005] 基于霍夫变换的方法：首先对参数空间（霍夫空间）离散化得到一系列子空间，即 累加器，然后对每个累加器投票，对票数进行统计，峰值累加器所对应的图像空间中的平面 即为检测到的最佳平面。然而累加器的设计存在共同的缺陷，即累加器的大小不一致和峰 值可辨性差，导致检测到的平面模型不准确，而且此方法耗占较大的存储空间。
[0006] 区域增长法：是根据同一物体区域内象素的相似性质来聚集象素点的方法，从初 始区域开始，将相邻的具有同样性质的象素或其它区域归并到目前的区域中从而逐步增长 区域，直至没有可以归并的点或其它小区域为止。区域生长方法是一种比较普遍的方法，在 没有先验知识可以利用时，能取得最佳的性能，可以用来分割比较复杂的图像，但是，区域 增长是一种迭代的方法，空间和时间开销都比较大。
[0007] 以上三种方法在大数据时内存开销大、计算时间长，对数据质量要求较高。


【发明内容】

[0008] 本发明要解决的技术问题在于针对现有技术中的上述缺陷，提供一种充分利用点 云空间分布特征来优化面片，降低错误率，提高精度的基于车载LiDAR点云数据的街景面 片提取及优化方法。
[0009] 本发明解决其技术问题所采用的技术方案是：
[0010] 提供一种基于车载LiDAR点云数据的街景面片提取及优化方法，其特征在于，包 括以下步骤：
[0011] 步骤I =LiDAR点云数据的面片提取，基于八叉树结构，将LiDAR点云数据划分成若 干个小立方体块，并将各个小立方体块拟合，得到符合预设要求的面片；
[0012] 步骤2 :对步骤1中得到的面片进行合并，对每一个面片，构建邻接表，计算所有相 邻面片法向量之间的夹角，若某相邻的面片P1和P 2的夹角Θ〈 Θ th，则将其进行合并，得到 新的面片P'，计算合并后面片的面精度的下降值，若符合事先给定的阈值范围，则合并面 片P 1和P2，得到新面片P，并更新邻接表；重复步骤2,直至没有可以再进行合并的面片；
[0013] 步骤3 :对步骤2中合并后的面片进行面片增长，将步骤2中得到的面片视为种子 面片，以最小单元为增长单元，26邻域进行增长；对增长后的面片再按照步骤2的方法进行 合并，得到最后合并的面片；
[0014] 步骤4 :对步骤3中得到的面片进行几何纠正，并消除噪声面片。
[0015] 本发明所述的方法中，步骤4中对面片进行几何纠正具体包括：
[0016] 根据面片包含的点，求得面片的最小包围矩形框的角点坐标；
[0017] 对交叉的最小包围矩形框，在交叉处对边界进行进一步的优化，具体为将所有矩 形框投影到XOY平面上，在二维空间中，对每条线段两端进行扩展，扩展长度为1 ;判断扩展 后的线段是否相交，若相交，则计算相交线段L1和L2对应的矩形框R1和R 2在Z轴方向的距 离h，公式如下：
[0018] h= (hlmax-h2max+hlmin_h2min)/2
[0019] 其中hl_，hlnin和h2_，I^niin分别为对应的矩形框R 1和R2中四个角点中z坐标 最大和最小值；若h小于给定的阈值hth，则求得相交点坐标，并执行下一步，否则，继续判 断下一对相交线段；分别计算交点坐标到两条线段的端点的水平距离，记为dl begin，dlend和 d2begin，d2end，并求得 dl- = min(dlbegin, dlend)和(12- = min(d2begin, d2end)。分别判断(Ilniin 和d2min是否小于给定距离阈值dth，若小于，则结合交点坐标信息，切除交点到相应端点之 间的矩形框。
[0020] 本发明所述的方法中，步骤4中，消除噪声面片包括以下步骤：
[0021] 获取每个站点对应的点云和面片信息；
[0022] 利用移动测量车记录下的每个站点在世界坐标系下的三维坐标，将预先获取的电 力线点云X wiM按照距离分配到每一个站点上，得到每个站点Si上对应的电力线点云XLre;
[0023] 对每一个站点Si上的面片集合Pi中的所有面片进行分析；如果面片= 1，2, 对应的点云分布较为稀疏，或者点云的数量小于某一阈值δ th，则认为该面片可 能为一个电力线片面；
[0024] 在Si对应的电力线点云；4α中进行搜索，得到出在面片If范围内且到If平面的 距离小于〇1米的点的个数NwiM，再计算出面片的单位宽度范围内平均的点云个数\ = \&/ w，w表示面片的宽度，如果NW>D2，则认为该面片是电力线面片，将其删除。
[0025] 本发明所述的方法中，消除噪声面片还包括对矩形框进行合并，针对之前的面片 合并做进一步的合并优化，基于面片的几何关系，优化合并结果，具体包括以下步骤：
[0026] 将所有的矩形框按照其面积进行排序得到排序之后的矩形框集合PSOTt，然后从 PSOTt中依次找出当前最大的矩形框Pirt;
[0027] 对吐⑴找到其对应的站点Si,将氏-3, SJ3]范围内的站点的矩形框作为当前 Ps'ort 附近的矩形框；
[0028] 将这些矩形框按照其与之间的位置关系分成水平方向和坚直方向的两类，对 与之间平面夹角小于Θ，且距离小于D的面片进行分析：首先对水平方向的矩形框进 行合并分析，先将待分析旋转到与当前矩形框巧。#平行，然后将旋转之后的矩形框投 影到Pirt平面上，得到一个二维矩形R hyp，计算出该矩形与PLrt四个顶点在该平面上构成的 矩形Rcm之间的重叠区域R tjl以及两者的最小包围矩形Rdiot，如果R。）λ Jhyp则将P/^删除； 如果R_> λ 2 (Rc^Rhyp)，则将以_合并到中，取矩形Rdiot作为合并之后的矩形框（λ P λ 2为设定的参数）；对于坚直方向的矩形框，只要其旋转之后在沿面片Pirt宽的方向上不 超过:APif5rt (λ为设定的参数），则将其与矩形框Picrt合并，取两者的最小包围矩形作为合 并之后的矩形框。
[0029] 本发明所述的方法中，步骤3中，面片增长具体包括步骤：
[0030] 计算邻接最小单元中所有点（个数为Μ)到种子面片的距离d，统计距离d小于给 定阈值d th的点的个数，记为N ;若N大于给定的阈值Nth并且大于rXM(0〈r〈l)，r为给定 的参数，则进行增长，同时重新估计面片参数，并更新邻接表；否则，不增长；并进行循环迭 代，直至所有面片处理完；
[0031] 逐步增大dth的值，增大幅度为d'，并进行循环迭代，直至dth小于给定的最大阈 值d' th，则停止迭代。
[0032] 本发明所述的方法中，步骤2中面精度的下降值计算和根据阈值范围进行面片合 并，具体包括以下步骤：
[0033] 若某相邻的面片P1和P2的夹角θ〈 Θ th，则将其进行合并，采用RASNAC算法重新 估计得到新的面P，并计算面P的方差，记为σ';分别记面片？1和己的方差为〇1和 σ 2，并记σ ^ = 2)。通过方差，计算合并后面精度的下降值Od = σ' -σ·。 若〇d〈max(3*〇Mx，〇th)，其中〇 th为事先给定的阈值，则合并面hjPP2，得到新面P'，并 更新邻接表；若。d>max(3· σ_，〇th)，则不进行合并操作。
[0034] 本发明产生的有益效果是：本发明通过基于车载LiDAR点云数据提取街景面片， 并优化结果，包括对面片矩形框的几何纠正，如交叉，倾斜等情况；对噪声面片进行消除，如 噪声点的干扰，电力线面片等，完成这些操作，实验证明即使点云数据质量较差也能得到较 好的结果。

【专利附图】

【附图说明】
[0035] 下面将结合附图及实施例对本发明作进一步说明，附图中：
[0036] 图1是本发明实施例点云建筑物面片提取及优化流程图；
[0037] 图2是本发明实施例基于八叉树的面片提取、面片增长及面片合并示意图；
[0038] 图3是本发明实施例交叉面片处理示意图。

【具体实施方式】
[0039] 为了使本发明的目的、技术方案及优点更加清楚明白，以下结合附图及实施例，对 本发明进行进一步详细说明。应当理解，此处所描述的具体实施例仅用以解释本发明，并不 用于限定本发明。
[0040] 本发明实施例的基于车载LiDAR点云数据的街景面片提取及优化方法，如图1所 示，包括以下步骤：
[0041] 步骤1，LiDAR点云数据的面片提取，共包括初始面片提取、面片合并、面片增长三 个步骤。
[0042] 步骤1. 1，初始面片提取，目前提取面片的常用方法包括，基于随机抽样一致性算 法（RANdom SAmple Consensus, RANSAC)、基于霍夫变换的算法、区域增长法，以及基于这三 种方法提出的很多改进算法。本发明面片提取过程如下：
[0043] 步骤1. 1. 1，基于八叉树结构来处理数据，假设三维点云数据最大包围立方盒为 B，根据该立方盒的大小将其切分成若干个大小为MXM(M = 2N)的小立方体，记为Bi，此处 即首先分成8个立方体，考虑到立方体太大无法拟合成面片，并且拟合的面片也不准确，故 根据先验知识将立方体继续划分η层，得到合适大小的立方体
[0044] 步骤I. 1. 2,针对每个小立方体Bi,采用最小二乘的方法拟合得到面Pi,并计算得 到每个点到面P i的距离Clij，同时求得最大距离dmax。
[0045] 若dmax〈dth，hth为给定的阈值，则面P i符合要求；
[0046] 若d_>dth，且d_〈 λ dth，λ为预先设置的参数，则采用RANSAC算法对该立方体进 行重新拟合，得到新的面Pi;
[0047] 若d_>dth，则对该立方体按照八叉树结构进行均匀分割，得到8块新的立方体，重 复迭代此步骤；
[0048] 步骤I. 1. 3,对于步骤I. 1. 2中的前两种情况，即dmax〈dth或dmax〈 λ dth，计算该立方 体中包含的点云的曲率平均值Rniean。
[0049] 若R_〈Rth，Rth给定的阈值，则保留该立方体；
[0050] 若R_>Rth，则对该立方体进行均匀分割，得到8块新的立方体，重复迭代步骤 I. 1. 2 和步骤 I. 1. 3。
[0051] 步骤I. 1. 4,直至每个立方体都被拟合成面或者已经为最小单元（大小为W)，则停 止迭代。
[0052] 步骤1.2,面片合并。
[0053] 步骤1. 2. 1，对每一个面片，构建邻接表。
[0054] 步骤1. 2. 2,根据邻接表，计算所有相邻面片法向量之间的夹角，并将所有夹 角从小到大进行排序，同时，将组成夹角的两个面片依次排序。设置阈值Q th，若某相 邻的面片P1和P2的夹角θ〈 Θ th，则将其进行合并，采用RASNAC算法重新估计得到新 的面W，并计算面W的方差，记为σ '。分别记面片？1和己的方差为〇:和σ2，并 记σ MX = max( σ σ 2)。通过方差，计算合并后面精度的下降值σ d = σ ' -σ _。若 σ d〈max(3 · σ _，σ th)，其中σ th为事先给定的阈值，则合并面P1和P2，得到新面并更 新邻接表；若。 d>max(3·〇th)，则不进行合并操作。
[0055] 步骤1. 2. 3,重复1. 2. 2步骤，直至没有可以再进行合并的面片。
[0056] 步骤I. 3,面片增长（如图2所示）
[0057] 步骤1. 3. 1，对步骤1. 2中得到的所有初始面片按照大小（包含最小单元的个数） 从大到小进行排序。
[0058] 步骤1.3. 2,针对每个初始面片，将其视为种子面片，以最小单元为增长单元，26 邻域进行增长。计算邻接最小单元中所有点（个数为M)到种子面片的距离d，统计距离d 小于给定阈值dth的点的个数，记为N。若N大于给定的阈值Nth并且大于rXM(0〈r〈l)，r 为给定的参数，则进行增长，同时重新估计面片参数，并更新邻接表；否则，不增长。对步骤 1. 3. 2进行循环迭代，直至所有面片处理完。
[0059] 步骤1.3. 3,逐步增大dth的值，增大幅度为d'，并重复步骤1.3. 2,直至dth小于 给定的最大阈值d' th，则停止迭代。
[0060] 步骤1. 3. 4,对增长后的面片执行步骤1. 2,如图2所示，即进行合并，得到合并后 的面片。
[0061] 步骤2,面片后处理，共包括面片粗滤除、面片矩形框获取及精细化、交叉矩形框优 化、面片噪声去除四个步骤。
[0062] 步骤2. 1，面片粗滤除。本发明所需的面片为街景坚立面，而步骤1所得的面片中， 存在较多的较小的面片以及水平面片，需要进行滤除。假设所有的面片集合为P i,针对每个 面片，计算其宽度W，高度H，面积S，点的个数N及法向量与Z轴夹角A，若符合下式条件，则 保留；否则，滤除。
[0063] W>Wth, H>Hth, S>Sth, A>Ath,
[0064] 其中Wth，Hth，Sth，R th及Ath为给定的合理阈值。
[0065] 步骤2. 2,面片矩形框获取。本发明实施例中，面片最终用矩形框进行表达，因此， 需要根据面片包含的点，求得其最小包围矩形框，具体过程如下：
[0066] 步骤2.2. 1，针对每一个面P，其法向量为（!^!^!^，假设其包含的点云集三维坐 标分别为X = {々於，Y = {yiW和Z= {Zjf，其中N为点的个数。首先将所有点投影到面p， 得到新的点集X'=以片，Y' = 和Z' = ?於，设点（Xp，ν Zp)为面p上任意一点，则公 式如下：
[0067] X' j = Xj+nxXD, y' t = Y^nyXD, z' i = Z^nzXD,
[0068] 其中D为点（Xi，yi，Zi)到面P的距离，其计算公式如下：
[0069] D = nx (Xp-Xi) +ny (yp-y) +nz (Zp-Zi)
[0070] 步骤2. 2. 2,绕Z轴旋转，使面法向量在XOY平面上的投影分量和Y轴平行，并求得 3X3旋转矩阵Rz。
[0071] 步骤2. 2. 3,基于该旋转矩阵，求取旋转后的点集y = ，Y7"=况垃:和 f = ，计算公式如下：
[0072] (x[,yf,z[) = (x-,yl，z·) x Rz
[0073] 步骤 2· 2· 4,求取点集（Xlr, Ylr, ZO 中的最大最小点 Oinax,yiiax, Zjnax)， (Xmin，ymin，zmin)，其中 ymax * Υη?η，设 Ymean - (ymax + ymin)/2? 故矩形框的四个 角点坐标分力lj 为（义min，ymean，Zmin) j (^-max'yinean? ^min) J Ymearp Zmax) J (-^max? yinean? ^max) °
[0074] 步骤2. 2. 5,将四个顶点旋转到原坐标系中，得到最终的矩形框角点坐标，旋转矩 阵为Rz的逆矩阵V z。
[0075] 步骤2. 3,矩形框精细化。在步骤2. 2中，由于所求的矩形框为最小包围盒，很多地 方没有点云，却处于矩形框的范围内，因此，需要进行进一步的精细化，即拆分矩形框，具体 过程如下：
[0076] 步骤2.3. 1，针对每一个面片，进行步骤2.2. 1，步骤2. 2. 2和步骤2. 2. 3,得到旋 转后的点集。沿X轴方向，对点集进行切块，切块宽度为Wx，并分别求得每一块中点最大的 z坐标。
[0077] 步骤2. 3. 2,对切块进行聚合，即考虑相邻切块最大z坐标差值的绝对值，若小于 给定的阈值zth，则合并切块。当合并完成后，将宽度小于W th的切块合并到邻接切块中。
[0078] 步骤2. 3. 3,针对每一个切块，执行步骤2. 2. 4和步骤2. 2. 5,完成矩形框的精细 化，得到所有矩形框i? = {&}?，η为矩形框个数。
[0079] 步骤2. 4,交叉矩形框优化。在获得的所有矩形框中，存在矩形框交叉的问题，为 了保证矩形框的紧密性以及合理性，在交叉处需要对边界进行进一步的优化，具体过程如 下：
[0080] 步骤2. 4. 1，如图3所示，将所有矩形框投影到XOY平面上，即将矩形框投影为一系 列的线段X =仏?}?，η为线段个数。在二维空间中，对每条线段两端进行扩展，扩展长度为 1〇
[0081] 步骤2. 4. 2,判断扩展后的线段是否相交，若相交，则计算相交线段L1和L2对应的 矩形框R 1和R2在Z轴方向的距离h，公式如下：
[0082] h= (hlmax-h2max+hlmin_h2min)/2
[0083] 其中hl_，hlnin和h2_，I^niin分别为对应的矩形框R 1和R2中四个角点中z坐标 最大和最小值。若h小于给定的阈值hth，则求得相交点坐标，并执行下一步，否则，继续判 断下一对相交线段。
[0084] 步骤2. 4. 3,分别计算交点坐标到两条线段（假设SLdPL2)的四个端点的 水平距离，记为 dlbegin，dlend 和 d2begin，d2end，并求得(Ilniin = min(dlbegin, dlend)和（12-= min(d2begin，d2end)。分别判断dlmin和d2 min是否小于给定距离阈值dth，若小于，则结合交点 坐标信息，切除交点到相应端点之间的矩形框。
[0085] 步骤3,面片噪声去除，包括电力线面去除和面片合并，具体过程如下：
[0086] 步骤3. 1，获取每个站点对应的点云和面片信息。
[0087] 利用移动测量车记录下的每个站点在世界坐标系下的三维坐标，将电力线点云 Xwire按照距离分配到每一个站点上，得到每个站点Si上对应的电力线点云；
[0088] 步骤3. 2,电力线面片去除。
[0089] 对每一个站点Si上的面片集合Pi中的所有面片进行分析。如果面片= 1，2, ..·，?}，对应的点云分布较为稀疏，或者点云的数量小于某一阈值δ th，则认为该面片可 能为一个电力线片面。接着在Si对应的电力线点云中进行搜索，得到出在面片If范 围内且到Pf平面的距离小于〇1米的点的个数Nwire，再计算出面片的单位宽度范围内平均的 点云个数\ = NwiM/w，w表示面片的宽度，如果NW>D2，则认为该面片是电力线面片，将其删 除。
[0090] 步骤3. 3,矩形框合并，针对步骤1. 2的面片合并做进一步的合并优化，基于面片 的几何关系，优化合并结果，具体过程如下：
[0091] 步骤3. 3. 1，将所有的矩形框按照其面积进行排序得到排序之后的矩形框集合 Pstjrt，然后从Pstjrt中依次找出当前最大的矩形框。
[0092] 步骤3. 3. 2,对巧^，找到其对应的站点Si,将氏-3, SJ3]范围内的站点的矩形框 作为当前附近的矩形框。
[0093] 步骤3. 3. 3,将这些矩形框按照其与:之间的位置关系分成水平方向和坚直方 向的两类。只需对与之间平面夹角小于Θ，且距离小于D的面片进行分析：首先对 水平方向的矩形框进行合并分析，先将待分析旋转到与当前矩形框Pj ert平行，然后将 旋转之后的矩形框投影到Pirt平面上，得到一个二维矩形Rhyp，计算出该矩形与Pia tI四个 顶点在该平面上构成的矩形Rcm之间的重叠区域Rtjl以及两者的最小包围矩形，如果 Rol> λ Ayp则将P/〇rt删除；如果D λ 2 (Rcur+Rhyp)，则将p/〇rt合并到Pjort中，取矩形R mer作为 合并之后的矩形框（λπ λ2为设定的参数）。对于坚直方向的矩形框，只要其旋转之后在 沿面片Piert宽的方向上不超过APi ert (λ为设定的参数），则将其与矩形框Piert合并，取两 者的最小包围矩形作为合并之后的矩形框。
[0094] 本发明综合现有技术，基于车载LiDAR点云数据提取街景面片，并优化结果，包括 对面片矩形框的几何纠正，如交叉，倾斜等情况；对噪声面片进行消除，如噪声点的干扰，电 力线面片等，完成这些操作，实验证明即使点云数据质量较差也能得到较好的结果。
[0095] 应当理解的是，对本领域普通技术人员来说，可以根据上述说明加以改进或变换， 而所有这些改进和变换都应属于本发明所附权利要求的保护范围。
【权利要求】
1. 一种基于车载LiDAR点云数据的街景面片提取及优化方法，其特征在于，包括以下 步骤： 步骤I=LiDAR点云数据的面片提取，将LiDAR点云数据划分成若干个小立方体块，并将 各个小立方体块拟合，得到符合预设要求的面片； 步骤2:对步骤1中得到的面片进行合并，对每一个面片，构建邻接表，计算所有相邻面 片法向量之间的夹角，若某相邻的面片P1和P2的夹角Θ〈Θth，则将其进行合并，得到新的 面片Pi，计算合并后面片的面精度的下降值，若符合事先给定的阈值范围，则合并面片P1 和己，得到新面片P，并更新邻接表；重复步骤2,直至没有可以再进行合并的面片； 步骤3 :对步骤2中合并后的面片进行面片增长，将步骤2中得到的面片视为种子面 片，以最小单元为增长单元，26邻域进行增长；对增长后的面片再按照步骤2的方法进行合 并，得到最后合并的面片； 步骤4 :对步骤3中得到的面片进行几何纠正，并消除噪声面片。
2. 根据权利要求1所述的方法，其特征在于，步骤4中对面片进行几何纠正具体包括： 根据面片包含的点，求得面片的最小包围矩形框的角点坐标； 对交叉的最小包围矩形框，在交叉处对边界进行进一步的优化，具体为将所有矩形框 投影到XOY平面上，在二维空间中，对每条线段两端进行扩展，扩展长度为1 ;判断扩展后的 线段是否相交，若相交，则计算相交线段L1和L2对应的矩形框R1和R2在Z轴方向的距离h， 公式如下： h= (hlmax-h2max+hlmin-h2min)/2 其中hlmax，hlmin和h2max，h2min分别为对应的矩形框R1和R2中四个角点中z坐标最大 和最小值；若h小于给定的阈值hth，则求得相交点坐标，并执行下一步，否则，继续判断下一 对相交线段；分别计算交点坐标到两条线段的端点的水平距离，记为dlbegin，dlend和d2begin， d2end，并求得dl- =min(dlbegin,dlend)和(12- =min(d2begin,d2end)。分别判断dl-和 (12- 是否小于给定距离阈值dth，若小于，则结合交点坐标信息，切除交点到相应端点之间的矩形 框。
3. 根据权利要求1所述的方法，其特征在于，步骤4中，消除噪声面片包括以下步骤： 获取每个站点对应的点云和面片信息； 利用移动测量车记录下的每个站点在世界坐标系下的三维坐标，将预先获取的电力线 点云XwiM按照距离分配到每一个站点上，得到每个站点Si上对应的电力线点云 对每一个站点Si上的面片集合Pi中的所有面片进行分析；如果面八{/f,/C= 1，2,…，i}，对应的点云分布较为稀疏，或者点云的数量小于某一阈值δth，则认为该面片可 能为一个电力线片面； 在Si对应的电力线点云中进行搜索，得到出在面片If范围内且到If平面的距离 小于D1米的点的个数Nwira，再计算出面片的单位宽度范围内平均的点云个数Nw =NwiM/w，w表示面片的宽度，如果NW>D2，则认为该面片是电力线面片，将其删除。
4. 根据权利要求1所述的方法，其特征在于，消除噪声面片还包括对矩形框进行合并， 具体包括以下步骤： 将所有的矩形框按照其面积进行排序得到排序之后的矩形框集合Pstjrt，然后从Pstjrt中 依次找出当前最大的矩形框Piat; 对'Pirt，找到其对应的站点Si,将[Si-Sji+S]范围内的站点的矩形框作为当前PirtB近的矩形框； 将这些矩形框按照其与之间的位置关系分成水平方向和坚直方向的两类，对与 之间平面夹角小于Θ，且距离小于D的面片进行分析：首先对水平方向的矩形框进行 合并分析，先将待分析Airt旋转到与当前矩形框Pirt平行，然后将旋转之后的矩形框投影 到Pirt平面上，得到一个二维矩形Rhyp，计算出该矩形与Pirt四个顶点在该平面上构成的 矩形Rcm之间的重叠区域Rtjl以及两者的最小包围矩形Rdiot，如果IVλJhyp则将PsU删除； 如果R_>λ2 (Rc^Rhyp)，则将合并到PsL中，取矩形Rdiot作为合并之后的矩形框（λi， λ2为设定的参数）；对于坚直方向的矩形框，只要其旋转之后在沿面片宽的方向上不 超过(λ为设定的参数），则将其与矩形框UPjfjrt合并，取两者的最小包围矩形作为合 并之后的矩形框。
5. 根据权利要求1所述的方法，其特征在于，步骤3中，面片增长具体包括步骤： 计算邻接最小单元中所有点（个数为Μ)到种子面片的距离d，统计距离d小于给定阈 值dth的点的个数，记为N;若N大于给定的阈值Nth并且大于rXM(0〈r〈I)，r为给定的参数， 则进行增长，同时重新估计面片参数，并更新邻接表；否则，不增长；并进行循环迭代，直至 所有面片处理完； 逐步增大dth的值，增大幅度为d'，并进行循环迭代，直至dth小于给定的最大阈值d'th，则停止迭代。
6. 根据权利要求1所述的方法，其特征在于，步骤2中面精度的下降值计算和根据阈值 范围进行面片合并，具体包括以下步骤： 若某相邻的面片P1和P2的夹角θ〈Θth，则将其进行合并，采用RASNAC算法重新估计 得到新的面P，并计算面P的方差，记为σ';分别记面片？1和己的方差为〇1和〇2， 并记σ_ =max(σσ2);通过方差，计算合并后面精度的下降值σd =σ'-σ_ ;若 σ/ηκχ(3·σ_，〇th)，其中〇th为事先给定的阈值，则合并面PjPP2，得到新面P，并更 新邻接表；若。d>max(3·〇th)，则不进行合并操作。
【文档编号】G06T5/00GK104463871SQ201410751955
【公开日】2015年3月25日 申请日期:2014年12月10日 优先权日:2014年12月10日
【发明者】姚剑, 李礼, 鲁小虎, 陈梦怡 申请人:武汉大学