[論文レビュー] Nonconvex optimization for optimum retrieval of the transmission matrix of a multimode fiber

Shengfu Cheng, Xuyu Zhang+6 more|arXiv (Cornell University)|Aug 2, 2023

Advanced X-ray Imaging Techniques被引用数 5

ひとこと要約

本論文は、強化勾配降下法と適応的ステップサイズを用いて、強化された非凸最適化手法を提案し、強度測定のみを用いてマルチモードファイバの複素数値伝送関数行列（TM）を高速かつ高精度に回復することを目的としている。低サンプリングレート（例：8）下でも、ホログラフィーを用いたゴールドスタンダードと比較して93.6％の最適焦点化効率を達成し、標準的なGPUアクセラレータ搭載システム上では、全8685×1024のTMをわずか42.3秒で回復可能である。

ABSTRACT

Transmission matrix (TM) allows light control through complex media such as multimode fibers (MMFs), gaining great attention in areas like biophotonics over the past decade. The measurement of a complex-valued TM is highly desired as it supports full modulation of the light field, yet demanding as the holographic setup is usually entailed. Efforts have been taken to retrieve a TM directly from intensity measurements with several representative phase retrieval algorithms, which still see limitations like slow or suboptimum recovery, especially under noisy environment. Here, a modified non-convex optimization approach is proposed. Through numerical evaluations, it shows that the nonconvex method offers an optimum efficiency of focusing with less running time or sampling rate. The comparative test under different signal-to-noise levels further indicates its improved robustness for TM retrieval. Experimentally, the optimum retrieval of the TM of a MMF is collectively validated by multiple groups of single-spot and multi-spot focusing demonstrations. Focus scanning on the working plane of the MMF is also conducted where our method achieves 93.6% efficiency of the gold standard holography method when the sampling rate is 8. Based on the recovered TM, image transmission through the MMF with high fidelity can be realized via another phase retrieval. Thanks to parallel operation and GPU acceleration, the nonconvex approach can retrieve an 8685$ imes$1024 TM (sampling rate=8) with 42.3 s on a regular computer. In brief, the proposed method provides optimum efficiency and fast implementation for TM retrieval, which will facilitate wide applications in deep-tissue optical imaging, manipulation and treatment.

研究の動機と目的

• マルチモードファイバにおける伝送関数行列（TM）回復のための既存の位相再構成アルゴリズムの限界、特に収束が遅く、ノイズ下で最適でない性能を示す問題を克服すること。
• 参照ビームを必要としないホログラフィー装置を避けるために、強度測定のみから正確かつ効率的な複素数値TM回復を可能とすること。
• ノイズ環境下でも高い焦点化効率と耐障害性を維持しながら、サンプリングレートと計算時間を最小限に抑えること。
• 実世界のマルチモードファイバーシステムにおいて、単一および複数スポット焦点化、ならびに画像伝送を通じて、実用的応用の有効性を実証すること。

提案手法

• 緩和平均交互反射（RAF）アルゴリズムに基づく修正非凸最適化フレームワークを採用し、2段階の反復的精錬プロセス（RAF 2-1）を組み込んだ。
• バジライ＝ボーベイン法による適応的ステップサイズを用いて、勾配降下最適化中の収束を加速する。
• 標準的な最急降下法と比較して、より高速な収束を実現するため、非線形共役勾配（NCG）およびL-BFGSといった高度な最適化ソルバを採用する。
• 2段階の反復戦略を導入：まず強度データから位相を推定し、次に勾配ベースの更新により推定値を精錬する。
• GPUアクセラレーションを適用して並列計算を実現し、大規模なTM（例：8685×1024）の処理にかかる実行時間を顕著に短縮する。
• 再構成精度の評価に、$ \| |\mathbf{X}^H \hat{\mathbf{a}}_i| - \mathbf{y}_i \|_2^2 $ で定義される正規化誤差指標を用いる。

実験結果

リサーチクエスチョン

• RQ1 非凸最適化アプローチは、TM回復における焦点化効率と収束速度の面で、既存の凸および非凸位相再構成手法を上回ることができるか？
• RQ2 提案されたRAF 2-1アルゴリズムは、ノイズ下での複素数値伝送関数行列の再構成において、元のRAFアルゴリズムに比べてどのように向上を達成するか？
• RQ3 マルチモードファイバにおける高精度なTM回復を達成するための、計算時間とサンプリングレートの最適なバランスは何か？
• RQ4 強度測定のみで、ホログラフィーを用いたゴールドスタンダード手法と同等の焦点化効率を達成できるか、その程度はどの程度か？
• RQ5 標準的なハードウェアにGPUアクセラレーションを適用した場合、大規模なTM（例：8685×1024）への適用において、本手法のスケーリング性能はどの程度か？

主な発見

• 提案された非凸手法は、サンプリングレート8の下で、ホログラフィーを用いたゴールドスタンダード手法と比較して93.6％の焦点化効率を達成し、従来の位相再構成技術を顕著に上回る。
• 8685×1024のTMに対しては、標準的なコンピュータにGPUアクセラレーションを組み合わせることで、わずか42.3秒で行列を回復可能であり、高い計算効率を示している。
• RAF 2-1はノイズのない条件下で再構成誤差を$10^{-4}$未満に低下させ、元のRAFアルゴリズムと比較して優れた精度を示している。
• NCGを勾配降下ソルバとして採用した場合、単位時間あたりの誤差低下が最も速く、最急降下法およびL-BFGSと比較しても収束速度が優れている。
• 本手法は、さまざまな信号対ノイズ比下でも耐障害性を維持しており、GGS、EKF、SDPベースの手法など、既存のアルゴリズムよりも優れた性能を示している。
• 実験的検証により、回復されたTMを用いたマルチモードファイバを介した高精度な画像伝送が確認され、光学的イメージングおよび波面整形における実用的有用性が裏付けられた。