[论文解读] 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

一句话总结

本文提出一种基于3D波方程的有限频率层析成像方法，用于超声 computed tomography (USCT)，通过考虑频率相关的时延敏感性和体积散射效应，实现了计算高效、高频的3D成像。通过使用从水中校准数据中导出的解析互相关时延敏感性，并采用雅可比矩阵的一维参数化，该方法实现了GPU加速的无矩阵反演，从而能够从2D逐片采集中实现高分辨率、低剖面外伪影的准确3D重建。

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.

研究动机与目标

• 通过考虑衍射和散射等有限频率波效应，克服射线基层析成像在高频USCT中的局限性。
• 通过有限频率核函数捕捉剖面外敏感性，实现从2D逐片USCT采集系统中真正实现3D重建。
• 通过敏感性核函数的一维参数化，开发一种计算高效、内存优化的反演框架，适用于大规模3D问题。
• 通过用物理上合理的、频率相关的敏感性核函数替代无限细射线近似，提高空间分辨率并减少反演伪影。
• 通过使用水中校准数据预先解析计算雅可比矩阵，提供一种临床可行、快速且准确的反演方法。

提出的方法

• 该方法采用Born近似线性化基于波方程的正向问题，从而将速度扰动建模为大规模线性反问题。
• 利用水中校准数据，通过解析方法计算互相关时延敏感性，避免了昂贵的数值波传播模拟。
• 一维参数化独立于发射器-接收器几何结构，编码了完整的雅可比算子，支持即时计算矩阵-向量乘积。
• 无矩阵方法高度可并行化，且针对GPU架构进行了优化，显著降低内存使用量，支持大规模3D应用。
• 采用伴随技术高效计算雅可比矩阵，确保敏感性核函数准确反映超声波的有限频率特性。
• 通过点扩散函数（PSF）估计评估分辨率，以量化空间分辨率和参数间权衡。

实验结果

研究问题

• RQ1 有限频率层析成像能否适配高频3D USCT，以实现超越射线基方法的分辨率？
• RQ2 有限频率敏感性核函数如何在2D逐片USCT系统中实现剖面外成像？
• RQ3 雅可比算子的一维参数化能否在大规模3D USCT问题中实现计算效率与可扩展性？
• RQ4 从水中校准数据中解析计算敏感性是否能显著缩短反演时间并提升临床可行性？
• RQ5 在分辨率和伪影抑制方面，有限频率核函数与射线理论近似相比表现如何？

主要发现

• 该方法通过在多个Fresnel区域分布的有限频率敏感性核函数，成功从2D逐片数据中重建出3D声学特性，捕捉了剖面外敏感性。
• 有限频率敏感性核函数分布在多个Fresnel区域，自然地将分辨率限制在物理可实现水平，减少了由无限细射线近似引起的伪影。
• 一维参数化实现了无矩阵、GPU加速的反演，显著降低内存使用量，使大规模3D区域应用成为可能。
• 通过水中校准数据预先解析计算雅可比矩阵，消除了重复波传播模拟的需求，大幅缩短求解时间。
• 在具有圆柱对称性的幻影实验数据上的验证，证实了该方法在3D重建中的准确性和鲁棒性。
• 通过PSF进行的分辨率分析表明，与传统射线基方法相比，该方法实现了更高的空间分辨率，并对速度异常的约束更优。