CN107240077B - 一种基于椭圆构像偏差迭代矫正的视觉测量方法 - Google Patents
一种基于椭圆构像偏差迭代矫正的视觉测量方法 Download PDFInfo
- Publication number
- CN107240077B CN107240077B CN201710408105.4A CN201710408105A CN107240077B CN 107240077 B CN107240077 B CN 107240077B CN 201710408105 A CN201710408105 A CN 201710408105A CN 107240077 B CN107240077 B CN 107240077B
- Authority
- CN
- China
- Prior art keywords
- ellipse
- pose
- correction
- iteration
- matrix
- 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
- 238000000691 measurement method Methods 0.000 title claims abstract description 14
- 230000000007 visual effect Effects 0.000 title description 4
- 239000011159 matrix material Substances 0.000 claims description 28
- 238000000034 method Methods 0.000 claims description 17
- 230000009466 transformation Effects 0.000 claims description 6
- 238000004364 calculation method Methods 0.000 claims description 5
- 238000001514 detection method Methods 0.000 claims description 3
- 239000000126 substance Substances 0.000 claims description 3
- 238000005259 measurement Methods 0.000 abstract description 20
- 238000010586 diagram Methods 0.000 description 3
- 230000003287 optical effect Effects 0.000 description 3
- 238000003384 imaging method Methods 0.000 description 2
- 238000004590 computer program Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 230000036544 posture Effects 0.000 description 1
- 238000000638 solvent extraction Methods 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T5/00—Image enhancement or restoration
- G06T5/80—Geometric correction
-
- 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
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/80—Analysis of captured images to determine intrinsic or extrinsic camera parameters, i.e. camera calibration
-
- 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/30—Subject of image; Context of image processing
- G06T2207/30244—Camera pose
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Length Measuring Devices By Optical Means (AREA)
Abstract
本发明公开一种基于椭圆构像偏差迭代矫正的视觉测量方法,该方法利用迭代实现了基于偏心差矫正的位姿估计,可以与单目相机标定、双目相机标定、位姿测量等应用结合,实现高精度的视觉测量;可以计算任意位置的圆形标志的偏心误差,因此对于具体测量应用中圆形标志的位置没有限制;通过调整迭代终止条件可以获得不同精度和速度要求的位姿估计结果。
Description
技术领域
本发明涉及计算机视觉测量领域,尤其涉及一种基于椭圆构像偏差迭代矫正的视觉测量方法。
背景技术
圆经过透视投影,当成像平面与圆平面不平行时投影为椭圆,然而投影椭圆的中心并非圆心在像平面上的投影,两者的偏差称为椭圆构像偏差。现有的视觉测量方法是将椭圆中心近似作为圆的中心进行位姿的计算,该近似过程所引入的误差在高精度测量要求下不能被忽略。文献“Eccentricity error identification and compensation forhigh-accuracy 3D optical measurement”(He D,Liu X,Peng X,et al.MeasurementScience&Technology,vol.24(7),075402,2013.)提出一种利用同心圆环进行摄像机标定的方法,通过获取的椭圆特征分析圆形标志投影中心,其优点是需要检测的特征少,缺点是在实际相机标定中精度较差,无法用于高精度的测量任务。文献“Eccentricity erroridentification and compensation for high-accuracy 3D optical measurement”(HeD,Liu X,Peng X,et al.Measurement Science&Technology,vol.24(7),075402,2013.)提出一种基于同心圆环的椭圆构像偏差矫正方法,但其中涉及近似计算,会对椭圆中心点的提取精度造成影响。文献“Eccentricity error compensation for geometric cameracalibration based on circular features”(Yang X,Fang S.Measurement Science&Technology,vol.25(2),pp.149-156,2014.)提出了椭圆构像偏差在相机标定中的影响,并给出了一种矫正的方式,但其主要在矫正相机畸变引起的椭圆中心偏差,没有很好的解决由于透视投影所造成的构像偏差问题。通过椭圆构像偏差模型可以计算出任意位姿下的椭圆构像偏差,但是前提是需要知道位姿信息,这一前提在实际工程测量应用中无法被满足。
发明内容
本发明的目的在于通过一种基于椭圆构像偏差迭代矫正的视觉测量方法,来解决以上背景技术部分提到的问题。
为达此目的,本发明采用以下技术方案:
一种基于椭圆构像偏差迭代矫正的视觉测量方法,其包括如下步骤:
S101、基于图像检测获得的椭圆中心作为特征点作出初始位姿估计,计算出基于初始位姿的椭圆构像偏差;
S102、对检测到的椭圆中心进行椭圆构像偏差矫正,获得矫正后的特征点,并作出新的位姿估计;
S103、进行迭代终止条件判断,若不满足终止条件,则基于新的位姿计算新的椭圆构像偏差值继续进行迭代计算;若满足终止条件,则得到最终结果。
特别地,所述步骤S101包括:对图像上的椭圆特征进行检测,获得椭圆中心点坐标pi,所有椭圆中心点的偏心误差初始化为μi0=0,i=0,1,2,L,位姿初始化为R=03×3,T=[00 0]T。
特别地,所述步骤S102包括:每次迭代对所有椭圆中心点进行椭圆构像偏差矫正,对于第k次矫正过程,矫正后的椭圆中心点坐标Pik=Pi-μik;由矫正后的椭圆中心点坐标Pik及标靶上圆形标志中心点坐标Qi计算出第k次矫正后的位姿Rk、Tk。
特别地,所述步骤S103中迭代终止条件包括:一、迭代次数达到设定上限值;二、所有标志点的椭圆偏心误差μik小于设定阈值;三、两次迭代所得位姿Rk-1、Tk-1、Rk、Tk的差值E小于设定阈值;当满足上述迭代终止条件的任一项时迭代均终止,得到最终估计位姿Rk、Tk,否则由Rk、Tk计算出偏心误差μi,k+1并进行第k+1次迭代矫正;其中,
为k-1次迭代计算得到的位姿变换矩阵,为k次迭代计算得到的位姿变换矩阵,为其逆矩阵;H矩阵为相邻两次迭代结果的差值;将H矩阵进行分块,D矩阵为其前三行前三列组成的3X3的矩阵,Xe,Ye,Ze分别为H矩阵第4列前三行的元素,而d11,d22,d33分别指代D矩阵中对角线上的三个元素。
本发明提出的基于椭圆构像偏差迭代矫正的视觉测量方法具有如下优点:一、利用迭代实现了基于偏心差矫正的位姿估计,可以与单目相机标定、双目相机标定、位姿测量等应用结合,实现高精度的视觉测量;二、可以计算任意位置的圆形标志的偏心误差,因此对于具体测量应用中圆形标志的位置没有限制;三、通过调整迭代终止条件可以获得不同精度和速度要求的位姿估计结果。
附图说明
图1为本发明实施例提供的基于椭圆构像偏差迭代矫正的视觉测量方法流程图;
图2为本发明实施例提供的单相机位姿测量示意图;
图3为本发明实施例提供的基于椭圆构像偏差迭代矫正的相机标定流程图;
图4为本发明实施例提供的单相机标定示意图;
图5A-图5D为本发明实施例提供的单相机标定时所拍摄图片的示意图。
具体实施方式
下面结合附图和实施例对本发明作进一步说明。可以理解的是,此处所描述的具体实施例仅仅用于解释本发明,而非对本发明的限定。另外还需要说明的是,为了便于描述,附图中仅示出了与本发明相关的部分而非全部内容,除非另有定义,本文所使用的所有技术和科学术语与属于本发明的技术领域的技术人员通常理解的含义相同。本文中所使用的术语只是为了描述具体的实施例,不是旨在于限制本发明。
请参照图1所示,图1为本发明实施例提供的基于椭圆构像偏差迭代矫正的视觉测量方法流程图。
本实施例中基于椭圆构像偏差迭代矫正的视觉测量方法包括如下步骤:
S101、基于图像检测获得的椭圆中心作为特征点作出初始位姿估计,计算出基于初始位姿的椭圆构像偏差。
S102、对检测到的椭圆中心进行椭圆构像偏差矫正,获得矫正后的特征点,并作出新的位姿估计。
S103、进行迭代终止条件判断,若不满足终止条件,则基于新的位姿计算新的椭圆构像偏差值继续进行迭代计算;若满足终止条件,则得到最终结果。
具体的,所述步骤S101包括:对图像上的椭圆特征进行检测,获得椭圆中心点坐标pi,所有椭圆中心点的偏心误差初始化为μi0=0,i=0,1,2,L,位姿初始化为R=03×3,T=[00 0]T。所述步骤S102包括:每次迭代对所有椭圆中心点进行椭圆构像偏差矫正,对于第k次矫正过程,矫正后的椭圆中心点坐标Pik=Pi-μik;由矫正后的椭圆中心点坐标Pik及标靶上圆形标志中心点坐标Qi计算出第k次矫正后的位姿Rk、Tk,需要说明的是,这里的位姿计算方法采用的是PnP方法,具体内容可参考文献:Hesch J A,Roumeliotis S I.A DirectLeast-Squares(DLS)method for PnP[C]//IEEE International Conference onComputer Vision.IEEE,2011:383-390。所述步骤S103中迭代终止条件包括:一、迭代次数达到设定上限值;二、所有标志点的椭圆偏心误差μik小于设定阈值;三、两次迭代所得位姿Rk-1、Tk-1、Rk、Tk的差值E小于设定阈值;当满足上述迭代终止条件的任一项时迭代均终止,得到最终估计位姿Rk、Tk,否则由Rk、Tk计算出偏心误差μi,k+1并进行第k+1次迭代矫正;其中,
为k-1次迭代计算得到的位姿变换矩阵,为k次迭代计算得到的位姿变换矩阵,为其逆矩阵;H矩阵为相邻两次迭代结果的差值;将H矩阵进行分块,D矩阵为其前三行前三列组成的3X3的矩阵,Xe,Ye,Ze分别为H矩阵第4列前三行的元素,而d11,d22,d33分别指代D矩阵中对角线上的三个元素。
本发明提出的基于椭圆构像偏差迭代矫正的视觉测量方法可以应用于位姿测量、单目相机标定、立体视觉标定等领域,下面分别对本发明应用于位姿测量、单目相机标定的具体实施方案进行详细说明:
当基于椭圆构像偏差迭代矫正的视觉测量方法应用于位姿测量时,实现过程如下:
一、如图2所示,相机11拍摄平面对象12上的圆形特征13。由于平面对象12可能不平行于成像平面,成像特征可能会是椭圆。在拍摄图像中提取椭圆特征,获取其中心点坐标。所有椭圆中心点的偏心误差初始化μi0=0,位姿初始化为R=03×3,T=[0 0 0]T。
二、对于所有的椭圆中心进行椭圆偏心误差进行第k次矫正,矫正之后的椭圆中心Pik=Pi-μik。由矫正后的椭圆中心点坐标Pik及标靶上圆形标志中心点坐标Qi可以计算出第k次矫正后的位姿Rk、Tk。
三、对于如下三项迭代终止条件:一、迭代次数达到设定上限值,在本实施例中设定上限值为10;二、所有标志点的椭圆偏心误差μik小于设定阈值,在本实施例中设定阈值为10-5;三、两次迭代所得位姿差值小于设定阈值,本实施例中设定阈值为10-5;当满足上述迭代终止条件的任一项时迭代均终止,得到最终估计位姿Rk、Tk,否则返回第二步,由Rk、Tk计算出偏心误差μi,k+1并进行第k+1次迭代矫正。
下表为椭圆构像偏差矫正前后的位姿估计结果,经过多次椭圆构像偏差的迭代矫正,估计的位姿结果相对于矫正前有了很大的提高。
当基于椭圆构像偏差迭代矫正的视觉测量方法应用于相机标定时,实现过程如图3所示:
一、如图4所示,相机21拍摄不同位置姿态的平面对象22上的圆形特征23,其中部分图像如图5A-图5D所示,这组图片是相机标定过程中拍摄的图片,图5A-图5D是4次在不同角度下拍摄到的图片。在拍摄图像中提取椭圆特征,获取其中心坐标Pji,i其中代表一副图像中椭圆中心点的编号,j代表所拍摄图片的编号。
二、对第一步中拍摄到的图像进行相机标定,标定得到相机畸变参数k1、k2、p1、p2。根据标定出的相机畸变参数进行去除畸变的处理,获得畸变矫正后的椭圆中心点坐标Dji。所有椭圆中心点的偏心误差初始化jμi0=0,位姿初始化为Rj=03×3,Tj=[0 0 0]T。
三、对于每幅图像中所有的椭圆中心进行椭圆偏心误差的第k次矫正,矫正之后的椭圆中心jDik=Dji-jμik。将相机畸变参数固定为0,由矫正后的椭圆中心点坐标jDik及标靶上圆形标志中心点坐标Qi可以标定出第k次矫正后的相机内参和对应第j幅图像的外参jRk、jTk。
四、对于如下三项迭代终止条件:对于如下三项迭代终止条件:一、迭代次数达到设定上限值,在本实施例中设定上限值为10;二、所有标志点的椭圆偏心误差jμik小于设定阈值,在本实施例中设定阈值为10-5;三、对于每幅图像,两次迭代所得位姿差值小于设定阈值,本实施例中设定阈值为10-5;当满足上述迭代终止条件的任一项时迭代均终止,得到最终估计内参与外参jRk、jTk,否则回到第三步进行第k+1次迭代标定。
五、基于第四步得到的相机内外参数对每一幅拍摄得到的图像进行椭圆构像偏差矫正,以矫正后的椭圆中心点坐标jDik及标靶上圆形标志中心点坐标Qi标定出相对于第二步更为精确的相机畸变参数k1、k2、p1、p2。
下表中三栏分别表示设置的相机内参数、未经过椭圆构像偏差矫正的标定结果和经过椭圆构像偏差矫正的标定结果。可以看出经过椭圆构像偏差矫正后相机的标定精度得到了很大的提高。
本发明提出的技术方案具有如下优点:一、利用迭代实现了基于偏心差矫正的位姿估计,可以与单目相机标定、双目相机标定、位姿测量等应用结合,实现高精度的视觉测量;二、可以计算任意位置的圆形标志的偏心误差,因此对于具体测量应用中圆形标志的位置没有限制;三、通过调整迭代终止条件可以获得不同精度和速度要求的位姿估计结果。
本领域普通技术人员可以理解实现上述实施例方法中的全部或部分流程,是可以通过计算机程序来指令相关的硬件来完成,所述的程序可存储于一计算机可读取存储介质中,该程序在执行时,可包括如上述各方法的实施例的流程。其中,所述的存储介质可为磁碟、光盘、只读存储记忆体(Read-Only Memory,ROM)或随机存储记忆体(Random AccessMemory,RAM)等。
以上结合具体实施例描述了本发明的技术原理。这些描述只是为了解释本发明的原理,而不能以任何方式解释为对本发明保护范围的限制。基于此处的解释,本领域的技术人员不需要付出创造性的劳动即可联想到本发明的其它具体实施方式,这些方式都将落入本发明的保护范围之内。
Claims (1)
1.一种基于椭圆构像偏差迭代矫正的视觉测量方法,其特征在于,包括如下步骤:
S101、基于图像检测获得的椭圆中心作为特征点作出初始位姿估计,计算出基于初始位姿的椭圆构像偏差;其中,所述步骤S101包括:对图像上的椭圆特征进行检测,获得椭圆中心点坐标pi,所有椭圆中心点的偏心误差初始化为μi0=0,i=0,1,2,…,位姿初始化为R=03×3,T=[0 0 0]T;
S102、对检测到的椭圆中心进行椭圆构像偏差矫正,获得矫正后的特征点,并作出新的位姿估计;其中,所述步骤S102包括:每次迭代对所有椭圆中心点进行椭圆构像偏差矫正,对于第k次矫正过程,矫正后的椭圆中心点坐标Pik=Pi-μik;由矫正后的椭圆中心点坐标Pik及标靶上圆形标志中心点坐标Qi计算出第k次矫正后的位姿Rk、Tk;
S103、进行迭代终止条件判断,若不满足终止条件,则基于新的位姿计算新的椭圆构像偏差值继续进行迭代计算;若满足终止条件,则得到最终结果;其中,所述步骤S103中迭代终止条件包括:一、迭代次数达到设定上限值;二、所有标志点的椭圆偏心误差μik小于设定阈值;三、两次迭代所得位姿Rk-1、Tk-1、Rk、Tk的差值E小于设定阈值;当满足上述迭代终止条件的任一项时迭代均终止,得到最终估计位姿Rk、Tk,否则由Rk、Tk计算出偏心误差μi,k+1并进行第k+1次迭代矫正;其中,
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710408105.4A CN107240077B (zh) | 2017-06-02 | 2017-06-02 | 一种基于椭圆构像偏差迭代矫正的视觉测量方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710408105.4A CN107240077B (zh) | 2017-06-02 | 2017-06-02 | 一种基于椭圆构像偏差迭代矫正的视觉测量方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN107240077A CN107240077A (zh) | 2017-10-10 |
CN107240077B true CN107240077B (zh) | 2021-02-02 |
Family
ID=59985972
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710408105.4A Active CN107240077B (zh) | 2017-06-02 | 2017-06-02 | 一种基于椭圆构像偏差迭代矫正的视觉测量方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107240077B (zh) |
Families Citing this family (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110009680B (zh) * | 2019-02-28 | 2022-04-22 | 中国人民解放军国防科技大学 | 基于圆特征及异面特征点的单目图像位置、姿态测量方法 |
CN112241944B (zh) * | 2019-07-16 | 2024-03-29 | 杭州海康机器人股份有限公司 | Agv调度系统、物料位姿矫正方法和物料位姿矫正装置 |
CN112489111B (zh) * | 2020-11-25 | 2024-01-30 | 深圳地平线机器人科技有限公司 | 相机外参标定方法、装置及相机外参标定系统 |
Family Cites Families (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US7159776B2 (en) * | 2005-01-05 | 2007-01-09 | Dell Products L.P. | System and method for optical medium label alignment |
CN102184545B (zh) * | 2011-06-02 | 2012-11-07 | 浙江大学 | 折反射全向相机镜面位姿的单图自标定方法 |
US10512441B2 (en) * | 2013-08-29 | 2019-12-24 | Siemens Healthcare Gmbh | Computed tomography having motion compensation |
CN103615998B (zh) * | 2013-12-13 | 2016-08-17 | 西安工业大学 | 齿轮测量中心工件装夹倾斜与偏心误差测量与补偿方法 |
CN104331896B (zh) * | 2014-11-21 | 2017-09-08 | 天津工业大学 | 一种基于深度信息的系统标定方法 |
CN106780576B (zh) * | 2016-11-23 | 2020-03-17 | 北京航空航天大学 | 一种面向rgbd数据流的相机位姿估计方法 |
-
2017
- 2017-06-02 CN CN201710408105.4A patent/CN107240077B/zh active Active
Non-Patent Citations (3)
Title |
---|
Detection of Concentric Circles for Camera Calibration;Guang JIANG 等;《Proceedings of the Tenth IEEE International Conference on Computer Vision (ICCV’05) 》;20051231;第1-8页 * |
Eccentricity error compensation for geometric camera calibration based oncircular features;Xieliu Yang等;《Measurement Science and Technology》;20141231;第25卷;第1-10页 * |
基于射影变换圆阵靶标中心像点的计算;陈天飞等;《仪器仪表学报》;20150430;第36 卷(第4 期);第895-902页 * |
Also Published As
Publication number | Publication date |
---|---|
CN107240077A (zh) | 2017-10-10 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107633536B (zh) | 一种基于二维平面模板的相机标定方法及系统 | |
CN109272570B (zh) | 一种基于立体视觉数学模型的空间点三维坐标求解方法 | |
CN105096329B (zh) | 一种精确校正超广角摄像头图像畸变的方法 | |
WO2018209968A1 (zh) | 摄像机标定方法及系统 | |
CN111524194B (zh) | 一种激光雷达和双目视觉相互融合的定位方法及终端 | |
CN107240077B (zh) | 一种基于椭圆构像偏差迭代矫正的视觉测量方法 | |
CN110099267A (zh) | 梯形校正系统、方法以及投影仪 | |
CN107633533B (zh) | 大畸变镜头下的高精度圆形标志点中心定位方法及装置 | |
CN113920205B (zh) | 一种非同轴相机的标定方法 | |
CN107871329B (zh) | 一种相机光学中心的快速标定方法及装置 | |
CN112686961A (zh) | 一种深度相机标定参数的修正方法、装置 | |
CN111062966A (zh) | 基于l-m算法和多项式插值对相机跟踪进行优化的方法 | |
CN112929626A (zh) | 一种基于智能手机影像的三维信息提取方法 | |
CN112529019A (zh) | 基于直线特征和关键点特征融合的图像拼接方法和系统 | |
CN111538029A (zh) | 一种视觉与雷达融合的测量方法及终端 | |
CN114998773A (zh) | 适用于无人机系统航拍图像的特征误匹配剔除方法及系统 | |
CN113962853A (zh) | 一种旋转线阵扫描图像位姿自动化精密解算方法 | |
CN108520541B (zh) | 一种广角摄像机的标定方法 | |
CN107993268B (zh) | 一种相机自标定的方法及系统 | |
CN110232715B (zh) | 一种多深度相机自校准的方法、装置及系统 | |
CN115311369A (zh) | 一种低成本适用少量重合区域的高精度多rgb-d相机的外参标定方法 | |
CN113920196A (zh) | 视觉定位方法、装置及计算机设备 | |
CN114897990A (zh) | 一种基于神经网络的相机畸变标定方法、系统和存储介质 | |
CN114511637A (zh) | 一种基于构造强特征的弱特征物体图像三维重建系统及方法 | |
CN112927299B (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 | ||
EE01 | Entry into force of recordation of patent licensing contract | ||
EE01 | Entry into force of recordation of patent licensing contract |
Application publication date: 20171010 Assignee: Jiangsu Jihui Huake Intelligent Equipment Technology Co., Ltd. Assignor: Wuxi research institute of the Central China University of Science and Technology Contract record no.: X2019980000288 Denomination of invention: Visual measurement method based on iterative correction of ellipse imaging deviation License type: Common License Record date: 20191023 |
|
GR01 | Patent grant | ||
GR01 | Patent grant |