Modeling of Errors Due to Uncertainties in Ultrasound Sensor Locations in Photoacoustic Tomography

Photoacoustic tomography is an imaging modality based on the photoacoustic effect caused by the absorption of an externally introduced light pulse. In the inverse problem of photoacoustic tomography, the initial pressure generated through the photoacoustic effect is estimated from a measured photoacoustic time-series utilizing a forward model for ultrasound propagation. Due to the ill-posedness of the inverse problem, errors in the forward model or measurements can result in significant errors in the solution of the inverse problem. In this work, we study modeling of errors caused by uncertainties in ultrasound sensor locations in photoacoustic tomography using a Bayesian framework. The approach is evaluated with simulated and experimental data. The results indicate that the inverse problem of photoacoustic tomography is sensitive even to small uncertainties in sensor locations. Furthermore, these uncertainties can lead to significant errors in the estimates and reduction of the quality of the photoacoustic images. In this work, we show that the errors due to uncertainties in ultrasound sensor locations can be modeled and compensated using Bayesian approximation error modeling.

Modelling of errors due to speed of sound variations in photoacoustic tomography using a Bayesian framework

Inverse problem of estimating initial pressure in photoacoustic tomography is ill-posed and thus sensitive to errors in modelling and measurements. In practical experiments, accurate knowledge of the speed of sound of the imaged target is commonly not available, and therefore an approximate speed of sound is used in the computational model. This can result in errors in the solution of the inverse problem that can appear as artefacts in the reconstructed images. In this paper, the inverse problem of photoacoustic tomography is approached in a Bayesian framework. Errors due to uncertainties in the speed of sound are modelled using Bayesian approximation error modelling. Estimation of the initial pressure distribution together with information on the reliability of these estimates are considered. The approach was studied using numerical simulations. The results show that uncertainties in the speed of sound can cause significant errors in the solution of the inverse problem. However, modelling of these uncertainties improves the accuracy of the solution.

Image Reconstruction with Reliability Assessment in Quantitative Photoacoustic Tomography

Quantitative photoacoustic tomography is a novel imaging method which aims to reconstruct optical parameters of an imaged target based on initial pressure distribution, which can be obtained from ultrasound measurements. In this paper, a method for reconstructing the optical parameters in a Bayesian framework is presented. In addition, evaluating the credibility of the estimates is studied. Furthermore, a Bayesian approximation error method is utilized to compensate the modeling errors caused by coarse discretization of the forward model. The reconstruction method and the reliability of the credibility estimates are investigated with two-dimensional numerical simulations. The results suggest that the Bayesian approach can be used to obtain accurate estimates of the optical parameters and the credibility estimates of these parameters. Furthermore, the Bayesian approximation error method can be used to compensate for the modeling errors caused by a coarse discretization, which can be used to reduce the computational costs of the reconstruction procedure. In addition, taking the modeling errors into account can increase the reliability of the credibility estimates.

Quantitative photoacoustic tomography augmented with surface light measurements

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.

Direct Estimation of Optical Parameters From Photoacoustic Time Series in Quantitative Photoacoustic Tomography

Estimation of optical absorption and scattering of a target is an inverse problem associated with quantitative photoacoustic tomography. Conventionally, the problem is expressed as two folded. First, images of initial pressure distribution created by absorption of a light pulse are formed based on acoustic boundary measurements. Then, the optical properties are determined based on these photoacoustic images. The optical stage of the inverse problem can thus suffer from, for example, artefacts caused by the acoustic stage. These could be caused by imperfections in the acoustic measurement setting, of which an example is a limited view acoustic measurement geometry. In this work, the forward model of quantitative photoacoustic tomography is treated as a coupled acoustic and optical model and the inverse problem is solved by using a Bayesian approach. Spatial distribution of the optical properties of the imaged target are estimated directly from the photoacoustic time series in varying acoustic detection and optical illumination configurations. It is numerically demonstrated, that estimation of optical properties of the imaged target is feasible in limited view acoustic detection setting.

Characterisation of a phantom for multiwavelength quantitative photoacoustic imaging

Quantitative photoacoustic imaging (qPAI) has the potential to provide high- resolution in vivo images of chromophore concentration, which may be indicative of tissue function and pathology. Many strategies have been proposed recently for extracting quantitative information, but many have not been experimentally verified. Experimental phantom-based validation studies can be used to test the robustness and accuracy of such algorithms in order to ensure reliable in vivo application is possible. The phantoms used in such studies must have well-characterised optical and acoustic properties similar to tissue, and be versatile and stable. Polyvinyl chloride plastisol (PVCP) has been suggested as a phantom for quality control and system evaluation. By characterising its multiwavelength optical properties, broadband acoustic properties and thermoelastic behaviour, this paper examines its potential as a phantom for qPAI studies too. PVCP's acoustic properties were assessed for various formulations, as well as its intrinsic optical absorption, and scattering with added TiO2, over a range of wavelengths from 400-2000 nm. To change the absorption coefficient, pigment-based chromophores that are stable during the phantom fabrication process, were used. These yielded unique spectra analogous to tissue chromophores and linear with concentration. At the high peak powers typically used in photoacoustic imaging, nonlinear optical absorption was observed. The Grüneisen parameter was measured to be [Formula: see text]  =  1.01  ±  0.05, larger than typically found in tissue, though useful for increased PA signal. Single and multiwavelength 3D PA imaging of various fabricated PVCP phantoms were demonstrated.