A real-time photoacoustic tomography system for small animals

Enabling high-speed wide-field dynamic imaging in multifocal photoacoustic computed microscopy: a simulation study

Photoacoustic-computed microscopy (PACM) is an emerging technology that employs thousands of optical foci to provide wide-field high-resolution images of tissue optical absorption. A major limitation of PACM is the slow imaging speed, limiting its usage in dynamic imaging. In this study, we improved the speed through a two-step approach. First, we employed compressed sensing with partially known support to reduce the transducer element number, which subsequently improved the imaging speed at each optical scanning step. Second, we use the high-speed low-resolution image acquired without microlens array to inform dynamic changes in the high-resolution PACM image. Combining both approaches, we achieved high-resolution dynamic imaging over a wide field.

Wearable 3-D Photoacoustic Tomography for Functional Brain Imaging in Behaving Rats

Understanding the relationship between brain function and behavior remains a major challenge in neuroscience. Photoacoustic tomography (PAT) is an emerging technique that allows for noninvasive in vivo brain imaging at micrometer-millisecond spatiotemporal resolution. In this article, a novel, miniaturized 3D wearable PAT (3D-wPAT) technique is described for brain imaging in behaving rats. 3D-wPAT has three layers of fully functional acoustic transducer arrays. Phantom imaging experiments revealed that the in-plane X-Y spatial resolutions were ~200 μm for each acoustic detection layer. The functional imaging capacity of 3D-wPAT was demonstrated by mapping the cerebral oxygen saturation via multi-wavelength irradiation in behaving hyperoxic rats. In addition, we demonstrated that 3D-wPAT could be used for monitoring sensory stimulus-evoked responses in behaving rats by measuring hemodynamic responses in the primary visual cortex during visual stimulation. Together, these results show the potential of 3D-wPAT for brain study in behaving rodents.

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.

Bond-selective photoacoustic imaging by converting molecular vibration into acoustic waves

The quantized vibration of chemical bonds provides a way of detecting specific molecules in a complex tissue environment. Unlike pure optical methods, for which imaging depth is limited to a few hundred micrometers by significant optical scattering, photoacoustic detection of vibrational absorption breaks through the optical diffusion limit by taking advantage of diffused photons and weak acoustic scattering. Key features of this method include both high scalability of imaging depth from a few millimeters to a few centimeters and chemical bond selectivity as a novel contrast mechanism for photoacoustic imaging. Its biomedical applications spans detection of white matter loss and regeneration, assessment of breast tumor margins, and diagnosis of vulnerable atherosclerotic plaques. This review provides an overview of the recent advances made in vibration-based photoacoustic imaging and various biomedical applications enabled by this new technology.

High-Throughput Sparsity-Based Inversion Scheme for Optoacoustic Tomography

The concept of sparsity is extensively exploited in the fields of data acquisition and image processing, contributing to better signal-to-noise and spatio-temporal performance of the various imaging methods. In the field of optoacoustic tomography, the image reconstruction problem is often characterized by computationally extensive inversion of very large datasets, for instance when acquiring volumetric multispectral data with high temporal resolution. In this article we seek to accelerate accurate model-based optoacoustic inversions by identifying various sources of sparsity in the forward and inverse models as well as in the single- and multi-frame representation of the projection data. These sources of sparsity are revealed through appropriate transformations in the signal, model and image domains and are subsequently exploited for expediting image reconstruction. The sparsity-based inversion scheme was tested with experimental data, offering reconstruction speed enhancement by a factor of 40 to 700 times as compared with the conventional iterative model-based inversions while preserving similar image quality. The demonstrated results pave the way for achieving real-time performance of model-based reconstruction in multi-dimensional optoacoustic imaging.

Micromachined Silicon Parallel Acoustic Delay Lines as Time Delayed Ultrasound Detector Array for Real-Time Photoacoustic Tomography

This paper reports the development of a new 16-channel parallel acoustic delay line (PADL) array for real-time photoacoustic tomography (PAT). The PADLs were directly fabricated from single-crystalline silicon substrates using deep reactive ion etching. Compared with other acoustic delay lines (e.g., optical fibers), the micromachined silicon PADLs offer higher acoustic transmission efficiency, smaller form factor, easier assembly, and mass production capability. To demonstrate its real-time photoacoustic imaging capability, the silicon PADL array was interfaced with one single-element ultrasonic transducer followed by one channel of DAQ electronics to receive 16 channels of photoacoustic signals simultaneously. A PAT image of an optically-absorbing target embedded in an optically-scattering phantom was reconstructed, which matched well with the actual size of the imaged target. Because the silicon PADL array allows a signal-to-channel reduction ratio of 16:1, it could significantly simplify the design and construction of ultrasonic receivers for real-time PAT.

Review on photoacoustic imaging of the brain using nanoprobes

Photoacoustic (PA) tomography (PAT) is a hybrid imaging modality that integrates rich optical contrasts with a high-ultrasonic spatial resolution in deep tissue. Among various imaging applications, PA neuroimaging is becoming increasingly important as it nicely complements the limitations of conventional neuroimaging modalities, such as the low-temporal resolution in magnetic resonance imaging and the low depth-to-resolution ratio in optical microscopy/tomography. In addition, the intrinsic hemoglobin contrast PA neuroimaging has also been greatly improved by recent developments in nanoparticles (NPs). For instance, near-infrared absorbing NPs greatly enhanced the vascular contrast in deep-brain PAT; tumor-targeting NPs allowed highly sensitive and highly specific delineation of brain tumors; and multifunctional NPs enabled comprehensive examination of the brain through multimodal imaging. We aim to give an overview of NPs used in PA neuroimaging. Classifications of various NPs used in PAT will be introduced at the beginning, followed by an overview of PA neuroimaging systems, and finally we will discuss major applications of NPs in PA neuroimaging and highlight representative studies.

Effects of the murine skull in optoacoustic brain microscopy

