A Comparative Study of Four Total Variational Regularization Reconstruction Algorithms for Sparse-View Photoacoustic Imaging.

Photoacoustic imaging (PAI) is a new nonionizing, noninvasive biomedical imaging technology that has been employed to reconstruct the light absorption characteristics of biological tissues. The latest developments in compressed sensing (CS) technology have shown that it is possible to accurately reconstruct PAI images from sparse data, which can greatly reduce scanning time. This study focuses on the comparative analysis of different CS-based total variation regularization reconstruction algorithms, aimed at finding a method suitable for PAI image reconstruction. The performance of four total variation regularization algorithms is evaluated through the reconstruction experiment of sparse numerical simulation signal and agar phantom signal data. The evaluation parameters include the signal-to-noise ratio and normalized mean absolute error of the PAI image and the CPU time. The comparative results demonstrate that the TVAL3 algorithm can well balance the quality and efficiency of the reconstruction. The results of this study can provide some useful guidance for the development of the PAI sparse reconstruction algorithm.

Evaluation of visible NIR-I and NIR-II light penetration for photoacoustic imaging in rat organs.

In this study, we evaluate the penetration capability of light in visible, near-infrared-I (NIR-I) and near-infrared-II (NIR-II) optical windows for photoacoustic macroscale imaging inside 9 biological tissues with three typical penetration depths. An acoustic resolution photoacoustic microscopy is designed to guarantee the consistent experiment conditions except excitation wavelength. Experimental results show that short NIR-II (1000-1150 nm) shows the best performance inside kidney, spleen and liver tissues at all depths, while NIR-I (700-1000 nm) works better for muscle, stomach, heart and brain tissues, especially in deep imaging. This study proposes the optimal selection of illumination wavelengths for photoacoustic macroscale imaging in rat organs, which enables the best signal-to-noise ratio (SNR) of the observed target.

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.

Spectral unmixing techniques for optoacoustic imaging of tissue pathophysiology.

A key feature of optoacoustic imaging is the ability to illuminate tissue at multiple wavelengths and therefore record images with a spectral dimension. While optoacoustic images at single wavelengths reveal morphological features, in analogy to ultrasound imaging or X-ray imaging, spectral imaging concedes sensing of intrinsic chromophores and externally administered agents that can reveal physiological, cellular and subcellular functions. Nevertheless, identification of spectral moieties within images obtained at multiple wavelengths requires spectral unmixing techniques, which present a unique mathematical problem given the three-dimensional nature of the optoacoustic images. Herein we discuss progress with spectral unmixing techniques developed for multispectral optoacoustic tomography. We explain how different techniques are required for accurate sensing of intrinsic tissue chromophores such as oxygenated and deoxygenated haemoglobin versus extrinsically administered photo-absorbing agents and nanoparticles. Finally, we review recent developments that allow accurate quantification of blood oxygen saturation (sO2) by transforming and solving the sO2 estimation problem from the spatial to the spectral domain.This article is part of the themed issue 'Challenges for chemistry in molecular imaging'.

Improved optoacoustic microscopy through three-dimensional spatial impulse response synthetic aperture focusing technique.

Synthetic aperture focusing technique (SAFT) is effective in restoring lateral resolution of ultrasonic images for scans with focusing-related distortions. Although successfully applied in pulse-echo ultrasonics, the physical nature of an optoacoustic modality requires a modified algorithm to return accurate results. The SIR-SAFT method reported here uses the spatial impulse response (SIR) of the transducer to weight the contributions to the SAFT and is tailored to provide significant resolution and signal gains for out-of-focus sources in scanning optoacoustic microscopy systems. Furthermore, the SIR-SAFT is implemented in full three dimensions, applicable to signals both far of and at the focus of the ultrasonic detector. The method has been further shown to outperform conventional SAFT algorithms for both simulated and experimental optoacoustic data.

Thermally modulated photoacoustic imaging with super-paramagnetic iron oxide nanoparticles.

Thermally modulated photoacoustic imaging (TMPI) is reported here for contrast enhancement when using nanoparticles as contrast agents. Exploiting the excellent sensitivity of the photoacoustic (PA) process on temperature and the highly selective heating capability of nanoparticles under electromagnetic field, the PA signals stemming from the nanoparticles labeled region can be efficiently modulated whereas those from highly light absorptive backgrounds are minimally affected. A coherent difference imaging procedure reduces the background signal and thus improves the imaging contrast. Phantom experiments with super-paramagnetic iron oxide nanoparticles (SPIONs) as contrast agents and alternating magnetic fields for heating are demonstrated. Further improvements toward clinical applications are also discussed.

Continuous scanning of a time-reversed ultrasonically encoded optical focus by reflection-mode digital phase conjugation.

