We designed, fabricated and tested a novel imaging system that fuses diffuse optical tomography (DOT) and photoacoustic tomography (PAT) in a single platform. This platform takes advantages of both DOT and PAT, and can potentially provide dual-modality two dimensional functional and cellular images of the breast quantitatively. Here we describe this integrated platform along with initial tissue phantom validations.
© 2011 OSA
Photoacoustic tomography (PAT) is a potentially powerful imaging technique for visualizing the internal structure of tissue with excellent spatial resolution and high optical contrast [1–7]. While absorption coefficient images of heterogeneous media can now be recovered by quantitative PAT [8,9], it still remains a major challenge to obtain scattering coefficient from photoacoustic measurements .
Diffuse optical tomography (DOT), on the other hand, is another rapidly growing modality due to its high contrast in both tissue absorption and scattering, and tissue functional information available from multispectral DOT [11–15]. However, the spatial resolution of DOT is relatively low, and the detection of small targets is often impossible or distorted significantly.
In this paper, we report a novel integrated PAT/DOT system that combines the advantages of both PAT and DOT. In our hybrid modality, high resolution absorption and its derived functional images are generated through PAT, while scattering images are produced by DOT. We designed and fabricated 32 PVDF transducers, and attached them to the source/detector fiber optic array of a multispectral DOT system . Diffused light and ultrasound signals are collected by DOT and PAT systems, respectively. Phantom experiments are used to validate the performance of the integrated system. Quantitative absorption and scattering images are obtained using our finite element (FE) based DOT and PAT reconstruction algorithms.
Our DOT system was previously described in detail in Ref 11. For DOT data acquisition (Fig. 1(a) , left), light from a diode laser at 775nm was transmitted sequentially to 16 source points at the phantom surface through an optical switch, and diffusing light was detected by 16 photodiodes. A set of 16 × 16 measured data was then input into the DOT reconstruction algorithm to generate 2D cross-sectional images of the absorption and scattering coefficients of the phantom.
In our PAT system (Fig. 1(a), right), an Nd:YAG laser generated a pulsed 532nm beam with a pulse repetition rate of 10Hz and a pulse width of <10ns, which was delivered to the top surface of the phantom via a reflection mirror and combined concave mirror/ground glass for beam extension (The incident optical fluence was controlled below 20mJ/cm2 which is the safety limit). 32 homemade transducers were used to detect the acoustic signals at 32 positions along a circular path around the phantom (Figs. 1(b) and 1(c)). Each transducer was fabricated in-house using a 110μm-thick Ag ink printed PVDF film as the sensing unit (5 × 30mm) and backing materials. A 3mm diameter hole (Fig. 1(b)) was drilled/eroded in the center of backing materials and the corresponding area of Ag ink on PVDF film for the transmission/receiving of DOT signals. The frequency response of each PVDF transducer had a bandwidth of up to 2MHz. These 32 transducers were grouped to 4 with each containing 8 transducers for signal acquisition. The complex acoustic wavefield signals were amplified by a multi-channel preamplifier and further amplified by an amplifier with a controllable gain from 30dB to 40dB. The data acquisition sampling rate was 50MHz.
We demonstrate the ability of our integrated imaging system through several tissue-like phantom experiments with different contrasts between the target and the background. Figure 2 depicts the geometrical configuration for the test cases which consists of a circular background region (radius = 35mm) with an embedded circular target (radius = 3mm) offsetting 10mm. Three tissue-like phantom experiments were conducted. Tissue absorption and scattering were simulated with India ink and Intralipid, respectively. Agar powder (2%) was used to solidify the mixed Intralipid–India ink solution. Table 1 lists the absorption and reduced scattering coefficients of the target and background for the three phantom experiments.
FE based dual-meshing reconstruction algorithms described in Refs. 8,9 were used to recover quantitative absorption images from PAT data. A fine mesh of triangular elements with 6285 nodes and a coarse mesh with 1604 nodes were applied in the reconstruction, and the images were converged within 50 iterations in a parallel computer.
DOT images were reconstructed using our FE based reconstruction algorithms detailed in Refs. 13,14. In the reconstruction, a single mesh of 700 nodes and 1334 triangular elements was applied, and the images were converged within 15 iterations in a 3GHz PC with 0.99-GB memory.
Results and Discussion
Reconstructed results are shown in Figs. 3 and 4 where the pink-colored region (at 3 o’clock) clearly indicates the location and size of the object (Figs. 3(a), 3(c), 3(e), 3(g), 3(i), 3(k), 4(a) and 4(c)). As can be seen, the images formed are shown to be qualitatively correct in visual content for all the cases. To provide a more quantitative assessment of these images, the reconstructed optical property distributions are displayed along one transect through the center of reconstructed target—the one perpendicular to the line through the centers of both the target and background regions—compared with the exact values, in Figs. 3(b), 3(d), 3(f), 3(h), 3(j), 3(l), 4(b) and 4(d). Table 2 provides the quantitative information about the center and the full-width-half-maximum (FWHM) of the object as well as the recovered average optical properties of the object estimated from a region-of-interest (ROI) determined by the center and the FWHM of each object.
From the recovered absorption coefficient images, while we can see that the value of absorption coefficient is well reconstructed both from DOT and PAT data for all the cases, the location and size of target are much better recovered by PAT than that by DOT. For example, the relative errors of the recovered absorption coefficient by DOT are, 2%, 14%, and 4%, for cases 1-3, respectively, while the relative errors of the reconstructed absorption coefficient by PAT are, 16%, 1%, and 3%, for cases 1-3, respectively. For the DOT images, while the target is detected, we note that the reconstructed target position is slightly shifted relative to the exact target position. The PAT recovered target size is found to be 7.1mm, 5.3mm, and 6.1mm for the three cases, respectively, whereas the DOT recovered target size (11-14mm) is significantly overestimated compared to the exact size.From Fig. 3, we noticed that the recovered target size and absorption value are more accurate when the contrast is increased for both DOT and PAT. The PAT image showed some target distortion (Fig. 3(a)) when only absorption contrast existed. This could be due to the use of a small number of transducers and the lack of target scattering contribution to the PAT reconstruction through the optical fluence/photon density. This can be improved by adding more transducers to the PAT system.
We also noted that some artifacts appear in the DOT images especially for case 1 (Fig. 3(b)). These artifacts generally are shown near the boundary source and detector positions where the measurement sensitivity is highest. However, such boundary artifacts are not seen in the PAT absorption images.
From Fig. 4, we see that the scattering images for both cases are recovered with high quality without boundary artifacts by DOT. In addition, we found that the recovered size of the target ranges from 5.9mm to 6.2mm, in good agreement with the actual object size of 6mm. The recovered off-center location of the target ranges from 9.4mm to 9.6mm, also in good agreement with the actual off-center location of 10mm. Finally, we see that the scattering coefficient values of target and background are well recovered quantitatively as seen from Table 2.
In sum, we have presented a hybrid PAT and DOT system that combines the advantages of both PAT and DOT. We showed that quantitative absorption and reduced scattering coefficient images can be obtained from this system using tissue-like phantom experiments. We are currently improving this hybrid modality by adding 64 transducers to the PAT system, and will apply it to the in vivo detection of breast cancer after its considerable phantom evaluation. In addition, we are implementing ideas that use the higher resolution absorption image by PAT as a priori knowledge for improved recovery of scattering image by DOT, or that consider the inhomogeneous scattering distribution provided by DOT to enhance the quantitative recovery of absorption coefficient by PAT, similar to the work described in Ref. 16 where we used the diffusing light measurements to assist in quantitative reconstruction of absorption coefficient by PAT. We plan to report the results from these studies in the near future.
This research was supported in part by the J. Crayton Pruitt Family Endowment and by the Bankhead Coley Florida Cancer Research Program.
References and links
1. S. A. Ermilov, T. Khamapirad, A. Conjusteau, M. H. Leonard, R. Lacewell, K. Mehta, T. Miller, and A. A. Oraevsky, “Laser optoacoustic imaging system for detection of breast cancer,” J. Biomed. Opt. 14(2), 024007 (2009). [CrossRef]
2. M. Pramanik, G. Ku, C. H. Li, and L. V. Wang, “Design and evaluation of a novel breast cancer detection system combining both thermoacoustic (TA) and photoacoustic (PA) tomography,” Med. Phys. 35(6), 2218–2223 (2008). [CrossRef]
3. G. Ku, B. D. Fornage, X. Jin, M. Xu, K. K. Hunt, and L. V. Wang, “Thermoacoustic and photoacoustic tomography of thick biological tissues toward breast imaging,” Technol. Cancer Res. Treat. 4(5), 559–566 (2005). [PubMed]
4. D. Piras, W. Xia, W. Steenbergen, T. G. Leeuwen, and S. Manohar, “Photoacoustic imaging of the breast using the twente photoacoustic mammoscope: present status and future perspectives,” IEEE J. Sel. Top. Quantum Electron. 16(4), 730–739 (2010). [CrossRef]
5. S. Manohar, S. E. Vaartjes, J. C. G. van Hespen, J. M. Klaase, F. M. van den Engh, W. Steenbergen, and T. G. van Leeuwen, “Initial results of in vivo non-invasive cancer imaging in the human breast using near-infrared photoacoustics,” Opt. Express 15(19), 12277–12285 (2007). [CrossRef]
6. S. Manohar, A. Kharine, J. C. G. van Hespen, W. Steenbergen, and T. G. van Leeuwen, “Photoacoustic mammography laboratory prototype: imaging of breast tissue phantoms,” J. Biomed. Opt. 9(6), 1172–1181 (2004). [CrossRef]
7. R. A. Kruger, R. B. Lam, D. R. Reinecke, S. P. Del Rio, and R. P. Doyle, “Target detection and quantification using a hybrid hand-held diffuse optical tomography and photoacoustic tomography system,” Med. Phys. 37(11), 6096–6100 (2010). [CrossRef]
9. Z. Yuan and H. Jiang, “Quantitative photoacoustic tomography: Recovery of optical absorption coefficient maps of heterogeneous media,” Appl. Phys. Lett. 88(23), 231101 (2006). [CrossRef]
10. B. T. Cox, J. G. Laufer, and P. C. Beard, “The challenges for quantitative photoacoustic imaging,” Proc. SPIE 7177, 717–713 (2009).
11. C. Li and H. Jiang, “Imaging of particle size and concentration in heterogeneous turbid media with multispectral diffuse optical tomography,” Opt. Express 12(25), 6313–6318 (2004). [CrossRef]
12. H. Jiang, N. V. Iftimia, Y. Xu, J. A. Eggert, L. L. Fajardo, and K. L. Klove, “Near-infrared optical imaging of the breast with model-based reconstruction,” Acad. Radiol. 9(2), 186–194 (2002). [CrossRef]
13. C. Li, H. Zhao, B. Anderson, and H. Jiang, “Multispectral breast imaging using a ten-wavelength, 64 x 64 source/detector channels silicon photodiode-based diffuse optical tomography system,” Med. Phys. 33(3), 627–636 (2006). [CrossRef]
14. H. Jiang, K. D. Paulsen, U. L. Osterberg, and M. S. Patterson, “Improved continuous light diffusion imaging in single- and multi-target tissue-like phantoms,” Phys. Med. Biol. 43(3), 675–693 (1998). [CrossRef]
16. L. Yin, Q. Wang, Q. Zhang, and H. Jiang, “Tomographic imaging of absolute optical absorption coefficient in turbid media using combined photoacoustic and diffusing light measurements,” Opt. Lett. 32(17), 2556–2558 (2007). [CrossRef]