Gpu-Based Acceleration And Mesh Optimization Of Finite-Element-Method-Based Quantitative Photoacoustic Tomography: A Step Towards Clinical Applications

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.

Publication Date

5-20-2017

Publication Title

Applied Optics

Volume

56

Issue

15

Number of Pages

4426-4432

Document Type

Article

Personal Identifier

scopus

DOI Link

https://doi.org/10.1364/AO.56.004426

Socpus ID

85019933664 (Scopus)

Source API URL

https://api.elsevier.com/content/abstract/scopus_id/85019933664

This document is currently not available here.

Share

COinS