Fast photoacoustic imaging system based on 320-element linear transducer array

Evaluation of 10 current image reconstruction algorithms for linear array photoacoustic imaging

Various reconstruction algorithms have been implemented for linear array photoacoustic imaging systems with the goal of accurately reconstructing the strength absorbers within the tissue being imaged. Since the existing algorithms have been introduced by different research groups and the context of performance evaluation was not consistent, it is difficult to make a fair comparison between them. In this study, we systematically compared the performance of 10 published image reconstruction algorithms (DAS, UBP, pDAS, DMAS, MV, EIGMV, SLSC, GSC, TR, and FD) using in-vitro phantom data. Evaluations were conducted based on lateral resolution of the reconstructed images, computational time, target detectability, and noise sensitivity. We anticipate the outcome of this study will assist researchers in selecting appropriate algorithms for their linear array PA imaging applications.

Carbon-Coated Magnetic Nanoparticle Dedicated to MRI/Photoacoustic Imaging of Tumor in Living Mice

Multimodality imaging can reveal complementary anatomic and functional information as they exploit different contrast mechanisms, which has broad clinical applications and promises to improve the accuracy of tumor diagnosis. Accordingly, to attain the particular goal, it is critical to exploit multimodal contrast agents. In the present work, we develop novel cobalt core/carbon shell-based nanoparticles (Cobalt at carbon NPs) with both magnetization and light absorption properties for dual-modality magnetic resonance imaging (MRI) and photoacoustic imaging (PAI). The nanoparticle consists of ferromagnetic cobalt particles coated with carbon for biocompatibility and optical absorption. In addition, the prepared Cobalt at carbon NPs are characterized by transmission electron microscope (TEM), visible-near-infrared spectra, Raman spectrum, and X-ray powder diffraction for structural analysis. Experiments verify that Cobalt at carbon NPs have been successfully constructed and the designed Cobalt at carbon NPs can be detected by both MRI and PAI in vitro and in vivo. Importantly, intravenous injection of Cobalt at carbon NPs into glioblastoma-bearing mice led to accumulation and retention of Cobalt at carbon NPs in the tumors. Using such a multifunctional probe, MRI can screen rapidly to identify potential lesion locations, whereas PAI can provide high-resolution morphological structure and quantitative information of the tumor. The Cobalt at carbon NPs are likely to become a promising candidate for dual-modality MRI/PAI of the tumor.

MEMS Ultrasound Transducers for Endoscopic Photoacoustic Imaging Applications

Photoacoustic imaging (PAI) is drawing extensive attention and gaining rapid development as an emerging biomedical imaging technology because of its high spatial resolution, large imaging depth, and rich optical contrast. PAI has great potential applications in endoscopy, but the progress of endoscopic PAI was hindered by the challenges of manufacturing and assembling miniature imaging components. Over the last decade, microelectromechanical systems (MEMS) technology has greatly facilitated the development of photoacoustic endoscopes and extended the realm of applicability of the PAI. As the key component of photoacoustic endoscopes, micromachined ultrasound transducers (MUTs), including piezoelectric MUTs (pMUTs) and capacitive MUTs (cMUTs), have been developed and explored for endoscopic PAI applications. In this article, the recent progress of pMUTs (thickness extension mode and flexural vibration mode) and cMUTs are reviewed and discussed with their applications in endoscopic PAI. Current PAI endoscopes based on pMUTs and cMUTs are also introduced and compared. Finally, the remaining challenges and future directions of MEMS ultrasound transducers for endoscopic PAI applications are given.

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.

Review of cost reduction methods in photoacoustic computed tomography

Photoacoustic Computed Tomography (PACT) is a major configuration of photoacoustic imaging, a hybrid noninvasive modality for both functional and molecular imaging. PACT has rapidly gained importance in the field of biomedical imaging due to superior performance as compared to conventional optical imaging counterparts. However, the overall cost of developing a PACT system is one of the challenges towards clinical translation of this novel technique. The cost of a typical commercial PACT system originates from optical source, ultrasound detector, and data acquisition unit. With growing applications of photoacoustic imaging, there is a tremendous demand towards reducing its cost. In this review article, we have discussed various approaches to reduce the overall cost of a PACT system, and provided a cost estimation to build a low-cost PACT system.

GPU-accelerated Double-stage Delay-multiply-and-sum Algorithm for Fast Photoacoustic Tomography Using LED Excitation and Linear Arrays

Double-stage delay-multiply-and-sum (DS-DMAS) is an algorithm proposed for photoacoustic image reconstruction. The DS-DMAS algorithm offers a higher contrast than conventional delay-and-sum and delay-multiply and-sum but at the expense of higher computational complexity. Here, we utilized a compute unified device architecture (CUDA) graphics processing unit (GPU) parallel computation approach to address the high complexity of the DS-DMAS for photoacoustic image reconstruction generated from a commercial light-emitting diode (LED)-based photoacoustic scanner. In comparison with a single-threaded central processing unit (CPU), the GPU approach increased speeds by nearly 140-fold for 1024 × 1024 pixel image; there was no decrease in accuracy. The proposed implementation makes it possible to reconstruct photoacoustic images with frame rates of 250, 125, and 83.3 when the images are 64 × 64, 128 × 128, and 256 × 256, respectively. Thus, DS-DMAS can be efficiently used in clinical devices when coupled with CUDA GPU parallel computation.

