[論文レビュー] 3D Wave-Equation-Based Finite-Frequency Tomography for Ultrasound Computed Tomography

Naiara Korta Martiartu, Christian Boehm+1 more|ArXiv.org|Aug 9, 2019

Seismic Imaging and Inversion Techniques参考文献 46被引用数 6

ひとこと要約

本論文は、周波数依存の伝搬時間感受性と体積散乱効果を考慮することで、計算的に効率的で高周波数の3次元イメージングを可能にする、3次元波動方程式に基づく有限周波数トモグラフィー手法を、超音波画像診断（USCT）に導入する。校正データを用いた解析的相互相関伝搬時間感受性と、ジャコビアンの1次元パrameterizationを用いることで、GPUアクセラレーションを実現したマトリクスフリーな逆問題解法が可能となり、2次元スライス単位の取得から高精度な3次元再構成が可能となり、解像度が向上し、縦方向のアーチファクトが低減される。

ABSTRACT

Ultrasound Computed Tomography (USCT) has great potential for 3D quantitative imaging of acoustic breast tissue properties. Typical devices include high-frequency transducers, which makes tomography techniques based on numerical wave propagation simulations computationally challenging, especially in 3D. Therefore, despite the finite-frequency nature of ultrasonic waves, ray-theoretical approaches to transmission tomography are still widely used. This work introduces finite-frequency traveltime tomography to medical ultrasound. In addition to being computationally tractable for 3D imaging at high frequencies, the method has two main advantages: (1) It correctly accounts for the frequency dependence and volumetric sensitivity of traveltime measurements, which are related to off-ray-path scattering and diffraction. (2) It naturally enables out-of-plane imaging and the construction of 3D images from 2D slice-by-slice acquisition systems. Our method rests on the availability of calibration data in water, used to linearize the forward problem and to provide analytical expressions of cross-correlation traveltime sensitivity. As a consequence of the finite frequency content, sensitivity is distributed in multiple Fresnel volumes, thereby providing out-of-plane sensitivity. To improve computational efficiency, we develop a memory-efficient implementation by encoding the Jacobian operator with a 1D parameterization, which allows us to extend the method to large-scale domains. We validate our tomographic approach using lab measurements collected with a 2D setup of transducers and using a cylindrically symmetric phantom. We then demonstrate its applicability for 3D reconstructions by simulating a slice-by-slice acquisition systems using the same dataset.

研究の動機と目的

• レイトラムトモグラフィーの限界を克服するため、回折や散乱を含む有限周波数波動効果を考慮する。
• 2次元スライス単位のUSCT取得システムにおいて、有限周波数カーネルを用いて縦方向の感受性を捉えることで、真正の3次元再構成を実現する。
• 大規模3次元問題における計算効率とメモリ最適化を実現するため、感受性カーネルの1次元パrameterizationを用いた、計算に高効率でスケーラブルな逆問題フレームワークを開発する。
• 無限に細いレイトラム近似を物理的に現実的な周波数依存感受性カーネルに置き換えることで、空間解像度を向上させ、逆問題のアーチファクトを低減する。
• 水中の校正データを用いて事前にジャコビアンを解析的に計算することで、臨床的実用性があり、高速かつ高精度な逆問題手法を提供する。

提案手法

• 波動方程式に基づく前方問題を線形化するためのボーン近似を採用し、速度摂動量のための大規模な線形逆問題を定式化する。
• 校正データを用いた解析的相互相関伝搬時間感受性の計算により、高価な数値的波動伝搬シミュレーションを回避する。
• 1次元パrameterizationにより、送信機受信機の幾何配置に依存しない完全なジャコビアン作用素を符号化し、マトリクス・ベクトル積の即時計算を可能にする。
• マトリクスフリーなアプローチは、並列処理に強く、GPUアーキテクチャに最適化されており、メモリ使用量を顕著に削減し、大規模3次元アプリケーションに適している。
• アドジョイント技法を用いてジャコビアンを効率的に計算し、有限周波数性を反映した正確な感受性カーネルを保証する。
• 空間解像度の評価には、点スプール関数（PSF）推定を用い、空間解像度とパラメータ間のトレードオフを定量的に評価する。

実験結果

リサーチクエスチョン

• RQ1 有限周波数トモグラフィーは、高周波数3次元USCTに適応可能であり、レイトラム法を上回る解像度を実現できるか？
• RQ2 有限周波数感受性カーネルは、2次元スライス単位のUSCTシステムにおいて、縦方向のイメージングをどのように可能にするか？
• RQ3 ジャコビアン作用素の1次元パrameterizationは、大規模3次元USCT問題における計算効率とスケーラビリティを達成できるか？
• RQ4 校正データからの解析的感受性計算は、逆問題の実行時間をどの程度短縮し、臨床的実用性を向上させるか？
• RQ5 有限周波数カーネルは、レイトラム理論近似と比較して、解像度とアーチファクト低減の点でどの程度優れているか？

主な発見

• 本手法は、体積的有限周波数カーネルを用いて縦方向の感受性を捉えることで、2次元スライス単位のデータから3次元音響特性を正確に再構成する。
• 複数のフレネルゾーンに分散する有限周波数感受性カーネルは、物理的に達成可能な解像度に自然に制限をかけるため、無限に細いレイトラム近似によるアーチファクトを低減する。
• 1次元パrameterizationにより、マトリクスフリーかつGPUアクセラレーションされた逆問題が実現され、メモリ使用量が顕著に削減され、大規模3次元ドメインへのスケーラブルな応用が可能になる。
• 水中の校正データを用いたジャコビアンの解析的計算により、繰り返し波動伝搬シミュレーションの必要がなくなり、解に至るまでの時間は著しく短縮される。
• 円対称なフィジカルフィジカルを用いたラボデータでの検証により、本手法の3次元再構成における正確性と頑健性が確認された。
• PSFを用いた解像度分析により、従来のレイトラム法と比較して、空間解像度が向上し、速度異常値の制約がより厳密に評価されていることが示された。