[論文レビュー] Imaging of buried obstacles in a two-layered medium with phaseless far-field data

Long Li, Jiansheng Yang+2 more|arXiv (Cornell University)|Jan 28, 2020

Numerical methods in inverse problems参考文献 55被引用数 7

ひとこと要約

本稿では、位相なし遠方場データのみを用いて二層媒体内に埋蔵された障害物を再構成する二段階の画像化手法を提案する。位相なし遠方場データにおける並進不変性を解消するため、二つの平面波の重ね合わせに基づく直接的画像化法を用いて小規模な異常を特定し、その後に再帰的ニュートン型反復アルゴリズムを適用して拡張された障害物の形状を再構成する。この手法により、k+ < k− の困難な状況下でも正確な位置特定と形状回復が達成される。

ABSTRACT

The inverse problem we consider is to reconstruct the location and shape of buried obstacles in the lower half-space of an unbounded two-layered medium in two dimensions from phaseless far-field data. A main difficulty of this problem is that the translation invariance property of the modulus of the far field pattern is unavoidable, which is similar to the homogenous background medium case. Based on the idea of using superpositions of two plane waves with different directions as the incident fields, we first develop a direct imaging method to locate the position of small anomalies and give a theoretical analysis of the algorithm. Then a recursive Newton-type iteration algorithm in frequencies is proposed to reconstruct extended obstacles. Finally, numerical experiments are presented to illustrate the feasibility of our algorithms.

研究の動機と目的

• 位相なし遠方場データのみを用いて二層媒体内に埋蔵された障害物の位置特定と再構成を行う逆散乱問題を解決すること。
• 位相なしデータに内在する並進不変性の問題を、二つの入射平面波の重ね合わせを用いることで克服すること。
• 小規模な異常の初期位置特定のための直接的画像化法と、拡張された障害物の高精度な形状再構成のための再帰的ニュートン型反復アルゴリズムを開発すること。
• 波動速度対比（k+ > k− および k+ < k−）の変動を想定した数値実験を通じて、提案手法の有効性を検証すること。

提案手法

• 位相なし遠方場データにおける並進不変性を解消するために、異なる方向を持つ二つの平面波の重ね合わせを入射場として用いる。
• 重ね合わせの原理に基づく直接的画像化法を開発し、サンプリング領域上で画像関数を最大化することで小規模な異常の位置を特定する。
• 周波数領域における再帰的ニュートン型反復アルゴリズムを用いて拡張された障害物を再構成し、直接的画像化の結果を初期推定値として用いる。
• 二層媒体内の散乱をモデル化するために、遠方パターンの漸近的挙動と界面を貫る境界条件を満たすヘルムホルツ方程式を用いる。
• 周波数を複数用いた継続的更新法を実装し、反復の各ステップで波数 k+ と k− を更新することで分解能を向上させる。
• 数値実験において、テスト障害物をパラメトリック曲線（例：楕円、リンゴ型、丸みを帯びた三角形／正方形）で定義する。

実験結果

リサーチクエスチョン

• RQ1 位相なし遠方場データから二層媒体内に埋蔵された小規模な障害物を、重ね合わせ平面波に基づく直接的画像化法で効果的に特定できるか？
• RQ2 位相なし遠方場データのみが利用可能な状況下で、再帰的ニュートン型反復アルゴリズムが拡張された障害物の形状と位置を再構成できるか？
• RQ3 波動速度対比（k+ > k− と k+ < k−）が、提案手法の再構成精度にどのように影響を与えるか？
• RQ4 直接的画像化法が、拡張障害物の反復的再構成に信頼できる初期推定値を提供できる範囲はどの程度か？
• RQ5 4%のノイズを含む遠方場データを用いた数値実験を通じて、提案手法がノイズに頑健であることが示されたか？

主な発見

• 直接的画像化法は、真の障害物位置の近くに局所的最大値を示し、小規模な異常を効果的に特定できた。例えば、リンゴ型障害物では (0.16, −3.69) の位置に明確なピークが観察された。
• 直接的画像化の結果を初期推定値として用いることで、再帰的ニュートン型反復アルゴリズムは、図4.4(c)および4.5(c)に示されるように、拡張された障害物の全形状を高精度に再構成した。
• k+ > k− の場合に比べ、k+ < k− の場合の再構成品質は劣り、特に下部領域の再構成が不正確であった。
• 複数の障害物（例：楕円、丸みを帯びた三角形、丸みを帯びた正方形）に対しては、上部および中央部の再構成が正確であったが、丸みを帯びた正方形の下部はやや不正確であった。
• 4%のノイズが加わった状況下でも、例4～7において一貫した位置特定と形状回復が達成されたため、本手法はノイズに頑健であることが示された。
• 本手法は、完全な位相データが入手できない二層媒体内の位相なし逆散乱問題に対して、実用的な解決策を提供する。