Photoacoustic Imaging with Capacitive Micromachined Ultrasound Transducers: Principles and Developments

Photoacoustic imaging (PAI) is an emerging imaging technique that bridges the gap between pure optical and acoustic techniques to provide images with optical contrast at the acoustic penetration depth. The two key components that have allowed PAI to attain high-resolution images at deeper penetration depths are the photoacoustic signal generator, which is typically implemented as a pulsed laser and the detector to receive the generated acoustic signals. Many types of acoustic sensors have been explored as a detector for the PAI including Fabry-Perot interferometers (FPIs), micro ring resonators (MRRs), piezoelectric transducers, and capacitive micromachined ultrasound transducers (CMUTs). The fabrication technique of CMUTs has given it an edge over the other detectors. First, CMUTs can be easily fabricated into given shapes and sizes to fit the design specifications. Moreover, they can be made into an array to increase the imaging speed and reduce motion artifacts. With a fabrication technique that is similar to complementary metal-oxide-semiconductor (CMOS), CMUTs can be integrated with electronics to reduce the parasitic capacitance and improve the signal to noise ratio. The numerous benefits of CMUTs have enticed researchers to develop it for various PAI purposes such as photoacoustic computed tomography (PACT) and photoacoustic endoscopy applications. For PACT applications, the main areas of research are in designing two-dimensional array, transparent, and multi-frequency CMUTs. Moving from the table top approach to endoscopes, some of the different configurations that are being investigated are phased and ring arrays. In this paper, an overview of the development of CMUTs for PAI is presented.

Fractional Regularization to Improve Photoacoustic Tomographic Image Reconstruction

Photoacoustic tomography involves reconstructing the initial pressure rise distribution from the measured acoustic boundary data. The recovery of the initial pressure rise distribution tends to be an ill-posed problem in the presence of noise and when limited independent data is available, necessitating regularization. The standard regularization schemes include Tikhonov, l₁ -norm, and total-variation. These regularization schemes weigh the singular values equally irrespective of the noise level present in the data. This paper introduces a fractional framework to weigh the singular values with respect to a fractional power. This fractional framework was implemented for Tikhonov, l₁ -norm, and total-variation regularization schemes. Moreover, an automated method for choosing the fractional power was also proposed. It was shown theoretically and with numerical experiments that the fractional power is inversely related to the data noise level for fractional Tikhonov scheme. The fractional framework outperforms the standard regularization schemes, Tikhonov, l₁ -norm, and total-variation by 54% in numerical simulations, experimental phantoms, and in vivo rat data in terms of observed contrast/signal-to-noise-ratio of the reconstructed images.

Wave front analysis for enhanced time-domain beamforming of point-like targets in optoacoustic imaging using a linear array

Using linear array transducers in combination with state-of-the-art multichannel electronics allows to perform optoacoustic imaging with frame rates only limited by the laser pulse repetition frequency and the acoustic time of flight. However, characteristic image artefacts resulting from the limited view and a lower SNR when compared to systems based on single-element focused transducers represent a burden for the clinical acceptance of the technology. In this paper, we present a new method for the improvement of image quality based on the analysis of the signal amplitudes along summation curves during the delay-and-sum beamforming process (DAS). The algorithm compares amplitude distributions along wave fronts with theoretical patterns from optoacoustic point sources. The method was validated on simulated and experimental phantom as well as in-vivo data. An improvement of the lateral resolution by more than a factor of two when comparing conventional DAS and our approach could be shown (numeric and experimental phantom data). For instance, on experimental data from a wire phantom, a PSF in the range of 0.18-0.22 mm was obtained with our approach against 0.48 mm for standard DAS. Furthermore, the SNR of a subcutaneous vessel 2.5 mm below the skin surface was improved by about 30 dB when compared to standard DAS.

Toward Intra-Operative Prostate Photoacoustic Imaging: Configuration Evaluation and Implementation Using the da Vinci Research Kit

We compare different possible scanning geometries for prostate photoacoustic tomography (PAT) while considering a realistic reconstruction scenario in which the limited view of the prostate and the directivity effect of the transducer are considered. Simulations and experiments confirm that an intra-operative configuration in which the photoacoustic signal is received by a pickup transducer from the anterior surface of the prostate provides the best approach. We propose a PAT acquisition system that includes a da Vinci system controlled by the da Vinci Research Kit, an illumination laser, and an ultrasound machine with parallel data acquisition. The robot maneuvers the pickup transducer to form a cylindrical detection surface around the prostate. The robot is programmed to acquire trajectories in which the transducer face is parallel to and oriented toward a rotational tomography axis, while the laser is fired and PAT data are collected at regular intervals. We present our initial images acquired with this novel system.

The double-stage delay-multiply-and-sum image reconstruction method improves imaging quality in a LED-based photoacoustic array scanner

