CN106846406A - 一种用于物体三维坐标测量的点云坐标系构建方法 - Google Patents
一种用于物体三维坐标测量的点云坐标系构建方法 Download PDFInfo
- Publication number
- CN106846406A CN106846406A CN201710102789.5A CN201710102789A CN106846406A CN 106846406 A CN106846406 A CN 106846406A CN 201710102789 A CN201710102789 A CN 201710102789A CN 106846406 A CN106846406 A CN 106846406A
- Authority
- CN
- China
- Prior art keywords
- coordinate system
- point cloud
- cloud
- construction method
- coordinates
- 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
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/10—Image acquisition modality
- G06T2207/10028—Range image; Depth image; 3D point clouds
Landscapes
- Length Measuring Devices By Optical Means (AREA)
Abstract
本发明涉及一种用于物体三维坐标测量的点云坐标系构建方法,包括以下步骤:1)获得某相机坐标系下编码标志群点云的三维坐标,标志点p=(xi,yi,zi),i=0,1,…,n,n为标志点总数;2)计算所述点云的重心po,以该重心作为点云坐标系的原点;3)计算协方差矩阵的特征值λ1、λ2、λ3,并获得与所述特征值相对应的三个单位特征向量η1、η2、η3,以η1、η2、η3作为点云坐标系的X、Y、Z轴。与现有技术相比,本发明具有可提高大型复杂物件摄影测量实际可操作性等优点。
Description
技术领域
本发明属于摄影测量领域,尤其是涉及一种用于物体三维坐标测量的点云坐标系构建方法。
背景技术
在摄影测量中,相机的外参数标定是实现点云拼接的重要前提,外参数标定的本质就是标定相机坐标系之间的转换矩阵或相机坐标系与物方坐标系之间的空间转换矩阵。摄影测量中传统相机标定法,根据不同的标定参照物可分为三维靶标定法、二维靶标定法二种。三维靶标定法需要已知靶标上的标志点在靶标坐标系的三维坐标,且加工、维护比较困难。二维靶标定法是三维靶标定法的简化,但也需要知道二维靶面内所有标志点在靶标坐标系的坐标。因此二维靶上标志点分布比较规则或不规则但需要通过其他方式先测得标定板上标志点的位置。如***在“数字近景工业摄影测量理论、方法与应用”(科学出版社,2016)提到的二维标定板,均匀分布9个编码标志和48个定向反光标定点,需要先使用V-STARS/S测量出所有标志点在标定靶坐标系的坐标。在实际测量中,通常需要粘贴大量信息唯一编码标志以实现大量照片的拼接,这两种二维标定靶都不实用。如何准确有效地描述位置关系随机、复杂的编码标志群内标志点云的位置关系,需要建立一种点云三维坐标系。
发明内容
本发明的目的就是为了克服上述现有技术存在的缺陷而提供一种可提高大型复杂物件摄影测量实际可操作性的用于物体三维坐标测量的点云坐标系构建方法。
本发明的目的可以通过以下技术方案来实现:
一种用于物体三维坐标测量的点云坐标系构建方法,包括以下步骤:
1)获得某相机坐标系下编码标志群点云的三维坐标,标志点p=(xi,yi,zi),i=0,1,…,n,n为标志点总数;
2)计算所述点云的重心po,以该重心作为点云坐标系的原点;
3)计算协方差矩阵的特征值λ1、λ2、λ3,并获得与所述特征值相对应的三个单位特征向量η1、η2、η3,以η1、η2、η3作为点云坐标系的X、Y、Z轴。
所述点云的重心po的计算公式为:
对于同一编码标志群,在不同相机坐标系获得的特征值顺序相同时,依次取三个特征值对应的特征向量为坐标系的X、Y、Z轴。
对于同一编码标志群,当利用不同相机坐标系前后两次获得的特征值顺序不同时,将两次获得的特征值顺序调整一致并相应调整对应的单位特征向量的顺序和方向后依次作为坐标系的X、Y、Z轴。
对于同一编码标志群,在不同相机坐标系下建立的点云坐标系为全局世界坐标系下的同一坐标系。
与现有技术相比,本发明具有以下优点:
1、在实际测量过程中往往需要灵活粘贴定位编码标志,通过本发明方法建立的三维点云坐标系,可以有效表达编码标志组上标志点云在局部世界坐标系下的位置关系,从而大大提高了对大型复杂物件摄影测量的实际可操作性。
2、本发明方法建立的三维点云坐标系为局部世界坐标系,它能准确表达位置关系随机、复杂的点云在该坐标系下的位置关系,从而方便地建立起换站前、换站后相机坐标系与局部世界坐标系的相互变换关系。
附图说明
图1为本发明点云三维坐标系与换站前、换站后相机坐标系的相对关系示意图;
图中:1为相机观测到的编码标志群点云;2为本发明建立的点云坐标系;3为换站前0号机位相机坐标系;4为换站后1号机位相机坐标系;
图2为本发明的方法的流程示意图。
具体实施方式
下面结合附图和具体实施例对本发明进行详细说明。本实施例以本发明技术方案为前提进行实施,给出了详细的实施方式和具体的操作过程,但本发明的保护范围不限于下述的实施例。
如图1所示,0号相机和1号相机有重叠的拍照,同时拍到了一群编码标志,可以从0号相机坐标系的角度建立该编码标志点云三维坐标系,也可以从1号相机坐标系角度建立该坐标系。点云坐标系在0号相机坐标系和1号相机坐标系下位置是不同的,但在全局世界坐标系的位置是唯一确定的。如图2所示,下面分别从0号相机坐标下和1号相机坐标下建立点云的三维坐标系:
1、从0号相机坐标下建立点云的三维坐标系
1)先求得0号相机坐标系下编码标志群中点云的三维坐标为P0=(x0i,y0i,z0i),i=0,1,…,n;
2)求得该点云重心的三维坐标作为该点云坐标系的原点;
3)求得该点云坐标的协方差矩阵求得cov0的三个特征值λ1,λ2,λ3对应的三个单位特征向量η1,η2,η3依次作为该点云三维坐标系的X、Y、Z轴。
2、从1号相机坐标下建立点云的三维坐标系
1)先求得1号相机坐标系下编码标志群点云的三维坐标为P1=(x1i,y1i,z1i),i=0,1,…,n;
2)求得该点云重心的三维坐标作为该点云坐标系的原点;
3)求得该点云坐标的协方差矩阵求得cov1的三个特征值顺序若也为λ1,λ2,λ3则对应的特征向量ξ1,ξ2,ξ3作为该点云三维坐标系的X、Y、Z轴,即当前后二次求得的特征值顺序刚好相同时,则依次取对应的特征向量为X、Y、Z轴即可)。若求得的三个特征值顺序为λ2,λ1,λ3,对应单位特征向量ζ1,ζ2,ζ3,应调整cov0或cov1特征值顺序使得同为λ1,λ2,λ3或者λ2,λ1,λ3,并相应调整特征值所对应单位特征向量的顺序和方向作为该点云三维坐标系的X、Y、Z轴,当前后两次求得特征值顺序不同时,不仅要调整二次的特征值顺序保持一致,再取对应特征向量作为XYZ轴即可,这时可能需要把第一次求得特征值顺序调整为第二次求得的特征值顺序,也可能是反之,在调整特征值顺序的同时,应同时调整对应的特征向量顺序,并且特征向量可能还应反号。
3、以上两次建立的点云三维坐标系为全局世界坐标系下的同一坐标系。
本发明建立的点云坐标系可应用于一种双目摄影测量方法中,该双目摄影测量方法用于标定换站前、换站后相机坐标系转换关系的小型编码标志是灵活粘贴在被测物上或周边的,编码标志之间的相对位置不定,实质构成一个编码标志点之间的相互位置关系不确定的整体三维标定靶。本发明方法建立的三维点云坐标系能准确表达位置关系随机、复杂的点云在该坐标系下的位置关系,从而建立起换站前、换站后相机坐标系与局部世界坐标系的相互变换关系。
以上详细描述了本发明的较佳具体实施例。应当理解,本领域的普通技术人员无需创造性劳动就可以根据本发明的构思作出诸多修改和变化。因此,凡本技术领域中技术人员依本发明的构思在现有技术的基础上通过逻辑分析、推理或者有限的实验可以得到的技术方案,皆应在由权利要求书所确定的保护范围内。
Claims (5)
1.一种用于物体三维坐标测量的点云坐标系构建方法,其特征在于,包括以下步骤:
1)获得某相机坐标系下编码标志群点云的三维坐标,标志点p=(xi,yi,zi),i=0,1,…,n,n为标志点总数;
2)计算所述点云的重心po,以该重心作为点云坐标系的原点;
3)计算协方差矩阵的特征值λ1、λ2、λ3,并获得与所述特征值相对应的三个单位特征向量η1、η2、η3,以η1、η2、η3作为点云坐标系的X、Y、Z轴。
2.根据权利要求1所述的用于物体三维坐标测量的点云坐标系构建方法,其特征在于,所述点云的重心po的计算公式为:
3.根据权利要求1所述的用于物体三维坐标测量的点云坐标系构建方法,其特征在于,对于同一编码标志群,在不同相机坐标系获得的特征值顺序相同时,依次取三个特征值对应的特征向量为坐标系的X、Y、Z轴。
4.根据权利要求1所述的用于物体三维坐标测量的点云坐标系构建方法,其特征在于,对于同一编码标志群,当利用不同相机坐标系前后两次获得的特征值顺序不同时,将两次获得的特征值顺序调整一致并相应调整对应的单位特征向量的顺序和方向后依次作为坐标系的X、Y、Z轴。
5.根据权利要求1所述的用于物体三维坐标测量的点云坐标系构建方法,其特征在于,对于同一编码标志群,在不同相机坐标系下建立的点云坐标系为全局世界坐标系下的同一坐标系。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710102789.5A CN106846406A (zh) | 2017-02-24 | 2017-02-24 | 一种用于物体三维坐标测量的点云坐标系构建方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710102789.5A CN106846406A (zh) | 2017-02-24 | 2017-02-24 | 一种用于物体三维坐标测量的点云坐标系构建方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN106846406A true CN106846406A (zh) | 2017-06-13 |
Family
ID=59133833
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710102789.5A Pending CN106846406A (zh) | 2017-02-24 | 2017-02-24 | 一种用于物体三维坐标测量的点云坐标系构建方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106846406A (zh) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2020237647A1 (zh) * | 2019-05-31 | 2020-12-03 | 深圳市大疆创新科技有限公司 | 点云数据的编码方法、编码装置、解码方法与解码装置 |
CN112731448A (zh) * | 2020-11-17 | 2021-04-30 | 南京恩瑞特实业有限公司 | 一种列车障碍物探测线路数据生成方法及*** |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105205858A (zh) * | 2015-09-18 | 2015-12-30 | 天津理工大学 | 一种基于单个深度视觉传感器的室内场景三维重建方法 |
CN105243637A (zh) * | 2015-09-21 | 2016-01-13 | 武汉海达数云技术有限公司 | 一种基于三维激光点云进行全景影像拼接方法 |
CN105931237A (zh) * | 2016-04-19 | 2016-09-07 | 北京理工大学 | 一种图像校准方法和*** |
CN105976353A (zh) * | 2016-04-14 | 2016-09-28 | 南京理工大学 | 基于模型和点云全局匹配的空间非合作目标位姿估计方法 |
CN106441200A (zh) * | 2016-07-20 | 2017-02-22 | 杭州先临三维科技股份有限公司 | 一种多测量模式的三维测量方法 |
-
2017
- 2017-02-24 CN CN201710102789.5A patent/CN106846406A/zh active Pending
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105205858A (zh) * | 2015-09-18 | 2015-12-30 | 天津理工大学 | 一种基于单个深度视觉传感器的室内场景三维重建方法 |
CN105243637A (zh) * | 2015-09-21 | 2016-01-13 | 武汉海达数云技术有限公司 | 一种基于三维激光点云进行全景影像拼接方法 |
CN105976353A (zh) * | 2016-04-14 | 2016-09-28 | 南京理工大学 | 基于模型和点云全局匹配的空间非合作目标位姿估计方法 |
CN105931237A (zh) * | 2016-04-19 | 2016-09-07 | 北京理工大学 | 一种图像校准方法和*** |
CN106441200A (zh) * | 2016-07-20 | 2017-02-22 | 杭州先临三维科技股份有限公司 | 一种多测量模式的三维测量方法 |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2020237647A1 (zh) * | 2019-05-31 | 2020-12-03 | 深圳市大疆创新科技有限公司 | 点云数据的编码方法、编码装置、解码方法与解码装置 |
CN112731448A (zh) * | 2020-11-17 | 2021-04-30 | 南京恩瑞特实业有限公司 | 一种列车障碍物探测线路数据生成方法及*** |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN108510551B (zh) | 一种远距离大视场条件下相机参数的标定方法及*** | |
CN102589530B (zh) | 基于二维相机和三维相机融合的非合作目标位姿测量方法 | |
CN100429476C (zh) | 一种双传感器激光视觉三维测量***校准方法 | |
CN108038902A (zh) | 一种面向深度相机的高精度三维重建方法和*** | |
CN110142805A (zh) | 一种基于激光雷达的机器人末端校准方法 | |
CN102646275A (zh) | 通过跟踪和定位算法实现虚拟三维叠加的方法 | |
CN106898025A (zh) | 一种基于八点编码标志的相机位移变换矩阵标定方法 | |
CN102645173A (zh) | 一种基于多目视觉的桥梁三维变形监测方法 | |
CN109099883A (zh) | 高精度大视场机器视觉测量与标定装置及方法 | |
CN103530880A (zh) | 基于投影高斯网格图案的摄像机标定方法 | |
CN103759669A (zh) | 一种大型零件的单目视觉测量方法 | |
CN103278138A (zh) | 一种复杂结构薄部件三维位置及姿态的测量方法 | |
CN104182969B (zh) | 一种单线阵相机内外参数标定方法 | |
CN101425185B (zh) | 一种基于复合型平面靶标的小尺寸视觉测量的摄像机标定方法 | |
CN110889873A (zh) | 一种目标定位方法、装置、电子设备及存储介质 | |
CN104807405B (zh) | 一种基于光线角度标定的三维坐标测量方法 | |
CN103852060A (zh) | 一种基于单目视觉的可见光图像测距方法 | |
CN105654547A (zh) | 三维重建方法 | |
CN102519434A (zh) | 一种用于立体视觉三维恢复数据精度测量的试验验证方法 | |
CN104123725B (zh) | 一种单线阵相机单应性矩阵h的计算方法 | |
CN102914295A (zh) | 基于计算机视觉立方体标定的三维测量方法 | |
CN105139355A (zh) | 一种深度图像的增强方法 | |
CN103606147A (zh) | 多台不共视场量测相机的坐标系转换标定方法 | |
CN101325724B (zh) | 一种立体图像对的极线校正方法 | |
Zheng et al. | Minimal solvers for 3d geometry from satellite imagery |
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 | ||
RJ01 | Rejection of invention patent application after publication | ||
RJ01 | Rejection of invention patent application after publication |
Application publication date: 20170613 |