3. 哈特曼波前传感器：微透镜阵列原理、质心计算、波前重构算法

哈特曼波前传感器，圈内人常叫它 SHWFS。这东西是自适应光学系统里最常用的波前探测工具，没有之一。我入行那会儿，第一次在实验室看到 CCD 上那密密麻麻的光斑阵列，说实话有点懵——这玩意儿怎么就能反推出波前畸变呢？后来亲手搭了一套系统，才真正搞明白它的脾气。

3.1 微透镜阵列：把波前切成小块

原理其实很朴素。一个理想平面波，经过微透镜阵列后，每个子透镜都会在焦平面上形成一个完美的艾里斑。但如果波前有畸变，局部波前的倾斜就会让光斑偏离理想位置。说白了，微透镜阵列就像一把尺子，把连续的波前离散成一个个小区域，每个区域的光斑偏移量，就反映了该区域的局部波前斜率。

微透镜阵列的制造工艺很关键。我见过一些廉价方案，直接用光刻胶热熔做透镜，但一致性差，质心标定能飘出半个像素。工业级产品通常用熔融石英刻蚀，或者用液晶空间光调制器模拟——后者灵活，但衍射效率是个坑。

参数	典型值	我的经验
子孔径数	10×10 ~ 40×40	20×20 是性价比平衡点
透镜焦距	5 ~ 50 mm	长焦提高灵敏度，但缩小动态范围
填充因子	> 95%	低于 90% 会漏光，重构误差大
阵列间距	100 ~ 500 μm	匹配 CCD 像素尺寸，避免欠采样

3.2 质心计算：从光斑到斜率

拿到光斑图像后，第一步就是找每个子孔径的光斑质心。最常用的方法是重心法，也叫质心法。公式很简单：

xc = Σ(x_i * I_i) / Σ(I_i)
yc = Σ(y_i * I_i) / Σ(I_i)

其中 I_i 是像素灰度值，x_i、y_i 是像素坐标。但实际用起来，坑不少。

还有一种改进算法叫加窗质心法。就是在计算前，先对光斑区域加一个高斯窗或矩形窗，把边缘噪声滤掉。我在一个激光通信项目里试过，加窗后质心抖动从 0.05 像素降到了 0.02 像素，效果很明显。

质心计算精度直接决定了波前重构的成败。你想想看，如果质心误差 0.1 像素，换算成波前斜率误差可能就有几个微弧度，再积分成波前，PV 值能差出 λ/10。所以，我每次做实验前，都会先做一次质心标定——用理想平面波打进去，记录每个子孔径的参考质心位置。

3.3 波前重构算法：从斜率到波前

有了每个子孔径的 x、y 方向斜率，接下来就是波前重构。主流方法分两类：区域法和模式法。我两种都用过，各有各的脾气。

3.3.1 区域法：逐点积分

区域法的思路很直接：把波前离散成网格点，每个网格点的相位值通过相邻子孔径的斜率积分得到。最经典的是 Fried 模型和 Hudgin 模型。

• Fried 模型：斜率定义在网格点之间，适合正方形排列的子孔径。
• Hudgin 模型：斜率定义在网格点上，适合六边形排列。

区域法的核心是解一个大型稀疏线性方程组。我一般用最小二乘法，配合迭代求解器（比如共轭梯度法）。优点是稳定，对噪声不敏感；缺点是计算量大，子孔径一多，矩阵维度就上去了。

# 伪代码：区域法波前重构
输入：斜率矩阵 Sx, Sy (M×N)
输出：波前相位矩阵 W (M×N)

1. 构建系数矩阵 A（稀疏，尺寸 2MN × MN）
2. 构建观测向量 b = [Sx(:); Sy(:)]
3. 求解最小二乘：W = (A^T A)^{-1} A^T b
4. 去除波前整体倾斜和活塞

我曾经在一个 40×40 子孔径的系统里用区域法，矩阵维度接近 3200×1600，直接求逆根本跑不动。后来改用预条件共轭梯度法，迭代 50 步就收敛了，速度提升了两个数量级。

3.3.2 模式法：用 Zernike 多项式拟合

模式法就不一样了。它假设波前可以用一组正交基函数展开，最常用的是 Zernike 多项式。每个 Zernike 模式对应一种像差：离焦、像散、彗差、球差……

做法是：先根据 Zernike 多项式的导数，建立斜率到模式系数的线性关系。然后解一个超定方程组，得到各阶 Zernike 系数。最后用这些系数合成波前。

Zernike 阶数	像差名称	典型应用
Z4	离焦	调焦补偿
Z5/Z6	像散	光学系统装调
Z7/Z8	彗差	离轴系统诊断
Z11	球差	大口径系统

模式法还有一个好处：可以直接得到像差系数，方便做闭环控制。我在自适应光学系统中，经常把 Zernike 系数直接送给变形镜控制器，省去了波前重构的中间步骤。

3.4 两种方法的对比

区域法和模式法没有绝对的好坏，看场景。我整理了一个对比表，方便你选型：

对比项	区域法	模式法
计算复杂度	高（解大型稀疏矩阵）	低（解小规模超定方程）
对噪声的鲁棒性	好	中等（高阶模式易受噪声干扰）
空间分辨率	高（每个子孔径一个点）	低（受限于 Zernike 阶数）
输出形式	波前相位图	Zernike 系数
适用场景	高频像差、不规则波前	低频像差、闭环控制

我个人习惯是：做诊断分析时用区域法，能看到波前的细节；做闭环控制时用模式法，直接驱动变形镜。如果系统资源允许，也可以两种都跑，互相验证。

3.5 哈特曼波前传感器的工作流程

下面这张 SVG 图，把整个流程串起来了。从光斑图像到波前输出，每一步都离不开前面讲的内容。

3.6 实战中的避坑指南

最后，分享几个我踩过的坑：

• 质心标定不能省：我曾经偷懒，直接用理论位置做参考，结果重构出来的波前全是假的倾斜。后来老老实实打平面波标定，误差立刻降下来。
• 注意子孔径的串扰：微透镜阵列的填充因子不是 100%，相邻子孔径的光斑可能会重叠。我一般会调整透镜焦距，让光斑尺寸控制在子孔径间距的 1/3 以内。
• 模式法别用太多阶：Zernike 阶数越高，对噪声越敏感。我试过取到 36 阶，结果高阶系数全是噪声，波前反而更差了。15~21 阶是个安全区间。
• 区域法别忘了边界条件：波前边缘的子孔径只有一侧有斜率数据，重构时容易产生边界畸变。我习惯在边界加虚拟子孔径，或者用平滑约束。