Light-emitting diode-based photoacoustic imaging is more compact and affordable than laser-based systems, but it has low power and hence a high number of replicates. Here, we describe double-stage delay-multiply-and-sum (DS-DMAS) to improve image quality collected on a LED-based scanner. DS-DMAS was evaluated experimentally using point targets (in different laterals and depths) as well as a hair and a rabbit eye. This algorithm can compensate for the low SNR of LED-based systems and offer better lateral resolution of about 60%, 25%, higher contrast ratio of about 97%, 34%, and better full-width-half-maximum of about 60%, 25%, versus delay-and-sum) and delay-multiply-and-sum, respectively. More importantly, DS-DMAS offers this using a smaller number of frames (only 2% of all the frames). These results indicate that DS-DMAS might be a valuable tool in the translation of LED-based and other low power PAI systems.

Polarized photoacoustic microscopy for vectorial-absorption-based anisotropy detection

We proposed polarized photoacoustic microscopy (PPAM) for quantitative detection of a target's microscopic anisotropy based on the vectorial optical absorption by applying four linearly polarized laser beams as excitation sources. Compared to conventional photoacoustic microscopy that treats targets as isotropic absorbers, PPAM allows us to quantitatively detect the target's anisotropic features beyond optical absorption with a newly proposed parameter valued between 0 and 1. The feasibility of the method was validated by dichroic phantoms. The dichroic compound eyes of mantis shrimps were imaged in situ to demonstrate the method's capability for quantitative three-dimensional biological imaging. The PPAM method provides an effective and straightforward strategy for tissue polarimetry, prefiguring great potential for biological imaging and material inspection.

3D photoacoustic imaging using Hadamard-bias encoding with a crossed electrode relaxor array

A novel 3D photoacoustic imaging technique is experimentally demonstrated using a 64×64 element bias-sensitive crossed-electrode relaxor array. This technique allows for large 2D arrays to receive across all elements while using minimal channel counts. Hadamard-bias patterns are applied to column electrodes while signals are measured from row electrodes. Photoacoustic signals are measured from a crossed-wire phantom in an intralipid scattering medium. The Hadamard-bias-encoded imaging scheme showed a signal-to-noise (SNR) of 25.3 dB, while the single-column biasing strategy (or identity-matrix-bias pattern) showed a SNR of 8.8 dB.

Internal-illumination photoacoustic computed tomography

We report a photoacoustic computed tomography (PACT) system using a customized optical fiber with a cylindrical diffuser to internally illuminate deep targets. The traditional external light illumination in PACT usually limits the penetration depth to a few centimeters from the tissue surface, mainly due to strong optical attenuation along the light propagation path from the outside in. By contrast, internal light illumination, with external ultrasound detection, can potentially detect much deeper targets. Different from previous internal illumination PACT implementations using forward-looking optical fibers, our internal-illumination PACT system uses a customized optical fiber with a 3-cm-long conoid needle diffuser attached to the fiber tip, which can homogeneously illuminate the surrounding space and substantially enlarge the field of view. We characterized the internal illumination distribution and PACT system performance. We performed tissue phantom and in vivo animal studies to further demonstrate the superior imaging depth using internal illumination over external illumination. We imaged a 7.5-cm-deep leaf target embedded in optically scattering medium and the beating heart of a mouse overlaid with 3.7-cm-thick chicken tissue. Our results have collectively demonstrated that the internal light illumination combined with external ultrasound detection might be a useful strategy to improve the penetration depth of PACT in imaging deep organs of large animals and humans.

Calibrated Linear Array-Driven Photoacoustic/Ultrasound Tomography

The anisotropic resolution of linear arrays, tools that are widely used in diagnostics, can be overcome by compounding approaches. We investigated the ability of a recently developed calibration and a novel algorithm to determine the actual radial transducer array distance and its misalignment (tilt) with respect to the center of rotation in a 2-D and 3-D tomographic setup. By increasing the time-of-flight accuracy, we force in-phase summation during the reconstruction. Our setup is composed of a linear transducer and a rotation and translation axis enabling multidimensional imaging in ultrasound and photoacoustic mode. Our approach is validated on phantoms and young mice ex vivo. The results indicate that application of the proposed analytical calibration algorithms prevents image artifacts. The spatial resolution achieved was 160 and 250 μm in photoacoustic mode of 2-D and 3-D tomography, respectively.

Synthetic-aperture based photoacoustic re-beamforming (SPARE) approach using beamformed ultrasound data

Photoacoustic (PA) imaging has been developed for various clinical and pre-clinical applications, and acquiring pre-beamformed channel data is necessary to reconstruct these images. However, accessing these pre-beamformed channel data requires custom hardware to enable parallel beamforming, and is available for a limited number of research ultrasound platforms. To broaden the impact of clinical PA imaging, our goal is to devise a new PA reconstruction approach that uses ultrasound post-beamformed radio frequency (RF) data rather than raw channel data, because this type of data is readily available in both clinical and research ultrasound systems. In our proposed Synthetic-aperture based photoacoustic re-beamforming (SPARE) approach, post-beamformed RF data from a clinical ultrasound scanner are considered as input data for an adaptive synthetic aperture beamforming algorithm. When receive focusing is applied prior to obtaining these data, the focal point is considered as a virtual element, and synthetic aperture beamforming is implemented assuming that the photoacoustic signals are received at the virtual element. The resolution and SNR obtained with the proposed method were compared to that obtained with conventional delay-and-sum beamforming with 99.87% and 91.56% agreement, respectively. In addition, we experimentally demonstrated feasibility with a pulsed laser diode setup. Results indicate that the post-beamformed RF data from any commercially available ultrasound platform can potentially be used to create PA images.