Despite the great promise behind the recent introduction of optoacoustic technology into the arsenal of small-animal neuroimaging methods, a variety of acoustic and light-related effects introduced by adult murine skull severely compromise the performance of optoacoustics in transcranial imaging. As a result, high-resolution noninvasive optoacoustic microscopy studies are still limited to a thin layer of pial microvasculature, which can be effectively resolved by tight focusing of the excitation light. We examined a range of distortions introduced by an adult murine skull in transcranial optoacoustic imaging under both acoustically- and optically-determined resolution scenarios. It is shown that strong low-pass filtering characteristics of the skull may significantly deteriorate the achievable spatial resolution in deep brain imaging where no light focusing is possible. While only brain vasculature with a diameter larger than 60 µm was effectively resolved via transcranial measurements with acoustic resolution, significant improvements are seen through cranial windows and thinned skull experiments.

Ultrasound in Radiology: From Anatomic, Functional, Molecular Imaging to Drug Delivery and Image-Guided Therapy

During the past decade, ultrasound has expanded medical imaging well beyond the "traditional" radiology setting: a combination of portability, low cost, and ease of use makes ultrasound imaging an indispensable tool for radiologists as well as for other medical professionals who need to obtain imaging diagnosis or guide a therapeutic intervention quickly and efficiently. Ultrasound combines excellent ability for deep penetration into soft tissues with very good spatial resolution, with only a few exceptions (ie, those involving overlying bone or gas). Real-time imaging (up to hundreds and thousands of frames per second) enables guidance of therapeutic procedures and biopsies; characterization of the mechanical properties of the tissues greatly aids with the accuracy of the procedures. The ability of ultrasound to deposit energy locally brings about the potential for localized intervention encompassing the following: tissue ablation, enhancing penetration through the natural barriers to drug delivery in the body and triggering drug release from carrier microparticles and nanoparticles. The use of microbubble contrast agents brings the ability to monitor and quantify tissue perfusion, and microbubble targeting with ligand-decorated microbubbles brings the ability to obtain molecular biomarker information, that is, ultrasound molecular imaging. Overall, ultrasound has become the most widely used imaging modality in modern medicine; it will continue to grow and expand.

Noninvasive high-speed photoacoustic tomography of cerebral hemodynamics in awake-moving rats

We present a noninvasive method of photoacoustic tomography (PAT) for imaging cerebral hemodynamics in awake-moving rats. The wearable PAT (wPAT) system has a size of 15 mm in height and 33 mm in diameter, and a weight of ~8 g (excluding cabling). The wPAT achieved an imaging rate of 3.33 frames/s with a lateral resolution of 243 μm. Animal experiments were designed to show wPAT feasibility for imaging cerebral hemodynamics on awake-moving animals. Results showed that the cerebral oxy-hemoglobin and deoxy-hemoglobin changed significantly in response to hyperoxia; and, after the injection of pentylenetetrazol (PTZ), cerebral blood volume changed faster over time and larger in amplitude for rats in awake-moving state compared with rats under anesthesia. By providing a light-weight, high-resolution technology for in vivo monitoring of cerebral hemodynamics in awake-behaving animals, it will be possible to develop a comprehensive understanding on how activity alters hemodynamics in normal and diseased states.

Imaging of small animal peripheral artery disease models: recent advancements and translational potential

Peripheral artery disease (PAD) is a broad disorder encompassing multiple forms of arterial disease outside of the heart. As such, PAD development is a multifactorial process with a variety of manifestations. For example, aneurysms are pathological expansions of an artery that can lead to rupture, while ischemic atherosclerosis reduces blood flow, increasing the risk of claudication, poor wound healing, limb amputation, and stroke. Current PAD treatment is often ineffective or associated with serious risks, largely because these disorders are commonly undiagnosed or misdiagnosed. Active areas of research are focused on detecting and characterizing deleterious arterial changes at early stages using non-invasive imaging strategies, such as ultrasound, as well as emerging technologies like photoacoustic imaging. Earlier disease detection and characterization could improve interventional strategies, leading to better prognosis in PAD patients. While rodents are being used to investigate PAD pathophysiology, imaging of these animal models has been underutilized. This review focuses on structural and molecular information and disease progression revealed by recent imaging efforts of aortic, cerebral, and peripheral vascular disease models in mice, rats, and rabbits. Effective translation to humans involves better understanding of underlying PAD pathophysiology to develop novel therapeutics and apply non-invasive imaging techniques in the clinic.

Imaging Approaches in Functional Assessment of Implantable Myogenic Biomaterials and Engineered Muscle Tissue

The fields of tissue engineering and regenerative medicine utilize implantable biomaterials and engineered tissues to regenerate damaged cells or replace lost tissues. There are distinct challenges in all facets of this research, but functional assessments and monitoring of such complex environments as muscle tissues present the current strategic priority. Many extant methods for addressing these questions result in the destruction or alteration of tissues or cell populations under investigation. Modern advances in non-invasive imaging modalities present opportunities to rethink some of the anachronistic methods, however, their standard employment may not be optimal when considering advancements in myology. New image analysis protocols and/or combinations of established modalities need to be addressed. This review focuses on efficacies and limitations of available imaging modalities to the functional assessment of implantable myogenic biomaterials and engineered muscle tissues.

Optical drug monitoring: photoacoustic imaging of nanosensors to monitor therapeutic lithium in vivo

Personalized medicine could revolutionize how primary care physicians treat chronic disease and how researchers study fundamental biological questions. To realize this goal, we need to develop more robust, modular tools and imaging approaches for in vivo monitoring of analytes. In this report, we demonstrate that synthetic nanosensors can measure physiologic parameters with photoacoustic contrast, and we apply that platform to continuously track lithium levels in vivo. Photoacoustic imaging achieves imaging depths that are unattainable with fluorescence or multiphoton microscopy. We validated the photoacoustic results that illustrate the superior imaging depth and quality of photoacoustic imaging with optical measurements. This powerful combination of techniques will unlock the ability to measure analyte changes in deep tissue and will open up photoacoustic imaging as a diagnostic tool for continuous physiological tracking of a wide range of analytes.

