Model-based optoacoustic inversions with incomplete projection data

被引:116
作者
Buehler, Andreas [1 ]
Rosenthal, Amir
Jetzfellner, Thomas
Dima, Alexander
Razansky, Daniel
Ntziachristos, Vasilis
机构
[1] Tech Univ Munich, Inst Biol & Med Imaging, D-85764 Neuherberg, Germany
基金
欧洲研究理事会;
关键词
photoacoustic imaging; limited-view-problem; reconstruction; PHOTOACOUSTIC TOMOGRAPHY; RECONSTRUCTION; REGULARIZATION; ALGORITHM;
D O I
10.1118/1.3556916
中图分类号
R8 [特种医学]; R445 [影像诊断学];
学科分类号
1002 ; 100207 ; 1009 ;
摘要
Purpose: Optoacoustic imaging is an emerging noninvasive imaging modality that can resolve optical contrast through several millimeters to centimeters of tissue with diffraction-limited resolution of ultrasound. Yet, quantified reconstruction of tissue absorption maps requires optoacoustic signals to be collected from as many locations around the object as possible. In many tomographic imaging scenarios, however, only limited-view or partial projection data are available, which has been shown to generate image artifacts and overall loss of quantification accuracy. Methods: In this article, the recently introduced interpolated-matrix-model optoacoustic inversion method is tested in different limited-view scenarios and compared to the standard backprojection algorithm. Both direct (TGSVD) and iterative (PLSQR) regularization methods are proposed to improve the accuracy of image reconstructions with their performance tested on simulated and experimental data. Results: While for full-view tomographic data the model-based inversion has been generally shown to attain higher reconstruction accuracy compared to backprojection algorithms, the incomplete tomographic datasets lead to ill-conditioned forward matrices and, consequently, to error-prone inversions, with strong artifacts following a distinct ripple-type pattern. The proposed regularization techniques are shown to stabilize the inversion and eliminate the artifacts. Conclusions: Overall, it has been determined that the regularized interpolated-matrix-model-based optoacoustic inversions show higher accuracy than reconstructions with the standard backprojection algorithm. Finally, the combination of model-based inversion with PLSQR or TGSVD regularization methods can lead to an accurate reconstruction of limited-projection angle optoacoustic data and practical systems for optoacoustic imaging in many realistic cases where the full-view dataset is unavailable. (C) 2011 American Association of Physicists in Medicine. [DOI: 10.1118/1.3556916]
引用
收藏
页码:1694 / 1704
页数:11
相关论文
共 32 条
[1]   deal. II - A general-purpose object-oriented finite element library [J].
Bangerth, W. ;
Hartmann, R. ;
Kanschat, G. .
ACM TRANSACTIONS ON MATHEMATICAL SOFTWARE, 2007, 33 (04)
[2]  
Golub G., 2013, J HOPKINS STUDIES MA
[3]  
Hansen P. C., 1994, Numerical Algorithms, V6, P1, DOI 10.1007/BF02149761
[4]  
Hansen P C., 1997, Rank-Deficient and Discrete Ill-Posed Problems. Numerical Aspects of Linear Inversion
[5]   TRUNCATED SINGULAR VALUE DECOMPOSITION SOLUTIONS TO DISCRETE ILL-POSED PROBLEMS WITH ILL-DETERMINED NUMERICAL RANK [J].
HANSEN, PC .
SIAM JOURNAL ON SCIENTIFIC AND STATISTICAL COMPUTING, 1990, 11 (03) :503-518
[6]   REGULARIZATION, GSVD AND TRUNCATED GSVD [J].
HANSEN, PC .
BIT, 1989, 29 (03) :491-504
[7]   Regularization tools version 4.0 for matlab 7.3 [J].
Hansen, Per Christian .
NUMERICAL ALGORITHMS, 2007, 46 (02) :189-194
[8]   Spatially varying optical and acoustic property reconstruction using finite-element-based photoacoustic tomography [J].
Jiang, HB ;
Yuan, Z ;
Gu, XJ .
JOURNAL OF THE OPTICAL SOCIETY OF AMERICA A-OPTICS IMAGE SCIENCE AND VISION, 2006, 23 (04) :878-888
[9]   Temporal backward projection of optoacoustic pressure transients using Fourier transform methods [J].
Köstli, KP ;
Frenz, M ;
Bebie, H ;
Weber, HP .
PHYSICS IN MEDICINE AND BIOLOGY, 2001, 46 (07) :1863-1872
[10]  
KRUGER A, PHOTOACOUSTIC