基于张正友与Tsai算法的摄像机标定系统

资 源 简 介

本项目在MATLAB环境下完整实现了两种主流的摄像机标定算法：张正友标定法（Zhang's Calibration Method）与Tsai两步标定法（Tsai's Two-Stage Calibration Method）。 针对张正友标定法，系统利用棋盘格标定板在不同位姿下的多幅图像，通过提取角点坐标并利用单应性矩阵约束关系，初步估算相机的内参数（主点坐标、焦距）和外参数（各图像对应的旋转矩阵与平移向量）。随后，系统采用极大似然准则建立目标函数，并利用勒文贝格-马夸特算法进行非线性全局优化，综合考虑径

详 情 说 明

基于MATLAB的摄像机标定两种经典算法实现系统

项目介绍

本项目是一个集成化的数学建模与仿真系统，旨在通过MATLAB语言复现计算机视觉领域最为经典的两种摄像机标定算法：张正友标定法（Zhang's Method）与Tsai两步标定法（Tsai's Two-Stage Method）。

系统涵盖了从原始数据生成、图像坐标模拟到最终模型优化的完整流程。通过对已知三维空间的物理点与其在二维图像平面上的投影点进行匹配，系统能够高精度地解算出相机的内部参数（包括焦距、主点、畸变系数）以及每一幅图像对应的外部参数（旋转矩阵和平移向量）。

功能特性

1. 全流程自动化：实现了从虚拟场景搭建、数据采集模拟、初值线性计算到非线性全局优化的自动化处理。
2. 算法多样化：同时提供了基于单应性矩阵的平面标定逻辑（张正友法）和基于径向排列约束的阶段式标定逻辑（Tsai法）。
3. 高保真模拟：数据生成环节集成了亚像素级别的随机噪声模拟以及径向畸变（k1, k2）注入，真实还原拍摄环境。
4. 可视化分析：内建误差评估模块，通过重投影误差趋势图、参数对比表等形式直观反映标定精度。
5. 稳健性能：利用勒文贝格-马夸特（LM）等高级优化算法，确保了在含噪声数据下解算的收敛度与准确性。

系统要求

1. 软件环境：MATLAB R2016b 或更高版本。
2. 必备工具箱：Optimization Toolbox（优化工具箱），用于执行非线性最小二乘优化。

实现逻辑与功能说明

1. 模拟数据生成逻辑

系统首先创建一个物理坐标系下的三维棋盘格点阵，默认规格为7x9个角点，每个方格步长为30mm。随后通过设定的相机内参真值（如1200px焦距）和径向畸变系数，对12组随机生成的位姿（旋转和平移）进行投影。在生成过程中，系统为每个2D像素坐标加入了0.05像素量级的正态随机噪声，以测试算法的鲁棒性。

2. 张正友标定法（Zhang's Method）实现步骤

这是本项目最核心的算法模块，遵循以下计算链条：

• 单应性矩阵计算：针对每张图像，利用直接线性变换（DLT）建立物理坐标与图像坐标的映射。
• 线性约束求解：利用旋转矩阵向量的正交性约束，构造线性方程组求解B矩阵，进而分解出相机的焦距、主点和倾斜因子。
• 外参初值估计：基于已获得的内参矩阵，快速推导每一幅图像的旋转向量（经罗德里格斯变换）和平移量。
• LM全局优化：使用勒文贝格-马夸特算法，以重投影误差最小化为目标，对内参、2个径向畸变系数以及所有外参进行联合优化。

3. Tsai两步标定法（Tsai's Method）实现步骤

针对该方法，系统采用了经典的两阶段策略：

• RAC约束阶段：基于径向排列约束原理，忽略由于畸变引起的径向位移方向上的变化，通过线性系统首先确定旋转矩阵的部分分量以及X、Y方向的平移。
• 迭代参数细化：在已知部分外参的基础上，通过非线性搜索（fminsearch）优化焦距、Z轴平移以及径向畸变系数。

4. 工具函数库详解

• 空间变换：系统内置了罗德里格斯（Rodrigues）变换函数，实现了旋转向量与旋转矩阵之间的相互转换，这是非线性优化的关键。
• 残差计算：专门编写了适用于不同算法的残差评估函数，通过计算投影点与实测点之间的欧氏距离来驱动优化器运行。
• 单应性提取：通过SVD奇点值分解处理超定方程组，确保从噪点数据中提取最稳健的单应性特征。

5. 结果评估与可视化逻辑

• 参数对比：系统将估计出的内参矩阵与模拟时使用的真值进行矩阵化对比，并输出畸变系数的误差。
• 重投影误差计算：标定完成后，系统将三维点重新投影至各个像平面，计算其与含有噪声的“观察值”之间的误差。
• 图形展示：自动生成重投影误差的柱状图，清晰展示每一张图像的平均误差水平，帮助用户判断标定的稳定性。

使用方法

1. 打开MATLAB，将当前工作目录定位至项目所在文件夹。
2. 在命令行窗口直接运行主入口函数。
3. 系统将先后执行模拟数据准备、张正友法计算和Tsai法计算。
4. 运行结束后，MATLAB将自动弹出评估图表，并在命令行输出两种算法的内参、畸变系数估计值与真值的对比结果。
5. 用户可以根据输出的重投影平均误差（单位：pixels）来判断当前模拟环境下的标定精度。