Tripling the detection view of high-frequency linear-array-based photoacoustic computed tomography by using two planar acoustic reflectors

BACKGROUND

Linear-array-based photoacoustic computed tomography (PACT) suffers from a limited view. Circular scanning does increase the detection view angle but is time-consuming. Therefore, it is desirable to increase the detection view angle of linear-array-based PACT without sacrificing imaging speed.

METHODS

Two planar acoustic reflectors placed at 120 degrees to each other were added to a linear-array-based PACT system. Each reflector redirects originally undetectable photoacoustic waves back to the transducer array elements, and together they triple the original detection view angle of the PACT system.

RESULTS

Adding two reflectors increased the detection view angle from 80 to 240 degrees. As a comparison, a single-reflector PACT has a detection view angle of only 160 degrees. A leaf skeleton phantom with a rich vascular network was imaged with the double-reflector PACT, and most of its features were recovered.

CONCLUSIONS

The two acoustic reflectors triple the detection view angle of a linear-array-based PACT without compromising the original imaging speed. This nearly full-view detection capability produces higher-quality images than single-reflector PACT or conventional PACT without reflectors.

Photoacoustic imaging platforms for multimodal imaging

Photoacoustic (PA) imaging is a hybrid biomedical imaging method that exploits both acoustical Epub ahead of print and optical properties and can provide both functional and structural information. Therefore, PA imaging can complement other imaging methods, such as ultrasound imaging, fluorescence imaging, optical coherence tomography, and multi-photon microscopy. This article reviews techniques that integrate PA with the above imaging methods and describes their applications.

In vivo deep brain imaging of rats using oral-cavity illuminated photoacoustic computed tomography

Using internal illumination with an optical fiber in the oral cavity, we demonstrate, for the first time, photoacoustic computed tomography (PACT) of the deep brain of rats in vivo. The experiment was performed on a full-ring-array PACT system, with the capability of providing high-speed cross-sectional imaging of the brain. Compared with external illumination through the cranial skull, internal illumination delivers more light to the base of the brain. Consequently, in vivo photoacoustic images clearly reveal deep brain structures such as the hypothalamus, brain stem, and cerebral medulla.

Expediting model-based optoacoustic reconstructions with tomographic symmetries

PURPOSE

Image quantification in optoacoustic tomography implies the use of accurate forward models of excitation, propagation, and detection of optoacoustic signals while inversions with high spatial resolution usually involve very large matrices, leading to unreasonably long computation times. The development of fast and memory efficient model-based approaches represents then an important challenge to advance on the quantitative and dynamic imaging capabilities of tomographic optoacoustic imaging.

METHODS

Herein, a method for simplification and acceleration of model-based inversions, relying on inherent symmetries present in common tomographic acquisition geometries, has been introduced. The method is showcased for the case of cylindrical symmetries by using polar image discretization of the time-domain optoacoustic forward model combined with efficient storage and inversion strategies.

RESULTS

The suggested methodology is shown to render fast and accurate model-based inversions in both numerical simulations and post mortem small animal experiments. In case of a full-view detection scheme, the memory requirements are reduced by one order of magnitude while high-resolution reconstructions are achieved at video rate.

CONCLUSIONS

By considering the rotational symmetry present in many tomographic optoacoustic imaging systems, the proposed methodology allows exploiting the advantages of model-based algorithms with feasible computational requirements and fast reconstruction times, so that its convenience and general applicability in optoacoustic imaging systems with tomographic symmetries is anticipated.

Detection of a novel mechanism of acousto-optic modulation of incoherent light

A novel form of acoustic modulation of light from an incoherent source has been detected in water as well as in turbid media. We demonstrate that patterns of modulated light intensity appear to propagate as the optical shadow of the density variations caused by ultrasound within an illuminated ultrasonic focal zone. This pattern differs from previous reports of acousto-optical interactions that produce diffraction effects that rely on phase shifts and changes in light directions caused by the acoustic modulation. Moreover, previous studies of acousto-optic interactions have mainly reported the effects of sound on coherent light sources via photon tagging, and/or the production of diffraction phenomena from phase effects that give rise to discrete sidebands. We aimed to assess whether the effects of ultrasound modulation of the intensity of light from an incoherent light source could be detected directly, and how the acoustically modulated (AOM) light signal depended on experimental parameters. Our observations suggest that ultrasound at moderate intensities can induce sufficiently large density variations within a uniform medium to cause measurable modulation of the intensity of an incoherent light source by absorption. Light passing through a region of high intensity ultrasound then produces a pattern that is the projection of the density variations within the region of their interaction. The patterns exhibit distinct maxima and minima that are observed at locations much different from those predicted by Raman-Nath, Bragg, or other diffraction theory. The observed patterns scaled appropriately with the geometrical magnification and sound wavelength. We conclude that these observed patterns are simple projections of the ultrasound induced density changes which cause spatial and temporal variations of the optical absorption within the illuminated sound field. These effects potentially provide a novel method for visualizing sound fields and may assist the interpretation of other hybrid imaging methods.

Handheld photoacoustic tomography probe built using optical-fiber parallel acoustic delay lines

