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.
©2013 Optical Society of America
Photoacoustic tomography (PAT) has been receiving growing attention in the last decade [1–17], since it combines good acoustic resolution in deep tissue with the optical absorption contrast. It has been widely applied to small animal imaging [1, 3], vasculature visualization [4, 5], osteoarthritis assessment [6–8], drug delivery monitoring , and so on. Recently, Sun et al. applied PAT to evaluate the tissue denaturation induced by high-intensity focused ultrasound treatment . Nie et al. imaged monkey brain  using PAT. Xia et al.  and Wang et al.  developed the three-dimensional PAT systems. PAT is usually limited to the biological tissues with relatively homogeneous acoustical properties . Many efforts, including ultrasound tomography , time reversal method [19, 20], statistical reconstruction , and coherence factor optimization  have been done to improve the PAT in inhomogeneous tissue. However, these studies usually need known some prior knowledge of tissue inhomogeneity [19,21] or only consider some relatively simple scenario of tissue inhomogeneous, such as speed inhomogeneity [18, 22] or single scattering .
Biological tissue is essentially a kind of acoustically stochastic medium. A large number of acoustically stochastic scatterers with the size comparable to wavelength could generate multiple scattering and randomly distort the propagation of photoacoustic (PA) signal. This random multiple scattered field is usually referred to as speckle noise . Speckle noise inherently exists in acoustical inhomogeneous tissue. It could bring about artifacts and distortion, and thus degrade image quality [19,21]. Therefore, speckle noises coming from acoustically inhomogeneous tissue are usually considered as nuisance and ignored in classical PAT. Speckle noise remains challenging the PAT.
In this study, an ultrasound-photoacoustics (US-PA) scheme is proposed to break the limitation of PAT in acoustically inhomogeneous tissue. An ultrasonic (US) method is used to determine the Green's function of inhomogeneous. Then, with the known Green’s function, we extend the PAT to speckle environment. Finally, we use the proposed method to improve the limited-view PAT by extracting the back-propagating information from the speckle noise.
Figure 1 represents the scenario considered in this study. The region of interest (ROI) is located below the surface of the tissue. A transducer array is attached on the surface to pick up the PA signal coming from the tissue. The randomly distributed circles represent scatterers deep in the tissue which scatter PA signals and give rise to speckle noise. In a PAT system, laser beam is expanded, diffused, and then irradiated on the tissue uniformly [1, 15]. After tissues are illuminated by pulsed laser, thermal expansion induced an initial acoustic sources q(r0, t0), which is positive proportional to optical energy absorption deposition A(r0), i.e., q(r0, t0) ~A(r0)∂H(t0)/∂t0 . H(t0) is the temporal profile of the laser pulse. The generated PA pressure p(rd, t) at position rd and time t (0 ≤ t ≤ T) isFig. 1) is propagating though the homogeneous tissue. Part of PA signal (pd2 shown in Fig. 1) is scattered and reflected (ps in Fig. 1) by the speckle tissue. Therefore, the Green's function g(rd, t | r0, t0) also consists of two terms g(rd, t | r0, t0) = gh(rd, t | r0, t0) + gs(rd, t | r0, t0), where gh and gs correspond to the homogeneous and speckle tissues, respectively. We can rewrite Eq. (1) as
Time reversing the detected PA signal pd(rd, t) to pd(rd, T – t) and reemitting the reversed signals pd(rd, T – t), we have the time reversal field pTR(r0, t),20, 24, 25], we have pTR(r0, t) ≈p(r0, T–t) and the time reversal wave will eventually focus toward the originally PA sources, i.e., pTR(r0, T) ~A(r0). The optical absorption distribution can be reconstruction from the detected PA signal by integrating Eq. (3). Therefore, Eq. (3) provides an algorithm to reconstruct PA image from the detected PA signal. For the homogeneous tissue with the known velocity of sound c, the Green’s function is given as g0(rd, t | r, t1) = δ(t–t1–|r0–rd|/c)/4π|r0–rd|. However, for the acoustically inhomogeneous tissue, the Green’s function gs(rd, t | r, t1) is unknown, the reconstruction of PAT become difficult. This is the reason that the applications of PAT are often prohibited from acoustically inhomogeneous tissue .
If the Green’s function gs(rd, t | r, t1) is determined, PAT can be extracted from speckle noise. Fortunately, US technique depends primarily on acoustic heterogeneity and it could provide the information of acoustically inhomogeneous properties of tissues [23, 26]. Therefore, US technique could determine the Green’s function of inhomogeneous tissue and thus make up the limitation of PAT in speckle tissue. This is the basic idea of our US-PA scheme:
Firstly, in order to estimate the Green’s function of speckle tissue, an N-element US transducer array is employed to emit and receive ultrasound, as shown in Fig. 1. A US pulse with a waveform u(t) is transmitted by the i-th (i = 1, 2, …, N) element of transducer array. The scattered field reflected by the speckle tissue is detected by the same transducer array, where the signals recorded by the j-th (j = 1, 2, …, N) element are denoted as h(j, i, t). Repeating the US pulse emitting and receiving for each element of the transducer array, we obtain N × N impulse responses h(j, i, t) (i, j, = 1, 2, …, N) between each pair of elements. With these impulse responses, the Green's function gs can be represented in frequency domain by 
For a discrete ultrasound array, the integral in Eq. (4) can be approximated by the summation,
Secondly, the ROI is illuminated by laser pulse and the excited PA signal is detected by the same US transducer. With the recorded PA signals and the estimated Green's function g(rd, t | r, t1) = gh(rd, t | r, t1) + gs(rd, t | r, t1), we can extract PAT images from the speckle noise by integrating Eq. (3).
The proposed US-PA scheme takes the advantages of both the US technique and PAT technique, which make it possible to extract PAT from acoustically inhomogeneous tissue. Moreover, the validities of the Green’s function estimation [Eq. (5)] and the image reconstruction [Eq. (3)] do not depend on the prior knowledge of tissue and the complexity of the medium perturbation. Therefore, this US-PA scheme could be generally applied to various scenarios, like scattering from stochastic tissue.
Numerical experiments are carried out to validate of the proposed method. As shown in Fig. 1, the randomly distributed white circles represent acoustic scatterers with a diameter of 1.6 mm. The two dark gray circles at the center of the ROI denote optical absorbers with a diameter of 2.0 mm. The relative optical energy absorption deposition A(r0) of optical absorbers is 1.0, while A(r0) of the medium is 0.0. The center of the shallow absorber is 3.0 mm distance from the US transducer and the deep one is 7.0 mm distance from the US transducer. In the following discussion, the deep absorber indicates the absorber, which is far from the US transducer, while the shallow absorber means the one, which is close to the transducer, as shown in Fig. 1. The acoustic parameters of the medium are c = 1500 m/s and ρ = 1000 kg/m3. An US transducer array with 16 elements (N = 16) is utilized to transmit and record US pulse for the Green’s function estimation. The same transducer array is also used to detect PA signals generated by laser pulse.
US pulse with a central frequency of 3.5 MHz is emitted by the transducer array. The response signals h(j, i, t) within time duration 50 µs are recorded. The Green's function Gs is estimated from the impulse response signals according to Eq. (5). PA signals from the optical absorbers are recorded, as shown in Fig. 2(a). With the known Green's function and detected photoacoustic signals, images are reconstructed according to Eq. (3).
For the sake of comparison, PA images are also reconstructed by only using the direct signal pd1 propagating through homogeneous tissue. In this situation, the speckle signal is ignored and the reconstruction equation [Eq. (3)] degraded to pTR(r0, t) = ∫dt1 ∫∫Σdrdpd(rd, T – t1)gh(r0, t | rd, t1). Using the summation instead of the integral in the above equation, we have27]. For the sake of simplicity, we let w(i, r0) = 1 in the following simulation. By ignoring speckle noise, the reconstruction equation is degraded to the widely used delay and sum method .
As shown in Fig. 2(b), when the speckle signal is ignored for the image reconstruction, the image of absorbers is seriously distorted due to the absence of the backside PA signals [29–31]. Moreover, although the two absorbers have the same optical absorber coefficient, they are presented significantly intensity difference in the reconstructed image . The reconstructed intensity of the deep absorber is only about 60% of that of the shallow one. This is because that the US linear array with a finite length does not enclose the ROI and can only capture the forward-propagating photoacoustic waves (e.g., pd1 in Fig. 1), which results in the limited-view PAT. The incomplete PA data will induce artifacts, distortion [20, 30], and false intensity contrast of absorbers located in different depth .
Figure 2(c) illustrates the PA image extracted from speckle noise by using the proposed method. When speckle noises are involved for imaging reconstruction, the distortion is restrained and the false intensity contrast between the two absorbers is decreased. The image quality is better than that without speckle noise. This phenomenon could be explained that the backside information carried by the speckle noise improves the image quality. As we have discussed, the US linear array only capture the forward-propagating photoacoustic waves (pd1) and the speckle signal (e.g., ps). The back-propagating PA wave (pd2) is missed. However, after the random scattering by stochastic tissue, the information carried by the back-propagating PA wave (pd2) is distorted, but it remains existence in the speckle noise (ps). They cannot be utilized by the classical PAT scheme due to the absence of the Green’s function of heterogeneous tissue. However, the US-PA scheme utilizes the US method to determine the Green’s function of speckle tissue and make up the PAT limitation. It allows us to extract the backside information of the ROI encoded in speckle noise for the image reconstruction. Therefore, the image quality is improved.
Multiple scattering induced by stochastic tissue expands a narrow pulse signal to a speckle noise with a long duration. The influence of signal length on the image quality is important. Figure 3 present the relationship between signal length and image quality, where the full width at half maximum (FWHM) and the maximum amplitude ratio γ between two absorbers are chosen to quantify distortion and false contrast, respectively. As shown in Fig. 3(a), when only direct wave (the first 5µs of signal pd) are employed for the image reconstruction, the FWHM of the deep absorber is only about 0.8 mm [Fig. 3(a)], which has far deviated from its actual diameter 2.0 mm. Moreover, the image intensity of the deeper absorber is much lower than that of the shallow one [Fig. 3(b)]. When more and more low-order scattering wave (about the first 20 µs of the detected PA signal) are included for imaging, the FWHM is approaching toward the actual value [Fig. 3(a)] and the false contrast between two absorbers are decreased [Fig. 3(b)]. Finally, the further increasing of signal length will only lead to relatively minor improvement. It could because the high-order scattering wave is generally much weaker than the direct wave and the low-order scattering wave.
The length of US array not only determines the estimation of Green’s function but also the recording of PA signal. Figure 4 presents the performance of the US-PA scheme with various transducer array lengths. In general, the image extracted from the signals including the speckle noise is better than those by ignoring speckle noise. Moreover, both FWHM and amplitude ratio are improved with the length of the transducer array. This change can be explained from two aspects. First, a longer US array will capture more direct wave, which decrease the distortion and false contrast induced by the incomplete data [20, 30, 31]. Second, a longer US array will increase the accuracy in approximating the Green's function. Additionally, it is noticed that when the transducer array is short (< 20 mm), the image size of the shallow absorber is better if speckle signal is not included. This could be the result of the estimation error of the Green’s function due to the short US transducer. A shorter US array will decrease the accuracy in approximating the Green's function. This error could bring additional distortion to the reconstructed image.
An US-PA imaging scheme is proposed to break the limitation PAT in acoustically inhomogeneous tissue. The US technique is carried out to determine the Green's function between each transducer element and the ROI. By utilizing US to make up the PAT limitation in acoustically heterogeneous tissue, PA images could be extracted from speckle noise with the determined Green’s function. Studies show that the method can effectively decrease the distortion and false intensity contrast in limited-view PA images. The information extracted from the speckle noise significantly improves the image quality. The US-PA method does not depend on the prior knowledge of tissue and the complexity of the medium perturbation. Moreover, the estimation of Green’s function and the detection of PA signal are performed by the same ultrasound transducer. Therefore, the proposed scheme could be easily integrated into a classical PAT system and improve PAT in acoustically inhomogeneous tissue.
This work was supported by the National Basic Research Program of China under Grant No. 2012CB921504, the Natural Science Foundation of China under Grant Nos. 11274167 and 11274171, and SRFDP Grant No. 20120091110001.
References and links
1. X. Wang, Y. Pang, G. Ku, X. Xie, G. Stoica, and L. V. Wang, “Noninvasive laser-induced photoacoustic tomography for structural and functional in vivo imaging of the brain,” Nat. Biotechnol. 21(7), 803–806 (2003). [CrossRef] [PubMed]
3. J. Gamelin, A. Maurudis, A. Aguirre, F. Huang, P. Guo, L. V. Wang, and Q. Zhu, “A real-time photoacoustic tomography system for small animals,” Opt. Express 17(13), 10489–10498 (2009). [CrossRef] [PubMed]
7. Y. Sun, E. S. Sobel, and H. Jiang, “First assessment of three-dimensional quantitative photoacoustic tomography for in vivo detection of osteoarthritis in the finger joints,” Med. Phys. 38(7), 4009–4017 (2011). [CrossRef] [PubMed]
8. J. Xiao, L. Yao, Y. Sun, E. S. Sobel, J. He, and H. Jiang, “Quantitative two-dimensional photoacoustic tomography of osteoarthritis in the finger joints,” Opt. Express 18(14), 14359–14365 (2010). [CrossRef] [PubMed]
9. J. R. Rajian, M. L. Fabiilli, J. B. Fowlkes, P. L. Carson, and X. Wang, “Drug delivery monitoring by photoacoustic tomography with an ICG encapsulated double emulsion,” Opt. Express 19(15), 14335–14347 (2011). [CrossRef] [PubMed]
10. Y. Sun, K. C. P. Li, and B. O'Neill, “Multispectral photoacoustic imaging of tissue denaturation induced by high-intensity focused ultrasound treatment,” Proc. SPIE 8581, 85813H, 85813H-5 (2013). [CrossRef]
12. J. Xia, Z. Guo, K. Maslov, A. Aguirre, Q. Zhu, C. Percival, and L. V. Wang, “Three-dimensional photoacoustic tomography based on the focal-line concept,” J. Biomed. Opt. 16(9), 090505 (2011). [CrossRef] [PubMed]
13. Y. Wang, T. N. Erpelding, L. Jankovic, Z. Guo, J. L. Robert, G. David, and L. V. Wang, “In vivo three-dimensional photoacoustic imaging based on a clinical matrix array ultrasound probe,” J. Biomed. Opt. 17(6), 061208 (2012). [CrossRef] [PubMed]
14. Z. Guo, S. Hu, and L. V. Wang, “Calibration-free absolute quantification of optical absorption coefficients using acoustic spectra in 3D photoacoustic microscopy of biological tissue,” Opt. Lett. 35, 2067–2069 (2010). [CrossRef] [PubMed]
15. C. Kim, T. N. Erpelding, L. Jankovic, M. D. Pashley, and L. V. Wang, “Deeply penetrating in vivo photoacoustic imaging using a clinical ultrasound array system,” Biomed. Opt. Express 1(1), 278–284 (2010). [CrossRef] [PubMed]
17. M. Xu and L. V. Wang, “Photoacoustic imaging in biomedicine,” Rev. Sci. Instrum. 77(4), 041101 (2006). [CrossRef]
19. D. Wu, C. Tao, and X. Liu, “Photoacoustic tomography in scattering biological tissue by using virtual time reversal mirror,” J. Appl. Phys. 109(8), 084702 (2011). [CrossRef]
20. D. Wu, X. Wang, C. Tao, and X. J. Liu, “Limited-view photoacoustic tomography utilizing backscatterers as virtual transducers,” Appl. Phys. Lett. 99(24), 244102 (2011). [CrossRef]
21. X. L. Deán-Ben, V. Ntziachristos, and D. Razansky, “Statistical optoacoustic image reconstruction using a-prior knowledge on the location of acoustic distortions,” Appl. Phys. Lett. 98(17), 171110 (2011). [CrossRef]
22. C. Yoon, J. Kang, S. Han, Y. Yoo, T.-K. Song, and J. H. Chang, “Enhancement of photoacoustic image quality by sound speed correction: ex vivo evaluation,” Opt. Express 20(3), 3082–3090 (2012). [CrossRef] [PubMed]
23. J.-L. Robert and M. Fink, “Green’s function estimation in speckle using the decomposition of the time reversal operator: Application to aberration correction in medical imaging,” J. Acoust. Soc. Am. 123(2), 866–877 (2008). [CrossRef] [PubMed]
25. E. Bossy, K. Daoudi, A.-C. Boccara, M. Tanter, J.-F. Aubry, G. Montaldo, and M. Fink, “Time reversal of photoacoustic waves,” Appl. Phys. Lett. 89(18), 184108 (2006). [CrossRef]
26. I. Vasconcelos, R. Snieder, and H. Douma, “Representation theorems and Green’s function retrieval for scattering in acoustic media,” Phys. Rev. E Stat. Nonlin. Soft Matter Phys. 80(3), 036605 (2009). [CrossRef] [PubMed]
27. F. L. Lizzi, M. Greenebaum, E. J. Feleppa, M. Elbaum, and D. J. Coleman, “Theoretical Framework for Spectrum Analysis in Ultrasonic Tissue Characterization,” J. Acoust. Soc. Am. 73(4), 1366–1373 (1983). [CrossRef] [PubMed]
28. 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(1), 89–95 (2005). [CrossRef] [PubMed]
31. D. Wu, C. Tao, X. Liu, and X. Wang, “Influence of limited-view scanning on depth imaging of photoacoustic tomography,” Chin. Phys. B 21(1), 014301 (2012). [CrossRef]