Comparison of Deconvolution Filters for Photoacoustic Tomography

In this work, we compare the merits of three temporal data deconvolution methods for use in the filtered backprojection algorithm for photoacoustic tomography (PAT). We evaluate the standard Fourier division technique, the Wiener deconvolution filter, and a Tikhonov L-2 norm regularized matrix inversion method. Our experiments were carried out on subjects of various appearances, namely a pencil lead, two man-made phantoms, an in vivo subcutaneous mouse tumor model, and a perfused and excised mouse brain. All subjects were scanned using an imaging system with a rotatable hemispherical bowl, into which 128 ultrasound transducer elements were embedded in a spiral pattern. We characterized the frequency response of each deconvolution method, compared the final image quality achieved by each deconvolution technique, and evaluated each method's robustness to noise. The frequency response was quantified by measuring the accuracy with which each filter recovered the ideal flat frequency spectrum of an experimentally measured impulse response. Image quality under the various scenarios was quantified by computing noise versus resolution curves for a point source phantom, as well as the full width at half maximum (FWHM) and contrast-to-noise ratio (CNR) of selected image features such as dots and linear structures in additional imaging subjects. It was found that the Tikhonov filter yielded the most accurate balance of lower and higher frequency content (as measured by comparing the spectra of deconvolved impulse response signals to the ideal flat frequency spectrum), achieved a competitive image resolution and contrast-to-noise ratio, and yielded the greatest robustness to noise. While the Wiener filter achieved a similar image resolution, it tended to underrepresent the lower frequency content of the deconvolved signals, and hence of the reconstructed images after backprojection. In addition, its robustness to noise was poorer than that of the Tikhonov filter. The performance of the Fourier filter was found to be the poorest of all three methods, based on the reconstructed images' lowest resolution (blurriest appearance), generally lowest contrast-to-noise ratio, and lowest robustness to noise. Overall, the Tikhonov filter was deemed to produce the most desirable image reconstructions.

Linear-array-based photoacoustic imaging of human microcirculation with a range of high frequency transducer probes

Photoacoustic imaging (PAI) with a linear-array-based probe can provide a convenient means of imaging the human microcirculation within its native structural context and adds functional information. PAI using a multielement linear transducer array combined with multichannel collecting system was used for in vivo volumetric imaging of the blood microcirculation, the total concentration of hemoglobin (HbT), and the hemoglobin oxygen saturation (sO₂) within human tissue. Three-dimensional (3-D) PA and ultrasound (US) volumetric scans were acquired from the forearm skin by linearly translating the transducer with a stepper motor over a region of interest, while capturing two-dimensional images using 15, 21, and 40 MHz frequency transducer probes. For the microvasculature imaging, PA images were acquired at 800- and 1064-nm wavelengths. For the HbT and sO₂ estimates, PA images were collected at 750- and 850-nm wavelengths. 3-D microcirculation, HbT, and sO₂ maps of the forearm skin were obtained from normal subjects. The linear-array-based PAI has been found promising in terms of resolution, imaging depth, and imaging speed for in vivo microcirculation imaging within human skin. We believe that a reflection type probe, similar to existing clinical US probes, is most likely to succeed in real clinical applications. Its advantages include ease of use, speed, and familiarity for radiographers and clinicians.

Basis pursuit deconvolution for improving model-based reconstructed images in photoacoustic tomography

The model-based image reconstruction approaches in photoacoustic tomography have a distinct advantage compared to traditional analytical methods for cases where limited data is available. These methods typically deploy Tikhonov based regularization scheme to reconstruct the initial pressure from the boundary acoustic data. The model-resolution for these cases represents the blur induced by the regularization scheme. A method that utilizes this blurring model and performs the basis pursuit deconvolution to improve the quantitative accuracy of the reconstructed photoacoustic image is proposed and shown to be superior compared to other traditional methods via three numerical experiments. Moreover, this deconvolution including the building of an approximate blur matrix is achieved via the Lanczos bidagonalization (least-squares QR) making this approach attractive in real-time.

Light in and sound out: emerging translational strategies for photoacoustic imaging

Photoacoustic imaging (PAI) has the potential for real-time molecular imaging at high resolution and deep inside the tissue, using nonionizing radiation and not necessarily depending on exogenous imaging agents, making this technique very promising for a range of clinical applications. The fact that PAI systems can be made portable and compatible with existing imaging technologies favors clinical translation even more. The breadth of clinical applications in which photoacoustics could play a valuable role include: noninvasive imaging of the breast, sentinel lymph nodes, skin, thyroid, eye, prostate (transrectal), and ovaries (transvaginal); minimally invasive endoscopic imaging of gastrointestinal tract, bladder, and circulating tumor cells (in vivo flow cytometry); and intraoperative imaging for assessment of tumor margins and (lymph node) metastases. In this review, we describe the basics of PAI and its recent advances in biomedical research, followed by a discussion of strategies for clinical translation of the technique.

