4、干涉检测原理（二）：相移干涉法（PSI）与波长调谐干涉法，相位提取算法（如四步、五步相移），以及相位解包裹技术
在干涉检测中，干涉条纹的强度变化本身并不直接给出面形误差的符号与数值。传统静态干涉图只能通过条纹弯曲程度做定性判断，难以达到纳米级定量精度。相移干涉法（Phase Shifting Interferometry, PSI）通过在参考光路中引入已知相位偏移，采集多帧干涉图，利用相位提取算法逐像素计算出包裹相位（wrapped phase），再通过相位解包裹（phase unwrapping）恢复连续的真实相位分布。波长调谐干涉法（Wavelength Tuning Interferometry, WTI）则是PSI的一种工程实现方式，特别适用于Fizeau型干涉腔，通过改变激光波长等效改变干涉腔相位，避免机械移相带来的误差。本章围绕PSI与WTI的原理、相位提取算法、解包裹技术展开，并结合实际检测场景给出完整案例。
4.1 相移干涉法的核心原理与数学模型
相移干涉法的基本思想是在参考光与测试光之间引入已知的、随时间变化的相位差，使得干涉场强度呈现余弦调制。设第i帧干涉图的强度分布为：
I_i(x,y) = I_0(x,y) + I_m(x,y)·cos[φ(x,y) + δ_i]其中，I₀为背景光强，I_m为调制幅度，φ(x,y)为待测相位，δ_i为第i帧的相移量。当相移量δ_i已知且等间隔分布时，可以通过多帧强度值联立求解φ(x,y)。最常用的相移方式包括：
- 时间相移（Temporal PSI）：在时间序列上依次采集N帧干涉图，每帧之间引入固定相移。适用于静态或准静态被测对象，是工业检测的主流方案。
- 空间相移（Spatial PSI）：在单帧图像中通过偏振或光栅分光同时获取多路相移干涉图，适用于动态测量，但空间分辨率会有所损失。
- 波长调谐相移（Wavelength Tuning）：通过改变激光波长改变干涉腔的光学相位差，等效于移相。特别适用于Fizeau型干涉仪，因为其参考面与测试面共光路，机械移相困难。
波长调谐的移相量由干涉腔的光程差（OPD）决定。对于Fizeau干涉腔，若参考面与测试面间距为L，则波长变化Δλ引起的相位变化为：
Δφ = 4π·L·Δλ / λ²因此，要实现2π的相移，所需的波长调谐量为Δλ = λ²/(2L)。例如，对于λ=632.8nm、L=10mm的干涉腔，完成2π相移仅需调谐约0.02nm，这对激光器的波长稳定性提出了较高要求。
4.2 相位提取算法：四步、五步相移及其误差特性
相位提取算法的核心是从N帧干涉图中求解φ(x,y)。不同帧数对应不同的抗噪能力和误差抑制特性。工程中最常用的是四步相移和五步相移算法。
4.2.1 四步相移算法
四步相移法采集四帧干涉图，相移量分别为0、π/2、π、3π/2。对应的强度表达式为：
I_1 = I_0 + I_m·cos(φ)
I_2 = I_0 + I_m·cos(φ + π/2) = I_0 - I_m·sin(φ)
I_3 = I_0 + I_m·cos(φ + π) = I_0 - I_m·cos(φ)
I_4 = I_0 + I_m·cos(φ + 3π/2) = I_0 + I_m·sin(φ)联立求解可得包裹相位：
φ_w = arctan[(I_4 - I_2) / (I_1 - I_3)]四步法的优点是计算简单、速度快，对探测器背景光强不敏感（I₀被差分消除）。但其对相移误差较为敏感——若实际相移量偏离π/2，会产生周期性误差，表现为面形上的“波纹”状伪影。
4.2.2 五步相移算法（Hariharan算法）
五步相移法采集五帧干涉图，相移量为-π、-π/2、0、π/2、π。其相位计算公式为：
φ_w = arctan[ 2·(I_2 - I_4) / (2·I_3 - I_1 - I_5) ]五步法对相移误差具有一阶不敏感性。当实际相移量存在偏差ε时，四步法的相位误差与ε成正比，而五步法的相位误差与ε²成正比。因此，在相移器标定不理想或环境扰动较大的场景下，五步法显著优于四步法。
| 算法 | 帧数 | 相移量 | 相位公式 | 对相移误差敏感度 | 适用场景 |
|---|---|---|---|---|---|
| 三步法 | 3 | -π/2, 0, π/2 | arctan[(I₁-I₃)/(2I₂-I₁-I₃)] | 高 | 快速测量，精度要求低 |
| 四步法 | 4 | 0, π/2, π, 3π/2 | arctan[(I₄-I₂)/(I₁-I₃)] | 中等 | 通用工业检测 |
| 五步法 | 5 | -π, -π/2, 0, π/2, π | arctan[2(I₂-I₄)/(2I₃-I₁-I₅)] | 低（一阶不敏感） | 高精度计量、环境扰动较大 |
| Carre算法 | 4 | 未知等间隔 | 需先估计相移量 | 对相移量不敏感 | 相移器未标定或非线性 |
4.3 相位解包裹技术：从包裹相位到连续相位
相位提取算法输出的φ_w(x,y)被截断在(-π, π]区间内，称为包裹相位。当被测面形起伏超过λ/2时，包裹相位图中会出现2π跳变。相位解包裹的目的就是恢复这些跳变处的整数倍2π，得到连续的真实相位分布。
4.3.1 空间解包裹（Spatial Unwrapping）
空间解包裹是最经典的方法，其基本思路是逐行或逐列比较相邻像素的相位差，当差值超过π时，自动加上或减去2π。具体流程为：
- 计算相邻像素的相位差：Δφ = φ(i+1,j) - φ(i,j)
- 若Δφ > π，则后续像素相位减去2π；若Δφ < -π，则后续像素相位加上2π
- 遍历整个图像，得到连续相位分布
空间解包裹的优点是实现简单、计算速度快。但其致命弱点是：当存在噪声、阴影、断点或高斜率区域时，误差会沿路径传播，形成“拉线”或“斑块”状错误。对于包含不连续面形（如台阶、沟槽）的样品，空间解包裹几乎无法正确恢复。
4.3.2 时间解包裹（Temporal Unwrapping）
时间解包裹是波长调谐干涉法特有的优势。其原理是：在波长调谐过程中，每个像素的相位随时间线性变化。通过采集多组不同波长下的干涉图，可以逐像素独立求解相位随时间的变化率，从而直接得到连续相位，无需空间邻域比较。
具体实现方式为：在波长调谐范围内采集M组干涉图，每组包含N帧相移图。对每个像素，利用最小二乘法拟合相位随波长的变化斜率，该斜率正比于该像素对应的光程差。由于每个像素独立处理，时间解包裹天然免疫于空间不连续问题。
4.3.3 质量图引导解包裹
对于复杂面形，可采用质量图引导的路径跟踪算法。质量图可以是调制幅度I_m(x,y)、相位梯度或局部一致性指标。算法优先从高质量区域开始解包裹，逐步向低质量区域扩展，避免误差沿低质量路径传播。常用的质量图引导算法包括：
- 洪水填充法（Flood Fill）：从质量最高的种子点开始，按质量降序扩展
- 枝切法（Branch Cut）：识别残差点并放置枝切线，阻止误差传播
- 最小二乘法（Least Squares）：全局优化，使解包裹结果与包裹相位的梯度差异最小
4.4 案例拆解：Fizeau干涉仪波长调谐检测平面镜面形
某光学工厂需要对口径为100mm的平面反射镜进行面形检测，要求PV值测量重复性优于λ/100（λ=632.8nm）。工厂原有静态干涉条纹判读法只能达到λ/20的重复性，无法满足需求。决定采用波长调谐干涉法进行升级。
初始条件与约束：
- 干涉仪类型：Fizeau型，参考面与测试面间距L=15mm
- 激光器：可调谐He-Ne激光器，中心波长632.8nm，调谐范围±0.5nm
- 相机：1024×1024像素，12位动态范围，帧率30fps
- 环境：恒温20±0.1℃，气浮隔振平台
关键参数计算：
完成2π相移所需波长调谐量：
Δλ = λ²/(2L) = (632.8×10⁻⁹)² / (2×15×10⁻³) = 0.0133 nm
若采用四步相移，需采集4帧，总调谐范围 = 3×0.0133 = 0.04 nm
若采用五步相移，需采集5帧，总调谐范围 = 4×0.0133 = 0.053 nm实施步骤：
- 将平面反射镜安装到Fizeau干涉仪的测试工装上，调整倾斜使条纹密度降至最低（零条纹模式）
- 设置波长调谐序列：从λ₀开始，每步增加0.0133nm，共采集5帧干涉图（采用五步相移算法）
- 对每帧图像进行暗场扣除和平场校正，消除探测器响应不均匀性
- 利用五步相移公式计算包裹相位φ_w(x,y)
- 由于平面镜面形连续，采用空间解包裹（质量图引导的洪水填充法）恢复连续相位
- 从连续相位中去除倾斜项和离焦项（通过Zernike多项式拟合），得到纯面形误差
- 计算PV值和RMS值，输出面形图
结果与对比：
| 指标 | 静态条纹判读法 | 波长调谐PSI（四步） | 波长调谐PSI（五步） |
|---|---|---|---|
| PV重复性 | ±λ/20 | ±λ/80 | ±λ/120 |
| RMS重复性 | ±λ/50 | ±λ/200 | ±λ/300 |
| 单次测量时间 | 约30秒（人工判读） | 0.5秒 | 0.7秒 |
| 对振动敏感度 | 低（静态） | 高 | 中 |
复盘与经验：
五步法相比四步法在重复性上提升了约50%，主要原因是波长调谐过程中激光器存在微小的非线性漂移，五步法对这类误差的抑制能力更强。但五步法对振动更敏感，因为采集时间更长。在实际生产中，若隔振条件良好（振动幅值小于λ/100），优先选择五步法；若隔振条件一般，则四步法配合多次平均更为稳妥。
4.5 常见误区与总结复盘
在PSI/WTI的实际应用中，工程师常陷入以下误区：
- 误区一：认为帧数越多精度越高。实际上，帧数增加会延长采集时间，增加振动和气流扰动引入的误差。在环境条件不佳时，三步法可能比五步法更可靠。
- 误区二：忽略探测器非线性。CCD/CMOS的光电响应并非理想线性，当干涉图调制深度较大时，非线性会产生高阶谐波误差，表现为面形上的“三瓣”或“四瓣”伪影。建议在算法中加入非线性校正项。
- 误区三：解包裹前不检查残差点。残差点（residue）是包裹相位图中梯度不守恒的位置，通常由噪声或欠采样引起。若不处理残差点直接解包裹，误差会沿路径传播。建议先计算残差点分布，再选择适当的枝切线或质量图引导策略。
- 误区四：忽略参考面误差。Fizeau干涉仪的参考面并非理想平面，其面形误差会直接叠加到测量结果中。高精度检测时，需先对参考面进行标定，并在数据处理中扣除。
总结：相移干涉法通过多帧强度采样实现了相位的高精度定量提取，其核心在于相移的准确性和相位提取算法的鲁棒性。四步法适合快速检测，五步法适合高精度计量，波长调谐则提供了无机械移相的替代方案。相位解包裹是PSI流程中最容易出错的环节，需要根据被测面形的连续性和噪声水平选择合适的策略。在实际工程中，建议结合质量图引导解包裹和多次测量平均，以获得稳定可靠的面形检测结果。