【视频开发】【计算机视觉】相机标定（Camera calibration）原理、步骤

相机标定（Camera calibration）原理、步骤

---

author@jason_ql（lql0716） 
http://blog.csdn.net/lql0716

---

• 在图像测量过程以及机器视觉应用中，为肯定空间物体表面某点的三维几何位置与其在图像中对应点之间的相互关系，必须创建相机成像的几何模型，这些几何模型参数就是相机参数。在大多数条件下这些参数必须经过实验与计算才能获得，这个求解参数的过程就称之为相机标定（或摄像机标定）。不管是在图像测量或者机器视觉应用中，相机参数的标定都是很是关键的环节，其标定结果的精度及算法的稳定性直接影响相机工做产生结果的准确性。所以，作好相机标定是作好后续工做的前提，提升标定精度是科研工做的重点所在。

• 经常使用术语 
  内参矩阵: Intrinsic Matrix 
  焦距: Focal Length 
  主点: Principal Point 
  径向畸变: Radial Distortion 
  切向畸变: Tangential Distortion 
  旋转矩阵: Rotation Matrices 
  平移向量: Translation Vectors 
  平均重投影偏差: Mean Reprojection Error 
  重投影偏差: Reprojection Errors 
  重投影点: Reprojected Points

一、坐标系的转换

1.1 世界坐标系

世界坐标系（world coordinate） (xw,yw,zw) ，也称为测量坐标系，是一个三维直角坐标系，以其为基准能够描述相机和待测物体的空间位置。世界坐标系的位置能够根据实际状况自由肯定。

1.2 相机坐标系

相机坐标系（camera coordinate） (xc,yc,zc) ，也是一个三维直角坐标系，原点位于镜头光心处，x、y轴分别与相面的两边平行，z轴为镜头光轴，与像平面垂直。

1.3 世界坐标系转换为相机坐标系

⎡⎣⎢⎢⎢xcyczc1⎤⎦⎥⎥⎥=[R0t1]⎡⎣⎢⎢⎢xwywzw1⎤⎦⎥⎥⎥

其中 R 为3*3的旋转矩阵， t 为3*1的平移矢量， (xc,yc,zc,1)T 为相机坐标系的齐次坐标， (xw,yw,zw,1)T 为世界坐标系的齐次坐标。

1.4 像素坐标系、图像坐标系

• 图1： 
  

• 像素坐标系（pixel coordinate） 
  如图1，像素坐标系 uov 是一个二维直角坐标系，反映了相机CCD/CMOS芯片中像素的排列状况。原点 o 位于图像的左上角， u 轴、 v 轴分别于像面的两边平行。像素坐标系中坐标轴的单位是像素（整数）。

  像素坐标系不利于坐标变换，所以须要创建图像坐标系 XOY ，其坐标轴的单位一般为毫米（mm），原点是相机光轴与相面的交点（称为主点），即图像的中心点， X 轴、 Y 轴分别与 u 轴、 v 轴平行。故两个坐标系实际是平移关系，便可以经过平移就可获得。

• 图像坐标系转换为像素坐标系

⎡⎣⎢uv1⎤⎦⎥=⎡⎣⎢1/dX0001/dY0u0v01⎤⎦⎥⎡⎣⎢XY1⎤⎦⎥

其中， dX 、 dY 分别为像素在 X 、 Y 轴方向上的物理尺寸， u0,v0 为主点（图像原点）坐标。

1.5 针孔成像原理

图2： 

如图，空间任意一点 P 与其图像点 p 之间的关系， P 与相机光心 o 的连线为 oP ， oP 与像面的交点 p 即为空间点 P 在图像平面上的投影。 
该过程为透视投影，以下矩阵表示：

s⎡⎣⎢XY1⎤⎦⎥=⎡⎣⎢f000f0001000⎤⎦⎥⎡⎣⎢⎢⎢xyz1⎤⎦⎥⎥⎥

其中， s 为比例因子（ s 不为0）， f 为有效焦距（光心到图像平面的距离）， (x,y,z,1)T 是空间点 P 在相机坐标系 oxyz 中的齐次坐标， (X,Y,1)T 是像点 p 在图像坐标系 OXY 中的齐次坐标。

1.6 世界坐标系转换为像素坐标系

s⎡⎣⎢uv1⎤⎦⎥=⎡⎣⎢1/dX0001/dY0u0v01⎤⎦⎥⎡⎣⎢f000f0001000⎤⎦⎥[R0t1]⎡⎣⎢⎢⎢xwywzw1⎤⎦⎥⎥⎥=⎡⎣⎢αx000αy0u0v01000⎤⎦⎥[R0t1]⎡⎣⎢⎢⎢xwywzw1⎤⎦⎥⎥⎥=M1M2Xw=MXw

