1 January 2007 Photoacoustic imaging of the microvasculature with a high-frequency ultrasound array transducer
Author Affiliations +
J. of Biomedical Optics, 12(1), 010501 (2007). doi:10.1117/1.2709850
Visualization of microvascular networks could provide new information about function and disease. We demonstrate the capabilities of a 30-MHz ultrasound array system for photoacoustic microscopy of small (≤300 μm) vessels in a rat. 3D images obtained by translating the array in the elevation direction are compared with photographs of excised skin. The system is shown to have 100-μm lateral resolution, 25-μm axial resolution, and 3-mm imaging depth. To our knowledge this is the first report on photoacoustic microscopy of the microvasculature with a high-frequency array transducer. It is anticipated that the system can be used for studying and diagnosing a number of diseases including cancer, atherosclerosis, dermatological disorders, and peripheral microvascular complications in diabetes.
Zemp, Bitton, Li, Shung, Stoica, and Wang: Photoacoustic imaging of the microvasculature with a high-frequency ultrasound array transducer

The microvasculature plays an important role in tissue function and pathogenesis,1 yet in many diseases the etiology of microvascular dysfunction is poorly understood. There is a significant need for technologies to visualize microvascular growth and remodeling in the pathogenic process. Early detection of abnormal microvascular dynamics and morphology could prove important in effective interventional strategies. Additionally, microvascular imaging technologies may provide a much needed tool for basic research and for monitoring therapeutic efficacy in vivo.

Photoacoustic imaging has the potential to address some of these needs. In photoacoustic imaging, short pulses of laser light are transmitted into the body. Upon light absorption, a rapid thermoelastic expansion launches acoustic waves that are detected and reconstructed to form images representing optical absorption. Since blood has a visible-light absorption coefficient one to two orders of magnitude larger than surrounding tissue, there is significant endogenous contrast available for imaging blood vessels. Recently, our group reported a novel scanning photoacoustic microscopy system with dark-field illumination and a high numerical aperture 50-MHz single-element ultrasound transducer.2, 3 This system produces clear images of microvessels to depths of about 3mm ; however, it requires mechanical scanning and has only 1-mm depth of field.

We present a photoacoustic microscopy system with a 30-MHz ultrasound array transducer. Advantages of this array over a mechanically scanned single element transducer include the ability to use dynamic receive focusing (greatly increasing the depth of field), electronic rather than mechanical B-scanning capability, and parallel data acquisition. Eventually the probe may be handheld and may provide real-time imaging capabilities. Other groups are beginning to use arrays for photoacoustic imaging,4 but this is the first report of high-frequency (above 20MHz ) photoacoustic imaging of the microvasculature. Due to small detection areas of individual elements, there may be some question regarding whether high-frequency arrays are sensitive enough for imaging in vivo. We show feasibility of the high-frequency array system for photoacoustic visualization of microvessels to depths of 3mm in rats in situ. Due to multiply scattered light, this depth is not attainable by present optical microscopy techniques such as confocal and two-photon microscopy, and is less invasive than intravital techniques. It is anticipated that systems similar to the one reported here may be used for studying and diagnosing a number of diseases where microvascular changes impact health, including cancer,1 dermatological disorders, cardiovascular disease, and peripheral microvascular complications in diabetics. High-frequency photoacoustic imaging may complement high-frequency Doppler ultrasound5, 6 to detect very small vessels as well as those with low flow velocities.

Our system, shown in Fig. 1 , consists of an ultrasound array, a laser system, and receiving electronics. The ultrasound transducer is a 30-MHz piezo-composite array of 48 elements7 with 50% pulse-echo 6-dB fractional bandwidth, 100-μm pitch (2 times the center wavelength), 2-mm elevation aperture and 8.2-mm elevational focus, and 19.1-dB insertion loss (41.2-dB uncompensated loss, with 18.37-, 0.6-, and 3.17-dB losses due to diffraction, reflector absorption, and attenuation, respectively). These pulse-echo measurements were tested with a steel reflector placed at the focal position in water. Depending on fabrication technology, comparable frequency single element transducers have been measured to have insertion losses of 20 to 50dB .8 Signals from the array were amplified by a custom rf board with 33-dB fixed gain and an additional 40dB of variable gain. The board multiplexed the 48 channels down to 16. Four 4-to-1 rf multiplexer boards (MAX4141) were then used to enable data acquisition with a four-channel oscilloscope (Tektronix, TDS5034). A PC running Labview controlled multiplexing via the parallel port, while also communicating with the oscilloscope via a general purpose interface bus (GPIB). The maximum data acquisition rate was limited by the 10-Hz pulse-repetition rate of the laser system. Acquired data were processed offline.

