Photoacoustic ultrasound (PAUS)--reconstruction tomography

Biomedical photoacoustic imaging

Photoacoustic (PA) imaging, also called optoacoustic imaging, is a new biomedical imaging modality based on the use of laser-generated ultrasound that has emerged over the last decade. It is a hybrid modality, combining the high-contrast and spectroscopic-based specificity of optical imaging with the high spatial resolution of ultrasound imaging. In essence, a PA image can be regarded as an ultrasound image in which the contrast depends not on the mechanical and elastic properties of the tissue, but its optical properties, specifically optical absorption. As a consequence, it offers greater specificity than conventional ultrasound imaging with the ability to detect haemoglobin, lipids, water and other light-absorbing chomophores, but with greater penetration depth than purely optical imaging modalities that rely on ballistic photons. As well as visualizing anatomical structures such as the microvasculature, it can also provide functional information in the form of blood oxygenation, blood flow and temperature. All of this can be achieved over a wide range of length scales from micrometres to centimetres with scalable spatial resolution. These attributes lend PA imaging to a wide variety of applications in clinical medicine, preclinical research and basic biology for studying cancer, cardiovascular disease, abnormalities of the microcirculation and other conditions. With the emergence of a variety of truly compelling in vivo images obtained by a number of groups around the world in the last 2-3 years, the technique has come of age and the promise of PA imaging is now beginning to be realized. Recent highlights include the demonstration of whole-body small-animal imaging, the first demonstrations of molecular imaging, the introduction of new microscopy modes and the first steps towards clinical breast imaging being taken as well as a myriad of in vivo preclinical imaging studies. In this article, the underlying physical principles of the technique, its practical implementation, and a range of clinical and preclinical applications are reviewed.

Biomedical photoacoustic imaging

Photoacoustic (PA) imaging, also called optoacoustic imaging, is a new biomedical imaging modality based on the use of laser-generated ultrasound that has emerged over the last decade. It is a hybrid modality, combining the high-contrast and spectroscopic-based specificity of optical imaging with the high spatial resolution of ultrasound imaging. In essence, a PA image can be regarded as an ultrasound image in which the contrast depends not on the mechanical and elastic properties of the tissue, but its optical properties, specifically optical absorption. As a consequence, it offers greater specificity than conventional ultrasound imaging with the ability to detect haemoglobin, lipids, water and other light-absorbing chomophores, but with greater penetration depth than purely optical imaging modalities that rely on ballistic photons. As well as visualizing anatomical structures such as the microvasculature, it can also provide functional information in the form of blood oxygenation, blood flow and temperature. All of this can be achieved over a wide range of length scales from micrometres to centimetres with scalable spatial resolution. These attributes lend PA imaging to a wide variety of applications in clinical medicine, preclinical research and basic biology for studying cancer, cardiovascular disease, abnormalities of the microcirculation and other conditions. With the emergence of a variety of truly compelling in vivo images obtained by a number of groups around the world in the last 2-3 years, the technique has come of age and the promise of PA imaging is now beginning to be realized. Recent highlights include the demonstration of whole-body small-animal imaging, the first demonstrations of molecular imaging, the introduction of new microscopy modes and the first steps towards clinical breast imaging being taken as well as a myriad of in vivo preclinical imaging studies. In this article, the underlying physical principles of the technique, its practical implementation, and a range of clinical and preclinical applications are reviewed.

Volumetric real-time multispectral optoacoustic tomography of biomarkers

Multispectral optoacoustic tomography (MSOT) has recently been developed to enable visualization of optical contrast and tissue biomarkers, with resolution and speed representative of ultrasound. In the implementation described here, MSOT enables operation in real-time mode by capturing single cross-sectional images in <1 ms from living small animals (e.g., mice) and other tissues of similar dimensions. At the core of the method is illumination of the object using multiple wavelengths in order to resolve spectrally distinct biomarkers over background tissue chromophores. The system allows horizontal placement of a mouse in the imaging chamber and three-dimensional scanning of the entire body without the need to immerse the mouse in water. Here we provide a detailed description of the MSOT scanner components, system calibration, selection of image reconstruction algorithms and animal handling. Overall, the entire protocol can be completed within 15-30 min for acquisition of a whole-body multispectral data set from a living mouse.

Image reconstruction in photoacoustic tomography involving layered acoustic media

Photoacoustic tomography (PAT), also known as thermoacoustic or optoacoustic tomography, is a rapidly emerging biomedical imaging technique that combines optical image contrast with ultrasound detection principles. Most existing reconstruction algorithms for PAT assume the object of interest possesses homogeneous acoustic properties. The images produced by such algorithms can contain significant distortions and artifacts when the object's acoustic properties are spatially variant. In this work, we establish an image reconstruction formula for PAT applications in which a planar detection surface is employed and the to-be-imaged optical absorber is embedded in a known planar layered acoustic medium. The reconstruction formula is exact in a mathematical sense and accounts for multiple acoustic reflections between the layers of the medium. Computer-simulation studies are conducted to demonstrate and investigate the proposed method.

Model-based optoacoustic inversions with incomplete projection data

PURPOSE

Optoacoustic imaging is an emerging noninvasive imaging modality that can resolve optical contrast through several millimeters to centimeters of tissue with diffraction-limited resolution of ultrasound. Yet, quantified reconstruction of tissue absorption maps requires optoacoustic signals to be collected from as many locations around the object as possible. In many tomographic imaging scenarios, however, only limited-view or partial projection data are available, which has been shown to generate image artifacts and overall loss of quantification accuracy.

METHODS

In this article, the recently introduced interpolated-matrix-model optoacoustic inversion method is tested in different limited-view scenarios and compared to the standard backprojection algorithm. Both direct (TGSVD) and iterative (PLSQR) regularization methods are proposed to improve the accuracy of image reconstructions with their performance tested on simulated and experimental data.

