[論文レビュー] H\"older stability of quantitative photoacoustic tomography based on partial data

Faouzi Triki, Qi Xue|arXiv (Cornell University)|Mar 30, 2021

Photoacoustic and Ultrasonic Imaging参考文献 35被引用数 3

ひとこと要約

本稿は、境界付近の部分的内部データから、定量的光超音波トモグラフィー（PAT）における拡散係数と吸収係数を再構成する際のホルダー安定性を確立する。これは、√D に対するコーシー問題の新規安定性解析に基づくもので、部分領域が境界に近づくにつれて安定性が向上し、ホルダー指数が部分領域に依存しない正の定数に収束する。数値再構成では、滑らかでない媒体や不連続な媒体に対しても相対誤差が低く（3–8.5％）、良好な結果が得られた。

ABSTRACT

We consider the reconstruction of the diffusion and absorption coefficients of the diffusion equation from the internal information of the solution obtained from the first step of the inverse photoacoustic tomography (PAT). In practice, the internal information is only partially provided near the boundary due to the high absorption property of the medium and the limitation of the equipment. Our main contribution is to prove a H\"older stability of the inverse problem in a subregion where the internal information is reliably provided based on the stability estimation of a Cauchy problem satisfied by the diffusion coefficient. The exponent of the H\"older stability converges to a positive constant independent of the subregion as the subregion contracts towards the boundary. Numerical experiments demonstrates that it is possible to locally reconstruct the diffusion and absorption coefficients for smooth and even discontinuous media.

研究の動機と目的

• 組織の吸収や測定制限のため、内部データ H(x) が境界付近でのみ利用可能な定量的 PAT における安定性解析の欠如に対処する。
• 測定境界に近い部分領域で、部分的内部データを用いて拡散係数 D(x) と吸収係数 µ(x) を再構成するホルダー安定性を証明する。
• 従来のコーシー問題の安定性推定が、部分領域が境界に近づくと発散するという限界を克服し、部分領域のサイズに依存しない境界を得る。
• 現実的な部分データ条件下で、滑らかでない媒体および不連続な媒体の局所的再構成の可能性を示す。

提案手法

• 内部データ H(x) = µ(x)u(x) を用いて、√D に対するコーシュー問題を定式化し、Γ ≡ 1 を仮定する。
• 解 uⱼ/u₁ の比を用いて、∇lnσ（σ = D|u₁|²）を Hⱼ/H₁ から回復する方程式系を導出する。
• ノイズ増幅を緩和するための数値微分およびフィルタリングを用いて、√D の楕円型方程式の係数および源項を再構成する。
• 未知領域の欠損データを背景値または平均値で補外することで、既知領域付近への影響を最小限に抑える。
• H(x) がノイズ依存のしきい値を超える部分領域でのみ、√D の楕円型 PDE を安定な数値スキームで解く。
• σ の直接微分を回避するため、係数項を Δ√σ / √σ = ½|∇lnσ|² + ½div(∇lnσ) の恒等式を用いて計算する。

実験結果

リサーチクエスチョン

• RQ1 内部データ H(x) が境界付近でのみ利用可能な場合、PAT における D と µ の逆問題の再構成に対してホルダー安定性を確立できるか？
• RQ2 再構成部分領域が測定境界に近づく際、安定性推定が有界かつ非退化のままであるか？
• RQ3 本手法は、部分データ条件下でも滑らかでない媒体および不連続な媒体を処理できるか？
• RQ4 部分領域が境界に収束する際、安定性推定のホルダー指数はどのように振る舞うか？
• RQ5 実際の応用において、数値的ノイズおよび微分の増幅は再構成品質にどの程度影響を及えるか？

主な発見

• 部分領域が境界に近づくにつれて、D と µ を再構成するホルダー安定性指数が、部分領域に依存しない正の定数に収束する。これは、照明源に近い領域で安定性が向上することを示唆する。
• 滑らかな媒体の数値再構成では、y > 0.2 の領域で D の相対誤差が 3.42%、µ の相対誤差が 3.19% にまで低下し、高い精度が達成された。
• 不連続な媒体に対しても、同領域で D の相対誤差が 8.56%、µ の相対誤差が 4.70% にとどまり、不連続性付近の微分増幅に対しても効果的に機能した。
• 血管モデルの再構成では、D の相対誤差が 4.74%、µ の相対誤差が 2.47% にとどまり、生物学的構造への頑健性が確認された。
• 不連続性の境界で真の極値がフィルタリングによって減少しても、再構成は安定かつ正確に保たれ、数値的ノイズに対して耐性があることが示された。
• 再構成は、H(x) がノイズしきい値を超える部分領域でのみ局在的かつ有効に機能する。これは、部分データ下での実用的妥当性を裏付けている。