CN102855663A - 一种根据激光雷达栅格点云构建csg模型的方法 - Google Patents

一种根据激光雷达栅格点云构建csg模型的方法 Download PDF

Info

Publication number
CN102855663A
CN102855663A CN201210138143XA CN201210138143A CN102855663A CN 102855663 A CN102855663 A CN 102855663A CN 201210138143X A CN201210138143X A CN 201210138143XA CN 201210138143 A CN201210138143 A CN 201210138143A CN 102855663 A CN102855663 A CN 102855663A
Authority
CN
China
Prior art keywords
point
cylinder
plane
cloud data
sphere
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.)
Granted
Application number
CN201210138143XA
Other languages
English (en)
Other versions
CN102855663B (zh
Inventor
王晏民
郭明
罗勇
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Beijing University of Civil Engineering and Architecture
Original Assignee
Beijing University of Civil Engineering and Architecture
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Beijing University of Civil Engineering and Architecture filed Critical Beijing University of Civil Engineering and Architecture
Priority to CN201210138143.XA priority Critical patent/CN102855663B/zh
Publication of CN102855663A publication Critical patent/CN102855663A/zh
Application granted granted Critical
Publication of CN102855663B publication Critical patent/CN102855663B/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

  • Image Analysis (AREA)

Abstract

本发明公开了一种根据激光雷达栅格点云构建CSG模型的方法,包括如下步骤:读取栅格点云数据,利用栅格点云数据构建空间索引结构树;按照随机抽样一致性算法在上述空间索引结构树的叶子节点数据中提取平面、球面、以及柱面的特征点云数据;并利用上述提取的特征点云数据根据最小二乘法分别拟合平面、球面、以及柱面的规则几何模型,并将同属一个规则几何模型的模型进行合并;根据平面、球面、柱面的规则几何模型,分别构造对应的长方体、球体、以及圆柱体三维模型;将长方体、球体、以及圆柱体的规则几何模型转换为三维空间体素模型,再将所述三维空间体素模型三角化为三角网模型,之后利用该三角网模型结合多面布尔运算算法生成CSG模型。

Description

一种根据激光雷达栅格点云构建CSG模型的方法
技术领域
本发明涉及一种构建CSG(constructive solid geometry,结构化实体)模型的方法,尤其涉及一种根据激光雷达栅格点云自动构建CSG模型的方法。
背景技术
激光雷达(LIDAR)是激光技术与现代光电探测技术结合的先进探测方式,用于三维空间信息获取的激光雷达测量技术为空间信息的数字化发展提供了高效的技术支撑。地面激光雷达扫描数据处理是近年来测绘领域最热门的研究方向之一。其中,地面激光雷达数据配准与建模是数据处理的核心内容。激光雷达扫描的最原始数据存储的是扫描步进角和测站点到目标点的距离,扫描得到的三维点与点之间具有栅格性,相邻点间的拓扑关系存在。现有扫描仪厂商使用的配套软件导出通用数据格式后往往转化成散乱点云,数据处理时在散乱点云的基础上通过计算 k-邻域或者建立索引的方式重建点与点或点云块之间的拓扑关系,这种方式显然没有充分利用激光雷达获取原始数据的拓扑信息。
实际工程项目中需要对激光雷达获取的多站散乱点云数据进行快速配准,配准前要从每站散乱点云数据进行点、线、面等特征的手动提取,然后根据人工提取到的这些特征对象进行两站两站点云之间的配准工作。即使是激光雷达扫描的单站原始数据,它的数据量也相对较大,人工提取配准需要的所有特征也是比较困难的,此时,通过计算k-领域或建立索引重建点云拓扑关系的方式在特征提取的效率上就不占据优势了。另外,在生成CSG模型和深度图像模型的过程中也需要首先获取规则几何模型特征,如平面、球面、圆柱面/体或者立方体等,生成CSG模型首先要提取组成其实体的各种体素,生成深度图像的前提是要首先计算各类深度图像的基准面。
不同的三维模型在描述不同实体时具有自己的独特的优势。不规则三角网模型具有表达细节详细的优势,对于精度要求不高的模型,彩色点云模型可以发挥出其表达真实的特点。但是这些三维模型的数据量都比较大,不易进行数据高效可视化操作。在三维场景中有很多不必要精细表达的规则物体,结构化实体模型(CSG)就可以表达这些物体,它数据量小,易操作,且不损失精度。所以,在数字城市建设、较为规则的大型建筑物建模等工程中要构建CSG模型,它生成的速度和精度便成为一个很有价值的研究方向。
发明内容
本发明针对现有技术的弊端,提供一种根据激光雷达栅格点云构建CSG模型的方法。
本发明所述的根据激光雷达栅格点云构建CSG模型的方法,包括如下步骤:
步骤一,读取地面激光雷达获取的栅格点云数据,同时利用所述栅格点云数据构建空间索引结构树;
步骤二,按照随机抽样一致性算法在上述空间索引结构树的叶子节点数据中提取平面特征点云数据、球面特征点云数据、以及柱面特征点云数据;并利用上述提取的特征点云数据根据最小二乘法分别拟合平面、球面、以及柱面的规则几何模型;
步骤三,根据上述平面、球面、以及柱面的规则几何模型,分别构造对应的长方体、球体、以及圆柱体三维模型;
步骤四,将长方体、球体、以及圆柱体的规则几何模型转换为三维空间体素模型,再将所述三维空间体素模型三角化为三角网模型,之后利用该三角网模型结合多面布尔运算算法生成CSG模型。
本发明所述的根据激光雷达栅格点云构建CSG模型的方法中,所述空间索引结构树的叶子节点的数据包括节点的三维最小外包盒信息、二维行列数外包矩形信息、节点名称标识信息、节点存储点坐标的标识信息、三维坐标偏移参数信息。
本发明所述的根据激光雷达栅格点云构建CSG模型的方法的步骤一中,依次读取栅格点云数据,通过读取栅格点云数据中点的标识信息或者直接赋予标识信息,将所读取的点标识依次插入空间索引树的叶子节点;
计算每个叶子节点的最小外包矩形体,并利用均匀采样的方法以自下而上的方式填充空间索引树的非叶子节点数据,直至根节点,从而构建出多细节层次多分辨率空间索引结构树。
本发明所述的根据激光雷达栅格点云构建CSG模型的方法中,在将所读取的点依次插入空间索引树的叶子节点时,判断空间索引树的叶子节点中点的数量,若点的数量大于设定阈值,则进行节点的分裂,直至所有的点都插入到空间索引树的叶子节点中。
本发明所述的根据激光雷达栅格点云构建CSG模型的方法的步骤二中,提取球面特征点云数据的步骤包括:
随机选取一点为中心点,以由该中心点向八个方向放射的星形结构上的预定数量的点为种子点;
将所选取的种子点拟合成球面,并计算各选定的小块点云中点与该球面的距离;
选取距离小于设定阈值的点为符合条件点,当所述符合条件点的数量达到预定值时,以这些符合条件点拟合成初始球面;
移动小块点云,计算各小块点云中点与该拟合的精确球面的距离,选取距离小于设定球面距离阈值的点为球面特征点,不断迭代拟合精确球面,直到获得全部球面特征点云数据。
本发明所述的根据激光雷达栅格点云构建CSG模型的方法的步骤二中,提取平面特征点云数据的步骤包括:
以随机选取的一点为基准点,以由该基准点自下而上、自左而右选取预定数量的点为种子点,以所选取的种子点拟合成平面;
对种子点进行数量扩展,并计算各扩展种子点与所述拟合的平面之间的距离,选取距离小于设定阈值的扩展种子点为新的种子点;
将此新的种子点加入前述种子点中并再次拟合成新的平面;
重复上述步骤直至无再进行扩展的种子点,所得新的平面上的点的数据即为平面特征点云数据。
本发明所述的根据激光雷达栅格点云构建CSG模型的方法中,计算所得的多个平面的法向,以第一个平面的法向为基准,若第二平面的法向与该第一平面的法向夹角小于设定的阈值,且两个平面间两个方向上的距离小于给定阈值,则此第二平面与第一平面为同一平面;
重复上述步骤,直至所有平面判断完毕;
将判断为同一平面的各平面进行合并,从而得到各平面所对应的平面特征点云数据。
本发明所述的根据激光雷达栅格点云构建CSG模型的方法的步骤二中,提取柱面特征点云数据的步骤包括:
以随机选取的一点为基准点,以由该基准点自下而上、自左而右选取预定数量的点为种子点,以所选取的种子点拟合成柱面;
对种子点进行数量扩展,并计算各扩展种子点与所述拟合的柱面之间的距离,选取距离小于设定阈值的扩展种子点为新的种子点;
将此新的种子点加入前述种子点中并迭代拟合成新的柱面;
重复上述步骤直至无再进行扩展的种子点,所得新的柱面上的点的数据即为柱面特征点云数据。
本发明所述的根据激光雷达栅格点云构建CSG模型的方法中,计算所得的多个柱面的轴向,以第一个柱面的轴向为基准,若第二个柱面的轴向与该第一个柱面的轴向夹角小于设定的阈值,且柱面间轴向距离小于给定阈值,则此第二个柱面与第一个柱面为同一柱面;
重复上述步骤,直至所有柱面判断完毕;
将判断为同一柱面的各柱面进行合并,从而得到各柱面所对应的柱面特征点云数据。
本发明所述的根据激光雷达栅格点云构建CSG模型的方法的步骤二中,若两个或多个叶子节点中的平面的法向夹角处于预定的阈值范围内、以及所述平面的边界点集连通,则将此两个或多个叶子节点中的平面特征点云数据进行合并;
若两个或多个叶子节点中的柱面的法向夹角处于预定的阈值范围内、以及柱面的边界点集连通,则将此两个或多个叶子节点中的柱面特征点云数据进行合并;
若两个或多个叶子节点中的球面的球心坐标的差值及球体半径的差值语出预设的阈值范围内,则经此两个或多个叶子节点中的球面特征点云数据进行合并。
本发明所述的根据激光雷达栅格点云构建CSG模型的方法的步骤三中,依据球面特征点云数据而确定球心坐标信息及球体半径信息,再依据球心三维坐标信息和球体半径信息构造球体三维模型。
本发明所述的根据激光雷达栅格点云构建CSG模型的方法的步骤三中,构造长方体三维模型的步骤包括:
提取预设几个平面的点集,计算各自的最小外包矩形体属性;
确定各个点集对应的拟合平面处于相互垂直状态,选取相邻的二到三个平面点集作为构成长方体所属的点集;
计算上述确定的长方体点集的有向包围盒属性(OBB),将该三维有向包围盒对象作为长方体三维模型。
本发明所述的根据激光雷达栅格点云构建CSG模型的方法中,若各个点集对应的拟合平面的法向夹角为90±5度,则各个点集对应的拟合平面视为处于相互垂直状态。
本发明所述的根据激光雷达栅格点云构建CSG模型的方法的步骤三中,以柱面特征点云数据确定圆柱体的半径及中心轴线,以柱面特征点云数据在中心轴线上的投影点的最大的距离差为圆柱体的长度而构造圆柱体三维模型。
本发明所述的根据激光雷达栅格点云构建CSG模型的方法中,充分利用原始数据的拓扑信息进行特征提取,在过程上避免了现有技术中耗时耗力的拓扑信息重建这一步骤,在效率上有很大的提高,且能够自动提取常用的几种规则几何模型。
附图说明
图1为本发明所述根据激光雷达栅格点云构建CSG模型的方法的流程示意图。
具体实施方式
下面结合附图对本发明做进一步的详细说明,以令本领域技术人员参照说明书文字能够据以实施。
如图1所示,本发明所述的根据激光雷达栅格点云构建CSG模型的方法,包括如下步骤:
步骤101,读取地面激光雷达获取的栅格点云数据,同时利用所述栅格点云数据构建空间索引结构树。
本发明中,所述空间索引结构树的叶子节点的数据包括节点的三维最小外包盒信息、二维行列数外包矩形信息、节点名称标识信息、节点存储点坐标的标识信息、三维坐标偏移参数信息。
本步骤中,依次读取栅格点云数据,通过读取栅格点云数据中点的标识信息或者直接赋予标识信息,将所读取的点标识依次插入空间索引树的叶子节点。因为文件中点的排序就是按时间乎序排列的,所以可依次读取点的属性信息,读取顺序即为该点的ID号并置于标志位中,ID号的首位为0。在将所读取的点依次插入空间索引树的叶子节点时,判断空间索引树的叶子节点中点的数量,若点的数量大于设定阈值,则进行节点的分裂,直至所有的点都插入到空间索引树的叶子节点中。例如,设定阈值为20000点,则当前叶子节点中点的数量大于20000时,就需要进行节点的分裂,以保证在每个叶子节点中都有不超过阈值范围的点数;所述剩余的点插入到分裂后的节点中,直至所有的点都插入到叶子节点中。
计算每个叶子节点的最小外包矩形体,并利用均匀采样的方法以自下而上的方式填充空间索引树的非叶子节点数据,直至根节点,从而构建出多细节层次多分辨率空间索引结构树。也就是说该分块后的树形结构的非叶子节点也存储多分辨率的点云数据,但不存储真实三维坐标数据,只是存储坐标数据的索引ID。这样一个以四叉树为基础的空间索引结构树(称为QMBB树),在数据读取的过程中一次性快速实时生成,栅格点云在逻辑上被分割成多分辨率细节层次数据块。基于四叉树分块索引的设计思路充分考虑到扫描的原始单站点云数据的栅格性,依扫描点云的行数和列数先以自上而下的方式构建基于四叉树的点云最小外包盒的空间索引,同时对栅格点云数据进行分割,数据分割的实质是每次按照地面激光雷达实际扫描的横向步进角度和纵向步进角度分别进行二分,完成球面或扇形分割。
步骤102,按照随机抽样一致性算法在上述空间索引结构树的叶子节点数据中提取平面特征点云数据、球面特征点云数据、以及柱面特征点云数据;并利用上述提取的特征点云数据根据最小二乘法分别拟合平面、球面、以及柱面的规则几何模型。
随机抽样一致性(RANSAC)算法的样本点是随机地选取数据点样本中所需要的种子点,这显然没有利用单站扫描数据中点之间已有的拓扑关系。如果利用它固有的点之间的拓扑关系的话可以很方便地选取出有利于进行下一步运算的种子点。根据要拟合的不同二次曲面选取不同拓扑关系的种子点。利用点之间的拓扑关系选点的另一个好处就是在选取下一组种子点的时候可以方便地移动选取点的位置而不会改变已经定好的种子点的拓扑关系。
本步骤中,提取球面特征点云数据的步骤包括:
随机选取一点为中心点,以由该中心点向八个方向放射的星形结构上的预定数量的点为种子点。由于球的特有性质,无论从球面的哪个方向扫描,得到的点云数量及形状都基本相同,把球投影到各个方向也都是圆。根据球的这些特有性质及栅格点云的分布情况要把种子点选为对称的形状,以一点为中心向八个方向放射的星形结构能较好地拟合球面,减小拟合球面所导致的误差。
将所选取的种子点拟合成球面,并计算各选定的小块点云中点与该球面的距离;选取距离小于设定阈值的点为符合条件点,当所述符合条件点的数量达到预定值时,以这些符合条件点拟合成初始球面。在选择好种子点后,将它们拟合成所需要的二次曲面。然后计算分好的每小块点云中点与该二次曲面的距离,当距离小于某一阈值时,可以认为该点是二次曲面上的点。当这些符合条件的点达到一定数量,便可以认为该符合条件的点集合能够拟合成我们需要的二次曲面,再用这些符合条件的点拟合成精确的二次曲面。距离阈值和符合条件点的个数(以下称合格点)是需要确定的两个自适应参数,参数选得理想与否直接关系到球面提取的正确性和精度。合格点是首先需要确定的。根据选择的种子点可以确定所选种子点附近的点密度a,又已知在点云中的球半径R,则球的投影为面积πR2,球面上合格点通过公式πR2÷a的值来判定。
距离阈值的选择也必须适当。太大的阈值会导致误差点的选入,影响球面提取的精度,失去了该方法剔除误差点的优势。由于每种三维激光扫描仪都有一定的点位精度,距离阈值太小会导致球面上的三维点无法全部提取,无论怎样选择种子点都不能提取合格点,不能从点云中提取出球面。
移动小块点云,计算各小块点云中点与该拟合的精确球面的距离,选取距离小于设定球面距离阈值的点为球面特征点,不断迭代拟合精确球面,直到获得全部球面特征点云数据。
球面方程
F ( x , y , z ) = ( x - x 0 ) 2 + ( y - y 0 ) 2 + ( z - z 0 ) 2 - R = 0 - - - ( 3 - 1 )
其中,x0,y0,z0,R为球面参数,o=(x0,y0,z0)表示球面中心,R表示球面半径。球面方程的一般形式
F ( x , y , z ) = x 2 + y 2 + z 2 + c 7 x + c 8 y + c 9 z + c 10 = 0 x 0 = - c 7 2 y 0 = - c 8 2 z 0 = - c 9 2 R = c 7 2 + c 8 2 + c 9 2 - 4 c 10 2 - - - ( 3 - 2 )
用种子点拟合球面,由球面方程(3-1)通过参数变换得到球面方程(3-2)球面拟合可以采用线性最小二乘法实现。
由于方程(3-4)中包含常数乍xi 2+yi 2+zi 2,故可以采用奇异值分解法求解。球面拟合方程组Ax=b,其中
A = x 1 y 1 z 1 1 x 2 y 2 z 2 1 . . . . . . . . . . . . x n y n z n 1 , x = c 7 c 8 c 9 c 10 , b = x 1 2 + y 1 2 + z 1 2 x 2 2 + y 2 2 + z 2 2 . . . x n 2 + y n 2 + z n 2 - - - ( 3 - 3 )
在求得最小二乘解x后,通过坐标转换公式(3-2)就可以得到球面参数:球心坐标(x0,y0,z0)和半径R。
根据公式 Dis = ( x - x 0 ) 2 + ( y - y 0 ) 2 + ( z - z 0 ) 2 - R - - - ( 3 - 4 )
得到样本点到拟合球面的距离。如果Dis小于阈值则判定该点为球面上的点,反之则留在样本点集中等待其他拟合对象的判定。
本步骤中,提取平面特征点云数据的步骤包括:
以随机选取的一点为基准点,以由该基准点自下而上、自左而右选取预定数量的点为种子点,以所选取的种子点拟合成平面。
平面的特征和球面特征完全不同,种子点的选取也不能相同。由于扫描仪具有点位误差,即使是一个平面上的点也会有偏移这个平面的可能,所以种子点个数不能太少,且根据点密度确定,否则不能很好地控制平面参数。栅格点云具有ID从下到上、从左向右递增的特点,根据前述所选择的基准点进行向上和向右扩展。
用最初的种子点拟合平面P,拟合平面г采用最小二乘拟合。设要拟合的平面方程为a0+a1x+a2y=-z,由点(x,y,z)组成的矛盾方程组为
a 0 + a 1 x 0 + a 2 y 0 = - z 0 a 0 + a 1 x 1 + a 2 y 1 = - z 1 . . . a 0 + a 1 x n + a 2 y n = - z n - - - ( 3 - 5 )
根据最小二乘法原理,要使给定点的偏差的平方和Q最小。
Q = Σ i = 0 n ( a 0 + a 1 x i + a 2 y i + z i ) 2 - - - ( 3 - 8 )
要使Q值最小,则当Q对a0,a1,a2的偏导数同时为零时,它将出现最小值。则:
∂ Q ∂ a 0 = Σ i = 0 n 2 ( a 0 + a 1 x i + a 2 y i + z i ) = 0 ∂ Q ∂ a 1 = Σ i = 0 n 2 ( a 0 + a 1 x i + a 2 y i + z i ) x i = 0 ∂ Q ∂ a 2 = Σ i = 0 n 2 ( a 0 + a 1 x i + a 2 + y i + z i ) y i = 0 - - - ( 3 - 9 )
用矩阵形式表示可得
1 1 . . . 1 x 0 x 1 . . . x n y 0 y 1 . . . y n ( 1 x 0 y 0 1 x 1 y 1 . . . . . . . . . 1 x n y n a 0 a 1 a 2 - - z 0 - z 1 . . . - z 2 ) = 0 0 . . . 0 - - - ( 3 - 10 )
M = 1 x 0 y 0 1 x 1 y 1 . . . . . . . . . 1 x n y n , A = a 0 a 1 a 2 , = - z 0 - z 1 . . . - z 2 , 则MT(MA-Z)=0,从而
A=(MTM)-1MTZ,代入点云的坐标x,y,z求得方程系数a0,a1,a2
对种子点进行数量扩展,并计算各扩展种子点与所述拟合的平面之间的距离,选取距离小于设定阈值的扩展种子点为新的种子点。
点到平面的距离可以用方程描述为
Figure BDA00001605102500104
其中A=a1,B=a2,C=-1,D=a0。判定每个扩展点的Dis与阈值的大小,Dis小于阈值的判定为平面上的点。
将此新的种子点加入前述种子点中并再次拟合成新的平面;重复上述步骤直至无再进行扩展的种子点,所得新的平面上的点的数据即为平面特征点云数据。直到加入种子点的个数为0时,第一个平面提取结束。重新选取第二个平面的种子点进行上述循环判定,最后提取出样本点内的所有平面。
由于现实世界中的平面不是绝对意义上的平面,而上述中拟合的平面都是绝对平面,这样就难免出现实际物体上同一平面的点云被划分为多个平面点集的情况。为避免这种误判断,本实施例采用平面法向一致性和边缘相邻的判断方法,即计算所得的多个平面的法向,以第一个平面的法向为基准,若第二平面的法向与该第一平面的法向夹角小于设定的阈值,且两个平面间两个方向上的距离小于给定阈值,则此第二平面与第一平面为同一平面;重复上述步骤,直至所有平面判断完毕;将判断为同一平面的各平面进行合并,从而得到各平面所对应的平面特征点云数据。
具体如下:
根据平面方程Ax+By+Cz+D=0,求得所有平面的法向(A,B,C)。以第一个平面为基准判断其它平面与第一个平面的法向夹角,如果小于预设阈值则判定为同一平面。剩下的平面也按此规则判定,最后把所有平面都进行合并,得到误差不超过预设阈值的平面点集,完成所有平面的提取。
以上提取平面的方法虽然有较高的精度和抗噪能力,但是该提取算法是串行算法,当样本点较多时该方法的效率会大大降低。考虑到以上算法的缺点,将法向量约束加入平面的提取算法中。
首先,要求得栅格点云中每个点的法向量,点的法向量是由该点和邻域点拟合平面的法向量确定的,并把该拟合平面的法向量当作该点的法向量。栅格点云不同于散乱点云,每个点之间都存在拓扑关系,计算法向量的时候不必再费较长时间计算邻域,直接根据点之间的拓扑关系就可以得到邻域点。
法向量一致的可以判定为在同一平面上的点,现实世界中没有绝对意义上的平面,所以在实际物体同一平面上的点就会有一定的偏差。将这个偏差设定一个阈值,如果两点之间的法向偏差小于阈值则判定两点在一个平面上。进一步的,要把法向一致但不是同一平面上的点分开聚类还需加入其他条件。如果两个平面点集之间的距离(平面之间的垂直距离)大于一定距离,则分别把它们划分成两个平面。
对于单视栅格点云模型数据,采用在QMBB树索引中的所有叶子节点数据分别进行提取的方式进行。通过遍历QMBB树所有叶子节点中存储的分割数据,将分割点云数据中满足特定条件的各个平面全部提取出来,最后利用所有提取出来的平面进行点云部分和平面部分的合并,合并算法与单块分割点云数据的多平面提取相似,在此不再赘述。
本步骤中,提取柱面特征点云数据的步骤包括:
以随机选取的一点为基准点,以由该基准点自下而上、自左而右选取预定数量的点为种子点,以所选取的种子点拟合成柱面。
对种子点进行数量扩展,并计算各扩展种子点与所述拟合的柱面之间的距离,选取距离小于设定阈值的扩展种子点为新的种子点;
将此新的种子点加入前述种子点中并迭代拟合成新的柱面;
重复上述步骤直至无再进行扩展的种子点,所得新的柱面上的点的数据即为柱面特征点云数据。
本实施例中,用种子点拟合圆柱面,因为圆柱面参数有个未知数,所以种子点的个数必须不少于7个。根据圆柱的特性,圆柱面上的点到某一直线的距离等于其半径Dis,
R = [ n ( y - y 0 ) - m ( z - z 0 ) ] 2 + [ l ( z - z 0 ) - n ( x - x 0 ) ] 2 + [ m ( x - x 0 ) - l ( y - y 0 ) ] 2 l 2 + m 2 + n 2 - - - ( 3 - 11 )
其中,(x0,y0,z0)为圆柱轴线上一点,(l,m,n)为圆柱轴线的方向向量。再根据非线性方程组的最小二乘解得到圆柱面的参数。
用种子点拟合柱面后,计算样本点与已拟合圆柱面的距离Dis,如果Dis小于预设阈值则判定为圆柱面上的点,否则将该点返回样本点等待再次判定。
平面、圆柱面的提取也会出现现实场景中一个圆柱面由于点位误差被分成多个圆柱的情况。此时需要把轴向夹角不超过阈值、边缘点相邻的柱面划分为一个柱面。
即计算所得的多个柱面的轴向,以第一个柱面的轴向为基准,若第二个柱面的轴向与该第一个柱面的轴向夹角小于设定的阈值,且柱面间轴向距离小于给定阈值,则此第二个柱面与第一个柱面为同一柱面;
重复上述步骤,直至所有柱面判断完毕;
将判断为同一柱面的各柱面进行合并,从而得到各柱面所对应的柱面特征点云数据。
上述样本点均来源于QMBB树空间索引的叶子节点中的三维空间点集,由于空间索引中节点划分的终止条件是叶子节点中包含的空间点数量达到某一阈值,而不是按对象划分,因此对于一个比较大的场景来说,一个物体可能被分到两个或多个叶子节点中。之前提取出的各种特征点云可能同属一个物体上,此时,我们需要将判定在同一物体上的同一种特征点并将其合并,重新拟合新的规则几何模型。
因此,本步骤二中,若两个或多个叶子节点中的平面的法向夹角处于预定的阈值范围内、以及所述平面的边界点集连通,则将此两个或多个叶子节点中的平面特征点云数据进行合并;
若两个或多个叶子节点中的柱面的法向夹角处于预定的阈值范围内、以及柱面的边界点集连通,则将此两个或多个叶子节点中的柱面特征点云数据进行合并;
若两个或多个叶子节点中的球面的球心坐标的差值及球体半径的差值语出预设的阈值范围内,则经此两个或多个叶子节点中的球面特征点云数据进行合并。
经过上述方法可以把小块栅格点云中的规则几何模型全部提取出来。有些规则几何模型可以根据其参数构成简单体素,如:球、圆柱,而立方体只能由两个或者多个平面的参数来组成简单体素。
步骤103,根据上述平面、球面、以及柱面的规则几何模型,分别构造对应的长方体、球体、以及圆柱体三维模型。
本步骤中,依据球面特征点云数据而确定球心坐标信息及球体半径信息,再依据球心坐标信息和球体半径信息构造球体三维模型。
长方体的构造要比球体的构造复杂的多,一个完整的长方体由六个面和八个顶点组成,而栅格点云数据对一个长方体的描述只能有最多三个面组成,所以从这有限的条件获取构造长方体所需的参数不像球体那么自动和直接,长方体构造的具体过程如下:
提取预设几个平面的点集,计算各自的最小外包矩形体属性。
确定各个点集对应的拟合平面处于相互垂直状态,选取相邻的二到三个平面点集作为构成长方体所属的点集。具体而言,若各个点集对应的拟合平面的法向夹角为90±5度的范围,则各个点集对应的拟合平面视为处于相互垂直状态。
计算上述确定的长方体点集的有向包围盒属性,将该三维有向包围盒对象作为长方体三维模型。
本步骤中,以柱面特征点云数据确定圆柱体的半径及中心轴线,以柱面特征点云数据在中心轴线上的投影点的最大的距离差为圆柱体的长度而构造圆柱体三维模型。
再经过选取柱面点云特征数据后,进行圆柱面的拟合。圆柱面可以理解为到某一条直线距离都等于固定值R的点的集合,R即为圆柱的半径,直线是圆柱的中心轴所在的直线,用公式表达圆柱为:
R = [ n ( y - y 0 ) - m ( z - z 0 ) ] 2 + [ l ( z - z 0 ) - n ( x - x 0 ) ] 2 + [ m ( x - x 0 ) - l ( y - y 0 ) ] 2 l 2 + m 2 + n 2 - - - ( 3 - 12 )
其中,(x0,y0,z0)为圆柱轴线上一点,(l,m,n)为圆柱中心轴线的方向向量。再根据非线性方程组的最小二乘解得到圆柱面的参数。
求出圆柱的中心轴向后构建一个直线方程:
x = x 0 + l * t y = y 0 + m * t z = z 0 + n * t - - - ( 3 - 13 )
求圆柱面上的点在中心轴线上的投影点,取最大的距离差值即为圆柱的长度。至此,为后面构建圆柱三角网所需的圆柱自身参数已经全部获取。
步骤104,将长方体、球体、以及圆柱体的规则几何模型转换为三维空间体素模型,再将所述三维空间体素模型三角化为三角网模型,之后利用该三角网模型结合多面布尔运算算法生成CSG模型。
构造实体几何(CSG)模型是一种通过各种简单体素进行布尔运算得到新的实体表达的方法。简单体素主要包括长方体、柱体、锥体、球、环或封闭的自由曲面等,其运算为几何变换或正则布尔运算,通过对简单体素的交、并、差等正则布尔运算计算新的实体。
CSG表示可以看成是一棵有序的二叉树,称为CSG树。CSG树的结构为:<CSG树>::=<体素>|
<CSG树><集合运算节点><CSG树>|
<CSG树><几何变换节点><几何变换参数>
CSG树的叶子节点是简单体素或者形体的变换参数,非叶子结点是正则的集合运算或几何变换操作。每棵子树表示其下两个节点组合及变换的正确结果,根节点表示了最终的形体。几何变换不一定为刚体变换,也可以是不同范围的剪切变换、比例变换和对称变换。CSG树无二义性的,但不具有唯一性,它的定义是有简单体素、几何变换和正则集合运算算子组成。如果体素是正则集,则还需要进行进一步的二分,只要CSG的叶子是正确的体素,正则集的性质能够保证任何CSG树都是正确的正则集。CSG表示的特点是:数据结构比较简单,数据量较小,修改比较容易,而且可以方便地转换成边界(Brep)表示;但是由于CSG表示受体素的种类和对体素操作的种类的限制,使得它表示形体的覆盖域有较大的局限性,而且对形体的局部操作(例如,倒角等等)不易实现,显示CSG表示的结果形体时需要的间也比较长。
由前述自动提取的规则几何模型生成CSG模型的一般过程如下:首先将规则几何模型向三位空间体素转换,由面模型构建体模型,再将生成的体素三角化为多面体表达,即将体模型用多个三角面片表示。然后在体素集中任意选定两个三角网表达的基本体素,选择组合两个基本体素的集合运算方式(包括并、交、差等算子)进行两两运算,生成一个同样是三角网表达的新物体,利用基本体素和每一步新创建的物体的算子组合,继续构造新物体,直至全部体素运算完成,最终构成一个以三角网模型表达的XSG模型。
前述谈及的几何模型就是简单的体素,如平面,圆球、圆柱及长方体。经过点的聚类、二次曲面及平面的拟合得到的还不是几何实体对象,利用上面得到的参数构造需要的三维几何模型。三维几何模型的构建用顶点和连接线的方式构建三角面片,最后组成几何实体。
1)球的三角化
我们从球面栅格点云中得到的参数是球心(x0,y0,z0)和球半径R,这两个参数足以构建球体模型。具体方法如下:
建立一个单位球的格网模型,球心坐标为(0,0,0),球半径为1,以经差α=10°、纬线差β=10°构建球的经线和纬线,经线和纬线的交点就是近似球面上的点(R×cosθ×cosα,R×cosθ×sinα,R×sinθ)。按照右手系将这些点按照顺序从上到下依次连接就构成了球体的规则三角网模型。
2)长方体的三角化
从栅格点云中通过平面特征提取可以获取到所有的属于平面的点集,这些点集可以拟合成标准的平面,通过相邻三个平面可以求出三个平面的交点即所要构建长方体的一个角点P,再根据这三个平面的点集计算出点集的最小外包盒,然而最小外包盒的中心也就是要构建长方体的中心,利用长方体的一个角点和中心便可以确定长方体的另外7个角点坐标,通过已知的八个角点按照右手系的原则连接成长方体的三角网模型。
3)圆柱体的三角化
从栅格点云中通过圆柱面特征提取得到属于该圆柱面的点集,将这些点集拟合成标准的圆柱面。通过拟合的圆柱面可以获取到圆柱的半径和长度,圆柱两底面的圆心也可以确定下来,以一定的经差将圆柱沿轴向进行剖分形成垂直经线,将圆柱以一定的间距进行横向剖分得到纬线,经线和纬线相交得到圆柱的经纬网格。以圆柱一端圆心(Point1)为起点按照右手系原则从上到下将经纬线的交点连接起来便构成圆柱体的三角网表达。
基本体素三角网表达是由体素的顶点集、颜色集及连接顶点的边集组成,顶点集由Vertex类型的顶点可变数组构成,颜色集由Color类型的颜色可变数组构成,边集由int类型标识动态数组类型表达,每个int类型标识顶点集中三维点的ID号,以ID号来确定构成三角形的顶点。在物理存储结构内,首先存储三角网体素中包含的所有三维点总数,然后按序号依次存储所有顶点的三维坐标值,物理存储结构的第二部分记录三角面片信息,首先记录生成的三角网总数,然后也按序号依次存储三角形的三个顶点ID,ID定义自三维点集中存储的各个ID值。
由已经建立好的球、圆柱和立方体的规则三角网模型进行布尔运算得到所要的实体模型。布尔运算的基本原理就是拓扑重构的过程。本发明中的多面体采用点-线-面的组织方式来存储数据。
在进行布尔运算之前要对多面体对象进行最小外包盒的计算,如果最小外包盒相交再计算多面体的相交,否则就不再进行后续计算,节省资源提高效率。
1、求多面体间的交点。输入两个用点-线-面表达的多面体A、B,用多面体B的每条边与A的每个面进行线面相交计算,运用前述的计算方法求出交点,并定义在三角形内的交点为有效交点,记录有效交点;同样用A的每条边与B的每个面进行线面相交计算,保存有效交点。
2、内、外点的判断。用多面体B的每个顶点与多面体A比较,判断点是否在多面体内,若在内部,依次记录顶点为内点,否则记为外点;再用多面体A的每个顶点与多面体B比较,判断点是否在多面体内,若在内部,依次记录顶点为内点,否则记为外点。分类结果置于相应的标识位中。
3、面面相交有效线的产生。A、B两实体有各自的面M1、M2,求M1的每个边界边与相应平面M2的交点,若交点既在M1的边界上,又在M2的边界内,则该交点为有效点,记录有效点产生的交线。
4、组合体的生成包括多面体的交、并、差运算。布尔运算的交运算:在多面体A的一个面上选取至少3个点,这三个点中应该包括有效交点和内点,以内点为出发点形成一个交面这个交面是按照右手系的原则把这些选出来的点连接起来构成的。对于A中每个平面重复以上步骤,分别求出其交面。同样地,对于多面体B也用上述方法求出其各个面的交面。布尔运算的并运算与交运算方法类似,不同之处在于并运算的出发点是两多面体外点。布尔运算的差运算与并运算类似区别在于形成交面的时候是从第一个多面体外点出发。
尽管本发明的实施方案已公开如上,但其并不仅仅限于说明书和实施方式中所列运用,它完全可以被适用于各种适合本发明的领域,对于熟悉本领域的人员而言,可容易地实现另外的修改,因此在不背离权利要求及等同范围所限定的一般概念下,本发明并不限于特定的细节和这里示出与描述的图例。

