CN104765915B - 三维激光扫描数据建模方法及系统 - Google Patents
三维激光扫描数据建模方法及系统 Download PDFInfo
- Publication number
- CN104765915B CN104765915B CN201510143817.9A CN201510143817A CN104765915B CN 104765915 B CN104765915 B CN 104765915B CN 201510143817 A CN201510143817 A CN 201510143817A CN 104765915 B CN104765915 B CN 104765915B
- Authority
- CN
- China
- Prior art keywords
- max
- cos
- laser scanning
- scanning data
- hemisphere face
- 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
- 238000000034 method Methods 0.000 title claims abstract description 36
- 238000013507 mapping Methods 0.000 claims abstract description 94
- ATJFFYVFTNAWJD-UHFFFAOYSA-N Tin Chemical compound [Sn] ATJFFYVFTNAWJD-UHFFFAOYSA-N 0.000 claims abstract description 22
- 239000000523 sample Substances 0.000 claims description 51
- 230000003287 optical effect Effects 0.000 claims 1
- 238000010586 diagram Methods 0.000 description 11
- 238000005516 engineering process Methods 0.000 description 3
- 230000001788 irregular Effects 0.000 description 2
- 239000002184 metal Substances 0.000 description 2
- 239000007787 solid Substances 0.000 description 2
- 238000010408 sweeping Methods 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000015572 biosynthetic process Effects 0.000 description 1
- 230000001419 dependent effect Effects 0.000 description 1
- 239000006185 dispersion Substances 0.000 description 1
- 235000013399 edible fruits Nutrition 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000012805 post-processing Methods 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
Landscapes
- Length Measuring Devices By Optical Means (AREA)
- Length Measuring Devices With Unspecified Measuring Means (AREA)
Abstract
本发明提供了一种基于半球面和平面联合投影的三维激光扫描数据建模方法及系统,所述方法包括将N个三维激光扫描数据点投影到半球面上,得到半球面映射点;将所述N个半球面映射点投影到平面上,得到N个三维激光扫描数据面映射点;将所述平面映射点转化为N个直角坐标,并进行任意轴坐标旋转,得到坐标系在水平面上的平面映射点;对水平面上的平面映射点进行三角剖分,得到平面映射点之间的不规则三角网拓扑关系;将所述平面映射点之间的不规则三角网拓扑关系对应地赋值给三维激光扫描数据点,得到三维激光扫描数据点之间的不规则三角网拓扑关系,完成对三维激光扫描数据点的建模。本发明能够实现三维激光扫描数据的快速准确建模。
Description
技术领域
本发明涉及激光扫描领域,具体涉及一种基于半球面和平面联合投影的三维激光扫描数据建模方法及系统。
背景技术
三维激光扫描仪是一种集成了多种高新技术的新型三维坐标测量仪器,采用非接触式高速激光测量方式,以点云形式获取被测对象表面的阵列式几何图形的三维数据,其激光测距的基本原理已发展成熟,数据获取的准确性和有效性也得以保障,然而扫描的结果只是待测实体表面的点云数据,在应用中需要更具有实际意义的三维实体表面模型,因此,对扫描数据的后处理一直面临着困难。
解决点云数据建模问题的方法可分为三类:基于隐式曲面法、区域生长法和基于Delaunay三角剖分法。Hoppe等最早提出的零集法属于基于隐式曲面法,该类方法不能在点云中插值点从而使其应用受到限制;Bernardini等提出的Ball-Pivoting法和CochenSteiner提出的贪婪算法属于区域生长法,该类方法构建的网格质量过分依赖于种子面片的选取和自定义参数的设置,采样不均匀的情况下还需要后续补孔;具有代表性的基于Delaunay三角剖分的算法有Guo等提出的α-shape法、Dey等提出的Cocone算法和Amenta等提出的Crust算法及Power Crust算法等,该类算法由于使用了中间数据结构,计算量大,时效性比较差。
针对上述现象,亟需一种快速准确的建模方法,解决三维激光扫描数据的建模问题。
发明内容
针对现有技术的不足,本发明提供一种基于半球面和平面联合投影的三维激光扫描数据建模方法,以实现三维激光扫描数据的快速准确建模。
为实现上述目的,本发明通过以下技术方案予以实现:
第一方面,本发明提供一种基于半球面和平面联合投影的三维激光扫描数据建模方法,包括:
S1.根据扫描头的坐标PO(pox,poy,poz)、扫描头的初始方位(θ,α)及N个三维激光扫描数据点的极坐标Poi(θi,αi,di),得到扫描数据点到所述扫描头的最长距离dmax;
S2.通过半球面投影法则将N个三维激光扫描数据点投影到半球面上,得到N个三维激光扫描数据半球面映射点Psi(θi/2,αi,dmax);
S3.通过平面投影法则将所述N个三维激光扫描数据半球面映射点投影到平面上,得到N个三维激光扫描数据平面映射点Ppi(θi/2,αi,dmax/cos(θi/2));
S4.将所述N个三维激光扫描数据平面映射点转化为N个直角坐标Pri(prix,priy,priz),并进行任意轴坐标旋转,得到坐标系在水平面上的平面映射点Pti(ptix,ptiy,ptiz);
S5.对水平面上的平面映射点Pti进行Delaunay三角剖分,得到Pti之间的不规则三角网拓扑关系;
S6.将所述Pti之间的不规则三角网拓扑关系对应地赋值给Poi,从而得到Poi之间的不规则三角网拓扑关系,完成对三维激光扫描数据点的建模。
其中,所述步骤S1根据扫描头的坐标PO(pox,poy,poz)、扫描头的初始方位(θ,α)及N个三维激光扫描数据点的极坐标Poi(θi,αi,di),得到扫描数据点到所述扫描头的最长距离dmax包括:
根据扫描头的坐标PO(pox,poy,poz)及扫描头的初始方位(θ,α)建立激光扫描仪坐标系,N个三维激光扫描数据点的极坐标Poi(θi,αi,di)为以激光扫描仪坐标系为基础的表示方式,将其转化为以大地坐标系为基础的表示方式Pdi(θi+θ,αi+α,di),比较N个三维激光扫描数据点Poi的di值,得到三维激光扫描数据点到扫描头的最长距离dmax。
其中,所述步骤S2通过半球面投影法则将N个三维激光扫描数据点投影到半球面上,得到N个三维激光扫描数据半球面映射点Psi(θi/2,αi,dmax)包括:
根据所述半球面投影法则Poi(θi,αi,di)→Psi(θi/2,αi,dmax)将极坐标(θi,αi,di)中的di值变换为dmax,dmax为一个固定值,数据点(θi,αi,dmax)落在一个球面上;且将倾角θi折半,此时数据点(θi,αi,dmax)在球面上收缩为一个半球面,得到半球面映射数据点Psi(θi/2,αi,dmax)。
其中,所述步骤S3通过平面投影法则将所述N个三维激光扫描数据半球面映射点投影到平面上,得到N个三维激光扫描数据平面映射点Ppi(θi/2,αi,dmax/cos(θi/2))包括:
根据所述平面投影法则Psi(θi/2,αi,dmax)→Ppi(θi/2,αi,dmax/cos(θi/2)),将半球面映射数据点(θi/2,αi,dmax)中的dmax值变换为dmax/cos(θi/2),此时数据点(θi/2,αi,dmax/cos(θi/2))便在一个平面上,从而得到平面映射数据点Ppi(θi/2,αi,dmax/cos(θi/2))。
其中,所述步骤S4将所述N个三维激光扫描数据平面映射点转化为N个直角坐标Pri(prix,priy,priz),并进行任意轴坐标旋转,得到坐标系在水平面上的平面映射点Pti(ptix,ptiy,ptiz)包括:
将以激光扫描仪坐标系为基础的平面映射点Ppi(θi/2,αi,dmax/cos(θi/2))的表示方式转化为以大地坐标系为基础的表示方式Pqi(θi/2+θ,αi+α,dmax/cos(θi/2));
按照下面公式将所述N个三维激光扫描数据平面映射点转化为N个直角坐标Pri(prix,priy,priz):
按照下面步骤进行任意轴坐标旋转:
A1.计算得到Ppi所在的平面法向量为(cosθcosα,cosθsinα,sinθ),Pti所在的水平面法向量为(0,0,1);
A2.计算得到旋转轴A为(cosθsinα,-cosθcosα,0);
A3.计算得到旋转角θr为arccos(sinθ);
A4.进行任意轴旋转:
其中,Ax表示cosθsinα,Ay表示-cosθcosα,Az表示0,C表示cosθr,S表示sinθr。
第二方面,本发明提供一种基于半球面和平面联合投影的三维激光扫描数据建模系统,包括预处理单元、半球面投影单元、平面投影单元、任意轴旋转单元、平面三角拓扑关系获取单元和三维激光数据建模单元;
其中,所述预处理单元用于根据扫描头的坐标PO(pox,poy,poz)、扫描头的初始方位(θ,α)及N个三维激光扫描数据点的极坐标Poi(θi,αi,di),得到扫描数据点到所述扫描头的最长距离dmax;
所述半球面投影单元用于通过半球面投影法则将N个三维激光扫描数据点投影到半球面上,得到N个三维激光扫描数据半球面映射点Psi(θi/2,αi,dmax);
所述平面投影单元用于通过平面投影法则将所述N个三维激光扫描数据半球面映射点投影到平面上,得到N个三维激光扫描数据平面映射点Ppi(θi/2,αi,dmax/cos(θi/2));
所述任意轴旋转单元用于将所述N个三维激光扫描数据平面映射点转化为N个直角坐标Pri(prix,priy,priz),并进行任意轴坐标旋转,得到坐标系在水平面上的平面映射点Pti(ptix,ptiy,ptiz);
所述平面三角拓扑关系获取单元用于对水平面上的平面映射点Pti进行Delaunay三角剖分,得到Pti之间的不规则三角网拓扑关系;
所述三维激光数据建模单元用于将所述Pti之间的不规则三角网拓扑关系对应地赋值给Poi,从而得到Poi之间的不规则三角网拓扑关系,完成对三维激光扫描数据点的建模。
其中,所述预处理单元具体用于:
根据扫描头的坐标PO(pox,poy,poz)及扫描头的初始方位(θ,α)建立激光扫描仪坐标系,N个三维激光扫描数据点的极坐标Poi(θi,αi,di)为以激光扫描仪坐标系为基础的表示方式,将其转化为以大地坐标系为基础的表示方式Pdi(θi+θ,αi+α,di),比较N个三维激光扫描数据点Poi的di值,得到三维激光扫描数据点到扫描头的最长距离dmax。
其中,所述半球面投影单元具体用于:
根据所述半球面投影法则Poi(θi,αi,di)→Psi(θi/2,αi,dmax),将极坐标(θi,αi,di)中的di值变换为dmax,dmax为一个固定值,数据点(θi,αi,dmax)在一个球面上;且将倾角θi折半,此时数据点(θi,αi,dmax)在球面上收缩为一个半球面,得到半球面映射数据点Psi(θi/2,αi,dmax)。
其中,所述平面投影单元具体用于:
根据所述平面投影法则Psi(θi/2,αi,dmax)→Ppi(θi/2,αi,dmax/cos(θi/2)),将半球面映射数据点(θi/2,αi,dmax)中的dmax值变换为dmax/cos(θi/2),此时数据点(θi/2,αi,dmax/cos(θi/2))便在一个平面上,从而得到平面映射数据点Ppi(θi/2,αi,dmax/cos(θi/2))。
其中,所述任意轴旋转单元具体用于:
将以激光扫描仪坐标系为基础的平面映射点Ppi(θi/2,αi,dmax/cos(θi/2))的表示方式转化为以大地坐标系为基础的表示方式Pqi(θi/2+θ,αi+α,dmax/cos(θi/2));
按照下面公式将所述N个三维激光扫描数据平面映射点转化为N个直角坐标Pri(prix,priy,priz):
按照下面步骤进行任意轴坐标旋转:
A1.计算得到Ppi所在的平面法向量为(cosθcosα,cosθsinα,sinθ),Pti所在的水平面法向量为(0,0,1);
A2.计算得到旋转轴A为(cosθsinα,-cosθcosα,0);
A3.计算得到旋转角θr为arccos(sinθ);
A4.进行任意轴旋转:
其中,Ax表示cosθsinα,Ay表示-cosθcosα,Az表示0,C表示cosθr,S表示sinθr。
通过以上描述可知,本发明所述的基于半球面和平面联合投影的三维激光扫描数据建模方法,充分利用了激光扫描数据可投影到一个假想的球面上这一特性,通过半球面和平面二次投影及任意轴旋转后,建模的算法从三维转换为二维,既减小了建模的空间复杂度、时间复杂度及算法难度,同时也保证了建模结果的准确性。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案,下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图是本发明的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图获得其他的附图。
图1示出了本发明实施例一提供的基于半球面和平面联合投影的三维激光扫描数据建模方法的流程图;
图2示出了一个金属矿山地下采空区的三维激光扫描数据点的示意图;
图3示出了三维激光扫描数据在半球面上的映射点的示意图;
图4示出了三维激光扫描数据在平面上的映射点的示意图;
图5示出了经过任意轴旋转后在水平面上的数据点的示意图;
图6示出了水平面上的数据点经过三角剖分后形成的不规则三角网的示意图;
图7示出了水平面上数据点的拓扑关系映射回原始数据点后形成的激光扫描数据建模结果示意图;
图8示出了本发明实施例二提供的基于半球面和平面联合投影的三维激光扫描数据建模系统的结构示意图。
具体实施方式
为使本发明实施例的目的、技术方案和优点更加清楚,下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整的描述,显然,所描述的实施例是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有作出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
三维激光扫描得到的点云数据相对于通常意义的点云数据建模方法而言具有其自身的空间几何特点:激光扫描数据可投影到一个假想的球面上,现有建模方法均未充分利用这一重要特性,现有的建模方法均是从散乱点云的角度重新探究数据点之间的拓扑关系,这既增加了建模的空间复杂度、时间复杂度及算法难度,同时也难以保证建模结果的准确性。针对上述描述,本发明提供了一种基于半球面和平面联合投影的三维激光扫描数据建模方法,以解决现有三维激光扫描数据建模复杂度高、建模所需时间较长、准确率较低的问题。
实施例一
图1示出了本发明实施例一提供的基于半球面和平面联合投影的三维激光扫描数据建模方法的流程图。参见图1,本发明实施例一提供的基于半球面和平面联合投影的三维激光扫描数据建模方法包括:
步骤101:根据扫描头的坐标PO(pox,poy,poz)、扫描头的初始方位(θ,α)及N个三维激光扫描数据点的极坐标Poi(θi,αi,di),得到扫描数据点到所述扫描头的最长距离dmax。
在本步骤中,根据扫描头的坐标PO(pox,poy,poz)及扫描头的初始方位(θ,α)建立激光扫描仪坐标系(其中θ为扫描头初始倾角,α为扫描头初始方位角),N个三维激光扫描数据点的极坐标Poi(θi,αi,di)为以激光扫描仪坐标系为基础的表示方式,将其转化为以大地坐标系为基础的表示方式Pdi(θi+θ,αi+α,di),比较N个三维激光扫描数据点Poi的di值,得到三维激光扫描数据点到扫描头的最长距离dmax。
步骤102:通过半球面投影法则将N个三维激光扫描数据点投影到半球面上,得到N个三维激光扫描数据半球面映射点Psi(θi/2,αi,dmax)。
在本步骤中,根据所述半球面投影法则Poi(θi,αi,di)→Psi(θi/2,αi,dmax)将极坐标(θi,αi,di)中的di值变换为dmax,dmax为一个固定值,数据点(θi,αi,dmax)在一个球面上;且将倾角θi折半,此时数据点(θi,αi,dmax)在球面上收缩为一个半球面,得到半球面映射数据点Psi(θi/2,αi,dmax)。
步骤103:通过平面投影法则将所述N个三维激光扫描数据半球面映射点投影到平面上,得到N个三维激光扫描数据平面映射点Ppi(θi/2,αi,dmax/cos(θi/2))。
在本步骤中,根据所述平面投影法则Psi(θi/2,αi,dmax)→Ppi(θi/2,αi,dmax/cos(θi/2)),将半球面映射数据点(θi/2,αi,dmax)中的dmax值变换为dmax/cos(θi/2),此时数据点(θi/2,αi,dmax/cos(θi/2))便在一个平面上,从而得到平面映射数据点Ppi(θi/2,αi,dmax/cos(θi/2))。
三维激光扫描仪器的性质及扫面原理决定θ的取值为-140°~140°,故0.342<cos(θi/2)≤1,因此dmax/cos(θi/2)有界。
步骤104:将所述N个三维激光扫描数据平面映射点转化为N个直角坐标Pri(prix,priy,priz),并进行任意轴坐标旋转,得到坐标系在水平面上的平面映射点Pti(ptix,ptiy,ptiz)。
在本步骤中,将以激光扫描仪坐标系为基础的平面映射点Ppi(θi/2,αi,dmax/cos(θi/2))的表示方式转化为以大地坐标系为基础的表示方式Pqi(θi/2+θ,αi+α,dmax/cos(θi/2));
按照下面公式将所述N个三维激光扫描数据平面映射点转化为N个直角坐标Pri(prix,priy,priz):
其中,上面将以激光扫描仪坐标系为基础的平面映射点Ppi转换为以大地坐标系为基础的表示方式Pqi是为了将极坐标表示的坐标点转化为直角坐标表示的坐标点。激光扫描仪坐标系下的表达方式无法直接转化为直角坐标系的坐标点。因为激光扫描仪坐标系是一个局部坐标系,而大地坐标系才是一个全局坐标系。
接着,按照下面步骤进行任意轴坐标旋转:
A1.计算得到Ppi所在的平面法向量为(cosθcosα,cosθsinα,sinθ),Pti所在的水平面法向量为(0,0,1);
A2.计算得到旋转轴A为(cosθsinα,-cosθcosα,0);
A3.计算得到旋转角θr为arccos(sinθ);
A4.进行任意轴旋转:
其中,Ax表示cosθsinα,Ay表示-cosθcosα,Az表示0,C表示cosθr,S表示sinθr;
由于Az=0,故任意轴旋转公式可简化为:
步骤105:对水平面上的平面映射点Pti进行Delaunay三角剖分,得到Pti之间的不规则三角网拓扑关系。
步骤106:将所述Pti之间的不规则三角网拓扑关系对应地赋值给Poi,从而得到Poi之间的不规则三角网拓扑关系,完成对三维激光扫描数据点的建模。
下面以一个具体实例来对上面的步骤进行详细介绍。图2示出了一个三维激光扫描数据点的示意图,扫描对象为一个金属矿山地下采空区,扫描方式为等距扫描,具体技术不作详细描述,可参考相关文献或相关专利。
步骤1:该实例中,扫描头的坐标PO为(X:84178.00m,Y:22508.68m,Z:-734.44m),扫描头的初始倾角θ为47.5°,初始方位角α为-12.5°。扫描结果共得到15756个数据点,均以激光扫描仪坐标系为基础记录各数据点的极坐标信息,分析得到扫描数据点到扫描头的最长距离dmax为67.157m。
步骤2:通过半球面投影法则将15756个三维激光扫描数据点投影到半球面上,得到三维激光扫描数据半球面映射点Psi(θi/2,αi,dmax),图3示出了三维激光扫描数据在半球面上的映射点的示意图。
在本步骤中,所述半球面投影法则如下:
Poi(θi,αi,di)→Psi(θi/2,αi,dmax);
首先将极坐标(θi,αi,di)中的di值变换为dmax,由于dmax是一个固定值,故数据点(θi,αi,dmax)便在一个球面上;其次将倾角θi折半,此时数据点(θi,αi,dmax)便在球面上收缩为一个半球面,从而得到半球面映射数据点Psi(θi/2,αi,dmax)。
步骤3:通过平面投影法则将15756个半球面映射点投影到平面上,得到三维激光扫描数据平面映射点Ppi(θi/2,αi,dmax/cos(θi/2)),图4示出了三维激光扫描数据在平面上的映射点的示意图。
在本步骤中,所述平面投影法则如下:
Psi(θi/2,αi,dmax)→Ppi(θi/2,αi,dmax/cos(θi/2));
将极坐标(θi/2,αi,dmax)中的dmax值变换为dmax/cos(θi/2),此时数据点(θi/2,αi,dmax/cos(θi/2))便在一个平面上,从而得到平面映射数据点Ppi(θi/2,αi,dmax/cos(θi/2))。
三维激光扫描仪器的性质及扫面原理决定θ的取值为-140°~140°,故0.342<cos(θi/2)≤1,因此dmax/cos(θi/2)有界。
步骤4:将15756个平面映射点转化为直角坐标Pri(prix,priy,priz),并进行任意轴坐标旋转,得到坐标系在水平面上的平面映射点Pti(ptix,ptiy,ptiz),图5示出了经过任意轴旋转后在水平面上的数据点的示意图。
在本步骤中,所述平面映射点Ppi(θi/2,αi,dmax/cos(θi/2))是以激光扫描仪坐标系为基础的表示方式,可转化为以大地坐标系为基础的表示方式Pqi(θi/2+47.5,αi-12.5,dmax/cos(θi/2))。
极坐标系点转化为直角坐标系点的公式如下:
任意轴旋转的步骤如下:
步骤41:计算得到Ppi所在的平面法向量为(0.660,-0.146,0.737),Pti所在的水平面法向量为(0,0,1);
步骤42:计算得到旋转轴A为(-0.146,-0.660,0);
步骤43:计算得到旋转角θr为42.5°;
步骤44:进行任意轴旋转:
步骤5:对水平面上的映射点Pti进行Delaunay三角剖分,得到Pti之间的不规则三角网拓扑关系,图6示出了水平面上的数据点经过三角剖分后形成的不规则三角网的示意图。
在本步骤中,所述Delaunay三角剖分形成的不规则三角网具有一些重要的性质:(1)空外接圆准则;(2)局部性准则;(3)最小内角最大化准则。
步骤6:将Pti之间的不规则三角网拓扑关系对应地赋值给Poi,从而得到Poi之间的不规则三角网拓扑关系,即完成对三维激光扫描数据点的建模,图7示出了水平面上数据点的拓扑关系映射回原始数据点后形成的激光扫描数据建模结果示意图。
通过本发明实施例描述的基于半球面和平面联合投影的三维激光扫描数据建模方法,可实现充分利用了激光扫描数据可投影到一个假想的球面上这一特性,在通过半球面和平面二次投影及任意轴旋转后,建模的算法从三维转换为二维,既减小了建模的空间复杂度、时间复杂度及算法难度,同时也保证了建模结果的准确性。
实施例二
图8示出了本发明实施例二提供的基于半球面和平面联合投影的三维激光扫描数据建模系统的结构示意图。参见图8,本发明实施例二提供的基于半球面和平面联合投影的三维激光扫描数据建模系统包括:包括预处理单元81、半球面投影单元82、平面投影单元83、任意轴旋转单元84、平面三角拓扑关系获取单元85和三维激光数据建模单元86;
其中,所述预处理单元81用于根据扫描头的坐标PO(pox,poy,poz)、扫描头的初始方位(θ,α)及N个三维激光扫描数据点的极坐标Poi(θi,αi,di),得到扫描数据点到所述扫描头的最长距离dmax;
所述半球面投影单元82用于通过半球面投影法则将N个三维激光扫描数据点投影到半球面上,得到N个三维激光扫描数据半球面映射点Psi(θi/2,αi,dmax);
所述平面投影单元83用于通过平面投影法则将所述N个三维激光扫描数据半球面映射点投影到平面上,得到N个三维激光扫描数据平面映射点Ppi(θi/2,αi,dmax/cos(θi/2));
所述任意轴旋转单元84用于将所述N个三维激光扫描数据平面映射点转化为N个直角坐标Pri(prix,priy,priz),并进行任意轴坐标旋转,得到坐标系在水平面上的平面映射点Pti(ptix,ptiy,ptiz);
所述平面三角拓扑关系获取单元85用于对水平面上的平面映射点Pti进行Delaunay三角剖分,得到Pti之间的不规则三角网拓扑关系;
所述三维激光数据建模单元86用于将所述Pti之间的不规则三角网拓扑关系对应地赋值给Poi,从而得到Poi之间的不规则三角网拓扑关系,完成对三维激光扫描数据点的建模。
其中,所述预处理单元81具体用于:
根据扫描头的坐标PO(pox,poy,poz)及扫描头的初始方位(θ,α)建立激光扫描仪坐标系,N个三维激光扫描数据点的极坐标Poi(θi,αi,di)为以激光扫描仪坐标系为基础的表示方式,将其转化为以大地坐标系为基础的表示方式Pdi(θi+θ,αi+α,di),比较N个三维激光扫描数据点Poi的di值,得到三维激光扫描数据点到扫描头的最长距离dmax。
其中,所述半球面投影单元82具体用于:
根据所述半球面投影法则Poi(θi,αi,di)→Psi(θi/2,αi,dmax),将极坐标(θi,αi,di)中的di值变换为dmax,dmax为一个固定值,数据点(θi,αi,dmax)在一个球面上;且将倾角θi折半,此时数据点(θi,αi,dmax)在球面上收缩为一个半球面,得到半球面映射数据点Psi(θi/2,αi,dmax)。
其中,所述平面投影单元83具体用于:
根据所述平面投影法则Psi(θi/2,αi,dmax)→Ppi(θi/2,αi,dmax/cos(θi/2)),将半球面映射数据点(θi/2,αi,dmax)中的dmax值变换为dmax/cos(θi/2),此时数据点(θi/2,αi,dmax/cos(θi/2))便在一个平面上,从而得到平面映射数据点Ppi(θi/2,αi,dmax/cos(θi/2))。
其中,所述任意轴旋转单元84具体用于:
将以激光扫描仪坐标系为基础的平面映射点Ppi(θi/2,αi,dmax/cos(θi/2))的表示方式转化为以大地坐标系为基础的表示方式Pqi(θi/2+θ,αi+α,dmax/cos(θi/2));
按照下面公式将所述N个三维激光扫描数据平面映射点转化为N个直角坐标Pri(prix,priy,priz):
按照下面步骤进行任意轴坐标旋转:
A1.计算得到Ppi所在的平面法向量为(cosθcosα,cosθsinα,sinθ),Pti所在的水平面法向量为(0,0,1);
A2.计算得到旋转轴A为(cosθsinα,-cosθcosα,0);
A3.计算得到旋转角θr为arccos(sinθ);
A4.进行任意轴旋转:
其中,Ax表示cosθsinα,Ay表示-cosθcosα,Az表示0,C表示cosθr,S表示sinθr。
本发明实施例所述的系统可以执行上述实施例所述的方法,其技术原理和有益效果相似,此处不再赘述。
以上实施例仅用于说明本发明的技术方案,而非对其限制;尽管参照前述实施例对本发明进行了详细的说明,本领域的普通技术人员应当理解:其依然可以对前述各实施例所记载的技术方案进行修改,或者对其中部分技术特征进行等同替换;而这些修改或替换,并不使相应技术方案的本质脱离本发明各实施例技术方案的精神和范围。
Claims (10)
1.一种基于半球面和平面联合投影的三维激光扫描数据建模方法,其特征在于,包括:
S1.根据扫描头的坐标PO(pox,poy,poz)、扫描头的初始方位(θ,α)及N个三维激光扫描数据点的极坐标Poi(θi,αi,di),得到扫描数据点到所述扫描头的最长距离dmax;其中,θ为扫描头初始倾角,α为扫描头初始方位角;
S2.通过半球面投影法则将N个三维激光扫描数据点投影到半球面上,得到N个三维激光扫描数据半球面映射点Psi(θi/2,αi,dmax);
S3.通过平面投影法则将所述N个三维激光扫描数据半球面映射点投影到平面上,得到N个三维激光扫描数据平面映射点Ppi(θi/2,αi,dmax/cos(θi/2));
S4.将所述N个三维激光扫描数据平面映射点转化为N个直角坐标Pri(prix,priy,priz),并进行任意轴坐标旋转,得到坐标系在水平面上的平面映射点Pti(ptix,ptiy,ptiz);
S5.对水平面上的平面映射点Pti进行Delaunay三角剖分,得到Pti之间的不规则三角网拓扑关系;
S6.将所述Pti之间的不规则三角网拓扑关系对应地赋值给Poi,从而得到Poi之间的不规则三角网拓扑关系,完成对三维激光扫描数据点的建模。
2.根据权利要求1所述的基于半球面和平面联合投影的三维激光扫描数据建模方法,其特征在于,所述步骤S1根据扫描头的坐标PO(pox,poy,poz)、扫描头的初始方位(θ,α)及N个三维激光扫描数据点的极坐标Poi(θi,αi,di),得到扫描数据点到所述扫描头的最长距离dmax包括:
根据扫描头的坐标PO(pox,poy,poz)及扫描头的初始方位(θ,α)建立激光扫描仪坐标系,N个三维激光扫描数据点的极坐标Poi(θi,αi,di)为以激光扫描仪坐标系为基础的表示方式,将其转化为以大地坐标系为基础的表示方式Pdi(θi+θ,αi+α,di),比较N个三维激光扫描数据点Poi的di值,得到三维激光扫描数据点到扫描头的最长距离dmax。
3.根据权利要求1所述的基于半球面和平面联合投影的三维激光扫描数据建模方法,其特征在于,所述步骤S2通过半球面投影法则将N个三维激光扫描数据点投影到半球面上,得到N个三维激光扫描数据半球面映射点Psi(θi/2,αi,dmax)包括:
根据所述半球面投影法则Poi(θi,αi,di)→Psi(θi/2,αi,dmax),将极坐标(θi,αi,di)中的di值变换为dmax,dmax为一个固定值,数据点(θi,αi,dmax)落在一个球面上;且将倾角θi折半,此时数据点(θi,αi,dmax)在球面上收缩为一个半球面,得到半球面映射数据点Psi(θi/2,αi,dmax)。
4.根据权利要求1所述的基于半球面和平面联合投影的三维激光扫描数据建模方法,其特征在于,所述步骤S3通过平面投影法则将所述N个三维激光扫描数据半球面映射点投影到平面上,得到N个三维激光扫描数据平面映射点Ppi(θi/2,αi,dmax/cos(θi/2))包括:
根据所述平面投影法则Psi(θi/2,αi,dmax)→Ppi(θi/2,αi,dmax/cos(θi/2)),将半球面映射数据点(θi/2,αi,dmax)中的dmax值变换为dmax/cos(θi/2),此时数据点(θi/2,αi,dmax/cos(θi/2))便在一个平面上,从而得到平面映射数据点Ppi(θi/2,αi,dmax/cos(θi/2))。
5.根据权利要求1所述的基于半球面和平面联合投影的三维激光扫描数据建模方法,其特征在于,所述步骤S4将所述N个三维激光扫描数据平面映射点转化为N个直角坐标Pri(prix,priy,priz),并进行任意轴坐标旋转,得到坐标系在水平面上的平面映射点Pti(ptix,ptiy,ptiz)包括:
将以激光扫描仪坐标系为基础的平面映射点Ppi(θi/2,αi,dmax/cos(θi/2))的表示方式转化为以大地坐标系为基础的表示方式Pqi(θi/2+θ,αi+α,dmax/cos(θi/2));
按照下面公式将所述N个三维激光扫描数据平面映射点转化为N个直角坐标Pri(prix,priy,priz):
按照下面步骤进行任意轴坐标旋转:
A1.计算得到Ppi所在的平面法向量为(cosθcosα,cosθsinα,sinθ),Pti所在的水平面法向量为(0,0,1);
A2.计算得到旋转轴A为(cosθsinα,-cosθcosα,0);
A3.计算得到旋转角θr为arccos(sinθ);
A4.进行任意轴旋转:
其中,Ax表示cosθsinα,Ay表示-cosθcosα,Az表示0,C表示cosθr,S表示sinθr。
6.一种基于半球面和平面联合投影的三维激光扫描数据建模系统,其特征在于,包括预处理单元、半球面投影单元、平面投影单元、任意轴旋转单元、平面三角拓扑关系获取单元和三维激光数据建模单元;
其中,所述预处理单元用于根据扫描头的坐标PO(pox,poy,poz)、扫描头的初始方位(θ,α)及N个三维激光扫描数据点的极坐标Poi(θi,αi,di),得到扫描数据点到所述扫描头的最长距离dmax;其中,θ为扫描头初始倾角,α为扫描头初始方位角;
所述半球面投影单元用于通过半球面投影法则将N个三维激光扫描数据点投影到半球面上,得到N个三维激光扫描数据半球面映射点Psi(θi/2,αi,dmax);
所述平面投影单元用于通过平面投影法则将所述N个三维激光扫描数据半球面映射点投影到平面上,得到N个三维激光扫描数据平面映射点Ppi(θi/2,αi,dmax/cos(θi/2));
所述任意轴旋转单元用于将所述N个三维激光扫描数据平面映射点转化为N个直角坐标Pri(prix,priy,priz),并进行任意轴坐标旋转,得到坐标系在水平面上的平面映射点Pti(ptix,ptiy,ptiz);
所述平面三角拓扑关系获取单元用于对水平面上的平面映射点Pti进行Delaunay三角剖分,得到Pti之间的不规则三角网拓扑关系;
所述三维激光数据建模单元用于将所述Pti之间的不规则三角网拓扑关系对应地赋值给Poi,从而得到Poi之间的不规则三角网拓扑关系,完成对三维激光扫描数据点的建模。
7.根据权利要求6所述的基于半球面和平面联合投影的三维激光扫描数据建模系统,其特征在于,所述预处理单元具体用于:
根据扫描头的坐标PO(pox,poy,poz)及扫描头的初始方位(θ,α)建立激光扫描仪坐标系,N个三维激光扫描数据点的极坐标Poi(θi,αi,di)为以激光扫描仪坐标系为基础的表示方式,将其转化为以大地坐标系为基础的表示方式Pdi(θi+θ,αi+α,di),比较N个三维激光扫描数据点Poi的di值,得到三维激光扫描数据点到扫描头的最长距离dmax。
8.根据权利要求6所述的基于半球面和平面联合投影的三维激光扫描数据建模系统,其特征在于,所述半球面投影单元具体用于:
根据所述半球面投影法则Poi(θi,αi,di)→Psi(θi/2,αi,dmax),将极坐标(θi,αi,di)中的di值变换为dmax,dmax为一个固定值,数据点(θi,αi,dmax)在一个球面上;且将倾角θi折半,此时数据点(θi,αi,dmax)在球面上收缩为一个半球面,得到半球面映射数据点Psi(θi/2,αi,dmax)。
9.根据权利要求6所述的基于半球面和平面联合投影的三维激光扫描数据建模系统,其特征在于,所述平面投影单元具体用于:
根据所述平面投影法则Psi(θi/2,αi,dmax)→Ppi(θi/2,αi,dmax/cos(θi/2)),将半球面映射数据点(θi/2,αi,dmax)中的dmax值变换为dmax/cos(θi/2),此时数据点(θi/2,αi,dmax/cos(θi/2))便在一个平面上,从而得到平面映射数据点Ppi(θi/2,αi,dmax/cos(θi/2))。
10.根据权利要求6所述的基于半球面和平面联合投影的三维激光扫描数据建模系统,其特征在于,所述任意轴旋转单元具体用于:
将以激光扫描仪坐标系为基础的平面映射点Ppi(θi/2,αi,dmax/cos(θi/2))的表示方式转化为以大地坐标系为基础的表示方式Pqi(θi/2+θ,αi+α,dmax/cos(θi/2));
按照下面公式将所述N个三维激光扫描数据平面映射点转化为N个直角坐标Pri(prix,priy,priz):
按照下面步骤进行任意轴坐标旋转:
A1.计算得到Ppi所在的平面法向量为(cosθcosα,cosθsinα,sinθ),Pti所在的水平面法向量为(0,0,1);
A2.计算得到旋转轴A为(cosθsinα,-cosθcosα,0);
A3.计算得到旋转角θr为arccos(sinθ);
A4.进行任意轴旋转:
其中,Ax表示cosθsinα,Ay表示-cosθcosα,Az表示0,C表示cosθr,S表示sinθr。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510143817.9A CN104765915B (zh) | 2015-03-30 | 2015-03-30 | 三维激光扫描数据建模方法及系统 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510143817.9A CN104765915B (zh) | 2015-03-30 | 2015-03-30 | 三维激光扫描数据建模方法及系统 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN104765915A CN104765915A (zh) | 2015-07-08 |
CN104765915B true CN104765915B (zh) | 2017-08-04 |
Family
ID=53647740
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510143817.9A Expired - Fee Related CN104765915B (zh) | 2015-03-30 | 2015-03-30 | 三维激光扫描数据建模方法及系统 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN104765915B (zh) |
Families Citing this family (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105069824B (zh) * | 2015-08-11 | 2017-12-08 | 中南大学 | 基于gps数据的露天矿山道路网自动构建方法及系统 |
JP6679746B2 (ja) * | 2016-03-11 | 2020-04-15 | サイバーオプティクス コーポレーション | 3次元非接触スキャニングシステムのフィールド校正 |
US10235797B1 (en) | 2017-11-27 | 2019-03-19 | Lowe's Companies, Inc. | Inverse rendering of visual material properties |
CN109253717B (zh) * | 2018-10-09 | 2020-11-27 | 安徽大学 | 一种矿区地表沉降三维激光扫描地表沉降监测设站方法 |
CN109514091B (zh) * | 2018-12-29 | 2020-07-31 | 广州新可激光设备有限公司 | 一种激光打标机的3d扫描建模方法、系统及其打标机 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102136155A (zh) * | 2010-01-27 | 2011-07-27 | 首都师范大学 | 基于三维激光扫描数据的物体立面矢量化方法和系统 |
CN103279989A (zh) * | 2013-05-30 | 2013-09-04 | 北京航天控制仪器研究所 | 一种三维激光成像系统平面点云数据三角化处理方法 |
CN104266605A (zh) * | 2014-06-27 | 2015-01-07 | 西北工业大学 | 一种三维激光扫描成像仪的成像方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US7340316B2 (en) * | 2004-06-28 | 2008-03-04 | Hanger Orthopedic Group, Inc. | System and method for producing medical devices |
-
2015
- 2015-03-30 CN CN201510143817.9A patent/CN104765915B/zh not_active Expired - Fee Related
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102136155A (zh) * | 2010-01-27 | 2011-07-27 | 首都师范大学 | 基于三维激光扫描数据的物体立面矢量化方法和系统 |
CN103279989A (zh) * | 2013-05-30 | 2013-09-04 | 北京航天控制仪器研究所 | 一种三维激光成像系统平面点云数据三角化处理方法 |
CN104266605A (zh) * | 2014-06-27 | 2015-01-07 | 西北工业大学 | 一种三维激光扫描成像仪的成像方法 |
Also Published As
Publication number | Publication date |
---|---|
CN104765915A (zh) | 2015-07-08 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN104765915B (zh) | 三维激光扫描数据建模方法及系统 | |
CN105136054B (zh) | 基于地面三维激光扫描的构筑物精细变形监测方法及系统 | |
CN106932784B (zh) | 基于二维激光雷达的车斗形容器三维扫描系统测量方法 | |
CN106960468B (zh) | 一种三维激光扫描点云精度评价方法 | |
CN104330073B (zh) | 一种烟囱倾斜度的测量方法 | |
Schwalbe et al. | 3D building model generation from airborne laser scanner data using 2D GIS data and orthogonal point cloud projections | |
CN104697502B (zh) | 基于最小二乘法的建筑物特征点坐标提取方法 | |
CN103196426A (zh) | 一种全站仪联合三维激光扫描仪用于建筑测量的方法 | |
CN102353678B (zh) | 一种测量文物病害的方法 | |
CN102044088B (zh) | 单站地面激光扫描海量散乱点云的lod模型快速构建方法 | |
CN102855663A (zh) | 一种根据激光雷达栅格点云构建csg模型的方法 | |
CN103617649B (zh) | 一种基于相机自标定技术的河工模型地形测量方法 | |
CN106885531A (zh) | 基于二维激光雷达的车斗形容器三维扫描系统标定方法 | |
CN107631700B (zh) | 三维扫描仪与全站仪相结合的三维数据测量方法 | |
CN111105496A (zh) | 一种基于机载激光雷达点云数据的高精dem构建方法 | |
CN105469404A (zh) | 一种基于三维点云数据的旋转体拟合方法及装置 | |
CN104751479A (zh) | 基于tin数据的建筑物提取方法和装置 | |
CN102298793A (zh) | 一种获取文物三维图形数据的方法 | |
CN103106632A (zh) | 一种基于均值漂移的不同精度三维点云数据的融合方法 | |
CN106097433A (zh) | 物体工业与影像模型的叠加方法及系统 | |
CN109146990B (zh) | 一种建筑轮廓的计算方法 | |
CN112068153A (zh) | 一种基于地基激光雷达点云的冠层间隙率估算方法 | |
CN103315770A (zh) | 基于三维超声图像的目标参数处理方法及系统 | |
CN103335606A (zh) | 三维激光扫描测绘方法 | |
CN108955520A (zh) | 一种结构光三维扫描可达性分析方法及分析系统 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
EXSB | Decision made by sipo to initiate 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: 20170804 |
|
CF01 | Termination of patent right due to non-payment of annual fee |