|
1.IntroductionPhotoacoustic tomography (PAT) is an emerging biomedical imaging modality, which detects optical absorbers, such as blood vessels, inside tissue. Based on the photoacoustic mechanism,1, 2 PAT uses an nonionizing illumination source and is noninvasive. Moreover, by using diffused light instead of ballistic light, PAT can image deeper into tissue than other pure high-resolution optical imaging methods, such as optical coherence tomography (OCT) and two-photon tomography. PAT has been successfully applied in imaging both small-animal and human tissues.3, 4, 5 Among many ultrasonic detection methods for PAT, scanning with a single finite-size flat transducer is widely used owing to its simplicity and high sensitivity. However, the finite-size flat transducer not only decreases lateral resolution as the object approaches the transducer (aperture effect),6 but also limits the detection region due to its narrow acceptance angle. In contrast, a real point detector with a very small active surface suffers from poor sensitivity but has a wide acceptance angle and a negligible aperture effect. Thus, a point detector with a high sensitivity is highly desirable for PAT. Moreover, both exact and approximate image reconstruction algorithms have been developed for PAT with point detectors.7, 8 Virtual point detection methods have already been introduced in photoacoustic tomography.9, 10, 11, 12 By taking advantage of the integrated ultrasonic detectors, detected ultrasonic waves are primarily equivalent to waves that come through a “virtual point,” such as the ring center or the focal point. This method uses the location of the virtual point to calculate the time delay of the ultrasonic signal. Three kinds of virtual point detectors have been developed: 3-D focused detector,9 2-D ring detector,10 and 2-D high numerical aperture (NA) transducer.12 The first has been proved to increase the image resolution for objects close to the focal zone in photoacoustic microscopy (PAM). The second and third virtual point detectors have similar construction. The ring detector has already been demonstrated by both phantom and animal experiments in PAT. However, the ring virtual point detector limits its field of view (FOV) within the ring. We previously introduced a 2-D high-NA focused ultrasonic transducer for PAT with phantom experiments. This virtual point detector has a large physical detection surface and a small effective detection volume. Phantom experiments demonstrated that this virtual point detector has a comparable sensitivity to a finite-size flat transducer, as well as a much wider acceptance angle than a flat transducer. Owing to the negligible aperture effect, this virtual point detector can generate images with uniform resolution. In the next section, we briefly describe the virtual point detector. In Sec. 3, we image the cerebral cortex of a mouse in situ using this virtual point detector. 2.A High-NA-Based Virtual Point DetectorThe transducer’s active material was a metallized polyvinylidene fluoride (PVDF) film (from Measurement Specialties, Inc.) with a thickness of . The 2-D high-NA focused transducer was constructed by cutting the film in a -wide strip and gluing it on an acrylic plastic surface, as seen in Fig. 1 . The transducer had a half-circular shape with a radius of , forming a positively focused transducer with . The estimated center frequency was around . From the calculation at its center frequency,12 the focal width of this high-NA transducer at the 2-D plane is about . The performance of this 2-D virtual point detector had been previously demonstrated by phantom experiments.12 The width of the PVDF strip constrains the elevation resolution to be about —in other words, the image slice practically has a thickness of about . For comparison, we further constructed a square flat PVDF transducer, by . 3.Animal ExperimentsThe cerebral cortex of a euthanized Swiss Webster mouse (Harlan Sprague Dawley, Inc., Indianapolis, Indiana, ) was imaged in situ by PAT. The hair on the mouse head was gently depilated by using a hair removal lotion, and the mouse head was then fixed on a homemade animal holder. The experimental setup is shown in Fig. 2 . The illumination source was a Nd:YAG laser (Brilliant B, Quantel), which generated , laser pulses with a repetition rate of . The illumination laser was homogenized by a ground glass so that the whole cortex region of the mouse could be illuminated. A hole at the bottom of the water tank was sealed by a thin transparent membrane. The mouse head, after a layer of water-based gelatin was applied to it, protruded up into the tank against that membrane. Ultrasonic detectors, either a flat or a high-NA positively focused transducer, were immersed in water at the same horizontal plane as the mouse cortex. Both detectors evenly scanned the cortex along a horizontal circle, stopping at 240 points, and the signals were averaged 20 times at each stop. The PA signals were first amplified by a preamplifier and were then recorded by an oscilloscope (Tektronix TDS640A) with a sampling rate of . Last, the recorded signals were sent to a PC for image reconstruction. Our image reconstruction is based on the solid-angle-weighted image reconstruction algorithm, as described in Ref. 7. We use the scanning center as the coordinate origin. Detectors, flat or virtual point detectors, lie at location . We further simplified the algorithm by approximately using the detected signal instead of the actual pressure and its temporal derivative term. The algorithm used in this paper is where is the reconstructed value at location , is the detected signal at time , is the unit vector of , and is the speed of the sound.Figure 3a shows the result of scanning with a high-NA-based virtual point detector; the scanning radius was about . Figure 3b shows the result of scanning with a flat transducer; the scanning radius was about . Fig. 3Reconstructed PAT images. (a) From data gathered by a positively focused PVDF transducer; the scanning radius of the virtual point was . (b) From data gathered by a flat transducer; the scanning radius was . (c) Photograph of the mouse cortex taken after photoacoustic data acquisition. ![]() Although the scanning radius of the virtual point detector was shorter than that of the flat transducer, the reconstructed image of Fig. 3a exhibits a uniform resolution. Moreover, the reconstructed images of peripheral blood vessels in Fig. 3b are blurred. For instance, we marked two cortex vessels, 1 and 2. Compared with the photograph of the mouse cortex in Fig. 3c, both marked vessels were successfully reconstructed by using a virtual point detector [Fig. 3a], and the image quality is significantly improved over the image obtained by using a flat transducer [Fig. 3b]. This result was consistent with the results from phantom experiments in Ref. 12. To alleviate the image blurring, the finite-size flat transducer had to be placed much farther away from the scanning center, at the expense of convenience and signal strength. However, the finite-size flat transducer has higher signal-to-noise ratio (SNR) in imaging regions close to the scanning center. This is because it has a larger detection surface than the size of the virtual point, and the ultrasound wave from regions close to the center approximately perpendicularly reach the surface of the finite-size flat transducer. By comparison, the high-NA 2-D detector effectively detects only ultrasonic waves passing through the virtual point. 4.ConclusionsIn summary, a 2-D virtual point detector was successfully applied in PAT for cortex imaging. Compared with a flat transducer, the virtual point detector, with a similar sensitivity to the flat transducer, can image the object at a much closer distance with more uniform resolution. In addition, the closer the detector to the source, the higher the SNR. Such a PAT system can also be built more compactly and will have a negligible aperture effect. This virtual point detector is also important to thermoacoustic tomography (TAT), where the illumination source is replaced by a microwave source. An important potential application is that the virtual point method can be applied in breast imaging by PAT or TAT. Due to the size of the breast, it is inconvenient to place ultrasonic detectors far away from the breast. Moreover, this 2-D virtual point detector can be used to construct the PAT array, by either evenly placing multiple detectors along the detection circular trajectory or by making a stack of detectors for multilayer detection. In addition to PAT and TAT, this virtual point detector can also be potentially implemented in other fields that use ultrasonic detectors, such as ultrasonic tomography. AcknowledgmentsThis project was sponsored in part by National Institutes of Health Grant Nos. R01 NS46214(BRP) and R01 EB000712. L.W. has a financial interest in Endra, Inc., which, however, did not support this work. ReferencesM. H. Xu and L. H. V. Wang,
“Photoacoustic imaging in biomedicine,”
Rev. Sci. Instrum., 77 41101
–41122
(2006). 0034-6748 Google Scholar
L. V. Wang and H.-I. Wu, Biomedical Optics: Principles and Imaging, Wiley, Hoboken, NJ
(2007). Google Scholar
C. G. A. Hoelen, F. F. M. de Mul, R. Pongers, and A. Dekker,
“Three-dimensional photoacoustic imaging of blood vessels in tissue,”
Opt. Lett., 23 648
–650
(1998). https://doi.org/10.1364/OL.23.000648 0146-9592 Google Scholar
X. D. Wang, Y. J. Pang, G. Ku, X. Y. Xie, G. Stoica, and L. H. V. Wang,
“Noninvasive laser-induced photoacoustic tomography for structural and functional in vivo imaging of the brain,”
Nat. Biotechnol., 21 803
–806
(2003). https://doi.org/10.1038/nbt839 1087-0156 Google Scholar
R. I. Siphanto, K. K. Thumma, R. G. M. Kolkman, T. G. van Leeuwen, F. F. M. de Mul, J. W. van Neck, L. N. A. van Adrichem, and W. Steenbergen,
“Serial noninvasive photoacoustic imaging of neovascularization in tumor angiogenesis,”
Opt. Express, 13 89
–95
(2005). https://doi.org/10.1364/OPEX.13.000089 1094-4087 Google Scholar
M. Xu and L. V. Wang,
“Analytical explanation of spatial resolution related to bandwidth and detector aperture size in thermoacoustic or photoacoustic reconstruction,”
Phys. Rev. E, 67 056605
(2003). https://doi.org/10.1103/PhysRevE.67.056605 1063-651X Google Scholar
M. H. Xu and L. H. V. Wang,
“Universal back-projection algorithm for photoacoustic computed tomography,”
Phys. Rev. E, 71 016706
(2005). https://doi.org/10.1103/PhysRevE.71.016706 1063-651X Google Scholar
R. A. Kruger, D. R. Reinecke, and G. A. Kruger,
“Thermoacoustic computed tomography-technical considerations,”
Med. Phys., 26 1832
–1837
(1999). https://doi.org/10.1118/1.598688 0094-2405 Google Scholar
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., 31 474
–476
(2006). https://doi.org/10.1364/OL.31.000474 0146-9592 Google Scholar
X. M. Yang, M. L. Li, and L. H. V. Wang,
“Ring-based ultrasonic virtual point detector with applications to photoacoustic tomography,”
Appl. Phys. Lett., 90 251103
(2007). https://doi.org/10.1063/1.2749856 0003-6951 Google Scholar
X. M. Yang and L. V. Wang,
“Photoacoustic tomography of a rat cerebral cortex with a ring-based ultrasonic virtual point detector,”
J. Biomed. Opt., 12 060507
(2007). https://doi.org/10.1117/1.2823076 1083-3668 Google Scholar
C. Li and L. H. V. Wang,
“High-numerical-aperture-based virtual point detectors for photoacoustic tomography,”
Appl. Phys. Lett., 3 033902
(2008). https://doi.org/10.1063/1.2963365 0003-6951 Google Scholar
|