Claims (10)

1.一种根据激光雷达栅格点云构建CSG模型的方法,其特征在于,包括如下步骤:
步骤一,读取地面激光雷达获取的栅格点云数据,同时利用所述栅格点云数据构建空间索引结构树;
步骤二,按照随机抽样一致性算法在上述空间索引结构树的叶子节点数据中提取平面特征点云数据、球面特征点云数据、以及柱面特征点云数据;并利用上述提取的特征点云数据根据最小二乘法分别拟合平面、球面、以及柱面的规则几何模型;
步骤三,根据上述平面、球面、以及柱面的规则几何模型,分别构造对应的长方体、球体、以及圆柱体三维模型;
步骤四,将长方体、球体、以及圆柱体的规则几何模型转换为三维空间体素模型,再将所述三维空间体素模型三角化为三角网模型,之后利用该三角网模型结合多面布尔运算算法生成CSG模型。
2.如权利要求1所述的方法,其特征在于,所述空间索引结构树的叶子节点的数据包括节点的三维最小外包盒信息、二维行列数外包矩形信息、节点名称标识信息、节点存储点坐标的标识信息、三维坐标偏移参数信息。
3.如权利要求2所述的方法,其特征在于,所述步骤一中,依次读取栅格点云数据,通过读取栅格点云数据中点的标识信息或者直接赋予标识信息,将所读取的点标识依次插入空间索引树的叶子节点;在将所读取的点依次插入空间索引树的叶子节点时,判断空间索引树的叶子节点中点的数量,若点的数量大于设定阈值,则进行节点的分裂,直至所有的点都插入到空间索引树的叶子节点中;
计算每个叶子节点的最小外包矩形体,并利用均匀采样的方法以自下而上的方式填充空间索引树的非叶子节点数据,直至根节点,从而构建出多细节层次多分辨率空间索引结构树。
4.如权利要求2所述的方法,其特征在于,所述步骤二中,提取球面特征点云数据的步骤包括:
随机选取一点为中心点,以由该中心点向八个方向放射的星形结构上的预定数量的点为种子点;
将所选取的种子点拟合成球面,并计算各选定的小块点云中点与该球面的距离;
选取距离小于设定阈值的点为符合条件点,当所述符合条件点的数量达到预定值时,以这些符合条件点拟合成初始球面;
移动小块点云,计算各小块点云中点与该拟合的精确球面的距离,选取距离小于设定球面距离阈值的点为球面特征点,不断迭代拟合精确球面,直到获得全部球面特征点云数据。
5.如权利要求2所述的方法,其特征在于,所述步骤二中,提取平面特征点云数据的步骤包括:
以随机选取的一点为基准点,以由该基准点自下而上、自左而右选取预定数量的点为种子点,以所选取的种子点拟合成平面;
对种子点进行数量扩展,并计算各扩展种子点与所述拟合的平面之间的距离,选取距离小于设定阈值的扩展种子点为新的种子点;
将此新的种子点加入前述种子点中并再次拟合成新的平面;
重复上述步骤直至无再进行扩展的种子点,所得新的平面上的点的数据即为平面特征点云数据;
并执行如下步骤:计算所得的多个平面的法向,以第一个平面的法向为基准,若第二平面的法向与该第一平面的法向夹角小于设定的阈值,且两个平面间两个方向上的距离小于给定阈值,则此第二平面与第一平面为同一平面;
重复上述步骤,直至所有平面判断完毕;
将判断为同一平面的各平面进行合并,从而得到各平面所对应的平面特征点云数据。
6.如权利要求2所述的方法,其特征在于,所述步骤二中,提取柱面特征点云数据的步骤包括:
以随机选取的一点为基准点,以由该基准点自下而上、自左而右选取预定数量的点为种子点,以所选取的种子点拟合成柱面;
对种子点进行数量扩展,并计算各扩展种子点与所述拟合的柱面之间的距离,选取距离小于设定阈值的扩展种子点为新的种子点;
将此新的种子点加入前述种子点中并迭代拟合成新的柱面;
重复上述步骤直至无再进行扩展的种子点,所得新的柱面上的点的数据即为柱面特征点云数据;
并执行如下步骤:计算所得的多个柱面的轴向,以第一个柱面的轴向为基准,若第二个柱面的轴向与该第一个柱面的轴向夹角小于设定的阈值,且柱面间轴向距离小于给定阈值,则此第二个柱面与第一个柱面为同一柱面;
重复上述步骤,直至所有柱面判断完毕;
将判断为同一柱面的各柱面进行合并,从而得到各柱面所对应的柱面特征点云数据。
7.如权利要求2所述的方法,其特征在于,所述步骤二中,若两个或多个叶子节点中的平面的法向夹角处于预定的阈值范围内、以及所述平面的边界点集连通,则将此两个或多个叶子节点中的平面特征点云数据进行合并;
若两个或多个叶子节点中的柱面的法向夹角处于预定的阈值范围内、以及柱面的边界点集连通,则将此两个或多个叶子节点中的柱面特征点云数据进行合并;
若两个或多个叶子节点中的球面的球心坐标的差值及球体半径的差值语出预设的阈值范围内,则经此两个或多个叶子节点中的球面特征点云数据进行合并。
8.如权利要求2所述的方法,其特征在于,所述步骤三中,依据球面特征点云数据而确定球心坐标信息及球体半径信息,再依据球心坐标信息和球体半径信息构造球体三维模型。
9.如权利要求2所述的方法,其特征在于,所述步骤三中,构造长方体三维模型的步骤包括:
提取预设几个平面的点集,计算各自的最小外包矩形体属性;
确定各个点集对应的拟合平面处于相互垂直状态,选取相邻的二到三个平面点集作为构成长方体所属的点集;其中,若各个点集对应的拟合平面的法向夹角为90±5度,则各个点集对应的拟合平面视为处于相互垂直状态;
计算上述确定的长方体点集的有向包围盒属性,将该三维有向包围盒对象作为长方体三维模型。
10.如权利要求2所述的方法,其特征在于,所述步骤三中,以柱面特征点云数据确定圆柱体的半径及中心轴线,以柱面特征点云数据在中心轴线上的投影点的最大的距离差为圆柱体的长度而构造圆柱体三维模型。
CN201210138143.XA 2012-05-04 2012-05-04 一种根据激光雷达栅格点云构建csg模型的方法 Expired - Fee Related CN102855663B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201210138143.XA CN102855663B (zh) 2012-05-04 2012-05-04 一种根据激光雷达栅格点云构建csg模型的方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201210138143.XA CN102855663B (zh) 2012-05-04 2012-05-04 一种根据激光雷达栅格点云构建csg模型的方法

