Consecutively reconstructing absorption and scattering distributions in turbid media with multiple-illumination photoacoustic tomography

Nonlinear iterative perturbation scheme with simplified spherical harmonics (SP3 ) light propagation model for quantitative photoacoustic tomography

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.

Fluence-matching technique using photoacoustic radiofrequency spectra for improving estimates of oxygen saturation

Photoacoustic (PA) signals encode information about the optical absorption and spatial distribution of absorbing chromophores as well as the light distribution in the medium. The wavelength dependence of the latter affects the accuracy in chromophore quantification, including estimations of oxygen saturation (sO₂) with depth. We propose the use of the ratio of the PA radiofrequency (RF) spectral slopes (SS) at different optical wavelengths to generate frequency filters which can be used to match the fluence profiles across separate images generated with different optical wavelengths. Proof-of-principle experiments were carried on a plastic tube with blood of a known oxygenation inserted into a porcine tissue. The algorithm was tested in-vivo in the hind leg of six CD1 mice, each under three different breathing conditions (100 % O₂, room air and 100 % CO₂). Imaging was done using the VevoLAZR system at the wavelengths 720 and 870 nm. The SS was calculated from the linear fit of the ratio of the photoacoustic RF power spectra at the two wavelengths. An ultrasound frequency filter was designed and applied to each segmented PA signal in the frequency domain and inversely transformed into the time domain to correct for the differences in the fluence profiles at both wavelengths. Linear spectral unmixing was used to estimate sO₂ before and after applying the ultrasound frequency filter. The estimated blood sO₂ in the plastic tube for the porcine tissue experiment improved by 10.3% after applying the frequency filter when compared to the sO₂ measured by a blood gas analyzer. For the in-vivo mouse experiments, the applied sO₂ correction was 2.67, 1.33 and -3.33% for every mm of muscle tissue for mice breathing 100% O₂, room air and 100% CO₂, respectively. The approach presented here provides a new approach for fluence matching that can potentially improve the accuracy of sO₂ estimates by removing the fluence depth dependence at different optical wavelengths.

Towards Quantitative Evaluation of Tissue Absorption Coefficients Using Light Fluence Correction in Optoacoustic Tomography

Optoacoustic tomography is a fast developing imaging modality, combining the high contrast available from optical excitation of tissue with the high resolution and penetration depth of ultrasound detection. Light is subject to both absorption and scattering when traveling through tissue; adequate knowledge of tissue optical properties and hence the spatial fluence distribution is required to create an optoacoustic image that is directly proportional to chromophore concentrations at all depths. Using data from a commercial multispectral optoacoustic tomography (MSOT) system, we implemented an iterative optimization for fluence correction based on a finite-element implementation of the delta-Eddington approximation to the Radiative Transfer Equation (RTE). We demonstrate a linear relationship between the image intensity and absorption coefficients across multiple wavelengths and depths in phantoms. We also demonstrate improved feature visibility and spectral recovery at depth in phantoms and with in vivo measurements, suggesting our approach could in the future enable quantitative extraction of tissue absorption coefficients in biological tissue.

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.