[论文解读] A practical preconditioner for wide-field continuum imaging of radio interferometric data

Hertzog L. Bester, Audrey Repetti+3 more|arXiv (Cornell University)|Jan 20, 2021

Soil Moisture and Remote Sensing被引用 4

一句话总结

本文提出 pfb-clean，一种预条件前向-后向邻近算法，通过将数据保真项的 Hessian 矩阵近似为与点扩散函数（PSF）的卷积，从而加速宽场射电干涉成像。该方法将梯度计算次数减少至仅 8 次，实现了快速、高动态范围的成像，且相比传统 CLEAN 方法在抗校准伪影方面具有更高的鲁棒性。

ABSTRACT

The celebrated CLEAN algorithm has been the cornerstone of deconvolution algorithms in radio interferometry almost since its conception in the 1970s. For all its faults, CLEAN is remarkably fast, robust to calibration artefacts and in its ability to model point sources. We demonstrate how the same assumptions that afford CLEAN its speed can be used to accelerate more sophisticated deconvolution algorithms.

研究动机与目标

• 为解决先进去卷积算法在宽场射电干涉测量中计算成本过高的问题，尽管其性能优越，但限制了实际应用。
• 开发一种实用的预条件器，以实现大规模干涉测量数据成像中邻近梯度方法的快速收敛。
• 减少现代成像算法中对昂贵的完整测量算子应用次数的需求。
• 在显著降低计算复杂度的同时，保持图像质量与鲁棒性，适用于宽场连续谱成像。

提出的方法

• 该方法将数据保真项的 Hessian 矩阵近似为与仪器点扩散函数（PSF）的卷积，从而实现高效的预条件处理。
• 采用预条件器 U = I^psf ⋆ + σId，其中 σ 用于正则化逆运算并支持高效的共轭梯度求解。
• 算法采用预条件前向-后向邻近方案，使用较大的步长（γ ≲ 1），从而最小化精确梯度计算的次数。
• 利用基于 dask 的 wgridder 和 dask-ms 实现，支持可扩展的、分布式测量算子计算。
• 邻近步骤采用 Condat（2013）提出的原始-对偶算法，以处理如 l2,1 范数等非光滑正则化项。
• 通过在正则化项中引入指示函数，强制执行非负性约束，确保图像重建结果具有物理意义。

实验结果

研究问题

• RQ1 基于 PSF 的 Hessian 近似能否有效预条件宽场射电成像中的邻近梯度算法？
• RQ2 该预条件器在多大程度上可减少成像算法中完整测量算子的应用次数？
• RQ3 与多尺度 CLEAN 相比，该预条件算法在图像质量与残差动态范围方面表现如何？
• RQ4 该方法在实现更快收敛的同时，能否保持对校准伪影的鲁棒性？

主要发现

• pfb-clean 仅通过 8 次梯度计算，便实现了与 wsclean 的多尺度 CLEAN 相当的图像质量，远少于先进算法的典型次数。
• 该方法生成了自然加权的残差，尊重通量的正性，而 ms-clean 可能产生负的残差。
• 尽管每次迭代的计算成本较高，但 pfb-clean 利用图像尺寸的 FFT 和小波变换，实现了高效的并行化与可扩展性。
• 当结合更具表达力的小波基正则化项（Ψ）时，该算法在抗校准伪影方面表现出更优的鲁棒性，但更强的正则化会增加对伪影的敏感性。
• 通过使用预条件器，算法实现了较大的步长（γ ≲ 1），从而大幅减少了达到收敛所需的完整梯度计算次数。