Time-reversed ultrasonically encoded (TRUE) optical focusing in turbid media was previously implemented using both analog and digital phase conjugation. The digital approach, in addition to its large energy gain, can improve the focal intensity and resolution by iterative focusing. However, performing iterative focusing at each focal position can be time-consuming. Here, we show that by gradually moving the focal position, the TRUE focal intensity is improved, as in iterative focusing at a fixed position, and can be continuously scanned to image fluorescent targets in a shorter time. In addition, our setup is, to the best of our knowledge, the first demonstration of TRUE focusing using a digital phase conjugate mirror in a reflection mode, which is more suitable for practical applications.

In vivo optically encoded photoacoustic flowgraphy.

We present an optically encoded photoacoustic (PA) flow imaging method based on optical-resolution PA microscopy. An intensity-modulated continuous-wave laser photothermally encodes the flowing medium, and a pulsed laser generates PA waves to image the encoded heat pattern. Flow speeds can be calculated by cross correlation. The method was validated in phantoms at flow speeds ranging from 0.23 to 11 mm/s. Venous blood flow speed in a mouse ear was also measured.

Comparison of photoacoustic image reconstruction algorithms using the channelized Hotelling observer.

We demonstrate the use of task-based image-quality metrics to compare various photoacoustic image-reconstruction algorithms, including a method based on the pseudoinverse of the system matrix, simple backprojection, filtered backprojection, and a method based on the Fourier transform. We use a three-dimensional forward model with a linear transducer array to simulate a photoacoustic imaging system. The reconstructed images correspond with two-dimensional slices of the object and are 128×128 pixels. In order to compare the algorithms, we use channelized Hotelling observers that predict the detection ability of human observers. We use two sets of channels: constant Q and difference of Gaussian spatial frequency channels. We look at three tasks, identification of a point source in a uniform background, identification of a 0.5-mm cube in a uniform background, and identification of a point source in a lumpy background. For the lumpy background task, which is the most realistic of the tasks, the method based on the pseudoinverse performs best according to both sets of channels.

Vessel orientation-dependent sensitivity of optoacoustic imaging using a linear array transducer.

For clinical optoacoustic imaging, linear probes are preferably used because they allow versatile imaging of the human body with real-time display and free-hand probe guidance. The two-dimensional (2-D) optoacoustic image obtained with this type of probe is generally interpreted as a 2-D cross-section of the tissue just as is common in echo ultrasound. We demonstrate in three-dimensional simulations, phantom experiments, and in vivo mouse experiments that for vascular imaging this interpretation is often inaccurate. The cylindrical blood vessels emit anisotropic acoustic transients, which can be sensitively detected only if the direction of acoustic radiation coincides with the probe aperture. Our results reveal for this reason that the signal amplitude of different blood vessels may differ even if the vessels have the same diameter and initial pressure distribution but different orientation relative to the imaging plane. This has important implications for the image interpretation, for the probe guidance technique, and especially in cases when a quantitative reconstruction of the optical tissue properties is required.

Parallel acoustic delay lines for photoacoustic tomography.

Achieving real-time photoacoustic (PA) tomography typically requires multi-element ultrasound transducer arrays and their associated multiple data acquisition (DAQ) electronics to receive PA waves simultaneously. We report the first demonstration of a photoacoustic tomography (PAT) system using optical fiber-based parallel acoustic delay lines (PADLs). By employing PADLs to introduce specific time delays, the PA signals (on the order of a few micro seconds) can be forced to arrive at the ultrasonic transducers at different times. As a result, time-delayed PA signals in multiple channels can be ultimately received and processed in a serial manner with a single-element transducer, followed by single-channel DAQ electronics. Our results show that an optically absorbing target in an optically scattering medium can be photoacoustically imaged using the newly developed PADL-based PAT system. Potentially, this approach could be adopted to significantly reduce the complexity and cost of ultrasonic array receiver systems.

Artefact reduction in optoacoustic tomographic imaging by estimating the distribution of acoustic scatterers.

The quality of optoacoustic tomographic reconstructions can be severely affected by acoustic reflections or scattering arising at interfaces of highly mismatched organs, such as bones, lungs, or other air-containing cavities. We present a procedure to reduce the associated artefacts based on estimation of the acoustic scatterers distribution within the imaged object. Signals generated by a strong optical absorber are processed and used in a weighted back-projection algorithm. Experimental results in a tissue-mimicking phantom clearly demonstrate improved performance as compared to the case in which no information on the distribution of acoustic scatterers is available.

Frequency-domain photoacoustic phased array probe for biomedical imaging applications.

We report the development of a frequency-domain biomedical photoacoustic imaging system that utilizes a continuous-wave laser source with a custom intensity modulation pattern, ultrasonic phased array for signal detection, and processing coupled with a beam-forming algorithm for reconstruction of photoacoustic correlation images. Sensitivity to optical contrast was demonstrated using tissue-mimicking phantoms and in-vivo tissue samples.