CN108389250B - 基于点云数据快速生成建筑物断面图的方法 - Google Patents
基于点云数据快速生成建筑物断面图的方法 Download PDFInfo
- Publication number
- CN108389250B CN108389250B CN201810191340.5A CN201810191340A CN108389250B CN 108389250 B CN108389250 B CN 108389250B CN 201810191340 A CN201810191340 A CN 201810191340A CN 108389250 B CN108389250 B CN 108389250B
- Authority
- CN
- China
- Prior art keywords
- point cloud
- point
- cloud data
- points
- plane
- 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.)
- Expired - Fee Related
Links
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T15/00—3D [Three Dimensional] image rendering
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01B—MEASURING LENGTH, THICKNESS OR SIMILAR LINEAR DIMENSIONS; MEASURING ANGLES; MEASURING AREAS; MEASURING IRREGULARITIES OF SURFACES OR CONTOURS
- G01B11/00—Measuring arrangements characterised by the use of optical techniques
- G01B11/002—Measuring arrangements characterised by the use of optical techniques for measuring two or more coordinates
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Computer Graphics (AREA)
- Theoretical Computer Science (AREA)
- Processing Or Creating Images (AREA)
- Image Analysis (AREA)
Abstract
本发明提供基于点云数据快速生成建筑物断面图的方法:测量采集建筑物的点云数据,计算坐标转换参数,进行点云配准和去噪;用裁剪立方体截取建筑物断面,得到边缘点云数据,形成二维平面数据P(xi,yi)T;对点云数据进行平移和放大,提取共线点云数据Pk1(xi,yi)T;利用角度θk和极径ρk剔除Pk1(xi,yi)T中的干扰点,得到共线点云数据Pk2(xi,yi)T;对其进行排序和分割,找出点数最多的一段保留下来,得到点云线段数据Pk3(xi,yi)T,将Pk3(xi,yi)T从P(xi,yi)T中删除;循环执行以上步骤提取点云线段,当点云线段的点数小于点数阈值t4时结束循环;将点云线段数据Pk3(xi,yi)T抽稀得到Pk4(xi,yi)T,再采用基于选权迭代的最小二乘法进行精确直线拟合,得到各线段的端点坐标;将各条线段的端点坐标首尾相连,得到建筑物断面图。
Description
技术领域
本发明属于三维激光扫描技术领域,具体涉及一种基于点云数据快速生成建筑物断面图的方法。
背景技术
工程应用中常需要获得建筑物的断面图,传统的建筑物断面图是利用全站仪等设备采集各个独立点的三维坐标,然后通过软件或人工绘制,该方法外业工作量大,效率低且内业处理过程繁琐,只能获得局部信息。
三维激光扫描是近年发展起来的一门新技术,以高精度、海量、快速、无接触测量等优势成为空间三维数据获取的有效手段,而获取的点云数据是没有直接拓扑关系的一系列离散点,要想利用点云数据获得建筑物断面图,首先需要将断面中各条线段检测出来,传统的直线检测方法主要是霍夫变换,而直接利用霍夫变换进行点云直线检测只能达到有限精度,特别是受噪声影响时,如果对精度要求过高可能导致检测失败。最小二乘法则可获得均方误差意义下的精确直线,但是容易受到噪声点干扰,如果直接进行最小二乘拟合会受到强噪声干扰,达不到获取高精度直线参数的目的。
发明内容
本发明是为了解决水光互补电站安全经济运行问题而进行的,目的在于提供一种基于点云数据快速生成建筑物断面图的方法,实现建筑物断面图快速精确绘制。
本发明为了实现上述目的,采用了以下方案:
本发明提供一种基于点云数据快速生成建筑物断面图的方法,其特征在于,包括以下步骤:步骤1.分站测量采集建筑物的点云数据,利用标靶坐标计算相邻两站之间的坐标转换参数,进而进行点云配准和去噪;步骤2.用裁剪立方体对建筑物断面进行截取,得到该断面的边缘点云数据,形成断面边缘点云的二维平面数据P(xi,yi)T;步骤3.对点云数据进行平移和放大,然后提取共线点云数据Pk1(xi,yi)T,其中1≤i≤N,N为集合P(xi,yi)T点数,包含:步骤3-1.对二维点云集合P(xi,yi)T,首先计算出中心位置(Xm,Ym),然后将点云数据整体平移,使中心位置与原点重合,再将点云坐标整体放大,利用转换后坐标算出向上取整后的极径最大值ρmax;步骤3-2.设置两个m×n累加器矩阵,m为将180度等分的数目,n=ρmax,计算每个θi下对应的所有点向上取整后的极径ρi,j=xjcosθi+yjsinθi,其中1≤j≤N,N为集合P(xi,yi)T点数,找到i和ρ所对应的单元进行投票,并将该单元的累加器加1;由于平移后计算出的极径ρ值有时为负,所以需要设置两个累加器矩阵;步骤3-3.找出步骤3-2中两个累加器中值最大的单元所对应的θk=180*i/m和ρk,其中1≤i≤m,1≤k≤R,R为提取出的线段总数,θk和ρk分别为极坐标系中的角度和极径;设置阈值t1,逐点计算角度θk下点云对应的极径值ρkj=xjcosθk+yjsinθk,其中1≤j≤N,N为集P(xi,yi)T点数;如果-t1≤ρkj-ρk≤t1,认为点在该条直线上,将其存储起来继续计算下一点,直至计算完毕提取出一条共线点云集合Pk1(xi,yi)T;步骤4.利用步骤3得到的直线参数角度θk和极径ρk剔除Pk1(xi,yi)T中的干扰点,得到共线点云数据Pk2(xi,yi)T;步骤5.对共线点云进行排序和分割,找出点数最多的一段保留下来,得到点云线段数据Pk3(xi,yi)T,将Pk3(xi,yi)T从P(xi,yi)T中删除;步骤6.循环执行以上步骤3至5提取点云线段,设定点数阈值t4作为点云线段的最少点数,当点云线段的点数小于该点数阈值t4时结束循环,至此提取出全部的点云线段;步骤7.对于每条点云线段,将点云线段数据Pk3(xi,yi)T按一定间隔进行抽稀得到Pk4(xi,yi)T,然后采用基于选权迭代的最小二乘法进行精确直线拟合,得到各条线段的端点坐标;步骤8.求得各条线段的端点坐标后,将其首尾相连,并可添加图廓注释等信息,得到建筑物断面图。
进一步地,本发明提供的基于点云数据快速生成建筑物断面图的方法,还可以具有以下特征:在步骤1中先利用全站仪布设控制网,获取控制点的绝对坐标,再将地面三维激光扫描仪分站架设在控制点上采集点云数据,用扫描仪可识别的标靶作为每两站之间的公共点,保证两站之间至少有四个共同标靶。
进一步地,本发明提供的基于点云数据快速生成建筑物断面图的方法,还可以具有以下特征:在步骤2中,裁剪立方体是由六个平面相交而成的实体,在计算各平面方程时规定法向量朝向立方体内部为正方向;由起点A(X1,Y1,Z1),终点B(X2,Y2,Z2)和边缘点C(X3,Y3,Z3)确定一个裁剪面ABC,AB⊥AC,由DEFG-HIJK确定该裁剪立方体,DHIE和GKJF分别表示裁剪立方体的前后平面且与XOY平面垂直,DHKG和EIJF分别表示裁剪立方体的左右平面且同时垂直于裁剪面和前后平面,DEFG和HIJK分别表示上下平面且分别平行于裁剪面;点C为DH中点,AC长为裁剪宽度B的一半,DH长为裁剪厚度T;
已知
裁剪面ABC的法向量为:
前后平面DHIE和GKJF分别与XOY平面垂直,因此法向量Z坐标为0;令dX=X2-X1,dY=Y2-Y1,则前后平面的法向量分别为:
由于前平面过点A(X1,Y1,Z1),后平面过点B(X2,Y2,Z2),根据平面点法式方程,可以求得各自的平面方程为:
左右平面DHKG和EIJF指距离裁剪面起点和终点连线半个宽度,且各自平行于该连线的两个平面。由于左右平面既垂直于前后面,又垂直于裁剪面,因此左右平面的法向量同时垂直于前后平面和裁剪平面的法向量,那么可以求得左右平面DHKG和EIJF的的法向量为:
设裁剪面起点和终点连线的中点坐标为(Xc,Yc,Zc),则:
令
则左右平面的方程分别为:
令
上下平面DEFG和HIJK的方程分别为:
假设一个点的坐标为(X1,Y1,Z1),所在平面的方程为
Ax+By+Cz+D=0 (10)
令
S=AX1+BY1+CZ1+D (11)
规定:
将待截取范围内的点云数据逐个带入以上公式(3)、(7)和(9)中,当点都在六个平面的正方向时,认为该点是要截取的点,存储下来继续计算下一点,直至计算完毕,得到断面边缘的二维点云集合P(xi,yi)T。从而实现建筑物断面点云数据的获取。
进一步地,本发明提供的基于点云数据快速生成建筑物断面图的方法,还可以具有以下特征:在步骤3中,是将点云坐标整体放大1000倍,这样在后续投票过程中可以精确到小数点后三位,可根据实际情况进行调整。
进一步地,本发明提供的基于点云数据快速生成建筑物断面图的方法,还可以具有以下特征:步骤4包含以下子步骤:
步骤4-1.设集合Pk1(xi,yi)T所在直线的斜截式方程为:
yi=akxi+bk (12)
计算集合Pk1(xi,yi)T中的点到公式(12)所确定的直线的距离dki:
式中(xi,yi)∈Pk1,1≤i≤Nk1,Nk1为集合Pk1(xi,yi)T点数;
步骤4-2.设定一个阈值t2作为点到直线的最远距离(0<t2<t1),如果dki≤t2,则有:
式中1≤i≤Nk2,Nk2≤Nk1,Nk2为集合Pk2(xi,yi)T点数;
步骤4-3.将距离直线不大于t2的点存储在集合Pk2(xi,yi)T中,从而有效去除线段连接处的干扰点。
进一步地,本发明提供的基于点云数据快速生成建筑物断面图的方法,还可以具有以下特征:步骤5包含以下子步骤:
步骤5-1.利用直线角度θk对Pk2(xi,yi)T坐标进行排序:如果0°≤θk≤45°或145°≤θk<180°,按照X坐标从小到大的顺序进行排序;如果45°<θk<145°,则按照Y坐标从小到大的顺序进行排序;
步骤5-2.逐点计算排序后的两点之间距离L1,设定距离阈值t3作为相邻两点间的最远距离;如果L1≤t3,认为两点在一条线段上,继续计算下一点;如果L1>t3,认为该点与前面点不在一条线段上,把前面点存在一个数组中,从该点继续计算,直到将共线点云分割完毕;
步骤5-3.比较步骤5-2中分割后各条点云线段的点数,取点数最多的一段作为提取出的第一条点云线段集合Pk3(xi,yi)T;
步骤5-4.将Pk3(xi,yi)T从集合P(xi,yi)T中删除。
进一步地,本发明提供的基于点云数据快速生成建筑物断面图的方法,还可以具有以下特征:在步骤7中,将集合Pk4(xi,yi)T的直线方程表示为:
yi=axi+b,
式中i=1,2…n,(xi,yi)为Pk4(xi,yi)T中的点,a为直线的斜率,b为y轴的截距,a、b为待估参数;在进行选权迭代之前应先利用如下公式计算出a、b的近似值a0、b0:
X=(BTB)-1BTL (15)
令a=a0+δa,b=b0+δb,以x为自变量,y为因变量,误差方程为:
可进一步表达为:
Vk=BδXk-lk,
按照最小二乘准则:
Vk TPkVk=min,
当k=1即初次进行最小二乘运算时,令
有
接下来根据单位权中误差σ给每个点重新定权,一般采用丹麦权,其作用是对中误差过大的点进行降权,以降低其对平差结果的影响,丹麦权有很多形式,这里采用如下一种:
得到重新定权后的权阵P2,代入公式(16)计算δX2、V2和σ,如此循环迭代,设定一个阈值t6,当|δak-δak-1|<t6且|δbk-δbk-1|<t6时结束迭代;令a=a0+δak,b=b0+δbk,此时求得的a,b即为最终的拟合参数。
进一步地,本发明提供的基于点云数据快速生成建筑物断面图的方法,还可以具有以下特征:在步骤7中,在线段角度小于45度或大于145度的情况下,对X和Y坐标采用最小二乘法拟合出直线方程,利用两端点的X坐标,代入直线方程求出Y坐标;如果线段角度大于45度且小于145度,对Y和X坐标采用最小二乘法拟合出直线方程,利用两端点的Y坐标,代入直线方程求出X坐标,最后将得到的端点坐标再平移缩放回去即可。
进一步地,本发明提供的基于点云数据快速生成建筑物断面图的方法,还可以具有以下特征:在步骤7中,可以每隔三个点或四个点等(根据点云密度情况),抽稀的目的是减少选权迭代时的次数,从而提高处理效率。由于提取出的点云线段Pk3(xi,yi)T已按坐标排好序,因此按一定间隔进行抽稀得到Pk4(xi,yi)T,不仅不会改变其形状,由于数据量的减少还会加快拟合速度。
进一步地,本发明提供的基于点云数据快速生成建筑物断面图的方法,还可以具有以下特征:在步骤7中,在抽稀得到Pk4(xi,yi)T后,设定长度阈值t5,用于对过长的线段进行分割,计算Pk4(xi,yi)T点云线段的长度L2,如果L2>t5,则以t5为尺度将L分成若干段,再对每一小段采用基于选权迭代的最小二乘法进行精确直线拟合,得到各条线段的端点坐标。分割后直线拟合效果更好。
发明的作用与效果
本发明所提出的基于点云数据快速生成建筑物断面图的方法,利用点云数据自动生成建筑物断面图,突破了传统断面图的绘制方法;改进了霍夫变换算法,使其适用于二维点云数据的直线提取,提高了直线提取的效率;采用选权迭代的最小二乘法拟合点云线段,有效消除了强噪声点的干扰,实现点云线段的精确拟合。
附图说明
图1为本发明实施例涉及的基于点云数据快速生成建筑物断面图的方法的流程图;
图2为本发明实施例涉及的裁剪立方体和裁剪面的结构示意图;
图3为本发明实施例涉及的截取变截面的示意图;
图4为本发明实施例涉及的变截面点云平面图;
图5为本发明实施例中提取出来的点云线段的示意图;
图6为本发明实施例涉及的隧道变截面纵断面图;
图7为本发明实施例涉及的叠加显示图;
图8为本发明实施例涉及的断面叠加局部放大图,其中(a)和(b)分别是左下角和右下角的断面叠加放大图。
具体实施方式
以下结合附图对本发明涉及的基于点云数据快速生成建筑物断面图的方法的具体实施方案进行详细地说明。
<实施例>
如图1所示,本实施例所提供的基于点云数据快速生成建筑物断面图的方法包括以下步骤:
步骤1.采用徕卡ScanStation P30三维激光扫描仪对武汉地铁2号线的一段隧道进行分段扫描后,利用标靶坐标计算两站之间的坐标转换参数,对获取的点云数据进行拼接、建立索引和点云可视化。
步骤2.取一个与隧道设计中心线平行的裁剪面,该裁剪面宽度为15m(比隧道宽度略大),厚度为0.02m,长度比变截面要略长,使变截面在该裁剪面之内,利用该裁剪面去截取该隧道相应位置的变截面,得到该断面边缘点云的二维平面数据P(xi,yi)T。图3给出了截取变截面示意图。图3中,中间平面部分为点云截取范围,包围了整个变截面,利用本文介绍的方法获取该部分的纵断面点云数据,共获得73272个点,图4给出了截取的变截面纵断面边缘点云数据的平面效果图。由图4可以看出,它与原始变截面轮廓相对应。
步骤3.提取共线点云数据Pk1(xi,yi)T。设置累加器矩阵,其大小为900×1000Rmax,即将180°划分成900份,Rmax表示原始变截面点云数据左下角到右上角距离向上取整后的值,找出中心位置(X0,Y0),将点云数据整体平移,使其中心与原点重合,再将点云坐标整体放大1000倍,作用是在投票过程中精确到小数点后三位。极径的阈值t1取0.03m。这里不同于图像的原点坐标一般在图像边缘,点云数据坐标通常距离原点较远,如果直接进行霍夫变化将占用大量内存,大大增加处理时间,有时会超出计算机的处理能力导致崩溃,又由于霍夫变换过程中计算极径时需要向上取整,势必会产生较大的误差,这样无法分别出差别较小的直线,因此本步骤中先对点云数据进行适当的平移和放大,再进行改进的霍夫变换。
步骤4.利用步骤3得到的直线参数角度θ和极径ρ剔除Pk1(xi,yi)T中的干扰点,得到共线点云数据点云Pk2(xi,yi)T。线段去噪时到直线距离的阈值t2取0.015m。
步骤5.对共线点云进行排序和分割,找出点数最多的一段保留下来,得到点云线段数据Pk3(xi,yi)T,将Pk3(xi,yi)T从P(xi,yi)T中删除。
步骤6.重复以上步骤3~5,提取下一条点云线段。设定一个阈值作为点云线段的最少点数,当点云线段的点数小于该阈值时结束循环。至此提取出全部的点云线段。点云分割时的阈值t3取0.15m,点云线段的最少个数的阈值t4取20个。图5是提取出来的点云线段,用不同的灰度将各条线段区分开来。
步骤7.将Pk3(xi,yi)T按一定间隔抽稀得到Pk4(xi,yi)T,设定一个长度阈值,用于对过长的线段进行分割,再对每一小段采用基于选权迭代的最小二乘法进行精确直线拟合,得到各条线段的端点坐标。分段拟合时的长度阈值t5取0.5m,因坐标放大了1000倍,t1、t2、t3和t5也应放大相应倍数。
步骤8.求得各条线段的端点坐标后,将其首尾相连,即可得到该变截面的纵断面图,如图6所示。
由图5、6可以看出:本方法成功提取出了大部分线段,各条线段之间的界限清晰,说明本方法既利用了霍夫变换在直线提取领域的优势,又有效解决了霍夫变换提取直线时精度不够的问题,同时该算法也兼顾了处理效率,对于数万个点云数据,可在数十秒内提取并拟合出大部分线段。自动生成的隧道变截面纵断面图与图5轮廓相对应。需要指出的是:平面上的一条直线在理想条件下对应参数空间中一个点,但是点云数据是一些不连续的离散点组成,又由于对参数空间进行了数字化,因此这些离散点通过霍夫变换在参数空间中对应的曲线也是由一些离散点组成,它们有时不会只相交于一点,从提取出来的点云线段可以看出,部分线段被分成了若干段,这就是所谓的峰值扩散现象,但并不影响最终的拟合效果,如有需要可根据两条线段的直线参数及最近距离等条件予以合并。
为了进一步考察断面图生成效果,将图6隧道变截面纵断面图与图4变截面点云平面图进行叠加显示,见图7,为了能清楚地显示效果,图7的左下角和右下角部分局部放大显示,其效果见图8。
由图7和8可以看出:直线拟合时采用的选权迭代算法具有良好的抗差效果,有效地消除了噪声点的干扰,拟合的直线与原始点云吻合良好。由于去除了部分干扰点,又对点云数据进行了抽稀处理,提高了选权迭代时的效率,可快速自动生成断面图。
以上实施例仅仅是对本发明技术方案所做的举例说明。本发明所涉及的基于点云数据快速生成建筑物断面图的方法并不限定于在以上实施例中所描述的内容,而是以权利要求所限定的范围为准。本发明所属领域技术人员在该实施例的基础上所做的任何修改或补充或等效替换,都在本发明的权利要求所要求保护的范围内。
Claims (9)
1.一种基于点云数据快速生成建筑物断面图的方法,其特征在于,包括以下步骤:
步骤1.分站测量采集建筑物的点云数据,利用标靶坐标计算相邻两站之间的坐标转换参数,进而进行点云配准和去噪;
步骤2.用裁剪立方体对建筑物断面进行截取,得到该断面的边缘点云数据,形成断面边缘点云的二维平面数据P(xi,yi)T;
步骤3.对点云数据进行平移和放大,然后提取共线点云数据Pk1(xi,yi)T,其中1≤i≤N,N为集合P(xi,yi)T点数,包含:
步骤3-1.对二维点云集合P(xi,yi)T,首先计算出中心位置(xm’,ym’),然后将点云数据整体平移,使中心位置与原点重合,再将点云坐标整体放大,利用转换后坐标算出向上取整后的极径最大值ρmax;
步骤3-2.设置两个m×n累加器矩阵,m为将180度等分的数目,n=ρmax,计算每个θe下对应的所有点向上取整后的极径ρe,j=xjcosθe+yjsinθe,其中1≤j≤N,N为集合P(xi,yi)T点数,θe=180°*e/m,1≤e≤m,找到e和ρe,j所对应的单元,并将该单元的累加器加1;
步骤3-3.找出步骤3-2中两个累加器中值最大的单元所对应的θk和ρk,其中1≤k≤R,R为提取出的线段总数,θk和ρk分别为极坐标系中的角度和极径;设置阈值t1,逐点计算角度θk下点云对应的极径值ρk,j=xjcosθk+yjsinθk,其中1≤j≤N,N为集合P(xi,yi)T点数;如果-t1≤ρk,j-ρk≤t1,认为点在直线上,将其存储起来继续计算下一点,直至计算完毕提取出一条共线点云集合Pk1(xi,yi)T;
步骤4.利用步骤3得到的直线参数角度θk和极径ρk剔除Pk1(xi,yi)T中的干扰点,得到共线点云数据Pk2(xi,yi)T;
步骤5.对共线点云进行排序和分割,找出点数最多的一段保留下来,得到点云线段数据Pk3(xi,yi)T,将Pk3(xi,yi)T从P(xi,yi)T中删除;
步骤6.循环执行以上步骤3至5提取点云线段,设定点数阈值t4作为点云线段的最少点数,当点云线段的点数小于该点数阈值t4时结束循环,至此提取出全部的点云线段;
步骤7.对于每条点云线段,将点云线段数据Pk3(xi,yi)T按一定间隔进行抽稀得到Pk4(xi,yi)T,然后采用基于选权迭代的最小二乘法进行精确直线拟合,得到各条线段的端点坐标;
步骤8.求得各条线段的端点坐标后,将其首尾相连,得到建筑物断面图。
2.根据权利要求1所述的基于点云数据快速生成建筑物断面图的方法,其特征在于:
其中,在步骤1中先利用全站仪布设控制网,获取控制点的绝对坐标,再将地面三维激光扫描仪分站架设在控制点上采集点云数据,用扫描仪可识别的标靶作为每两站之间的公共点,保证两站之间至少有四个共同标靶。
3.根据权利要求1所述的基于点云数据快速生成建筑物断面图的方法,其特征在于:
其中,在步骤2中,裁剪立方体是由六个平面相交而成的实体,在计算各平面方程时规定法向量朝向立方体内部为正方向;由起点A(x1,y1,z1),终点B(x2,y2,z2)和边缘点C(x3,y3,z3)确定一个裁剪面ABC,AB⊥AC,由DEFG-HIJK确定该裁剪立方体,DHIE和GKJF分别表示裁剪立方体的前后平面且与XOY平面垂直,DHKG和EIJF分别表示裁剪立方体的左右平面且同时垂直于裁剪面和前后平面,DEFG和HIJK分别表示上下平面且分别平行于裁剪面;点C为DH中点,AC长为裁剪宽度W的一半,DH长为裁剪厚度T;
裁剪面ABC的法向量为:
令dX=x2-x1,dY=y2-y1,则前后平面的法向量分别为:
由于前平面过点A(x1,y1,z1),后平面过点B(x2,y2,z2),根据平面点法式方程,求得各自的平面方程为:
左右平面DHKG和EIJF的的法向量为:
设裁剪面起点和终点连线的中点坐标为(xc,yc,zc),则
令
则左右平面的方程分别为:
令
上下平面DEFG和HIJK的方程分别为:
假设平面的方程为
A′x+B′y+C′z+D′=0 (10)
空间中一个点的坐标为(xp,yp,zp),
令
S=A′xp+B′yp+C′zp+D′ (11)
规定:S>0表示点在正面,S=0表示点平面上,S<0表示点在背面;
将待截取范围内的点云数据逐个代入以上公式(3)、(7)和(9)中,当点都在六个平面的正方向时,认为该点是要截取的点,存储下来继续计算下一点,直至计算完毕,得到断面边缘的二维点云集合P(xi,yi)T。
4.根据权利要求1所述的基于点云数据快速生成建筑物断面图的方法,其特征在于:
其中,在步骤3中,是将点云坐标整体放大1000倍。
6.根据权利要求1所述的基于点云数据快速生成建筑物断面图的方法,其特征在于:
其中,步骤5包含以下子步骤:
步骤5-1.利用直线角度θk对Pk2(xi,yi)T坐标进行排序:如果0°≤θk≤45°或145°≤θk<180°,按照X坐标从小到大的顺序进行排序;如果45°<θk<145°,则按照Y坐标从小到大的顺序进行排序;
步骤5-2.逐点计算排序后的两点之间距离L1,设定距离阈值t3作为相邻两点间的最远距离;如果L1≤t3,认为两点在一条线段上,继续计算下一点;如果L1>t3,认为当前计算点与前面点不在一条线段上,把前面点存在一个数组中,从当前计算点开始继续计算,直到将共线点云分割完毕;
步骤5-3.比较步骤5-2中分割后各条点云线段的点数,取点数最多的一段作为提取出的第一条点云线段集合Pk3(xi,yi)T;
步骤5-4.将Pk3(xi,yi)T从集合P(xi,yi)T中删除。
7.根据权利要求1所述的基于点云数据快速生成建筑物断面图的方法,其特征在于:
其中,在步骤7中,将集合Pk4(xi,yi)T的直线方程表示为:
yi=axi+b (14)
式中i=1,2…n,(xi,yi)为Pk4(xi,yi)T中的点,a为直线的斜率,b为y轴的截距,a、b为待估参数;在进行选权迭代之前应先利用如下公式计算出a、b的近似值a0、b0:
X=(BTB)-1BTL (15)
令a=a0+δa,b=b0+δb,以x为自变量,y为因变量,误差方程为:
进一步表达为:
Vk=BδXk-lk,
按照最小二乘准则:
Vk TPkVk=min,
当k=1即初次进行最小二乘运算时,令
有
根据单位权中误差σ给每个点重新定权,得到重新定权后的权阵P2,代入公式(16)计算δX2、V2和σ,如此循环迭代,设定一个阈值t6,当|δak-δak-1|<t6且|δbk-δbk-1|<t6时结束迭代;令a=a0+δak,b=b0+δbk,此时求得的a,b即为最终的拟合参数。
8.根据权利要求7所述的基于点云数据快速生成建筑物断面图的方法,其特征在于:
其中,在步骤7中,在线段角度小于45度或大于145度的情况下,对X和Y坐标采用最小二乘法拟合出直线方程,利用两端点的X坐标,代入直线方程求出Y坐标;如果线段角度大于45度且小于145度,对Y和X坐标采用最小二乘法拟合出直线方程,利用两端点的Y坐标,代入直线方程求出X坐标。
9.根据权利要求1所述的基于点云数据快速生成建筑物断面图的方法,其特征在于:
其中,在步骤7中,在抽稀得到Pk4(xi,yi)T后,以设定的长度阈值t5为尺度对过长的线段进行分割,再对每一小段采用基于选权迭代的最小二乘法进行精确直线拟合,得到各条线段的端点坐标。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810191340.5A CN108389250B (zh) | 2018-03-08 | 2018-03-08 | 基于点云数据快速生成建筑物断面图的方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810191340.5A CN108389250B (zh) | 2018-03-08 | 2018-03-08 | 基于点云数据快速生成建筑物断面图的方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN108389250A CN108389250A (zh) | 2018-08-10 |
CN108389250B true CN108389250B (zh) | 2020-05-22 |
Family
ID=63067075
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201810191340.5A Expired - Fee Related CN108389250B (zh) | 2018-03-08 | 2018-03-08 | 基于点云数据快速生成建筑物断面图的方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN108389250B (zh) |
Families Citing this family (12)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110148144B (zh) * | 2018-08-27 | 2024-02-13 | 腾讯大地通途(北京)科技有限公司 | 点云数据的分割方法和装置、存储介质、电子装置 |
CN108931983B (zh) * | 2018-09-07 | 2020-04-24 | 深圳市银星智能科技股份有限公司 | 地图构建方法及其机器人 |
CN109668543A (zh) * | 2019-01-22 | 2019-04-23 | 南京理工大学 | 基于激光雷达的倾斜度测量方法 |
CN109949349B (zh) * | 2019-01-24 | 2021-09-21 | 北京大学第三医院(北京大学第三临床医学院) | 一种多模态三维图像的配准及融合显示方法 |
CN110095062B (zh) * | 2019-04-17 | 2021-01-05 | 北京华捷艾米科技有限公司 | 一种物体体积参数测量方法、装置及设备 |
CN111563892B (zh) * | 2020-05-11 | 2024-01-26 | 江苏易图地理信息科技股份有限公司 | 基于点云数据立面采集精度优化的拟合算法 |
CN112232248B (zh) * | 2020-10-22 | 2023-04-07 | 中国人民解放军战略支援部队信息工程大学 | 一种多线LiDAR点云数据平面特征的提取方法及装置 |
CN112164145B (zh) * | 2020-10-30 | 2022-04-29 | 武汉大学 | 基于点云数据快速提取室内三维线段结构的方法 |
CN114880332B (zh) * | 2022-07-08 | 2022-09-16 | 深圳市信润富联数字科技有限公司 | 点云数据的存储方法、装置、电子设备及存储介质 |
CN115201785B (zh) * | 2022-09-15 | 2022-12-13 | 建设综合勘察研究设计院有限公司 | 一种建筑物数字断面提取方法、系统、终端及存储介质 |
CN115661220B (zh) * | 2022-12-28 | 2023-03-17 | 深圳煜炜光学科技有限公司 | 点云数据配准方法、装置、设备及存储介质 |
CN116626706B (zh) * | 2023-05-12 | 2024-01-16 | 北京交通大学 | 轨道交通隧道侵限检测方法及系统 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105184855A (zh) * | 2015-08-25 | 2015-12-23 | 广州市城市规划勘测设计研究院 | 基于三维点云的特征面构建方法及装置 |
CN106548479A (zh) * | 2016-12-06 | 2017-03-29 | 武汉大学 | 一种多层次激光点云建筑物边界规则化方法 |
CN107462204A (zh) * | 2017-09-21 | 2017-12-12 | 武汉武大卓越科技有限责任公司 | 一种三维路面标准轮廓提取方法及系统 |
-
2018
- 2018-03-08 CN CN201810191340.5A patent/CN108389250B/zh not_active Expired - Fee Related
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105184855A (zh) * | 2015-08-25 | 2015-12-23 | 广州市城市规划勘测设计研究院 | 基于三维点云的特征面构建方法及装置 |
CN106548479A (zh) * | 2016-12-06 | 2017-03-29 | 武汉大学 | 一种多层次激光点云建筑物边界规则化方法 |
CN107462204A (zh) * | 2017-09-21 | 2017-12-12 | 武汉武大卓越科技有限责任公司 | 一种三维路面标准轮廓提取方法及系统 |
Non-Patent Citations (1)
Title |
---|
A New Progressive Simplification Method for Point Cloud Using Local Entropy of Normal Angle;Wei Xuan等;《Journal of the Indian Society of Remote Sensing》;20171231;第581-589页 * |
Also Published As
Publication number | Publication date |
---|---|
CN108389250A (zh) | 2018-08-10 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN108389250B (zh) | 基于点云数据快速生成建筑物断面图的方法 | |
CN110443836B (zh) | 一种基于平面特征的点云数据自动配准方法及装置 | |
CN109903327B (zh) | 一种稀疏点云的目标物尺寸测量方法 | |
CN108107444B (zh) | 基于激光数据的变电站异物识别方法 | |
CN109767463B (zh) | 一种三维点云自动配准方法 | |
CN105809651B (zh) | 基于边缘非相似性对比的图像显著性检测方法 | |
CN104778688A (zh) | 点云数据的配准方法及装置 | |
CN103727930A (zh) | 一种基于边缘匹配的激光测距仪与相机相对位姿标定方法 | |
CN111179321B (zh) | 一种基于模板匹配的点云配准方法 | |
CN112115774A (zh) | 结合rpa和ai的文字识别方法、装置、电子设备和存储介质 | |
CN112164145B (zh) | 基于点云数据快速提取室内三维线段结构的方法 | |
CN113192179A (zh) | 一种基于双目立体视觉的三维重建方法 | |
CN111145129A (zh) | 基于超体素的点云去噪方法 | |
CN114241018B (zh) | 一种牙齿点云配准方法、系统及可读存储介质 | |
CN110942102B (zh) | 一种概率松弛核线匹配方法及系统 | |
CN114463396B (zh) | 一种利用平面形状和拓扑图投票的点云配准方法 | |
Wang | Automatic extraction of building outline from high resolution aerial imagery | |
CN111664845B (zh) | 交通标志定位、视觉地图制作方法及装置、定位系统 | |
CN117237428B (zh) | 一种面向三维点云的数据配准方法、装置及介质 | |
CN116763295B (zh) | 牲畜体尺测量方法、电子设备及存储介质 | |
CN111127667B (zh) | 基于区域曲率二进制描述符的点云初始配准方法 | |
CN116721228B (zh) | 一种基于低密度点云的建筑物高程提取方法及系统 | |
CN113409332A (zh) | 一种基于三维点云的建筑物平面分割方法 | |
CN111091614B (zh) | 一种变电站三维实景模型建立方法及系统 | |
CN117115390A (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 | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20200522 Termination date: 20210308 |
|
CF01 | Termination of patent right due to non-payment of annual fee |