CN110517193B - 一种海底声呐点云数据处理方法 - Google Patents
一种海底声呐点云数据处理方法 Download PDFInfo
- Publication number
- CN110517193B CN110517193B CN201910576810.4A CN201910576810A CN110517193B CN 110517193 B CN110517193 B CN 110517193B CN 201910576810 A CN201910576810 A CN 201910576810A CN 110517193 B CN110517193 B CN 110517193B
- Authority
- CN
- China
- Prior art keywords
- point
- points
- point cloud
- neighborhood
- normal vector
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
- 238000003672 processing method Methods 0.000 title claims abstract description 8
- 238000000034 method Methods 0.000 claims abstract description 51
- 238000001914 filtration Methods 0.000 claims abstract description 18
- 239000013598 vector Substances 0.000 claims description 63
- 238000005070 sampling Methods 0.000 claims description 39
- 238000012545 processing Methods 0.000 claims description 17
- 230000005484 gravity Effects 0.000 claims description 10
- 230000008859 change Effects 0.000 claims description 8
- 239000011159 matrix material Substances 0.000 claims description 8
- 230000008569 process Effects 0.000 claims description 8
- 230000008439 repair process Effects 0.000 claims description 6
- 238000004458 analytical method Methods 0.000 claims description 5
- 230000000694 effects Effects 0.000 claims description 5
- 230000002146 bilateral effect Effects 0.000 claims description 4
- 238000004364 calculation method Methods 0.000 claims description 4
- 238000010276 construction Methods 0.000 claims description 4
- 230000009191 jumping Effects 0.000 claims description 4
- 238000009499 grossing Methods 0.000 claims description 3
- 238000000638 solvent extraction Methods 0.000 claims description 3
- 241000764238 Isis Species 0.000 claims description 2
- 238000012937 correction Methods 0.000 claims description 2
- 230000007547 defect Effects 0.000 claims description 2
- 238000005457 optimization Methods 0.000 claims description 2
- 238000000354 decomposition reaction Methods 0.000 claims 1
- 238000000265 homogenisation Methods 0.000 claims 1
- 238000001035 drying Methods 0.000 abstract description 4
- 230000004927 fusion Effects 0.000 abstract description 2
- 238000005192 partition Methods 0.000 description 5
- 230000014759 maintenance of location Effects 0.000 description 3
- 238000005259 measurement Methods 0.000 description 3
- 230000009467 reduction Effects 0.000 description 3
- 230000011218 segmentation Effects 0.000 description 3
- 238000012876 topography Methods 0.000 description 3
- 238000012423 maintenance Methods 0.000 description 2
- 238000012935 Averaging Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000005056 compaction Methods 0.000 description 1
- 238000013500 data storage Methods 0.000 description 1
- 238000001514 detection method Methods 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 238000000926 separation method Methods 0.000 description 1
- 238000010408 sweeping Methods 0.000 description 1
- 230000001131 transforming effect Effects 0.000 description 1
- 230000000007 visual effect Effects 0.000 description 1
- 238000013316 zoning Methods 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T5/00—Image enhancement or restoration
- G06T5/70—Denoising; Smoothing
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/10—Image acquisition modality
- G06T2207/10028—Range image; Depth image; 3D point clouds
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/20—Special algorithmic details
- G06T2207/20024—Filtering details
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Theoretical Computer Science (AREA)
- Measurement Of Velocity Or Position Using Acoustic Or Ultrasonic Waves (AREA)
Abstract
本发明公开了一种海底声呐点云数据处理方法,首先根据声呐点云的特性进行点云数据的融合与去燥,然后针对声呐点云的几何特性进行分区的点云精简,最后利用基于同心圆局部扩张的点云孔洞修补算法进行四维海底声呐点云孔洞的修补的海底点云数据处理方法;经本发明相关方法处理后的点云数据便于滤波算法对声呐点云数据中的人工目标进行识别与检测,既节省了人工对声呐图像数据判读所消耗的时间,又提高了人工目标识别算法的精确性,具有十分重大的意义。
Description
技术领域
本发明属于点云图像处理技术领域,具体涉及一种海底声呐点云数据处理方法。
背景技术
我国目前用于探测海底的装置为海底侧扫声呐,通过声呐侧扫的方式获得了大量的侧扫声呐扫描图像。侧扫声呐装置不但可以利用影像直观地描述海底人工目标的几何特性与海底形貌状况,而且通过对声呐数据三维坐标的分析可以对海底人工目标进行准确的定位,但如何有效利用声呐数据实现对海底人工目标准确和快速的判读是首先要解决的难题。当前针对海底声呐数据的判读主要依靠判读人员通过目视进行,然而随着大规模扫海工作的推进,会产生海量的待判读声呐数据。判读人员通常要对这类海量声呐数据进行回放识别,这极大地增加了判读人员的工作量。目前海底声呐图像主要包括两类,一类是经侧扫声呐对海底进行扫描探测所获得的包含经度、纬度和回声强度的三维点云图像,另一类是经原始点云图像进行解析变换形成的二维影像数据。二维声呐图像由于受水下环境和声呐自身特性的影响,生成的图像质量较差,噪声干扰严重,细节模糊,与光学图像存在着很大的区别,单纯针对声呐二维影像处理不易于海底目标的检测。
发明内容
本发明的目的是提供一种海底声呐点云数据处理方法,解决了现有技术中存在的人工对声呐图像数据判读,准确性不高的问题。
本发明所采用的技术方案是,首先,针对侧扫声呐数据与多波束测深***测深数据的特点进行了分析,将两种数据进行融合,形成四维点云数据。侧扫声呐在测量过程中由于多种因素会存在误差,受误差影响令原始点云数据存在很多噪声点。本发明将噪声分为离群的高频噪声和靠近点云主体的低频噪声两类,针对不同的噪声分别应用不同方法进行处理:应用基于近邻点平均距离的滤波算法去除高频噪声,对低频噪声应用基于双边滤波的各向异性光顺去噪算法对低频噪声进行滤波。
其次,融合后海底声呐点云数据点总量为百万级,且点与点之间密度分布不均匀。庞大且分布不均匀的采样点集增加了数据存储和处理的难度。因此,原始海底声呐点云数据的简化处理亟待解决。整块点云数据由平坦区域和高低起伏的几何特征不同的区域组成,并且几何特征明显的区域点云数据极有可能蕴含着海底人工目标。针对海底声呐点云数据的采样点数量庞大且分布不均匀、特征点多的特性,本发明提出了一种基于海底声呐点云特征的分区精简方法。首先,使用kd-tree进行点的k邻域查找。其次,针对采样点与其相邻点的曲率与法向量差和局部平均距离等定义了采样点的特征值,应用基于区域生长的方式进行点云区域的分割。最后,依据分割区域的标准对平坦区域和非平坦区域应用不同的点云精简方式,自动调整精简策略。该算法提高了平坦区域点云精简率,降低了非平坦区域的精简率,实现了对原始点云特征区域的高度保持,在提高精简率、精简效率和特征保留度方面与传统全局精简算法相比具有明显优势。
最后,由于受声呐测量原理和海底海貌特征的影响,经过点云精简后的数据会存在细小的孔洞,且孔洞的特点为孔洞边界凹凸不平分布不均匀。与传统三维点云不同的是,海底声呐数据除包含x、y和z三维坐标外,存在第四维回声强度坐标i。为此,本发明将海底声呐精简后点云作为待修补对象,针对海底声呐点云孔洞特点改进了基于局部扩张同心圆的点云孔洞修补算法,并以此为基础加入了第四维声强信息的修补,实现对海底声呐点云数据的最终处理。
本发明的有益效果为:本发明根据声呐点云的特性进行点云数据的融合与去燥,然后针对声呐点云的几何特性进行分区的点云精简,最后利用基于同心圆局部扩张的点云孔洞修补算法进行四维海底声呐点云孔洞的修补的海底点云数据处理方法;经本发明相关方法处理后的点云数据便于滤波算法对声呐点云数据中的人工目标进行识别与检测,既节省了人工对声呐图像数据判读所消耗的时间,又提高了人工目标识别算法的精确性,具有十分重大的意义。
附图说明
图1是本发明一种海底声呐点云数据处理方法的全局技术路线图
图2是本发明一种海底声呐点云数据处理方法实施例海底声呐点云孔洞修补前后效果对比图
具体实施方式
下面结合附图和具体实施方式对本发明进行详细说明。
一种海底声呐点云数据处理方法,如图1,的实现步骤为:
步骤1:进行原始点云数据的融合,步骤如下:
(1)处理测深数据,计算每个y值对应的x的个数n;
(2)处理侧扫数据,计算每个y值对应的x的个数m,每个y值对应的强度i;
(3)计算具有相同y值的测深数据并对应m/n个侧扫数据;
(4)遍历测深数据,对每个测深数据找到其对应的m/n个侧扫数据,对m/n个侧扫数据的强度求平均值并将其赋值给测深数据,这样每个测深数据形成了(x,y,z,i)。
步骤2:对融合后的点云数据进行高频滤波,高频噪声主要指漂浮在主点群上方的稀疏点和远离主点群中心的具有高频率和大幅值等特点的点群集合。针对高频噪声的特点,采用一种基于近邻点平均距离的滤波算法,将明显离群点进行剔除。去噪后得到的新的点云数据集S'为:
算法具体步骤如下:
(1)针对原始点云中的某一点Si,计算其k邻域内的所有点;
步骤3:对高频去燥后点云进行点云法向量修正,由于低频噪声的去噪双边滤波算法无法直接对法向散乱的无序点云直接进行滤波,所以需要针对该区域进行法向量修正。具体步骤如下:
将原始点云经高频滤波处理后进行读入:
S={Sj},j∈J,SJ=A3 (5)
存在S中的某一点Sj可以完成协方差矩阵C的构造,其中J代表原始点云。
其中;Ns为Sj的邻域;Sc则代表Ns内所有点质量的中心。
协方差分析法其自身虽具有抗噪声效果,但是当面对点云的杂乱噪声时,该方法需要进一步修改法向以满足需求。本发明采用高斯权函数进行法矢nj的平滑,求得平滑后的法矢ni'。
步骤4:进行声呐点云低频去燥,该算法表达式为:
s'i=si-ε·n (8)
其中:s'i为滤波后的点;si为初始采样点;ε为原始双边滤波因子;n为初始采样点si的法向量方向。
其中:
光顺滤波权函数Wc(x)定义为:
特征保持权函数Wp(y)定义为:
其中σc和σs为高斯滤波参数,通常情况下,σc等于采样点si的邻域半径,σs等于采样点si的邻域标准差。当σc为定值时,σs的大小变化则十分明显地影响着点云特征变化。当σs为定值时,σc的大小变化则十分明显地影响着点云的光顺程度。
算法具体步骤如下:
(1)针对某采样点si,计算其k个近邻点zij,j=1,2,...,k;针对每个k邻近点,计算Wc(x)函数的x=||si-zij||,其中x代表采样点si到其k近邻点zij的距离;计算函数Wp(y)的y=<nj,si-zij>,其中y代表采样点si与其k近邻点zij距离向量si-zij与采样点si法向量方向的内积;
(2)由(10)和(11)计算出Wc(x)和Ws(x);
(3)将Wc(x)和Ws(x)代入式(9),计算出ε;
(4)计算去噪处理后的新增点集;
(5)在所有目标点计算完毕后,结束算法,输出处理后的点云数据。
步骤5:针对去燥后的点云进行点云分区精简。针对海底声呐点云数据的采样点数量庞大且分布不均匀、特征点多的特性,本发明提出一种基于海底声呐点云特征的分区精简方法。首先,使用kd-tree进行点的k邻域查找。其次,针对采样点与其相邻点的曲率与法向量差和局部平均距离等定义了采样点的特征值,应用基于区域生长的方式进行点云区域的分割。最后,依据分割区域的标准对平坦区域和非平坦区域应用不同的点云精简方式,自动调整精简策略。该算法提高了平坦区域点云精简率,降低了非平坦区域的精简率,实现了对原始点云特征区域的高度保持,在提高精简率、精简效率和特征保留度方面与传统全局精简算法相比具有明显优势。具体实现步骤如下:
(1)计算法向量差
对局部平面S的法向量n进行标准化,使其||n||=1,因此得到的协方差矩阵M可表示为:
如果稠密且光滑点云中存在某两个点xi和xj,且它们的距离足够接近,则可近似地认为两点的法向量ni和nj平行,即ni·nj≈±1,如果数据点的法向量方向相同,则ni·nj≈1,否则,ni和nj中的一个需要反转过来。具体步骤如下:
(1)构建无向连通图G。对于点xi和xj距离越近,就越有可能法向量同向,所以,以1-|ni·nj|作为边(i,j)的权重。
(2)确定初始传播点。将极值点作为法向估计的初始传播点,该点的法向方向应与坐标方向一致。
(3)确定传播方向。根据无向图的最小生成树的传播方向,最小化传播路径的权重以确保全局优化。
(4)通过传播方向。将目标点云中所有点进行重定向,便可得到所有点集的法向量估计。
计算采样点xi的法向量ni,设其邻域点集为N(xi),则点xi与邻域点的法向差可表示为:
法向量差综合考虑了目标点与其邻域点之间的法向量偏差。该值越大,该点起伏越大,该点是特征点的可能性越大;该值越小,该点周围曲面较为平坦,该点处于点云平滑区域,是特征点的可能性较小。
(2)计算区域平均距离
区域平均距离是指点云中的点与其相邻区域中的所有点的欧式距离均值。通过上文采用的kd-tree法计算点云中的某一点x的邻域点集可表示为:ω(xi)(i=1,2,3,...,k)。对任意一点x的k领域计算如下:
区域平均距离反映了点云的密集程度,d(x)的值越大,该点成为特征点的几率越高。
(3)计算曲率
曲率反映目标点周围弯曲程度,由(14)求得的协方差矩阵M的特征值λj(j=0,1,2)均为非负实数,假设:(0<λ0<λ1<λ2),与之相对应的三个特征向量vj(j=0,1,2)组成一个正交基。定义目标点的曲率cui为:
(4)计算特征阈值及参数
其中α是特征系数,N(xi)表示点xi的邻域集合,cui表示该点的最终曲率,d(x)表示该点的局部平均距离。
精简的特征阈值可以由点云密度和点云尺寸进行定义,点云密度可以由八叉树划分后的叶子节点对角线长度l进行定义,而目标点与点云中心点的最大距离dmax可以对点云尺寸进行衡量,所以有:
γ=l/dmax (19)
经过对特征阈值和特征参数的定义与设定,就可以有效地制定针对点云非平坦区域进行精简的策略。
步骤6:点云区域分割
(1)选取初始点P1,将该点与其k邻域内欧式距离最小的点进行连接得到点P2,形成初始边P1P2,在P1和P2的k邻域中找到符合P2P3+P1P3最短的点P3,从而构造初始三角形P1P2P3。
(2)对于初始ΔP1P2P3,将三角形的边界边作为初始边,在曲面S上进行ΔABC的扩展,由于D点的夹角是ΔABC和ΔDBC大于规定范围,选取夹角范围内的点E作为网格点。在实验过程中,对C值选定一个合适的阈值范围对点云的区域生长效果尤为重要。计算相邻三角网格法向量的夹角α,并设定C值对点云表面变化进行表征。假设存在两点M(x1,x2,x3)和N(y1,y2,y3),根据两点之间距离公式有:
对M和N的k邻域点进行计算,并采用同样方式寻找构造点进行扩展:
其中和分别表示一条重叠边关联的两个平面法向量。由(21)分析可得,C值的取值范围是[0-2]且随法向量夹角α单调递增,如果C值范围较小,则说明新加入的三角形较为平缓,将该点加入初始点所属分片中,否则停止该边进行扩展。
(3)在扩展过程中,遇见已扩展的点集时,停止该边的扩展。如果所有边都已扩展完成,将原始数据点与已扩展的点进行分离,跳转至(1)直到所有点云数据均分离完毕。
应用区域生长法进行点云分割,可以基本保证平缓区域的点在同一分组内,所以该方法可以基本保证平坦区域和非平坦区域的分割。
步骤7:点云分区域精简
针对海底声呐点云这类大场景复杂点云的精简,不能单纯应用传统全局简化算法。全局精简算法虽然在精简率、精简时间和特征点保留方面具有优势,但点云场景复杂时无法达到兼顾。点云的平坦区域只需取部分点就可以做完整描述,而几何特征丰富的区域,点云特征点保留的越多,越能更好的描述区域细节。为了平衡海底声呐点云数据集上精简率和精简精度两大指标,本发明提出了基于声呐点云特征的分区精简方法。
针对平坦区域主要采用包围盒法进行精简,算法步骤如下:
(1)计算点面平均距离,将超过点面平均距离的点提取并保存。
(2)应用包围盒对该区域的点进行划分,每个包围盒里面的点用其重心点代替。重心点是距离中心最近的那个点,通过对中心点位置的确定,便可计算重心点。
通过计算得到距离中心点最近的重心点Sg,并将包围盒内的重心点进行保留。
(3)提取平坦区域边界点,并与步骤(1)和(2)的点进行点云的拼接形成精简后点云。
针对非平坦区域具体步骤为:
(1)计算点的k邻域,计算特征参数(法向差、局部平均距离和曲率),并对点进行特征判断,提取特征点进行保留。
(2)将阈值内的点进行曲率判断,对于曲率较小的点做较大程度精简。
(3)将两类点云进行拼接,形成精简后点云。
步骤8:受声呐测量原理和海底海貌特征的影响,经过点云精简后的数据会存在细小的孔洞,且孔洞的特点为孔洞边界凹凸不平分布不均匀。与传统三维点云不同的是,海底声呐数据除包含x、y和z三维坐标外,存在第四维回声强度坐标i。所以选择合适的点云修补算法并改进应用于第四维坐标修补尤为重要。
基于点对领域关系的方法是通过对数据点的邻域关系从而进行判断某点是否为边界点。先计算每个数据点的k近邻,然后根据求该点的k邻域点的最小二乘拟合曲面,将该点和其邻域点进行投影,再根据某种规则进行边界识别,识别后的边界点包含内外边界,而外边界属于点云边缘点而非孔洞边界点。
步骤9:由上一步骤获得了孔洞的边界特征点,使用质心做圆来解决海底声呐点云四维的孔洞修补。该方法可以有效修补海底声呐点云数据产生的三维孔洞,并可以针对第四维声强信息进行合理有效的修补。
算法主要思路是:
假设获取到的边界点的集合为T={Ti(xi,yi,zi,ii),i=1,2...,n}计算边界的质心如下:
设边界点到质心的最短距离为dmin则有:
其中dave是上文求得的点云全局平均距离。将l向下取整求得待扩张的同心圆数目,将所有孔洞点与质心连线,连线与同心圆的交点便是孔洞修补的填充点。设Ti与孔洞边界点的质心连线与同心圆的交点为Tadd,则有:
其中:N对应着向外扩的第N个圆,N=0,1,...,l。
关于第四维数据强度I,这里规定位于同一半径连接线上的所有点的I值相等,即从质心到最外边界点Ti的连线所有与同心圆的交点的强度值相等。质心强度值为所有孔洞边界强度平均值。
由质心法扩圆的具体实现步骤如下:
(1)计算孔洞边界数据点的质心;
(2)计算(24)中的l并向下取整作为待扩张的同心圆数目;
(3)进行同心圆扩充,扩充个数为全局平均距离dave的N倍;
(4)将边界点与所有边界点的质心点进行连接,连接线与步骤(3)中扩充的同心圆的交点作为新增的孔洞三维填充点;
(5)将同一质心与边界线连线上的所有点集的强度值应用连线最外边界强度值进行赋值。得到修补后的新四维点集数据集合A。
由于同心圆的特征可知,修补后的点云数据集具有一定的缺点,越靠近内圈的时候点云越密集,这样相对于原始数据会显得不真实。所以要对修补后的孔洞进行均匀化处理。具体实现步骤如下:
(1)读取孔洞边界及填补后的点云集Pn;
(2)计算点与点之间相邻数据之间线段的长度,如果所有线段之间的长度均大于全局平均点距离dave,进行(4),否则继续进行(3);
(3)对每个点的相邻线段长度和进行计算,剔除最短线段上的点,保存现有点云,并跳转至(1);
(4)保存所有的点集得到新的集合Pm。
一种海底声呐点云数据处理方法的实施例
下面通过某真实海底四维声呐点云作为实例来说明一种海底声呐点云数据处理方法的实现过程。
首先,进行原始点云数据的融合,。
其次,执行步骤2和步骤3得到去燥后的点云图像。然后,执行步骤5~步骤7,完成对该点云数据的精简过程,实现基于声呐点云特征的分区精简,
最后,执行步骤8,对精简后带有孔洞的区域采用四维局部扩张同心圆算法进行点云孔洞修补,其结果如图2所示。
Claims (7)
1.一种海底声呐点云数据处理方法,其特征在于,具体包括以下步骤:
步骤1:进行原始点云数据的融合;
步骤2:对步骤1融合后的点云数据进行高频滤波,高频噪声主要指漂浮在主点群上方的稀疏点和远离主点群中心的具有高频率和大幅值特点的点群集合,针对高频噪声的特点,采用一种基于近邻点平均距离的滤波算法,将明显离群点进行剔除,去噪后得到的新的点云数据集S'为:
步骤3:对步骤2中高频滤波后点云进行点云法向量修正;
步骤4:进行声呐点云低频去噪,该算法表达式为:
s’i=si-ε·n (5)
其中:s’i为滤波后的点;si为初始采样点;ε为原始双边滤波因子;n为初始采样点si的法向量方向;
光顺滤波权函数Wc(x)定义为:
特征保持权函数Wp(y)定义为:
其中σc和σs为高斯滤波参数,σc等于采样点si的邻域半径,σs等于采样点si的邻域标准差;当σc为定值时,σs的大小变化影响点云特征变化;当σs为定值时,σc的大小变化影响点云的光顺程度;
步骤5:针对步骤4中低频去噪后的点云进行全局精简;
首先,使用kd-tree进行点的k邻域查找;
其次,针对采样点与其相邻点的曲率与法向量差和局部平均距离定义了采样点的特征值,应用基于区域生长的方式进行点云区域的分割;
最后,依据分割区域的标准对平坦区域和非平坦区域应用不同的点云精简方式,自动调整精简策略;
步骤6:将步骤5中精简过后的点云区域进行分割
a.选取初始点P1,将该点与其k邻域内欧式距离最小的点进行连接得到点P2,形成初始边P1P2,在P1和P2的k邻域中找到符合P2P3+P1P3最短的点P3,从而构造初始三角形P1P2P3;
b.对于初始ΔP1P2P3,将三角形的边界边作为初始边,在曲面S上进行ΔABC的扩展,由于D点的夹角是ΔABC和ΔDBC大于规定范围,选取夹角范围内的点E作为网格点;在实验过程中,对C值选定一个合适的阈值范围对点云的区域生长效果尤为重要;计算相邻三角网格法向量的夹角α,并设定C值对点云表面变化进行表征;假设存在两点M(x1,x2,x3)和N(y1,y2,y3),对M和N的k邻域点进行计算,并采用同样方式寻找构造点进行扩展:
其中和分别表示一条重叠边关联的两个平面法向量;由(9)分析可得,C值的取值范围是[0-2]且随法向量夹角α单调递增,如果C值范围较小,则说明新加入的三角形较为平缓,将该点加入初始点所属分片中,否则停止该边进行扩展;
c.在扩展过程中,遇见已扩展的点集时,停止该边的扩展;如果所有边都已扩展完成,将原始数据点与已扩展的点进行分离,跳转至a直到所有点云数据均分离完毕;
步骤7:将点云分为平坦与非平坦区域进行局部精简,确保复杂场景下的精简精度
针对平坦区域采用包围盒法进行精简,取部分点就可以做完整描述
针对非平坦区域采用阈值分析法进行精简,采用点云特征点保留的越多,越能更好的描述区域细节;
步骤8:海底声呐数据除包含x、y和z三维坐标外,存在第四维回声强度坐标i,采用点云修补算法对第四维坐标进行修补;
通过对数据点的邻域关系从而进行判断某点是否为边界点,先计算每个数据点的k近邻,然后根据求该点的k邻域点的最小二乘拟合曲面,将该点和其邻域点进行投影,再进行边界识别,识别后的边界点包含内外边界,外边界属于点云边缘点,内边界属于孔洞边界点,k点指代某点的邻域点个数为k,是一种点云邻域关系的指代方式;
步骤9:由上一步骤获得了孔洞的边界特征点,使用质心做圆来解决海底声呐点云四维的孔洞修补,该方法可以有效修补海底声呐点云数据产生的三维孔洞,并可以针对第四维声强信息进行合理有效的修补,经本方法修补后的点云孔洞分布均匀,填充点与孔洞点平滑相接,孔洞填充点周围点云声呐分布合理,实现了对海底声呐点云数据的有效处理。
4.根据权利要求1所述的一种海底声呐点云数据处理方法,其特征在于,所述步骤4中低频去噪的具体步骤如下:
1)针对某采样点si,计算其k个近邻点zij,j=1,2,...,k;针对每个k邻近点,计算Wc(x)函数的x=||si-zij||,其中x代表采样点si到其k近邻点zij的距离;计算函数Wp(y)的y=<nj,si-zij>,其中y代表采样点si与其k近邻点zij距离向量si-zij与采样点si法向量方向的内积;
2)由(7)和(8)计算出Wc(x)和Ws(x);
3)将Wc(x)和Ws(x)代入式(6),计算出ε;
4)计算去噪处理后的新增点集;
5)在所有目标点计算完毕后,结束算法,输出处理后的点云数据。
5.根据权利要求1所述的一种海底声呐点云数据处理方法,其特征在于,所述步骤5中全局精简的具体实现步骤如下:
a.计算法向量差
法向量差主要指某一点xi与其邻域点N(xi)的法向量偏差,它反映了某点的弯曲程度;假设存在某局部曲面S,某点xi与坐标中心点的距离为di,则:
对局部平面S的法向量n进行标准化,使其||n||=1,因此得到的协方差矩阵M可表示为:
如果稠密且光滑点云中存在某两个点xi和xj,且它们的距离足够接近,则可近似地认为两点的法向量ni和nj平行,即ni·nj≈±1,如果数据点的法向量方向相同,则ni·nj≈1,否则,ni和nj中的一个需要反转过来;具体步骤如下:
1)构建无向连通图G;对于点xi和xj距离越近,就越有可能法向量同向,所以,以1-|ni·nj|作为边(i,j)的权重;
2)确定初始传播点;将极值点作为法向估计的初始传播点,该点的法向方向应与坐标方向一致;
3)确定传播方向;根据无向图的最小生成树的传播方向,最小化传播路径的权重以确保全局优化;
4)通过传播方向;将目标点云中所有点进行重定向,便可得到所有点集的法向量估计;
计算采样点xi的法向量ni,设其邻域点集为N(xi),则点xi与邻域点的法向差可表示为:
法向量差综合考虑了目标点与其邻域点之间的法向量偏差;该值越大,该点起伏越大,该点是特征点的可能性越大;该值越小,该点周围曲面较为平坦,该点处于点云平滑区域,是特征点的可能性较小;
b.计算区域平均距离
区域平均距离是指点云中的点与其相邻区域中的所有点的欧式距离均值;通过上文采用的kd-tree法计算点云中的某一点x的邻域点集可表示为:ω(xi)(i=1,2,3,...,k);对任意一点x的k领域计算如下:
区域平均距离反映了点云的密集程度,d(x)的值越大,该点成为特征点的几率越高;
c.计算曲率
曲率反映目标点周围弯曲程度,由(15)求得的协方差矩阵M的特征值λj(j=0,1,2)均为非负实数,假设:(0<λ0<λ1<λ2),与之相对应的三个特征向量vj(j=0,1,2)组成一个正交基;定义目标点的曲率cui为:
d.计算特征阈值及参数
其中α是特征系数,N(xi)表示点xi的邻域集合,cui表示该点的最终曲率,d(x)表示该点的局部平均距离;
精简的特征阈值可以由点云密度和点云尺寸进行定义,点云密度可以由八叉树划分后的叶子节点对角线长度l进行定义,而目标点与点云中心点的最大距离dmax可以对点云尺寸进行衡量,所以有:
γ=l/dmax (20)
经过对特征阈值和特征参数的定义与设定,就可以有效地制定针对点云非平坦区域进行精简的策略。
6.根据权利要求1所述的一种海底声呐点云数据处理方法,其特征在于,所述步骤7中
a.针对平坦区域主要采用包围盒法进行精简,步骤如下:
1)计算点面平均距离,将超过点面平均距离的点提取并保存;
2)应用包围盒对该区域的点进行划分,每个包围盒里面的点用其重心点代替;重心点是距离中心最近的那个点,通过对中心点位置的确定,便可计算重心点;
通过计算得到距离中心点最近的重心点Sg,并将包围盒内的重心点进行保留;
3)提取平坦区域边界点,并与步骤(1)和(2)的点进行点云的拼接形成精简后点云;
b.针对非平坦区域,步骤如下:
1)计算点的k邻域,计算特征参数,即法向差、局部平均距离和曲率,并对点进行特征判断,提取特征点进行保留;
2)将阈值内的点进行曲率判断,对于曲率较小的点做较大程度精简;
3)将两类点云进行拼接,形成精简后点云。
7.根据权利要求1所述的一种海底声呐点云数据处理方法,其特征在于,所述步骤9中四维的孔洞修补具体步骤为:
假设获取到的边界点的集合为T={Ti(xi,yi,zi,ii),i=1,2...,n}计算边界的质心如下:
设边界点到质心的最短距离为dmin则有:
其中dave是上文求得的点云全局平均距离;将l向下取整求得待扩张的同心圆数目,将所有孔洞点与质心连线,连线与同心圆的交点便是孔洞修补的填充点;设Ti与孔洞边界点的质心连线与同心圆的交点为Tadd,则有:
其中:N对应着向外扩的第N个圆,N=0,1,...,l;
关于第四维数据强度I,这里规定位于同一半径连接线上的所有点的I值相等,即从质心到最外边界点Ti的连线所有与同心圆的交点的强度值相等;质心强度值为所有孔洞边界强度平均值;
由质心法扩圆的具体实现步骤如下:
1)计算孔洞边界数据点的质心;
2)计算(23)中的l并向下取整作为待扩张的同心圆数目;
3)进行同心圆扩充,扩充个数为全局平均距离dave的N倍;
4)将边界点与所有边界点的质心点进行连接,连接线与步骤(3)中扩充的同心圆的交点作为新增的孔洞三维填充点;
5)将同一质心与边界线连线上的所有点集的强度值应用连线最外边界强度值进行赋值;得到修补后的新四维点集数据集合A;
由于同心圆的特征可知,修补后的点云数据集具有一定的缺点,越靠近内圈的时候点云越密集,这样相对于原始数据会显得不真实;所以要对修补后的孔洞进行均匀化处理;具体实现步骤如下:
1)读取孔洞边界及填补后的点云集Pn;
2)计算点与点之间相邻数据之间线段的长度,如果所有线段之间的长度均大于全局平均点距离dave,进行(4),否则继续进行(3);
3)对每个点的相邻线段长度和进行计算,剔除最短线段上的点,保存现有点云,并跳转至(1);
4)保存所有的点集得到新的集合Pm。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910576810.4A CN110517193B (zh) | 2019-06-28 | 2019-06-28 | 一种海底声呐点云数据处理方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910576810.4A CN110517193B (zh) | 2019-06-28 | 2019-06-28 | 一种海底声呐点云数据处理方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN110517193A CN110517193A (zh) | 2019-11-29 |
CN110517193B true CN110517193B (zh) | 2022-04-12 |
Family
ID=68623554
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910576810.4A Active CN110517193B (zh) | 2019-06-28 | 2019-06-28 | 一种海底声呐点云数据处理方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN110517193B (zh) |
Families Citing this family (13)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111722201A (zh) * | 2020-06-23 | 2020-09-29 | 常州市贝叶斯智能科技有限公司 | 一种用于室内机器人激光雷达数据去噪方法 |
CN111861941B (zh) * | 2020-07-31 | 2024-03-22 | 富德康(北京)科技股份有限公司 | 一种三维空间测量结果数据的补偿算法 |
CN112147619B (zh) * | 2020-09-24 | 2021-12-10 | 中交第三航务工程局有限公司江苏分公司 | 基于声呐点云数据的桩体间距离迭代确定方法 |
CN112529891B (zh) * | 2020-12-21 | 2024-03-19 | 深圳辰视智能科技有限公司 | 一种基于点云的空心孔洞识别与轮廓检测方法、装置及存储介质 |
CN112950622B (zh) * | 2021-03-29 | 2023-04-18 | 上海商汤临港智能科技有限公司 | 一种目标检测方法、装置、计算机设备和存储介质 |
CN113421290A (zh) * | 2021-07-05 | 2021-09-21 | 山西大学 | 一种基于无人机图像采集的电厂锅炉内三维重建方法 |
CN113640829A (zh) * | 2021-08-18 | 2021-11-12 | 上海同陆云交通科技有限公司 | 一种基于LiDAR的无人机桥底检测*** |
CN114677454B (zh) * | 2022-03-25 | 2022-10-04 | 杭州睿影科技有限公司 | 一种图像生成方法和装置 |
CN115239625B (zh) * | 2022-06-21 | 2023-05-09 | 厦门微图软件科技有限公司 | 顶盖焊点云缺陷检测方法、装置、设备及存储介质 |
CN115578524B (zh) * | 2022-12-09 | 2023-03-31 | 华东交通大学 | 红外三维重构方法、***、存储介质及计算机设备 |
CN116152454A (zh) * | 2023-02-15 | 2023-05-23 | 中铁水利信息科技有限公司 | 基于gis和三维建模的水利实时监测管理*** |
CN116071274A (zh) * | 2023-03-23 | 2023-05-05 | 成都工业学院 | 一种基于点云区域生长的跨面孔洞修补方法及*** |
CN117235506B (zh) * | 2023-11-10 | 2024-02-13 | 四川大学 | 一种基于相空间重构的信号提取方法及装置 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103197308A (zh) * | 2013-03-15 | 2013-07-10 | 浙江大学 | 基于多波束相控阵声纳***的三维声纳可视化处理方法 |
CN109241978A (zh) * | 2018-08-23 | 2019-01-18 | 中科光绘(上海)科技有限公司 | 地基三维激光点云中平面片的快速提取方法 |
CN109801223A (zh) * | 2017-11-17 | 2019-05-24 | 宁波盈芯信息科技有限公司 | 一种点云去噪及空洞修补方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US8670591B2 (en) * | 2012-03-01 | 2014-03-11 | Exelis, Inc. | Foliage penetration based on 4D LIDAR datasets |
-
2019
- 2019-06-28 CN CN201910576810.4A patent/CN110517193B/zh active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103197308A (zh) * | 2013-03-15 | 2013-07-10 | 浙江大学 | 基于多波束相控阵声纳***的三维声纳可视化处理方法 |
CN109801223A (zh) * | 2017-11-17 | 2019-05-24 | 宁波盈芯信息科技有限公司 | 一种点云去噪及空洞修补方法 |
CN109241978A (zh) * | 2018-08-23 | 2019-01-18 | 中科光绘(上海)科技有限公司 | 地基三维激光点云中平面片的快速提取方法 |
Non-Patent Citations (1)
Title |
---|
Accurate Image Super-Resolution Using Very Deep Convolutional Networks;Jiwon Kim 等;《2016 IEEE Conference on Computer Vision and Pattern Recognition (CVPR)》;20160630;第1063-6919页 * |
Also Published As
Publication number | Publication date |
---|---|
CN110517193A (zh) | 2019-11-29 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN110517193B (zh) | 一种海底声呐点云数据处理方法 | |
CN112581457B (zh) | 一种基于三维点云的管道内表面检测方法及装置 | |
CN107123164B (zh) | 保持锐利特征的三维重建方法及*** | |
CN112070769B (zh) | 一种基于dbscan的分层点云分割方法 | |
CN106600622B (zh) | 一种基于超体素的点云数据分割方法 | |
CN106780524B (zh) | 一种三维点云道路边界自动提取方法 | |
CN110807781B (zh) | 一种保留细节与边界特征的点云精简方法 | |
CN111986115A (zh) | 激光点云噪声和冗余数据的精准剔除方法 | |
CN107358636B (zh) | 一种基于纹理合成的疏松缺陷图像生成方法 | |
CN112347550B (zh) | 耦合式室内三维语义建图及建模方法 | |
CN111612728B (zh) | 一种基于双目rgb图像的3d点云稠密化方法和装置 | |
CN107845073B (zh) | 一种基于深度图的局部自适应三维点云去噪方法 | |
CN107680168B (zh) | 三维重建中基于平面拟合的网格简化方法 | |
CN113920134B (zh) | 一种基于多线激光雷达的斜坡地面点云分割方法及*** | |
CN111524127B (zh) | 面向低空机载激光雷达数据的城市道路面提取方法 | |
CN113640826B (zh) | 一种基于3d激光点云的障碍物识别方法及*** | |
KR101549155B1 (ko) | 라이다 자료를 활용한 구조물의 직선경계 추출방법 | |
CN111968224A (zh) | 船舶3d扫描点云数据处理方法 | |
CN111783722B (zh) | 一种激光点云的车道线提取方法和电子设备 | |
CN108765445B (zh) | 一种肺气管分割方法及装置 | |
CN114332291A (zh) | 一种倾斜摄影模型建筑物外轮廓规则提取方法 | |
CN111986308A (zh) | 基于点云法线与曲率变化双约束曲面误差显著点识别方法 | |
CN107610174B (zh) | 一种鲁棒的基于深度信息的平面检测方法及*** | |
CN109359653A (zh) | 一种棉花叶部粘连病斑图像分割方法和*** | |
CN114898118A (zh) | 基于多源点云的输电线路房屋拆迁量自动统计方法及*** |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |