【Optica】抗湍流宏观傅里叶叠层成像：复杂大气环境下的长距离合成孔径超分辨率突破

摘要：针对长距离光学成像中大气湍流引发的极端相位畸变与分辨率退化物理瓶颈，中国科学院光电技术研究所联合南京理工大学团队提出了一种抗湍流宏观傅里叶叠层成像（TRMFP）计算框架。该研究摒弃了传统自适应光学硬件，将散斑干涉统计学先验与闭环叠层相位恢复算法深度融合。通过短曝光图像序列的系综平均（Ensemble Averaging）提取衍射极限信息，并引入全变分（TV）与奇异值分解（SVD）低秩约束联合优化空变残余像差。TRMFP在无需额外波前传感器的情况下，成功在极端动态对流湍流（高达7.5）与光学粗糙表面目标环境中，实现了1.6至2倍的超分辨增益，为复杂真实大气条件下的远距离合成孔径成像确立了严谨的纯计算物理范式。

1. 研究背景与挑战

实际光学系统的空间分辨率受制于两大基础物理屏障：物理孔径尺寸与大气湍流。由瑞利判据可知，系统的极限角分辨率为 ；然而，当成像孔径 超过大气相干长度（Fried参数 ）时，系统分辨率从“孔径受限”跌落至“湍流主导”域。 传统干涉合成孔径（SA）在射电天文中应用成熟，但向光学波段移植时面临致命挑战：光电探测器仅能响应光强而丢失相位，迫使光学SA必须依赖极端精密的机械姿态控制以维持亚波长级的共焦（Confocality）与共相（Cophase）对齐。 宏观傅里叶叠层成像（Macroscopic FP）曾被视为一种优雅的计算破局方案，通过空间相机扫描与交替投影算法实现数值共相。但现有FP技术体系隐含了一个脆弱的理想化假设——静态且确定性的波前。一旦推向长距离真实场景，大气湍流引发的折射率起伏（服从 Kolmogorov 幂律）会带来灾难性破坏：

1. 相干重叠失效：动态时变的空变相位误差直接撕裂了相邻扫描孔径间频域重叠区所需的高度相位一致性。
2. 正向物理模型崩塌：每单次曝光实质上是畸变波前的一个随机实现，彻底违背了FP经典的确定性正向传播模型。传统的静态像差校正算法（如EPRY）或相位多样性（Phase Diversity）技术，在面对毫秒级演化的随机相屏时均无法收敛。

图1：TRMFP系统架构与动态重构流程（对应原论文 Fig. 1）。揭示了系统如何在各个扫描位置采集短曝光图像序列，并通过两阶段的“预处理去畸变+闭环联合像差优化”突破湍流物理限制。

2. 系统、方法与物理原理深度拆解

为在强烈且动态的波前畸变中提取高频信息，该研究严格设定曝光时间小于大气相干时间 （通常小于10 ms），以物理“冻结”瞬间湍流状态。系统的等晕区正向物理退化模型可表述为：

其中， 为合成孔径扫描位置标号， 为当前位置下的第 次独立短曝光， 代表中心相对频移， 为目标复数频域分布， 囊括了系统孔径与该瞬间的随机湍流波前畸变。为破解这一极端病态的逆问题，团队设计了物理逻辑严密的“两阶段降维”策略。

阶段一：基于散斑干涉学的湍流统计抑制（Preliminary Turbulence Mitigation）该步骤的底层物理直觉源自散斑干涉学：尽管畸变严重，短曝光图像仍携带有直至衍射极限的高空间频率能量。通过构建空间卷积退化模型 ，团队巧妙规避了相干FP与非相干散斑算法的数学冲突。

1. 相位恢复：将随机生成的湍流PSF拆解为理想PSF与随机偏移权重的卷积：。通过互相关（Cross-correlation）操作对原始畸变序列进行对齐补偿（Shift-correction），获得 。此时的PSF系综平均 将严格退化为 叠加一个实数常数背景。由于该结果是纯实数且空间对称的，其傅里叶变换具有零相位特性。因此，目标的傅里叶相位可直接由相干平均图像的频域导出：
2. 振幅恢复：高频振幅则通过统计各独立短曝光图像能谱的均方根提取：。
3. 参数抉择：团队在物理边界上进行了权衡，强湍流模拟（至10）采用 帧以保证收敛，而实验考量采集效率则锁定 帧。两者重新组合逆傅里叶变换后，输出了初步剥离严重湍流相位的衍射极限子图。

图2：湍流预处理理论验证与去畸变对比（对应原论文 Fig. 2）。相较于传统单帧FP的伪影崩溃与直接长曝光平均的不可逆高频抹除，统计学重构成功找回了高频频谱边缘信息。

3. 核心创新点深挖

阶段二：空变残余像差校正与联合低秩投影（Residual Aberration Correction, RAC）第一阶段虽然大幅压制了宏观畸变，但针对强湍流，子图像仍残留未完全熨平的空变相位误差（Spatially Varying Residual Aberrations）。若直接拼接，微弱的相位失配仍会引发FP频域更新的连锁崩溃。 团队创新性地构建了一个内嵌正交模式的交替投影联合优化网络。在每个子迭代循环中：

1. 全变分（TV）平滑极值：物理上，光瞳残差应具备空间连续性。采用快速梯度投影（FGP）算法施加TV近端算子：，随后施加二进制支持域掩膜（CTF）强行切断孔径外散逸能量。
2. EPRY空变更新：基于校正后的光瞳 ，沿数据保真度梯度反向传播，同步更新目标频带 与局部光瞳估计。
3. 最核心的物理先验——SVD低秩约束矩阵投影：这是该算法抑制病态发散的“定海神针”。物理学表明，消除主体湍流后的残余畸变在相邻重叠率为85%以上的扫描点之间高度相关，由少数低阶共同正交模式（如离焦、像散等）主导。算法将全孔径光瞳集堆叠为矩阵 ，执行截断奇异值分解（SVD）：团队根据奇异值能谱将秩硬性截断于 （保留95%主导能量）。这一降维打击粗暴且有效，从数学底层剥离了高频随机无序噪声，迫使所有空变光瞳共享同一个低维物理子空间。

图3：联合像差校正算法中光瞳模式的空频演化（对应原论文 Fig. 5）。实验揭示，在D/r0=3时RAC可单独完美收敛；但在D/r0=5时RAC失效，这反证了“预先抑制大尺度湍流（阶段一）”这一前置条件的不可或缺。

4. 性能验证与多维实验结果

文献提供了极度详尽的多尺度验证，囊括数值模拟与两种物理实验形态，确保了数据的严密性。

仿真验证：对光学粗糙表面的鲁棒性真实物体表面通常具有粗糙度，会导致严重的散斑散射。模拟中，团队向目标注入 到 的随机极值相位，构建光学粗糙体。在此恶劣设定及极强湍流（）下，辅以散斑去噪模块的TRMFP展现了压倒性优势，其皮尔逊相关系数（PCC）仍维持在 0.8334，而传统FP已断崖式跌落至 0.5701，彻底失去分辨能力。

1.5米静态相位屏严谨对照实验：采用 639 nm 单模激光照明， 孔径（直径 6.25 mm），机械扫描 网格形成 21.35 mm 的合成孔径，相邻重叠率达 85%。利用旋转相位屏捕获 2420 张短曝光图像。 定量的对比度（）计算表明（以瑞利判据的理论对比度阈值 0.1 为基准）：传统单帧FP极限分辨率骤降至 Group 1 Element 1；而TRMFP逆势将分辨率拉升至 Group 2 Element 2，将原本受湍流限制的角分辨率从 rad 精确恢复至 rad，逼近无湍流理论衍射极限（ rad）的 80%。在具有复杂灰度梯度的纸币样本测试中，传统FP的均方根对比度为 0.2247，完全模糊；TRMFP提升至 0.2857，微缩汉字拼音笔画毫厘毕现。

5米极端实验室动态湍流实验：这是系统工程上的极限施压。靶距增至5米，团队利用两台大功率空调机组制造跨度达12°C（30°C vs 18°C）的室温对流，同时通过两台鼓风机向光路横向注入热风，重现真实大气的强闪烁与剧烈相变。 成像系统改用 （孔径 13.39 mm），曝光时间缩短至极限的 10 ms。TRMFP 依然在分辨率板上实现了由 Group -1 Element 6 到 Group 0 Element 4 的跨越式攀升，绝对分辨力增益达 1.6倍。更重要的是，在对非合作材质（强吸收帆布、高反射硬币压痕、透明塑料标签）的泛化测试中，TRMFP 呈现出极高保真度的重构，结构相似性指数（SSIM）全面压制对照组。

图4：5米对流强湍流环境及多材质泛化验证（对应原论文 Fig. 7）。定量对比度及SSIM数值验证了TRMFP在面对帆布纤维、金属压痕等真实散射体时的极高空间解析稳定性。

5. 总结与展望

1. 跨学科的物理模型贯通：首度在底层算法维度跨越了非相干领域（散斑干涉）与相干合成孔径体系（FP）的鸿沟。
2. 纯计算波前补偿范式：不依赖任何额外的自适应光学（AO）硬件，无需基于深度学习的海量先验数据训练，实现了完全基于物理过程可解释性的全频段湍流抑制。
3. “解耦-降维”联合重构拓扑：利用统计平均剥离强畸变，继而通过SVD低秩约束将高频空变残差束缚于物理可行的低维子空间，彻底解决了病态逆问题的发散难题。
   
   局限性剖析与演进预测（Discussion）：需客观指出，当前TRMFP严格绑定于“等晕区（Isoplanatic）”成像假设。这意味着它在小视场（FOV）内表现卓越，但在面临长基线大视场的严重各向异性（Anisoplanatism）湍流时，单一全局退化核将不再适用。未来的技术演进势必向构建“视场依赖的光瞳函数（Field-dependent pupil functions）”的多等晕拼片模型发展。同时，海量短曝光导致的庞大数据流及其引发的2-3倍算力开销，需通过异构GPU并行运算优化。更长远看，将TRMFP的精细高频波前解算能力与低阶AO系统融合，极有可能在深空观测与光通信领域，触达分辨率超越大气相干长度 物理禁锢的全新边界。

6. 参考文献

Junhao Zhang, Sheng Li, Weilong Wei, Bowen Wang, Kaiyuan Yang, Hongmei   Wang, Qiang Zhou, Haotong Ma, Ge Ren, Chao Zuo, and Zongliang Xie,   "Turbulence-resilient macroscopic Fourier ptychography," Optica 13, 771-785 (2026)