Fig. 1

System diagram. GG: ground-glass diffuser; L: lens; MUX: multiplexer.


A Q-switched Nd:YAG laser (Quantel, Brilliant B) was used to pump a dye laser (Continuum ND6000) to obtain 6-ns pulses of wavelength 584nm , an isosbestic point. Laser pulses were coupled into a 0.37NA , 600-μm optical fiber using a fiber coupler consisting of a ground-glass diffuser, then a lens to direct a fraction of the diffused light directly into the core of the fiber. Approximately 0.3mJ of light was delivered through the fiber, which was positioned obliquely at the side of the array, forming an elliptical illumination pattern (6×3mm) with the major axis parallel to the array scanning plane. Approximately 6mJcm2 of light was delivered to the skin surface, approximately 30% of the permissible ANSI visible light exposure limit of 20mJcm2 per pulse.

Channel data were acquired at 250MSs and processed using a sector-scan beam-forming algorithm modified from ultrasound beam forming. To reconstruct an A-scan line along a line of sight defined by the angle θ , we apply dynamic delays τn(R,θ) to the received channel data S(n,t),n=1,2,,48 , where R is the distance from the image point to the beam-forming origin (taken here as the center of the array), then add the delayed data. The delays are given as


where xn is the distance to the array element from the array center, and c is the speed of sound. This expression is a second-order Taylor approximation in xn to the transit time
between a field point and element xn . The first term of Eq. 1 represents steering, while the second represents focusing. The transit time between a point of interest and the array for pulse-echo ultrasound is t=2Rc , while in photoacoustic imaging propagation is one way, and t=Rc . Hence the steering term xnsinθc is identical to pulse-echo ultrasound, while the focusing term xn2cos2θ2c2t is different. In beam forming, we used all 48 elements to leverage the maximum possible aperture. We also found it advantageous to apply a coherence factor (CF) weighting in the beam-forming procedure.9 A coherence factor for an A-scan line with direction θm is given as


where N is the number of elements. This time-varying function is applied multiplicatively to each beam-formed A-scan and tends to emphasize highly in-phase signals while rejecting low coherence side lobes, grating lobes, and noise.

To demonstrate the spatial resolution and dynamic focusing capability of the system, we imaged 6-μm-diam carbon fibers in water at various spatial locations to produce a composite image, shown in Fig. 2a . Importantly, the point-spread functions are well focused throughout the region of interest. Figure 2b shows the projection of a single point-spread function on the x axis. The 6-dB axial and lateral spatial resolution for points near the elevation focus were measured as 25 and 100±5μm , respectively. Clutter and noise levels (taken as the mean level of the envelope-detected signal relative to the peak) were approximately 33dB when no coherence-factor weighting was used in beam forming. When coherence-factor weighting was applied, mean clutter and noise levels reduced to 60dB .

Fig. 2

(a) Composite B-scan of 6μm carbon fibers in water. (b) B-scan of a single 6-μm carbon fiber in water, projected on the x axis. (c) Photoacoustic B-scan of subcutaneous microvessels in a young rat. In (a) and (c) the z axis is depth from the transducer surface in millimeters. (d) Surface-rendered photoacoustic image of volumetric data acquired by translating the array in the y direction. The z axis is depth from the skin surface in millimeters, with the skin surface at the bottom of the plot. (e) Photograph of the underside of excised skin in the area of photoacoustic imaging.


To demonstrate the ability to see microvessels in animals, we imaged a young Sprague Dawley rat approximately 100g in weight. The hair was removed from the back right side of the rat while the rat was anesthetized, and the animal was sacrificed prior to imaging to reduce motion artifacts. All animal procedures were performed in accordance with approved protocols.

The array and optical fibers were mounted on a 3-axis translation stage, and lowered into a water tank with an acoustically and optically transparent aperture for the animal below the tank. Figure 2c shows a photoacoustic B-scan image at one imaging location with several small vessels visible at 2- to 3-mm depths below the surface (the skin is at approximately the 6-mm level). With 40 times averaging, the mean clutter and noise in this image is 61-dB below the peak vessel signal. A sequence of 51 B-scans was acquired at 0.254-mm elevational intervals ( 0.01 ) to form a volumetric dataset. Presently, elevational scanning was accomplished by manually translating the array on a linear stage. The elevation resolution defined here as the 6-dB one-way receive beam waist was measured as 260±10μm . Figure 2d shows a 3-D surface rendering of the photoacoustic data compared with a photograph of the underside of the excised skin [Fig. 2e]. By eye, these vessels were not visible from the skin surface. The largest vessel diameter as measured in the photograph was approximately 300μm .