RESULTS

While for full-view tomographic data the model-based inversion has been generally shown to attain higher reconstruction accuracy compared to backprojection algorithms, the incomplete tomographic datasets lead to ill-conditioned forward matrices and, consequently, to error-prone inversions, with strong artifacts following a distinct ripple-type pattern. The proposed regularization techniques are shown to stabilize the inversion and eliminate the artifacts.

CONCLUSIONS

Overall, it has been determined that the regularized interpolated-matrix-model-based optoacoustic inversions show higher accuracy than reconstructions with the standard backprojection algorithm. Finally, the combination of model-based inversion with PLSQR or TGSVD regularization methods can lead to an accurate reconstruction of limited-projection angle optoacoustic data and practical systems for optoacoustic imaging in many realistic cases where the full-view dataset is unavailable.

Fast, limited-data photoacoustic imaging for multiplexed systems using a frequency-domain estimation technique

PURPOSE

A new frequency-domain estimation algorithm has been developed that uses a priori information to simultaneously improve imaging quality and time resolution in photoacoustic tomography with incomplete data sets.

METHODS

The method involves application of a single-stage Wiener optimal filter to augment data sets by interpolation between measurement locations using relationships determined in a reference scan. The filter can be applied in real-time using FFT methods using either fixed or dynamic references and used with any imaging algorithm. The performance of the method is compared to a modified version of constrained backprojection algorithms using simulations and experimental investigations.

RESULTS

Simulations demonstrate the effectiveness of the approach for tracking dynamic photoacoustic activity for data sets with limited views (90 degrees) or tomographic views with a reduced number of acquisition angles at any given time (< or = 32). Experimental data of contrast uptake and washout using a 512-element curved transducer with 8:1 electronic multiplexing with the algorithm demonstrate full two-dimensional tomographic imaging with a temporal resolution better than 130 ms.

CONCLUSIONS

The estimation algorithm enables high spatial resolution, real-time imaging of dynamic physiological events or volumetric regions for photoacoustic systems employing multiplexing or scanning.

Target detection and quantification using a hybrid hand-held diffuse optical tomography and photoacoustic tomography system

We present a photoacoustic tomography-guided diffuse optical tomography approach using a hand-held probe for detection and characterization of deeply-seated targets embedded in a turbid medium. Diffuse optical tomography guided by coregistered ultrasound, MRI, and x ray has demonstrated a great clinical potential to overcome lesion location uncertainty and to improve light quantification accuracy. However, due to the different contrast mechanisms, some lesions may not be detectable by a nonoptical modality but yet have high optical contrast. Photoacoustic tomography utilizes a short-pulsed laser beam to diffusively penetrate into tissue. Upon absorption of the light by the target, photoacoustic waves are generated and used to reconstruct, at ultrasound resolution, the optical absorption distribution that reveals optical contrast. However, the robustness of optical property quantification of targets by photoacoustic tomography is complicated because of the wide range of ultrasound transducer sensitivity, the orientation and shape of the targets relative to the ultrasound array, and the uniformity of the laser beam. We show in this paper that the relative optical absorption map provided by photoacoustic tomography can potentially guide the diffuse optical tomography to accurately reconstruct target absorption maps.

Target detection and quantification using a hybrid hand-held diffuse optical tomography and photoacoustic tomography system

We present a photoacoustic tomography-guided diffuse optical tomography approach using a hand-held probe for detection and characterization of deeply-seated targets embedded in a turbid medium. Diffuse optical tomography guided by coregistered ultrasound, MRI, and x ray has demonstrated a great clinical potential to overcome lesion location uncertainty and to improve light quantification accuracy. However, due to the different contrast mechanisms, some lesions may not be detectable by a nonoptical modality but yet have high optical contrast. Photoacoustic tomography utilizes a short-pulsed laser beam to diffusively penetrate into tissue. Upon absorption of the light by the target, photoacoustic waves are generated and used to reconstruct, at ultrasound resolution, the optical absorption distribution that reveals optical contrast. However, the robustness of optical property quantification of targets by photoacoustic tomography is complicated because of the wide range of ultrasound transducer sensitivity, the orientation and shape of the targets relative to the ultrasound array, and the uniformity of the laser beam. We show in this paper that the relative optical absorption map provided by photoacoustic tomography can potentially guide the diffuse optical tomography to accurately reconstruct target absorption maps.

Photoacoustic tomography of ex vivo mouse hearts with myocardial infarction

In the present study, we evaluated the applicability of ex vivo photoacoustic imaging (PAI) on small animal organs. We used photoacoustic tomography (PAT) to visualize infarcted areas within murine hearts and compared these data to other imaging techniques [magnetic resonance imaging (MRI), micro-computed tomography] and histological slices. In order to induce ischemia, an in vivo ligation of the left anterior descending artery was performed on nine wild-type mice. After varying survival periods, the hearts were excised and fixed in formaldehyde. Samples were illuminated with nanosecond laser pulses delivered by a Nd:YAG pumped optical parametric oscillator. Ultrasound detection was achieved using a Mach-Zehnder interferometer (MZI) working as an integrating line detector. The voxel data were computed using a Fourier-domain based reconstruction algorithm, followed by inverse Radon transforms. The results clearly showed the capability of PAI to visualize myocardial infarction and to produce three-dimensional images with a spatial resolution of approximately 120 μm. Regions of affected muscle tissue in PAI corresponded well with the results of MRI and histology. Photoacoustic tomography utilizing a MZI for ultrasound detection allows for imaging of small tissue samples. Due to its high spatial resolution, good soft tissue contrast and comparatively low cost, PAT offers great potentials for imaging.

Optoacoustic imaging for clinical applications: devices and methods

INTRODUCTION

Optoacoustic (photoacoustic) imaging offers visualization of optical contrast in tissues, within several millimeters to centimeters, with resolutions that are typical of ultrasound imaging. This performance can offer a natural extension to widespread optical microscopy approaches, for applications from small animals to humans.

AREAS COVERED

An increasing number of optoacoustic approaches are considered for biomedical imaging. Implementations range from handheld and endoscopic operations to fixed scanner set-ups that can address a wide range of preclinical and clinical needs. This article illuminates aspects of the underlying principles of optoacoustic imaging operation and critically reviews the major system trends developed for clinical application. In addition to anatomical imaging, typically performed using single wavelength illumination, multispectral methods are also reviewed as they pertain to functional and molecular imaging. This article also highlights the advantages and limitations as well as the potential of this technology for clinical practice.

EXPERT OPINION

Optoacoustic imaging is an emerging and highly promising area of the imaging sciences that can offer high-resolution optical visualization deep within tissues. Therefore, it offers a promising alternative to existing optical systems developed for clinical use, which are generally limited to superficial or low-resolution imaging. These up-and-coming features offer a wider variety of optoacoustic approaches that are likely to be clinically deployed in the near future.

Nanometer-scale optical imaging of collagen fibers using gold nanoparticles

We describe 3D single particle tracking of gold nanoparticles (AuNPs) moving along collagen fibers in aqueous environment with two-photon excitation conditions. The photoacoustic effect at the collagen fiber caused by the irradiation with ultrashort, near-infrared laser pulses propels the particles adsorbed to the surface of the collagen fibers. We report the tracking of individual AuNPs in three dimensions with high spatial and temporal resolution, of few nanometers and milliseconds, respectively. Due to the emission signal caused by the interaction between the AuNPs and the weak chromophores in the collagen fiber, the trajectories of individual AuNPs reveal the fiber topography with nanometric resolution. The intensity along the trajectory shows that we are sensitive to the distribution of the weak chromophores on the fiber.

Statistical approach for optoacoustic image reconstruction in the presence of strong acoustic heterogeneities

A method is presented to reduce artefacts produced in optoacoustic tomography images due to internal reflection or scattering of the acoustic waves. It is based on weighting the tomographic contribution of each detector with the probability that a signal affected by acoustic mismatches is measured at that position. The correction method does not require a priori knowledge of the acoustic or optical properties of the imaged sample. Performance tests were made with agar phantoms that included air gaps for mimicking strong acoustic reflections as well as with an acoustically heterogeneous adult Zebrafish. The results obtained with the method proposed show a clear reduction of the artefacts with respect to the original images reconstructed with filtered back-projection algorithm. This performance is directly related to in vivo small animal imaging applications involving imaging in the presence of bones, lungs, and other highly mismatched organs.

Optoacoustic methods for frequency calibration of ultrasonic sensors

The frequency response of ultrasonic detectors is commonly calibrated by finding their sensitivity to incident plane waves at discrete frequencies. For certain applications, such as the emerging field of optoacoustic tomography, it is the response to point sources emitting broadband spectra that needs to be found instead. Although these two distinct sensitivity characteristics are interchangeable in the case of a flat detector and a point source at infinity, it is not the case for detectors with size considerably larger than the acoustic wavelength of interest or those having a focused aperture. Such geometries, which are common in optoacoustics, require direct calibration of the acoustic detector using a point source placed in the relevant position. In this paper, we report on novel cross-validating optoacoustic methods for measuring the frequency response of wideband acoustic sensors. The approach developed does not require pre-calibrated hydrophones and therefore can be readily adopted in any existing optoacoustic measurement configuration. The methods are successfully confirmed experimentally by measuring the frequency response of a common piezoelectric detector having a cylindrically focused shape.

Passive element enriched photoacoustic computed tomography (PER PACT) for simultaneous imaging of acoustic propagation properties and light absorption

We present a 'hybrid' imaging approach which can image both light absorption properties and acoustic transmission properties of an object in a two-dimensional slice using a computed tomography (CT) photoacoustic imager. The ultrasound transmission measurement method uses a strong optical absorber of small cross-section placed in the path of the light illuminating the sample. This absorber, which we call a passive element acts as a source of ultrasound. The interaction of ultrasound with the sample can be measured in transmission, using the same ultrasound detector used for photoacoustics. Such measurements are made at various angles around the sample in a CT approach. Images of the ultrasound propagation parameters, attenuation and speed of sound, can be reconstructed by inversion of a measurement model. We validate the method on specially designed phantoms and biological specimens. The obtained images are quantitative in terms of the shape, size, location, and acoustic properties of the examined heterogeneities.

Photoacoustic angiography of the breast

PURPOSE

The authors report a noninvasive technique and instrumentation for visualizing vasculature in the breast in three dimensions without using either ionizing radiation or exogenous contrast agents, such as iodine or gadolinium. Vasculature is visualized by virtue of its high hemoglobin content compared to surrounding breast parenchyma. The technique is compatible with dynamic contrast-enhanced studies.

METHODS

Photoacoustic sonic waves were stimulated in the breast with a pulsed laser operating at 800 nm and a mean exposure of 20 mJ/pulse over an area of approximately 20 cm2. These waves were subsequently detected by a hemispherical array of piezoelectric transducers, the temporal signals from which were filtered and backprojected to form three-dimensional images with nearly uniform k-space sampling.

RESULTS

Three-dimensional vascular images of a human volunteer demonstrated a clear visualization of vascular anatomy with submillimeter spatial resolution to a maximum depth of 40 mm using a 24 s image acquisition protocol. Spatial resolution was nearly isotropic and approached 250 microm over a 64 x 64 x 50 mm field of view.

CONCLUSIONS

The authors have successfully visualized submillimeter breast vasculature to a depth of 40 mm using an illumination intensity that is 32 times less than the maximum permissible exposure according to the American National Standard for Safe Use of Lasers. Clearly, the authors can achieve greater penetration depth in the breast by increasing the intensity and the cross-sectional area of the illumination beam. Given the 24 s image acquisition time without contrast agent, dynamic, contrast-enhanced, photoacoustic breast imaging using optically absorbing contrast agents is conceivable in the future.

Optoacoustic tomography with varying illumination and non-uniform detection patterns

Quantification of tissue morphology and biomarker distribution by means of optoacoustic tomography is an important and longstanding challenge, mainly caused by the complex heterogeneous structure of biological tissues as well as the lack of accurate and robust reconstruction algorithms. The recently introduced model-based inversion approaches were shown to mitigate some of reconstruction artifacts associated with the commonly used back-projection schemes, while providing an excellent platform for obtaining quantified maps of optical energy deposition in experimental configurations of various complexity. In this work, we introduce a weighted model-based approach, capable of overcoming reconstruction challenges caused by per-projection variations of object's illumination and other partial illumination effects. The universal weighting procedure is equally shown to reduce reconstruction artifacts associated with other experimental imperfections, such as non-uniform transducer sensitivity fields. Significant improvements in image fidelity and quantification are showcased both numerically and experimentally on tissue phantoms and mice.

Comparison of intensity-modulated continuous-wave lasers with a chirped modulation frequency to pulsed lasers for photoacoustic imaging applications

Using a Green's function solution to the photoacoustic wave equation, we compare intensity-modulated continuous-wave (CW) lasers with a chirped modulation frequency to pulsed lasers for photoacoustic imaging applications. Assuming the same transducer is used in both cases, we show that the axial resolution is identical and is determined by the transducer and material properties of the object. We derive a simple formula relating the signal-to-noise ratios (SNRs) of the two imaging systems that only depends on the fluence of each pulse and the time-bandwidth product of the chirp pulse. We also compare the SNR of the two systems assuming the fluence is limited by the American National Standards Institute (ANSI) laser safety guidelines for skin. We find that the SNR is about 20 dB to 30 dB larger for pulsed laser systems for reasonable values of the parameters. However, CW diode lasers have the advantage of being compact and relatively inexpensive, which may outweigh the lower SNR in many applications.

Integration of photoacoustic imaging and high-intensity focused ultrasound

We have developed an integrated photoacoustic imaging (PAI) and high-intensity focused ultrasound (HIFU) system for solid tumor treatments. A single-element, spherically focused ultrasonic transducer, with a central frequency of 5 MHz, was used to induce HIFU lesions in soft tissue. The same ultrasonic transducer was also used as a detector during PAI to guide HIFU ablation. The use of same transducer for PAI and HIFU can reduce the requirement on acoustic windows during the imaging-guided therapy, as well as ensuring the correct alignment between the therapeutic beam and the planned treatment volume. During an experiment, targeted soft tissue was first imaged by PAI. The resulted image was used to plan the subsequent HIFU ablation. After the HIFU ablation, targeted soft tissue was imaged again by PAI to evaluate the effectiveness of treatments. Good contrast was obtained between photoacoustic images before and after HIFU ablation. In conclusion, our results demonstrated that PAI technology may potentially be integrated with HIFU ablation for image-guided therapy.

Fast semi-analytical model-based acoustic inversion for quantitative optoacoustic tomography

We present a fast model-based inversion algorithm for quantitative 2-D and 3-D optoacoustic tomography. The algorithm is based on an accurate and efficient forward model, which eliminates the need for regularization in the inversion process while providing modeling flexibility essential for quantitative image formation. The resulting image-reconstruction method eliminates stability problems encountered in previously published model-based techniques and, thus, enables performing image reconstruction in real time. Our model-based framework offers a generalization of the forward solution to more comprehensive optoacoustic propagation models, such as including detector frequency response, without changing the inversion procedure. The reconstruction speed and other algorithmic performances are demonstrated using numerical simulation studies and experimentally on tissue-mimicking optically heterogeneous phantoms and small animals. In the experimental examples, the model-based reconstructions manifested correctly the effect of light attenuation through the objects and did not suffer from the artifacts which usually afflict the commonly used filtered backprojection algorithms, such as negative absorption values.

Image reconstruction in photoacoustic tomography with variable speed of sound using a higher-order geometrical acoustics approximation

Previous research correcting for variable speed of sound in photoacoustic tomography (PAT) based on a generalized radon transform (GRT) model assumes first-order geometrical acoustics (GA) approximation. In the GRT model, the pressure is related to the optical absorption, in an acoustically inhomogeneous medium, through integration over nonspherical isochronous surfaces. Previous research based on the GRT model assumes that the path taken by acoustic rays is linear and neglects amplitude perturbations to the measured pressure. We have derived a higher-order GA expression that takes into account the first-order effect in the amplitude of the measured signal and higher-order perturbation to the travel times. The higher-order perturbation to travel time incorporates the effect of ray bending. Incorrect travel times can lead to image distortion and blurring. These corrections are expected to impact image quality and quantitative PAT. We have previously shown that travel-time corrections in 2-D suggest that perceivable differences in the isochronous surfaces can be seen when the second-order travel-time perturbations are taken into account with a 10% speed-of-sound variation. In this work, we develop iterative image reconstruction algorithms that incorporate this higher-order GA approximation assuming that the speed of sound map is known. We evaluate the effect of higher-order GA approximation on image quality and accuracy.

Silver nanosystems for photoacoustic imaging and image-guided therapy

Due to their optical absorption properties, metallic nanoparticles are excellent photoacoustic imaging contrast agents. A silver nanosystem is presented here as a potential contrast agent for photoacoustic imaging and image-guided therapy. Currently, the nanosystem consists of a porous silver layer deposited on the surface of spherical silica cores ranging in diameter from 180 to 520 nm. The porous nature of the silver layer will allow for release of drugs or other therapeutic agents encapsulated in the core in future applications. In their current PEGylated form, the silver nanosystem is shown to be nontoxic in vitro at concentrations of silver up to 2 mgml. Furthermore, the near-infrared absorbance properties of the nanosystem are demonstrated by measuring strong, concentration-dependent photoacoustic signal from the silver nanosystem embedded in an ex vivo tissue sample. Our study suggests that silver nanosystems can be used as multifunctional agents capable of augmenting image-guided therapy techniques.

Real-time optoacoustic monitoring and three-dimensional mapping of a human arm vasculature

We present our findings from a real-time laser optoacoustic imaging system (LOIS). The system utilizes a Q-switched Nd:YAG laser; a standard 128-channel ultrasonic linear array probe; custom electronics and custom software to collect, process, and display optoacoustic (OA) images at 10 Hz. We propose that this system be used during preoperative mapping of forearm vessels for hemodialysis treatment. To demonstrate the real-time imaging capabilities of the system, we show OA images of forearm vessels in a volunteer and compare our results to ultrasound images of the same region. Our OA images show blood vessels in high contrast. Manipulations with the probe enable us to locate and track arteries and veins of a forearm in real time. We also demonstrate the ability to combine a series of OA image slices into a volume for spatial representation of the vascular network. Finally, we use frame-by-frame analysis of the recorded OA video to measure dynamic changes of the crossection of the ulnar artery.

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.

Fast and robust deconvolution-based image reconstruction for photoacoustic tomography in circular geometry: experimental validation

Photoacoustic tomography (PAT) is a fast-developing biomedical imaging technology suitable for in vivo imaging. PAT in spherical or circular geometry gives good image resolution yet is slow or expensive in signal acquisition and image formation. Reducing the number of detection angles can ameliorate such issues usually at the expense of image quality. This paper introduces a deconvolution-based algorithm that models the imaging process as a linear and shift-invariant system. As demonstrated by the in vivo experiment, this algorithm not only runs much faster than the back-projection algorithm, but also shows stronger robustness in that it provides better image quality when detection angles are sparse. Therefore, this algorithm promises to enable real-time PAT in circular geometry.

Photoacoustic tomography for imaging nanoparticles

Nanotechnology is the key to a new, noninvasive photoacoustic imaging technique that could detect early stages of disease tissues. The combination of photoacoustic imaging with nanotechnology holds promise for determining the structural and functional properties of tissues with enhanced sensitivity and specificity and for monitoring the treatment of diseases. In this chapter, we described in detail photoacoustic reconstruction methods and imaging systems. We also review the recent advances in nanoparticles and their in vivo applications in the field of photoacoustic imaging.

Multispectral optoacoustic tomography (MSOT) scanner for whole-body small animal imaging

A major difficulty arising from whole-body optoacoustic imaging is the long acquisition times associated with recording signals from multiple spatial projections. The acquired signals are also generally weak and the signal-to-noise-ratio is low, problems often solved by signal averaging, which complicates acquisition and increases acquisition times to an extent that makes many in vivo applications challenging or even impossible. Herein we present a fast acquisition multispectral optoacoustic tomography (MSOT) scanner for whole-body visualization of molecular markers in small animals. Multi-wavelength illumination offers the possibility to resolve exogenously administered fluorescent probes, biomarkers, and other intrinsic and exogenous chromophores. The system performance is determined in phantom experiments involving molecular probes and validated by imaging of small animals of various scales.

Quantitative three-dimensional photoacoustic tomography of the finger joints: an in vivo study

We present for the first time in vivo full three-dimensional (3-D) photoacoustic tomography (PAT) of the distal interphalangeal joint in a human subject. Both absorbed energy density and absorption coefficient images of the joint are quantitatively obtained using our finite-element-based photoacoustic image reconstruction algorithm coupled with the photon diffusion equation. The results show that major anatomical features in the joint along with the side arteries can be imaged with a 1-MHz transducer in a spherical scanning geometry. In addition, the cartilages associated with the joint can be quantitatively differentiated from the phalanx. This in vivo study suggests that the 3-D PAT method described has the potential to be used for early diagnosis of joint diseases such as osteoarthritis and rheumatoid arthritis.

Effects of different imaging models on least-squares image reconstruction accuracy in photoacoustic tomography

In the classic formulation of photoacoustic tomography (PAT), two distinct descriptions of the imaging model have been employed for developing reconstruction algorithms. We demonstrate that the numerical and statistical properties of unweighted least-squares reconstruction algorithms associated with each imaging model are generally very different. Specifically, some PAT reconstruction algorithms, including many of the iterative algorithms previously explored, do not work directly with the raw measured pressure wavefields, but rather with an integrated data function that is obtained by temporally integrating the photoacoustic wavefield. The integration modifies the statistical distribution of the data, introducing statistical correlations among samples. This change is highly significant for iterative algorithms, many of which explicitly or implicitly seek to minimize a statistical cost function. In this work, we demonstrate that iterative reconstruction by least-squares minimization yields better resolution-noise tradeoffs when working with the raw pressure data than with the integrated data commonly employed. In addition, we demonstrate that the raw-data based approach is less sensitive to certain deterministic errors, such as dc offset errors.

Measurement of the spectral directivity of optoacoustic and ultrasonic transducers with a laser ultrasonic source

Comprehensive characterization of wideband ultrasonic transducers and specifically optoacoustic detectors is achieved through the analysis of their frequency response as a function of the incident angle. The tests are performed under well-defined, repeatable operating conditions. Backillumination of a blackened, acoustically matched planar surface with a short laser pulse creates an acoustic impulse which is used as a wideband ultrasonic source. Upon illumination with a short laser pulse, the bandwidth of our source shows a -6 dB point of 12 MHz and a low-frequency roll-off around 300 kHz. Using proprietary software, we examine thoroughly the planarity of the emitted wave front within a specified amplitude cutoff and phase incoherence. Analysis of the angular dependence of the frequency response yields invaluable directivity information about the detector under study: a necessary component toward accurate optoacoustic image reconstruction and quantitative tomography. The laser ultrasonic source we developed is the main feature of our directivity measurement setup. Due to its simplicity, it can easily be adapted to various calibration devices. This paper focuses on the development and characterization of the flatness and the bandwidth of our wideband ultrasonic source.

Coregistered three-dimensional ultrasound and photoacoustic imaging system for ovarian tissue characterization

Ovarian cancer has the highest mortality of all gynecologic cancers, with a five-year survival rate of only 30% or less. Current imaging techniques are limited in sensitivity and specificity in detecting early stage ovarian cancer prior to its widespread metastasis. New imaging techniques that can provide functional and molecular contrasts are needed to reduce the high mortality of this disease. One such promising technique is photoacoustic imaging. We develop a 1280-element coregistered 3-D ultrasound and photoacoustic imaging system based on a 1.75-D acoustic array. Volumetric images over a scan range of 80 deg in azimuth and 20 deg in elevation can be achieved in minutes. The system has been used to image normal porcine ovarian tissue. This is an important step toward better understanding of ovarian cancer optical properties obtained with photoacoustic techniques. To the best of our knowledge, such data are not available in the literature. We present characterization measurements of the system and compare coregistered ultrasound and photoacoustic images of ovarian tissue to histological images. The results show excellent coregistration of ultrasound and photoacoustic images. Strong optical absorption from vasculature, especially highly vascularized corpora lutea and low absorption from follicles, is demonstrated.

Multiscale photoacoustic microscopy and computed tomography

Photoacoustic tomography (PAT) is probably the fastest growing biomedical imaging technology owing to its capability of high-resolution sensing of rich optical contrast in vivo at depths beyond the optical transport mean free path (~1 mm in the skin). Existing high-resolution optical imaging technologies, such as confocal microscopy and two-photon microscopy, have fundamentally impacted biomedicine but cannot reach such depths. Taking advantage of low ultrasonic scattering, PAT indirectly improves tissue transparency by 100 to 1000 fold and consequently enables deeply penetrating functional and molecular imaging at high spatial resolution. Further, PAT holds the promise of in vivo imaging at multiple length scales ranging from subcellular organelles to organs with the same contrast origin, an important application in multiscale systems biology research.

Photoacoustics, thermoacoustics, and acousto-optics for biomedical imaging

Recently there have been significant advances in developing hybrid techniques combining electromagnetic waves with ultrasound for biomedical imaging, namely photoacoustic, thermoacoustic, and acousto-optic (or ultrasound modulated optical) tomography. All three techniques take advantage of tissue contrast offered by electromagnetic (EM) waves, while achieving good spatial resolution in deeper tissue facilitated by ultrasound. In this review the principles of the three techniques are introduced. A description of existing experimental and image reconstruction techniques is provided. Some recent key developments are highlighted and current issues in each of the areas are discussed.

Implementation and comparison of reconstruction algorithms for two-dimensional optoacoustic tomography using a linear array

Our goal is to compare and contrast various image reconstruction algorithms for optoacoustic tomography (OAT) assuming a finite linear aperture of the kind that arises when using a linear-array transducer. Because such transducers generally have tall, narrow elements, they are essentially insensitive to out-of-plane acoustic waves, and the usually 3-D OAT problem reduces to a 2-D problem. Algorithms developed for the 3-D problem may not perform optimally in 2-D. We have implemented and evaluated a number of previously described OAT algorithms, including an exact (in 3-D) Fourier-based algorithm and a synthetic-aperture-based algorithm. We have also implemented a 2-D algorithm developed by Norton for reflection mode tomography that has not, to the best of our knowledge, been applied to OAT before. Our simulation studies of resolution, contrast, noise properties, and signal detectability measures suggest that Norton's approach-based algorithm has the best contrast, resolution, and signal detectability.

RF testbed for thermoacoustic tomography

Thermoacoustic signal excitation is a function of intrinsic tissue properties and illuminating electric field. De-ionized (DI) water is a preferred acoustic coupling medium for thermoacoustics because acoustic and electromagnetic waves propagate in DI water with very little loss. We have designed a water-filled testbed propagating a controlled electric field with respect to pulse shape, power, and polarization. Directional coupler line sections permit measurement of incident, reflected, and transmitted powers. Both S-parameters and E(y) measurement show that the electric-field distribution is relatively uniform in testbed. Comparing baseline power measurements to those taken with a test object in place yields power loss in the object, which should correlate to thermoacoustic signal strength. Moreover, power loss--and therefore thermoacoustic computerized tomography signal strength--is sensitive to the orientation of the object to the polarization of the electric field. This testbed will enable quantitative characterization of the thermoacoustic contrast mechanism in ex vivo tissue specimens.

Acoustic fingerprints of dye-labeled protein submicrosphere photoacoustic contrast agents