The development of the first miniaturized parallel acoustic delay line (PADL) probe for handheld photoacoustic tomography (PAT) is reported. Using fused-silica optical fibers with low acoustic attenuation, we constructed two arrays of eight PADLs. Precision laser micromachining was conducted to produce robust and accurate mechanical support and alignment structures for the PADLs, with minimal acoustic distortion and interchannel coupling. The 16 optical-fiber PADLs, each with a different time delay, were arranged to form one input port and two output ports. A handheld PADL probe was constructed using two single-element transducers and two data acquisition channels (equal to a channel reduction ratio of 8∶1). Photoacoustic (PA) images of a black-ink target embedded in an optically scattering phantom were successfully acquired. After traveling through the PADLs, the eight channels of differently time-delayed PA signals reached each single-element ultrasonic transducer in a designated nonoverlapping time series, allowing clear signal separation for PA image reconstruction. Our results show that the PADL technique and the handheld probe can potentially enable real-time PAT, while significantly reducing the complexity and cost of the ultrasound receiver system.

Photoacoustic radiation force on a microbubble

We investigate the radiation force on a microbubble due to the photoacoustic wave which is generated by using a pulsed laser. In particular, we focus on the dependence of pulsed laser parameters on the radiation force. In order to do so, we first obtain a new and comprehensive analytical solution to the photoacoustic wave equation based on the Fourier transform for various absorption profiles. Then, we write an expression of the radiation force containing explicit laser parameters, pulse duration, and beamwidth of the laser. Furthermore, we calculate the primary radiation force acting on a microbubble. We show that laser parameters and the position of the microbubble relative to a photoacoustic source have a considerable effect on the primary radiation force. By means of recent developments in laser technologies that render tunability of pulse duration and repetition frequency possible, an adjustable radiation force can be applied to microbubbles. High spatial control of applied force is ensured on account of smaller focal spots achievable by focused optics. In this context, conventional piezoelectric acoustic source applications could be surpassed. In addition, it is possible to increase the radiation force by making source wavelength with the absorption peak of absorber concurrent. The application of photoacoustic radiation force can open a cache of opportunities such as manipulation of microbubbles used as contrast agents and as carrier vehicles for drugs and genes with a desired force along with in vivo applications.

Photoacoustic microscopy and computed tomography: from bench to bedside

Photoacoustic imaging (PAI) of biological tissue has seen immense growth in the past decade, providing unprecedented spatial resolution and functional information at depths in the optical diffusive regime. PAI uniquely combines the advantages of optical excitation and those of acoustic detection. The hybrid imaging modality features high sensitivity to optical absorption and wide scalability of spatial resolution with the desired imaging depth. Here we first summarize the fundamental principles underpinning the technology, then highlight its practical implementation, and finally discuss recent advances toward clinical translation.

Fast spatiotemporal image reconstruction based on low-rank matrix estimation for dynamic photoacoustic computed tomography

In order to monitor dynamic physiological events in near-real time, a variety of photoacoustic computed tomography (PACT) systems have been developed that can rapidly acquire data. Previously reported studies of dynamic PACT have employed conventional static methods to reconstruct a temporally ordered sequence of images on a frame-by-frame basis. Frame-by-frame image reconstruction (FBFIR) methods fail to exploit correlations between data frames and are known to be statistically and computationally suboptimal. In this study, a low-rank matrix estimation-based spatiotemporal image reconstruction (LRME-STIR) method is investigated for dynamic PACT applications. The LRME-STIR method is based on the observation that, in many PACT applications, the number of frames is much greater than the rank of the ideal noiseless data matrix. Using both computer-simulated and experimentally measured photoacoustic data, the performance of the LRME-STIR method is compared with that of conventional FBFIR method followed by image-domain filtering. The results demonstrate that the LRME-STIR method is not only computationally more efficient but also produces more accurate dynamic PACT images than a conventional FBFIR method followed by image-domain filtering.

Small-animal whole-body photoacoustic tomography: a review

With the wide use of small animals for biomedical studies, in vivo small-animal whole-body imaging plays an increasingly important role. Photoacoustic tomography (PAT) is an emerging whole-body imaging modality that shows great potential for preclinical research. As a hybrid technique, PAT is based on the acoustic detection of optical absorption from either endogenous tissue chromophores, such as oxyhemoglobin and deoxyhemoglobin, or exogenous contrast agents. Because ultrasound scatters much less than light in tissue, PAT generates high-resolution images in both the optical ballistic and diffusive regimes. Using near-infrared light, which has relatively low blood absorption, PAT can image through the whole body of small animals with acoustically defined spatial resolution. Anatomical and vascular structures are imaged with endogenous hemoglobin contrast, while functional and molecular images are enabled by the wide choice of exogenous optical contrasts. This paper reviews the rapidly growing field of small-animal whole-body PAT and highlights studies done in the past decade.

Volumetric optoacoustic imaging with multi-bandwidth deconvolution

Optoacoustic (photoacoustic) imaging based on cylindrically focused 1-D transducer arrays comes with powerful characteristics in visualizing optical contrast. Parallel reading of multiple detectors arranged around a tissue cross section enables capturing data for generating images of this plane within micro-seconds. Dedicated small animals scanners and handheld systems using 1-D cylindrically focused ultrasound transducer arrays have demonstrated real-time cross-sectional imaging and high in-plane resolution. Yet, the resolution achieved along the axis perpendicular to the focal plane, i.e., the elevation resolution, is determined by the focusing capacities of the detector and is typically lower than the in-plane resolution. Herein, we investigated whether deconvolution of the sensitivity field of the transducer could lead to tangible image improvements. We showcase the findings on experimental measurements from phantoms and animals and discuss the features and the limitations of the approach in improving resolution along the elevation dimension.

Multicontrast photoacoustic in vivo imaging using near-infrared fluorescent proteins

