Speed-of-sound compensated photoacoustic tomography for accurate imaging

被引:119
作者
Jose, Jithin [1 ]
Willemink, Rene G. H. [2 ]
Steenbergen, Wiendelt [1 ]
Slump, C. H. [2 ]
van Leeuwen, Ton G. [1 ,3 ]
Manohar, Srirang [1 ]
机构
[1] Univ Twente, Biomed Photon Imaging Grp, MIRA Inst Biomed Technol & Tech Med, NL-7500 AE Enschede, Netherlands
[2] Univ Twente, Signals & Syst Grp, MIRA Inst Biomed Technol & Tech Med, NL-7500 AE Enschede, Netherlands
[3] Univ Amsterdam, Acad Med Ctr, NL-1100 DE Amsterdam, Netherlands
关键词
photoacoustic tomography; photoacoustic imaging; acoustic velocity; acoustic attenuation; ultrasound computed tomography; ULTRASOUND COMPUTED-TOMOGRAPHY; THERMOACOUSTIC TOMOGRAPHY; RECONSTRUCTION; VASCULATURE; MAMMOSCOPE; ALGORITHM;
D O I
10.1118/1.4764911
中图分类号
R8 [特种医学]; R445 [影像诊断学];
学科分类号
1002 ; 100207 ; 1009 ;
摘要
Purpose: In most photoacoustic (PA) tomographic reconstructions, variations in speed-of-sound (SOS) of the subject are neglected under the assumption of acoustic homogeneity. Biological tissue with spatially heterogeneous SOS cannot be accurately reconstructed under this assumption. The authors present experimental and image reconstruction methods with which 2D SOS distributions can be accurately acquired and reconstructed, and with which the SOS map can be used subsequently to reconstruct highly accurate PA tomograms. Methods: The authors begin with a 2D iterative reconstruction approach in an ultrasound transmission tomography setting, which uses ray refracted paths instead of straight ray paths to recover accurate SOS images of the subject. Subsequently, they use the SOS distribution in a new 2D iterative PA reconstruction approach, where refraction of rays originating from PA sources is accounted for in accurately retrieving the distribution of these sources. Both the SOS reconstruction and SOS-compensated PA reconstruction methods utilize the Eikonal equation to model acoustic wavefront propagation. The equation is solved using a high accuracy fast marching method. Results: The authors validated the new reconstruction algorithms using numerical phantoms. For experiments they utilized the recently introduced PER-PACT method which can be used to simultaneously acquire SOS and PA data from subjects. Conclusions: It is first confirmed that it is important to take SOS inhomogeneities into account in high resolution PA tomography. The iterative reconstruction algorithms, that model acoustic refractive effects, in reconstructing SOS distributions, and subsequently using these distributions to correct PA tomograms, yield artifact-free highly accurate images. The approach of using the hybrid measurement method and the new reconstruction algorithms is successful in substantially improving the quality of PA images with a minimization of blurring and artifacts. (C) 2012 American Association of Physicists in Medicine. [http://dx.doi.org/10.1118/1.4764911]
引用
收藏
页码:7262 / 7271
页数:10
相关论文
共 54 条
[1]   ANALYSIS OF BOUNDED VARIATION PENALTY METHODS FOR ILL-POSED PROBLEMS [J].
ACAR, R ;
VOGEL, CR .
INVERSE PROBLEMS, 1994, 10 (06) :1217-1229
[2]   Image reconstruction in thermoacoustic tomography with compensation for acoustic heterogeneities [J].
Anastasio, MA ;
Zhang, J ;
Pan, XC .
MEDICAL IMAGING 2005: ULTRASONIC IMAGING AND SIGNAL PROCESSING, 2005, 5750 :298-304
[3]  
Bamber J. C., 2007, ENCY ACOUSTICS, V4
[4]   Biomedical photoacoustic imaging [J].
Beard, Paul .
INTERFACE FOCUS, 2011, 1 (04) :602-631
[5]   Statistical optoacoustic image reconstruction using a-priori knowledge on the location of acoustic distortions [J].
Dean-Ben, X. Luis ;
Ntziachristos, Vasilis ;
Razansky, Daniel .
APPLIED PHYSICS LETTERS, 2011, 98 (17)
[6]   Laser optoacoustic imaging system for detection of breast cancer [J].
Ermilov, Sergey A. ;
Khamapirad, Tuenchit ;
Conjusteau, Andre ;
Leonard, Morton H. ;
Lacewell, Ron ;
Mehta, Ketan ;
Miller, Tom ;
Oraevsky, Alexander A. .
JOURNAL OF BIOMEDICAL OPTICS, 2009, 14 (02)
[7]   A real-time photoacoustic tomography system for small animals [J].
Gamelin, John ;
Maurudis, Anastasios ;
Aguirre, Andres ;
Huang, Fei ;
Guo, Puyun ;
Wang, Lihong V. ;
Zhu, Quing .
OPTICS EXPRESS, 2009, 17 (13) :10489-10498
[8]  
Hansen P. C., 1998, RANK DEFICIENT DISCR, P243
[9]   THE USE OF THE L-CURVE IN THE REGULARIZATION OF DISCRETE III-POSED PROBLEMS [J].
HANSEN, PC ;
OLEARY, DP .
SIAM JOURNAL ON SCIENTIFIC COMPUTING, 1993, 14 (06) :1487-1503
[10]   Visualizing breast cancer using the Twente photoacoustic mammoscope: What do we learn from twelve new patient measurements? [J].
Heijblom, M. ;
Piras, D. ;
Xia, W. ;
van Hespen, J. C. G. ;
Klaase, J. M. ;
van den Engh, F. M. ;
van Leeuwen, T. G. ;
Steenbergen, W. ;
Manohar, S. .
OPTICS EXPRESS, 2012, 20 (11) :11582-11597