[论文解读] Full-Wave Iterative Image Reconstruction in Photoacoustic Tomography with Acoustically Inhomogeneous Media

Chao Huang, Kun Wang+3 more|arXiv (Cornell University)|Mar 22, 2013

Photoacoustic and Ultrasonic Imaging参考文献 39被引用 7

一句话总结

本文提出了一种用于声学非均匀介质中光声 computed tomography (PACT) 的全波迭代图像重建方法，基于精确的光声波方程构建离散成像模型。通过采用 k-space 伪谱法实现匹配的前向和反向投影算子，该方法在噪声、数据不完整性和声学非均匀性条件下均能实现鲁棒的重建，显著减少伪影，其有效性在模拟与实验数据上均得到验证。

ABSTRACT

Existing approaches to image reconstruction in photoacoustic computed tomography (PACT) with acoustically heterogeneous media are limited to weakly varying media, are computationally burdensome, and/or cannot effectively mitigate the effects of measurement data incompleteness and noise. In this work, we develop and investigate a discrete imaging model for PACT that is based on the exact photoacoustic (PA) wave equation and facilitates the circumvention of these limitations. A key contribution of the work is the establishment of a procedure to implement a matched forward and backprojection operator pair associated with the discrete imaging model, which permits application of a wide-range of modern image reconstruction algorithms that can mitigate the effects of data incompleteness and noise. The forward and backprojection operators are based on the k-space pseudospectral method for computing numerical solutions to the PA wave equation in the time domain. The developed reconstruction methodology is investigated by use of both computer-simulated and experimental PACT measurement data.

研究动机与目标

• 解决现有 PACT 重建方法在声学非均匀介质中的局限性，包括对弱变化介质的假设以及高计算成本。
• 克服实际 PACT 应用中因数据不完整、噪声和建模误差引起的伪影。
• 构建基于精确光声波方程的离散成像模型，以准确表征声速、密度和衰减参数变化的非均质介质中的波传播。
• 建立匹配的前向与反向投影算子对，以实现复杂介质中迭代重建算法的高效实现。
• 使现代迭代重建技术（如基于原始-对偶混合梯度法的方法）能够应用于具有真实声学非均匀性的全波 PACT 模型。

提出的方法

• 基于精确光声波方程构建离散成像模型，整合空间变化的声速、质量密度和频率幂律衰减。
• 采用 k-space 伪谱法实现前向与反向投影算子，以实现光声波方程的时间域数值求解。
• 开发匹配的算子对（前向、反向投影及其伴随运算），以高效应用迭代重建算法。
• 将算子集成至原始-对偶混合梯度（PDHG）算法中，以在噪声和数据不完整条件下稳健求解逆问题。
• 定义关键算子：$\mathcal{P}_l$ 用于计算离散散度类运算，$\mathcal{P}_c$ 用于非负性强制，$\mathcal{P}_l^T$ 用于伴随运算，$\mathcal{P}_p$ 用于梯度更新中的向量归一化。
• 将该算法应用于计算机模拟与实验 PACT 数据，验证其在真实非均匀介质中的性能。

实验结果

研究问题

• RQ1 基于精确光声波方程的全波迭代重建框架，能否有效减轻声学非均匀性存在时 PACT 中的伪影？
• RQ2 所提出的匹配前向与反向投影算子对在数据不完整与噪声存在时，如何提升重建精度与收敛性？
• RQ3 k-space 伪谱法在多大程度上能够实现 3D 非均匀介质中高效且精确的波传播建模，以支持 PACT？
• RQ4 该迭代重建方法在具有强声学畸变（如颅骨引起的波前畸变）的实验数据上表现如何？
• RQ5 在真实测量条件下，该方法能否在伪影抑制与空间分辨率方面优于传统的时间反演或射线基方法？

主要发现

• 所提出的全波迭代重建方法显著减少了由声学非均匀性（包括颅骨引起的波前畸变）导致的图像伪影，实验结果在经颅 PACT 数据中得到验证。
• 与传统的时间反演和射线基方法相比，该方法在数据不完整或存在噪声的场景下，实现了更高的图像质量与空间分辨率。
• k-space 伪谱法能够实现 3D 中即使存在强声学非均匀性时，光声波方程的精确且稳定的数值求解。
• 匹配的前向与反向投影算子对实现了迭代算法的高效实现，在保持高重建保真度的同时降低了计算负担。
• 基于 PDHG 的算法在稀疏或噪声测量数据条件下仍能稳定收敛，并生成高精度重建结果，其有效性在模拟与实验数据集上均得到验证。
• 由于其基于精确波方程的理论基础，该方法对建模误差与物理参数不确定性（如声速与衰减的变化）表现出强鲁棒性。