Dye-labeled protein microspheres, submicron in size and capable of producing thermoelastically generated ultrasound in response to laser stimulation, are presented as contrast agents for photoacoustic imaging. Incident laser energy absorbed by fluorescein isothiocyanate (FITC)-labeled elastin submicrospheres results in thermoelastically generated sound production. Plotted A-line graphs reveal a distinctive morphology and a greater than two orders of magnitude increase in signal amplitude subsequent to converting FITC elastin into submicrospheres (despite a four orders of magnitude decrease in concentration). Evidence of nonlinearity and enhancement of ultrasound backscatter indicate a potential use in contrast-enhanced harmonic imaging. Photoacoustic and ultrasound imaging of FITC-elastin submicrospheres in a water-filled phantom vessel shows enhanced contrast at low concentration and clear delineation of the phantom vessel wall.

The use of pulse synthesis for optimization of photoacoustic measurements

In this paper the use of pulse shaping in photoacoustic (PA) measurements is presented. The benefits of this approach are demonstrated by utilizing it for optimization of either the responsivity or the sensitivity of PA measurements. The optimization is based on the observation that the temporal properties of the PA effect can be represented as a linear system which can be fully characterized by its impulse response. Accordingly, the response of the PA system to an input optical pulse, whose instantaneous power is arbitrarily shaped, can be analytically predicted via a convolution between the pulse envelope and the PA impulse response. Additionally, the same formalism can be used to show that the response of the PA system to a pulse whose instantaneous power is a reversed version of the impulse response, i.e. a matched pulse, would exhibit optimal peak amplitude when compared with all other pulses with the same energy. Pulses can also be designed to optimize the sensitivity of the measurement to a variation in a specific system parameter. The use of the matched pulses can improve SNR and enable a reduction in the total optical energy required for obtaining a detectable signal. This may be important for applications where the optical power is restricted or for dynamical measurements where long integration times are prohibited. To implement this new approach, a novel PA optical setup which enabled synthesis of excitation waveforms with arbitrary temporal envelopes was constructed. The setup was based on a tunable laser source, operating in the near-IR range, and an external electro-optic modulator. Using this setup, our approach for system characterization and response prediction was tested and the superiority of the matched pulses over other common types of pulses of equal energy was demonstrated.

The application of compressed sensing for photo-acoustic tomography

Photo-acoustic (PA) imaging has been developed for different purposes, but recently, the modality has gained interest with applications to small animal imaging. As a technique it is sensitive to endogenous optical contrast present in tissues and, contrary to diffuse optical imaging, it promises to bring high resolution imaging for in vivo studies at midrange depths (3-10 mm). Because of the limited amount of radiation tissues can be exposed to, existing reconstruction algorithms for circular tomography require a great number of measurements and averaging, implying long acquisition times. Time-resolved PA imaging is therefore possible only at the cost of complex and expensive electronics. This paper suggests a new reconstruction strategy using the compressed sensing formalism which states that a small number of linear projections of a compressible image contain enough information for reconstruction. By directly sampling the image to recover in a sparse representation, it is possible to dramatically reduce the number of measurements needed for a given quality of reconstruction.

Quantitative photoacoustic measurement of tissue optical absorption spectrum aided by an optical contrast agent

In photoacoustic imaging, the intensity of photoacoustic signal induced by optical absorption in biological tissue is proportional to light energy deposition, which is the product of the absorption coefficient and the local light fluence. Because tissue optical properties are highly dependent on the wavelength, the spectrum of the local light fluence at a target tissue beneath the sample surface is different than the spectrum of the incident light fluence. Therefore, quantifying the tissue optical absorption spectrum by using a photoacoustic technique is not feasible without the knowledge of the local light fluence. In this work, a highly accurate photoacoustic measurement of the subsurface tissue optical absorption spectrum has been achieved for the first time by introducing an extrinsic optical contrast agent with known optical properties. From the photoacoustic measurements with and without the contrast agent, a quantified measurement of the chromophore absorption spectrum can be realized in a strongly scattering medium. Experiments on micro-flow vessels containing fresh canine blood buried in phantoms and chicken breast tissues were carried out in a wavelength range from 680 nm to 950 nm. Spectroscopic photoacoustic measurements of both oxygenated and deoxygenated blood specimens presented an improved match with the references when employing this technique.

Tangential resolution improvement in thermoacoustic and photoacoustic tomography using a negative acoustic lens

We developed a novel concept of using a negative acoustic lens to increase the acceptance angle of an unfocused large-area ultrasonic transducer (detector), leading to more than twofold improvement of the tangential resolution in both thermoacoustic and photoacoustic tomography. In both thermoacoustic and photoacoustic tomography, for a given transducer bandwidth, the aperture size of the detector affects the tangential resolution greatly when the object of interest is near the detector surface. We were able to overcome such tangential resolution deterioration by attaching an acoustic concave lens, made of acrylic in front of the flat detector surface. We then quantified the tangential resolution improvement using phantom images. We also showed that the use of the negative lens preserves the shape of an object after the image is reconstructed.

Photoacoustic tomography with a single detector in a reverberant cavity

In conventional biomedical photoacoustic tomography (PAT), ultrasonic pulses generated through the absorption of nanosecond pulses of near-infrared light are recorded over an array of detectors and used to recover an image of the initial acoustic pressure distribution within soft tissue. This image is related to the tissue optical coefficients and therefore carries information about the tissue physiology. For high resolution imaging, a large-area detector array with a high density of small, sensitive elements is required. Such arrays can be expensive, so reverberant-field PAT has been suggested as a means of obtaining PAT images using arrays with a smaller number of detectors. By recording the reflections from an acoustically reverberant cavity surrounding the sample, in addition to the primary acoustic pulse, sufficient information may be captured to allow an image to be reconstructed without the need for a large-area array. An initial study using two-dimensional simulations was performed to assess the feasibility of using a single detector for PAT. It is shown that reverberant-field data recorded at a single detector are sufficient to reconstruct the initial pressure distribution accurately, so long as the shape of the reverberant cavity makes it ray-chaotic. The practicalities of such an approach to photoacoustic imaging are discussed.