Publications (2)

Publication Number Publication Date
CN102855663A true CN102855663A (zh) 2013-01-02
CN102855663B CN102855663B (zh) 2015-04-01

Family

ID=47402222

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201210138143.XA Expired - Fee Related CN102855663B (zh) 2012-05-04 2012-05-04 一种根据激光雷达栅格点云构建csg模型的方法

Country Status (1)

Country Link
CN (1) CN102855663B (zh)

Cited By (36)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103366006A (zh) * 2013-07-19 2013-10-23 武汉睿数信息技术有限公司 一种基于三维激光点云建立地理空间索引的方法
CN104102557A (zh) * 2014-06-27 2014-10-15 武汉理工大学 一种基于聚类的云计算平台数据备份方法
CN104732544A (zh) * 2015-04-01 2015-06-24 郑州辰维科技股份有限公司 一种快速查找形状目标点的方法
CN105160344A (zh) * 2015-06-18 2015-12-16 北京大学深圳研究生院 一种三维点云的局部特征提取方法及装置
CN105513051A (zh) * 2015-11-26 2016-04-20 福州华鹰重工机械有限公司 一种点云数据处理方法和设备
CN105809738A (zh) * 2016-03-09 2016-07-27 同济大学 一种基于生态系统和csg的复杂产品模型建立及更新方法
CN105809615A (zh) * 2016-03-10 2016-07-27 广州欧科信息技术股份有限公司 点云数据影像化方法
CN105844600A (zh) * 2016-04-27 2016-08-10 北京航空航天大学 一种空间目标三维点云光顺去噪方法
CN105893658A (zh) * 2016-03-28 2016-08-24 同济大学 一种基于复杂网络的复杂产品生态网络建模方法
CN106500594A (zh) * 2016-10-31 2017-03-15 武汉大学 融合反射强度和几何特征的铁路轨道半自动检测方法
CN106600684A (zh) * 2016-11-29 2017-04-26 浙江科澜信息技术有限公司 一种倾斜模型组织构建方法
CN106875439A (zh) * 2017-02-20 2017-06-20 天津大学 基于三维点云模型的单晶硅棒外形尺寸测量方法
CN107301277A (zh) * 2017-06-05 2017-10-27 中国空间技术研究院 一种基于航天器数字化模型的产品结构树重构方法
CN107958209A (zh) * 2017-11-16 2018-04-24 深圳天眼激光科技有限公司 一种违建识别方法、系统及电子设备
CN109141266A (zh) * 2018-07-19 2019-01-04 北京建筑大学 一种钢结构测量方法及系统
CN109359614A (zh) * 2018-10-30 2019-02-19 百度在线网络技术(北京)有限公司 一种激光点云的平面识别方法、装置、设备和介质
CN109410260A (zh) * 2018-09-27 2019-03-01 先临三维科技股份有限公司 点云数据网格化方法、装置、计算机设备和存储介质
CN109461207A (zh) * 2018-11-05 2019-03-12 胡翰 一种点云数据建筑物单体化方法及装置
CN109919955A (zh) * 2019-03-11 2019-06-21 南京林业大学 地基式激光雷达点云的隧道轴线提取和分割方法
CN109947884A (zh) * 2019-03-06 2019-06-28 武汉大学 一种全球ICESat/GLAS点云的高效组织查询方法
CN110095062A (zh) * 2019-04-17 2019-08-06 北京华捷艾米科技有限公司 一种物体体积参数测量方法、装置及设备
CN110223390A (zh) * 2019-06-12 2019-09-10 辽宁工程技术大学 一种基于直线体素遍历原理的多段线嵌入tin算法
WO2020073483A1 (zh) * 2018-10-09 2020-04-16 上海莉莉丝科技股份有限公司 在编辑器中显示物体的方法、系统、设备和介质
CN111322950A (zh) * 2020-04-17 2020-06-23 易思维(杭州)科技有限公司 利用线结构光传感器定位圆柱体位置的方法及其用途
CN111350214A (zh) * 2020-03-23 2020-06-30 中交第三航务工程局有限公司江苏分公司 多波束水下钢管桩桩位测量方法
CN111445472A (zh) * 2020-03-26 2020-07-24 达闼科技成都有限公司 激光点云地面分割方法、装置、计算设备及存储介质
CN111680747A (zh) * 2020-06-08 2020-09-18 北京百度网讯科技有限公司 用于占据栅格子图的闭环检测的方法和装置
CN112712596A (zh) * 2021-03-29 2021-04-27 深圳大学 一种密集匹配点云建筑物结构化模型精细重建方法
CN112767551A (zh) * 2021-01-18 2021-05-07 北京房江湖科技有限公司 三维模型构建方法、装置以及电子设备、存储介质
CN112902844A (zh) * 2021-02-24 2021-06-04 三门核电有限公司 一种基于激光跟踪的设备圆形端口建模方法
CN112991547A (zh) * 2021-01-18 2021-06-18 北京房江湖科技有限公司 模型简化处理方法、装置以及电子设备、介质
CN113160398A (zh) * 2020-12-25 2021-07-23 中国人民解放军国防科技大学 一种快速三维栅格构建系统、方法、介质、设备、无人车
CN113412505A (zh) * 2019-03-22 2021-09-17 华为技术有限公司 用于对通过探测和测距传感器获得的点云进行有序表示和特征提取的系统和方法
CN113569856A (zh) * 2021-07-13 2021-10-29 盎锐(上海)信息科技有限公司 用于实测实量的模型语义分割方法及激光雷达
CN114782508A (zh) * 2022-06-21 2022-07-22 杭州三坛医疗科技有限公司 髋关节配准点筛选方法及装置
CN117934732A (zh) * 2024-03-22 2024-04-26 北京华航唯实机器人科技股份有限公司 点云模型的重建方法、装置及电子设备

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110110687B (zh) * 2019-05-15 2020-11-17 江南大学 基于颜色信息和三维轮廓信息的树上水果自动识别方法

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5537519A (en) * 1993-06-29 1996-07-16 Electronic Data Systems Corporation System and method for converting boundary representations to constructive solid geometry representations for three-dimensional solid object modeling
CN101976467A (zh) * 2010-09-13 2011-02-16 天津市星际空间地理信息工程有限公司 一种集成机载lidar和车载移动激光扫描技术的高精细三维城市场景构建方法

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5537519A (en) * 1993-06-29 1996-07-16 Electronic Data Systems Corporation System and method for converting boundary representations to constructive solid geometry representations for three-dimensional solid object modeling
CN101976467A (zh) * 2010-09-13 2011-02-16 天津市星际空间地理信息工程有限公司 一种集成机载lidar和车载移动激光扫描技术的高精细三维城市场景构建方法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
JYUN-YUAN CHEN ET AL.: "Point Cloud Modeling Using Algebraic Template", 《INTERNATIONAL JOURNAL OF INNOVATIVE COMPUTING, INFORMATION AND CONTROL》 *
陈传法等: "基于近景目标3维测量系统和CSG技术构建3维景观", 《测绘与空间地理信息》 *

