[論文レビュー] Fourier Phase Retrieval with a Single Mask by Douglas-Rachford Algorithm

Pengwen Chen, Albert Fannjiang|arXiv (Cornell University)|Sep 2, 2015

Advanced X-ray Imaging Techniques参考文献 39被引用数 9

ひとこと要約

本稿では、1つのランダムな位相マスクを用いたフーリエ位相再構成にDouglas-Rachford（DR）アルゴリズムの適用を提案し、一意な固定点への局所的幾何的収束を証明するとともに、数値的にグローバル収束を示している。この手法は、特にX線自由電子レーザー画像法における「破壊前に回折」の応用において、物体構造の再構成における一意性と安定性を保証する。

ABSTRACT

Douglas-Rachford (DR) algorithm is analyzed for Fourier phase retrieval with a single random phase mask. Local, geometric convergence to a unique fixed point is proved with numerical demonstration of global convergence.

研究の動機と目的

• 1つのランダム位相マスクを用いたフーリエ位相再構成におけるDouglas-Rachfordアルゴリズムの収束保証を確立すること。
• 放射線損傷のため複数回の測定が制限される「破壊前に回折」X線画像法における位相再構成問題に取り組むこと。
• 現実的なスキャッタリングポテンシャルのセクター制約下で、1つの符号化回折パターンを用いて解の一意性を回復すること。
• HIO や ER のようなヒューリスティックな反復アルゴリズムの代わりに、理論的裏付けが強く数値的に安定した代替手法を提供すること。
• 最小限の測定オーバーヘッドで、結晶性でない単粒子イメージングへの位相再構成の適用範囲を拡張すること。

提案手法

• Douーリングス・ラッハフォードアルゴリズムを位相再構成に適用するため、2つの射影作用素を定義する：1つはフーリエ変換を介した回折場の集合への射影、もう1つは測定された強度と一致する制約への射影。
• アルゴリズムは次の式に従って反復的に解を更新する：$ y^{(k+1)} = y^{(k)} + P_1(2P_2 - I)y^{(k)} - P_2y^{(k)} $。計算効率を高めるために高速フーリエ変換（FFTs）を活用する。
• 理論的解析により、1つのランダム位相マスクとスキャッタリングポテンシャルのセクター制約の下で、局所的かつ幾何的収束が一意な固定点へと成立することを証明した。
• 解の一意性をセクター制約なしに保証するために、2パターンの測定（符号化済みおよび通常の回折）に拡張した。
• 数値実験により、グローバル収束の挙動が確認され、理論的予測を支持した。
• 代数的数論に基づくアプローチを採用し、z変換と既約性の議論を用いて、ランダムマスキング下での一意性を確立した。

実験結果

リサーチクエスチョン

• RQ1 1つのランダム位相マスクを用いたフーリエ位相再構成において、Douglas-Rachfordアルゴリズムが局所的幾何的収束を達成できるか？
• RQ2 一般にグローバル収束保証がないにもかかわらず、実際の応用ではアルゴリズムが真の解へグローバルに収束するか？
• RQ3 どのような条件下で、1つの符号化回折パターンが位相再構成解の一意性を保証するか？
• RQ4 スキャッタリングポテンシャルのセクター制約が、DRアルゴリズムの解の一意性および収束性に与える影響は何か？
• RQ5 2パターンの測定（符号化済み＋通常）にDRアルゴリズムを拡張することで、セクター制約の必要性を排除できるか？

主な発見

• 1つのランダム位相マスクを用いた場合、スキャッタリングポテンシャルのセクター制約の下で、Douglas-Rachfordアルゴリズムは一意な固定点への局所的幾何的収束を示す。
• 数値的結果により、初期値が解に近いとは限らない状況でも、DRアルゴリズムが真の解へグローバルに収束することが確認された。
• スキャッタリングポテンシャルが複素平面の既知のセクター内にある限り、1つのランダムマスクを用いることで、解の一意性が高確率で保証される。
• 理論的解析により、マスキングされたオブジェクトのz変換がほとんど確実に既約であり、共役対称でないことが示され、偽の解の発生を防ぐ。
• 2つの回折パターン（符号化済みおよび通常）を用いることで、セクター制約なしに解の一意性がほぼ確実に回復され、より高いロバストネスが得られる。
• FFTと明示的な射影を用いることで、計算が効率的であり、XFELを用いた単粒子イメージングなどの大規模な画像処理応用に適している。