Laser optoacoustic imaging system for detection of breast cancer

We designed, fabricated and tested the laser optoacoustic imaging system for breast cancer detection (LOIS-64), which fuses optical and acoustic imaging techniques in one modality by utilizing pulsed optical illumination and ultrawide-band ultrasonic detection of resulting optoacoustic (OA) signals. The system was designed to image a single breast slice in craniocaudal or mediolateral projection with an arc-shaped array of 64 ultrawide-band acoustic transducers. The system resolution on breast phantoms was at least 0.5 mm. The single-channel sensitivity of 1.66 mVPa was estimated to be sufficient for single-pulse imaging of 6 to 11 mm tumors through the whole imaging slice of the breast. The implemented signal processing using the wavelet transform allowed significant reduction of the low-frequency (LF) acoustic noise, allowed localization of the optoacoustic signals from tumors, and enhanced the contrast and sharpened the boundaries of the optoacoustic images of the tumors. During the preliminary clinical studies on 27 patients, the LOIS-64 was able to visualize 18 out of 20 malignant lesions suspected from mammography and ultrasound images and confirmed by the biopsy performed after the optoacoustic tomography (OAT) procedure.

Picomolar sensitivity MRI and photoacoustic imaging of cobalt nanoparticles

Multimodality imaging based on complementary detection principles has broad clinical applications and promises to improve the accuracy of medical diagnosis. This means that a tracer particle advantageously incorporates multiple functionalities into a single delivery vehicle. In the present work, we explore a unique combination of MRI and photoacoustic tomography (PAT) to detect picomolar concentrations of nanoparticles. The nanoconstruct consists of ferromagnetic (Co) particles coated with gold (Au) for biocompatibility and a unique shape that enables optical absorption over a broad range of frequencies. The end result is a dual-modality probe useful for the detection of trace amounts of nanoparticles in biological tissues, in which MRI provides volume detection, whereas PAT performs edge detection.

Estimating chromophore distributions from multiwavelength photoacoustic images

Biomedical photoacoustic tomography (PAT) can provide qualitative images of biomedical soft tissue with high spatial resolution. However, whether it is possible to give accurate quantitative estimates of the spatially varying concentrations of the sources of photoacoustic contrast-endogenous or exogenous chromophores-remains an open question. Even if the chromophores' absorption spectra are known, the problem is nonlinear and ill-posed. We describe a framework for obtaining such quantitative estimates. When the optical scattering distribution is known, adjoint and gradient-based optimization techniques can be used to recover the concentration distributions of the individual chromophores that contribute to the overall tissue absorption. When the scattering distribution is unknown, prior knowledge of the wavelength dependence of the scattering is shown to be sufficient to overcome the absorption-scattering nonuniqueness and allow both distributions of chromophore concentrations and scattering to be recovered from multiwavelength photoacoustic images.

Poly(lactic-co-glycolic) acid as a carrier for imaging contrast agents

PURPOSE

With the broadening field of nanomedicine poised for future molecular level therapeutics, nano- and microparticles intended for the augmentation of either single- or multimodal imaging are created with PLGA as the chief constituent and carrier.

METHODS

Emulsion techniques were used to encapsulate hydrophilic and hydrophobic imaging contrast agents in PLGA particles. The imaging contrast properties of these PLGA particles were further enhanced by reducing silver onto the PLGA surface, creating a silver cage around the polymeric core.

RESULTS

The MRI contrast agent Gd-DTPA and the exogenous dye rhodamine 6G were both encapsulated in PLGA and shown to enhance MR and fluorescence contrast, respectively. The silver nanocage built around PLGA nanoparticles exhibited strong near infrared light absorbance properties, making it a suitable contrast agent for optical imaging strategies such as photoacoustic imaging.

CONCLUSIONS

The biodegradable polymer PLGA is an extremely versatile nano- and micro-carrier for several imaging contrast agents with the possibility of targeting diseased states at a molecular level.

Photoacoustic recovery of an absolute optical absorption coefficient with an exact solution of a wave equation

We describe a new approach to recover an absolute optical absorption coefficient from measured photoacoustic signals based on an exact solution of the wave equation. We present a numerical model to describe the detection of photoacoustic waves by introducing a point spread function of a photoacoustic measuring system. This model is employed in an iterative fitting procedure in order to recover an absolute optical absorption coefficient. In addition, the fitting procedure utilizes the ratios of the amplitudes of the characteristic peaks of photoacoustic signals as the input, as opposed to using the total measured photoacoustic data. Such an arrangement increases computational speed and reduces the influence of signals from neighboring objects. The proposed method is then verified both numerically and experimentally on the gel samples, from which the accuracy is found to be in the range -5.2% to 4.8%.

Deconvolution reconstruction of full-view and limited-view photoacoustic tomography: a simulation study

Although many algorithms are available for full-view photoacoustic tomography (PAT), no exact and stable algorithm for limited-view PAT has been proposed. In this paper the deconvolution reconstruction (DR) algorithm is proposed for both full-view and limited-view PAT. In the DR algorithm, first a new function is constructed from detected photoacoustic signals and approximately simplified, and then the tissue's electromagnetic absorption is derived from this function on the basis of Fourier-based deconvolution. Computer simulations are carried out to compare the DR algorithm with two popular PAT algorithms, the time-domain reconstruction (TDR) and the filtered back projection (FBP). Although the error of the DR algorithm increases with the size of the detected object, it is shown that the DR algorithm has good precision and strong robustness to noise in the full-view PAT, nearly equivalent to the TDR and FBP. Yet the DR algorithm is more than ten times faster in computation speed. In the limited-view PAT, the DR is superior to the TDR and FBP in terms of both accuracy and robustness to noise.