Cited By (56)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103366006B (zh) * 2013-07-19 2016-06-22 武汉睿数信息技术有限公司 一种基于三维激光点云建立地理空间索引的方法
CN103366006A (zh) * 2013-07-19 2013-10-23 武汉睿数信息技术有限公司 一种基于三维激光点云建立地理空间索引的方法
CN104102557A (zh) * 2014-06-27 2014-10-15 武汉理工大学 一种基于聚类的云计算平台数据备份方法
CN104102557B (zh) * 2014-06-27 2017-11-10 武汉理工大学 一种基于聚类的云计算平台数据备份方法
CN104732544B (zh) * 2015-04-01 2017-07-11 郑州辰维科技股份有限公司 一种快速查找形状目标点的方法
CN104732544A (zh) * 2015-04-01 2015-06-24 郑州辰维科技股份有限公司 一种快速查找形状目标点的方法
CN105160344A (zh) * 2015-06-18 2015-12-16 北京大学深圳研究生院 一种三维点云的局部特征提取方法及装置
CN105513051A (zh) * 2015-11-26 2016-04-20 福州华鹰重工机械有限公司 一种点云数据处理方法和设备
CN105809738A (zh) * 2016-03-09 2016-07-27 同济大学 一种基于生态系统和csg的复杂产品模型建立及更新方法
CN105809615A (zh) * 2016-03-10 2016-07-27 广州欧科信息技术股份有限公司 点云数据影像化方法
CN105893658A (zh) * 2016-03-28 2016-08-24 同济大学 一种基于复杂网络的复杂产品生态网络建模方法
CN105893658B (zh) * 2016-03-28 2019-07-05 同济大学 一种基于复杂网络的复杂产品生态网络建模方法
CN105844600A (zh) * 2016-04-27 2016-08-10 北京航空航天大学 一种空间目标三维点云光顺去噪方法
CN105844600B (zh) * 2016-04-27 2018-03-16 北京航空航天大学 一种空间目标三维点云光顺去噪方法
CN106500594B (zh) * 2016-10-31 2017-10-03 武汉大学 融合反射强度和几何特征的铁路轨道半自动检测方法
CN106500594A (zh) * 2016-10-31 2017-03-15 武汉大学 融合反射强度和几何特征的铁路轨道半自动检测方法
CN106600684A (zh) * 2016-11-29 2017-04-26 浙江科澜信息技术有限公司 一种倾斜模型组织构建方法
CN106875439B (zh) * 2017-02-20 2019-06-07 天津大学 基于三维点云模型的单晶硅棒外形尺寸测量方法
CN106875439A (zh) * 2017-02-20 2017-06-20 天津大学 基于三维点云模型的单晶硅棒外形尺寸测量方法
CN107301277A (zh) * 2017-06-05 2017-10-27 中国空间技术研究院 一种基于航天器数字化模型的产品结构树重构方法
CN107301277B (zh) * 2017-06-05 2020-07-14 中国空间技术研究院 一种基于航天器数字化模型的产品结构树重构方法
CN107958209A (zh) * 2017-11-16 2018-04-24 深圳天眼激光科技有限公司 一种违建识别方法、系统及电子设备
CN109141266A (zh) * 2018-07-19 2019-01-04 北京建筑大学 一种钢结构测量方法及系统
CN109410260A (zh) * 2018-09-27 2019-03-01 先临三维科技股份有限公司 点云数据网格化方法、装置、计算机设备和存储介质
CN109410260B (zh) * 2018-09-27 2020-12-29 先临三维科技股份有限公司 点云数据网格化方法、装置、计算机设备和存储介质
WO2020073483A1 (zh) * 2018-10-09 2020-04-16 上海莉莉丝科技股份有限公司 在编辑器中显示物体的方法、系统、设备和介质
US11294541B2 (en) 2018-10-09 2022-04-05 Shanghai Lilith Technology Corporation Method, system, and device for displaying geometry in editor, and medium
CN109359614A (zh) * 2018-10-30 2019-02-19 百度在线网络技术(北京)有限公司 一种激光点云的平面识别方法、装置、设备和介质
CN109461207A (zh) * 2018-11-05 2019-03-12 胡翰 一种点云数据建筑物单体化方法及装置
CN109947884A (zh) * 2019-03-06 2019-06-28 武汉大学 一种全球ICESat/GLAS点云的高效组织查询方法
CN109947884B (zh) * 2019-03-06 2021-04-02 武汉大学 一种全球ICESat/GLAS点云的高效组织查询方法
CN109919955A (zh) * 2019-03-11 2019-06-21 南京林业大学 地基式激光雷达点云的隧道轴线提取和分割方法
CN113412505B (zh) * 2019-03-22 2022-09-09 华为技术有限公司 用于对通过探测和测距传感器获得的点云进行有序表示和特征提取的处理单元和方法
US11556745B2 (en) 2019-03-22 2023-01-17 Huawei Technologies Co., Ltd. System and method for ordered representation and feature extraction for point clouds obtained by detection and ranging sensor
CN113412505A (zh) * 2019-03-22 2021-09-17 华为技术有限公司 用于对通过探测和测距传感器获得的点云进行有序表示和特征提取的系统和方法
CN110095062A (zh) * 2019-04-17 2019-08-06 北京华捷艾米科技有限公司 一种物体体积参数测量方法、装置及设备
CN110223390A (zh) * 2019-06-12 2019-09-10 辽宁工程技术大学 一种基于直线体素遍历原理的多段线嵌入tin算法
CN110223390B (zh) * 2019-06-12 2023-02-17 辽宁工程技术大学 一种基于直线体素遍历原理的多段线嵌入tin算法
CN111350214A (zh) * 2020-03-23 2020-06-30 中交第三航务工程局有限公司江苏分公司 多波束水下钢管桩桩位测量方法
CN111350214B (zh) * 2020-03-23 2021-07-30 中交第三航务工程局有限公司江苏分公司 多波束水下钢管桩桩位测量方法
CN111445472A (zh) * 2020-03-26 2020-07-24 达闼科技成都有限公司 激光点云地面分割方法、装置、计算设备及存储介质
CN111445472B (zh) * 2020-03-26 2023-05-12 达闼机器人股份有限公司 激光点云地面分割方法、装置、计算设备及存储介质
CN111322950A (zh) * 2020-04-17 2020-06-23 易思维(杭州)科技有限公司 利用线结构光传感器定位圆柱体位置的方法及其用途
CN111322950B (zh) * 2020-04-17 2021-08-17 易思维(杭州)科技有限公司 利用线结构光传感器定位圆柱体位置的方法及其用途
CN111680747B (zh) * 2020-06-08 2023-09-01 北京百度网讯科技有限公司 用于占据栅格子图的闭环检测的方法和装置
CN111680747A (zh) * 2020-06-08 2020-09-18 北京百度网讯科技有限公司 用于占据栅格子图的闭环检测的方法和装置
CN113160398A (zh) * 2020-12-25 2021-07-23 中国人民解放军国防科技大学 一种快速三维栅格构建系统、方法、介质、设备、无人车
CN112991547B (zh) * 2021-01-18 2022-05-17 贝壳找房(北京)科技有限公司 模型简化处理方法、装置以及电子设备、介质
CN112767551A (zh) * 2021-01-18 2021-05-07 北京房江湖科技有限公司 三维模型构建方法、装置以及电子设备、存储介质
CN112991547A (zh) * 2021-01-18 2021-06-18 北京房江湖科技有限公司 模型简化处理方法、装置以及电子设备、介质
CN112902844A (zh) * 2021-02-24 2021-06-04 三门核电有限公司 一种基于激光跟踪的设备圆形端口建模方法
CN112712596A (zh) * 2021-03-29 2021-04-27 深圳大学 一种密集匹配点云建筑物结构化模型精细重建方法
CN113569856A (zh) * 2021-07-13 2021-10-29 盎锐(上海)信息科技有限公司 用于实测实量的模型语义分割方法及激光雷达
CN113569856B (zh) * 2021-07-13 2024-06-04 盎锐(杭州)信息科技有限公司 用于实测实量的模型语义分割方法及激光雷达
CN114782508A (zh) * 2022-06-21 2022-07-22 杭州三坛医疗科技有限公司 髋关节配准点筛选方法及装置
CN117934732A (zh) * 2024-03-22 2024-04-26 北京华航唯实机器人科技股份有限公司 点云模型的重建方法、装置及电子设备

