Four-dimensional photoacoustic imaging of moving targets

Photoacoustic Imaging in Biomedicine and Life Sciences

Photo-acoustic imaging, also known as opto-acoustic imaging, has become a widely popular modality for biomedical applications. This hybrid technique possesses the advantages of high optical contrast and high ultrasonic resolution. Due to the distinct optical absorption properties of tissue compartments and main chromophores, photo-acoustics is able to non-invasively observe structural and functional variations within biological tissues including oxygenation and deoxygenation, blood vessels and spatial melanin distribution. The detection of acoustic waves produced by a pulsed laser source yields a high scaling range, from organ level photo-acoustic tomography to sub-cellular or even molecular imaging. This review discusses significant novel technical solutions utilising photo-acoustics and their applications in the fields of biomedicine and life sciences.

Development of a Stationary 3D Photoacoustic Imaging System Using Sparse Single-Element Transducers: Phantom Study

Photoacoustic imaging (PAI) is an emerging label-free and non-invasive modality for imaging biological tissues. PAI has been implemented in different configurations, one of which is photoacoustic computed tomography (PACT) with a potential wide range of applications, including brain and breast imaging. Hemispherical Array PACT (HA-PACT) is a variation of PACT that has solved the limited detection-view problem. Here, we designed an HA-PACT system consisting of 50 single element transducers. For implementation, we initially performed a simulation study, with parameters close to those in practice, to determine the relationship between the number of transducers and the quality of the reconstructed image. We then used the greatest number of transducers possible on the hemisphere and imaged copper wire phantoms coated with a light absorbing material to evaluate the performance of the system. Several practical issues such as light illumination, arrangement of the transducers, and an image reconstruction algorithm have been comprehensively studied.

Analysis of Negatively Focused Ultrasound Detectors in Optoacoustic Tomography

In optoacoustic tomography, negatively focused transducers may be used for improving the tangential image resolution while preserving a high signal-to-noise ratio. Commonly, image reconstruction in such scenarios is facilitated by the use of the virtual-detector approach. Although the validity of this approach has been experimentally verified, it is based on an approximation whose effect on optoacoustic image reconstruction has not yet been studied. In this paper, we analyze the response of negatively focused acoustic detectors in 2D in both time and frequency domains. Based on this analysis, tradeoffs between the detector size, curvature, and sensitivity are formulated. In addition, our analysis reveals the geometrical underpinning of the virtual-detector approximation and quantifies its deviation from the exact solution. The error involved in the virtual-detector approximation is studied in image reconstruction simulations and its effect on image quality is shown. The theoretical tools developed in this work may be used in the design of new optoacoustic detection geometries as well as for improved image reconstruction.

Objective assessment and design improvement of a staring, sparse transducer array by the spatial crosstalk matrix for 3D photoacoustic tomography

Accurate reconstruction of 3D photoacoustic (PA) images requires detection of photoacoustic signals from many angles. Several groups have adopted staring ultrasound arrays, but assessment of array performance has been limited. We previously reported on a method to calibrate a 3D PA tomography (PAT) staring array system and analyze system performance using singular value decomposition (SVD). The developed SVD metric, however, was impractical for large system matrices, which are typical of 3D PAT problems. The present study consisted of two main objectives. The first objective aimed to introduce the crosstalk matrix concept to the field of PAT for system design. Figures-of-merit utilized in this study were root mean square error, peak signal-to-noise ratio, mean absolute error, and a three dimensional structural similarity index, which were derived between the normalized spatial crosstalk matrix and the identity matrix. The applicability of this approach for 3D PAT was validated by observing the response of the figures-of-merit in relation to well-understood PAT sampling characteristics (i.e. spatial and temporal sampling rate). The second objective aimed to utilize the figures-of-merit to characterize and improve the performance of a near-spherical staring array design. Transducer arrangement, array radius, and array angular coverage were the design parameters examined. We observed that the performance of a 129-element staring transducer array for 3D PAT could be improved by selection of optimal values of the design parameters. The results suggested that this formulation could be used to objectively characterize 3D PAT system performance and would enable the development of efficient strategies for system design optimization.

Acoustic Inversion in Optoacoustic Tomography: A Review