Non-invasive imaging of biological processes in vivo is invaluable in advancing biology. Photoacoustic tomography is a scalable imaging technique that provides higher resolution at greater depths in tissue than achievable by purely optical methods. Here we report the application of two spectrally distinct near-infrared fluorescent proteins, iRFP670 and iRFP720, engineered from bacterial phytochromes, as photoacoustic contrast agents. iRFPs provide tissue-specific contrast without the need for delivery of any additional substances. Compared to conventional GFP-like red-shifted fluorescent proteins, iRFP670 and iRFP720 demonstrate stronger photoacoustic signals at longer wavelengths, and can be spectrally resolved from each other and hemoglobin. We simultaneously visualized two differently labeled tumors, one with iRFP670 and the other with iRFP720, as well as blood vessels. We acquired images of a mouse as 2D sections of a whole animal, and as localized 3D volumetric images with high contrast and sub-millimeter resolution at depths up to 8 mm. Our results suggest iRFPs are genetically-encoded probes of choice for simultaneous photoacoustic imaging of several tissues or processes in vivo.

Artifact removal in photoacoustic section imaging by combining an integrating cylindrical detector with model-based reconstruction

Photoacoustic section imaging reveals optically absorbing structures within a thin slice of an object. It requires measuring acoustic waves excited by absorption of short laser pulses with a cylindrical acoustic lens detector rotating around the object. Owing to the finite detector size and its limited depth of focus, various artifacts arise, seen as distortions within the imaging slice and cross-talk from neighboring areas of the object. The presented solution aims at avoiding these artifacts by a special design of the sensor and by use of a model-based reconstruction algorithm that improves section images by incorporating information from neighboring sections. The integrating property of the cylindrical detector, which exceeds in direction of the cylinder axis the size of the imaged object, avoids the lateral blurring that normally results from the finite width of a small detector. Applying a maximum likelihood reconstruction method for the inversion of the imaging system matrix to the temporal pressure signals yields line projections of the initial energy distribution, from which section images are obtained by applying the inverse Radon transform. By using data from few sections, a significant reduction of artifacts related to the imperfections of the sensor is demonstrated both in simulations and in phantom experiments.

Acoustic Inversion in Optoacoustic Tomography: A Review

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

High-resolution photoacoustic tomography of resting-state functional connectivity in the mouse brain

The increasing use of mouse models for human brain disease studies presents an emerging need for a new functional imaging modality. Using optical excitation and acoustic detection, we developed a functional connectivity photoacoustic tomography system, which allows noninvasive imaging of resting-state functional connectivity in the mouse brain, with a large field of view and a high spatial resolution. Bilateral correlations were observed in eight functional regions, including the olfactory bulb, limbic, parietal, somatosensory, retrosplenial, visual, motor, and temporal regions, as well as in several subregions. The borders and locations of these regions agreed well with the Paxinos mouse brain atlas. By subjecting the mouse to alternating hyperoxic and hypoxic conditions, strong and weak functional connectivities were observed, respectively. In addition to connectivity images, vascular images were simultaneously acquired. These studies show that functional connectivity photoacoustic tomography is a promising, noninvasive technique for functional imaging of the mouse brain.

Retrospective respiration-gated whole-body photoacoustic computed tomography of mice

Photoacoustic tomography (PAT) is an emerging technique that has a great potential for preclinical whole-body imaging. To date, most whole-body PAT systems require multiple laser shots to generate one cross-sectional image, yielding a frame rate of <1  Hz. Because a mouse breathes at up to 3 Hz, without proper gating mechanisms, acquired images are susceptible to motion artifacts. Here, we introduce, for the first time to our knowledge, retrospective respiratory gating for whole-body photoacoustic computed tomography. This new method involves simultaneous capturing of the animal's respiratory waveform during photoacoustic data acquisition. The recorded photoacoustic signals are sorted and clustered according to the respiratory phase, and an image of the animal at each respiratory phase is reconstructed subsequently from the corresponding cluster. The new method was tested in a ring-shaped confocal photoacoustic computed tomography system with a hardware-limited frame rate of 0.625 Hz. After respiratory gating, we observed sharper vascular and anatomical images at different positions of the animal body. The entire breathing cycle can also be visualized at 20  frames/cycle.

Wide-field two-dimensional multifocal optical-resolution photoacoustic-computed microscopy

Optical-resolution photoacoustic microscopy (OR-PAM) is an emerging technique that directly images optical absorption in tissue at high spatial resolution. To date, the majority of OR-PAM systems are based on single-focused optical excitation and ultrasonic detection, limiting the wide-field imaging speed. While 1D multifocal OR-PAM (1D-MFOR-PAM) has been developed, the potential of microlens and transducer arrays has not been fully realized. Here we present the development of 2D multifocal optical-resolution photoacoustic-computed microscopy (2D-MFOR-PACM), using a 2D microlens array and a full-ring ultrasonic transducer array. The 10 mm×10 mm microlens array generates 1800 optical foci within the focal plane of the 512-element transducer array, and raster scanning the microlens array yields optical-resolution photoacoustic images. The system has improved the in-plane resolution of a full-ring transducer array from ≥100 to 29 μm and achieved an imaging time of 36 s over a 10 mm×10 mm field of view. In comparison, the 1D-MFOR-PAM would take more than 4 min to image over the same field of view. The imaging capability of the system was demonstrated on phantoms and animals both ex vivo and in vivo.

Volumetric real-time tracking of peripheral human vasculature with GPU-accelerated three-dimensional optoacoustic tomography

Optoacoustic tomography provides a unique possibility for ultra-high-speed 3-D imaging by acquiring complete volumetric datasets from interrogation of tissue by a single nanosecond-duration laser pulse. Yet, similarly to ultrasound, optoacoustics is a time-resolved imaging method, thus, fast 3-D imaging implies real-time acquisition and processing of high speed data from hundreds of detectors simultaneously, which presents significant technological challenges. Herein we present a highly efficient graphical processing unit (GPU) framework for real-time reconstruction and visualization of 3-D tomographic optoacoustic data. By utilizing a newly developed 3-D optoacoustic scanner, which simultaneously acquires signals with a handheld 256-element spherical ultrasonic array system, we further demonstrate tracking of deep tissue human vasculature rendered at a rate of 10 volumetric frames per second. The flexibility provided by the handheld hardware design, combined with the real-time operation, makes the developed platform highly usable for both clinical imaging practice and small animal research applications.