其中， αx=f/dX 、 αy=f/dY ，称为 u 、 v 轴的尺度因子， M1 称为相机的内部参数矩阵， M2 称为相机的外部参数矩阵， M 称为投影矩阵。

2 相机内参与畸变参数

2.1 相机内参

• 参看1.6节所述

2.2 畸变参数

• 畸变参数(distortion parameters)

  畸变（distortion）是对直线投影（rectilinear projection）的一种偏移。

  简单来讲直线投影是场景内的一条直线投影到图片上也保持为一条直线。

  畸变简单来讲就是一条直线投影到图片上不能保持为一条直线了，这是一种光学畸变（optical aberration）,可能因为摄像机镜头的缘由。

• 畸变通常能够分为：径向畸变、切向畸变

  一、径向畸变来自于透镜形状 
  二、切向畸变来自于整个摄像机的组装过程

  畸变还有其余类型的畸变，可是没有径向畸变、切向畸变显著

• 畸变图示 
  

• 径向畸变

  实际摄像机的透镜老是在成像仪的边缘产生显著的畸变，这种现象来源于“筒形”或“鱼眼”的影响。

  以下图，光线在原理透镜中心的地方比靠近中心的地方更加弯曲。对于经常使用的普通透镜来讲，这种现象更加严重。筒形畸变在便宜的网络摄像机中很是厉害，但在高端摄像机中不明显，由于这些透镜系统作了不少消除径向畸变的工做。 
  

  对于径向畸变，成像仪中心（光学中心）的畸变为0，随着向边缘移动，畸变愈来愈严重。

• 径向畸变包括：枕形畸变、桶形畸变

• 切向畸变

  • 切向畸变是因为透镜制造上的缺陷使得透镜自己与图像平面不平行而产生的。

  • 切向畸变可分为：薄透镜畸变、离心畸变

  • 切向畸变图示： 
    

2.2.1 opencv中的畸变模型

• 径向畸变模型：

如下公式由泰勒展式得出，在opencv中 K=1 ， r2=x2+y2 ,  (x,y) 为真实坐标（发生畸变）， (x′,y′) 为理想坐标。

δxr=x(k1r2+k2r4+k3r6+K)

δyr=y(k1r2+k2r4+k3r6+K)

• 切向畸变模型：

δxd=2p1xy+p2(r2+2x2)+K

δyd=2p1(r2+2y2)+2p2xy+K

• 理想坐标 (x′,y′) 与真实坐标 (x,y) ：

x′=x+δxr+δxd

y′=y+δyr+δyd

也即：

[x′y′]=(1+k1r2+k2r4+k3r6)[xy]+[2p1xy+p2(r2+2x2)2p1(r2+2y2)+2p2xy]

实际计算过程当中，若是考虑太多高阶的畸变参数，会致使标定求解的不稳定。

2.2.2 张氏标定中的畸变模型

• 张正友的方法只考虑了径向畸变，没有考虑切向畸变

• 模型：

  • 理想状况下(没有畸变)图片的像素坐标为：

    (u,v)

  • 真实的像素坐标为：

    (u¯,v¯)

  • 真实坐标与理想坐标的关系式（泰勒展开）： 
    (x,y)  and  (x¯,y¯)  are the ideal (distortion-free) and real (distorted) normalized image coordinates. 
    

    x¯=x+x[k1(x2+y2)+k2(x2+y2)2]

  y¯=y+y[k1(x2+y2)+k2(x2+y2)2]

  其中 k1 ， k2 是径向畸变系数(coefficients of the radial distortion)。

  像素坐标表示为： 
  

  u¯=u+(u−u0)[k1(x2+y2)+k2(x2+y2)2]

  v¯=v+(v−v0)[k1(x2+y2)+k2(x2+y2)2]

  即： 
  

  [(u−u0)(x2+y2)(v−v0)(x2+y2)(u−u0)(x2+y2)2(v−v0)(x2+y2)][k1k2]=[u¯−uv¯−v]

2.2.3 Tasi畸变模型

• Tasi畸变模型同张氏畸变模型相似，也只考虑了径向畸变

2.3 相机标定

• 相机标定步骤： 
  

  一、打印一张棋盘格，把它贴在一个平面上，做为标定物。 
  二、经过调整标定物或摄像机的方向，为标定物拍摄一些不一样方向的照片。 
  三、从照片中提取棋盘格角点。 
  四、估算理想无畸变的状况下，五个内参和六个外参。 
  五、应用最小二乘法估算实际存在径向畸变下的畸变系数。 
  六、极大似然法，优化估计，提高估计精度。

3 Matlab相机标定

