Wavelength-dependent error minimization for quantitative spectroscopic photoacoustic tomography with a ring-array system

PURPOSE

Photoacoustic tomography (PAT) is a non-invasive and high-resolution imaging technique that can provide functional and molecular information from the optical properties of pathological tissues, such as cancer. Spectroscopic PAT (sPAT) is capable of supplying information such as oxygen saturation (sO2), which is an important biological indicator for diseases such as cancer. However, the wavelength dependent nature of sPAT makes it challenging to provide accurate quantitative measurements of tissue oxygenation beyond shallow depths. We have previously reported the utility of combined ultrasound tomography and PAT to achieve optical and acoustic compensated PAT images at a single wavelength and for enhanced PAT images at larger depths. In this work we further explore the utility of the optical and acoustic compensation PAT algorithm to minimize the wavelength dependency in sPAT by showcasing improvements in spectral unmixing.

MATERIALS AND METHODS

Two optically and acoustically characterized heterogenous phantoms were manufactured to test the ability of the system and developed algorithm to minimize the wavelength-dependence driven error in sPAT spectral unmixing. The PA inclusions within each phantom were composed of a mixture of two sulfate dyes, copper sulfate (CuSO4) and nickel sulfate (NiSO4), with known optical spectra. Improvements between uncompensated and optically and acoustically compensated PAT (OAcPAT) were quantified as the relative percent error between the measured results and the ground truth.

RESULTS

The results of our phantom studies demonstrate that OAcPAT can significantly improve the accuracy of sPAT measurements in a heterogenous medium and especially at larger inclusions depths which can reach to up to 12% improvement in measurement errors. This significant improvement can play a vital role in reliability of future in-vivo biomarker quantifications.

CONCLUSIONS

Utilizing UST for model-based optical and acoustic compensation of PAT images was proposed by our group previously. In this work, we further demonstrated the efficacy of the developed algorithm in sPAT by minimizing the error caused by the tissue's optical heterogeneity on improving spectral unmixing, which is a major limiting factor in reliability of sPAT measurements. Such synergistic combination of UST and PAT provides a window of opportunity to achieve bias-free quantitative sPAT measurements, which plays an important role in future pre-clinical and clinical utility of PAT.

Model-based optical and acoustical compensation for photoacoustic tomography of heterogeneous mediums

Photoacoustic tomography (PAT) is a non-invasive, high-resolution imaging modality, capable of providing functional and molecular information of various pathologies, such as cancer. One limitation of PAT is the depth and wavelength dependent optical fluence, which results in reduced PA signal amplitude from deeper tissue regions. These factors can therefore introduce errors into quantitative measurements such as oxygen saturation (sO₂) or the localization and concentration of various chromophores. The variation in the speed-of-sound between different tissues can also lead to distortions in object location and shape. Compensating for these effects allows PAT to be used more quantitatively. We have developed a proof-of-concept algorithm capable of compensating for the heterogeneity in speed-of-sound and depth dependent optical fluence. Speed-of-sound correction was done by using a straight ray-based algorithm for calculating the family of iso-time-of-flight contours between the transducers and every pixel in the imaging grid, while fluence compensation was done by utilizing the graphics processing unit (GPU) accelerated software MCXCL for Monte Carlo modeling of optical fluence variation. This algorithm was tested on a polyvinyl chloride plastisol (PVCP) phantom, which contained cyst mimics and blood inclusions to test the algorithm under relatively heterogeneous conditions. Our results indicate that our PAT algorithm can compensate for the speed-of-sound variation and depth dependent fluence effects within a heterogeneous phantom. The results of this study will pave the way for further development and evaluation of the proposed method in more complex in-vitro and ex-vivo phantoms, as well as compensating for the wavelength-dependent optical fluence in spectroscopic PAT.