|
1.IntroductionThe mechanical properties of tissue, such as stiffness and viscosity, are often correlated to pathological changes. For example, during the growth of malignant tumors, healthy tissues are often replaced by stiffer fibrotic tissues.1 In fatty liver disease, large vacuoles of fat accumulate in liver cells, leading to an increase of liver viscosity.2 For decades, researchers have been developing tissue viscoelasticity imaging to measure tissue mechanical properties noninvasively and quantitatively. One approach is to push tissue with either a static or dynamic force and to inversely estimate elastic modulus from the resulting tissue response.3 An alternative method is to generate and track shear-wave propagation in tissue and to calculate the elastic modulus by measuring the shear-wave velocity.4 This second method suffers fewer boundary effects,5 and in addition to the elasticity, tissue viscosity can be evaluated by measuring the shear-wave dispersion, e.g., measuring shear-wave attenuation or velocity at various shear-wave frequencies, and the viscoelasticity can be obtained from the Voigt’s model.6 Magnetic resonance imaging and ultrasound have shown a good capacity for shear-wave tracking in tissue.7,8 However, tracking shear waves with optical methods may bring additional advantages—it is sensitive at the optical wavelength scale and offers the opportunity to combine optical imaging with viscoelasticity imaging, providing complementary optical information for clinical diagnosis. In our previous study,9 we proposed an optical method to track transient shear waves at centimeter depth within tissue-mimicking phantoms by shear-wave laser speckle contrast analysis (SW-LASCA). Based on a time-of-flight method, local shear-wave speed was calculated by finding the time shift of a time-resolved contrast difference signal due to the transient shear waves generated at two locations separated by a certain distance. Elastography of a heterogeneous phantom was then achieved based on the measured shear-wave speeds. In this article, with the same setup, we proposed a method to track continuous shear-wave propagation and to quantify the viscosity of tissue-mimicking phantoms by measuring the dispersion of the shear-wave attenuation. First, the principle of SW-LASCA is briefly introduced, including the methodology and experimental setup. Then, the method to measure shear-wave attenuation dispersion is described, and experiments are conducted on phantoms with various viscosities. The viscosity is then calculated from the measured shear-wave dispersion and compared with the values in literature. Finally, both advantages and limitations for this method are discussed and summarized. 2.MethodThe method of SW-LASCA was described in Ref. 9, and a brief summary is provided here. Speckle pattern images are formed when coherent light is transmitted trough a turbid medium and detected by a charged-coupled device (CCD) camera. Contrast of speckle pattern images is calculated by where and are the standard deviation and mean of CCD pixel intensities, respectively. is close to unity if there is no movement within the medium, and decreases otherwise. In SW-LASCA, shear waves introduce movement of the optical scatterers in an optical turbid medium, and the time-resolved speckle contrast difference () is calculated to reveal the shear-wave propagation. The contrast difference () is the reduction of the speckle contrast with and without () shear waves:Figure 1(a) shows our SW-LASCA setup. The tissue-mimicking phantom was fixed in a perspex tank filled with degased water. A 532-nm Nd:YAG long-coherence laser (Excelsior 532, Newport Inc., Irvine, CA) illuminated the phantom, and the transmitted laser speckle patterns were recorded on a CCD camera (Qimaging Retiga EXi, Surrey, BC, Canada). The CCD camera was aligned with the laser axis to allow maximized light intensity to be received. Continuous shear waves were generated by an amplitude-modulated ultrasound using a focused ultrasound transducer (Panametric NDT Videoscan 307, 5-MHz central frequency, 5-cm focal length, Olympus, Essex, UK) tens of millimeters away from the optical axis. The full-width–at-half-maximum of the ultrasound focus was 1 cm in the axial direction and 0.1 cm in the lateral direction. Shear waves propagating toward the scattered light region [the green region in Fig. 1(a)] modulate the light and thus the CCD speckle contrast difference signal. Figure 1(b) is a typical time-resolved contrast difference signal due to the continuous shear-wave propagation, where 400-Hz shear waves are generated. The signal was produced by recording CCD speckle images with different time delays (-axis) after the ultrasound has been launched. As shown, the contrast difference signal varies with the shear-wave frequency over time. The distance between laser and ultrasound axes can be modified by moving the transducer with a translation stage, as shown in Fig. 1(a). Therefore, local shear-wave information, such as speed and attenuation coefficient, could be obtained from the time-resolved contrast difference signals produced by the shear waves originating from two ultrasound locations.9 Fig. 1(a) Shear-wave laser speckle contrast analysis (SW-LASCA) experimental setup. (b) Time-resolved speckle contrast difference to track shear wave. ![]() In a viscoelastic medium, such as tissue, shear-wave attenuation may change at different frequencies. This phenomenon is called shear-wave attenuation dispersion. The Voigt model10 may be used to describe the shear attenuation dispersion as below where () is the shear-wave attenuation coefficient, () is the angular frequency of the shear wave, and (), (Pa), (Pa s) are tissue density, elasticity, and viscosity, respectively. If is measured at different shear-wave frequencies, elasticity and viscosity can be estimated with Eq. (3).To measure the shear-wave attenuation with the time-resolved speckle contrast difference signals, the following relationships are deduced. First, the relationship between shear-wave attenuation coefficient and amplitude is obtained as where is the shear-wave amplitude at position , is the initial amplitude of shear wave at the ultrasound focus, is the attenuation coefficient of the shear wave at frequency , and is the shear-wave path length between the focus and the position . According to Eq. (4), the shear-wave attenuation coefficient could be derived if , , and are known.In the experiments, , which is the distance between the laser and ultrasound axes, is known. is assumed to be linear with the magnitude of the acoustic radiation force (ARF), and the ARF is proportional to the time-averaged ultrasound intensity .11 Therefore, the relationship between shear-wave initial amplitude and can be related as where is a constant. On the other hand, the time-averaged ultrasound intensity, , is proportional to the mean square of the ultrasound pressure. With a hydrophone measurement, could be further related to the input voltage of ultrasound transducer. Figure 2 shows the hydrophone measurement result (red dots), where the -axis is the ultrasound input voltage and the -axis is the ultrasound intensity (mean square of acoustic pressure). With linear regression (blue solid curve), the relationship between the time-averaged ultrasound intensity and the transducer input voltage (between 300 and 600 mVpp) was found to beBy substituting Eqs. (5) and (6) into Eq. (4), we can obtain In Eq. (7), and are unknown, and is the averaged attenuation coefficient over . In order to eliminate and to calculate the local shear-wave attenuation coefficient, shear waves were generated at two locations with path lengths of and . The ratio of the corresponding shear-wave amplitudes is calculated as In Eq. (8), the unknown constant is eliminated. To calculate the left side of Eq. (8), we assume that the identical shear-wave amplitudes at the detection position result in identical contrast differences, which means if , then , despite of the original position of the shear waves. Consequently, Eq. (8) can be rewritten as and the shear-wave attenuation coefficient at frequency can be calculated byBelow, we describe how we found () pair in Eq. (10). Figure 3 shows the relationship between ultrasound input voltage and contrast difference induced by a 400-Hz shear wave, in which the markers are the mean values of the contrast difference signals over time [e.g., averaged value from Fig. 1(b)]. P1 and P2 are the two positions where shear waves are generated. Because P2 is located further away from the detection area than P1, the corresponding mean contrast difference values are generally smaller because of the shear-wave attenuation. By changing the ultrasound input voltage, the mean contrast differences increase for both P1 and P2. Solid lines in the figure indicate linear regressions of the data. As indicated by the dashed lines in the figure, for a given , and can be obtained from the fitted curves to calculate the shear-wave attenuation coefficient using Eq. (10). Fig. 3Relationship between input voltage of ultrasound transducer and mean speckle contrast difference. Shear-wave frequency is 400 Hz. ![]() Shear waves can be generated by the ARF resulting from the momentum transfer from the ultrasound field.12 In the experiment, continuous shear waves were generated by an oscillatory ARF via modulating the amplitude of ultrasound over time. Shear-wave frequencies were determined by the amplitude modulation frequency ranging from 100 to 600 Hz. Figure 4 is the schematic of the experimental setup. The green area is the photon probability density mapping of scattered light in the phantom predicted by Monte Carlo simulation13 using a million photons to get a good statistics of the photon probability density mapping, where a 6-mm radius detection area was defined to match the experiment. The parameters in the Monte Carlo simulation were determined based on Ref. 14 including optical absorption coefficient , optical scattering coefficient , and anisotropic coefficient . Shear waves were generated at P1 and P2 with amplitude-modulated ultrasound and tracked by the time-resolved contrast difference signal [e.g., Fig. 1(b)]. The attenuation coefficients of shear wave at each frequency were then calculated with Eq. (10). The tissue-mimicking phantoms used in the experiment were made of 4% gelatin and 1% agar following the recipe in Ref. 15. A 4% intralipid was added in the solution to obtain a reduced optical scattering coefficient of .14 In order to alter viscosity of the phantom, 0.3% Xanthan gum was added in one of the phantom (#1). The size of phantoms was , and both phantoms were made to be homogeneous. 3.ResultTable 1 shows the shear-wave attenuation coefficients measured in the experiment. Standard deviations were produced from three repeated measurements. As expected, the shear-wave attenuation dispersion is observed on both phantoms. Shear-wave attenuation increases with the increase of shear-wave frequency. The phantom with higher viscosity (#1) has a larger attenuation coefficients compared with the less viscous phantom (#2). Furthermore, according to Table 1, the standard deviation of attenuation coefficient measured with our method ranges from 2% to 23%. Table 1Attenuation coefficients over shear-wave frequency.
Equation (3) was used to calculate viscosity of phantoms with the shear-wave attenuation coefficients in Table 1 and a constrained elasticity range (within 6 kPa) estimated with the method in our previous article.9 Table 2 shows that the viscosities measured in the experiment (Row 2) agree with the reference values in Ref. 15 (Row 3) and have a maximum standard deviation of 9%. Table 2Comparison of our result with reference value.
4.Discussion and ConclusionThis study shows promising results of viscosity measurements in tissue-mimicking phantoms by SW-LASCA. The prior information on elasticity can be obtained by measuring the shear-wave speed described in Ref. 9. These two methods are based on the same setup; therefore, it is easy to switch between them. Because of the small wavelength of light, the sensitivity of shear-wave tracking with SW-LASCA could be better than ultrasound-based methods. In addition, optical detection inherently enables the detected signal to contain optical property information. In the future, we will combine the SW-LASCA and ultrasound-modulated optical tomography16 to build a dual-mode imaging system sensitive to both mechanical and optical contrasts. With the complementary information provided by optical detection, it could potentially be used as a diagnosis tool for cancerous tissues. The previously reported ultrasound-based method had a bias of on viscosity measurement.10,15 Our method shows a bias of . One possible reason is the variation of phantoms even though we followed the recipe used in previous reports. Furthermore, there are two additional factors. First, shear-wave reflection at phantom boundaries may affect the measured contrast difference values and, thus, the results in the bias on the attenuation coefficient estimation. In the experiment, the phantom was fixed within a water tank made of Perspex. Due to large shear impendence mismatch between the Perspex and the phantom, shear waves were almost totally reflected at the Perspex–phantom interface. Therefore, both the original and reflected shear waves contributed to the contrast difference signal, and this may bias the estimation of the shear-wave attenuation coefficient. However, this artificially strong shear-wave reflection in the experiments may not occur in real applications. Second, in Eq. (9), we assumed that the identical shear-wave amplitudes at detection point result in identical contrast differences, even if the shear waves originate from two positions along the propagation path. This assumption may not be totally true due to shear-wave diffraction: the shear wave with a longer propagation distance will overlap more with the scattered light and so generates a different contrast difference amplitude compared with the shear wave propagated a shorter distance. This effect could be minimized either by generating plane shear waves, so that less shear wave-spreading will be encountered, or by decreasing the distance between the two locations, so that the change of shear-wave geometry may be ignored. In this feasibility study, relatively high intensities of ultrasound were used in the experiments to generate shear waves. A mechanical index (MI) of 1.9 is usually recognized as the upper limit for clinical ultrasound use. Therefore, a MI of 1.28 was applied to prove the feasibility of reducing ultrasound intensity for the measurement. The same setup and system parameters were used to obtain Fig. 5, except for a lower ultrasound intensity. In Fig. 5, we used an ultrasound pressure equivalent to a MI of 1.28 at 5 MHz to generate shear waves, which is well below the ultrasound safety limit. Results show a good signal-to-noise ratio at . The standard deviation is produced from the three repeated measurements. Fig. 5Speckle contrast difference due to continuous shear wave generated by an ultrasound mechanical index (MI) of 1.28. Standard deviation is produced with three repeated measurements. ![]() In laser speckle contrast measurements, the CCD exposure time is important. On one hand, we hope to use a short enough CCD exposure time to minimize the speckle decorrelation effect due to blood flow and other tissue motions. On the other hand, we want to use a long enough CCD exposure time to maximize the light intensity received on the CCD in order to get a better signal-to-noise ratio. Generally, an exposure time less than 0.5 ms leads to a noisy signal, and an exposure time larger than 10 ms leads to 50% reduction of contrast in our measurements. Therefore, 2-ms exposure time is usually used in the experiments. However, this exposure time could be further reduced for in vivo measurements, where speckle decorrelation time is much shorter than 10 ms. In the future, it would also be interesting to evaluate to what extent the shear-wave signal will be affected by the speckle decorrelation by using a flow phantom. Phantoms in the study simulated the scattering properties of tissue with a thickness of 24 mm and a reduced optical scattering coefficient of about for the 532-nm laser.17 The reduced optical scattering coefficient was relatively weak compared with that of some biological tissues like breast.18 Therefore, further investigations are needed to assess the feasibility of our method on biological tissues. In addition, the green laser used in the current study will be replaced by a near-infrared laser for a better penetration depth in biological tissues. Currently, we used one single element transducer to generate shear waves at two positions by translating the transducer to different positions. In the future, we will program a linear array transducer to replace the “mechanical” translation to “electrical” translation by dynamic ultrasound beam focusing. Dynamic beam focusing will provide a faster and more flexible measurement. In summary, we proposed an optical method to quantitatively measure the viscosity at centimeter depth within tissue-mimicking phantoms. It is potentially more sensitive to viscosity changes than ultrasound-based methods, because the light wavelength is much smaller than that of the ultrasound. It may also complement viscoelastography techniques based on ultrasound or magnetic resonance imaging by providing optical information of tissues. By combining SW-LASCA with ultrasound-modulated optical tomography, we aim to build a dual-mode imaging system that is sensitive to both mechanical and optical contrasts. However, further investigation is needed to assess the feasibility of this method in vivo, where the speckle decorrelation time is short and the optical scattering coefficient is large. AcknowledgmentsThe authors want to thank Professor David O. Cosgrove for fruitful discussion, and EPSRC (EP/H02316X/1), the Royal Society, and China scholarship council for financial support. The correspondence should be made to M.X.T and D.S.E (mengxing.tang/daniel.elson@imperial.ac.uk). ReferencesD. Cosgroveet al.,
“Shear wave elastography for breast masses is highly reproducible,”
Eur. Radiol., 22
(5), 1023
–1032
(2012). http://dx.doi.org/10.1007/s00330-011-2340-y EURAE3 1432-1084 Google Scholar
N. Salamehet al.,
“Early detection of steatohepatitis in fatty rat liver by using MR elastography,”
Radiology, 253
(1), 90
–97
(2009). http://dx.doi.org/10.1148/radiol.2523081817 RADLAX 0033-8419 Google Scholar
J. Bamberet al.,
“EFSUMB guidelines and recommendations on the clinical use of ultrasound elastography. Part 1: basic principles and technology,”
Ultraschall Med., 34
(02), 169
–184
(2013). http://dx.doi.org/10.1055/s-00000089 0172-4614 Google Scholar
K. NightingaleS. McAleaveyG. Trahey,
“Shear-wave generation using acoustic radiation force: in vivo and ex vivo results,”
Ultrasound Med. Biol., 29
(12), 1715
–1723
(2003). http://dx.doi.org/10.1016/j.ultrasmedbio.2003.08.008 USMBA3 0301-5629 Google Scholar
M. L. PalmeriK. R. Nightingale,
“Acoustic radiation force-based elasticity imaging methods,”
Interface Focus, 1
(4), 553
–564
(2011). Google Scholar
C. Shigaoet al.,
“Shearwave dispersion ultrasound vibrometry (SDUV) for measuring tissue elasticity and viscosity,”
IEEE Trans. Ultrason. Ferroelectr. Freq. Control, 56
(1), 55
–62
(2009). http://dx.doi.org/10.1109/TUFFC.2009.1005 ITUCER 0885-3010 Google Scholar
A. E. Bohteet al.,
“MR elastography of the liver: defining thresholds for detecting viscoelastic changes,”
Radiology, 269
(3), 768
–776
(2013). RADLAX 0033-8419 Google Scholar
J. M. Hudsonet al.,
“Inter- and intra-operator reliability and repeatability of shear wave elastography in the liver: a study in healthy volunteers,”
Ultrasound Med. Biol., 39
(6), 950
–955
(2013). http://dx.doi.org/10.1016/j.ultrasmedbio.2012.12.011 USMBA3 0301-5629 Google Scholar
Y. Chenget al.,
“Shear wave elasticity imaging based on acoustic radiation force and optical detection,”
Ultrasound Med. Biol., 38
(9), 1637
–1645
(2012). http://dx.doi.org/10.1016/j.ultrasmedbio.2012.04.022 USMBA3 0301-5629 Google Scholar
S. ChenM. FatemiJ. F. Greenleaf,
“Quantifying elasticity and viscosity from measurement of shear wave speed dispersion,”
J. Acoust. Soc. Am., 115
(6), 2781
–2785
(2004). http://dx.doi.org/10.1121/1.1739480 JASMAN 0001-4966 Google Scholar
G. R. Torr,
“The acoustic radiation force,”
Am. J. Phys., 52
(5), 402
–408
(1984). http://dx.doi.org/10.1119/1.13625 AJPIAS 0002-9505 Google Scholar
V. G. Andreevet al.,
“Observation of shear waves excited by focused ultrasound in rubber-like medium,”
Acoust. Phys., 43
(2), 123
–128
(1997). AOUSEK 1063-7710 Google Scholar
L. WangS. L. JacquesL. Zheng,
“MCML—Monte Carlo modeling of light transport in multi-layered tissues,”
Comput. Meth. Prog. Biomed., 47
(2), 131
–146
(1995). http://dx.doi.org/10.1016/0169-2607(95)01640-F CMPBEK 0169-2607 Google Scholar
H. J. van Staverenet al.,
“Light scattering in Intralipid-10% in the wavelength range of 400–1100 nm,”
Appl. Opt., 30
(31), 4507
–4514
(1991). http://dx.doi.org/10.1364/AO.30.004507 APOPAI 0003-6935 Google Scholar
J. Bercoffet al.,
“The role of viscosity in the impulse diffraction field of elastic waves induced by the acoustic radiation force,”
IEEE Trans. Ultrason. Ferroelectr. Freq. Control, 51
(11), 1523
–1536
(2004). http://dx.doi.org/10.1109/TUFFC.2004.1367494 ITUCER 0885-3010 Google Scholar
R. Liet al.,
“Effects of acoustic radiation force and shear waves for absorption and stiffness sensing in ultrasound modulated optical tomography,”
Opt. Express, 19
(8), 7299
–7311
(2011). http://dx.doi.org/10.1364/OE.19.007299 OPEXFF 1094-4087 Google Scholar
K. DaoudiA.-C. BoccaraE. Bossy,
“Detection and discrimination of optical absorption and shear stiffness at depth in tissue-mimicking phantoms by transient optoelastography,”
Appl. Phys. Lett., 94
(15), 154103
(2009). http://dx.doi.org/10.1063/1.3120552 APPLAB 0003-6951 Google Scholar
T. Durduranet al.,
“Bulk optical properties of healthy female breast tissue,”
Phys. Med. Biol., 47
(16), 2847
–2861
(2002). http://dx.doi.org/10.1088/0031-9155/47/16/302 PHMBA7 0031-9155 Google Scholar
BiographyYi Cheng received her BEng degree from Tianjin University, China, with a specialization in opto-electronic engineering in 2010. In 2009 she visited University of Dundee, UK, and worked on the development of an ultrasonic needle for surgery. She is now achieving her PhD degree at Imperial College London, UK, working on optical shear wave elastography for early cancer detection and ultrasound modulated optical tomography for deep tissue imaging. She has developed a dual mode optical system that is sensitive on both tissue mechanical and optical contrast. Sinan Li joined Imperial College London as a PhD student in biomedical engineering in 2011. He is working on ultrasound mediated optical tomography, optical shear wave elastography, optical sensing of acoustic waves in turbid medium, and phase-transition contrast agents and their applications in ultrasound and optical imaging. He also carries out research on multiply scattering of light and laser speckle statistical analysis. He has developed a Monte Carlo simulation for laser speckle contrast detection of shear wave in an optically turbid medium. Daniel Elson is a reader in surgical imaging in the Hamlyn Centre for Robotic Surgery, Institute of Global Health Innovation and Department of Surgery and Cancer at Imperial College London and St. Mary's Hospital. His research interests are based around the development and application of photonics technology to medical imaging, including multispectral endoscopy, ultrasound mediated optical tomography, endoscopic structured lighting, light sources in endoscopy, and polarized endoscopic imaging. Mengxing Tang joined the Department of Bioengineering, Imperial College as a lecturer in 2006 and became a senior lecturer in 2011. Before moving to London, he was a departmental lecturer in the Department of Engineering Science, University of Oxford. He founded the Ultrasound Imaging Group at Imperial College in 2006 in close collaboration with Dr. Robert Eckersley from KCL. His main research interest is to develop new imaging, image analysis, and sensing techniques with ultrasound, microbubbles and light and apply them to a range of applications, including the diagnosis and treatment of cardiovascular diseases and cancer. He has authored 50 peer-reviewed journal papers. |