Optoacoustic tomography enables volumetric imaging with optical contrast in biological tissue at depths beyond the optical mean free path by the use of optical excitation and acoustic detection. The hybrid nature of optoacoustic tomography gives rise to two distinct inverse problems: The optical inverse problem, related to the propagation of the excitation light in tissue, and the acoustic inverse problem, which deals with the propagation and detection of the generated acoustic waves. Since the two inverse problems have different physical underpinnings and are governed by different types of equations, they are often treated independently as unrelated problems. From an imaging standpoint, the acoustic inverse problem relates to forming an image from the measured acoustic data, whereas the optical inverse problem relates to quantifying the formed image. This review focuses on the acoustic aspects of optoacoustic tomography, specifically acoustic reconstruction algorithms and imaging-system practicalities. As these two aspects are intimately linked, and no silver bullet exists in the path towards high-performance imaging, we adopt a holistic approach in our review and discuss the many links between the two aspects. Four classes of reconstruction algorithms are reviewed: time-domain (so called back-projection) formulae, frequency-domain formulae, time-reversal algorithms, and model-based algorithms. These algorithms are discussed in the context of the various acoustic detectors and detection surfaces which are commonly used in experimental studies. We further discuss the effects of non-ideal imaging scenarios on the quality of reconstruction and review methods that can mitigate these effects. Namely, we consider the cases of finite detector aperture, limited-view tomography, spatial under-sampling of the acoustic signals, and acoustic heterogeneities and losses.

4-D photoacoustic tomography

Photoacoustic tomography (PAT) offers three-dimensional (3D) structural and functional imaging of living biological tissue with label-free, optical absorption contrast. These attributes lend PAT imaging to a wide variety of applications in clinical medicine and preclinical research. Despite advances in live animal imaging with PAT, there is still a need for 3D imaging at centimeter depths in real-time. We report the development of four dimensional (4D) PAT, which integrates time resolutions with 3D spatial resolution, obtained using spherical arrays of ultrasonic detectors. The 4D PAT technique generates motion pictures of imaged tissue, enabling real time tracking of dynamic physiological and pathological processes at hundred micrometer-millisecond resolutions. The 4D PAT technique is used here to image needle-based drug delivery and pharmacokinetics. We also use this technique to monitor 1) fast hemodynamic changes during inter-ictal epileptic seizures and 2) temperature variations during tumor thermal therapy.

Three-dimensional optoacoustic tomography at video rate

Using optoacoustic excitation, a complete volumetric tomographic data sets from the imaged object can in principle be generated with a single interrogating laser pulse. Thus, optoacoustic imaging intrinsically has the potential for fast three-dimensional imaging. We have developed a system capable of acquiring volumetric optoacoustic data in real time and showcase in this work the undocumented capacity to generate high resolution three-dimensional optoacoustic images at a rate of 10 Hz, currently mainly limited by the pulse repetition rate of the excitation laser.

Model-based optoacoustic inversion with arbitrary-shape detectors

PURPOSE

Optoacoustic imaging enables mapping the optical absorption of biological tissue using optical excitation and acoustic detection. Although most image-reconstruction algorithms are based on the assumption of a detector with an isotropic sensitivity, the geometry of the detector often leads to a response with spatially dependent magnitude and bandwidth. This effect may lead to attenuation or distortion in the recorded signal and, consequently, in the reconstructed image.

METHODS

Herein, an accurate numerical method for simulating the spatially dependent response of an arbitrary-shape acoustic transducer is presented. The method is based on an analytical solution obtained for a two-dimensional line detector. The calculated response is incorporated in the forward model matrix of an optoacoustic imaging setup using temporal convolution, and image reconstruction is performed by inverting the matrix relation.

RESULTS

The method was numerically and experimentally demonstrated in two dimensions for both flat and focused transducers and compared to the spatial-convolution method. In forward simulations, the developed method did not suffer from the numerical errors exhibited by the spatial-convolution method. In reconstruction simulations and experiments, the use of both temporal-convolution and spatial-convolution methods lead to an enhancement in resolution compared to a reconstruction with a point detector model. However, because of its higher modeling accuracy, the temporal-convolution method achieved a noise figure approximated three times lower than the spatial-convolution method.

CONCLUSIONS

The demonstrated performance of the spatial-convolution method shows it is a powerful tool for reducing reconstruction artifacts originating from the detector finite size and improving the quality of optoacoustic reconstructions. Furthermore, the method may be used for assessing new system designs. Specifically, detectors with nonstandard shapes may be investigated.

Singular value decomposition analysis of a photoacoustic imaging system and 3D imaging at 0.7 FPS

Photoacoustic imaging is a non-ionizing imaging modality that provides contrast consistent with optical imaging techniques while the resolution and penetration depth is similar to ultrasound techniques. In a previous publication [Opt. Express 18, 11406 (2010)], a technique was introduced to experimentally acquire the imaging operator for a photoacoustic imaging system. While this was an important foundation for future work, we have recently improved the experimental procedure allowing for a more densely populated imaging operator to be acquired. Subsets of the imaging operator were produced by varying the transducer count as well as the measurement space temporal sampling rate. Examination of the matrix rank and the effect of contributing object space singular vectors to image reconstruction were performed. For a PAI system collecting only limited data projections, matrix rank increased linearly with transducer count and measurement space temporal sampling rate. Image reconstruction using a regularized pseudoinverse of the imaging operator was performed on photoacoustic signals from a point source, line source, and an array of point sources derived from the imaging operator. As expected, image quality increased for each object with increasing transducer count and measurement space temporal sampling rate. Using the same approach, but on experimentally sampled photoacoustic signals from a moving point-like source, acquisition, data transfer, reconstruction and image display took 1.4 s using one laser pulse per 3D frame. With relatively simple hardware improvements to data transfer and computation speed, our current imaging results imply that acquisition and display of 3D photoacoustic images at laser repetition rates of 10Hz is easily achieved.

