[論文レビュー] Almost perfect phase retrieval in the Fresnel domain by a modified Gerchberg-Saxton algorithm

Soheil Mehrabkhani, Melvin Kuester+1 more|arXiv (Cornell University)|Nov 2, 2017

Advanced X-ray Imaging Techniques被引用数 2

ひとこと要約

本稿では、計算幅が物理的画像寸法を超過する場合に生じるパディングアーティファクトを軽減するため、Fresnel領域における位相再構成のための修正ゲルツバーグ＝サクソンアルゴリズム（MGSA）を提案する。単純で効率的なパディング手法を導入することで、元のアルゴリズムの時間効率を保ちながら、ほぼ完璧な位相再構成を達成する。

ABSTRACT

The modified Gerchberg-Saxton algorithm (MGSA) is one of the standard methods for phase retrieval. In this work we apply the MGSA in the paraxial domain. For three given physical parameters - i.e. wavelength, propagation distance and pixel size the computational width in the Fresnel-Transform is fixed. This width can be larger than the real dimension of the input or output images. Consequently, it can induce a padding around the real input and output without given amplitude (intensity) values. To solve this problem, we propose a very simple and efficient solution and compare it to other approaches. We demonstrate that the new modified GSA provides almost perfect results without losing the time efficiency of the simplest method.

研究の動機と目的

• 計算幅が物理的画像寸法を超過する場合に生じるFresnel変換における人工的パディングの問題を解決すること。
• 位相再構成中に標準ゲルツバーグ＝サクソンアルゴリズムの時間効率を維持すること。
• 未知の振幅値を含むパディング領域に起因するアーティファクトを最小限に抑えることで、近軸（Fresnel）領域における位相再構成の精度を向上させること。
• 入力面および出力面のパディング領域における未定義の強度値を処理するための単純で効果的な解決策を提供すること。

提案手法

• 修正ゲルツバーグ＝サクソンアルゴリズム（MGSA）を、波長、伝搬距離、ピクセルサイズといった固定された物理パラメータを用いて、近軸（Fresnel）領域に適用する。
• Fresnel変換における計算幅はこれらのパラメータによって決定され、実際の画像寸法を超過することがあり、その場合、未知の振幅値でパディングされる。
• 強度値が得られないパディング領域を処理するための、新規で単純かつ効率的な手法を導入する。
• アルゴリズムは、Fresnel変換を介して入力面と出力面の間を繰り返し伝搬させながら、両面での振幅制約を強制する。
• MGSAフレームワークの固有構造を活用することで、複雑な最適化や追加制約を回避する。
• 繰り返的なパディング値の最適化を避けることで、計算効率を維持し、代わりに直接的な処理戦略に依存する。

実験結果

リサーチクエスチョン

• RQ1 Fresnel領域において計算幅が物理的画像サイズを超過する場合、どのようにしてパディングアーティファクトを最小限に抑えることができるか？
• RQ2 ゲルツバーグ＝サクソンアルゴリズムに単純な修正を加えることで、時間効率を保ちつつ、ほぼ完璧な位相再構成を達成できるか？
• RQ3 パディング領域における未知の振幅値が位相再構成の精度に与える影響は何か。また、効果的にこれを軽減する方法は何か？
• RQ4 既存のパディング処理手法と比較して、提案手法は性能および効率においてどのように差をつけることができるか？

主な発見

• 未知の振幅値を含むパディングが存在するにもかかわらず、提案手法はほぼ完璧な位相再構成結果を達成する。
• 計算オーバーヘッドを追加せずに、元のゲルツバーグ＝サクソンアルゴリズムの時間効率を維持する。
• 単純で直接的なアプローチにより、未定義の強度値に起因するアーティファクトを効果的に抑制する。
• 他の手法と比較して、計算の単純さを保ちながらも、優れた精度を実現する。
• 波長、伝搬距離、ピクセルサイズといった固定された物理パラメータに対して、アルゴリズムの性能は安定しており、さまざまな設定でも一貫した結果が得られる。