MATLAB相机标定实战指南：从棋盘格到参数优化的完整工作流

计算机视觉领域的基础工作之一就是相机标定，它直接影响着后续三维重建、目标检测等任务的精度。作为行业标准工具，MATLAB Camera Calibrator提供了一套直观的图形化界面，但真正高效使用它需要掌握从数据采集到参数优化的全链路技巧。本文将带您深入每个关键环节，避开那些新手常踩的坑。

1. 标定前的准备工作：棋盘格与拍摄环境

棋盘格的质量直接影响标定结果的可靠性。建议使用激光打印的高精度棋盘格图案，粘贴在平整的刚性表面上。棋盘格尺寸的选择需要权衡：

• 小尺寸（如15mm方格）：适合近距离标定，但远距离拍摄时角点检测困难
• 大尺寸（如50mm方格）：远距离拍摄效果好，但需要更大的物理空间
• 推荐折中方案：25-30mm方格适用于大多数室内场景

拍摄时需要注意以下关键点：

1. 光照均匀性：避免反光和阴影，可使用漫反射光源
2. 拍摄角度多样性：建议覆盖以下典型角度组合：
   • 正视角度（3-5张）
   • ±30°倾斜（各3张）
   • ±45°倾斜（各2张）
   • 俯视/仰视（各2张）
3. 距离变化：从最近对焦距离到3米外分3-4个层次拍摄
4. 总数量控制：初学者建议20-30张，熟练后可减少到15-20张

常见错误：所有照片都在同一角度拍摄，这会导致某些畸变参数无法准确估计。

2. Camera Calibrator工具的核心参数解析

启动Camera Calibrator后，在Options面板中有几个关键选项需要理解其背后的数学含义：

为什么建议勾选3 Coefficients？
即使使用普通镜头，保留k₃参数也能更好地建模边缘畸变。现代计算资源充足，多一个参数不会显著增加计算负担，却可以提高标定鲁棒性。

勾选Skew的数学意义是考虑图像传感器像素可能不是严格正交的情况，其内参矩阵形式变为：

code复制[f_x  s  c_x
  0  f_y c_y
  0   0   1 ]

其中s即为skew参数。对于大多数消费级相机，这个值很小，但精确建模时应当保留。

3. 标定执行与误差诊断技巧

点击Calibrate按钮后，工具会显示两个关键可视化结果：

1. 重投影误差直方图：理想情况下应呈窄峰分布（<0.5像素）
2. 误差空间分布图：检查误差是否在图像边缘集中

图片剔除策略：

• 绝对阈值法：删除误差>2倍平均误差的图片
• 相对排序法：保留误差最小的15-20张
• 空间覆盖法：确保剩余图片仍能覆盖各种角度

实际操作中，可以结合MATLAB命令行进行更精细的控制：

matlab复制% 查看每张图片的误差
errors = cameraParams.ReprojectionErrors;
meanErrors = mean(sqrt(sum(errors.^2,2)),1);

% 找出误差最大的3张图片
[~,idx] = maxk(meanErrors,3);

% 交互式查看问题图片
figure;
imshow(imageFiles(idx(1)).name);
title(['问题图片: ' num2str(meanErrors(idx(1))) '像素误差']);

4. 参数导出与实际应用

导出相机参数后，关键数据结构包含：

matlab复制cameraParams.IntrinsicMatrix  % 内参矩阵
cameraParams.RadialDistortion % 径向畸变[k1 k2 k3]
cameraParams.TangentialDistortion % 切向畸变[p1 p2]

内参矩阵的实际应用示例：

假设我们要将图像坐标(u,v)转换为归一化相机坐标(x,y)：

matlab复制% 内参矩阵示例
K = [fx 0 cx; 0 fy cy; 0 0 1];

% 图像点（假设已经校正畸变）
uv_point = [u; v; 1];

% 转换为相机坐标
camera_point = K \ uv_point;  % 等价于inv(K)*uv_point
normalized_xy = camera_point(1:2)/camera_point(3);

对于无人机应用，还需要考虑镜头到IMU的变换。完整的传感器标定流程包括：

1. 相机内参标定（本文内容）
2. 相机-IMU外参标定
3. 时间同步校准
4. 多传感器联合优化

5. 高级技巧与疑难解答

标定结果不稳定的可能原因：

• 棋盘格平面度不足（使用亚克力板衬底）
• 角点检测不准确（尝试调整检测阈值）
• 拍摄角度分布不均（增加大角度样本）
• 镜头存在明显渐晕（考虑光照补偿）

标定后的验证方法：

1. 重投影测试：使用未参与标定的新图片
2. 三维重建测试：检查尺度一致性
3. 运动恢复结构(SfM)：验证多视图几何关系

MATLAB与OpenCV标定的差异：

当需要在不同平台间迁移参数时，需要注意这些细微差别。一个实用的转换函数示例：

matlab复制function ocam = matlabToOpenCV(cameraParams)
    ocam.K = cameraParams.IntrinsicMatrix';
    ocam.d = [cameraParams.RadialDistortion(1:2),...
              cameraParams.TangentialDistortion,...
              cameraParams.RadialDistortion(3)];
end

在实际项目中，我发现标定质量对后续算法性能的影响往往被低估。特别是在SLAM系统中，一个好的标定可以减少后端优化的负担。经过多次实践，我总结出一个经验法则：当重投影误差超过0.3像素时，就应该考虑重新标定或优化拍摄方案。