Structural and functional photoacoustic molecular tomography aided by emerging contrast agents

Photoacoustic tomography (PAT) can offer structural, functional and molecular contrasts at scalable observation level. By ultrasonically overcoming the strong optical scattering, this imaging technology can reach centimeters penetration depth while retaining high spatial resolution in biological tissue. Recent extensive research has been focused on developing new contrast agents to improve the imaging sensitivity, specificity and efficiency. These emerging materials have substantially accelerated PAT applications in signal sensing, functional imaging, biomarker labeling and therapy monitoring etc. Here, the potentials of different optical probes as PAT contrast agents were elucidated. We first describe the instrumental embodiments and the measured functional parameters, then focus on emerging contrast agent-based PAT applications, and finally discuss the challenges and prospects.

Biomedical photoacoustics in China

During the last decade, along with its explosive growth globally, biomedical photoacoustics has become a rapidly growing research field in China as well. In particular, photoacoustic tomography (PAT), capable of imaging intact biological tissue in vivo at great depths, has generated intense interest among Chinese researchers. This review briefly summarizes the current status and recent progress of the research in PAT in China. The focus is on the technology development and biomedical applications of three representative embodiments of PAT: photoacoustic microscopy, photoacoustic computed tomography, and photoacoustic endoscopy. In addition, recent development and studies in other related areas are also reviewed shortly.

Impact of inhomogeneous optical scattering coefficient distribution on recovery of optical absorption coefficient maps using tomographic photoacoustic data

We present a study through extensive simulation that considers the impact of inhomogeneous optical scattering coefficient distribution on recovery of optical absorption coefficient maps using tomographic photoacoustic data collected from media mimicking breast tissue. We found that while the impact of scattering heterogeneities/targets is modest on photoacoustic recovery of optical absorption coefficients, the impact of scattering contrast caused by adipose tissue, a layer of normal tissue along the boundary of the breast, is dramatic on reconstruction of optical absorption coefficients using photoacoustic data-up to 25.8% relative error in recovering the absorption coefficient is estimated in such cases. To overcome this problem, we propose a new method to enhance photoacoustic recovery of the optical absorption coefficient in heterogeneous media by considering inhomogeneous scattering coefficient distribution provided by diffuse optical tomography (DOT). Results from extensive simulations show that photoacoustic recovery of absorption coefficient maps can be improved considerably with a priori scattering information from DOT.

Three-dimensional thermoacoustic imaging for early breast cancer detection

PURPOSE

Microwave-induced thermoacoustic tomography (TAT) is a noninvasive modality based on the differences in microwave absorption of various biological tissues. In this paper, the feasibility of the early breast tumor detection by TAT system has been discussed and validated experimentally.

METHODS

A fast TAT system, which based on three 128-elements transducers, a 384-64ch switch and a parallel data acquisition system (DAS), was developed to reconstruct the three-dimensional (3D) image of a breast model with similar microwave absorption coefficient to breast tissue. A novel method to explore the ability of TAT system to distinguish absorption coefficient was introduced and the minimum absorption coefficient difference that can be distinguished clearly by our TAT system is 12 m(-1).

RESULTS

The potential applications of the TAT system were clearly demonstrated by successfully mapping breast model with mimicked tumors and microcalcification. An imaging experiment of human breast tumor embedding in the breast model was performed and the tumor was visualized by the 3D thermoacoustic volume.

CONCLUSIONS

The thermoacoustic images match well with the samples and achieve penetration depth of 6 cm. The experimental results indicate that TAT has a great potential to be used for detecting early-stage breast cancers with high contrast and high resolution.

Improving image quality by accounting for changes in water temperature during a photoacoustic tomography scan

The emerging field of photoacoustic tomography is rapidly evolving with many new system designs and reconstruction algorithms being published. Many systems use water as a coupling medium between the scanned object and the ultrasound transducers. Prior to a scan, the water is heated to body temperature to enable small animal imaging. During the scan, the water heating system of some systems is switched off to minimize the risk of bubble formation, which leads to a gradual decrease in water temperature and hence the speed of sound. In this work, we use a commercially available scanner that follows this procedure, and show that a failure to model intra-scan temperature decreases as small as 1.5°C leads to image artifacts that may be difficult to distinguish from true structures, particularly in complex scenes. We then improve image quality by continuously monitoring the water temperature during the scan and applying variable speed of sound corrections in the image reconstruction algorithm. While upgrading to an air bubble-free heating pump and keeping it running during the scan could also solve the changing temperature problem, we show that a software correction for the temperature changes provides a cost-effective alternative to a hardware upgrade. The efficacy of the software corrections was shown to be consistent across objects of widely varying appearances, namely physical phantoms, ex vivo tissue, and in vivo mouse imaging. To the best of our knowledge, this is the first study to demonstrate the efficacy of modeling temporal variations in the speed of sound during photoacoustic scans, as opposed to spatial variations as focused on by previous studies. Since air bubbles pose a common problem in ultrasonic and photoacoustic imaging systems, our results will be useful to future small animal imaging studies that use scanners with similarly limited heating units.

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.

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.

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.

