1
|
Wang Z, Tao W, Zhao H. Extractor-attention-predictor network for quantitative photoacoustic tomography. PHOTOACOUSTICS 2024; 38:100609. [PMID: 38745884 PMCID: PMC11091525 DOI: 10.1016/j.pacs.2024.100609] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 11/22/2023] [Revised: 04/18/2024] [Accepted: 04/20/2024] [Indexed: 05/16/2024]
Abstract
Quantitative photoacoustic tomography (qPAT) holds great potential in estimating chromophore concentrations, whereas the involved optical inverse problem, aiming to recover absorption coefficient distributions from photoacoustic images, remains challenging. To address this problem, we propose an extractor-attention-predictor network architecture (EAPNet), which employs a contracting-expanding structure to capture contextual information alongside a multilayer perceptron to enhance nonlinear modeling capability. A spatial attention module is introduced to facilitate the utilization of important information. We also use a balanced loss function to prevent network parameter updates from being biased towards specific regions. Our method obtains satisfactory quantitative metrics in simulated and real-world validations. Moreover, it demonstrates superior robustness to target properties and yields reliable results for targets with small size, deep location, or relatively low absorption intensity, indicating its broader applicability. The EAPNet, compared to the conventional UNet, exhibits improved efficiency, which significantly enhances performance while maintaining similar network size and computational complexity.
Collapse
Affiliation(s)
- Zeqi Wang
- School of Sensing Science and Engineering, Shanghai Jiao Tong University, Shanghai, 200240, China
| | - Wei Tao
- School of Sensing Science and Engineering, Shanghai Jiao Tong University, Shanghai, 200240, China
| | - Hui Zhao
- School of Sensing Science and Engineering, Shanghai Jiao Tong University, Shanghai, 200240, China
| |
Collapse
|
2
|
Liang Z, Zhang S, Liang Z, Mo Z, Zhang X, Zhong Y, Chen W, Qi L. Deep learning acceleration of iterative model-based light fluence correction for photoacoustic tomography. PHOTOACOUSTICS 2024; 37:100601. [PMID: 38516295 PMCID: PMC10955667 DOI: 10.1016/j.pacs.2024.100601] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Indexed: 03/23/2024]
Abstract
Photoacoustic tomography (PAT) is a promising imaging technique that can visualize the distribution of chromophores within biological tissue. However, the accuracy of PAT imaging is compromised by light fluence (LF), which hinders the quantification of light absorbers. Currently, model-based iterative methods are used for LF correction, but they require extensive computational resources due to repeated LF estimation based on differential light transport models. To improve LF correction efficiency, we propose to use Fourier neural operator (FNO), a neural network specially designed for estimating partial differential equations, to learn the forward projection of light transport in PAT. Trained using paired finite-element-based LF simulation data, our FNO model replaces the traditional computational heavy LF estimator during iterative correction, such that the correction procedure is considerably accelerated. Simulation and experimental results demonstrate that our method achieves comparable LF correction quality to traditional iterative methods while reducing the correction time by over 30 times.
Collapse
Affiliation(s)
- Zhaoyong Liang
- School of Biomedical Engineering, Southern Medical University, 1023 Shatai Rd., Baiyun District, Guangzhou, Guangdong 510515, China
- Guangdong Provincial Key Laboratory of Medical Image Processing, Southern Medical University, 1023 Shatai Rd., Baiyun District, Guangzhou, Guangdong 510515, China
- Guangdong Province Engineering Laboratory for Medical Imaging and Diagnostic Technology, Southern Medical University, 1023 Shatai Rd., Baiyun District, Guangzhou, Guangdong 510515, China
| | - Shuangyang Zhang
- School of Biomedical Engineering, Southern Medical University, 1023 Shatai Rd., Baiyun District, Guangzhou, Guangdong 510515, China
- Guangdong Provincial Key Laboratory of Medical Image Processing, Southern Medical University, 1023 Shatai Rd., Baiyun District, Guangzhou, Guangdong 510515, China
- Guangdong Province Engineering Laboratory for Medical Imaging and Diagnostic Technology, Southern Medical University, 1023 Shatai Rd., Baiyun District, Guangzhou, Guangdong 510515, China
| | - Zhichao Liang
- School of Biomedical Engineering, Southern Medical University, 1023 Shatai Rd., Baiyun District, Guangzhou, Guangdong 510515, China
- Guangdong Provincial Key Laboratory of Medical Image Processing, Southern Medical University, 1023 Shatai Rd., Baiyun District, Guangzhou, Guangdong 510515, China
- Guangdong Province Engineering Laboratory for Medical Imaging and Diagnostic Technology, Southern Medical University, 1023 Shatai Rd., Baiyun District, Guangzhou, Guangdong 510515, China
| | - Zongxin Mo
- School of Biomedical Engineering, Southern Medical University, 1023 Shatai Rd., Baiyun District, Guangzhou, Guangdong 510515, China
- Guangdong Provincial Key Laboratory of Medical Image Processing, Southern Medical University, 1023 Shatai Rd., Baiyun District, Guangzhou, Guangdong 510515, China
- Guangdong Province Engineering Laboratory for Medical Imaging and Diagnostic Technology, Southern Medical University, 1023 Shatai Rd., Baiyun District, Guangzhou, Guangdong 510515, China
| | - Xiaoming Zhang
- School of Biomedical Engineering, Southern Medical University, 1023 Shatai Rd., Baiyun District, Guangzhou, Guangdong 510515, China
- Guangdong Provincial Key Laboratory of Medical Image Processing, Southern Medical University, 1023 Shatai Rd., Baiyun District, Guangzhou, Guangdong 510515, China
- Guangdong Province Engineering Laboratory for Medical Imaging and Diagnostic Technology, Southern Medical University, 1023 Shatai Rd., Baiyun District, Guangzhou, Guangdong 510515, China
| | - Yutian Zhong
- School of Biomedical Engineering, Southern Medical University, 1023 Shatai Rd., Baiyun District, Guangzhou, Guangdong 510515, China
- Guangdong Provincial Key Laboratory of Medical Image Processing, Southern Medical University, 1023 Shatai Rd., Baiyun District, Guangzhou, Guangdong 510515, China
- Guangdong Province Engineering Laboratory for Medical Imaging and Diagnostic Technology, Southern Medical University, 1023 Shatai Rd., Baiyun District, Guangzhou, Guangdong 510515, China
| | - Wufan Chen
- School of Biomedical Engineering, Southern Medical University, 1023 Shatai Rd., Baiyun District, Guangzhou, Guangdong 510515, China
- Guangdong Provincial Key Laboratory of Medical Image Processing, Southern Medical University, 1023 Shatai Rd., Baiyun District, Guangzhou, Guangdong 510515, China
- Guangdong Province Engineering Laboratory for Medical Imaging and Diagnostic Technology, Southern Medical University, 1023 Shatai Rd., Baiyun District, Guangzhou, Guangdong 510515, China
| | - Li Qi
- School of Biomedical Engineering, Southern Medical University, 1023 Shatai Rd., Baiyun District, Guangzhou, Guangdong 510515, China
- Guangdong Provincial Key Laboratory of Medical Image Processing, Southern Medical University, 1023 Shatai Rd., Baiyun District, Guangzhou, Guangdong 510515, China
- Guangdong Province Engineering Laboratory for Medical Imaging and Diagnostic Technology, Southern Medical University, 1023 Shatai Rd., Baiyun District, Guangzhou, Guangdong 510515, China
| |
Collapse
|
3
|
Zheng S, Yingsa H, Meichen S, Qi M. Quantitative photoacoustic tomography with light fluence compensation based on radiance Monte Carlo model. Phys Med Biol 2023; 68. [PMID: 36821863 DOI: 10.1088/1361-6560/acbe90] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/23/2022] [Accepted: 02/23/2023] [Indexed: 02/25/2023]
Abstract
Objective. Photoacoustic tomography (PAT) is a rapidly evolving imaging modality that provides images with high contrast and spatial resolution showing the optical properties of biological tissues. The photoacoustic pressure is proportional to the product of the optical absorption coefficient and the local light fluence. The essential challenge in reconstructing quantitative images representing spatially varying absorption coefficients is the unknown light fluence. In addition, optical attenuation induces spatial variations in the light fluence, and the heterogeneity of the fluence determines the limits of reconstruction quality and depth.Approach.In this work, a reconstruction enhancement scheme is proposed to compensate for the variation in the light fluence in the absorption coefficient recovery. The inverse problem of the radiance Monte Carlo model describing light transport through the tissue is solved by using an alternating optimization strategy. In the iteration, the absorption coefficients and photon weights are alternately updated.Main results.The method provides highly accurate quantitative images of absorption coefficients in simulations, phantoms, andin vivostudies. The results show that the method has great potential for improving the accuracy of absorption coefficient recovery compared to conventional reconstruction methods that ignore light fluence variations. Comparison with state-of-the-art fluence compensation methods shows significant improvements in root mean square error, normalized mean square absolute distance, and structural similarity metrics.Significance.This method achieves high precision quantitative imaging by compensating for nonuniform light fluence without increasing the complexity and operation of the imaging system.
Collapse
Affiliation(s)
- Sun Zheng
- Department of Electronic and Communication Engineering, North China Electric Power University, Baoding 071003, Hebei, People's Republic of China
- Hebei Key Laboratory of Power Internet of Things Technology, North China Electric Power University, Baoding 071003, Hebei, People's Republic of China
| | - Hou Yingsa
- Department of Electronic and Communication Engineering, North China Electric Power University, Baoding 071003, Hebei, People's Republic of China
| | - Sun Meichen
- Department of Electronic and Communication Engineering, North China Electric Power University, Baoding 071003, Hebei, People's Republic of China
| | - Meng Qi
- Department of Electronic and Communication Engineering, North China Electric Power University, Baoding 071003, Hebei, People's Republic of China
| |
Collapse
|
4
|
Zhang S, Liu J, Liang Z, Ge J, Feng Y, Chen W, Qi L. Pixel-wise reconstruction of tissue absorption coefficients in photoacoustic tomography using a non-segmentation iterative method. PHOTOACOUSTICS 2022; 28:100390. [PMID: 36051488 PMCID: PMC9424605 DOI: 10.1016/j.pacs.2022.100390] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 04/24/2022] [Revised: 06/30/2022] [Accepted: 08/14/2022] [Indexed: 06/15/2023]
Abstract
In Photoacoustic Tomography (PAT), the acquired image represents a light energy deposition map of the imaging object. For quantitative imaging, the PAT image is converted into an absorption coefficient (μ a ) map by dividing the light fluence (LF). Previous methods usually assume a uniform tissueμ a distribution, and consequently degrade the LF correction results. Here, we propose a simple method to reconstruct the pixel-wiseμ a map. Our method is based on a non-segmentation-based iterative algorithm, which alternately optimizes the LF distribution and theμ a map. Using simulation data, as well as phantom and animal data, we implemented our algorithm and compared it to segmentation-based correction methods. The results show that our method can obtain accurate estimation of the LF distribution and therefore improve the image quality and feature visibility of theμ a map. Our method may facilitate efficient calculation of the concentration distributions of endogenous and exogenous agents in vivo.
Collapse
Affiliation(s)
- Shuangyang Zhang
- School of Biomedical Engineering, Southern Medical University, Guangzhou, Guangdong, China
- Guangdong Provincial Key Laboratory of Medical Image Processing, Southern Medical University, Guangzhou, Guangdong, China
- Guangdong Province Engineering Laboratory for Medical Imaging and Diagnostic Technology, Southern Medical University, Guangzhou, Guangdong, China
| | - Jiaming Liu
- School of Biomedical Engineering, Southern Medical University, Guangzhou, Guangdong, China
- Guangdong Provincial Key Laboratory of Medical Image Processing, Southern Medical University, Guangzhou, Guangdong, China
- Guangdong Province Engineering Laboratory for Medical Imaging and Diagnostic Technology, Southern Medical University, Guangzhou, Guangdong, China
| | - Zhichao Liang
- School of Biomedical Engineering, Southern Medical University, Guangzhou, Guangdong, China
- Guangdong Provincial Key Laboratory of Medical Image Processing, Southern Medical University, Guangzhou, Guangdong, China
- Guangdong Province Engineering Laboratory for Medical Imaging and Diagnostic Technology, Southern Medical University, Guangzhou, Guangdong, China
| | - Jia Ge
- School of Biomedical Engineering, Southern Medical University, Guangzhou, Guangdong, China
- Guangdong Provincial Key Laboratory of Medical Image Processing, Southern Medical University, Guangzhou, Guangdong, China
- Guangdong Province Engineering Laboratory for Medical Imaging and Diagnostic Technology, Southern Medical University, Guangzhou, Guangdong, China
| | - Yanqiu Feng
- School of Biomedical Engineering, Southern Medical University, Guangzhou, Guangdong, China
- Guangdong Provincial Key Laboratory of Medical Image Processing, Southern Medical University, Guangzhou, Guangdong, China
- Guangdong Province Engineering Laboratory for Medical Imaging and Diagnostic Technology, Southern Medical University, Guangzhou, Guangdong, China
| | - Wufan Chen
- School of Biomedical Engineering, Southern Medical University, Guangzhou, Guangdong, China
- Guangdong Provincial Key Laboratory of Medical Image Processing, Southern Medical University, Guangzhou, Guangdong, China
- Guangdong Province Engineering Laboratory for Medical Imaging and Diagnostic Technology, Southern Medical University, Guangzhou, Guangdong, China
| | - Li Qi
- School of Biomedical Engineering, Southern Medical University, Guangzhou, Guangdong, China
- Guangdong Provincial Key Laboratory of Medical Image Processing, Southern Medical University, Guangzhou, Guangdong, China
- Guangdong Province Engineering Laboratory for Medical Imaging and Diagnostic Technology, Southern Medical University, Guangzhou, Guangdong, China
| |
Collapse
|
5
|
Hänninen N, Pulkkinen A, Arridge S, Tarvainen T. Adaptive stochastic Gauss-Newton method with optical Monte Carlo for quantitative photoacoustic tomography. JOURNAL OF BIOMEDICAL OPTICS 2022; 27:083013. [PMID: 35396833 PMCID: PMC8993421 DOI: 10.1117/1.jbo.27.8.083013] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Figures] [Subscribe] [Scholar Register] [Received: 12/10/2021] [Accepted: 03/18/2022] [Indexed: 06/14/2023]
Abstract
SIGNIFICANCE The image reconstruction problem in quantitative photoacoustic tomography (QPAT) is an ill-posed inverse problem. Monte Carlo method for light transport can be utilized in solving this image reconstruction problem. AIM The aim was to develop an adaptive image reconstruction method where the number of photon packets in Monte Carlo simulation is varied to achieve a sufficient accuracy with reduced computational burden. APPROACH The image reconstruction problem was formulated as a minimization problem. An adaptive stochastic Gauss-Newton (A-SGN) method combined with Monte Carlo method for light transport was developed. In the algorithm, the number of photon packets used on Gauss-Newton (GN) iteration was varied utilizing a so-called norm test. RESULTS The approach was evaluated with numerical simulations. With the proposed approach, the number of photon packets needed for solving the inverse problem was significantly smaller than in a conventional approach where the number of photon packets was fixed for each GN iteration. CONCLUSIONS The A-SGN method with a norm test can be utilized in QPAT to provide accurate and computationally efficient solutions.
Collapse
Affiliation(s)
- Niko Hänninen
- University of Eastern Finland, Department of Applied Physics, Kuopio, Finland
| | - Aki Pulkkinen
- University of Eastern Finland, Department of Applied Physics, Kuopio, Finland
| | - Simon Arridge
- University College London, Department of Computer Science, London, United Kingdom
| | - Tanja Tarvainen
- University of Eastern Finland, Department of Applied Physics, Kuopio, Finland
- University College London, Department of Computer Science, London, United Kingdom
| |
Collapse
|
6
|
Quantitative Photoacoustic Reconstruction of the Optical Properties of Intervertebral Discs Using a Gradient Descent Scheme. PHOTONICS 2022. [DOI: 10.3390/photonics9020116] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 11/16/2022]
Abstract
The intervertebral discs (IVD) are among the essential organs of the human body, ensuring the mobility of the spine. These organs possess a high proportion of water. However, as the discs age, this content decreases, which can potentially lead to various diseases called degenerative disc diseases. This water content is therefore an important indicator of the well-being of the disc. In this paper, we propose photoacoustic imaging as a means of probing a disc and quantitatively recovering its molecular composition, which should allow concluding on its state. An adjoint-assisted gradient descent scheme is implemented to recover the optical absorption coefficient in the disc, from which, if spectroscopic measurements are performed, the molecular composition can be deduced. The algorithm was tested on synthetic measurements. A realistic numerical phantom was built from magnetic resonance imaging of an actual IVD of a pig. A simplified experiment, with a single laser source, was performed. Results show the feasibility of using photoacoustics imaging to probe IVDs. The influences of exact and approximate formulations of the gradient are studied. The impact of noise on the reconstructions is also evaluated.
Collapse
|
7
|
In-silico investigation towards the non-invasive optical detection of blood lactate. Sci Rep 2021; 11:14274. [PMID: 34253775 PMCID: PMC8275594 DOI: 10.1038/s41598-021-92803-x] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/28/2021] [Accepted: 06/08/2021] [Indexed: 02/06/2023] Open
Abstract
This paper uses Monte Carlo simulations to investigate the interaction of short-wave infrared (SWIR) light with vascular tissue as a step toward the development of a non-invasive optical sensor for measuring blood lactate in humans. The primary focus of this work was to determine the optimal source-detector separation, penetration depth of light at SWIR wavelengths in tissue, and the optimal light power required for reliable detection of lactate. The investigation also focused on determining the non-linear variations in absorbance of lactate at a few select SWIR wavelengths. SWIR photons only penetrated 1.3 mm and did not travel beyond the hypodermal fat layer. The maximum output power was only 2.51% of the input power, demonstrating the need for a highly sensitive detection system. Simulations optimized a source-detector separation of 1 mm at 1684 nm for accurate measurement of lactate in blood.
Collapse
|
8
|
Wang Y, Xu M, Gao F, Kang F, Zhu S. Nonlinear iterative perturbation scheme with simplified spherical harmonics (SP 3 ) light propagation model for quantitative photoacoustic tomography. JOURNAL OF BIOPHOTONICS 2021; 14:e202000446. [PMID: 33576563 DOI: 10.1002/jbio.202000446] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 11/06/2020] [Revised: 01/31/2021] [Accepted: 02/06/2021] [Indexed: 06/12/2023]
Abstract
When using quantitative photoacoustic tomography (q-PAT) reconstruction to recover the optical absorption coefficients of tissue, the commonly used diffusion equation has several limitations in the case of the objects that have small geometries and high-absorption or low-scattering areas. Furthermore, the conventional perturbation reconstruction strategy is unsatisfactory when the target tissue containing large heterogeneous features. We herein present a modified q-PAT implementation that employs the higher-order photon migration model achieving the tradeoff between mathematical rigidity and computational efficiency. Besides, a nonlinear iterative method is proposed to obtain the perturbations of optical absorption considering the updating of the sensitivity matrix in calculating the fluence perturbations. Consequently, the distribution of tissue optical properties can be recovered in a robust way even if the targets with high absorption are included. The proposed approach has been validated by simulation, phantom and in vivo experiments, exhibiting promising performances in image fidelity and quantitative feasibility for practical applications.
Collapse
Affiliation(s)
- Yihan Wang
- School of Life Science and Technology, Xidian University, Xi'an, China
- Engineering Research Center of Molecular and Neuro Imaging of Ministry of Education, Xi'an, China
| | - Menglu Xu
- School of Life Science and Technology, Xidian University, Xi'an, China
- Engineering Research Center of Molecular and Neuro Imaging of Ministry of Education, Xi'an, China
| | - Feng Gao
- College of Precision Instrument and Optoelectronics Engineering, Tianjin University, Tianjin, China
- Tianjin Key Laboratory of Biomedical Detecting Techniques and Instruments, Tianjin University, Tianjin, China
| | - Fei Kang
- Department of Nuclear Medicine, Xijing Hospital, Fourth Military Medical University, Xi'an, China
| | - Shouping Zhu
- School of Life Science and Technology, Xidian University, Xi'an, China
- Engineering Research Center of Molecular and Neuro Imaging of Ministry of Education, Xi'an, China
| |
Collapse
|
9
|
Dantuma M, Kruitwagen S, Ortega-Julia J, Pompe van Meerdervoort RP, Manohar S. Tunable blood oxygenation in the vascular anatomy of a semi-anthropomorphic photoacoustic breast phantom. JOURNAL OF BIOMEDICAL OPTICS 2021; 26:JBO-200370RR. [PMID: 33728828 PMCID: PMC7961914 DOI: 10.1117/1.jbo.26.3.036003] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 11/25/2020] [Accepted: 02/19/2021] [Indexed: 05/21/2023]
Abstract
SIGNIFICANCE Recovering accurate oxygenation estimations in the breast with quantitative photoacoustic tomography (QPAT) is not straightforward. Accurate light fluence models are required, but the unknown ground truth of the breast makes it difficult to validate them. Phantoms are often used for the validation, but most reported phantoms have a simple architecture. Fluence models developed in these simplistic objects are not accurate for application on the complex tissues of the breast. AIM We present a sophisticated breast phantom platform for photoacoustic (PA) and ultrasound (US) imaging in general, and specifically for QPAT. The breast phantom is semi-anthropomorphic in distribution of optical and acoustic properties and contains wall-less channels with blood. APPROACH 3D printing approaches are used to develop the solid 3D breast phantom from custom polyvinyl chloride plastisol formulations and additives for replicating the tissue optical and acoustic properties. A flow circuit was developed to flush the channels with bovine blood with a controlled oxygen saturation level. To showcase the phantom's functionality, PA measurements were performed on the phantom with two oxygenation levels. Image reconstructions with and without fluence compensation from Monte Carlo simulations were analyzed for the accuracy of oxygen saturation estimations. RESULTS We present design aspects of the phantom, demonstrate how it is developed, and present its breast-like appearance in PA and US imaging. The oxygen saturations were estimated in two regions of interest with and without using the fluence models. The fluence compensation positively influenced the SO2 estimations in all cases and confirmed that highly accurate fluence models are required to minimize estimation errors. CONCLUSIONS This phantom allows studies to be performed in PA in carefully controlled laboratory settings to validate approaches to recover both qualitative and quantitative features sought after in in-vivo studies. We believe that testing with phantoms of this complexity can streamline the transition of new PA technologies from the laboratory to studies in the clinic.
Collapse
Affiliation(s)
- Maura Dantuma
- University of Twente, Multi-Modality Medical Imaging, Techmed Centre, Enschede, The Netherlands
- Address all correspondence to Maura Dantuma,
| | - Saskia Kruitwagen
- University of Twente, Multi-Modality Medical Imaging, Techmed Centre, Enschede, The Netherlands
- Medisch Spectrum Twente, Enschede, The Netherlands
| | - Javier Ortega-Julia
- University of Twente, Multi-Modality Medical Imaging, Techmed Centre, Enschede, The Netherlands
| | | | - Srirang Manohar
- University of Twente, Multi-Modality Medical Imaging, Techmed Centre, Enschede, The Netherlands
| |
Collapse
|
10
|
Leino AA, Lunttila T, Mozumder M, Pulkkinen A, Tarvainen T. Perturbation Monte Carlo Method for Quantitative Photoacoustic Tomography. IEEE TRANSACTIONS ON MEDICAL IMAGING 2020; 39:2985-2995. [PMID: 32217473 DOI: 10.1109/tmi.2020.2983129] [Citation(s) in RCA: 8] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/10/2023]
Abstract
Quantitative photoacoustic tomography aims at estimating optical parameters from photoacoustic images that are formed utilizing the photoacoustic effect caused by the absorption of an externally introduced light pulse. This optical parameter estimation is an ill-posed inverse problem, and thus it is sensitive to measurement and modeling errors. In this work, we propose a novel way to solve the inverse problem of quantitative photoacoustic tomography based on the perturbation Monte Carlo method. Monte Carlo method for light propagation is a stochastic approach for simulating photon trajectories in a medium with scattering particles. It is widely accepted as an accurate method to simulate light propagation in tissues. Furthermore, it is numerically robust and easy to implement. Perturbation Monte Carlo maintains this robustness and enables forming gradients for the solution of the inverse problem. We validate the method and apply it in the framework of Bayesian inverse problems. The simulations show that the perturbation Monte Carlo method can be used to estimate spatial distributions of both absorption and scattering parameters simultaneously. These estimates are qualitatively good and quantitatively accurate also in parameter scales that are realistic for biological tissues.
Collapse
|
11
|
Zhou X, Akhlaghi N, Wear KA, Garra BS, Pfefer TJ, Vogt WC. Evaluation of Fluence Correction Algorithms in Multispectral Photoacoustic Imaging. PHOTOACOUSTICS 2020; 19:100181. [PMID: 32405456 PMCID: PMC7210453 DOI: 10.1016/j.pacs.2020.100181] [Citation(s) in RCA: 16] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Received: 01/10/2020] [Revised: 04/03/2020] [Accepted: 04/07/2020] [Indexed: 05/07/2023]
Abstract
Multispectral photoacoustic imaging (MPAI) is a promising emerging diagnostic technology, but fluence artifacts can degrade device performance. Our goal was to develop well-validated phantom-based test methods for evaluating and comparing MPAI fluence correction algorithms, including a heuristic diffusion approximation, Monte Carlo simulations, and an algorithm we developed based on novel application of the diffusion dipole model (DDM). Phantoms simulated a range of breast-mimicking optical properties and contained channels filled with chromophore solutions (ink, hemoglobin, or copper sulfate) or connected to a previously developed blood flow circuit providing tunable oxygen saturation (SO2). The DDM algorithm achieved similar spectral recovery and SO2 measurement accuracy to Monte Carlo-based corrections with lower computational cost, potentially providing an accurate, real-time correction approach. Algorithms were sensitive to optical property uncertainty, but error was minimized by matching phantom albedo. The developed test methods may provide a foundation for standardized assessment of MPAI fluence correction algorithm performance.
Collapse
Affiliation(s)
- Xuewen Zhou
- Fischell Department of Bioengineering, University of Maryland, College Park, MD, 02742, United States
- Center for Devices and Radiological Health, U.S. Food and Drug Administration, Silver Spring, MD 20993, United States
| | - Nima Akhlaghi
- Center for Devices and Radiological Health, U.S. Food and Drug Administration, Silver Spring, MD 20993, United States
| | - Keith A. Wear
- Center for Devices and Radiological Health, U.S. Food and Drug Administration, Silver Spring, MD 20993, United States
| | - Brian S. Garra
- Center for Devices and Radiological Health, U.S. Food and Drug Administration, Silver Spring, MD 20993, United States
| | - T. Joshua Pfefer
- Center for Devices and Radiological Health, U.S. Food and Drug Administration, Silver Spring, MD 20993, United States
| | - William C. Vogt
- Center for Devices and Radiological Health, U.S. Food and Drug Administration, Silver Spring, MD 20993, United States
- Corresponding author.
| |
Collapse
|
12
|
Macdonald CM, Arridge S, Powell S. Efficient inversion strategies for estimating optical properties with Monte Carlo radiative transport models. JOURNAL OF BIOMEDICAL OPTICS 2020; 25:JBO-200101R. [PMID: 32798354 PMCID: PMC7426481 DOI: 10.1117/1.jbo.25.8.085002] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Figures] [Subscribe] [Scholar Register] [Received: 04/14/2020] [Accepted: 07/23/2020] [Indexed: 06/11/2023]
Abstract
SIGNIFICANCE Indirect imaging problems in biomedical optics generally require repeated evaluation of forward models of radiative transport, for which Monte Carlo is accurate yet computationally costly. We develop an approach to reduce this bottleneck, which has significant implications for quantitative tomographic imaging in a variety of medical and industrial applications. AIM Our aim is to enable computationally efficient image reconstruction in (hybrid) diffuse optical modalities using stochastic forward models. APPROACH Using Monte Carlo, we compute a fully stochastic gradient of an objective function for a given imaging problem. Leveraging techniques from the machine learning community, we then adaptively control the accuracy of this gradient throughout the iterative inversion scheme to substantially reduce computational resources at each step. RESULTS For example problems of quantitative photoacoustic tomography and ultrasound-modulated optical tomography, we demonstrate that solutions are attainable using a total computational expense that is comparable to (or less than) that which is required for a single high-accuracy forward run of the same Monte Carlo model. CONCLUSIONS This approach demonstrates significant computational savings when approaching the full nonlinear inverse problem of optical property estimation using stochastic methods.
Collapse
Affiliation(s)
- Callum M. Macdonald
- University College London, Department of Medical Physics and Biomedical Engineering, London, United Kingdom
| | - Simon Arridge
- University College London, Department of Computer Science, London, United Kingdom
| | - Samuel Powell
- University of Nottingham, Faculty of Engineering, Nottingham, United Kingdom
| |
Collapse
|
13
|
Mahmoodkalayeh S, Zarei M, Ansari MA, Kratkiewicz K, Ranjbaran M, Manwar R, Avanaki K. Improving vascular imaging with co-planar mutually guided photoacoustic and diffuse optical tomography: a simulation study. BIOMEDICAL OPTICS EXPRESS 2020; 11:4333-4347. [PMID: 32923047 PMCID: PMC7449743 DOI: 10.1364/boe.385017] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 12/06/2019] [Revised: 07/01/2020] [Accepted: 07/01/2020] [Indexed: 05/04/2023]
Abstract
Diffuse optical tomography (DOT) and photoacoustic tomography (PAT) are functional imaging modalities that provide absorption coefficient maps of the tissue. Spatial resolution of DOT is relatively low due to light scattering characteristics of the tissue. On the other hand, although PAT can resolve regions of different absorptions with a high spatial resolution, measuring the absolute value of optical absorptions using PAT is challenging due to unknown light fluence distribution in the tissue. Development of image guidance techniques using a priori information of imaging target structure has been shown to increase the accuracy of DOT. PAT is one such method that can be used as a complementary modality to serve as a guide for DOT image reconstruction. On the other hand, estimated fluence map provided by DOT can be used to quantitatively correct PAT images. In this study we introduce a mutually-guided imaging system for fast and simultaneous optical and photoacoustic measurements of tissue absorption map, where DOT is guided by the PAT image and vice versa. Using the obtained absorption map of the tissue, we then estimate the tissue scattering map. We conducted this study using a series of simulations on digital phantoms and demonstrated the effectiveness of the proposed method.
Collapse
Affiliation(s)
| | - Mehrdad Zarei
- Laser and Plasma Research Institute, Shahid Beheshti University, Tehran 19839 69411, Iran
| | - Mohammad Ali Ansari
- Laser and Plasma Research Institute, Shahid Beheshti University, Tehran 19839 69411, Iran
| | - Karl Kratkiewicz
- Wayne State University, Bioengineering Department, Detroit, Michigan 48201, USA
| | - Mohsen Ranjbaran
- Department of Physics, University of Isfahan, Isfahan 81746-73441, Iran
| | - Rayyan Manwar
- Wayne State University, Bioengineering Department, Detroit, Michigan 48201, USA
| | - Kamran Avanaki
- Richard and Loan Hill Department of Bioengineering, University of Illinois at Chicago, Chicago, Illinois 60607, USA
- Department of Dermatology, University of Illinois at Chicago, Chicago, Illinois 60607, USA
| |
Collapse
|
14
|
Pandey PK, Bharadwaj J, Naik N, Aggrawal HO. One-step fluorescence photoacoustic tomography with the optical radiative transfer model. JOURNAL OF THE OPTICAL SOCIETY OF AMERICA. A, OPTICS, IMAGE SCIENCE, AND VISION 2020; 37:1175-1192. [PMID: 32609678 DOI: 10.1364/josaa.389476] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 01/29/2020] [Accepted: 05/13/2020] [Indexed: 06/11/2023]
Abstract
We present adjoint-based Jacobian as well as gradient evaluations and corresponding reconstruction schemes to solve the fully nonlinear, optical radiative transfer modeled one-step fluorescence photoacoustic tomographic (FPAT) problem, which aims to reconstruct the map of absorption coefficient of the exogenous fluorophore from boundary photoacoustic data. The radiative transport equation (RTE) and frequency-domain photoacoustic equation have been employed to model light and photoacoustic wave propagation, respectively. Levenberg-Marquardt and Broyden-Fletcher-Goldfarb-Shanno reconstruction schemes have been used corresponding to the evaluated Jacobians and gradients, respectively. Numerical reconstructions obtained from the two schemes have been validated for scattering-dominant as well as nonscattering-dominant media in 2D. To the best of our knowledge, these are the first one-step FPAT reconstruction results in literature based on the optical RTE model.
Collapse
|
15
|
Asllanaj F, Addoum A. Simultaneous reconstruction of absorption, scattering and anisotropy factor distributions in quantitative photoacoustic tomography. Biomed Phys Eng Express 2020; 6:045010. [DOI: 10.1088/2057-1976/ab90a0] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/12/2022]
|
16
|
Pandey PK, Gottam O, Naik N, Pradhan A. Gradient-based one-step fluorescence photoacoustic tomography. APPLIED OPTICS 2020; 59:4357-4366. [PMID: 32400412 DOI: 10.1364/ao.382879] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 11/08/2019] [Accepted: 04/11/2020] [Indexed: 06/11/2023]
Abstract
We propose a gradient-based scheme to solve the fluorescence photoacoustic tomographic (FPAT) problem in a fully nonlinear one-step setting, which aims to reconstruct the map of the absorption coefficient of an exogenous fluorophore from boundary photoacoustic pressure data. Adjoint-based gradient evaluation is presented for the FPAT problem in a frequency-domain photoacoustic equation setting. Numerical validations of the resulting Broyden-Fletcher-Goldfarb-Shanno (BFGS) reconstruction scheme are carried out in two dimensions for full- as well as limited-data test cases, and the results are compared with existing Jacobian-based one-step FPAT reconstructions. The reasonably comparable results of the one-step gradient- and Jacobian-based FPAT reconstruction schemes, coupled with the significant computational savings of the former, potentially set up the one-step gradient-based schemes as an advantageous method of choice for FPAT reconstructions. Further reconstruction studies carried out using quantitative photoacoustic tomography (QPAT)-based chromophore reconstructions as inputs to the FPAT inversions show a robustness of fluorophore absorption coefficient reconstructions to the QPAT-obtained inputs.
Collapse
|
17
|
Buchmann J, Kaplan B, Powell S, Prohaska S, Laufer J. Quantitative PA tomography of high resolution 3-D images: Experimental validation in a tissue phantom. PHOTOACOUSTICS 2020; 17:100157. [PMID: 31956487 PMCID: PMC6961715 DOI: 10.1016/j.pacs.2019.100157] [Citation(s) in RCA: 20] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Received: 10/23/2019] [Revised: 11/29/2019] [Accepted: 12/05/2019] [Indexed: 05/18/2023]
Abstract
Quantitative photoacoustic tomography aims to recover the spatial distribution of absolute chromophore concentrations and their ratios from deep tissue, high-resolution images. In this study, a model-based inversion scheme based on a Monte-Carlo light transport model is experimentally validated on 3-D multispectral images of a tissue phantom acquired using an all-optical scanner with a planar detection geometry. A calibrated absorber allowed scaling of the measured data during the inversion, while an acoustic correction method was employed to compensate the effects of limited view detection. Chromophore- and fluence-dependent step sizes and Adam optimization were implemented to achieve rapid convergence. High resolution 3-D maps of absolute concentrations and their ratios were recovered with high accuracy. Potential applications of this method include quantitative functional and molecular photoacoustic tomography of deep tissue in preclinical and clinical studies.
Collapse
Affiliation(s)
- Jens Buchmann
- Institut für Physik, Martin-Luther-Universität Halle-Wittenberg, von-Danckelmann-Platz 3, 06120 Halle (Saale), Germany
- Institut für Optik und Atomare Physik, Technische Universität Berlin, Straße des 17, Juni 135, 10623 Berlin, Germany
| | - Bernhard Kaplan
- Visual Data Analysis, Zuse Institute Berlin, Takustr. 7, 14195 Berlin, Germany
| | - Samuel Powell
- Optics and Photonics Group, Faculty of Engineering, University of Nottingham, University Park, Nottingham NG7 2RD, United Kingdom
| | - Steffen Prohaska
- Visual Data Analysis, Zuse Institute Berlin, Takustr. 7, 14195 Berlin, Germany
| | - Jan Laufer
- Institut für Physik, Martin-Luther-Universität Halle-Wittenberg, von-Danckelmann-Platz 3, 06120 Halle (Saale), Germany
- Corresponding author.
| |
Collapse
|
18
|
Buchmann J, Kaplan BA, Powell S, Prohaska S, Laufer J. Three-dimensional quantitative photoacoustic tomography using an adjoint radiance Monte Carlo model and gradient descent. JOURNAL OF BIOMEDICAL OPTICS 2019; 24:1-13. [PMID: 31172727 PMCID: PMC6977014 DOI: 10.1117/1.jbo.24.6.066001] [Citation(s) in RCA: 16] [Impact Index Per Article: 3.2] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 01/09/2019] [Accepted: 04/24/2019] [Indexed: 05/18/2023]
Abstract
Quantitative photoacoustic tomography aims to recover maps of the local concentrations of tissue chromophores from multispectral images. While model-based inversion schemes are promising approaches, major challenges to their practical implementation include the unknown fluence distribution and the scale of the inverse problem. We describe an inversion scheme based on a radiance Monte Carlo model and an adjoint-assisted gradient optimization that incorporates fluence-dependent step sizes and adaptive moment estimation. The inversion is shown to recover absolute chromophore concentrations, blood oxygen saturation, and the Grüneisen parameter from in silico three-dimensional phantom images for different radiance approximations. The scattering coefficient is assumed to be homogeneous and known a priori.
Collapse
Affiliation(s)
- Jens Buchmann
- Technische Universität Berlin, Institut für Optik und Atomare Physik, Berlin, Germany
| | | | - Samuel Powell
- King’s College London, Biomedical Engineering and Imaging Sciences, Becket House, London, United Kingdom
| | | | - Jan Laufer
- Martin-Luther-Universität Halle-Wittenberg, Institut für Physik, Halle (Saale), Germany
| |
Collapse
|
19
|
Yao J, Wang LV. Recent progress in photoacoustic molecular imaging. Curr Opin Chem Biol 2018; 45:104-112. [PMID: 29631120 PMCID: PMC6076847 DOI: 10.1016/j.cbpa.2018.03.016] [Citation(s) in RCA: 57] [Impact Index Per Article: 9.5] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/13/2018] [Revised: 03/24/2018] [Accepted: 03/27/2018] [Indexed: 01/08/2023]
Abstract
By acoustically detecting the optical absorption contrast, photoacoustic (PA) tomography (PAT) has broken the penetration limits of traditional high-resolution optical imaging. Through spectroscopic analysis of the target's optical absorption, PAT can identify a wealth of endogenous and exogenous molecules and thus is inherently capable of molecular imaging with high sensitivity. PAT's molecular sensitivity is uniquely accompanied by non-ionizing radiation, high spatial resolution, and deep penetration in biological tissues, which other optical imaging modalities cannot achieve yet. In this concise review, we summarize the most recent technological advancements in PA molecular imaging and highlight the novel molecular probes specifically made for PAT in deep tissues. We conclude with a brief discussion of the opportunities for future advancements.
Collapse
Affiliation(s)
- Junjie Yao
- Department of Biomedical Engineering, Duke University, Durham, NC 27708, USA.
| | - Lihong V Wang
- Caltech Optical Imaging Laboratory, Andrew and Peggy Cherng Department of Medical Engineering, Department of Electrical Engineering, California Institute of Technology, Pasadena, CA 91125, USA.
| |
Collapse
|
20
|
Nykänen O, Pulkkinen A, Tarvainen T. Quantitative photoacoustic tomography augmented with surface light measurements. BIOMEDICAL OPTICS EXPRESS 2017; 8:4380-4395. [PMID: 29082072 PMCID: PMC5654787 DOI: 10.1364/boe.8.004380] [Citation(s) in RCA: 6] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Received: 07/06/2017] [Revised: 08/18/2017] [Accepted: 08/26/2017] [Indexed: 05/11/2023]
Abstract
Quantitative photoacoustic tomography is an imaging modality in which distributions of optical parameters inside tissue are estimated from photoacoustic images. This optical parameter estimation is an ill-posed problem and it needs to be approached in the framework of inverse problems. In this work, utilising surface light measurements in quantitative photoacoustic tomography is studied. Estimation of absorption and scattering is formulated as a minimisation problem utilising both internal quantitative photoacoustic data and surface light data. The image reconstruction problem is studied with two-dimensional numerical simulations in various imaging situations using the diffusion approximation as the model for light propagation. The results show that quantitative photoacoustic tomography augmented with surface light data can improve both absorption and scattering estimates when compared to the conventional quantitative photoacoustic tomography.
Collapse
Affiliation(s)
- Olli Nykänen
- Department of Applied Physics, University of Eastern Finland, P.O. Box 1627, 70211 Kuopio,
Finland
| | - Aki Pulkkinen
- Department of Applied Physics, University of Eastern Finland, P.O. Box 1627, 70211 Kuopio,
Finland
| | - Tanja Tarvainen
- Department of Applied Physics, University of Eastern Finland, P.O. Box 1627, 70211 Kuopio,
Finland
- Department of Computer Science, University College London, Gower Street, London WC1E 6BT,
UK
| |
Collapse
|
21
|
Carbone NA, Iriarte DI, Pomarico JA. GPU accelerated Monte Carlo simulation of light propagation in inhomogeneous fluorescent turbid media: application to whole field CW imaging. Biomed Phys Eng Express 2017. [DOI: 10.1088/2057-1976/aa7b8f] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/03/2023]
|