Deep-tissue photoacoustic tomography of Förster resonance energy transfer

Förster resonance energy transfer (FRET) is a distance-dependent process that transfers excited state energy from a donor molecule to an acceptor molecule without the emission of a photon. The FRET rate is determined by the proximity between the donor and the acceptor molecules; it becomes significant only when the proximity is within several nanometers. Therefore, FRET has been applied to visualize interactions and conformational changes of biomolecules, such as proteins, lipids, and nucleic acids that cannot be resolved by optical microscopy. Here, we report photoacoustic tomography of FRET efficiency at a 1-cm depth in chicken breast tissue, whereas conventional high-resolution fluorescence imaging is limited to <0.1  cm. Photoacoustic tomography is expected to facilitate the examination of FRET phenomena in living organisms.

An optimized ultrasound detector for photoacoustic breast tomography

PURPOSE

Photoacoustic imaging has proven to be able to detect vascularization-driven optical absorption contrast associated with tumors. In order to detect breast tumors located a few centimeter deep in tissue, a sensitive ultrasound detector is of crucial importance for photoacoustic mammography. Further, because the expected photoacoustic frequency bandwidth (a few MHz to tens of kHz) is inversely proportional to the dimensions of light absorbing structures (0.5-10+ mm), proper choices of materials and their geometries and proper considerations in design have to be made to implement optimal photoacoustic detectors. In this study, we design and evaluate a specialized ultrasound detector for photoacoustic mammography.

METHODS

Based on the required detector sensitivity and its frequency response, a selection of active material and matching layers and their geometries is made leading to functional detector models. By iteration between simulation of detector performances, fabrication and experimental characterization of functional models an optimized implementation is made and evaluated. For computer simulation, we use 1D Krimholtz-Leedom-Matthaei and 3D finite-element based models.

RESULTS

The experimental results of the designed first and second functional detectors matched with the simulations. In subsequent bare piezoelectric samples the effect of lateral resonances was addressed and their influence minimized by subdicing the samples. Consequently, using simulations, a final optimized detector was designed, with a center frequency of 1 MHz and a -6 dB bandwidth of 0.4-1.25 MHz (fractional bandwidth of ~80%). The minimum detectable pressure was measured to be 0.5 Pa.

CONCLUSIONS

A single-element, large-aperture, sensitive, and broadband detector is designed and developed for photoacoustic tomography of the breast. The detector should be capable of detecting vascularized tumors with 1-2 mm resolution. The minimum detectable pressure is 0.5 Pa, which will facilitate deeper imaging compared to the current systems. Further improvements by proper electrical grounding and shielding and implementation of this design into an arrayed detector will pave the way for clinical applications of photoacoustic mammography.

Imaging challenges in biomaterials and tissue engineering

Biomaterials are employed in the fields of tissue engineering and regenerative medicine (TERM) in order to enhance the regeneration or replacement of tissue function and/or structure. The unique environments resulting from the presence of biomaterials, cells, and tissues result in distinct challenges in regards to monitoring and assessing the results of these interventions. Imaging technologies for three-dimensional (3D) analysis have been identified as a strategic priority in TERM research. Traditionally, histological and immunohistochemical techniques have been used to evaluate engineered tissues. However, these methods do not allow for an accurate volume assessment, are invasive, and do not provide information on functional status. Imaging techniques are needed that enable non-destructive, longitudinal, quantitative, and three-dimensional analysis of TERM strategies. This review focuses on evaluating the application of available imaging modalities for assessment of biomaterials and tissue in TERM applications. Included is a discussion of limitations of these techniques and identification of areas for further development.

Enhancement of photoacoustic tomography by ultrasonic computed tomography based on optical excitation of elements of a full-ring transducer array

Photoacoustic computed tomography (PACT) is a hybrid technique that combines optical excitation and ultrasonic detection to provide high-resolution images in deep tissues. In the image reconstruction, a constant speed of sound (SOS) is normally assumed. This assumption, however, is often not strictly satisfied in deep tissue imaging, due to acoustic heterogeneities within the object and between the object and the coupling medium. If these heterogeneities are not accounted for, they will cause distortions and artifacts in the reconstructed images. In this Letter, we incorporated ultrasonic computed tomography (USCT), which measures the SOS distribution within the object, into our full-ring array PACT system. Without the need for ultrasonic transmitting electronics, USCT was performed using the same laser beam as for PACT measurement. By scanning the laser beam on the array surface, we can sequentially fire different elements. As a first demonstration of the system, we studied the effect of acoustic heterogeneities on photoacoustic vascular imaging. We verified that constant SOS is a reasonable approximation when the SOS variation is small. When the variation is large, distortion will be observed in the periphery of the object, especially in the tangential direction.

Calibration-free quantification of absolute oxygen saturation based on the dynamics of photoacoustic signals

Photoacoustic tomography (PAT) is a hybrid imaging technique that has broad preclinical and clinical applications. Based on the photoacoustic effect, PAT directly measures specific optical absorption, which is the product of the tissue-intrinsic optical absorption coefficient and the local optical fluence. Therefore, quantitative PAT, such as absolute oxygen saturation (sO₂) quantification, requires knowledge of the local optical fluence, which can only be estimated through invasive measurements or sophisticated modeling of light transportation. In this Letter, we circumvent this requirement by taking advantage of the dynamics in sO₂. The new method works when the sO₂ transition can be simultaneously monitored with multiple wavelengths. For each wavelength, the ratio of photoacoustic amplitudes measured at different sO₂ states is utilized. Using the ratio cancels the contribution from optical fluence and allows calibration-free quantification of absolute sO₂. The new method was validated through both phantom and in vivo experiments.

Photoacoustic tomography extracted from speckle noise in acoustically inhomogeneous tissue