Photoacoustic imaging based on MEMS mirror scanning

A microelectromechanical systems (MEMS)-based photoacoustic imaging system is reported for the first time. In this system, the MEMS-based light scanning subsystem and a ring-shaped polyvinylidene fluoride (PVDF) transducer are integrated into a miniaturized probe that is capable of three-dimensional (3D) photoacoustic imaging. It is demonstrated that the imaging system is able to image small objects embedded in phantom materials and in chicken and to in vivo visualize blood vessels under the skin of a human hand.

Video rate optoacoustic tomography of mouse kidney perfusion

Optoacoustic tomography can visualize optical contrast in tissues while capitalizing on the advantages of ultrasound, such as high spatial resolution and fast imaging capabilities. We report a novel multispectral optoacoustic tomography system for deep tissue small animal imaging. The previously undocumented capacity of whole-body optoacoustic tomography at a video rate has been demonstrated by visualizing mouse kidney perfusion using Indocyanine Green in vivo.

Photoacoustic monitoring of granulation tissue grown in a grafted artificial dermis on rat skin

In this study, we investigated the validity of photoacoustic (PA) measurement for monitoring granulation tissue and hence adhesion of grafted artificial dermis (AD). A 2.5 cm x 2.5 cm, 3-mm-thick AD composed of an atelocollagen sponge sheet and a silicone film was grafted on a full-thickness open wound in rat dorsal skin. The grafted AD was irradiated with low-energy, 532-nm nanosecond laser pulses to photoacoustically excite blood in neovascularities, and the PA signals induced were measured using a piezoelectric transducer as a function of postgrafting time. The PA signals were compared with results of laser Doppler imaging and histological analysis. We found a significant correlation between the depths of the first or shallowest PA signal peaks and the depths of granulation tissues estimated from histology with hematoxylin & eosin staining (R=0.951, p<0.05). There was also a significant correlation between the amplitudes of the first PA signal peaks and densities of CD31-positive cells evaluated from histology with immunohistochemical staining (R=0.859, p<0.05). With laser Doppler imaging, no clear perfusion signals were observed, which is attributable to a high light scattering loss in ADs. These findings suggest the validity of PA measurement for monitoring the adhesion of grafted ADs.

Photoacoustic imaging and characterization of the microvasculature

Photoacoustic (optoacoustic) tomography, combining optical absorption contrast and highly scalable spatial resolution (from micrometer optical resolution to millimeter acoustic resolution), has broken through the fundamental penetration limit of optical ballistic imaging modalities-including confocal microscopy, two-photon microscopy, and optical coherence tomography-and has achieved high spatial resolution at depths down to the diffusive regime. Optical absorption contrast is highly desirable for microvascular imaging and characterization because of the presence of endogenous strongly light-absorbing hemoglobin. We focus on the current state of microvascular imaging and characterization based on photoacoustics. We first review the three major embodiments of photoacoustic tomography: microscopy, computed tomography, and endoscopy. We then discuss the methods used to characterize important functional parameters, such as total hemoglobin concentration, hemoglobin oxygen saturation, and blood flow. Next, we highlight a few representative applications in microvascular-related physiological and pathophysiological research, including hemodynamic monitoring, chronic imaging, tumor-vascular interaction, and neurovascular coupling. Finally, several potential technical advances toward clinical applications are suggested, and a few technical challenges in contrast enhancement and fluence compensation are summarized.

In vivo detection and imaging of low-density foreign body with microwave-induced thermoacoustic tomography

Radiography or computed tomography is the most widely available imaging tool for foreign body detection. However, the detectability of low-density substances by x ray is very poor when located in soft tissues. Various dielectric loss factors of foreign bodies contribute great microwave absorption heterogeneity compared with the surrounding tissue. A fast thermoacoustic tomography system at 1.2 and 6 GHz was developed to detect foreign targets in small animals. The 6 GHz system had a much higher signal-to-noise ratio in near-surface imaging but smaller imaging depth than the 1.2 GHz system. The effects of microwave distribution inhomogeneity on nonuniform excitation of acoustic pressure were studied and a corresponding calibration algorithm for image distortion was provided and experimentally examined. Thermoacoustic images of radiolucent objects including glass fiber, wood, and bamboo hidden in phantom and residual in living mice were compared with radiography and ultrasonography. Good contrast was obtained between the foreign bodies and the tissue surrounding it, and the location and size of the lesion targets in thermoacoustic images were in good agreement with the actual sample. The experimental results demonstrate that thermoacoustic tomography may become the ideal modality for radiolucent foreign body detection and imaging in animals and human.

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.

Photoacoustic tomography of small animal brain with a curved array transducer

