4. 图像畸变基础
在HUD（Head-Up Display，平视显示）系统中，图像从像源（PGU）发出，经折转镜、自由曲面镜、风挡玻璃等多级光学面传递，最终投射到驾驶员前方。理想情况下，像源上的直线在虚像面上仍应是直线，像源上的矩形应保持矩形比例。但真实光学系统与风挡玻璃的曲面特性，会使图像产生位置、形状、比例上的偏差，这类偏差统称为图像畸变。理解畸变的分类与成因，是后续建立校正模型、设计标定流程、评估校正效果的基础。
一、几何畸变的定义与数学描述
几何畸变（Geometric Distortion）指图像中像素点的实际成像位置与理想针孔成像位置之间的偏差。它不改变像素灰度值，只改变像素的空间坐标。设理想成像点为 (x_u, y_u)，实际成像点为 (x_d, y_d)，则畸变可表示为：
x_d = x_u + δx(x_u, y_u)
y_d = y_u + δy(x_u, y_u)
其中 δx、δy 为畸变量，是理想坐标的非线性函数。在HUD工程中，畸变量通常以像素、毫米或角分（arcmin）为单位度量。若畸变随位置线性变化，表现为整体缩放或平移，称为线性畸变；若随位置非线性变化，则表现为直线弯曲、局部拉伸或压缩，称为非线性畸变。HUD校正的核心任务，就是估计 δx、δy 并构造反向映射，使实际成像逼近理想位置。
二、径向畸变：沿半径方向的拉伸与压缩
径向畸变（Radial Distortion）是最常见的畸变类型，由透镜或曲面镜的球面特性引起，畸变量沿图像中心到边缘的半径方向变化。其经典模型为：
x_d = x_u (1 + k1 r² + k2 r⁴ + k3 r⁶)
y_d = y_u (1 + k1 r² + k2 r⁴ + k3 r⁶)
r² = x_u² + y_u²
当 k1 > 0 时，边缘点向外扩张，形成枕形畸变（Pincushion）；当 k1 < 0 时，边缘点向内收缩，形成桶形畸变（Barrel）。HUD中自由曲面镜和风挡玻璃的曲率组合，常使虚像边缘出现明显的枕形或桶形趋势。径向畸变的特点是：中心区域畸变小，越靠近边缘畸变越大；且畸变方向沿半径，不改变图像的旋转对称性。
三、切向畸变：装配偏差引起的非对称偏移
切向畸变（Tangential Distortion）由光学元件装配时透镜平面与成像平面不平行引起，表现为图像在切向方向上的偏移。其模型为：
x_d = x_u + [2 p1 x_u y_u + p2 (r² + 2 x_u²)]
y_d = y_u + [p1 (r² + 2 y_u²) + 2 p2 x_u y_u]
其中 p1、p2 为切向畸变系数。切向畸变通常比径向畸变小一个数量级，但在HUD这种多镜面、长光路系统中，折转镜与自由曲面镜的微小角度偏差会累积成可观测的切向偏移。其典型表现是图像出现“倾斜”或“单侧拉伸”，且畸变分布不再关于中心对称。在标定中，切向畸变常与径向畸变联合求解，但若装配偏差过大，仅靠多项式模型难以完全补偿。
四、梯形畸变：投影几何引起的透视变形
梯形畸变（Keystone Distortion）属于透视畸变，源于像源平面与虚像平面不平行，或观察方向与投影光轴不重合。其几何本质是投影变换：
[x_d, y_d, w]ᵀ = H · [x_u, y_u, 1]ᵀ
x_d' = x_d / w, y_d' = y_d / w
其中 H 为3×3单应矩阵。梯形畸变使矩形图像变为梯形，表现为上宽下窄或左宽右窄。在HUD中，风挡玻璃的楔形角、安装倾角以及驾驶员眼盒位置偏移，都会引入梯形畸变。与径向畸变不同，梯形畸变是全局性的线性透视变形，可通过单应变换精确校正；但若与径向畸变耦合，则需先分离再分别补偿。
五、非线性畸变分类与工程对照
实际HUD系统的畸变往往是多种类型的叠加。为便于工程处理，可按数学特性将畸变分为以下几类，并对应不同的校正策略：
| 畸变类型 | 数学特性 | 典型成因 | HUD中的表现 | 校正方法 |
|---|---|---|---|---|
| 径向畸变 | 沿半径非线性 | 球面镜、自由曲面曲率 | 边缘枕形/桶形弯曲 | 多项式反向映射 |
| 切向畸变 | 非对称切向偏移 | 镜片装配不平行 | 图像倾斜、单侧拉伸 | 切向系数联合标定 |
| 梯形畸变 | 全局透视线性 | 风挡楔角、安装倾角 | 矩形变梯形 | 单应矩阵校正 |
| 非线性局部畸变 | 局部高阶非线性 | 风挡局部曲率突变 | 局部波纹、扭曲 | 网格标定+插值 |
| 色差畸变 | 波长相关偏移 | 折射率随波长变化 | RGB边缘错位 | 通道独立校正 |
工程提示：在HUD畸变建模中，建议先通过单应变换消除梯形畸变，再拟合径向与切向系数，最后用残差网格处理局部非线性畸变。这种“先线性后非线性”的分层策略，可显著降低高阶多项式拟合的阶数，提高数值稳定性。
六、案例场景：某型HUD的畸变实测与分解
某型车载HUD在标定中发现：虚像四角向外拉伸约12像素，中心区域基本无偏移；同时图像整体向右下倾斜约0.8°。工程师按以下流程分解畸变：
- 用棋盘格标定板采集虚像，提取角点实际坐标与理想坐标。
- 计算单应矩阵
H，发现梯形分量占总畸变的35%，主要来自风挡安装倾角。 - 去除梯形分量后，剩余畸变呈明显径向对称，拟合得
k1 = -0.021，判定为桶形畸变。 - 切向系数
p1 = 0.0007、p2 = -0.0004，量级较小，与折转镜装配偏差吻合。 - 残差网格显示右上角存在局部波纹，对应风挡玻璃局部曲率变化，采用双三次插值补偿。
最终校正后，全视场畸变残差从初始的18像素降至2.3像素，满足HUD虚像畸变小于3像素的验收要求。
七、常见误区与复盘
误区一：把梯形畸变当作径向畸变处理。梯形畸变是全局透视变形，若强行用径向多项式拟合，会导致中心区域过校正、边缘欠校正。应先分离单应分量。
误区二：忽略切向畸变的累积效应。单镜片切向畸变虽小，但HUD光路中多个镜面串联后，切向偏差会线性累积，不可简单忽略。
误区三：用单一多项式拟合全视场。风挡玻璃边缘曲率变化剧烈，高阶多项式易产生振荡。建议采用“全局低阶+局部网格插值”的混合模型。
误区四：只关注静态畸变。眼盒位置变化会改变观察视角，引入动态梯形畸变。校正模型需在眼盒范围内进行多位置标定与验证。
掌握几何畸变的定义与分类，是建立HUD畸变校正模型的第一步。径向、切向、梯形与非线性畸变各有明确的数学形式和物理成因，工程中需先分解、再建模、后验证，才能实现高精度、高鲁棒性的风挡校正。