[论文解读] Mathematics of thermoacoustic tomography

Peter Kuchment, Leonid Kunyansky|ArXiv.org|Apr 2, 2007

Photoacoustic and Ultrasonic Imaging参考文献 82被引用 4

一句话总结

本文对热声计算机断层扫描（TAT）进行了全面的数学综述，将成像问题表述为波动方程或球均值算子的反问题。文中建立了唯一性、反演公式及稳定性条件，并针对各种几何结构和函数空间，得出了关于可注入集与范围条件的关键结果。

ABSTRACT

The paper presents a survey of mathematical problems, techniques, and challenges arising in the Thermoacoustic and Photoacoustic Tomography.

研究动机与目标

• 系统分析热声计算机断层扫描（TAT）中的数学挑战，TAT是一种结合电磁与声学信号的混合医学成像方式。
• 解决从测量的声压数据中重建初始能量吸收分布的反问题。
• 研究重建唯一性、稳定性与可逆性的理论条件。
• 对数学模型进行详细处理，包括球均值变换与波动方程表述。
• 识别TAT重建中的开放问题与未解决难题，特别是关于非紧支集与非球形观测表面的问题。

提出的方法

• 将TAT表述为初始条件由能量吸收函数 $ f(x) $ 给出的波动方程反问题。
• 当声速恒定时，将问题简化为球均值变换（球Radon变换），从而通过积分几何进行分析。
• 应用积分几何与微局部分析工具，包括 $ \kappa $-算子与傅里叶分析，研究反演与范围条件。
• 利用时间域的傅里叶变换与空间变换推导重建公式与稳定性估计。
• 分析观测表面几何（平面、柱面、球面、一般闭曲面）对反演公式可行性与结构的影响。
• 考虑部分数据、非紧支集与非恒定声速对反问题可解性与稳定性的影响。

实验结果

研究问题

• RQ1 在何种条件下，初始能量吸收函数 $ f(x) $ 可从表面 $ S $ 上的声压测量中唯一恢复？
• RQ2 当声速恒定时，能否为非球形观测表面推导出显式反演公式？
• RQ3 在一般几何结构下，球均值变换的数据（即范围）的必要与充分条件是什么？
• RQ4 重建的稳定性如何依赖于观测表面的几何形状以及 $ f(x) $ 支集的位置？
• RQ5 位于观测表面之外的 $ f(x) $ 的奇点在多大程度上仍能稳定恢复？

主要发现

• 对于紧支集的 $ f \in L^p(\mathbb{R}^d) $ 且 $ p \leq 2d/(d-1) $，重建唯一性成立，且集合 $ S $ 是此类函数的可注入集。
• 当 $ p > 2d/(d-1) $ 时，唯一性不成立，表明在函数空间可积性上存在一个尖锐的阈值，以保证可注入性。
• 当 $ S $ 为平面、柱面或球面时，存在显式反演公式，但当 $ f $ 的支集超出 $ S $ 时，这些公式的有效性会失效。
• 当 $ f $ 的支集部分位于观测表面 $ S $ 外时，重建的稳定性会丧失，如定理9所示。
• 球均值变换的范围条件对一般闭曲面 $ S $ 是必要的，但预计并非充分，尤其当 $ f $ 的支集位于 $ S $ 外时。
• 球面 $ S $ 存在显式反演公式的事实与 $ \kappa $-算子理论的预期相矛盾，导致一个尚未解决的理论悖论。