Maximum data acquisition rate to form a single B-scan is presently limited by the 10-Hz laser pulse-repetition rate and hardware communication ( 2sec per B-scan). Future work will include the development of highly parallel data acquisition electronics. This, along with higher pulse-repetition rate lasers, should allow real-time imaging capabilities. Present imaging depths are largely limited by light delivery. With improved light delivery, imaging to greater depths should be likely. The present dynamic-receive focusing capabilities of our system offer a significant advantage over single-element mechanically scanned systems, which have a much smaller, predefined depth of field.

High-frequency photoacoustic imaging with arrays should enable visualization of key features of microvascular function and pathology. This technology may also be used to study microvessel constriction and dilation in response to physiological stimuli. Multiple optical wavelength approaches will enable assessment of microvascular oxygenation.10 It should be of significant interest to basic scientists studying angiogenesis in cancer, as well as clinicians such as dermatologists for assessing early stages of peripheral microvascular disease.

In summary, photoacoustic microscopy with high-frequency arrays should provide unique opportunities for clinical and biological impact, and the promising data in this letter should motivate future developments.


This project was sponsored by NIH grants R01 EB000712 and R01 NS46214.


1.  P. Carmeliet, “Angiogenesis in life, disease and medicine,” Nature (London)0028-0836 438, 932–936 (2005). Google Scholar

2.  K. Maslov, G. Stoica, and L. Wang, “In vivo dark-field reflection-mode photoacoustic microscopy,” Opt. Lett.0146-9592 10.1364/OL.30.000625, 30(6), 625–627 (2005). Google Scholar

3.  H. F. Zhang, K. Maslov, G. Stoica, and L. H. Wang, “Functional photoacoustic microscopy for high-resolution and noninvasive in vivo imaging,” Nat. Biotechnol.1087-0156 10.1038/nbt1220 24, 848–851 (2006). Google Scholar

4.  J. J. Niederhauser, M. Jaeger, R. Lemor, P. Weber, and M. Frenz, “Combined ultrasound and optoacoustic system for real-time high-contrast vascular imaging in vivo,” IEEE Trans. Med. Imaging0278-0062 10.1109/TMI.2004.843199 24(4), 436–440 (2005). Google Scholar

5.  D. E. Kruse and K. W. Ferrara, “A new high resolution color flow system using an eigendecomposition-based adaptive filter for clutter rejection,” IEEE Trans. Ultrason. Ferroelectr. Freq. Control0885-3010 49(12), 1739–1754 (2002). Google Scholar

6.  D. E. Goertz, D. A. Christopher, J. L. Yu, R. S. Kerbel, P. N. Burns, and F. S. Foster, “High frequency color flow imaging of the microcirculation,” Ultrasound Med. Biol.0301-5629 26, 6371 (2000). Google Scholar

7.  T. A. Ritter, T. R. Shrout, R. Tutwiler, and K. K. Shung, “A 30-MHz piezo-composite ultrasound array for medical imaging applications,” IEEE Trans. Ultrason. Ferroelectr. Freq. Control0885-3010 10.1109/58.985706 49(2), 217–230 (2002). Google Scholar

8.  K. A. Snook, J. Z. Zhao, C. H. Alves, J. M. Cannata, W. H. Chen, R. J. Meyer, T. A. Ritter, K. K. Shung, and M. F. Insana, “High frequency transducers for medical ultrasonic imaging,” Proc. SPIE0277-786X 3982, 92–99 (2000). Google Scholar

9.  M. L. Li, H. F. Zhang, K. Maslov, G. Stoica, and L. H. V. Wang, “Improved in vivo photoacoustic microscopy based on a virtual-detector concept,” Opt. Lett.0146-9592 10.1364/OL.31.000474 31(4), 474–476 (2006). Google Scholar

10.  K. Maslov, M. Sivaramakrishnan, H. F. Zhang, G. Stoica, and L. V. Wang, “Technical considerations in quantitative blood oxygenation measurement using photoacoustic microscopy in vivo,” Proc. SPIE0277-786X 6086, R1–11 (2006). Google Scholar

Roger J. Zemp, Rachel Bitton, K. Kirk Shung, Meng-Lin Li, George Stoica, Lihong V. Wang, "Photoacoustic imaging of the microvasculature with a high-frequency ultrasound array transducer," Journal of Biomedical Optics 12(1), 010501 (1 January 2007). https://doi.org/10.1117/1.2709850

Back to Top