一种基于分段模型的摄像机标定方法技术

一种基于分段模型的摄像机标定方法属于摄像机标定方法领域，该标定方法计算标定模型参数的图像特征点分别在像素平面中划分的三个区域内，标定模型是分段函数；其用畸变较少的圆形区域内的特征点计算线性模型，求解并在后续标定模型中延用这些准确性高的内部参数矩阵A和外部参数R、t；仅需对圆环区域模型中的畸变系数求解，并对畸变系数的优化初始值进行求解，降低了待求解的维数，避免了计算结果的不收敛情况并大幅缩短了计算时间；其将第三区域所建立的畸变模型转化为连续的样条平滑函数以修正第三区域内的全部像素点畸变，使整幅图像中距图像中线较远位置的特殊畸变都能够被精确修正，大幅度地提高了标定精度。

【技术实现步骤摘要】

本专利技术属于摄像机标定方法领域，具体涉及。
技术介绍
摄像机标定是由已知特征点的像素坐标和世界坐标去求解标定模型参数的过程，从而建立图像像素位置与三维场景点坐标之间的对应投影关系。闻精度摄像机标定是数字图像测量和二维重构技术的基础，摄像机两步标定法是目前传统标定方法中应用最为广泛的方法，具有较高的标定精度。当应用场合所要求的精度很高且摄像机的参数不经常变化时，传统标定方法中的两步标定法通常会作为首选方法。两步标定法的第一步是先利用直接线性变换方法求解摄像机参数，第二步考虑畸变因素，将第一步中求得的参数作为初始值，利用优化算法重新对标定模型中的参数进行非线性求解，从而获得更准确的标定模型。目如，两步标定法中有代表性的是Tsai, Heikkila和张正友的标定方法,这二种方法都米用了统一形式的畸变函数来描述整个镜头的畸变，其中Tsai和张正友的模型中考虑了径向畸变的影响，而Heikkila在此基础上又增加了切向畸变。但是，由于镜头制造误差存在随机性，各种畸变误差在图像中的分布规律并不完全一致，因此，在实际应用中很难找到能够精准体现图像扭曲程度的畸变模型。如果想更准确和更详尽的描述镜头畸变，模型中的参数就会增加，复杂的模型反而会引起数值计算的不稳定性。基于上述因素，现有的标定方法只能选择一些具体的畸变形式。也就是说，图像中会存在一些畸变类型，并未被选定的畸变模型所体现，最终会影响摄像机标定的精度。而且，现有的标定方法中待求解的维数往往较大，导致在优化过程中容易出现结果不收敛或计算时间过长的情况。此外，现有方法中的畸变模型都试图利用一个统一的函数来描述整幅图像中的畸变，最终，标定模型中的参数是利用图像中所有特征点计算获得的。然而，位于距离图像中心较远的区域中的点通常带有不确定型畸变，这些点通常位于距离图像中心较远的区域中，在优化求解过程中会将大量误差带入模型中，最终降低和损害整幅图像的标定精度。Tsai R Y在《A versatile camera calibration technique for high-accuracy3D machine vision metrology using off-the-shelf TV cameras and lenses》[J].IEEEJournal of Robotics and Automation，1987，3 (4): 323-344 中提出了一种摄像机标定法使用的转换坐标系，即:建立世界坐标系、摄像机坐标系、图像物理坐标系和图像像素平面坐标系。张正友在Zhang Z.(2000).A flexible new technique forcamera calibration.1EEE Trans, on Pattern Analysis and MachineIntelligence, 22 (11): 1330-1334.提出了一种相机标定方法，其中包括一种具体的线性模型求解和标定方法。
技术实现思路
为了解决现有摄像机两步标定法标定模型中的参数都是利用图像中所有特征点计算获得并都是采用一个统一的函数来描述整幅图像中的畸变，图像中会存在一些畸变类型并未被选定的畸变模型所体现，以及现有的标定方法虽然大都能获得较好的优化初值，但是由于待求解的维数较大，在优化过程中的大量非线性计算容易出现计算结果不收敛或计算时间过长的情况的技术问题，本专利技术提供。本专利技术解决技术 问题所采取的技术方案如下:，包括如下步骤:步骤一:固定摄像机并将棋盘格标定板的中心置于摄像机光轴附近，使棋盘格标定板在一个相对固定的成像位置围绕摄像机镜头变换不同的倾角和姿态，用摄像机对标定板的每一种姿态拍摄一幅照片图像，最终取得多幅照片图像；步骤二:按照Tsai R Y的坐标转换方法确立转换坐标系，即:建立世界坐标系、摄像机坐标系、图像物理坐标系和图像像素平面坐标系，世界坐标系内的三维坐标点坐标(Xw, Yw, Zw)被依次转换为摄像机点坐标(X。，Y。, Z。)、理想图像点坐标(xu, yu)、真实图像点坐标( yd)和像素点坐标(xp，yp);另外，参照张正友的摄像机标定方法将世界坐标系的Zdi坐标设为O,即Zw=O ;此时,世界坐标系内任意点的坐标表示为(xw, Yw, O),或仅表示为世界坐标系XY轴平面的二维点坐标(Xw, Yw)；步骤三:利用亚像素角点检测方法分别对步骤一所取得的各幅照片图像提取对应的角点像素坐标；步骤四:将步骤二所述的图像像素坐标平面划分成三个相邻的区域:圆形区域、圆环形区域和第三区域；步骤五:将步骤三所取得的各幅照片图像的角点像素坐标按照其在步骤四所属的三个不同区域分类，分别得到属于圆形区域内的图像角点像素检测坐标、属于圆环形区域内的图像角点像素检测坐标和属于第三区域内的图像角点像素检测坐标；步骤六:按照步骤二所建立的转换坐标系在步骤四所取得的圆形区域建立线性标定模型，将棋盘格标定板上方格角点在世界坐标系对应的三维点坐标(Xw，Yw，0)转换为属于圆形区域内的图像像素坐标(Xp，yp)，忽略畸变的影响，所建立的圆形区域内对应平面的标定线性模型如下:本文档来自技高网...

