[論文レビュー] Mathematical Challenges Arising in Thermoacoustic Tomography with Line Detectors

Markus Haltmeier, T C FIDLER|ArXiv.org|Oct 4, 2006

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

ひとこと要約

本論文は、線形検出器を用いた熱音響トモグラフィのための数学的枠組みを確立し、3次元画像再構成問題を、直線上で測定された境界データを用いた2次元波動方程式への還元する。本稿では、データの可積分性に課題があるにもかかわらず、高速かつ安定な数値実装が可能なフーリエに基づく解析的再構成公式を導出し、厳密に分析する。

ABSTRACT

Thermoacoustic computed tomography (thermoacoustic CT) has the potential to become a mayor non-invasive medical imaging method. In this paper we derive a general mathematical framework of a novel measuring setup introduced in [P. Burgholzer, C. Hofer, G. Paltauf, M. Haltmeier, and O. Scherzer, "Thermoacoustic tomography with integrating area and line detectors", IEEE Transactions on Ultrasonics, Ferroelectrics, and Frequency Control, 52 (2005)], that uses line shaped detectors instead of the usual point like ones. We show that the three dimensional thermoacoustic imaging problem reduces to the mathematical problem of reconstructing the initial data of the two dimensional wave equation from boundary measurements of its solution. We derive and analyze an analytic reconstruction formula which allows for fast numerical implementation.

研究の動機と目的

• 点検出器の代わりに線形検出器を用いた熱音響トモグラフィのための厳密な数学的モデルの構築。
• 高帯域幅・小型化実装において技術的制約を受ける従来の点検出器近似の限界を克服すること。
• 線形検出器の測定値を用いて、3次元熱音響イメージング問題を2次元波動方程式の初期データ再構成問題に還元すること。
• 実用的実装に適した数学的に妥当で数値的に効率的な再構成公式の提供。
• 非絶対可積分測定データを含む現実的なデータ条件下での再構成公式の収束性と安定性の分析。

提案手法

• 熱音響プロセスを、吸収エネルギーを表す初期圧力分布 f(x) を持つ3次元波動方程式として定式化する。
• 線形検出器を線分に沿った圧力場の積分としてモデル化し、記録曲線 C_rec における境界測定値を導出する。
• 検出器線に垂直な断面で領域をスライスすることで、3次元前方演算子を2次元波動方程式問題の族に分解する。
• 線形測定値から2次元波動方程式初期データを再構成するためのフーリエに基づく解析的再構成公式を導出する。
• C変換およびその逆変換を適用し、フーリエ解析を介して測定データと初期圧力分布を関連付ける。
• C変換が対称的緩和分布上で位相的同型写像であることを証明することで、測定演算子の安定な逆写像が可能であることを分析する。

実験結果

リサーチクエスチョン

• RQ1 線形検出器を用いた場合、3次元熱音響イメージング問題はどのように2次元波動方程式問題に還元できるか？
• RQ2 線形検出器を用いた熱音響トモグラフィにおける前方演算子の数学的構造は何か？また、点検出器モデルとはどのように異なるか？
• RQ3 測定データが絶対可積分でない場合でも、安定的かつ高速な解析的再構成公式を導出できるか？
• RQ4 C変換は2次元波動方程式フレームワークにおいて、測定データと初期圧力分布を結ぶ役割を果たすか？
• RQ5 提案手法は、現実的な物理的測定条件下で、数値的安定性と収束性をどのように確保するか？

主な発見

• 線形検出器を用いた3次元熱音響イメージング問題は、解を直線上に制限した2次元波動方程式の初期データ再構成問題に還元される。
• フーリエに基づく解析的再構成公式が導出され、厳密に分析されており、測定データが絶対可積分でない場合でも数学的妥当性が保証される。
• C変換が対称的緩和分布上で位相的同型写像であることが証明され、測定演算子の安定な逆写像が可能となる。
• 逆C変換は、測定データのフーリエ変換とコサイン変調を含む二重積分を用いて明示的に与えられる。
• 再構成公式の解析的構造と周波数領域および空間領域における分離可能性のおかげで、高速な数値実装が可能である。
• 本手法は、線形検出器アレイを単一軸の周りに回転させることで3次元画像再構成を達成し、3次元直接逆問題に比べて低い数値的複雑性を維持する。