[論文レビュー] A New Numerical Algorithm for Thermoacoustic and Photoacoustic Tomography with Variable Sound Speed

Jianliang QianPlamen Stefanov, Stefanov, Plamen+2 more|arXiv (Cornell University)|Jan 19, 2011

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

ひとこと要約

本稿では、近年の理論的進展に基づいて導出された明示的なノイマン級数の公式を用いて、可変または不連続な音速を伴う熱音響および光音響トモグラフィーのための新しい数値アルゴリズムを提示する。この手法により、非捕捉および捕捉音速、部分境界データ、および区分的連続な媒体といった困難な状況においても、安定的かつ高精度な再構成が可能となり、数値実験でその安定性と正確性が確認された。

ABSTRACT

We present a new algorithm for reconstructing an unknown source in Thermoacoustic and Photoacoustic Tomography based on the recent advances in understanding the theoretical nature of the problem. We work with variable sound speeds that might be also discontinuous across some surface. The latter problem arises in brain imaging. The new algorithm is based on an explicit formula in the form of a Neumann series. We present numerical examples with non-trapping, trapping and piecewise smooth speeds, as well as examples with data on a part of the boundary. These numerical examples demonstrate the robust performance of the new algorithm.

研究の動機と目的

• 可変音速条件下における数値的に安定的かつ効率的な再構成アルゴリズムの開発を目的とする。
• 特に、脳画像診断において軟部組織とくびれ組織の間で顕著に異なる音響特性を示す頭蓋骨組織を含む、不連続な音速の問題に対処することを目的とする。
• 従来の定常音速仮定を超えて、一般の、捕捉または非捕捉の可能性がある音速分布に対しても正確または高精度な再構成を可能にする、既存の逆問題解法の拡張を目的とする。
• 境界の一部でのみ測定が得られるような不完全データの状況（実用的画像診断で一般的な設定）を扱うことを目的とする。
• 非捕捉、捕捉、および区分的連続な音速モデルにおける包括的な数値実験を通じて、アルゴリズムの堅牢性と正確性を検証することを目的とする。

提案手法

• 理論的進展に基づき、可変音速を伴う波動伝搬の逆問題から導出された収束するノイマン級数の形をとる明示的逆公式に基づくアルゴリズムである。
• ノイマン級数の定式化により、境界における音圧測定値から初期源分布の反復的数値再構成が可能となる。
• 滑らかでない音速、たとえば脳画像診断における頭蓋骨-軟部組織界面のような区分的定数または区分的連続な媒体に対しても適用可能である。
• アルゴリズムは、境界距離マップおよび測地線に沿った特異点の伝播を通じて、領域の幾何構造と音速プロファイルを組み込む。
• 数値実装では、有限項数（k）までの級数展開を切り捨て、誤差推定を用いて収束を監視・制御する。
• 全境界および部分境界データ、非捕捉および捕捉音速を含むさまざまな音速モデルにおける設定で、アプローチがテストされた。

実験結果

リサーチクエスチョン

• RQ1 音速が可変または不連続である場合に、熱音響および光音響トモグラフィーに対して正確または高精度な再構成公式を導出可能か？
• RQ2 一部の測定境界がある場合、特に捕捉や観測不完全性により一部の測地線が観測不能となる状況下で、ノイマン級数に基づくアルゴリズムの性能はいかがなものか？
• RQ3 捕捉音速が再構成品質に与える影響は何か？また、十分に大きな時間Tを用いることで、安定な結果を得られるか？
• RQ4 音速の急激な変化が生じる場所（例：人間の頭蓋骨内）のような不連続音速をどのように処理できるか？
• RQ5 音速が定常音速から著しく逸脱しても、不安定または非捕捉な構成において、アルゴリズムはどれほど堅牢かつ正確に保たれるか？

主な発見

• 非捕捉音速と3辺での部分データに対してk=10で再構成誤差6.11%を達成し、不安定性のリスクがあるにもかかわらず高い正確性を示した。
• 非捕捉音速と2辺での部分データに対してk=10で誤差は10.6%に増加し、観測幾何構造が欠落していることによる不可視特異点の存在と整合的であった。
• 捕捉音速と3辺でのデータに対して再構成が著しく劣化し、特に観測されていない境界付近で顕著であった。これは、境界に到達しないまま捕捉される測地線が原因であった。
• 観測時間Tを8.62（安定性閾値T₀ ≈ 2.6）に引き上げたことで、捕捉音速c₂の再構成が著しく改善され、2辺でのデータでも安定した結果が得られた。
• 非捕捉および捕捉条件の下で、シェップ・ローグランのホワイトフィルムおよび変更版「ゼブラ」フィルムの再構成に成功し、特異点が幾何的に不可視な領域でのみ明確なアーチファクトが観察された。
• ノイマン級数は、不連続および捕捉音速を含むすべてのテスト設定で堅牢に収束した。これは、本手法の理論的基盤の妥当性を確認するものであった。