Abstract
Finite element method (FEM)-based time-domain quantitative photoacoustic tomography (TD-qPAT) is a powerful approach, as it provides highly accurate quantitative imaging capability by recovering absolute tissue absorption coefficients for functional imaging. However, this approach is extremely computationally demanding, and requires days for the reconstruction of one set of images, making it impractical to be used in clinical applications, where a large amount of data needs to be processed in a limited time scale. To address this challenge, here we present a graphic processing unit (GPU)-based parallelization method to accelerate the image reconstruction using FEM-based TD-qPAT. In addition, to further optimize FEM-based TD-qPAT reconstruction, an adaptive meshing technique, along with mesh density optimization, is adopted. Phantom experimental data are used in our study to evaluate the GPU-based TD-qPAT algorithm, as well as the adaptive meshing technique. The results show that our new approach can considerably reduce the computation time by at least 136-fold over the current central processing unit (CPU)-based algorithm. The quality of image reconstruction is also improved significantly when adaptive meshing and mesh density optimization are applied.
© 2017 Optical Society of America
Full Article | PDF ArticleMore Like This
Kuan Peng, Ling He, Ziqiang Zhu, Jingtian Tang, and Jiaying Xiao
Appl. Opt. 52(34) 8270-8279 (2013)
Bo Zhang, Xiang Yang, Fei Yang, Xin Yang, Chenghu Qin, Dong Han, Xibo Ma, Kai Liu, and Jie Tian
Opt. Express 18(19) 20201-20214 (2010)
Lei Yao and Huabei Jiang
Appl. Opt. 50(25) 5031-5041 (2011)