Iterative image reconstruction in transcranial photoacoustic tomography based on the elastic wave equation.

Photoacoustic computed tomography (PACT) is an emerging computed imaging modality that exploits optical contrast and ultrasonic detection principles to form images of the photoacoustically induced initial pressure distribution within tissue. The PACT reconstruction problem corresponds to a time-domain inverse source problem, where the initial pressure distribution is recovered from the measurements recorded on an aperture outside the support of the source. A major challenge in transcranial PACT of the brain is to compensate for aberrations and attenuation in the measured data due to the propagation of the photoacoustic wavefields through the skull. To properly account for these effects, a wave equation-based inversion method can be employed that can model the heterogeneous elastic properties of the medium. In this study, an optimization-based image reconstruction method for 3D transcranial PACT is developed based on the elastic wave equation. To accomplish this, a forward-adjoint operator pair based on a finite-difference time-domain discretization of the 3D elastic wave equation is utilized to compute penalized least squares estimates of the initial pressure distribution. Computer-simulation and experimental studies are conducted to investigate the robustness of the reconstruction method to model mismatch and its ability to effectively resolve cortical and superficial brain structures.

A survey of computational frameworks for solving the acoustic inverse problem in three-dimensional photoacoustic computed tomography.

Photoacoustic computed tomography (PACT), also known as optoacoustic tomography, is an emerging imaging technique that holds great promise for biomedical imaging. PACT is a hybrid imaging method that can exploit the strong endogenous contrast of optical methods along with the high spatial resolution of ultrasound methods. In its canonical form that is addressed in this article, PACT seeks to estimate the photoacoustically-induced initial pressure distribution within the object. Image reconstruction methods are employed to solve the acoustic inverse problem associated with the image formation process. When an idealized imaging scenario is considered, analytic solutions to the PACT inverse problem are available; however, in practice, numerous challenges exist that are more readily addressed within an optimization-based, or iterative, image reconstruction framework. In this article, the PACT image reconstruction problem is reviewed within the context of modern optimization-based image reconstruction methodologies. Imaging models that relate the measured photoacoustic wavefields to the sought-after object function are described in their continuous and discrete forms. The basic principles of optimization-based image reconstruction from discrete PACT measurement data are presented, which includes a review of methods for modeling the PACT measurement system response and other important physical factors. Non-conventional formulations of the PACT image reconstruction problem, in which acoustic parameters of the medium are concurrently estimated along with the PACT image, are also introduced and reviewed.

All-reflective ring illumination system for photoacoustic tomography.

Given that breast cancer is the second leading cause of cancer-related deaths among women in the United States, it is necessary to continue improving the sensitivity and specificity of breast imaging systems that diagnose breast lesions. Photoacoustic (PA) imaging can provide functional information during in vivo studies and can augment the structural information provided by ultrasound (US) imaging. A full-ring, all-reflective, illumination system for photoacoustic tomography (PAT) coupled to a full-ring US receiver is developed and tested. The US/PA tomography system utilizes a cone mirror and conical reflectors to optimize light delivery for PAT imaging and has the potential to image objects that are placed within the ring US transducer. The conical reflector used in this system distributes the laser energy over a circular cross-sectional area, thereby reducing the overall fluence. This, in turn, allows the operator to increase the laser energy achieving better cross-sectional penetration depth. A proof-of-concept design utilizing a single cone mirror and a parabolic reflector is used for imaging cylindrical phantoms with light-absorbing objects. For the given phantoms, it has been shown that there was no restriction in imaging a given targeted cross-sectional area irrespective of vertical depth, demonstrating the potential of mirror-based, ring-illuminated PAT system. In addition, the all-reflective ring illumination method shows a uniform PA signal across the scanned cross-sectional area.

Parameterized joint reconstruction of the initial pressure and sound speed distributions for photoacoustic computed tomography.

