Fast Eikonal Phase Retrieval for High-Throughput Beamlines

技術要約：高スループット・ビームラインのための高速エイコナル位相回復

問題提起
第四世代放射光光源における回折ベースの位相コントラスト顕微鏡法（PPC-µCT）は、試料から検出器までの伝搬距離を長く取ることで、弱屈折特性を持つ構造の高感度なイメージングを可能にする。しかし、この領域では再構成において重大な課題が生じる。標準的な単一距離の位相回復法（例：Paganinフィルタ）は、線形な強度輸送の仮定に依存している。軟組織と強い界面（例：骨と空気、あるいは骨と流体）が混在する不均質な試料では、急峻な位相勾配が強い強度変調を生じさせ、これらの線形仮定を破綻させる。これは、画像品質を低下させ、ダウンストリームの解析を損なう特徴的な非線形アーティファクト（長距離のストリークなど）を引き起こす。

フル・フレネル波動光学の定式化を用いればこれらの効果をモデル化できるが、それには複素波面のピクセル単位での明示的な表現が必要となる。これは厳格なサンプリング制約を課すことになる。位相変化が急速な場合（鋭い界面によって引き起こされる場合）、標準的なフレネル伝搬はエイリアシングや膨大なメモリ要件により、計算上極めて困難になる。一方、以前に導入されたエイコナル位相回復（EPR）アルゴリズムは、エイコナル近似を用いて強い屈折をモデル化することに成功したが、大規模なデータセットに対して数十時間を要するという高い計算コストに悩まされており、現代のビームラインにおける高スループットな取得速度とは相容れないものであった。

手法
著者らは、元のEPRアルゴリズムを、制御された精度を維持しつつ2桁以上高速化する「高速エイコナル位相回復（Fast EPR）」の定式化を提案する。この手法は、伝搬距離 $L$ に関するフレネル伝搬子の二次の漸近展開と、補完的な主要レント・クレンメル・ブリーレン（WKB）波動光学補正に基づいている。

本アプローチは、異なる物理的領域を扱うための2つの補完的なソルバーを用いた統一フレームワークを採用している：

1. ローカル・ソルバー（サブピクセル・シフト）： エイコナルによる横方向のレイ・シフト $\Delta(r)$ が1ピクセル未満（$\eta < 1$）である領域では、二次的な $O(L^2)$ クローズ（closure）を使用する。これには、二次的なエイコナル項（WKB0）と主要な回折圧力項（WKB1）を含む、強度の輸送に関する二次微分展開が含まれる。このローカルモデルは、収束を加速させるために、FFT対角のエネルギー依存逆演算子（Paganin型フィルタ）によって前処理される。
2. ノンローカル・ソルバー（マルチピクセル・シフト）： シフトが1ピクセルを超える（$\eta > 1$）場合、ローカルなテイラー展開は破綻する。このとき、強度は明示的なエイコナル・レイ・マッピングを通じて伝搬されるノンローカルな定式化へと切り替わる。これは、検出器グリッド上での質量保存的な双線形再分配として実装される。検出器の残差は、対応する随伴（転置）演算子（双線形補間として実装）を用いて物体グリッドへと転送された後、FFT対角の前処理器が適用される。

多色性サポート
本フレームワークは、コンパクトなスペクトル離散化を通じて多波長データ（polychromatic data）をサポートする。スペクトルは $N_s$ 個のスペクトル成分の非干渉和としてモデル化される。すべての線ごとの量（位相/吸収パラメータおよびシフト）はエネルギー依存的に扱われ、FFTの効率性を維持しながら、エネルギー依存の輸送と反転を可能にする。

主な貢献

• アルゴリズムの加速： 新しい定式化は、同じGPUリソースに対して正規化した場合、元のEPR実装と比較して、位相回復の計算コストを約580倍（2.8桁）削減する。
• ハイブリッド・ソルバー戦略： ローカル（サブピクセル・シフト用の $O(L^2)$）とノンローカル（マルチピクセル・シフト用のレイ・マッピング）のデュアル・ソルバー戦略の導入により、軟組織のような弱い屈折から強い界面に至るまで、幅広い実験条件における堅牢性を確保している。
• 二次的な物理学： WKB展開の完全な $O(L^2)$ の内容（回折圧力項を含む）を保持することで、一次の線形モデルに残留する非線形ストリーク・アーティファクトを抑制する。
• 多波長効率： 本手法は、高スループットなワークフローに必要な計算効率を犠牲にすることなく、多波長データの正確な再構成を可能にする。

結果
著者らは、羊頭部のHiP-CTデータセットと竹のサンプルによるストレス・テストを用いて手法を検証した：

• アーティファクトの低減： 羊頭部データセットを用いた比較により、二次（$L^2$）モデルが一次（$L^1$）モデルと比較して残留非線形伝搬アーティファクトを有意に低減し、強い勾配付近での軟組織コントラストの描画をより鮮明にすることを示した。
• 収束性： アルゴリズムは迅速に収束する。二次的な物理学を有効にした後、改善の大部分は最初の反復（イテレーション）で達成される。
• 堅牢性： レイ・シフトが検出器のピクセルサイズを超えた竹のサンプルにおいて、ローカル・ソルバーは失敗（第1反復後に発散）したが、ノンローカル・ソルバーは安定かつ正確に動作した。
• パフォーマンス： 8,000枚の放射線透過画像（$256 \times 3104$ ピクセル）のデータセットに対し、Fast EPR法は、単一のGPU上で単色データの場合に約1.1〜1.9分、多波長データ（5つのスペクトル線）の場合に1.47〜1.9分で位相回復を完了した。これは、同等のボリュームに対して約870分を要した元のEPRのベンチマークとは対照的である。

意義
本論文は、この統一されたアプローチが、高スループットなPPC-µCT実験に関連する伝搬条件下において、正確かつ計算効率の高い位相回復を可能にすると主張している。再構成時間を、高スループットなワークフロー（数時間に一度のサンプル取得）に適合するレベルまで短縮することで、本手法は非線形位相回復を日常的なワークフローに組み込むことを容易にする。さらに、明示的な波動光学補正を備えた高速な二次近傍場モデルを利用することで、研究者は、標準的な線形アプローチが失敗する領域（低エネルギーまたは高解像度など、より強い位相勾配が存在する領域）に近い近傍場境界で、法外なコストをかけることなく操作することが可能になる。著者らは、新しい手法のアルゴリズムの複雑さが、並列ハードウェア構成におけるトモグラフィック・バックプロジェクションの背後に隠蔽可能であることを指摘しており、実用的な再構成パイプラインへのさらなる統合を可能にしている。