[論文レビュー] Explicit Reconstructions in QPAT, QTAT, TE, and MRE

Guillaume Bal|arXiv (Cornell University)|Feb 14, 2012

Medical Imaging Techniques and Applications参考文献 36被引用数 5

ひとこと要約

本稿では、初期波再構成から導かれる内部関数型を活用することで、定量的光音響（QPAT）、熱音響（QTAT）、一時的エラストグラフィ（TE）、磁気共鳴エラストグラフィ（MRE）の明示的で安定した再構成手順を提示する。第二階テンソル（例：拡散の異方性）および複素係数（例：吸収および減衰）が、データの滑らかさと導関数の損失に依存する安定性推定のもとで一意に再構成可能であることが示され、特にスカラ成分が全テンソルよりも安定に再構成されることを明らかにした。

ABSTRACT

Photo-acoustic Tomography (PAT) and Thermo-acoustic Tomography (TAT) are medical imaging modalities that combine the high contrast of radiative properties of tissues with the high resolution of ultrasound. In both modalities, a first step concerns the reconstruction of the radiation-induced source of ultrasound. Transient Elastography (TE) and Magnetic Resonance Elastography (MRE) combine the high elastic contrast of tissues with the high resolution of ultrasound and magnetic resonance, respectively. In both modalities, a first step concerns the reconstruction of the elastic displacement. The result of this first step, which is not considered in this paper, is the availability of internal functionals of the unknown tissue properties. All imaging modalities are recast as the reconstruction of parameters in elliptic equations from knowledge of solutions to such equations. This paper provides a characterization of the parameters that may or may not be reconstructed from such internal functionals. We provide explicit reconstruction procedures and indicate how stable they are with respect to errors in the available measurements. The modalities PAT, TAT, TE, and MRE allow us to reconstruct high-contrast optical and elastic properties of tissues with the high resolution of ultrasound or magnetic resonance imaging. They provide a means to reconstruct second-order tensors modeling tissue anisotropy as well as complex-valued coefficients modeling absorbing and dissipative effects.

研究の動機と目的

• 初期逆波問題を解いた後に得られる内部関数型から、QPAT、QTAT、TE、MREにおけるどのパラメータが一意に再構成可能であるかを特定すること。
• 内部関数型の比を用いて、光学的および弾性組織特性の明示的再構成アルゴリズムを開発すること。
• 測定誤差に対する再構成の安定性を定量化すること、特に再構成プロセスにおける導関数の損失を分析すること。
• スカラおよびテンソル値の係数（複素数値および異方的パrameterを含む）に一般化可能なフレームワークを提供することで、先行研究を拡張すること。
• 再構成安定性が未知係数の構造（スカラ対全テンソル）に依存することを示し、スカラ成分がより安定である理由を明らかにすること。

提案手法

• QPAT、QTAT、TE、MREを、波解の比から得られる内部関数型を有する楕円型方程式の逆問題として再定式化する。
• 未知解を消去し、未知係数を分離するために、内部関数型の比 $ H_j / H_1 $ を用いる。
• テンソル係数を再構成する場合には関数型比の二階微分を、スカラ成分を再構成する場合には一階微分を用いて再構成式を構築する。
• SobolevおよびH"olderノルムを用いて安定性解析を行い、データから再構成係数への誤差伝搬を定量化する。
• ノイズのあるデータ処理のため、反復的固定点法およびローパスフィルタリングを用い、関数型の微分処理における誤差増幅を抑制する。
• 拡散テンソル $ a $ のスカラ成分 $ B $ の再構成において、積分のおかげで1階分の導関数が得られ、全テンソル再構成よりも安定性が向上することを示した。

実験結果

リサーチクエスチョン

• RQ1 初期逆波問題を解いた後に得られる内部関数型から、QPAT、QTAT、TE、MREにおけるどのパラメータが一意に再構成可能であるか。
• RQ2 再構成の安定性がデータの滑らかさおよび再構成プロセスにおける導関数の損失数にどのように依存するか。
• RQ3 複素数値および異方的係数（例：拡散テンソル）が、同一の関数型枠組みから明示的かつ安定に再構成可能か。
• RQ4 係数構造（スカラ対全テンソル）が再構成安定性に与える影響は何か。なぜスカラ再構成がより安定なのか。
• RQ5 ノイズのあるデータをどのように処理すれば、誤差を増幅させることなく内部関数型の正確な微分が可能になるか。

主な発見

• QPATにおける吸収係数 $ \sigma $ の再構成は、$ \|\delta\sigma\|_{\infty} \leq C\|\delta H_1\|_{\infty} $ と表される一様安定性を示し、測定誤差下でも安定である。
• スカラ拡散係数 $ a $ の再構成では1階分の導関数損失が生じるが、全テンソル再構成では2階分の導関数損失が生じるため、スカラ成分の方が安定性が高い。
• $ \Gamma $ が既知であるとき、安定性推定 $ \|\delta\hat{a}\|_{C^0} + \|\delta\sigma\|_{C^0} + \|\delta B\|_{C^1} \leq C\|\delta H\|_{C^2} $ が成り立ち、スカラ部 $ B $ の安定性が向上していることが示された。
• 拡散テンソルのスカラ成分 $ B $ の再構成は、積分のおかげで1階分の導関数が得られ、$ \sigma $ よりも安定性が向上する。
• 数値実装により、ノイズのあるデータの微分処理がローパスフィルタリングを用いることで安定化され、再構成精度が維持されることを確認した。
• 本手法により、複素数値係数および異方的テンソルの明示的かつ安定な再構成が可能となり、先行研究をより広い種類の組織特性に拡張した。