2、光学基础：透镜成像原理、理想成像模型、畸变产生的物理根源
光学畸变测量的一切方法论，都建立在对“光如何通过透镜形成图像”这一基本过程的深刻理解之上。若将畸变测量比作诊断疾病，那么本章要解决的便是“健康人体的生理构造”——只有先弄清理想成像的数学描述与物理边界，才能准确识别病理性的像差偏移。本章从几何光学的核心定律出发，逐步构建理想成像模型，并深入剖析畸变区别于其他像差的独特物理根源，为后续的数学模型标定与测量方案设计奠定不可动摇的理论基石。
2.1 透镜成像的核心原理：从费马原理到高斯光学
透镜成像的物理本质，是光波在穿越不同折射率介质时，遵循费马原理（Fermat's Principle）——光总是沿着光程取极值的路径传播。当一束发自物点P的光线经过透镜折射后，若所有光线都能汇聚于同一点P'，则称P'为P的像点。这一过程在数学上可由折射定律（Snell's Law）逐点追踪，但在工程应用中，我们更常采用近轴近似下的高斯光学（Gaussian Optics）来建立简洁的成像关系。
高斯光学将透镜视为理想薄透镜，假设光线与光轴的夹角极小（sinθ≈θ），且透镜厚度可忽略。在此条件下，物距l、像距l'与焦距f之间满足高斯公式：
1/f = 1/l' - 1/l其中符号约定遵循笛卡尔坐标系：光线从左向右传播，物距l为负值（实物在透镜左侧），像距l'为正值（实像在透镜右侧）。横向放大率β定义为像高y'与物高y之比：
β = y'/y = l'/l这一公式揭示了理想成像的核心特征：同一物平面上的所有点，其像点均位于同一像平面，且放大率恒定。这意味着一个正方形物体，其像必然是完美的正方形——这正是畸变测量的“零基准”。
2.2 理想成像模型的数学框架与坐标映射
在实际光学系统中，我们通常用矩阵光学或光线追迹来描述成像过程。对于理想成像模型，物点坐标(x, y, z)与像点坐标(x', y', z')之间的映射关系可简化为：
x' = β·x
y' = β·y
z' = β²·z（纵向放大率）在二维图像传感器平面上，我们只关心横向映射。若以图像中心为原点，理想像高r₀与物高r的关系为：
r₀ = β·r其中r = √(x² + y²)为物点的径向距离。这一线性关系是畸变测量的基准线——实际光学系统由于像差的存在，实际像高r'往往偏离r₀，偏离量Δr = r' - r₀即为畸变量。
为了工程化描述，我们引入归一化视场坐标。设最大像高为r_max，归一化视场ρ = r/r_max（0≤ρ≤1）。理想像高可写为：
r₀(ρ) = β·ρ·r_max而实际像高通常可用奇次多项式展开（因畸变具有旋转对称性，偶次项相互抵消）：
r'(ρ) = k₁·ρ + k₃·ρ³ + k₅·ρ⁵ + k₇·ρ⁷ + ...其中k₁为线性放大系数（理想情况下k₁=β·r_max），k₃、k₅、k₇为畸变系数。这一多项式模型是后续畸变标定的数学基础，也是工业相机标定中Brown模型的核心构成。
| 模型参数 | 物理含义 | 理想值 | 实际典型值 |
|---|---|---|---|
| k₁ | 线性放大系数 | β·r_max | 与设计值偏差<0.5% |
| k₃ | 三阶畸变系数 | 0 | 10⁻⁴ ~ 10⁻² 量级 |
| k₅ | 五阶畸变系数 | 0 | 10⁻⁶ ~ 10⁻⁴ 量级 |
| k₇ | 七阶畸变系数 | 0 | 通常可忽略 |
2.3 畸变产生的物理根源：主光线与光阑的博弈
畸变与其他像差（球差、彗差、像散等）的本质区别在于：畸变不影响成像的清晰度，只改变像的几何形状。其物理根源在于不同视场的主光线（Chief Ray，即通过孔径光阑中心的光线）经过光学系统后，其像平面交点高度与理想像高不成正比。
具体而言，畸变的大小取决于两个关键因素：
- 孔径光阑（Aperture Stop）的位置：光阑位于透镜前方或后方，会显著改变主光线的入射高度，从而影响畸变的符号与大小。例如，单透镜加光阑位于前方时，通常产生负畸变（桶形）；光阑位于后方时，产生正畸变（枕形）。
- 透镜球面折射的非线性：即使对于薄透镜，大视场角下光线在球面上的入射角较大，折射定律的线性近似失效，导致实际像高偏离理想值。这种偏离随视场角的增大而加速增长，表现为高次畸变项。
从像差理论的角度，畸变属于赛德尔像差（Seidel Aberrations）中的一种，其波前表达式为：
W_distortion = S_V · r³ · cosθ其中S_V为第五赛德尔和，r为归一化孔径坐标，θ为方位角。值得注意的是，畸变的波前误差与视场坐标的三次方成正比，这意味着视场越大，畸变增长越剧烈。
2.4 畸变的分类与工程表征
根据畸变的方向性，可将其分为两大类：
- 正畸变（枕形畸变，Pincushion Distortion）：实际像高大于理想像高，图像边缘向外膨胀，直线向外弯曲。常见于长焦镜头、望远系统。
- 负畸变（桶形畸变，Barrel Distortion）：实际像高小于理想像高，图像边缘向内收缩，直线向内弯曲。常见于广角镜头、鱼眼镜头、监控摄像头。
在工程测量中，畸变量通常用相对畸变（Relative Distortion）表示：
D = (Y' - Y₀) / Y₀ × 100%其中Y'为实际像高，Y₀为理想像高。工业镜头通常要求D<1%，高精度测量镜头要求D<0.1%，而机器视觉测量系统往往需要D<0.01%才能保证亚像素精度。
2.5 案例拆解：广角镜头畸变对测量精度的影响
背景与约束：某自动化产线需要安装一台视觉测量系统，用于检测直径200mm工件的圆度。选用焦距8mm的广角工业镜头，搭配1/1.8" CMOS传感器（像面尺寸7.2mm×5.4mm），工作距离300mm。系统要求测量重复性±0.02mm。
初始现象：在系统调试阶段，工程师发现当工件位于视场边缘时，测量直径比中心位置小约0.35mm，且边缘区域的直线边缘呈现明显的桶形弯曲。
关键数据分析：
镜头标称畸变：-2.5%（边缘视场）
传感器半对角线：√(3.6²+2.7²) = 4.5mm
工件半径100mm，在300mm工作距离下，物方视场角：
半视场角 = arctan(100/300) ≈ 18.4°
对应像高 = 8mm × tan(18.4°) ≈ 2.66mm
归一化视场 ρ = 2.66/4.5 ≈ 0.59
查镜头畸变曲线，ρ=0.59处畸变约-0.8%分析判断：在ρ=0.59视场处，-0.8%的畸变意味着实际像高比理想像高小0.8%。对于工件半径100mm，边缘测量误差为100mm×0.8% = 0.8mm，远超±0.02mm的重复性要求。若不校正，系统完全无法满足测量精度。
实际操作步骤：
- 使用标准棋盘格标定板（格距10mm，精度±0.001mm），在同样工作距离下采集15张不同姿态的图像。
- 利用OpenCV的calibrateCamera函数求解内参，获得畸变系数k₁、k₃、k₅。
- 对实时图像应用undistort操作，将像素坐标映射到理想坐标。
- 在校正后的图像上重新测量工件直径。
结果对比：校正前边缘测量误差-0.8mm，校正后误差降至+0.008mm，满足±0.02mm要求。但代价是边缘视场图像被重采样，有效分辨率下降约15%。
复盘经验：广角镜头的畸变是系统性的、可预测的，通过标定完全可以补偿。但需注意：畸变标定精度受限于标定板精度、特征点提取精度（通常为0.1像素）以及温度变化引起的镜头参数漂移。在温度变化超过±5℃的环境中，建议增加温度补偿系数。
2.6 常见误区与总结复盘
在理解畸变物理根源时，工程师常陷入以下误区：
- 误区一：将畸变与模糊混为一谈。畸变是几何位移，不降低MTF（调制传递函数）；模糊是分辨率下降。两者可独立存在，也可同时存在。测量时必须区分：畸变用几何标定评估，清晰度用MTF或SFR评估。
- 误区二：认为薄透镜没有畸变。即使理想薄透镜，当视场角较大时，由于球面折射的非线性，仍会产生畸变。薄透镜只是消除了厚度影响，并未消除视场相关的像差。
- 误区三：忽略光阑位置的影响。同一透镜，光阑位置不同，畸变符号可能反转。在镜头选型时，必须确认光阑位置与畸变方向的关系。
- 误区四：认为畸变是旋转对称的。实际镜头由于装配误差、镜片偏心，畸变往往包含非对称分量（如切向畸变）。高精度测量必须同时考虑径向畸变和切向畸变。
总结复盘：本章的核心逻辑链为：费马原理→折射定律→高斯光学→理想成像模型（线性映射）→实际系统偏离线性→畸变产生。畸变的物理根源在于主光线通过光学系统时，不同视场的横向放大率不一致，其数学表征为奇次多项式展开。理解这一链条，是后续进行畸变标定、图像校正和测量精度分析的前提。在工程实践中，务必牢记：畸变是系统误差，可通过标定消除；但标定本身的精度决定了最终测量精度的上限。