We present the application of a curved array photoacoustic tomographic imaging system that can provide rapid, high-resolution photoacoustic imaging of small animal brains. The system is optimized to produce a B-mode, 90-deg field-of-view image at sub-200-microm resolution at a frame rate of approximately 1 frame/second when a 10-Hz pulse repetition rate laser is employed. By rotating samples, a complete 360-deg scan can be achieved within 15 s. In previous work, two-dimensional (2-D) ex vivo mouse brain cortex imaging has been reported. We report three-dimensional (3-D) small animal brain imaging obtained with the curved array system. The results are presented as a series of 2-D cross-sectional images. Besides structural imaging, the blood oxygen saturation of the animal brain cortex is also measured in vivo. In addition, the system can measure the time-resolved relative changes in blood oxygen saturation level in the small animal brain cortex. Last, ultrasonic gel coupling, instead of the previously adopted water coupling, is conveniently used in near-real-time 2-D imaging.

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.

Quantitative photoacoustic tomography

In this paper, several algorithms that allow for quantitative photoacoustic reconstruction of tissue optical, acoustic and physiological properties are described in a finite-element method based framework. These quantitative reconstruction algorithms are compared, and the merits and limitations associated with these methods are discussed. In addition, a multispectral approach is presented for concurrent reconstructions of multiple parameters including deoxyhaemoglobin, oxyhaemoglobin and water concentrations as well as acoustic speed. Simulation and in vivo experiments are used to demonstrate the effectiveness of the reconstruction algorithms presented.

A real-time photoacoustic tomography system for small animals

A real-time 512-element photoacoustic tomography system for small animal imaging using a ring ultrasound array has been developed. The system, based upon a 5 MHz transducer array formed along a 50 mm circular aperture, achieves sub-200 micron lateral resolution over a 2 cm disk-shaped region. Corresponding elevation resolutions of 0.6 to 2.5 mm over the central volume enable depth-resolved 3D tomographic imaging with linear translation. Using 8:1 electronic multiplexing, imaging at up to 8 frame/sec is demonstrated for both dynamic phantoms and in vivo mouse and brain samples. The real-time, full 2D tomographic capability of the system paves the way for functional photoacoustic tomographic imaging studies in small animals with sub-second time frame.

In vivo carbon nanotube-enhanced non-invasive photoacoustic mapping of the sentinel lymph node

Sentinel lymph node biopsy (SLNB), a less invasive alternative to axillary lymph node dissection (ALND), has become the standard of care for patients with clinically node-negative breast cancer. In SLNB, lymphatic mapping with radio-labeled sulfur colloid and/or blue dye helps identify the sentinel lymph node (SLN), which is most likely to contain metastatic breast cancer. Even though SLNB, using both methylene blue and radioactive tracers, has a high identification rate, it still relies on an invasive surgical procedure, with associated morbidity. In this study, we have demonstrated a non-invasive single-walled carbon nanotube (SWNT)-enhanced photoacoustic (PA) identification of SLN in a rat model. We have successfully imaged the SLN in vivo by PA imaging (793 nm laser source, 5 MHz ultrasonic detector) with high contrast-to-noise ratio (=89) and good resolution ( approximately 500 microm). The SWNTs also show a wideband optical absorption, generating PA signals over an excitation wavelength range of 740-820 nm. Thus, by varying the incident light wavelength to the near infrared region, where biological tissues (hemoglobin, tissue pigments, lipids and water) show low light absorption, the imaging depth is maximized. In the future, functionalization of the SWNTs with targeting groups should allow the molecular imaging of breast cancer.

Photoacoustic tomography imaging system based on digital B-mode ultrasound diagnosis equipment

A photoacoustic tomography imaging system with B-mode digital ultrasound diagnosis equipment was designed, fabricated, and tested for biomedical imaging. The system for biomedical diagnostics was designed with four parts consisting of an irradiating source, digital B-mode ultrasound diagnostic equipment with a liner transducer array, data acquisition and transmission system, and rotation devices for imaging scanning. Multielement phase-controlled focus technique was applied to reduce the data-acquisition time and enhance signal-to-noise ratio. Under the self-designed working sequence, data-acquisition time of an image can be reduced to 12 s. A series of experiments were performed to estimate the imaging quality of the system. The lateral and elevational resolutions of the system were detected to be 0.25 and 1.4 mm, respectively. The results show that the system has the feasibility to be used as a medical imaging method, and it may provide a new rapid and noninvasive imaging modality for clinical applications.

Microwave-induced thermoacoustic scanning CT for high-contrast and noninvasive breast cancer imaging

A fast thermoacoustic computed tomography system with a multielement linear transducer array was developed to image biological tissues with circular scanning. The spatial resolution of the imaging system and the spectra of the thermoacoustic signals were analyzed. A modified integration backprojection algorithm using velocity potential was employed to recover the direct energy deposition distribution, signal processing methods, and reconstruction algorithms were validated by imaging a phantom. The differences of the microwave-frequency dielectric properties between malignant and normal adipose-dominated tissues in the breast are considerable, and the absorption contrast can reach as large as 6:1 at 1.2 GHz. An experiment of human breast tissue with a tumor was performed with this system; the thermoacoustic images reconstructed by a limited-field-filtered backprojection algorithm and a modified integration backprojection algorithm were also compared with a mammogram. Our results show that the system can provide a rapid and noninvasive approach for high-contrast breast cancer imaging.