Ultrafast photoacoustic imaging and its application to real-time 3D imaging with improved focusing

The restricted temporal resolution ofphotoacoustic imaging due to limited frame rates often prohibits its applications in areas such as real-time 3D imaging. This paper presents an ultrasound/photoacoustic multimodality imaging system that provides an ultrafast frame rate and consists of an ultrasound transducer array with plane wave excitation and a laser with pulse repetition frequency up to 2000 Hz. Its application to real-time 3D photoacoustic imaging is demonstrated and a synthetic-aperture focusing technique is applied to improve the elevational focusing quality of the mechanically-scanned 1D array. A 3D frame rate of 12 Hz in a volume covering a 19.2 mm x 19.2 mm scanning surface is demonstrated.

Localization of spherical lesions in tumor-mimicking phantoms by 3D sparse array photoacoustic imaging

PURPOSE

The authors have developed a sparse-array photoacoustic imaging (SPAI) system that is capable of mapping 3D distributions of optical absorption using a small number of laser pulses with no mechanical scanning needed. In previous studies, the authors have shown the localization accuracy and the high frame-rate image acquisition on simple phantoms with limited medical relevance. The purpose of this study was to test the imaging capabilities of SPAI in the context of breast tumor detection and localization.

METHODS

The authors constructed an array of phantoms that include spherical lesions of sizes 1.5-9 mm, buried in highly scattering tissue phantoms at depths of 3-30 mm. The authors investigated both homogeneous lesions made of blood at varying concentrations and heterogeneous lesions containing vessel-like structures. Volumetric images of the deeply buried lesions were taken at increasingly shallower depths and image-based localization was compared to measured depth.

RESULTS

The authors were able to detect and accurately localize homogeneous lesions having a realistic absorption coefficient of 0.2 cm(-1) down to depths of 9-20 mm, and heterogeneous lesions containing 0.5 mm diameter vessel-like structures down to depths of 13-20 mm. Image acquisition required 2.5 s for each volumetric lesion image.

CONCLUSIONS

These results suggest that 3D SPAI can detect highly vascularized lesions well below 1 cm in diameter and can overcome optical scatter of tissue to depths of 1-2 cm. With further improvement in the sensitivity and noise characteristics of the imaging system, similar imaging depths should be within reach in real breast tissue. The method, due to its optical contrast, 3D imaging, and fast acquisition, may prove useful in the clinic as an adjunct to existing breast screening tools.

Sentinel lymph nodes and lymphatic vessels: noninvasive dual-modality in vivo mapping by using indocyanine green in rats--volumetric spectroscopic photoacoustic imaging and planar fluorescence imaging

PURPOSE

To noninvasively map sentinel lymph nodes (SLNs) and lymphatic vessels in rats in vivo by using dual-modality nonionizing imaging-volumetric spectroscopic photoacoustic imaging, which measures optical absorption, and planar fluorescence imaging, which measures fluorescent emission-of indocyanine green (ICG).

MATERIALS AND METHODS

Institutional animal care and use committee approval was obtained. Healthy Sprague-Dawley rats weighing 250-420 g (age range, 60-120 days) were imaged by using volumetric photoacoustic imaging (n = 5) and planar fluorescence imaging (n = 3) before and after injection of 1 mmol/L ICG. Student paired t tests based on a logarithmic scale were performed to evaluate the change in photoacoustic signal enhancement of SLNs and lymphatic vessels before and after ICG injection. The spatial resolutions of both imaging systems were compared at various imaging depths (2-8 mm) by layering additional biologic tissues on top of the rats in vivo. Spectroscopic photoacoustic imaging was applied to identify ICG-dyed SLNs.

RESULTS

In all five rats examined with photoacoustic imaging, SLNs were clearly visible, with a mean signal enhancement of 5.9 arbitrary units (AU) + or - 1.8 (standard error of the mean) (P < .002) at 0.2 hour after injection, while lymphatic vessels were seen in four of the five rats, with a signal enhancement of 4.3 AU + or - 0.6 (P = .001). In all three rats examined with fluorescence imaging, SLNs and lymphatic vessels were seen. The average full width at half maximum (FWHM) of the SLNs in the photoacoustic images at three imaging depths (2, 6, and 8 mm) was 2.0 mm + or - 0.2 (standard deviation), comparable to the size of a dissected lymph node as measured with a caliper. However, the FWHM of the SLNs in fluorescence images widened from 8 to 22 mm as the imaging depth increased, owing to strong light scattering. SLNs were identified spectroscopically in photoacoustic images.