Also Published As

Publication number Publication date
CN102855663B (zh) 2015-04-01

Similar Documents

Publication Publication Date Title
CN102855663B (zh) 一种根据激光雷达栅格点云构建csg模型的方法
CN111008422B (zh) 一种建筑物实景地图制作方法及系统
CN111932671A (zh) 基于密集点云数据的三维实体模型重建方法
CN103258345B (zh) 一种基于地面激光雷达三维扫描的树木枝干参数提取方法
CN103247041B (zh) 一种基于局部采样的多几何特征点云数据的分割方法
CN103500467B (zh) 基于图像的三维模型构成方法
LU102117B1 (en) Method and system for measuring mountain view visible area in city
CN102044088B (zh) 单站地面激光扫描海量散乱点云的lod模型快速构建方法
CN104778744B (zh) 基于Lidar数据的大规模三维森林可视化场景建立方法
CN103065353B (zh) 三维模型特征提取方法和系统、三维模型检索方法和系统
JP2011512606A (ja) 樹冠植生ポイントのグループに対しシミュレーションによる木の幹と枝とを与えるジオスペースモデリングシステム及び関連方法
CN111737790B (zh) 仿真城市模型的构建方法和设备
CN107657659A (zh) 基于长方体拟合扫描三维点云的曼哈顿结构建筑物自动建模方法
CN108776993A (zh) 带有孔洞的三维点云的建模方法及地下电缆工井建模方法
CN106127857A (zh) 综合数据驱动与模型驱动的机载LiDAR数据建模方法
Zhao Application of 3D CAD in landscape architecture design and optimization of hierarchical details
CN107767453A (zh) 一种基于规则约束的建筑物lidar点云重构优化方法
US20120265494A1 (en) Method of Online Building-Model Reconstruction Using Photogrammetric Mapping System
CN113916130B (zh) 一种基于最小二乘法的建筑物位置测量方法
CN104751479A (zh) 基于tin数据的建筑物提取方法和装置
CN107545602A (zh) 基于LiDAR点云的空间拓扑关系约束下的建筑物建模方法
CN102567465A (zh) 复眼视觉基础上的三维模型几何相似性搜索方法
CN106067160A (zh) 大屏幕融合投影方法
CN109472862A (zh) 一种变电站三维建模系统
CN114972672A (zh) 输电线路实景三维模型构建方法、装置、设备和存储介质

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
C14 Grant of patent or utility model
GR01 Patent grant
C56 Change in the name or address of the patentee
CB03 Change of inventor or designer information

Inventor after: Guo Ming

Inventor after: Wang Yanmin

Inventor after: Luo Yong

Inventor before: Wang Yanmin

Inventor before: Guo Ming

Inventor before: Luo Yong

COR Change of bibliographic data
CP01 Change in the name or title of a patent holder

Address after: 100044 Xicheng District Exhibition Hall Road, Beijing, No. 1

Patentee after: Beijing University of Civil Engineering and Architecture

Address before: 100044 Xicheng District Exhibition Hall Road, Beijing, No. 1

Patentee before: Beijing Engineering and Architecture School

CF01 Termination of patent right due to non-payment of annual fee
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20150401

Termination date: 20210504