【相机标定】OpenCV 相机标定中的重投影误差与角点三维坐标计算详解

摘要：
本文将从以下几个方面展开，结合典型代码深入解析 OpenCV 中的相机标定过程，重点阐述重投影误差的计算方法与实际意义，并通过一个 calcBoardCornerPositions() 函数详细讲解棋盘格角点三维坐标的构建逻辑。

在计算机视觉领域，相机标定（Camera Calibration）是获取相机内参数和畸变参数的关键步骤。而重投影误差（Reprojection Error）则是衡量标定精度的重要指标。在使用 OpenCV 进行标定时，我们经常会接触到一个名为 computeReprojectionErrors() 的函数，它用于计算重投影误差，帮助我们评估标定结果的准确性。

↓↓↓↓↓↓ 以下正文 ↓↓↓↓↓↓

一、重投影误差是什么？

在相机标定中，我们使用一个带有已知几何尺寸的标定板（如棋盘格）进行拍摄，通过提取图像中的角点并与其在世界坐标系中的真实三维位置进行比较，来拟合相机的内参（焦距、主点）与畸变参数。

重投影误差定义为：

将三维点通过估计得到的相机内参、外参投影回图像平面后，与实际检测到的二维图像点之间的距离。

重投影误差越小，说明标定模型与实际相机系统越吻合。

在 OpenCV 中，这个误差通常用像素为单位来表示，计算结果常用于衡量整个标定质量的好坏。

OpenCV 中的计算方式

computeReprojectionErrors(objectPoints, rvecs, tvecs, reprojectionErrors);

objectPoints: 世界坐标系中的三维点数组
rvecs: 每一张标定图像的旋转向量
tvecs: 每一张标定图像的平移向量
reprojectionErrors: 输出每张图像的重投影误差

重投影误差一般多大算合理？

一般经验值如下：

小于 0.5 像素：标定质量非常好，适用于高精度需求。
0.5 ~ 1.0 像素：精度良好，适用于大部分应用。
1.0 ~ 2.0 像素：仍可接受，但存在优化空间。
大于 5 像素：标定可能存在错误或数据异常。

因此，如果你得到的误差值为 20.5 像素，就需要对标定流程和输入数据进行彻底排查。

二、重投影误差异常的原因分析

当我们遇到过大的重投影误差（比如 10 像素以上）时，很可能是以下几个方面出现了问题：

1. 棋盘格角点检测不准确

照片模糊、曝光不均、反光严重导致角点提取失败或偏移。
建议使用 drawChessboardCorners() 显示检测结果，手动验证角点匹配质量。

2. 标定板尺寸设置错误

如果你的标定板上一个格子的实际物理大小是 25mm，而你误填为 1mm 或 100mm，都会导致估计的相机矩阵出现异常。

3. 数据单位不统一

确保所有物理坐标（如角点位置）单位一致，且以米或毫米为准。

4. 输入图像质量差或数量太少

用于标定的图片最好在 10 张以上，覆盖不同视角与姿态。
图像需保证清晰、没有严重的畸变或遮挡。

5. 相机模型选择不当

OpenCV 支持多种畸变模型，如 CV_CALIB_RATIONAL_MODEL。错误的模型可能造成拟合能力下降。

三、异常相机矩阵的常见原因

相机矩阵的一般形式如下：

[ fx   0   cx ]
[  0  fy   cy ]
[  0   0    1 ]

fx, fy: 焦距，单位为像素
cx, cy: 图像主点（通常为图像中心）

如果你得到了一个如下的相机矩阵：

[42880.11, 0, 959.5;0, 42880.11, 539.5;0, 0, 1]

说明焦距异常大，很可能是以下问题导致的：

1. 标定板单位或格子大小设置错误

如果设定了极小的格子大小，比如 0.01（但实际是 10mm），那么焦距会被放大一千倍。

2. 视角变化不足

如果所有标定图片角度过于一致，优化过程无法正确拟合真实焦距。

3. 初始估计参数错误

有些标定代码会使用初始猜测值进行非线性优化。如果初始估计离实际值相差太远，会导致最终估计错误。

建议检查棋盘格实际物理尺寸，并可尝试使用 OpenCV 的 calibrateCamera() 函数中的 CALIB_USE_INTRINSIC_GUESS 标志手动输入初值。

四、棋盘格角点三维坐标构造

在相机标定时，我们需要构造真实世界中棋盘格角点的位置（objectPoints），这组点是在一个统一世界坐标系中的固定值，是整个标定过程的关键输入。

下面是一段典型的构造函数：

private void calcBoardCornerPositions(Mat corners) {final int cn = 3;float[] positions = new float[mCornersSize * cn];for (int i = 0; i < mPatternSize.height; i++) {for (int j = 0; j < mPatternSize.width * cn; j += cn) {positions[(int) (i * mPatternSize.width * cn + j + 0)] =(2 * (j / cn) + i % 2) * (float) mSquareSize;positions[(int) (i * mPatternSize.width * cn + j + 1)] =i * (float) mSquareSize;positions[(int) (i * mPatternSize.width * cn + j + 2)] = 0;}}corners.create(mCornersSize, 1, CvType.CV_32FC3);corners.put(0, 0, positions);
}