【技术保护点】
一种基于分段模型的摄像机标定方法，其特征在于，该方法包括如下步骤：步骤一：固定摄像机并将棋盘格标定板的中心置于摄像机光轴附近，使棋盘格标定板在一个相对固定的成像位置围绕摄像机镜头变换不同的倾角和姿态，用摄像机对标定板的每一种姿态拍摄一幅照片图像，最终取得多幅照片图像；步骤二：按照Tsai?R?Y的坐标转换方法确立转换坐标系，即：建立世界坐标系、摄像机坐标系、图像物理坐标系和图像像素平面坐标系，世界坐标系内的三维坐标点坐标(Xw,Yw,Zw)被依次转换为摄像机点坐标(Xc,Yc,Zc)、理想图像点坐标(xu,yu)、真实图像点坐标(xd,yd)和像素点坐标(xp,yp)；另外，参照张正友的摄像机标定方法将世界坐标系的Zw轴坐标设为0，即Zw=0；此时，世界坐标系内任意点的坐标表示为(Xw,Yw,0)，或仅表示为世界坐标系XY轴平面的二维点坐标(Xw,Yw)；步骤三：利用亚像素角点检测方法分别对步骤一所取得的各幅照片图像提取对应的角点像素坐标；步骤四：将步骤二所述的图像像素坐标平面划分成三个相邻的区域：圆形区域（1）、圆环形区域（2）和第三区域（3）；步骤五：将步骤三所取得的各幅照片图像的角点像素坐标按照其在步骤四所属的三个不同区域分类，分别得到属于圆形区域（1）内的图像角点像素检测坐标、属于圆环形区域（2）内的图像角点像素检测坐标和属于第三区域（3）内的图像角点像素检测坐标；步骤六：按照步骤二所建立的转换坐标系在步骤四所取得的圆形区域（1）建立线性标定模型，将棋盘格标定板上方格角点在世界坐标系对应的三维点坐 标(Xw,Yw,0)转换为属于圆形区域（1）内的图像像素坐标(xp,yp)，忽略畸变的影响，所建立的圆形区域（1）内对应平面的标定线性模型如下：sxpyp1=Ar1r2r3tXwYw01=A[R,t]XwYw01...(1)当世界坐标系的Zw轴坐标设为0时，即世界坐标系内任意点的坐标表示为(Xw,Yw)时，式（1）可变换为下式sxpyp1=Ar1r2tXwYw1=HXwYw1...(1.1)式（1）和式（1.1）中，s表示比例因子，R=[r1,r2,r3]是世界坐标的旋转矩阵，ri（i=1,2,3）表示旋转矩阵R的第i列，t表示平移向量，R和t统称为标定模型的外部参数；A=αγu00βv0001是内部参数矩阵，H表示3×3的单应性矩阵；α和β表示像素平面中U轴和V轴的比例因子，γ表示像素平面两坐标轴的不垂直因子，(u0,v0)表示相机光轴与图像平面的交点在像素坐标系上的坐标；步骤七：利用步骤五所述圆形区域（1）内的角点像素检测坐标(xp,yp)以及这些角点在步骤二所述世界坐标系中对应的二维点坐标值(Xw,Yw)，并使用张正友的相机标定方法中线性模型求解方法，求解其对应在步骤六中所述的平面标定线性模型，完成对步骤四所述圆形区域（1）的线性标定；在上述求解过程中，同时求得步骤六所述的内部参数A，并且同时求得各组步骤一所述每幅照片图像对应在步骤六所述的外部参数R和t；步骤八：利用步骤七求出的内部参数A以及外部参数R和t，对步骤四所述属于圆环形区域（2）建立非线性畸变模型如下：sxuyu1=r1r2tXwYw1...(2)xuyu=(1+k1r2+k2r4)xdyd+2p1xdyd+p2(r2+2xd2)p1(r2+2yd2)+2p2xdyd...(3)xpyp1=αγu00βv0001xdyd1...(4)式（3）中的r=xu2+yu2;k1,k2,p1,p2表示图像物理坐标系中径向与切向畸变函数的系数；步骤九：求解步骤八中的畸变系数k1,k2,p1,p2，并根据已求得的畸变系数k1,k2,p1,p2代入步骤八所述的由式（2）、式（3）和式（4）构成的整体非线性畸变模型，以完成对步骤四所述圆环形区域（2）的标定；步骤十：获取步骤四所述属于第三区域（3）内的图像的角点的理想像素坐标(ud,vd)；将步骤一所取得的每一幅照片图像中属于第三区域（3）内的图像的角点在世界坐标系中对应的二维点坐标值(Xw,Yw)作为转换目标，利用在步骤七获得的圆形区域（1）内对应平面的标定线性模型的内部参数矩阵A和每幅照片图像对应的外部参数R和t对前述的二维点坐标值(Xw,Yw)投影到像素坐标系平面，以获得它们对应的理想参考像素坐标(ud,vd)；步骤十一：将步骤十所述理想参考像素坐标(ud,vd)与步骤五属于第三区域（3）内的图像角点像素检测坐标(up,vp)作差比较，获得属于第三区域...