Accurate estimation of the initial pressure distribution in photoacoustic computed tomography (PACT) depends on knowledge of the sound speed distribution. However, the sound speed distribution is typically unknown. Further, the initial pressure and sound speed distributions cannot both, in general, be stably recovered from PACT measurements alone. In this work, a joint reconstruction (JR) method for the initial pressure distribution and a low-dimensional parameterized model of the sound speed distribution is proposed. By employing a priori information about the structure of the sound speed distribution, both the initial pressure and sound speed can be accurately recovered. The JR problem is solved by use of a proximal optimization method that allows constraints and non-smooth regularization functions for the initial pressure distribution. The gradients of the cost function with respect to the initial pressure and sound speed distributions are calculated by use of an adjoint state method that has the same per-iteration computational cost as calculating the gradient with respect to the initial pressure distribution alone. This approach is evaluated through 2D computer-simulation studies for a small animal imaging model and by application to experimental in vivo measurements of a mouse.

Mitigation of artifacts due to isolated acoustic heterogeneities in photoacoustic computed tomography using a variable data truncation-based reconstruction method.

Photoacoustic computed tomography (PACT) is an emerging computed imaging modality that exploits optical contrast and ultrasonic detection principles to form images of the absorbed optical energy density within tissue. If the object possesses spatially variant acoustic properties that are unaccounted for by the reconstruction method, the estimated image can contain distortions. While reconstruction methods have recently been developed to compensate for this effect, they generally require the object's acoustic properties to be known a priori. To circumvent the need for detailed information regarding an object's acoustic properties, we previously proposed a half-time reconstruction method for PACT. A half-time reconstruction method estimates the PACT image from a data set that has been temporally truncated to exclude the data components that have been strongly aberrated. However, this method can be improved upon when the approximate sizes and locations of isolated heterogeneous structures, such as bones or gas pockets, are known. To address this, we investigate PACT reconstruction methods that are based on a variable data truncation (VDT) approach. The VDT approach represents a generalization of the half-time approach, in which the degree of temporal truncation for each measurement is determined by the distance between the corresponding ultrasonic transducer location and the nearest known bone or gas void location. Computer-simulated and experimental data are employed to demonstrate the effectiveness of the approach in mitigating artifacts due to acoustic heterogeneities.

A forward-adjoint operator pair based on the elastic wave equation for use in transcranial photoacoustic computed tomography.

Photoacoustic computed tomography (PACT) is an emerging imaging modality that exploits optical contrast and ultrasonic detection principles to form images of the photoacoustically induced initial pressure distribution within tissue. The PACT reconstruction problem corresponds to an inverse source problem in which the initial pressure distribution is recovered from measurements of the radiated wavefield. A major challenge in transcranial PACT brain imaging is compensation for aberrations in the measured data due to the presence of the skull. Ultrasonic waves undergo absorption, scattering and longitudinal-to-shear wave mode conversion as they propagate through the skull. To properly account for these effects, a wave-equation-based inversion method should be employed that can model the heterogeneous elastic properties of the skull. In this work, a forward model based on a finite-difference time-domain discretization of the three-dimensional elastic wave equation is established and a procedure for computing the corresponding adjoint of the forward operator is presented. Massively parallel implementations of these operators employing multiple graphics processing units (GPUs) are also developed. The developed numerical framework is validated and investigated in computer19 simulation and experimental phantom studies whose designs are motivated by transcranial PACT applications.

In vivo photoacoustic flowmetry at depths of the diffusive regime based on saline injection.

We propose a saline injection-based method to quantify blood flow velocity in vivo with acoustic-resolution photoacoustic tomography. By monitoring the saline­blood interface propagating in the blood vessel, the flow velocity can be resolved. We first demonstrated our method in phantom experiments, where a root mean square error of prediction of 0.29 mm/s was achieved. By injecting saline into a mouse tail vein covered with 1 mm chicken tissue, we showed that the flow velocity in the tail vein could be measured at depths, which is especially pertinent to monitoring blood flow velocity in patients undergoing intravenous infusion.