CONCLUSION

These two modalities, when used together with ICG, have the potential to help map SLNs in axillary staging and to help evaluate tumor metastasis in patients with breast cancer.

Analysis of a photoacoustic imaging system by the crosstalk matrix and singular value decomposition

Photoacoustic imaging is a hybrid imaging modality capable of producing contrast similar to optical imaging techniques but with increased penetration depth and resolution in turbid media by encoding the information as acoustic waves. In general, it is important to characterize the performance of a photoacoustic imaging system by parameters such as sensitivity, resolution, and contrast. However, system characterization can extend beyond these metrics by implementing advanced analysis via the crosstalk matrix and singular value decomposition. A method was developed to experimentally measure a matrix that represented the imaging operator for a photoacoustic imaging system. Computations to produce the crosstalk matrix were completed to provide insight into the spatially dependent sensitivity and aliasing for the photoacoustic imaging system. Further analysis of the imaging operator was done via singular value decomposition to estimate the capability of the imaging system to reconstruct objects and the inherent sensitivity to those objects. The results provided by singular value decomposition were compared to SVD results from a de-noised imaging operator to estimate the number of measurable singular vectors for the system. These characterization techniques can be broadly applied to any photoacoustic system and, with regards to the studied system, could be used as a basis for improvements to future iterations.

Photoacoustic imaging velocimetry for flow-field measurement

We present the photoacoustic imaging velocimetry (PAIV) method for flow-field measurement based on a linear transducer array. The PAIV method is realized by using a Q-switched pulsed laser, a linear transducer array, a parallel data-acquisition equipment and dynamic focusing reconstruction. Tracers used to track liquid flow field were real-timely detected, two-dimensional (2-D) flow visualization was successfully reached, and flow parameters were acquired by measuring the movement of the tracer. Experimental results revealed that the PAIV method would be developed into 3-D imaging velocimetry for flow-field measurement, and potentially applied to research the security and targeting efficiency of optical nano-material probes.

Compressed sensing in photoacoustic tomography in vivo

The data acquisition speed in photoacoustic computed tomography (PACT) is limited by the laser repetition rate and the number of parallel ultrasound detecting channels. Reconstructing an image with fewer measurements can effectively accelerate the data acquisition and reduce the system cost. We adapt compressed sensing (CS) for the reconstruction in PACT. CS-based PACT is implemented as a nonlinear conjugate gradient descent algorithm and tested with both phantom and in vivo experiments.

Reconstruction of high quality photoacoustic tomography with a limited-view scanning

The goal of this work is to resolve the limited-view problem of photoacoustic tomography (PAT). We report a two-loop iteration method to inverse the photoacoustic sources from the measured photoacoustic signals. PAT reconstruction with this method does not depend on the detection path. Therefore, the proposed method can provide recognizable image even when the detector only scans a small angle (about 20 degrees approximately 30 degrees). The comparison with the delay-and-sum method shows the advantage of the proposed method in reconstructing image from incomplete data.

Development and characterization of an omnidirectional photoacoustic point source for calibration of a staring 3D photoacoustic imaging system

Photoacoustic imaging is a modality which makes use of the contrast provided by optical imaging techniques and the spatial resolution and penetration depth similar to acoustic imaging modalities. We have developed a method for fast 3D photoacoustic imaging using a sparse hemispherical array of transducers. Such a system requires characterization of the transducer's response to an ideal point source in order to accurately reconstruct objects in the imaging volume. First, an attempt was made to design an ideal photoacoustic point source via a combination of liquids which would appropriately scatter and absorb the light such that a spherical distribution was achieved. Methylene blue (MB(+)) was used as the primary optical absorber while Intralipid (IL) was used as the liquid responsible for the optical scatter. A multitude of combinations were tested and the signal uniformity was characterized. The combination of 200 microM MB(+) and 0.09% IL was found to produce the most uniform signal over the range of transducers in the hemispherical array. The liquid source was then characterized over a broader range of azimuthal and zenith angles where it was shown the azimuthal consistency was much greater than the stability seen in different zenith elevations. The source was then used in a calibration scan for an imaging volume of 40 x 40 x 40 mm(3). At 216 points evenly spaced in the imaging volume, parameters were recorded for signal amplitude, width, and time-of-flight. These calibration parameters could then be applied to an iterative reconstruction algorithm in an attempt to more accurately produce images.