CN113390394B - 一种带有摄影比例尺的光束法平差算法 - Google Patents
一种带有摄影比例尺的光束法平差算法 Download PDFInfo
- Publication number
- CN113390394B CN113390394B CN202110648377.8A CN202110648377A CN113390394B CN 113390394 B CN113390394 B CN 113390394B CN 202110648377 A CN202110648377 A CN 202110648377A CN 113390394 B CN113390394 B CN 113390394B
- Authority
- CN
- China
- Prior art keywords
- camera
- object point
- scale
- initial value
- point
- 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.)
- Active
Links
- 238000004422 calculation algorithm Methods 0.000 title claims abstract description 16
- 238000000034 method Methods 0.000 claims abstract description 6
- 238000003384 imaging method Methods 0.000 claims description 13
- 238000005457 optimization Methods 0.000 claims description 7
- 238000012897 Levenberg–Marquardt algorithm Methods 0.000 claims description 2
- 238000004364 calculation method Methods 0.000 abstract 1
- 239000003550 marker Substances 0.000 description 19
- 238000005259 measurement Methods 0.000 description 3
- 230000007547 defect Effects 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 239000000523 sample Substances 0.000 description 1
- 230000000007 visual effect Effects 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01C—MEASURING DISTANCES, LEVELS OR BEARINGS; SURVEYING; NAVIGATION; GYROSCOPIC INSTRUMENTS; PHOTOGRAMMETRY OR VIDEOGRAMMETRY
- G01C11/00—Photogrammetry or videogrammetry, e.g. stereogrammetry; Photographic surveying
Landscapes
- Engineering & Computer Science (AREA)
- Multimedia (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Length Measuring Devices By Optical Means (AREA)
Abstract
本发明公开了一种带有摄影比例尺的光束法平差算法。首先利用常规三维重建算法得到被测物点的三维坐标初值,然后根据各被测物点的三维坐标初值计算出各物点在对应相机中成像的摄影比例尺,并建立带有摄影比例尺的目标函数,优化相机参数和被测物点三维坐标初值,得到最终的相机参数和三维坐标。该方法适用于工业摄影测量中相机与标志点的距离变化较大时的三维坐标重建,且计算结果可靠。
Description
技术领域
本发明涉及视觉测量技术领域,尤其涉及一种带有摄影比例尺的光束法平差算法。
背景技术
在航空航天、汽车制造等领域中,工业摄影测量技术常常需要在测量现场布置一定数量的编码标志点,利用相机从多个方位采集标志点图像,再对标志点进行三维重建。为了提升重建精度,需要利用光束法平差对重建初值进行非线性优化,获得高精度的三维坐标值。但常规光束法平差在系统的实际应用中存在以下限制:
在实际的系统应用中,相机光心在各方位与同一个标志点之间的距离不同,该标志点在不同成像平面上相同大小的重投影误差所对应的空间定位误差大小也会不同,因此各成像面重投影误差所表征的空间定位误差无法统一在相同尺度上。常规光束法平差以成像平面上的重投影误差作为误差方程,在相机与标志点距离变化较大的情况下,难以对最终重建结果的空间定位误差进行评价。
发明内容
本发明目的就是为了弥补已有技术的缺陷,提供一种带有摄影比例尺的光束法平差算法,以增强光束法平差在相机与标志点距离变化较大情况下对标志点空间定位误差的优化,实现对标志点的高精度测量。
本发明是通过以下技术方案实现的:
一种带有摄影比例尺的光束法平差算法,具体包括以下步骤:
(1)、利用相机标定算法得到相机的内、外参数初值,利用三维重建算法得到n个被测物点的三维坐标初值;
(2)、根据各被测物点的三维坐标初值计算出各物点在对应相机中成像的摄影比例尺,并建立带有摄影比例尺的误差方程及目标函数;
(3)、根据建立的误差方程和目标函数,对相机的内、外参数初值和被测物点的三维坐标初值进行整体的非线性优化,得到最终优化后的相机内、外参数和三维坐标。
步骤(2)中所述的建立带有摄影比例尺的误差方程及目标函数,具体步骤如下:
(2a)设空间中第i个物点Pi的初值为(Xi,Yi,Zi),第j个相机的主距为(fx,fy),计算得到物点Pi在第j个相机中成像的摄影比例尺为mx=Zi/fx,my=Zi/fy;
(2b)设空间中第i个物点Pi在第j个相机的成像平面上的像点坐标为pij(uij,vij),根据共线方程计算得到第i个物点Pi在第j个相机的成像平面上的重投影点pij′(uij′,vij′),第i个物点Pi在第j个相机中成像的带有摄影比例尺的误差方程为得到在m台相机和n个被测物点情况下的目标函数
利用Levenberg-Marquardt算法对相机的内、外参数初值和被测物点的三维坐标初值进行整体的非线性优化,得到最终优化后的相机内、外参数和三维坐标。
本发明的优点是:本发明可以对标志点的空间定位误差进行优化;本发明不受相机与标志点距离变化的影响,可以更好地评价标志点空间定位误差。
附图说明
图1为本发明的方法流程图。
图2为摄影比例尺示意图。
具体实施方式
如图1、2所示,一种带有摄影比例尺的光束法平差算法,具体包括以下步骤:
(1)、利用相机标定算法得到相机的内、外参数初值,利用三维重建算法得到n个被测物点的三维坐标初值;
(2)、根据各被测物点的三维坐标初值计算出各物点在对应相机中成像的摄影比例尺,并建立带有摄影比例尺的误差方程及目标函数;
(3)、根据建立的误差方程和目标函数,对相机的内、外参数初值和被测物点的三维坐标初值进行整体的非线性优化,得到最终优化后的相机内、外参数和三维坐标。
现结合附图和具体实施例对本发明进一步详细的说明。
实施例一
本实施例是对三维扫描测头上的标志点进行三维重建,在4台相机,6个标志点的情况下,利用算法对相机参数和物点三维坐标进行优化,包括以下步骤:
步骤(1):对相机进行标定得到相机的内、外参数初值(C1,C2,C3,C4),对标志点进行三维重建,得到各点的三维坐标初值(P1,P2,P3,P4);
步骤(2):根据各标志点的三维坐标初值计算出各物点在对应相机中成像的摄影比例尺,并建立带有摄影比例尺的误差方程以及目标函数,包括如下具体步骤:
(2a)设空间中第i个标志点Pi的初值为(Xi,Yi,Zi),第j个相机的主距为(fx,fy),计算得到标志点Pi在第j个相机中成像的摄影比例尺为mx=Zi/fn,my=Zi/fy;
(2b)设空间中第i个物点Pi在第j个相机的成像平面上的像点坐标为pij(uij,vij),根据共线方程计算得到第i个物点Pi在第j个相机的成像平面上的重投影点pij′(uij′,vij′),第i个物点Pi在第j个相机中成像的带有摄影比例尺的误差方程为得到在m台相机和n个被测物点情况下的目标函数
(3)、根据目标函数,利用LM(Levenberg-Marquardt)算法对4个相机参数(C1,C2,...,Cm)和6个标志点三维坐标初值(P1,P2,...,Pn)进行整体的非线性优化,得到最终优化后的相机参数和三维坐标。
以上实施方式仅用于说明本发明,而非对本发明的限制,有关技术领域的普通技术人员,在不脱离本发明的精神和范围的情况下,还可以做出各种变化和变型,因此所有等同的技术方案也属于本发明的保护范畴。
Claims (1)
1.一种带有摄影比例尺的光束法平差算法,其特征在于:具体包括以下步骤:
(1)、利用相机标定算法得到相机的内、外参数初值,利用三维重建算法得到n个被测物点的三维坐标初值;
(2)、根据各被测物点的三维坐标初值计算出各物点在对应相机中成像的摄影比例尺,并建立带有摄影比例尺的误差方程及目标函数;
(3)、根据建立的误差方程和目标函数,对相机的内、外参数初值和被测物点的三维坐标初值进行整体的非线性优化,得到最终优化后的相机内、外参数和三维坐标;
步骤(2)中所述的建立带有摄影比例尺的误差方程及目标函数,具体步骤如下:
(2a)设空间中第i个物点Pi的初值为(Xi,Yi,Zi),第j个相机的主距为(fx,fy),计算得到物点Pi在第j个相机中成像的摄影比例尺为mx=Zi/fx,my=Zi/fy;
(2b)设空间中第i个物点Pi在第j个相机的成像平面上的像点坐标为pij(uij,vij),根据共线方程计算得到第i个物点Pi在第j个相机的成像平面上的重投影点pij′(uij′,vij′),第i个物点Pi在第j个相机中成像的带有摄影比例尺的误差方程为得到在m台相机和n个被测物点情况下的目标函数
利用Levenberg-Marquardt算法对相机的内、外参数初值和被测物点的三维坐标初值进行整体的非线性优化,得到最终优化后的相机内、外参数和三维坐标。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110648377.8A CN113390394B (zh) | 2021-06-10 | 2021-06-10 | 一种带有摄影比例尺的光束法平差算法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110648377.8A CN113390394B (zh) | 2021-06-10 | 2021-06-10 | 一种带有摄影比例尺的光束法平差算法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN113390394A CN113390394A (zh) | 2021-09-14 |
CN113390394B true CN113390394B (zh) | 2022-12-06 |
Family
ID=77620213
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202110648377.8A Active CN113390394B (zh) | 2021-06-10 | 2021-06-10 | 一种带有摄影比例尺的光束法平差算法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN113390394B (zh) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN114596341B (zh) * | 2022-02-16 | 2024-07-23 | 合肥工业大学 | 面向大视场运动目标的多相机高精度三维位姿跟踪方法 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103759670A (zh) * | 2014-01-06 | 2014-04-30 | 四川虹微技术有限公司 | 一种基于数字近景摄影的物体三维信息获取方法 |
CN104931022A (zh) * | 2015-04-21 | 2015-09-23 | 国家测绘地理信息局卫星测绘应用中心 | 基于星载激光测高数据的卫星影像立体区域网平差方法 |
CN105627934A (zh) * | 2014-10-30 | 2016-06-01 | 宁波舜宇光电信息有限公司 | 一种基于机器视觉的视觉比例系数获取方法 |
CN108871207A (zh) * | 2018-09-14 | 2018-11-23 | 佛山市天目工业测量科技有限公司 | 一种摄影测量基准尺长度标定装置及使用方法 |
WO2019179200A1 (zh) * | 2018-03-22 | 2019-09-26 | 深圳岚锋创视网络科技有限公司 | 多目相机的三维重建方法、vr相机和全景相机 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20150070468A1 (en) * | 2013-09-10 | 2015-03-12 | Faro Technologies, Inc. | Use of a three-dimensional imager's point cloud data to set the scale for photogrammetry |
-
2021
- 2021-06-10 CN CN202110648377.8A patent/CN113390394B/zh active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103759670A (zh) * | 2014-01-06 | 2014-04-30 | 四川虹微技术有限公司 | 一种基于数字近景摄影的物体三维信息获取方法 |
CN105627934A (zh) * | 2014-10-30 | 2016-06-01 | 宁波舜宇光电信息有限公司 | 一种基于机器视觉的视觉比例系数获取方法 |
CN104931022A (zh) * | 2015-04-21 | 2015-09-23 | 国家测绘地理信息局卫星测绘应用中心 | 基于星载激光测高数据的卫星影像立体区域网平差方法 |
WO2019179200A1 (zh) * | 2018-03-22 | 2019-09-26 | 深圳岚锋创视网络科技有限公司 | 多目相机的三维重建方法、vr相机和全景相机 |
CN108871207A (zh) * | 2018-09-14 | 2018-11-23 | 佛山市天目工业测量科技有限公司 | 一种摄影测量基准尺长度标定装置及使用方法 |
Non-Patent Citations (1)
Title |
---|
多目视觉测量系统的光束法平差改进;别梓钒等;《计测技术》;20210831;23-27 * |
Also Published As
Publication number | Publication date |
---|---|
CN113390394A (zh) | 2021-09-14 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN102034238B (zh) | 基于光学成像测头和视觉图结构的多摄像机系统标定方法 | |
CN109242908B (zh) | 用于水下双目视觉测量系统的标定方法 | |
CN105913439B (zh) | 一种基于激光跟踪仪的大视场摄像机标定方法 | |
CN101943563B (zh) | 基于空间平面约束的线结构光视觉传感器的快速标定方法 | |
CN107576281B (zh) | 一种测量管路弯曲半径的方法和装置 | |
CN109163657B (zh) | 一种基于双目视觉三维重建的圆形目标位姿检测方法 | |
CN116051659B (zh) | 一种线阵相机与2d激光扫描仪联合标定方法 | |
CN104266608B (zh) | 视觉传感器现场标定装置和标定方法 | |
CN101231750A (zh) | 一种双目立体测量系统的标定方法 | |
CN110398208A (zh) | 基于摄影全站仪系统的大数据变形监测方法 | |
CN114705122B (zh) | 一种大视场立体视觉标定方法 | |
CN108648242B (zh) | 基于激光测距仪辅助无公共视场的两相机标定方法及装置 | |
CN103606147B (zh) | 多台不共视场量测相机的坐标系转换标定方法 | |
CN105716542A (zh) | 一种基于柔性特征点的三维数据拼接方法 | |
CN103559707B (zh) | 基于运动方靶标定物的工业定焦相机参数标定方法 | |
CN109900205B (zh) | 一种高精度的单线激光器和光学相机的快速标定方法 | |
CN113269673B (zh) | 一种基于标准球架的三维点云拼接方法 | |
CN105486289A (zh) | 一种激光摄影测量系统及相机标定方法 | |
CN110345921A (zh) | 立体视场视觉测量及垂轴像差和轴向像差校正方法及系统 | |
CN104697463A (zh) | 一种双目视觉传感器的消隐特征约束标定方法及装置 | |
CN114018167A (zh) | 一种基于单目三维视觉的桥梁挠度测量方法 | |
CN113390394B (zh) | 一种带有摄影比例尺的光束法平差算法 | |
CN109342008B (zh) | 基于单应性矩阵的风洞试验模型迎角单相机视频测量方法 | |
CN110827359A (zh) | 基于棋盘格三面体的相机与激光器外参检校方法及装置 | |
CN113160301B (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 |