[論文レビュー] 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

ひとこと要約

本稿では、データ適合項のヘッセ行列をポイント spread 関数（PSF）との畳み込みとして近似することで、広視野の電波干渉計画像化を高速化する、pfb-clean と呼ばれるプリコンディショニング付き前向き後向きプロキシマルアルゴリズムを提案する。この手法により勾配評価回数をわずか 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.

研究の動機と目的

• 広視野電波干渉計測定における高度なデコンボリューションアルゴリズムの高い計算コストが、優れた性能にもかかわらずその導入を制限している問題に対処すること。
• 大規模な干渉計測定データの画像化に適した、高速収束を実現する実用的なプリコンディショナーを開発すること。
• 現代の画像化アルゴリズムが要請する高価な全測定演算子の適用回数を削減すること。
• 広視野連続スペクトル画像化において、計算複雑度を著しく低減しつつも、画像品質と耐性を維持すること。

提案手法

• 本手法は、データ適合項のヘッセ行列を装置のポイントスプレッド関数（PSF）との畳み込みとして近似することで、効率的なプリコンディショニングを可能にする。
• σ が逆行列の正則化と効率的な共役勾配解法の実現を可能にする、U = I^psf ⋆ + σId というプリコンディショナーを用いる。
• 大きなステップサイズ（γ ≲ 1）を用いるプリコンディショニング付き前向き後向きプロキシマルスキームを採用し、正確な勾配評価回数を最小限に抑える。
• スケーラブルかつ分散処理可能な測定演算子評価のため、dask に基づく wgridder および dask-ms の実装を活用する。
• 非滑らか正則化子（例：l2,1 範囲）を扱うために、Condat（2013）の原双対アルゴリズムをプロキシマルステップに用いる。
• 物理的に意味のある画像再構成を保証するため、正の制約を正則化子内のインジケータ関数によって強制する。

実験結果

リサーチクエスチョン

• RQ1 PSF を用いたヘッセ行列近似は、広視野電波画像化におけるプロキシマル勾配アルゴリズムのプリコンディショニングに有効に機能するか？
• RQ2 このプリコンディショナーは、画像化アルゴリズムにおける全測定演算子の適用回数をどの程度削減できるか？
• RQ3 画像品質および残差ダイナミクスの観点から、このプリコンディショニング付きアルゴリズムはマルチスケール CLEAN と比べてどの程度の性能を示すか？
• RQ4 高速収束を達成しつつ、キャリブレーションアーティファクトに対して耐性を維持できるか？

主な発見

• pfb-clean は、通常の高度なアルゴリズムと比較してはるかに少ない 8 回の勾配評価で、wsclean のマルチスケール CLEAN と同等の画像品質を達成した。
• 本手法は自然重み付けされた残差を生成し、ms-clean が負の残差を生じるのとは異なり、フラックスの正の性質を尊重した。
• 1 回の反復あたりのコストが高価であるものの、画像サイズの FFT およびウェーブレット変換を活用することで、並列処理とスケーラビリティが実現された。
• より表現力の高いウェーブレットベースの正則化子（Ψ）を組み合わせた場合、キャリブレーションアーティファクトに対する耐性が向上したが、正則化が強いほどアーティファクトへの感受性が増加した。
• プリコンディショナーの使用により、大きなステップサイズ（γ ≲ 1）が可能となり、収束に必要な正確な勾配評価回数が著しく削減された。