Real-Time Model-Based Inversion in Cross-Sectional Optoacoustic Tomography
Real-Time Model-Based Inversion in Cross-Sectional Optoacoustic Tomography
复制标题
DOI：
10.1109/tmi.2016.2536779
复制
发表时间：
2016-08-01
影响因子：
10.6
通讯作者：
Razansky, Daniel
中科院分区：
文献类型：
--
作者：
Ding, Lu;Dean-Ben, Xose Luis;Razansky, Daniel
分析（封闭形式）反演方案由于其快速重建能力和低内存需求而成为光声断层扫描中图像重建的标准方法。然而，对定量成像和伪影减少的需求导致了更精确的反演方法的发展，这些方法依赖于光声波生成和传播的精确正向建模。通过这种方式，可以合并多个实验因素，例如精确的检测几何形状、换能器的时空响应以及声学异质性。基于模型的反演通常会产生非常大的稀疏矩阵公式，需要大量计算和内存要求的正则化方案来进行图像重建，从而阻碍了它们在实时成像应用中的有效实现。在这里，我们引入了一种新的离散化过程，用于二维光声断层扫描中基于模型的高效重建，该过程允许在图形处理单元（GPU）上以相对较低的数值复杂度并行实现。通过在反演过程的每次迭代中动态计算模型矩阵，新方法导致成像帧速率超过 10 Hz，从而使用基于模型的方法实现实时图像渲染。
Analytical (closed-form) inversion schemes have been the standard approach for image reconstruction in optoacoustic tomography due to their fast reconstruction abilities and low memory requirements. Yet, the need for quantitative imaging and artifact reduction has led to the development of more accurate inversion approaches, which rely on accurate forward modeling of the optoacoustic wave generation and propagation. In this way, multiple experimental factors can be incorporated, such as the exact detection geometry, spatio-temporal response of the transducers, and acoustic heterogeneities. The model-based inversion commonly results in very large sparse matrix formulations that require computationally extensive and memory demanding regularization schemes for image reconstruction, hindering their effective implementation in real-time imaging applications. Herein, we introduce a new discretization procedure for efficient model-based reconstructions in two-dimensional optoacoustic tomography that allows for parallel implementation on a graphics processing unit (GPU) with a relatively low numerical complexity. By on-the-fly calculation of the model matrix in each iteration of the inversion procedure, the new approach results in imaging frame rates exceeding 10 Hz, thus enabling real-time image rendering using the model-based approach.