CN1620153A - 利用平面控制点场进行非量测数码相机标定的方法 - Google Patents
利用平面控制点场进行非量测数码相机标定的方法 Download PDFInfo
- Publication number
- CN1620153A CN1620153A CN 200410061262 CN200410061262A CN1620153A CN 1620153 A CN1620153 A CN 1620153A CN 200410061262 CN200410061262 CN 200410061262 CN 200410061262 A CN200410061262 A CN 200410061262A CN 1620153 A CN1620153 A CN 1620153A
- Authority
- CN
- China
- Prior art keywords
- digital camera
- plane
- image
- coordinate
- points
- 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
Images
Landscapes
- Studio Devices (AREA)
Abstract
利用平面控制点场进行非量测数码相机标定的方法,包括以下步骤:利用待标定的数码相机获取控制场的影像;利用边缘提取和图像匹配技术从图像中提取控制点的影像坐标;通过二维直接线性变换和共线方程间的对应关系求出数码相机的内外方位元素初值;按照严密光束法平差原理进行数码相机的高精度标定,从而求得数码相机内外方位元素的精确值;实现非量测数码相机标定。本发明的优点是不需要传统的很难建立的三维控制场,代之以比较容易实现的二维平面控制点场,实现非量测数码相机的高精度标定。
Description
技术领域
本发明涉及一种利用平面控制点场进行非量测数码相机标定的方法,尤其是应用于数字近景摄影测量等领域。
背景技术
检查和较准相机的内方位元素和光学畸变参数的过程称为相机标定。对于日益广泛使用的非量测数码相机来说,标定是从二维图像获取三维信息必不可少的步骤。长期以来,相机的标定都是在实验室采用物理方法进行,成本高,自动化程度低,周期长,无法做到实用化。传统的准直管和精密测角仪器设备昂贵,相机标定必须在专门的实验室进行,无法在作业现场进行。而利用三维控制场进行相机标定时,高精度三维控制场的建立很困难,控制场的移动很不方便,存放也占用较大空间。计算机视觉界提出了一些不需要控制场的相机自标定方法,但这些自标定方法共同的缺点是标定精度较低,无法满足摄影测量的需要。可以说,简便高精度的相机标定方法是人们一直以来的追求,也是数码相机能够真正广泛应用于数字近景摄影测量等领域的关键因素之一。
发明内容
本发明的目的是提供一种利用平面控制点场进行非量测数码相机标定的方法,该方法不仅简单且能够满足数字近景摄影测量的精度要求。
本发明提供的技术方案是,一种利用平面控制点场进行非量测数码相机标定的方法,包括以下步骤:
1.利用待标定的数码相机按照下述方法获取4张以上控制场的影像;
以其上标记有若干坐标已知的特征点(如十字丝的交点、格网点)的平面作为平面控制点场,以垂直于平面控制场的直线作为Z轴,数码相机的主光线与Z轴的夹角不大于30度;沿Z轴与平面控制点场的交点拍照,对于每张像片,数码相机的位置及三个旋转角则互不相同。
2.利用现有技术中的边缘提取和图像匹配技术从图像中提取控制点的影像坐标。
3.利用已知的控制点平面坐标和图像中提取的影像坐标,通过二维直接线性变换和共线方程问的对应关系求出数码相机的内外方位元素初值。
二维直接线性变换表示的是空间平面和像平面间的映射关系
其中H=(h1,h2,h3,h4,h5,h6,h7,h8)T为二维直接线性变换的八个变换参数,X,Y为平面控制点空间坐标(Z坐标为零),x,y为相应的像坐标。当像片点数大于等于4个时,即可通过解线性方程AH=0求得二维直接线性变换参数。
摄影测量中最常用的共线方程为:
其中x0,y0,f为摄像机的内方位元素(主点和焦距);(XS,YS,ZS)为摄站坐标;(X,Y,Z)为物方空间坐标;(x,y)为相应的像点坐标;R={ai,bi,ci,i=1,2,3}为摄影测量中常用的旋转角,ω,κ(Y为主轴)构成的旋转矩阵。对于平面场而言,一般将坐标系建立在控制场中心,则X,Y轴在平面内,Z轴与平面垂直,此时所有控制点的Z坐标为零。
将共线方程进行变换即可得到与二维直接线性变换相同的表达形式,因而每一项都可由二维直接线性变换参数表示。根据旋转矩阵的性质
和
可得:
因此如果有两张以上像片,就可以利用上式求解主点(x0,y0)的初值。
根据旋转矩阵的性质a1b1+a2b2+a3b3=0即可导出焦距的求解方程:
求出主点及焦距的初值后,即可进一步求得外方位元素的初值。在Y为主轴的转角***下,
sinω=-b3,
因而三个旋转角的值都可以唯一确定。XS,YS,ZS的初值时可通过解如下线性方程组获得:
4.利用已知的控制点平面坐标和图像中提取的影像坐标,以及上述求得的数码相机内外方位元素初值,按照下述公式所表达的严密光束法平差原理进行数码相机的高精度标定,从而求得数码相机内外方位元素的精确值。
其中数码相机的镜头畸变差为:
Δx=(x-x0)(K1r2+K2r4)+P1(r2+2(x-x0)2)+2P2(x-x0)·(y-y0)
Δy=(y-y0)(K1r2+K2r4)+P2(r2+2(y-y0)2)+2P1(x-x0)·(y-y0)
本发明所涉及的平面控制点场是指在某一平面上布设4个以上坐标已知的特征点,用来作为标定或三维重建的绝对控制。控制场的载体理论上可以是任何大小的薄板,为了移动方便,实用上一般以1米见方为宜。薄板可以用有机玻璃等平整物体,薄板本身为浅色,如白色或浅灰色。在薄板上面布设一定数量的深色控制点,如黑色。控制点场可以用一般机械工厂都具有的加工设备刻划和加工,控制点的坐标可以用精密坐标量测仪器量测得到。
在本发明的方法中数码相机的移动可以用平面控制场的反移动来代替,两者是等价的,增加了应用的灵活性,这在数码相机必须保持固定的场合具有得天独厚的优势。
本发明的优点是不需要传统的很难建立的三维控制场,代之以制造简单,移动非常方便的二维平面控制点场进行非量测数码相机的高精度离线(off line)和在线(on line)标定该方法不仅简单且能够满足数字近景摄影测量的精度要求。如果硬件由计算机控制,则整个标定过程只需少量人工干预,总耗时小于5分钟。本发明的优点是可以利用利用本发明方法进行标定后的数码相机,三维重建时目标物的重建精度与控制场本身的精度相当。
本发明对工作环境的要求较低,只需正常的照明条件(如自然光或普通日光灯),图像拍摄过程中控制场或数码相机其中一个保持固定不动即可。对温度,湿度,压力等无特殊要求。
附图说明
图1为本发明的流程图;
图2为本发明方法中非量测数码相机相对于二维控制点的位置示意图;
图3为本发明方法中格网点式二维控制点场的示意图。
具体实施方式
参见图1、图2和图3,本发明包括以下步骤:
1.利用待标定的数码相机按照下述方法获取4张以上控制场的影像;
以其上标记有4个以上(为了获得稳定可靠的结果,实用上以每平方米300个左右为最佳)坐标已知的特征点(如十字丝的交点2、格网点4)的平面作为平面控制点场3,以垂直于平面控制场的直线作为Z轴,数码相机1的主光线与Z轴的夹角不大于30度;沿Z轴与平面控制点场的交点拍照;对于每张像片,数码相机的位置及三个旋转角则互不相同;即每拍一张只需移动一下相机位置,并旋转相机的朝向。
2.利用现有技术中的边缘提取和图像匹配技术从图像中提取控制点的影像坐标。具体实施方法为先利用图像处理中的Canny算子提取影像边缘,然后将这些边缘拟合成线段,再利用图像匹配即可得到精确的直线段,然后将构成控制点的十字丝或格网线的两条直线段相交即可获得控制点的精确影像坐标。
3.利用已知的控制点平面坐标和图像中提取的影像坐标,通过二维直接线性变换和共线方程间的对应关系求出数码相机的内外方位元素初值。
二维直接线性变换表示的是空间平面和像平面间的映射关系
当像片点数大于等于4个时,即可通过解线性方程AH=0求得二维直接线性变换参数。
摄影测量中最常用的共线方程为:
其中x0,y0,f为摄像机的内方位元素(主点和焦距);(XS,YS,ZS)为摄站坐标;(X,Y,Z)为物方空间坐标;(x,y)为相应的像点坐标;R={ai,bi,ci,i=1,2,3}为摄影测量中常用的旋转角,ω,κ(Y为主轴)构成的旋转矩阵。对于平面场而言,一般将坐标系建立在控制场中心,则X,Y轴在平面内,Z轴与平面垂直,此时所有控制点的Z坐标为零。
将共线方程进行变换即可得到与二维直接线性变换相同的表达形式,因而每一项都可由二维直接线性变换参数表示。根据旋转矩阵的性质
和
可得:
因此如果有两张以上像片,就可以利用上式求解主点(x0,y0)的初值。
根据旋转矩阵的性质a1b1+a2b2+a3b3=0即可导出焦距的求解方程:
4.利用已知的控制点平面坐标和图像中提取的影像坐标,以及上述求得的数码相机内外方位元素初值,即可按照下述公式表达的严密光束法平差原理进行数码相机的高精度标定,从而求得数码相机内外方位元素的精确值。
其中镜头畸变差为:
Δx=(x-x0)(K1r2+K2r4)+P1(r2+2(x-x0)2)+2P2(x-x0)·(y-y0)
Δy=(y-y0)(K1r2+K2r4)+P2(r2+2(y-y0)2)+2P1(x-x0)·(y-y0)
Claims (2)
1.一种利用平面控制点场进行非量测数码相机标定的方法,其特征是:包括以下步骤
一、利用待标定的数码相机按照下述方法获取4张以上控制场的影像;
以其上标记有4个以上坐标已知的特征点的平面作为平面控制点场,以垂直于平面控制场的直线作为Z轴,数码相机的主光线与Z轴的夹角不大于30度;沿Z轴与平面控制点场的交点拍照,对于每张像片,数码相机的位置及三个旋转角则互不相同;
二、利用现有技术中的边缘提取和图像匹配技术从图像中提取控制点的影像坐标;
三、利用已知的控制点平面坐标和图像中提取的影像坐标,通过二维直接线性变换和共线方程间的对应关系求出数码相机的内外方位元素初值;
二维直接线性变换表示的是空间平面和像平面间的映射关系
其中H=(h1,h2,h3,h4,h5,h6,h7,h8)T为二维直接线性变换的八个变换参数,X,Y为平面控制点空间坐标,x,y为相应的像坐标。当像片点数大于4个时,可将上式进行适当变换,通过解超定方程AH=0求得二维直接线性变换参数;
摄影测量中最常用的共线方程为:
其中x0,y0,f为摄像机的内方位元素,即主点和焦距;(XS,YS,ZS)为摄站坐标;(X,Y,Z)为物方空间坐标;(x,y)为相应的像点坐标;R={ai,bi,ci,i=1,2,3}为摄影测量中常用的旋转角,ω,κ构成的旋转矩阵;对于平面场而言,坐标系建立在控制场中心,则X,Y轴在平面内,Z轴与平面垂直,此时所有控制点的Z坐标为零;
将共线方程进行变换即可得到与二维直接线性变换相同的表达形式,因而每一项都可由二维直接线性变换参数表示;根据旋转矩阵的性质
和 可得:
对于两张以上像片,利用上式求解主点(x0,y0)的初值;
根据旋转矩阵的性质a1b1+a2b1+a3b3=0即可导出焦距的求解方程:
求出主点及焦距的初值后,进一步求得外方位元素的初值;在Y为主轴的转角***下,
sinω=-b3,
因而三个旋转角的值都可以唯一确定;解如下线性方程组获得XS,YS,ZS的初值:
四、利用已知的控制点平面坐标和图像中提取的影像坐标,以及上述求得的数码相机内外方位元素初值,按照下述公式所表达的严密光束法平差原理进行数码相机的高精度标定,从而求得数码相机内外方位元素的精确值;实现非量测数码相机标定
其中镜头畸变差为:
Δx=(x-x0)(K1r2+K2r4)+P1(r2+2(x-x0)2)+2P2(x-x0)·(y-y0)
Δy=(y-y0)(K1r2+K2r4)+P2(r2+2(y-y0)2)+2P1(x-x0)·(y-y0)。
2.根据权利要求1所述的方法,其特征是:平面控制点场的特征点为十字丝的交点或格网点。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN 200410061262 CN1620153A (zh) | 2004-12-02 | 2004-12-02 | 利用平面控制点场进行非量测数码相机标定的方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN 200410061262 CN1620153A (zh) | 2004-12-02 | 2004-12-02 | 利用平面控制点场进行非量测数码相机标定的方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN1620153A true CN1620153A (zh) | 2005-05-25 |
Family
ID=34764469
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN 200410061262 Pending CN1620153A (zh) | 2004-12-02 | 2004-12-02 | 利用平面控制点场进行非量测数码相机标定的方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN1620153A (zh) |
Cited By (15)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN100409658C (zh) * | 2006-09-08 | 2008-08-06 | 南京大学 | 图像缩放器中抗边缘锯齿噪声的方法 |
CN101226638B (zh) * | 2007-01-18 | 2010-05-19 | 中国科学院自动化研究所 | 一种对多相机***的标定方法及装置 |
CN102012236A (zh) * | 2010-09-26 | 2011-04-13 | 郑州辰维科技股份有限公司 | 一种月球车双目视觉避障***标定方法 |
CN101551907B (zh) * | 2009-04-28 | 2011-05-04 | 浙江大学 | 一种多照相机自动化高精度标定方法 |
CN101561269B (zh) * | 2009-05-26 | 2011-08-17 | 空气动力学国家重点实验室 | 近景摄影测量的特征线自动匹配方法 |
CN101661617B (zh) * | 2008-08-30 | 2011-11-02 | 华为终端有限公司 | 摄像机标定的方法及装置 |
CN101876555B (zh) * | 2009-11-04 | 2011-12-21 | 北京控制工程研究所 | 一种月球车双目视觉导航***标定方法 |
CN102865857A (zh) * | 2012-09-04 | 2013-01-09 | 北京信息科技大学 | 一种摄影测量图像匹配方法 |
CN102901490A (zh) * | 2012-09-04 | 2013-01-30 | 北京信息科技大学 | 一种基于动态阈值的图像匹配方法及*** |
CN103175544A (zh) * | 2013-03-04 | 2013-06-26 | 中国科学院长春光学精密机械与物理研究所 | 用地面控制点计算和标定航天tdiccd相机偏流角的方法 |
CN104729532A (zh) * | 2015-03-02 | 2015-06-24 | 山东科技大学 | 一种全景相机的严密标定方法 |
CN105424058A (zh) * | 2015-11-06 | 2016-03-23 | 中国人民解放军信息工程大学 | 基于摄影测量技术的数码相机投影中心位置精确标定方法 |
CN106023142A (zh) * | 2016-05-06 | 2016-10-12 | 北京信息科技大学 | 利用共面直线阵列标定摄影测量相机的算法 |
CN109087339A (zh) * | 2018-06-13 | 2018-12-25 | 武汉朗视软件有限公司 | 一种激光扫描点与影像配准方法 |
CN112308894A (zh) * | 2019-08-01 | 2021-02-02 | 司法鉴定科学研究院 | 一种基于近景摄影测量的图像处理***及处理方法 |
-
2004
- 2004-12-02 CN CN 200410061262 patent/CN1620153A/zh active Pending
Cited By (20)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN100409658C (zh) * | 2006-09-08 | 2008-08-06 | 南京大学 | 图像缩放器中抗边缘锯齿噪声的方法 |
CN101226638B (zh) * | 2007-01-18 | 2010-05-19 | 中国科学院自动化研究所 | 一种对多相机***的标定方法及装置 |
CN101661617B (zh) * | 2008-08-30 | 2011-11-02 | 华为终端有限公司 | 摄像机标定的方法及装置 |
CN101551907B (zh) * | 2009-04-28 | 2011-05-04 | 浙江大学 | 一种多照相机自动化高精度标定方法 |
CN101561269B (zh) * | 2009-05-26 | 2011-08-17 | 空气动力学国家重点实验室 | 近景摄影测量的特征线自动匹配方法 |
CN101876555B (zh) * | 2009-11-04 | 2011-12-21 | 北京控制工程研究所 | 一种月球车双目视觉导航***标定方法 |
CN102012236A (zh) * | 2010-09-26 | 2011-04-13 | 郑州辰维科技股份有限公司 | 一种月球车双目视觉避障***标定方法 |
CN102901490B (zh) * | 2012-09-04 | 2014-08-06 | 北京信息科技大学 | 一种基于动态阈值的图像匹配方法及*** |
CN102901490A (zh) * | 2012-09-04 | 2013-01-30 | 北京信息科技大学 | 一种基于动态阈值的图像匹配方法及*** |
CN102865857A (zh) * | 2012-09-04 | 2013-01-09 | 北京信息科技大学 | 一种摄影测量图像匹配方法 |
CN103175544A (zh) * | 2013-03-04 | 2013-06-26 | 中国科学院长春光学精密机械与物理研究所 | 用地面控制点计算和标定航天tdiccd相机偏流角的方法 |
CN103175544B (zh) * | 2013-03-04 | 2015-09-09 | 中国科学院长春光学精密机械与物理研究所 | 用地面控制点计算和标定航天tdiccd相机偏流角的方法 |
CN104729532A (zh) * | 2015-03-02 | 2015-06-24 | 山东科技大学 | 一种全景相机的严密标定方法 |
CN104729532B (zh) * | 2015-03-02 | 2018-05-01 | 山东科技大学 | 一种全景相机的严密标定方法 |
CN105424058A (zh) * | 2015-11-06 | 2016-03-23 | 中国人民解放军信息工程大学 | 基于摄影测量技术的数码相机投影中心位置精确标定方法 |
CN105424058B (zh) * | 2015-11-06 | 2017-11-14 | 中国人民解放军信息工程大学 | 基于摄影测量技术的数码相机投影中心位置精确标定方法 |
CN106023142A (zh) * | 2016-05-06 | 2016-10-12 | 北京信息科技大学 | 利用共面直线阵列标定摄影测量相机的算法 |
CN106023142B (zh) * | 2016-05-06 | 2019-03-01 | 北京信息科技大学 | 利用共面直线阵列标定摄影测量相机的算法 |
CN109087339A (zh) * | 2018-06-13 | 2018-12-25 | 武汉朗视软件有限公司 | 一种激光扫描点与影像配准方法 |
CN112308894A (zh) * | 2019-08-01 | 2021-02-02 | 司法鉴定科学研究院 | 一种基于近景摄影测量的图像处理***及处理方法 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN105716542B (zh) | 一种基于柔性特征点的三维数据拼接方法 | |
CN1620153A (zh) | 利用平面控制点场进行非量测数码相机标定的方法 | |
CN109919911B (zh) | 基于多视角光度立体的移动三维重建方法 | |
CN103292710B (zh) | 一种应用双目视觉视差测距原理的距离测量方法 | |
CN103743352B (zh) | 一种基于多相机匹配的三维变形测量方法 | |
CN1294533C (zh) | 摄像机或照相机几何畸变的标定方法 | |
CN111192235A (zh) | 一种基于单目视觉模型和透视变换的图像测量方法 | |
WO1991015924A1 (en) | Computer assisted video surveying and method therefor | |
CN101487702A (zh) | 基于双目视觉的交通事故现场摄影测量方法 | |
CN115880344B (zh) | 一种双目立体匹配数据集视差真值获取方法 | |
CN112325788B (zh) | 基于数字照相的隧道内壁变形特征原位量测装置及方法 | |
CN111047586B (zh) | 一种基于机器视觉的像素当量测量方法 | |
CN101033950A (zh) | 岩石表面三维形貌的测量方法 | |
CN108550171B (zh) | 基于交比不变性的含有八卦编码信息的线阵相机标定方法 | |
CN109191527A (zh) | 一种基于最小化距离偏差的对位方法及装置 | |
CN111145269A (zh) | 一种鱼眼相机与单线激光雷达的外方位元素的标定方法 | |
CN109003308B (zh) | 一种基于相位编码的特殊成像范围相机标定***及方法 | |
CN105374067A (zh) | 一种基于pal相机的三维重建方法及其重建*** | |
CN111768451A (zh) | 一种基于移动显示屏的大尺寸双目视觉离焦标定方法 | |
CN108510546A (zh) | 一种适用于图谱及结构信息同步探测***的相机标定方法 | |
KR0184039B1 (ko) | 화질 검사 장치 및 그 화상 합성 방법 | |
CN102750698B (zh) | 纹理摄像机标定装置与方法及其纹理图像几何校正方法 | |
CN104167001A (zh) | 基于正交补偿的大视场摄像机标定方法 | |
CN114972447A (zh) | 一种基于无人机拍照的水体表面流迹线测量方法 | |
Schneider et al. | Geometric modelling and calibration of a high resolution panoramic camera |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C02 | Deemed withdrawal of patent application after publication (patent law 2001) | ||
WD01 | Invention patent application deemed withdrawn after publication |