Photoacoustic tomography is usually limited to acoustically homogeneous tissue. A hybrid scheme is developed to break this limitation by utilizing ultrasound to determine the unknown Green's function of inhomogeneous tissue. The method can effectively decrease the distortion and false contrast in images by extracting information from speckle noise. The method does not depend on the prior knowledge of tissue and the medium complexity. Moreover, the estimation of Green's function and the photoacoustic detection are performed by the same transducer. Therefore, the scheme could be easily integrated into a classical photoacoustic tomography system and extend its application in speckle environment.

Spectroscopic Imaging of Deep Tissue through Photoacoustic Detection of Molecular Vibration

The quantized vibration of chemical bonds provides a way of imaging target molecules in a complex tissue environment. Photoacoustic detection of harmonic vibrational transitions provides an approach to visualize tissue content beyond the ballistic photon regime. This method involves pulsed laser excitation of overtone transitions in target molecules inside a tissue. Fast relaxation of the vibrational energy into heat results in a local temperature rise on the order of mK and a subsequent generation of acoustic waves detectable with an ultrasonic transducer. In this perspective, we review recent advances that demonstrate the advantages of vibration-based photoacoustic imaging and illustrate its potential in diagnosing cardiovascular plaques. An outlook into future development of vibrational photoacoustic endoscopy and tomography is provided.

Optoacoustic imaging and tomography: reconstruction approaches and outstanding challenges in image performance and quantification

This paper comprehensively reviews the emerging topic of optoacoustic imaging from the image reconstruction and quantification perspective. Optoacoustic imaging combines highly attractive features, including rich contrast and high versatility in sensing diverse biological targets, excellent spatial resolution not compromised by light scattering, and relatively low cost of implementation. Yet, living objects present a complex target for optoacoustic imaging due to the presence of a highly heterogeneous tissue background in the form of strong spatial variations of scattering and absorption. Extracting quantified information on the actual distribution of tissue chromophores and other biomarkers constitutes therefore a challenging problem. Image quantification is further compromised by some frequently-used approximated inversion formulae. In this review, the currently available optoacoustic image reconstruction and quantification approaches are assessed, including back-projection and model-based inversion algorithms, sparse signal representation, wavelet-based approaches, methods for reduction of acoustic artifacts as well as multi-spectral methods for visualization of tissue bio-markers. Applicability of the different methodologies is further analyzed in the context of real-life performance in small animal and clinical in-vivo imaging scenarios.

Hybrid photoacoustic and ultrasound section imaging with optical ultrasound detection

A setup is proposed that provides perfectly co-registered photoacoustic (PA) and ultrasound (US) section images. Photoacoustic and ultrasound backscatter signals are generated by laser pulses coming from the same laser system, the latter by absorption of some of the laser energy on an optically absorbing target near the imaged object. By measuring both signals with the same optical detector, which is focused into the selected section by use of a cylindrical acoustic mirror, the information for both images is acquired simultaneously. Co-registered PA and US images are obtained after applying the inverse Radon transform to the data, which are gathered while rotating the object relative to the detector. Phantom experiments demonstrate a resolution of 1.1 mm between the sections of both imaging modalities and a in-plane resolution of about 60 µm and 120 µm for the US and PA modes, respectively. The complementary contrast mechanisms of the two modalities are shown by images of a zebrafish.

Optoacoustic imaging of blood perfusion: techniques for intraoperative tissue viability assessment

Reliably assessing tissue viability during surgery is of major importance in surgical procedures. The most basic requirement for viability is sufficient oxygen supply to the tissue. Therefore it is highly desirable to visualize in real-time the dynamic process of blood perfusion up to and within the microvasculature. A modality sensitive to structures in the range of few hundred micrometers and offering high contrast to the embedding tissue is then needed. To this end, a number of methods have been developed, but have had no significant impact on the clinical routine due to various deficiencies. In this paper we demonstrate the applicability of optoacoustic imaging, which combines ultrasonic resolution with strong optical contrast. A method for optoacoustic perfusion assessment, based on a local and repeatable injection of saline, was proposed and assessed ex-vivo on large pig bowels and in-vivo in mouse tails. The obtained dynamic perfusion images highlight the method's potential to enable immediate and quantitative assessment of tissue viability during surgery.

Real-time handheld multispectral optoacoustic imaging

Multispectral optoacoustic tomography (MSOT) of functional and molecular contrast has the potential to find broad deployment in clinical practice. We have developed the first handheld MSOT imaging device with fast wavelength tuning achieving a frame rate of 50 Hz. In this Letter, we demonstrate its clinical potential by dynamically resolving multiple disease-relevant tissue chromophores, including oxy-/deoxyhemoglobin, and melanin, in human volunteers.

A new acoustic lens material for large area detectors in photoacoustic breast tomography

OBJECTIVES

We introduce a new acoustic lens material for photoacoustic tomography (PAT) to improve lateral resolution while possessing excellent acoustic acoustic impedance matching with tissue to minimize lens induced image artifacts.

BACKGROUND

A large surface area detector due to its high sensitivity is preferable to detect weak signals in photoacoustic mammography. The lateral resolution is then limited by the narrow acceptance angle of such detectors. Acoustic lenses made of acrylic plastic (PMMA) have been used to enlarge the acceptance angle of such detectors and improve lateral resolution. However, such PMMA lenses introduce image artifacts due to internal reflections of ultrasound within the lenses, the result of acoustic impedance mismatch with the coupling medium or tissue.

METHODS

A new lens is proposed based on the 2-component resin Stycast 1090SI. We characterized the acoustic properties of the proposed lens material in comparison with commonly used PMMA, inspecting the speed of sound, acoustic attenuation and density. We fabricated acoustic lenses based on the new material and PMMA, and studied the effect of the acoustic lenses on detector performance comparing finite element (FEM) simulations and measurements of directional sensitivity, pulse-echo response and frequency response. We further investigated the effect of using the acoustic lenses on the image quality of a photoacoustic breast tomography system using k-Wave simulations and experiments.