Study of cerebrovascular reserve capacity by magnetic resonance perfusion weighted imaging and photoacoustic imaging

The aim of this work was to assess the feasibility of photoacoustic imaging (PAI) and MR imaging for evaluating the cerebrovascular reserve capacity (CVRC) in animal models. Wistar-Kyoto (WKY) rats and spontaneous hypertensive rats (SHR) were used for MRI. BALB/c mice were used for PAI. MR perfusion weighted imaging (PWI) was performed on a 1.5-T whole-body MR system before and after oral administration of acetazolamide (ACZ). The region of interest (ROI) was chosen in the bilateral frontal lobe for measuring regional cerebral blood flow (rCBF), regional cerebral blood volume (rCBV) and mean transit time (MTT). The vessel diameters of the superficial layer of the cortex were measured by PAI in the resting and ACZ-activated mice. The results showed that there was a statistical difference between the resting and ACZ-activated animals in vessel diameter, rCBV and rCBF values. The increments in rCBV and rCBF of WKY rats between resting and ACZ test states were significantly higher than that of SHR. The pathological findings of small arterial walls and lumen of the brain were also different between WKY and SHR rats. The diameters of blood vessels in the superficial layer of the brain measured by PAI were enlarged after the ACZ tolerance test. This result was also observed in the MRI CBV map, where the signal of the vessel in the superficial layer of the cortex became redder after the ACZ stimulation, suggesting the increase of blood flow. It can be concluded that MR PWI and PAI combined with the ACZ test might be useful in evaluating the CVRC and revealing the pathologic changes in cerebral vessels.

Curved array photoacoustic tomographic system for small animal imaging

We present systematic characterization of a photoacoustic imaging system optimized for rapid, high-resolution tomographic imaging of small animals. The system is based on a 128-element ultrasonic transducer array with a 5-MHz center frequency and 80% bandwidth shaped to a quarter circle of 25 mm radius. A 16-channel data-acquisition module and dedicated channel detection electronics enable capture of a 90-deg field-of-view image in less than 1 s and a complete 360-deg scan using sample rotation within 15 s. Measurements on cylindrical phantom targets demonstrate a resolution of better than 200 microm and high-sensitivity detection of 580-microm blood tubing to depths greater than 3 cm in a turbid medium with reduced scattering coefficient mu(s) (')=7.8 cm(-1). The system is used to systematically investigate the effects of target size, orientation, and geometry on tomographic imaging. As a demonstration of these effects and the system imaging capabilities, we present tomographic photoacoustic images of the brain vasculature of an ex vivo mouse with varying measurement aperture. For the first time, according to our knowledge, resolution of sub-200-microm vessels with an overlying turbid medium of greater than 2 cm depth is demonstrated using only intrinsic biological contrast.

Photoacoustic Imaging Application in Tumor Diagnosis and Treatment Monitoring

Photoacoustic imaging (also called optoacoustic or thermoacoustic imaging) can image vascularity clearly with simultaneous high contrast and high spatial resolution, and has the potential to be an application for tumor diagnosis and treatment monitoring. In a unique photoacoustic system, a single pulse laser beam was used as the light source for both cancer treatment and for concurrently generating ultrasound signals for photoacoustic imaging. The photoacoustic system was used to detect early tumor on the rat back, and the vascular structure around the tumor could be imaged clearly with optimal contrast. This system was also used to monitoring damage of the vascular structures before, during and after photodynamic therapy of tumor. This work demonstrates that photoacoustic imaging can potentially be used to guide photodynamic therapy and other phototherapies using vascular changes during treatment. Prospective application of photoacoustic imaging is to characterize and monitor the accumulation of gold nanoshells in vivo to guide nanoshell-based thermal tumor therapy.

High antinoise photoacoustic tomography based on a modified filtered backprojection algorithm with combination wavelet

How to extract the weak photoacoustic signals from the collected signals with high noise is the key to photoacoustic signal processing. We have developed a modified filtered backprojection algorithm based on combination wavelet for high antinoise photoacoustic tomography. A Q-switched-Nd: yttrium-aluminum-garnet laser operating at 532 nm is used as light source. The laser has a pulse width of 7 ns and a repetition frequency of 20 Hz. A needle polyvinylidene fluoride hydrophone with diameter of 1 mm is used to capture photoacoustic signals. The modified algorithm is successfully applied to imaging vascular network of a chick embryo chorioallantoic membrane in situ and brain structure of a rat brain in vivo, respectively. In the reconstructed images, almost all of the capillary vessels and the vascular ramifications of the chick embryo chorioallantoic membrane are accurately resolved, and the detailed brain structures of the rat brain organization are clearly identified with the skull and scalp intact. The experimental results demonstrate that the modified algorithm has much higher antinoise capacity, and can greatly improve the reconstruction image quality. The spatial resolution of the reconstructed images can reach 204 microm. The modified filtered back-projection algorithm based on the combination wavelet has the potential in the practical high-noise signal processing for deeply penetrating photoacoustic tomography.