【技术特征摘要】
1.一种基于分段模型的摄像机标定方法，其特征在于，该方法包括如下步骤: 步骤一:固定摄像机并将棋盘格标定板的中心置于摄像机光轴附近，使棋盘格标定板在一个相对固定的成像位置围绕摄像机镜头变换不同的倾角和姿态，用摄像机对标定板的每一种姿态拍摄一幅照片图像，最终取得多幅照片图像； 步骤二:按照Tsai R Y的坐标转换方法确立转换坐标系，即:建立世界坐标系、摄像机坐标系、图像物理坐标系和图像像素平面坐标系，世界坐标系内的三维坐标点坐标(Xw, Yw, Zw)被依次转换为摄像机点坐标(X。，Y。, Z。)、理想图像点坐标(xu, yu)、真实图像点坐标( yd)和像素点坐标(xp，yp);另外，参照张正友的摄像机标定方法将世界坐标系的Zdi坐标设为O,即Zw=O ;此时,世界坐标系内任意点的坐标表示为(xw, Yw, O),或仅表示为世界坐标系XY轴平面的二维点坐标(Xw, Yw)； 步骤三:利用亚像素角点检测方法分别对步骤一所取得的各幅照片图像提取对应的角点像素坐标； 步骤四:将步骤二所述的图像像素坐标平面划分成三个相邻的区域:圆形区域(I)、圆环形区域(2)和第三区域(3); 步骤五:将步骤三所取得的各幅照片图像的角点像素坐标按照其在步骤四所属的三个不同区域分类，分别得到属于圆形区域(I)内...

【专利技术属性】
技术研发人员：孙秋成，周晓东，李纯净，
申请(专利权)人：长春工业大学，
类型：发明
国别省市：