RESULTS

Our acoustic characterization shows that Stycast 1090SI has tissue-like acoustic impedance, high speed of sound and low acoustic attenuation. These acoustic properties ensure an excellent acoustic lens material to minimize the acoustic insertion loss. Both acoustic lenses show significant enlargement of detector acceptance angle and lateral resolution improvement from modeling and experiments. However, the image artifacts induced by the presence of an acoustic lens are reduced using the proposed lens compared to PMMA lens, due to the minimization of internal reflections.

CONCLUSIONS

The proposed Stycast 1090SI acoustic lens improves the lateral resolution of photoacoustic tomography systems while not suffering from internal reflection-induced image artifacts compared a lens made of PMMA.

Speed-of-sound compensated photoacoustic tomography for accurate imaging

PURPOSE

In most photoacoustic (PA) tomographic reconstructions, variations in speed-of-sound (SOS) of the subject are neglected under the assumption of acoustic homogeneity. Biological tissue with spatially heterogeneous SOS cannot be accurately reconstructed under this assumption. The authors present experimental and image reconstruction methods with which 2D SOS distributions can be accurately acquired and reconstructed, and with which the SOS map can be used subsequently to reconstruct highly accurate PA tomograms.

METHODS

The authors begin with a 2D iterative reconstruction approach in an ultrasound transmission tomography setting, which uses ray refracted paths instead of straight ray paths to recover accurate SOS images of the subject. Subsequently, they use the SOS distribution in a new 2D iterative PA reconstruction approach, where refraction of rays originating from PA sources is accounted for in accurately retrieving the distribution of these sources. Both the SOS reconstruction and SOS-compensated PA reconstruction methods utilize the Eikonal equation to model acoustic wavefront propagation. The equation is solved using a high accuracy fast marching method.

RESULTS

The authors validated the new reconstruction algorithms using numerical phantoms. For experiments they utilized the recently introduced PER-PACT method which can be used to simultaneously acquire SOS and PA data from subjects.

CONCLUSIONS

It is first confirmed that it is important to take SOS inhomogeneities into account in high resolution PA tomography. The iterative reconstruction algorithms, that model acoustic refractive effects, in reconstructing SOS distributions, and subsequently using these distributions to correct PA tomograms, yield artifact-free highly accurate images. The approach of using the hybrid measurement method and the new reconstruction algorithms is successful in substantially improving the quality of PA images with a minimization of blurring and artifacts.

Noninvasive photoacoustic detecting intraocular foreign bodies with an annular transducer array

We present a fast photoacoustic imaging system based on an annular transducer array for detection of intraocular foreign bodies. An eight-channel data acquisition system is applied to capture the photoacoustic signals using multiplexing and the total time of data acquisition and transferring is within 3 s. A limited-view filtered back projection algorithm is used to reconstruct the photoacoustic images. Experimental models of intraocular metal and glass foreign bodies were constructed on ex vivo pig's eyes and clear photoacoustic images of intraocular foreign bodies were obtained. Experimental results demonstrate the photoacoustic imaging system holds the potential for in clinic detecting the intraocular foreign bodies.

Micromachined silicon acoustic delay lines for ultrasound applications

In this paper, we present the design, fabrication and testing of novel micromachined silicon-based acoustic delay lines. The acoustic properties of different silicon delay line structures have been characterized. Based on the experiment results, two different acoustic delay line systems (parallel and serial) have been successfully demonstrated to create controlled time delays in multiple channels of ultrasound signals. The time-delayed ultrasound signals are received with a single-element ultrasound transducer in a time-serial manner. This unique capability could be used to merge multiple signal channels, thereby enabling new ultrasound receiver designs with potentially less complexity and lower cost.

Parallel acoustic delay lines for photoacoustic tomography

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

Video-rate functional photoacoustic microscopy at depths

We report the development of functional photoacoustic microscopy capable of video-rate high-resolution in vivo imaging in deep tissue. A lightweight photoacoustic probe is made of a single-element broadband ultrasound transducer, a compact photoacoustic beam combiner, and a bright-field light delivery system. Focused broadband ultrasound detection provides a 44-μm lateral resolution and a 28-μm axial resolution based on the envelope (a 15-μm axial resolution based on the raw RF signal). Due to the efficient bright-field light delivery, the system can image as deep as 4.8 mm in vivo using low excitation pulse energy (28 μJ per pulse, 0.35  mJ/cm² on the skin surface). The photoacoustic probe is mounted on a fast-scanning voice-coil scanner to acquire 40 two-dimensional (2-D) B-scan images per second over a 9-mm range. High-resolution anatomical imaging is demonstrated in the mouse ear and brain. Via fast dual-wavelength switching, oxygen dynamics of mouse cardio-vasculature is imaged in realtime as well.

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.

Noninvasive photoacoustic computed tomography of mouse brain metabolism in vivo

We have demonstrated the feasibility of imaging mouse brain metabolism using photoacoustic computed tomography (PACT), a fast, noninvasive and functional imaging modality with optical contrast and acoustic resolution. Brain responses to forepaw stimulations were imaged transdermally and transcranially. 2-NBDG, which diffuses well across the blood-brain-barrier, provided exogenous contrast for photoacoustic imaging of glucose response. Concurrently, hemoglobin provided endogenous contrast for photoacoustic imaging of hemodynamic response. Glucose and hemodynamic responses were quantitatively decoupled by using two-wavelength measurements. We found that glucose uptake and blood perfusion around the somatosensory region of the contralateral hemisphere were both increased by stimulations, indicating elevated neuron activity. While the glucose response area was more homogenous and confined within the somatosensory region, the hemodynamic response area had a clear vascular pattern and spread wider than the somatosensory region. Our results demonstrate that 2-NBDG-enhanced PACT is a promising tool for noninvasive studies of brain metabolism.