[论文解读] A Reconstruction Procedure for Thermoacoustic Tomography in the Case of Limited Boundary Data

Dustin Steinhauer|ArXiv.org|May 18, 2009

Photoacoustic and Ultrasonic Imaging参考文献 9被引用 8

一句话总结

本文提出了一种用于有限边界数据条件下热声计算机层析成像的显式重构方法，利用傅里叶积分算子和微局部分析，从平面上或曲边界上的测量数据中恢复初始数据的奇异性。关键贡献是一个公式，该公式在区域中无共轭点的假设下，可重构初始数据的波前集，当所有射线在有限时间内到达测量表面时，可实现对所有奇异性的完整恢复。

ABSTRACT

We derive an explicit method for reconstructing singularities of the initial data in a thermoacoustic tomography problem, in the case of variable sound speed and limited boundary data. In order to obtain this explicit formula we assume the metric induced by the sound speed does not have conjugate points inside the region to be observed.

研究动机与目标

• 开发一种当仅有部分边界数据可用时，热声计算机层析成像的显式重构方法。
• 从平面上或边界上的相对开子集上的测量数据中，恢复初始数据的波前集，而非完整的边界数据。
• 建立所有初始数据奇异性均可从有限测量中重构的条件。
• 为可变声速介质中奇异性恢复提供严格的微局部框架。

提出的方法

• 使用具有可变声速 $ c(x) $ 的波动方程，将热声计算机层析成像建模为初始速度为零的柯西问题。
• 应用傅里叶积分算子 $ E_+ $ 和 $ E_- $ 表示解，其典范关系由双特征线带定义。
• 通过求解特征方程 $ \phi_t^2 - c^2(x)|\nabla_x \phi|^2 = 0 $ 的解构造拟解 $ S_{\pm} $，该解在局部时间内有效。
• 通过双特征线传播定义可见性集 $ U $，确保到达测量表面的奇异性被捕捉。
• 使用形式伴随算子 $ S^{*}_{\pm} $ 和伪微分算子 $ R_{j,\pm} $，在光滑算子模下实现测量过程的逆运算。
• 通过单位分解 $ \theta_j $ 的局部重构组合，实现对可见性集 $ U $ 上波前集的全局重构。

实验结果

研究问题

• RQ1 在热声计算机层析成像中，能否从有限边界测量数据中重构初始数据的奇异性？
• RQ2 在区域和声速的何种几何条件下，可见性集可覆盖整个波前集？
• RQ3 如何利用微局部分析和傅里叶积分算子推导有限数据下的显式重构公式？
• RQ4 无共轭点在确保稳定且完整重构中起何种作用？

主要发现

• 本文推导出仅使用平面上或边界上相对开子集上的测量数据，对初始数据 $ f $ 的波前集的显式重构公式。
• 该重构在由 $ c(x) $ 诱导的黎曼度量在感兴趣区域内无共轭点的假设下成立。
• 当所有来自 $ f $ 的波前集的双特征线在时间 $ T $ 内到达测量表面 $ \tilde{\Gamma} $ 时，$ f $ 的整个波前集均可被重构。
• 该方法通过拟解、伴随算子和伪微分算子的组合，实现了波前集的精确重构（模光滑算子）。
• 可见性集 $ U $，定义为 $ U_{j,\pm} $ 的并集，决定了可从数据中恢复的波前集区域。
• 即使测量表面不是平面，而是 $ \partial\Omega $ 的相对开子集，只要区域严格凸且 $ c(x) = 1 $ 在 $ \Omega $ 外部，该结果依然成立。