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

一句话总结

本文为使用线型探测器的光声计算机断层扫描建立了一个数学框架，将三维图像重建问题简化为在一条线上测量边界数据的二维波动方程求解问题。它推导并严格分析了一种基于傅里叶的解析重建公式，即使在数据可积性方面存在挑战，也能实现快速且稳定的数值实现。

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.

研究动机与目标

• 开发一种使用线型探测器而非点型探测器的光声计算机断层扫描的严格数学模型。
• 解决传统点型探测器近似方法在高带宽、小尺寸实现中所受的技术限制。
• 通过线型探测器测量，将三维光声成像问题简化为二维波动方程初始数据重建问题。
• 提供一种数学上严谨且数值高效的重建公式，以实现实际应用。
• 在真实数据条件下（包括非绝对可积测量）分析重建公式的收敛性和稳定性。

提出的方法

• 将光声过程建模为具有初始压力分布 f(x)（表示吸收能量）的三维波动方程。
• 将线型探测器建模为对压力场沿线段的积分，从而在记录曲线 C_rec 上产生边界测量。
• 通过垂直于探测器线的截面分解三维正向算子，形成一族二维波动方程问题。
• 从线测量中推导出基于傅里叶的二维波动方程初始数据的解析重建公式。
• 应用 C 变换及其逆变换，通过傅里叶分析将测量数据与初始压力分布关联起来。
• 通过证明 C 变换在对称温和分布上是拓扑同构，分析重建公式的收敛性，从而确保其稳定性和可逆性。

实验结果

研究问题

• RQ1 当使用线型探测器时，如何将三维光声成像问题简化为二维波动方程问题？
• RQ2 线型探测器光声计算机断层扫描中正向算子的数学结构是什么？它与点型探测器模型有何不同？
• RQ3 即使测量数据不是绝对可积的，是否可以为线型探测器数据推导出稳定且快速的解析重建公式？
• RQ4 在二维波动方程框架下，C 变换在连接测量数据与初始压力分布方面起什么作用？
• RQ5 所提出的方法如何在真实物理测量条件下确保数值稳定性和收敛性？

主要发现

• 使用线型探测器的三维光声成像问题，可简化为从解在一条线上的限制中重建二维波动方程的初始数据。
• 推导并严格分析了一种基于傅里叶的解析重建公式，即使在测量数据非绝对可积的情况下，也能保证数学上的有效性。
• 证明了 C 变换在对称温和分布上是拓扑同构，从而实现了测量算子的稳定逆运算。
• 通过涉及测量数据傅里叶变换与余弦调制的双重积分，显式给出了逆 C 变换。
• 由于其解析结构以及在频率域与空间域的可分性，该重建公式支持快速的数值实现。
• 通过绕单一轴旋转线型探测器阵列实现三维图像重建，相比三维直接反演，保持了更低的数值复杂度。