• 打开matlab，找到“Camera Calibrator”并打开

• 在新窗口中选择添加图片“Add Images”

• 添加图片以后，会有以下提示，设置棋盘格的实际大小以后，点击ok

• 点击“Calibrate”开始计算

• 点击“Export Camera Parameters”，输出到matlab命令窗口

• 导出的数据以下：

  
  
  
  

   
   
   
   
1
   
   
   
   
2
   
   
   
   
3
   
   
   
   
4
   
   
   
   
5
   
   
   
   
6
   
   
   
   
7
   
   
   
   
8
   
   
   
   
9
   
   
   
   
10
   
   
   
   
11
   
   
   
   
12
   
   
   
   
13
   
   
   
   
14
   
   
   
   
15
   
   
   
   
16
   
   
   
   
17
   
   
   
   
18
   
   
   
   
19
   
   
   
   
20
   
   
   
   
21
   
   
   
   
22
   
   
   
   
23
   
   
   
   
24
   
   
   
   
25
   
   
   
   
26
   
   
   
   
27
   
   
   
   
28
   
   
   
   
29
   
   
   
   
30
   
   
   
   
31
   
   
   
   
32
   
   
   
   
33
   
   
   
   
34
   
   
   
   
35
   
   
   
   
36
   
   
   
   
37
   
   
   
   
38
  
  
  
  
cameraParams = 

  cameraParameters (具备属性):

   Camera Intrinsics
    IntrinsicMatrix: [3x3 double]
        FocalLength: [510.6720 511.9534]
     PrincipalPoint: [332.3800 238.4849]
               Skew: 0

   Lens Distortion
        RadialDistortion: [0.0618 -0.4877]
    TangentialDistortion: [0 0]

   Camera Extrinsics
      RotationMatrices: [3x3x50 double]
    TranslationVectors: [50x3 double]

   Accuracy of Estimation
    MeanReprojectionError: 0.1923
       ReprojectionErrors: [77x2x50 double]
        ReprojectedPoints: [77x2x50 double]

   Calibration Settings
                        NumPatterns: 50
                        WorldPoints: [77x2 double]
                         WorldUnits: 'mm'
                       EstimateSkew: 0
    NumRadialDistortionCoefficients: 2
       EstimateTangentialDistortion: 0


estimationErrors = 

  cameraCalibrationErrors (具备属性):

    IntrinsicsErrors: [1x1 intrinsicsEstimationErrors]
    ExtrinsicsErrors: [1x1 extrinsicsEstimationErrors]

• 在命令窗口输入如下命令便可得到内参矩阵和径向畸变：

  
  
  
  

   
   
   
   
1
   
   
   
   
2
   
   
   
   
3
   
   
   
   
4
   
   
   
   
5
   
   
   
   
6
   
   
   
   
7
   
   
   
   
8
   
   
   
   
9
   
   
   
   
10
   
   
   
   
11
   
   
   
   
12
   
   
   
   
13
  
  
  
  
>> cameraParams.IntrinsicMatrix

ans =

  510.6720         0         0
         0  511.9534         0
  332.3800  238.4849    1.0000

>> cameraParams.RadialDistortion

ans =

    0.0618   -0.4877

参考资料

• 相关文章

  • 机器学习、深度学习、计算机视觉、天然语言处理及应用案例——干货分享（持续更新……）
  • 旋转矩阵、欧拉角、四元数理论及其转换关系

• 基于生长的棋盘格角点检测方法–（1）原理介绍

• opencv角点检测、棋盘格检测、亚像素cvFindCornerSubPix（）

• 相机标定(Camera calibration)

• OpenCV相机标定和姿态更新

• Camera Calibration and 3D Reconstruction

• Camera Calibration

• 相机姿态估计（二）–单目POSIT算法

• 相机位姿估计1：根据四个特征点估计相机姿态

• Head Pose Estimation using OpenCV and Dlib

• opencv角点检测、棋盘格检测、亚像素cvFindCornerSubPix（）

• 关于OpenCV的那些事——相机标定

• 计算机视觉-相机内参数和外参数

• 张正友相机标定算法 【2】

• Python+OpenCV学习（17）—摄像机标定

• 相机标定Opencv实现以及标定流程&&标定结果评价&&图像矫正流程解析（附标定程序和棋盘图）

• 张正友相机标定Opencv实现以及标定流程&&标定结果评价&&图像矫正流程解析（附标定程序和棋盘图）

• MATLAB–相机标定教程

• 《学习opencv》

• 张正友相机标定论文《A Flexible New Technique for Camera Calibration》

• OpenCV相机标定

• 畸变模型

• 相机的那些事儿 - 概念、模型及标定

• 相机标定（二）之相机畸变模型