CN106056614A - 一种地面激光点云数据的建筑物分割与轮廓线提取方法 - Google Patents
一种地面激光点云数据的建筑物分割与轮廓线提取方法 Download PDFInfo
- Publication number
- CN106056614A CN106056614A CN201610393443.0A CN201610393443A CN106056614A CN 106056614 A CN106056614 A CN 106056614A CN 201610393443 A CN201610393443 A CN 201610393443A CN 106056614 A CN106056614 A CN 106056614A
- Authority
- CN
- China
- Prior art keywords
- point
- point cloud
- building
- neighborhood
- segmentation
- 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.)
- Pending
Links
- 230000011218 segmentation Effects 0.000 title claims abstract description 52
- 238000000605 extraction Methods 0.000 title claims abstract description 16
- 238000004422 calculation algorithm Methods 0.000 claims abstract description 33
- 238000000513 principal component analysis Methods 0.000 claims abstract description 21
- 238000004364 calculation method Methods 0.000 claims abstract description 6
- 238000000034 method Methods 0.000 claims description 58
- 239000011159 matrix material Substances 0.000 claims description 20
- VMXUWOKSQNHOCA-UKTHLTGXSA-N ranitidine Chemical compound [O-][N+](=O)\C=C(/NC)NCCSCC1=CC=C(CN(C)C)O1 VMXUWOKSQNHOCA-UKTHLTGXSA-N 0.000 claims description 8
- 238000009826 distribution Methods 0.000 claims description 7
- 239000000284 extract Substances 0.000 claims description 6
- 238000002372 labelling Methods 0.000 claims description 6
- 238000012217 deletion Methods 0.000 claims description 5
- 230000037430 deletion Effects 0.000 claims description 5
- 238000003491 array Methods 0.000 claims description 3
- 238000006243 chemical reaction Methods 0.000 claims description 3
- 238000003709 image segmentation Methods 0.000 claims description 3
- 230000000717 retained effect Effects 0.000 claims description 3
- 238000005070 sampling Methods 0.000 claims description 3
- 230000005484 gravity Effects 0.000 claims 1
- 239000006185 dispersion Substances 0.000 description 3
- 230000007812 deficiency Effects 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 238000002474 experimental method Methods 0.000 description 2
- 238000000638 solvent extraction Methods 0.000 description 2
- 241001269238 Data Species 0.000 description 1
- 230000003044 adaptive effect Effects 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 238000005192 partition Methods 0.000 description 1
- 238000004088 simulation Methods 0.000 description 1
- 238000005303 weighing Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T3/00—Geometric image transformations in the plane of the image
- G06T3/06—Topological mapping of higher dimensional structures onto lower dimensional surfaces
- G06T3/067—Reshaping or unfolding 3D tree structures onto 2D planes
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F18/00—Pattern recognition
- G06F18/20—Analysing
- G06F18/24—Classification techniques
- G06F18/241—Classification techniques relating to the classification model, e.g. parametric or non-parametric approaches
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T3/00—Geometric image transformations in the plane of the image
- G06T3/06—Topological mapping of higher dimensional structures onto lower dimensional surfaces
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06V—IMAGE OR VIDEO RECOGNITION OR UNDERSTANDING
- G06V10/00—Arrangements for image or video recognition or understanding
- G06V10/40—Extraction of image or video features
-
- 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/10004—Still image; Photographic image
- G06T2207/10012—Stereo images
Landscapes
- Engineering & Computer Science (AREA)
- Theoretical Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Data Mining & Analysis (AREA)
- Bioinformatics & Cheminformatics (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Evolutionary Biology (AREA)
- Evolutionary Computation (AREA)
- Bioinformatics & Computational Biology (AREA)
- General Engineering & Computer Science (AREA)
- Artificial Intelligence (AREA)
- Life Sciences & Earth Sciences (AREA)
- Multimedia (AREA)
- Image Analysis (AREA)
- Image Processing (AREA)
Abstract
本发明公开了一种地面激光点云数据的建筑物分割与轮廓线提取方法,包括对点云进行竖直投影、生成二维灰度图像、利用Otsu算法分割建筑物点云、主成分分析、最佳邻域计算、轮廓线提取等6个步骤;本发明实现了从原始点云中快速、准确地分割出建筑物点云,并且全自动地完成建筑物轮廓线点云的提取。
Description
技术领域
本发明属于地面激光点云数据处理应用技术领域,涉及一种对原始三维点云数据进行建筑物的精确分割和自动的轮廓线提取方法。
背景技术
地面激光扫描技术是一项新兴的三维测量扫描技术,其获取的海量高精度的点云数据,可以提供大规模精细的三维信息,更好地反映建筑物的细节和层次信息,对构建高精度数字模型而言,是一种崭新的方法。地面激光为我们提供高质量点云数据的同时,数据量也成倍的增长,对我们处理它产生了一定的难度。造成这一问题的原因主要有两个:一是在扫描目标建筑物点云的同时,不可避免的引入噪声和非目标地物点云数据;二是大多数建筑物立面可看成由许多大小不等的平面组成,每个平面除去少量外轮廓的点云对后续建模有帮助,大部分平面内部点云是可以被删除的。本发明旨在解决上述两个问题,并将发明分为建筑物点云分割和轮廓线提取两个部分。
目前,建筑物点云分割的方法有区域增长法、基于边缘的方法、基于特征聚类的方法、基于二维影像的方法、基于回波强度的方法等。上述分割算法中,有的精度高但计算复杂,有的耗时短但存在误分割或适应性的问题。目前,还没有一种完善的算法,在分割的计算中可以很好地兼顾精度和速度。本发明结合基于二维影像的方法和区域增长法,并做出一定的改进,完成了一种兼顾精度和速度的分割算法。
目前,建筑物点云轮廓线提取的方法主要有基于三角格网的方法、基于格网划分的方法、基于点的特征信息的方法等。这些方法多是根据点云的细节信息,利用特征提取或者条件约束完成轮廓线点云的提取,它们有着共同的问题就是计算较为复杂,时间复杂度高,并且针对不同的点云数据实验效果好坏不一。本发明重点研究了点云特征描述方法和最佳尺度确定的方法,充分利用点云的邻域信息和空间几何关系,并引入主成分分析与香农熵的计算,不但提高了轮廓线提取算法的精度和自适应性,而且针对不同的地面激光点云数据,无需反复调整阈值,实现了轮廓线的自动提取。
发明内容
为了解决上述技术问题,本发明提供了一种地面激光点云数据的建筑物分割与轮廓线提取方法,本发明主要涉及两部分内容,即建筑物点云分割模块与建筑物点云轮廓线提取模块。本发明的目的在于能够实现从原始三维点云中快速、准确地分割出建筑物点云、并且全自动地完成建筑物轮廓线点云的提取。
本发明所采用的技术方案是:一种地面激光点云数据的建筑物分割与轮廓线提取方法,其特征在于,包括以下步骤:
步骤1:对点云进行竖直投影;
输入原始三维点云,对原始三维点云进行竖直投影,凸显不同地物的点云分布特点;
步骤2:生成二维灰度图像;
根据点云的X、Y坐标范围建立水平格网,统计每个格网内投影点个数,并以此生成二维灰度图像,建立三维点云与二维影像之间的联系;
步骤3:利用Otsu算法分割建筑物点云;
对原始三维点云进行建筑物分割;根据建筑物在二维灰度图像中的灰度特点,使用Otsu算法对图像进行分割,图像分割后再将影像转换为三维点云,即可完成建筑物分割;
步骤4:主成分分析;
对分割后的建筑物点云进行主成分分析,计算出每个点的维度特征值;
步骤5:最佳邻域计算;
在不同尺度下分别进行主成分分析,计算不同尺度下维度特征值,再根据熵函数,自动计算出每个点的最佳邻域,得出每个点的主维度,从而将点云分为线状点云、面状点云、散乱点云三类;
步骤6:轮廓线提取;
根据几何位置关系,对面状点云进行轮廓线提取处理;先为每个点确定其所在的平面,再判断此点是在平面的边缘还是平面内部,根据平面中的点与其邻域点形成的线段夹角和约为360°来滤除平面中的点,从而保留每个平面的轮廓线点云。
作为优选,步骤1中所述的对原始三维点云进行竖直投影,是将原始三维点云沿Z轴投影至X、Y平面上,仅保留每个点的X、Y坐标;其投影公式为:
作为优选,步骤2中所述的统计每个格网内投影点个数,记为投影点密度DoPP,调整所有格网的DoPP的取值范围,将其线性拉伸到0-255,用拉伸后的值模拟每个格网的灰度,生成二维灰度图像。
作为优选,步骤3中所述的Otsu算法对图像进行分割,是计算出灰度阈值T,以阈值T将影像分为目标和背景两类,使它们的类间方差g为最大,g由公式(2)计算得出;
g=ω0ω1(μ0-μ1)2 (2);
其中ω0、ω1为两类各占比例,μ0、μ1为两类的平均灰度。
作为优选,所述的Otsu算法的包括以下两点改进:
(1)分块求解;在进行Otsu算法的时候,并不是对所有的点云数据仅进行一次计算,而是依据建筑物立面在整个场景内的分布情况,对场景范围进行划分,在每个子场景内分别进行一次Otsu算法;划分原则为,尽可能多的满足,在每个子场景中,同时含有建筑物立面数据和其他地物数据,尤其是保证建筑物点云中DoPP值相对较小的部分满足上述要求,这样做可以显著提升点云分割的正确率;
将图像L在X和Y方向上分别划分为m和n等份,得到每个子场景的范围:
Lij=xi*yi(i=1,2,…,m;j=1,2,…,n) (3);
在每个Lij内进行一次Otsu算法,最后取每个子场景分割出的建筑物立面数据的合集,得到最终的分割结果;
(2)结合建筑物立面点云在竖直投影后可形成连续的线段这一特点,对于漏分割的点云部分,利用区域增长算法进行补偿;其具体的计算流程包括以下子步骤:
步骤3.1:选取种子点;设原始三维点云投影生成灰度图像的像素灰度为g[i][j],经过Otsu分块求解算法处理后像素灰度为G[i][j],设定标记数组B[i][j],记录每个像素的属性,三个数组的相同的下标对应同一个像素;遍历每个像素,若有像素p满足G[px][py]>0且其四邻域或八邻域有任意一个像素q符合G[qx][qy]=0且g[qx][qy]>0,那么B[px][py]=1,否则B[px][py]=0,其中像素p、q的关系满足公式(4);选取所有标记为1的像素作为种子点;
步骤3.2:设定区域增长条件,即阈值T;设种子点像素为p、其邻域点像素为q,将种子点G[px][py]与其邻域点g[qx][qy]作比较,若其商满足G[px][py]/g[qx][qy]<T,则将g[qx][qy]值赋予G[qx][qy],并将B[qx][qy]标记为1,若其商满足G[px][py]/g[qx][qy]≥T,则将B[qx][qy]标记为0,当p与其对应的所有q均完成上述操作时,将B[px][py]标记为0;
步骤3.3:区域增长;对B=1的像素看做新的种子点,并返回执行步骤3.2;
步骤3.4:结束条件;数组B中的所有元素均为0,即所有点的标记均为0,区域增长结束。
作为优选,步骤4中所述的主成分分析是通过对每个点及其邻域点的三维坐标构建协方差矩阵,计算出矩阵的三个特征值,以此分析出每个点属于三个维度特征的概率a1D、a2D、a3D;其具体实现过程包括以下子步骤:
步骤4.1:根据点云密度和精度要求,设定邻域半径r的取值范围[rmin,rmax],以及取值间隔;
步骤4.2:在不同邻域半径下,分别进行主成分分析;
首先利用每个扫描点Xi及其邻域点νr的三维坐标{Xi=(xi,yi,zi)|i∈vr}构造协方差矩阵:
矩阵C为一个3*3的矩阵,其中为点集νr的重心坐标,矩阵M的具体形式为:
然后计算矩阵C的三个特征值λ1、λ2、λ3,并按照λ1≥λ2≥λ3的规则进行排列;
步骤4.3:计算每个点的三个维度特征;
令根据公式(7)计算每个点的三个维度特征,即一维线状特征a1D、二维面状特征a2D和三维散乱点特征a3D;其中a1D、a2D、a3D的和为1,换言之,a1D、a2D、a3D分别代表了扫描点属于三个维度特征的概率;
作为优选,步骤5中所述的最佳邻域计算,是根据熵函数计算每个点的最佳邻域,见公式(8):
Ef(vr)=-a1Dln(a1D)-a2Dln(a2D)-a3Dln(a3D) (8);
其中a1D、a2D、a3D是主成分分析中得到的每个点属于三个维度特征的概率。
根据公式(8)计算不同r值下的邻域点集的熵值;当熵取最小值时,表示在此邻域下该点的主维度特征最突出,此时对应的半径r*即为最佳邻域半径。
作为优选,步骤6中所述的根据几何位置关系,对面状点云进行轮廓线提取处理;其具体实现过程包括以下子步骤:
步骤6.1:根据最佳邻域r*确定X点的最佳邻域点集
步骤6.2:确定点集所在平面。根据RANCAC随机抽样一致算法,从中随机抽取两个不同点并重复多次(重复次数与的点个数相同),计算X、Xi、Xj三点所在平面方程并统计中其余点到平面的欧氏距离和,选取使距离和最小的平面Sxpq的两邻域点Xp、Xq,将Xp、Xq从中删除;
步骤6.3:删除平面外部点。遍历计算每个Xi与Xp、Xq三点确定的平面Sipq方程,设定角度阈值θ1,当平面Sipq与平面Sxpq的夹角小于θ1时,认定Xi在平面Sxpq上,否则将Xi从中删除。
步骤6.4:设XXp为起始边,遍历计算每个∠XPXXi,设定角度阈值θ2,选取小于θ2中使∠XPXXi最小的邻域点X1';
步骤6.5:判断X1'与Xp的方向,如果X1'在Xp的顺时针方向,就按照顺时针方向去搜索;否则,按照逆时针搜索下一个邻域点;
步骤6.6:将步骤6.4得到的邻域点与X构成的线段设为起始边,重复步骤6.4寻求满足条件的下一个邻域点X2',X3',…,Xn',直到搜索不到满足条件的邻域点;
步骤6.7:计算的角度和,若和约等于360°,将点X视为平面中的点,将其滤除,否则视为边缘点,将其保留。
由于传统的建筑物点云分割方法和轮廓线提取方法在精度、速度、适用性、自动化等方面有着不同程度的不足,并且无法做到同时兼顾,因此本发明针对以上问题进行了改进,本发明的优点在于:
1、本发明提出的建筑物点云分割方法兼顾速度和精度。
(1)本发明在建筑物点云分割中,把三维点云转换成二维影像,使用整体分割的方法,避免了针对每一个点的计算,从而提高了计算速度。
(2)针对Otsu在建筑物点云分割中精度的不足,提出了分块求解和区域增长的改进方法,既提高了建筑物点云的分割精度,又不增加过多的计算,保证了分割速度快的优点。
2、本发明提出的建筑物点云轮廓线提取方法是一种全自动的方法,无需人工干预,并且适用性强。本发明通过主成分分析和熵函数可自动计算出每个点的最佳邻域半径,并将点云按照维度特征进行分类,最后根据面状点云内部点固有的几何关系自动提取轮廓线点云,在算法设计中特地将距离阈值转换为角度阈值,用每个点与邻域点的相对位置关系替代绝对位置关系,这样针对不同点云数据无需修改实验中的阈值,达到自动提取的目的。
附图说明
图1:本发明实施例的流程图。
具体实施方式
为了便于本领域普通技术人员理解和实施本发明,下面结合附图及实施例对本发明作进一步的详细描述,应当理解,此处所描述的实施示例仅用于说明和解释本发明,并不用于限定本发明。
请见图1,本发明提供的一种地面激光点云数据的建筑物分割与轮廓线提取方法,包括以下步骤:
步骤1:对点云进行竖直投影;
输入原始三维点云,对原始三维点云进行竖直投影,凸显不同地物的点云分布特点;
对原始三维点云进行竖直投影,是将原始三维点云沿Z轴投影至X、Y平面上,仅保留每个点的X、Y坐标;其投影公式为:
步骤2:生成二维灰度图像;
根据点云的X、Y坐标范围建立水平格网,统计每个格网内投影点个数,记为投影点密度DoPP,调整所有格网的DoPP的取值范围,将其线性拉伸到0-255,用拉伸后的值模拟每个格网的灰度,并以此生成二维灰度图像,建立三维点云与二维影像之间的联系;
步骤3:对点云进行初始的建筑物分割;
根据建筑物部分在二维灰度图像中的灰度特点,使用Otsu算法对图像进行分割,Otsu算法的主要原理是计算出灰度阈值T,以阈值T将影像分为目标和背景两类,使它们的类间方差g为最大,g由公式(2)计算得出。其中ω0、ω1为两类各占比例,μ0、μ1为两类的平均灰度。
g=ω0ω1(μ0-μ1)2 (2);
图像分割后再将影像转换为三维点云,即可完成建筑物初始分割。一般情况下,初始分割可得到建筑物主体部分,但存在部分错误分割区域。
步骤4:对点云进行精确的建筑物分割;
结合初始分割结果中,错误分割区域的分布特点,使用一种改进的Otsu算法进行精确分割。改进算法包含两个部分,分块求解和区域增长补偿漏分割区域。经过此步骤,可以准确、完成地分割出建筑物点云,剔除大量地面、树木等非目标地物点云数据。
(1)分块求解;在进行Otsu算法的时候,并不是对所有的点云数据仅进行一次计算,而是依据建筑物立面在整个场景内的分布情况,对场景范围进行划分,在每个子场景内分别进行一次Otsu算法;划分原则为,尽可能多的满足,在每个子场景中,同时含有建筑物立面数据和其他地物数据,尤其是保证建筑物点云中DoPP值相对较小的部分满足上述要求,这样做可以显著提升点云分割的正确率;
将图像L在X和Y方向上分别划分为m和n等份,得到每个子场景的范围:
Lij=xi*yi(i=1,2,…,m;j=1,2,…,n) (3);
在每个Lij内进行一次Otsu算法,最后取每个子场景分割出的建筑物立面数据的合集,得到最终的分割结果;
(2)结合建筑物立面点云在竖直投影后可形成连续的线段这一特点,对于漏分割的点云部分,利用区域增长算法进行补偿;其具体的计算流程包括以下子步骤:
步骤4.1:选取种子点;设原始三维点云投影生成灰度图像的像素灰度为g[i][j],经过Otsu分块求解算法处理后像素灰度为G[i][j],设定标记数组B[i][j],记录每个像素的属性,三个数组的相同的下标对应同一个像素;遍历每个像素,若有像素p满足G[px][py]>0且其四邻域或八邻域有任意一个像素q符合G[qx][qy]=0且g[qx][qy]>0,那么B[px][py]=1,否则B[px][py]=0,其中像素p、q的关系满足公式(4);选取所有标记为1的像素作为种子点;
步骤4.2:设定区域增长条件,即阈值T;设种子点像素为p、其邻域点像素为q,将种子点G[px][py]与其邻域点g[qx][qy]作比较,若其商满足G[px][py]/g[qx][qy]<T,则将g[qx][qy]值赋予G[qx][qy],并将B[qx][qy]标记为1,若其商满足G[px][py]/g[qx][qy]≥T,则将B[qx][qy]标记为0,当p与其对应的所有q均完成上述操作时,将B[px][py]标记为0;
步骤4.3:区域增长;对B=1的像素看做新的种子点,并返回执行步骤3.2;
步骤4.4:结束条件;数组B中的所有元素均为0,即所有点的标记均为0,区域增长结束。
步骤5:主成分分析;
轮廓线提取的首要步骤是按照点与其邻域点的位置关系,将点云分为线状点云、面状点云、散乱点云三类。主成分分析可以对复杂无序的点云数据进行分析计算,进而得到三个维度X、Y、Z之间的关系,确定每个点的类型。主成分分析的主要计算流程是通过对每个点及其邻域点的三维坐标构建协方差矩阵,计算出矩阵的三个特征值,以此分析出每个点属于三个维度特征的概率a1D、a2D、a3D;其具体实现过程包括以下子步骤:
步骤5.1:根据点云密度和精度要求,设定邻域半径r的取值范围[rmin,rmax],以及取值间隔;
步骤5.2:在不同邻域半径下,分别进行主成分分析;
首先利用每个扫描点Xi及其邻域点νr的三维坐标{Xi=(xi,yi,zi)|i∈vr}构造协方差矩阵:
矩阵C为一个3*3的矩阵,其中为点集νr的重心坐标,矩阵M的具体形式为:
然后计算矩阵C的三个特征值λ1、λ2、λ3,并按照λ1≥λ2≥λ3的规则进行排列;
步骤5.3:计算每个点的三个维度特征;
令根据公式(7)计算每个点的三个维度特征,即一维线状特征a1D、二维面状特征a2D和三维散乱点特征a3D;其中a1D、a2D、a3D的和为1,换言之,a1D、a2D、a3D分别代表了扫描点属于三个维度特征的概率;
步骤6:最佳邻域计算;
在不同尺度下分别进行主成分分析,计算不同尺度下维度特征值,再根据熵函数,自动计算出每个点的最佳邻域,得出每个点的主维度,从而将点云分为线状点云、面状点云、散乱点云三类;
在主成分分析计算每个点的维度特征时,点的维度特征会随着邻域点集的变化而改变,针对不同区域的点云数据,由于每个点的深度的不同,点云密度和邻域点数量的不确定性,我们很难在某一固定的尺度下得到最正确的特征点集。为此,根据熵函数计算每个点的最佳邻域,见公式(8);
Ef(vr)=-a1Dln(a1D)-a2Dln(a2D)-a3Dln(a3D) (8);
根据公式(8)计算不同r值下的邻域点集的熵值;当熵取最小值时,表示在此邻域下该点的主维度特征最突出,此时对应的半径r*即为最佳邻域半径。
步骤7:轮廓线提取;
根据几何位置关系,对面状点云进行轮廓线提取处理;先为每个点确定其所在的平面,再判断此点是在平面的边缘还是平面内部,根据平面中的点与其邻域点形成的线段夹角和约为360°来滤除平面中的点,从而保留每个平面的轮廓线点云。
根据几何位置关系,对面状点云进行轮廓线提取处理;其具体实现过程包括以下子步骤:
步骤7.1:根据最佳邻域r*确定X点的最佳邻域点集
步骤7.2:确定点集所在平面。根据RANCAC随机抽样一致算法,从中随机抽取两个不同点并重复多次(重复次数与的点个数相同),计算X、Xi、Xj三点所在平面方程并统计中其余点到平面的欧氏距离和,选取使距离和最小的平面Sxpq的两邻域点Xp、Xq,将Xp、Xq从中删除;
步骤7.3:删除平面外部点。遍历计算每个Xi与Xp、Xq三点确定的平面Sipq方程,设定角度阈值θ1,当平面Sipq与平面Sxpq的夹角小于θ1时,认定Xi在平面Sxpq上,否则将Xi从中删除。
步骤7.4:设XXp为起始边,遍历计算每个∠XPXXi,设定角度阈值θ2,选取小于θ2中使∠XPXXi最小的邻域点X1';
步骤7.5:判断X1'与Xp的方向,如果X1'在Xp的顺时针方向,就按照顺时针方向去搜索;否则,按照逆时针搜索下一个邻域点;
步骤7.6:将步骤7.4得到的邻域点与X构成的线段设为起始边,重复步骤7.4寻求满足条件的下一个邻域点X2',X3',…,Xn',直到搜索不到满足条件的邻域点;
步骤7.7:计算的角度和,若和约等于360°,将点X视为平面中的点,将其滤除,否则视为边缘点,将其保留。
应当理解的是,本说明书未详细阐述的部分均属于现有技术。
应当理解的是,上述针对较佳实施例的描述较为详细,并不能因此而认为是对本发明专利保护范围的限制,本领域的普通技术人员在本发明的启示下,在不脱离本发明权利要求所保护的范围情况下,还可以做出替换或变形,均落入本发明的保护范围之内,本发明的请求保护范围应以所附权利要求为准。
Claims (8)
1.一种地面激光点云数据的建筑物分割与轮廓线提取方法,其特征在于,包括以下步骤:
步骤1:对点云进行竖直投影;
输入原始三维点云,对原始三维点云进行竖直投影,凸显不同地物的点云分布特点;
步骤2:生成二维灰度图像;
根据点云的X、Y坐标范围建立水平格网,统计每个格网内投影点个数,并以此生成二维灰度图像,建立三维点云与二维影像之间的联系;
步骤3:利用Otsu算法分割建筑物点云;
对原始三维点云进行建筑物分割;根据建筑物在二维灰度图像中的灰度特点,使用Otsu算法对图像进行分割,图像分割后再将影像转换为三维点云,即可完成建筑物分割;
步骤4:主成分分析;
对分割后的建筑物点云进行主成分分析,计算出每个点的维度特征值;
步骤5:最佳邻域计算;
在不同尺度下分别进行主成分分析,计算不同尺度下维度特征值,再根据熵函数,自动计算出每个点的最佳邻域,得出每个点的主维度,从而将点云分为线状点云、面状点云、散乱点云三类;
步骤6:轮廓线提取;
根据几何位置关系,对面状点云进行轮廓线提取处理;先为每个点确定其所在的平面,再判断此点是在平面的边缘还是平面内部,根据平面中的点与其邻域点形成的线段夹角和约为360°来滤除平面中的点,从而保留每个平面的轮廓线点云。
2.根据权利要求1所述的地面激光点云数据的建筑物分割与轮廓线提取方法,其特征在于:步骤1中所述的对原始三维点云进行竖直投影,是将原始三维点云沿Z轴投影至X、Y平面上,仅保留每个点的X、Y坐标;其投影公式为:
3.根据权利要求1所述的地面激光点云数据的建筑物分割与轮廓线提取方法,其特征在于:步骤2中所述的统计每个格网内投影点个数,记为投影点密度DoPP,调整所有格网的DoPP的取值范围,将其线性拉伸到0-255,用拉伸后的值模拟每个格网的灰度,生成二维灰度图像。
4.根据权利要求1所述的地面激光点云数据的建筑物分割与轮廓线提取方法,其特征在于:步骤3中所述的Otsu算法对图像进行分割,计算出灰度阈值T,以阈值T将影像分为目标和背景两类,使它们的类间方差g为最大,g由公式(2)计算得出;
g=ω0ω1(μ0-μ1)2 (2);
其中ω0、ω1为两类各占比例,μ0、μ1为两类的平均灰度。
5.根据权利要求4所述的地面激光点云数据的建筑物分割与轮廓线提取方法,其特征在于,所述的Otsu算法的包括以下两点改进:
(1)分块求解;在进行Otsu算法的时候,并不是对所有的点云数据仅进行一次计算,而是依据建筑物立面在整个场景内的分布情况,对场景范围进行划分,在每个子场景内分别进行一次Otsu算法;划分原则为,尽可能多的满足,在每个子场景中,同时含有建筑物立面数据和其他地物数据,尤其是保证建筑物点云中DoPP值相对较小的部分满足上述要求,这样做可以显著提升点云分割的正确率;
将图像L在X和Y方向上分别划分为m和n等份,得到每个子场景的范围:
Lij=xi*yi(i=1,2,…,m;j=1,2,…,n) (3);
在每个Lij内进行一次Otsu算法,最后取每个子场景分割出的建筑物立面数据的合集,得到最终的分割结果;
(2)结合建筑物立面点云在竖直投影后可形成连续的线段这一特点,对于漏分割的点云部分,利用区域增长算法进行补偿;其具体的计算流程包括以下子步骤:
步骤3.1:选取种子点;设原始三维点云投影生成灰度图像的像素灰度为g[i][j],经过Otsu分块求解算法处理后像素灰度为G[i][j],设定标记数组B[i][j],记录每个像素的属性,三个数组的相同的下标对应同一个像素;遍历每个像素,若有像素p满足G[px][py]>0且其四邻域或八邻域有任意一个像素q符合G[qx][qy]=0且g[qx][qy]>0,那么B[px][py]=1,否则B[px][py]=0,其中像素p、q的关系满足公式(4);选取所有标记为1的像素作为种子点;
步骤3.2:设定区域增长条件,即阈值T;设种子点像素为p、其邻域点像素为q,将种子点G[px][py]与其邻域点g[qx][qy]作比较,若其商满足G[px][py]/g[qx][qy]<T,则将g[qx][qy]值赋予G[qx][qy],并将B[qx][qy]标记为1,若其商满足G[px][py]/g[qx][qy]≥T,则将B[qx][qy]标记为0,当p与其对应的所有q均完成上述操作时,将B[px][py]标记为0;
步骤3.3:区域增长;对B=1的像素看做新的种子点,并返回执行步骤3.2;
步骤3.4:结束条件;数组B中的所有元素均为0,即所有点的标记均为0,区域增长结束。
6.根据权利要求1所述的地面激光点云数据的建筑物分割与轮廓线提取方法,其特征在于:步骤4中所述的主成分分析是通过对每个点及其邻域点的三维坐标构建协方差矩阵,计算出矩阵的三个特征值,以此分析出每个点属于三个维度特征的概率a1D、a2D、a3D;其具体实现过程包括以下子步骤:
步骤4.1:根据点云密度和精度要求,设定邻域半径r的取值范围[rmin,rmax],以及取值间隔;
步骤4.2:在不同邻域半径下,分别进行主成分分析;
首先利用每个扫描点Xi及其邻域点νr的三维坐标{Xi=(xi,yi,zi)|i∈vr}构造协方差矩阵:
矩阵C为一个3*3的矩阵,其中 为点集νr的重心坐标,矩阵M的具体形式为:
然后计算矩阵C的三个特征值λ1、λ2、λ3,并按照λ1≥λ2≥λ3的规则进行排列;
步骤4.3:计算每个点的三个维度特征;
令根据公式(7)计算每个点的三个维度特征,即一维线状特征a1D、二维面状特征a2D和三维散乱点特征a3D;其中a1D、a2D、a3D的和为1,换言之,a1D、a2D、a3D分别代表了扫描点属于三个维度特征的概率;
7.根据权利要求6所述的地面激光点云数据的建筑物分割与轮廓线提取方法,其特征在于:步骤5中所述的最佳邻域计算,是根据熵函数计算每个点的最佳邻域,见公式(8):
Ef(vr)=-a1Dln(a1D)-a2Dln(a2D)-a3Dln(a3D) (8);
其中a1D、a2D、a3D是主成分分析中得到的每个点属于三个维度特征的概率;
根据公式(8)计算不同r值下的邻域点集的熵值;当熵取最小值时,表示在此邻域下该点的主维度特征最突出,此时对应的半径r*即为最佳邻域半径。
8.根据权利要求7所述的地面激光点云数据的建筑物分割与轮廓线提取方法,其特征在于:步骤6中所述的根据几何位置关系,对面状点云进行轮廓线提取处理;其具体实现过程包括以下子步骤:
步骤6.1:根据最佳邻域r*确定X点的最佳邻域点集
步骤6.2:确定点集所在平面;根据RANCAC随机抽样一致算法,从中随机抽取两个不同点并重复多次,重复次数与的点个数相同,计算X、Xi、Xj三点所在平面方程并统计中其余点到平面的欧氏距离和,选取使距离和最小的平面Sxpq的两邻域点Xp、Xq,将Xp、Xq从中删除;
步骤6.3:删除平面外部点;遍历计算每个Xi与Xp、Xq三点确定的平面Sipq方程,设定角度阈值θ1,当平面Sipq与平面Sxpq的夹角小于θ1时,认定Xi在平面Sxpq上,否则将Xi从中删除;
步骤6.4:设XXp为起始边,遍历计算每个∠XPXXi,设定角度阈值θ2,选取小于θ2中使∠XPXXi最小的邻域点X1';
步骤6.5:判断X1'与Xp的方向,如果X1'在Xp的顺时针方向,就按照顺时针方向去搜索;否则,按照逆时针搜索下一个邻域点;
步骤6.6:将步骤6.4得到的邻域点与X构成的线段设为起始边,重复步骤6.4寻求满足条件的下一个邻域点X2',X3',…,Xn',直到搜索不到满足条件的邻域点;
步骤6.7:计算的角度和,若和约等于360°,将点X视为平面中的点,将其滤除,否则视为边缘点,将其保留。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610393443.0A CN106056614A (zh) | 2016-06-03 | 2016-06-03 | 一种地面激光点云数据的建筑物分割与轮廓线提取方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610393443.0A CN106056614A (zh) | 2016-06-03 | 2016-06-03 | 一种地面激光点云数据的建筑物分割与轮廓线提取方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN106056614A true CN106056614A (zh) | 2016-10-26 |
Family
ID=57170323
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201610393443.0A Pending CN106056614A (zh) | 2016-06-03 | 2016-06-03 | 一种地面激光点云数据的建筑物分割与轮廓线提取方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106056614A (zh) |
Cited By (32)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106407925A (zh) * | 2016-09-09 | 2017-02-15 | 厦门大学 | 基于局部区间极大值的激光扫描点云树木自动提取方法 |
CN106846494A (zh) * | 2017-01-16 | 2017-06-13 | 青岛海大新星软件咨询有限公司 | 倾斜摄影三维建筑物模型自动单体化算法 |
CN107862738A (zh) * | 2017-11-28 | 2018-03-30 | 武汉大学 | 一种基于移动激光测量点云进行室内结构化三维重建方法 |
CN108241871A (zh) * | 2017-12-27 | 2018-07-03 | 华北水利水电大学 | 基于多特征的激光点云与影像融合数据分类方法 |
CN108734769A (zh) * | 2017-04-17 | 2018-11-02 | 宏达国际电子股份有限公司 | 三维模型分析方法、电子装置及非暂态电脑可读取媒体 |
CN108984599A (zh) * | 2018-06-01 | 2018-12-11 | 青岛秀山移动测量有限公司 | 一种利用行驶轨迹参考的车载激光点云路面提取方法 |
CN108986024A (zh) * | 2017-06-03 | 2018-12-11 | 西南大学 | 一种基于网格的激光点云规则排列处理方法 |
CN109118500A (zh) * | 2018-07-16 | 2019-01-01 | 重庆大学产业技术研究院 | 一种基于图像的三维激光扫描点云数据的分割方法 |
CN109325953A (zh) * | 2018-02-05 | 2019-02-12 | 黑龙江科技大学 | 一种大规模稠密点云法线的确定方法 |
CN110047036A (zh) * | 2019-04-22 | 2019-07-23 | 重庆交通大学 | 基于极坐标格网的地面激光扫描数据建筑物立面提取方法 |
CN110223267A (zh) * | 2018-03-14 | 2019-09-10 | 浙江大学山东工业技术研究院 | 基于高度直方图分割的耐火砖深度缺陷的识别方法 |
CN110377776A (zh) * | 2018-07-23 | 2019-10-25 | 北京京东尚科信息技术有限公司 | 生成点云数据的方法和装置 |
CN110969624A (zh) * | 2019-11-07 | 2020-04-07 | 哈尔滨工程大学 | 一种激光雷达三维点云分割方法 |
CN111179297A (zh) * | 2018-11-13 | 2020-05-19 | 北京航空航天大学 | 点云的多重轮廓生成方法、装置及*** |
CN111754385A (zh) * | 2019-03-26 | 2020-10-09 | 深圳中科飞测科技有限公司 | 数据点模型处理方法及***、检测方法及***、可读介质 |
CN111783721A (zh) * | 2020-07-13 | 2020-10-16 | 湖北亿咖通科技有限公司 | 一种激光点云的车道线提取方法及电子设备 |
CN111783722A (zh) * | 2020-07-13 | 2020-10-16 | 湖北亿咖通科技有限公司 | 一种激光点云的车道线提取方法和电子设备 |
CN111833286A (zh) * | 2019-03-26 | 2020-10-27 | 深圳中科飞测科技有限公司 | 点云处理方法及其***、检测方法及其***、可读介质 |
WO2020234678A1 (en) * | 2019-05-21 | 2020-11-26 | International Business Machines Corporation | Progressive 3d point cloud segmentation into object and background from tracking sessions |
CN112132969A (zh) * | 2020-09-01 | 2020-12-25 | 济南市房产测绘研究院(济南市房屋安全检测鉴定中心) | 一种车载激光点云建筑物目标分类方法 |
CN112232248A (zh) * | 2020-10-22 | 2021-01-15 | 中国人民解放军战略支援部队信息工程大学 | 一种多线LiDAR点云数据平面特征的提取方法及装置 |
CN112488010A (zh) * | 2020-12-05 | 2021-03-12 | 武汉中海庭数据技术有限公司 | 基于无人机点云数据的高精度目标提取方法及*** |
CN112529952A (zh) * | 2020-12-15 | 2021-03-19 | 武汉万集信息技术有限公司 | 物体体积测量方法、装置及电子设备 |
CN112595258A (zh) * | 2020-11-23 | 2021-04-02 | 扆亮海 | 基于地面激光点云的地物轮廓提取方法 |
CN112633657A (zh) * | 2020-12-16 | 2021-04-09 | 中冶建筑研究总院有限公司 | 一种施工质量管理方法、装置、设备及存储介质 |
CN113409332A (zh) * | 2021-06-11 | 2021-09-17 | 电子科技大学 | 一种基于三维点云的建筑物平面分割方法 |
CN113591597A (zh) * | 2021-07-07 | 2021-11-02 | 东莞市鑫泰仪器仪表有限公司 | 一种基于热成像的智慧治安信息*** |
CN113838072A (zh) * | 2021-11-01 | 2021-12-24 | 江苏集萃智能光电***研究所有限公司 | 一种高动态星图图像分割方法 |
CN114820986A (zh) * | 2022-05-13 | 2022-07-29 | 广西微车检智能科技有限公司 | 一种基于激光雷达的挂车外廓参数测量方法 |
WO2022198745A1 (en) * | 2021-03-24 | 2022-09-29 | Shanghai Huizi Cosmetics Co., Ltd. | Method and system for extracting nail contours |
US11532120B2 (en) | 2017-10-06 | 2022-12-20 | Interdigital Vc Holdings, Inc. | Method and device for hole filling of a point cloud |
CN115880325A (zh) * | 2022-12-07 | 2023-03-31 | 重庆市地理信息和遥感应用中心 | 基于点云维度和空间距离聚类的建筑物轮廓自动提取方法 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103020637A (zh) * | 2012-11-22 | 2013-04-03 | 北京航空航天大学 | 一种基于K-plane算法的建筑物顶面点云数据分割方法 |
CN104484668A (zh) * | 2015-01-19 | 2015-04-01 | 武汉大学 | 一种无人机多重叠遥感影像的建筑物轮廓线提取方法 |
-
2016
- 2016-06-03 CN CN201610393443.0A patent/CN106056614A/zh active Pending
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103020637A (zh) * | 2012-11-22 | 2013-04-03 | 北京航空航天大学 | 一种基于K-plane算法的建筑物顶面点云数据分割方法 |
CN104484668A (zh) * | 2015-01-19 | 2015-04-01 | 武汉大学 | 一种无人机多重叠遥感影像的建筑物轮廓线提取方法 |
Non-Patent Citations (2)
Title |
---|
秦家鑫等: "一种建筑物点云轮廓线的自动提取方法", 《遥感信息》 * |
秦家鑫等: "基于Otsu的建筑物点云分割改进算法", 《地理空间信息》 * |
Cited By (48)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106407925B (zh) * | 2016-09-09 | 2019-09-27 | 厦门大学 | 基于局部区间极大值的激光扫描点云树木自动提取方法 |
CN106407925A (zh) * | 2016-09-09 | 2017-02-15 | 厦门大学 | 基于局部区间极大值的激光扫描点云树木自动提取方法 |
CN106846494A (zh) * | 2017-01-16 | 2017-06-13 | 青岛海大新星软件咨询有限公司 | 倾斜摄影三维建筑物模型自动单体化算法 |
CN108734769A (zh) * | 2017-04-17 | 2018-11-02 | 宏达国际电子股份有限公司 | 三维模型分析方法、电子装置及非暂态电脑可读取媒体 |
CN108986024A (zh) * | 2017-06-03 | 2018-12-11 | 西南大学 | 一种基于网格的激光点云规则排列处理方法 |
CN108986024B (zh) * | 2017-06-03 | 2024-01-23 | 西南大学 | 一种基于网格的激光点云规则排列处理方法 |
US11532120B2 (en) | 2017-10-06 | 2022-12-20 | Interdigital Vc Holdings, Inc. | Method and device for hole filling of a point cloud |
CN107862738A (zh) * | 2017-11-28 | 2018-03-30 | 武汉大学 | 一种基于移动激光测量点云进行室内结构化三维重建方法 |
CN107862738B (zh) * | 2017-11-28 | 2019-10-11 | 武汉大学 | 一种基于移动激光测量点云进行室内结构化三维重建方法 |
CN108241871A (zh) * | 2017-12-27 | 2018-07-03 | 华北水利水电大学 | 基于多特征的激光点云与影像融合数据分类方法 |
CN109325953B (zh) * | 2018-02-05 | 2021-09-21 | 黑龙江科技大学 | 一种大规模稠密点云法线的确定方法 |
CN109325953A (zh) * | 2018-02-05 | 2019-02-12 | 黑龙江科技大学 | 一种大规模稠密点云法线的确定方法 |
CN110223267A (zh) * | 2018-03-14 | 2019-09-10 | 浙江大学山东工业技术研究院 | 基于高度直方图分割的耐火砖深度缺陷的识别方法 |
CN108984599A (zh) * | 2018-06-01 | 2018-12-11 | 青岛秀山移动测量有限公司 | 一种利用行驶轨迹参考的车载激光点云路面提取方法 |
CN108984599B (zh) * | 2018-06-01 | 2021-08-20 | 青岛秀山移动测量有限公司 | 一种利用行驶轨迹参考的车载激光点云路面提取方法 |
CN109118500B (zh) * | 2018-07-16 | 2022-05-10 | 重庆大学产业技术研究院 | 一种基于图像的三维激光扫描点云数据的分割方法 |
CN109118500A (zh) * | 2018-07-16 | 2019-01-01 | 重庆大学产业技术研究院 | 一种基于图像的三维激光扫描点云数据的分割方法 |
CN110377776A (zh) * | 2018-07-23 | 2019-10-25 | 北京京东尚科信息技术有限公司 | 生成点云数据的方法和装置 |
CN111179297A (zh) * | 2018-11-13 | 2020-05-19 | 北京航空航天大学 | 点云的多重轮廓生成方法、装置及*** |
CN111179297B (zh) * | 2018-11-13 | 2023-09-19 | 北京航空航天大学 | 点云的多重轮廓生成方法、装置及*** |
CN111754385A (zh) * | 2019-03-26 | 2020-10-09 | 深圳中科飞测科技有限公司 | 数据点模型处理方法及***、检测方法及***、可读介质 |
CN111833286A (zh) * | 2019-03-26 | 2020-10-27 | 深圳中科飞测科技有限公司 | 点云处理方法及其***、检测方法及其***、可读介质 |
CN111754385B (zh) * | 2019-03-26 | 2024-06-11 | 深圳中科飞测科技股份有限公司 | 数据点模型处理方法及***、检测方法及***、可读介质 |
CN110047036A (zh) * | 2019-04-22 | 2019-07-23 | 重庆交通大学 | 基于极坐标格网的地面激光扫描数据建筑物立面提取方法 |
GB2598512B (en) * | 2019-05-21 | 2022-10-05 | Ibm | Progressive 3D point cloud segmentation into object and background from tracking sessions |
WO2020234678A1 (en) * | 2019-05-21 | 2020-11-26 | International Business Machines Corporation | Progressive 3d point cloud segmentation into object and background from tracking sessions |
GB2598512A (en) * | 2019-05-21 | 2022-03-02 | Ibm | Progressive 3D point cloud segmentation into object and background from tracking sessions |
CN110969624A (zh) * | 2019-11-07 | 2020-04-07 | 哈尔滨工程大学 | 一种激光雷达三维点云分割方法 |
CN111783722A (zh) * | 2020-07-13 | 2020-10-16 | 湖北亿咖通科技有限公司 | 一种激光点云的车道线提取方法和电子设备 |
CN111783722B (zh) * | 2020-07-13 | 2021-07-06 | 湖北亿咖通科技有限公司 | 一种激光点云的车道线提取方法和电子设备 |
CN111783721A (zh) * | 2020-07-13 | 2020-10-16 | 湖北亿咖通科技有限公司 | 一种激光点云的车道线提取方法及电子设备 |
CN112132969B (zh) * | 2020-09-01 | 2023-10-10 | 济南市房产测绘研究院(济南市房屋安全检测鉴定中心) | 一种车载激光点云建筑物目标分类方法 |
CN112132969A (zh) * | 2020-09-01 | 2020-12-25 | 济南市房产测绘研究院(济南市房屋安全检测鉴定中心) | 一种车载激光点云建筑物目标分类方法 |
CN112232248A (zh) * | 2020-10-22 | 2021-01-15 | 中国人民解放军战略支援部队信息工程大学 | 一种多线LiDAR点云数据平面特征的提取方法及装置 |
CN112595258A (zh) * | 2020-11-23 | 2021-04-02 | 扆亮海 | 基于地面激光点云的地物轮廓提取方法 |
CN112595258B (zh) * | 2020-11-23 | 2022-04-22 | 湖南航天智远科技有限公司 | 基于地面激光点云的地物轮廓提取方法 |
CN112488010A (zh) * | 2020-12-05 | 2021-03-12 | 武汉中海庭数据技术有限公司 | 基于无人机点云数据的高精度目标提取方法及*** |
CN112529952A (zh) * | 2020-12-15 | 2021-03-19 | 武汉万集信息技术有限公司 | 物体体积测量方法、装置及电子设备 |
CN112529952B (zh) * | 2020-12-15 | 2023-11-14 | 武汉万集光电技术有限公司 | 物体体积测量方法、装置及电子设备 |
CN112633657B (zh) * | 2020-12-16 | 2024-05-14 | 中冶建筑研究总院有限公司 | 一种施工质量管理方法、装置、设备及存储介质 |
CN112633657A (zh) * | 2020-12-16 | 2021-04-09 | 中冶建筑研究总院有限公司 | 一种施工质量管理方法、装置、设备及存储介质 |
WO2022198745A1 (en) * | 2021-03-24 | 2022-09-29 | Shanghai Huizi Cosmetics Co., Ltd. | Method and system for extracting nail contours |
CN113409332A (zh) * | 2021-06-11 | 2021-09-17 | 电子科技大学 | 一种基于三维点云的建筑物平面分割方法 |
CN113591597A (zh) * | 2021-07-07 | 2021-11-02 | 东莞市鑫泰仪器仪表有限公司 | 一种基于热成像的智慧治安信息*** |
CN113838072A (zh) * | 2021-11-01 | 2021-12-24 | 江苏集萃智能光电***研究所有限公司 | 一种高动态星图图像分割方法 |
CN114820986A (zh) * | 2022-05-13 | 2022-07-29 | 广西微车检智能科技有限公司 | 一种基于激光雷达的挂车外廓参数测量方法 |
CN114820986B (zh) * | 2022-05-13 | 2024-04-09 | 广西微车检智能科技有限公司 | 一种基于激光雷达的挂车外廓参数测量方法 |
CN115880325A (zh) * | 2022-12-07 | 2023-03-31 | 重庆市地理信息和遥感应用中心 | 基于点云维度和空间距离聚类的建筑物轮廓自动提取方法 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106056614A (zh) | 一种地面激光点云数据的建筑物分割与轮廓线提取方法 | |
Jin et al. | Deep learning: individual maize segmentation from terrestrial lidar data using faster R-CNN and regional growth algorithms | |
CN104049245B (zh) | 基于LiDAR点云空间差异分析的城市建筑物变化检测方法 | |
CN106372277B (zh) | 森林立地指数时空估测中的变异函数模型优化方法 | |
CN109711410A (zh) | 一种三维物体快速分割和识别方法、装置及*** | |
CN111340875B (zh) | 一种基于三维激光雷达的空间移动目标检测方法 | |
CN105469098A (zh) | 一种自适应特征权重合成的lidar数据地物精确分类方法 | |
CN107316048A (zh) | 点云分类方法及装置 | |
CN106199557A (zh) | 一种机载激光雷达数据植被提取方法 | |
Marinelli et al. | A novel approach to 3-D change detection in multitemporal LiDAR data acquired in forest areas | |
CN104091321A (zh) | 适用于地面激光雷达点云分类的多层次点集特征的提取方法 | |
CN103679655A (zh) | 一种基于坡度与区域生长的LiDAR点云滤波方法 | |
CN103745441A (zh) | 一种机载LiDAR点云滤波方法 | |
Stone et al. | Determining an optimal model for processing lidar data at the plot level: results for a Pinus radiata plantation in New South Wales, Australia | |
CN110532963B (zh) | 一种车载激光雷达点云驱动的道路标线精准提取方法 | |
CN105809194A (zh) | 一种sar影像翻译为光学影像的方法 | |
CN108109139A (zh) | 基于灰度体元模型的机载lidar三维建筑物检测方法 | |
CN112257605A (zh) | 基于自标注训练样本的三维目标检测方法、***及装置 | |
CN110363771B (zh) | 一种基于三维点云数据的隔离护栏形点提取方法及装置 | |
CN114037836A (zh) | 一种将人工智能识别技术应用于三维输变电工程测算的方法 | |
CN107292039B (zh) | 一种基于小波聚类的uuv巡岸轮廓构建方法 | |
Ayrey et al. | Ecologically-based metrics for assessing structure in developing area-based, enhanced forest inventories from LiDAR | |
Pawłuszek et al. | Landslides identification using airborne laser scanning data derived topographic terrain attributes and support vector machine classification | |
CN117765006A (zh) | 基于无人机影像与激光点云的多层次密集树冠分割方法 | |
Cherlinka | Using geostatistics, DEM and remote sensing to clarify soil cover maps of Ukraine |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
RJ01 | Rejection of invention patent application after publication |
Application publication date: 20161026 |
|
RJ01 | Rejection of invention patent application after publication |