Research Papers: Imaging

Shear wave pulse compression for dynamic elastography using phase-sensitive optical coherence tomography

[+] Author Affiliations
Thu-Mai Nguyen, Bastien Arnal, Emily Y. Wong, Matthew O’Donnell

University of Washington, Department of Bioengineering, 3720 15th Avenue NE, Seattle, Washington 98195

Shaozhen Song

University of Washington, Department of Bioengineering, 3720 15th Avenue NE, Seattle, Washington 98195

University of Dundee, School of Engineering, Physics and Mathematics, Scotland, United Kingdom

Zhihong Huang

University of Dundee, School of Engineering, Physics and Mathematics, Scotland, United Kingdom

Ruikang K. Wang

University of Washington, Department of Bioengineering, 3720 15th Avenue NE, Seattle, Washington 98195

University of Washington, Department of Ophthalmology, 325 9th Avenue, Seattle, Washington 98104

J. Biomed. Opt. 19(1), 016013 (Jan 17, 2014). doi:10.1117/1.JBO.19.1.016013
History: Received October 17, 2013; Revised December 4, 2013; Accepted December 9, 2013
Text Size: A A A

Open Access Open Access

Abstract.  Assessing the biomechanical properties of soft tissue provides clinically valuable information to supplement conventional structural imaging. In the previous studies, we introduced a dynamic elastography technique based on phase-sensitive optical coherence tomography (PhS-OCT) to characterize submillimetric structures such as skin layers or ocular tissues. Here, we propose to implement a pulse compression technique for shear wave elastography. We performed shear wave pulse compression in tissue-mimicking phantoms. Using a mechanical actuator to generate broadband frequency-modulated vibrations (1 to 5 kHz), induced displacements were detected at an equivalent frame rate of 47 kHz using a PhS-OCT. The recorded signal was digitally compressed to a broadband pulse. Stiffness maps were then reconstructed from spatially localized estimates of the local shear wave speed. We demonstrate that a simple pulse compression scheme can increase shear wave detection signal-to-noise ratio (>12dB gain) and reduce artifacts in reconstructing stiffness maps of heterogeneous media.

Figures in this Article

Elastography is a recently developed field to evaluate the biomechanical properties of soft tissue for medical diagnosis. Stiffness variations have indeed proven to be good indicators of pathological state, providing valuable additional information to structural images in many applications: breast tumor detection,1 liver fibrosis diagnosis,2,3 and vascular pathologies characterization.4,5 Several techniques have been developed in the past 20 years, most of them involving magnetic resonance imaging (MRI)6 or ultrasound7,8 modalities. More recently, optical coherence elastography (OCE) has emerged9 based on the optical coherence tomography (OCT) imaging. OCT provides high spatial resolution (micron scale) and is, therefore, an ideal technology to examine small organs such as skin layers, ocular tissues (cornea, intraocular lens), and peripheral vessels. The earliest OCE techniques used imaging of tissue displacements following a static deformation.10,11 Static elastography, however, suffers from the following limitations: operator dependence (variability of the applied deformation) and nonquantitative estimates of the stiffness due to high sensitivity to boundary conditions (uneven stress distribution).

Phase-sensitive OCT (PhS-OCT)1215 is one of the latest improvements enabling displacement detection with much higher sensitivity (nanometer scale) than the previous optical intensity-based techniques. Dynamic elastography (i.e., based on dynamic excitations1618 and later on transient elastography methods1921) have been implemented using PhS-OCT. Transient elastography has been introduced to address the abovementioned limitations of static elastography. In the transient regime, propagating shear waves are induced in the sample and their propagation is tracked with PhS-OCT. The stiffness of the sample is then characterized by retrieving the shear modulus from the shear wave propagation speed. PhS-OCT has also been combined with multifrequency excitations to investigate the frequency dependence of the mechanical response of viscoelastic media.22

In this work, we propose to combine dynamic elastography with a digital pulse compression method to improve shear modulus estimation. Pulse compression was first introduced for radar detection in the early 1960s23,24 and later extended to ultrasound applications.25 The main idea is to spread the instantaneous peak energy level over a longer period while simultaneously improving the axial resolution and increasing the signal-to-noise ratio (SNR) through pulse compression. Here, we propose to implement pulse compression with shear waves to: (1) improve the SNR of detected displacements induced in the sample and (2) compress spatially and temporally the shear wave for better reconstruction of shear modulus maps. For that purpose, we propose to generate a broadband frequency-modulated (chirp) mechanical excitation and apply digital pulse compression to the resulting displacement field detected with PhS-OCT. We evaluate the performance of this pulse compression scheme by comparing the results to those from monochromatic excitation as employed in our previous studies.26

We performed shear wave elastography experiments on agar phantoms. Our setup consists of a PhS-OCT imaging system combined with a mechanical shear source.

Sample Preparation

Agar phantoms were made to mimic the elastic properties of biological soft tissues. Two types of tissue-mimicking phantoms were constructed using the same methods detailed in the previous studies:26

  • A homogeneous phantom was made from a 0.5%-agar (w/v) solution. This simple medium was used to assess the performance of our pulse compression algorithm.
  • A heterogeneous phantom was made from a 0.5%-agar background and a 1%-agar, 1-mm diameter cylindrical inclusion. This phantom constitutes a crude model of a diseased tissue (e.g., tumor) that we used to demonstrate the benefits of pulse compression for elastic modulus reconstruction in complex media.

A few drops of 10% latex microspheres suspension (0.3-μm diameter, Duke Scientific Corporation, California) were added to both phantoms as optical scatterers. Such scatterers were used in the previous studies27 and found to provide adequate backscattering to reach imaging ranges up to 2-mm deep.

Shear Wave Generation

Shear waves were induced in the sample using a piezoelectric actuator (AE0505, Thorlabs, New Jersey) driven by a function generator (Tektronix, Oregon) amplified by a power amplifier with controllable gain (AE Techron, Indiana). Typical voltages of 80 V were used to induce 8μm amplitude displacements of the actuator. For the purpose of this study, we used a linear frequency modulation to generate an 8-ms long chirp with a frequency content ranging from 1 to 5 kHz as the driving signal for the actuator. In addition, to compare with the previous techniques, we also repeated all the measurements using an 8-ms toneburst at a single frequency of 3 kHz. This excitation method is similar to that used in our previous work26 and serves as a reference to evaluate the performances of the pulse compression method.

The coupling between the actuator and the sample was done through a polished stainless steel tip attached to the actuator. The tip was wedge-like in shape, with a minimal width of 1.7 mm. The waves generated by this tip can therefore be considered approximately planar. The tip was positioned in such a manner that it touches the sample side without deforming it (Fig. 1). Accurate positioning was performed using a precision translation stage and monitored using real-time OCT B-mode images. The vertical vibration of the tip generates displacements that are mainly polarized axially (i.e., along the OCT probing beam axis) and propagate laterally across the sample, creating a propagating shear wave. The tip was placed on the side of the sample so that significant displacements were induced across the whole sample depth, generating bulk shear waves.

Graphic Jump LocationF1 :

Setup for dynamic elastography combining a phase-sensitive optical coherence tomography system for detection and a mechanical shear source. (x,z) are the lateral (perpendicular to the incident light beam) and axial (along the incident light beam) coordinates, respectively.

Such a shear wave has a propagation speed that depends on the sample elastic moduli. For a locally homogeneous, incompressible, and isotropic medium that does not exhibit strong shear viscosity, the shear wave speed cS is28Display Formula

cS=μ/ρ,(1)
where μ is the shear modulus of the sample and ρ is its density.

PhS-OCT Imaging System

Figure 1 illustrates the experimental setup. Our PhS-OCT system is a spectral-domain OCT setup, based on a fiber Michelson interferometer. The light source is a super-luminescent diode [1310-nm central wavelength, 56-nm spectral bandwidth (full-width half-maximum), DenseLight Semiconductors Ltd, Singapore]. Light in the reference arm is delivered to a stationary mirror while that in the sample arm is focused into a tissue-mimicking phantom. Backscattered light from both arms is recombined through an optical coupler. The resulting spectral interferogram is detected by a line-scan camera (14 bit, 1024 pixel, InGaAs detector, Sensors Ltd, New Jersey) capable of 47,000 A-scans per second. The system has a measured dynamic range of 105dB at 0.5-mm depth and a measured phase noise of 3.0 mrad.

Shear waves were detected by operating the PhS-OCT system in M-B mode, as detailed in our previous publications.26 In short, each M-scan is a set of 512 successive acquisitions performed at one given lateral location at a 47 kHz repetition rate, yielding a 10ms recording time. The incident light beam is then swept across the sample to perform a B-scan while the mechanical excitation is repeated. Each B-scan consists of 128 A-lines with a 23.4 μm spacing, covering a lateral imaging range of 2.5 mm. The axial pixel size is 5μm. A full M-B scan therefore consists of 512×128 A-lines and is acquired in 1.4s. The actuator excitation was triggered at the beginning of each M-scan using custom LabVIEW code.

Each M-scan was used to retrieve the displacements induced by shear wave propagation. Phase differences Δφ between two successive A-lines yield the local axial displacement uz(x,z,t) according to Display Formula

uz(x,z,t)=Δφ(x,z,t)4πnλ,(2)
where (x,z) are the spatial coordinates of a given pixel, t is a given time, λ=1310nm is the central wavelength of the light source, and n1.33 is the refractive index of the sample. Thus, the axial displacement field was computed as a function of time for each pixel of the region of interest, providing a “movie” of shear wave propagation (512 frames at 47 kHz).

Postprocessing
Pulse compression

Digital pulse compression is theoretically achieved by performing the autocorrelation of the displacement field. In practice, we used an inverse filter method to optimize the autocorrelation, mainly by reducing side lobes of the autocorrelation function. This approach closely parallels that discussed detailed in 25 for the conventional ultrasound imaging, where the filter was designed using the following steps:

  • A reference signal uzref(x=0,t) is computed by averaging along the depth of the axial displacements at x=0 (nearest point to the shear source).
  • An ideal output pulse pulse(τ) is derived from the autocorrelation of the reference signal Display Formula
    pulse(τ)=[uzref(t)uzref(t)]*w(τ),(3)
    where w(τ) is a weighting function designed to suppress side lobes of the autocorrelation function, and * denote respectively the convolution and the multiplication.
  • An inverse filter f(τ) is then designed by inverting Display Formula
    [uzref(t)uzref(t)]f(τ)=pulse(τ)Uzref(τ,τ)·f(τ)=pulse(τ),(4)
    where Uzref(τ,τ) is the convolution matrix of the reference signal. The inverse filter is then applied (i.e., convolved) to the autocorrelation of the axial displacement field at each location in the region of interest.

Shear modulus reconstruction

The local shear wave speed was determined by applying a custom time-of-flight algorithm to the axial displacements after pulse compression. The details on the principle of time-of-flight estimation for shear wave elastography can be found in Refs. 1 and 29. Briefly, correlation along the time of two adjacent A-lines of the axial displacement field provides the travel time of the shear wave between these two adjacent locations, yielding the local shear wave speed (group velocity) at each location within the region of interest. Maps of the local shear modulus are then retrieved using Eq. (1).

A “time-gating” of the displacement field can be applied prior to the time-of-flight algorithm to reduce the effects of reflected shear waves on the reconstruction in heterogeneous media. Time-gating selects only the earliest incident shear wavefront that is less likely to be affected by reflections (on elastic heterogeneities) occurring at later times (i.e., “near-ballistic propagation”). For that purpose, we developed an algorithm to detect the earliest arriving shear wave (at each location of the imaging plane) and applied a Gaussian temporal window to eliminate later wavefronts.

Homogeneous Phantom

Figure 2 shows the broadband displacement field detected in a 0.5%-agar phantom resulting from an 8-ms long chirp excitation. As illustrated in Fig. 2(a), the shear wave propagates laterally (from left to right) with lower frequencies (i.e., larger wavelengths) at early times and higher frequencies (i.e., smaller wavelengths) at later times. The SNR was computed from the spectrum of the displacements at one location [Fig. 2(b)] as the ratio between the energy contained in the [1 to 5] kHz range and the energy out of that range. For this example, the SNR was estimated to be 10.2 dB.

Graphic Jump LocationF2 :

Displacement field in a 0.5%-agar phantom resulting from a chirp excitation. (a) Axial displacements (color scale) superimposed on the B-mode image (gray scale) at different instants. (b) Normalized axial displacements profile (upper) and spectrum (lower) at one location.

The distortion of the displacement field at the phantom surface is probably due to interferences between the surface and bulk shear waves, which propagate with slightly different speeds. Thus, the surface region was not considered for the pulse compression.

Figures 3(a) and 3(b) show, respectively, the ideal output pulse [as defined in Eq. (3)] and the inverse filter [as defined in Eq. (4)] derived from the displacement field shown in Fig. 2. The result of that inverse filter is shown in Fig. 4.

Graphic Jump LocationF3 :

Pulse compression applied to a 0.5%-agar phantom. (a) Ideal output pulse (red thick line) derived from the autocorrelation function (black line) as defined in Eq. (3). (b) Inverse filter obtained from Eq. (4).

Graphic Jump LocationF4 :

Displacement field in a 0.5%-agar phantom after compression. (a) Normalized axial displacements (color scale) superimposed on the B-mode image (gray scale) at different instants. The white dash line delineates the region of interest where the compression has been applied. (b) Axial displacements profile (upper) and spectrum (lower) at one location.

Figure 4 shows the result of pulse compression obtained by applying the filter shown in Fig. 3(b) to the data shown in Fig. 2. The initial displacement field has been spatially and temporally compressed to a short broadband pulse of 1-ms duration, using the procedure described in Sec. 2.4.1. Figure 4(b) is an example of the axial displacements at one given location. The SNR was estimated to be 23.0 dB from the spectrum shown in Fig. 4(b), which corresponds to a 12.8 dB increase compared to the initial chirp.

The compressed displacement field at a given depth can also be represented in the (t,x) domain as displayed in Fig. 5(a). In the (t,x) domain, the local slope of the signal represents the local shear wave propagation speed. The local shear wave speed values are converted to the local shear modulus values, leading to the shear modulus map shown in Fig. 5(b). For this phantom, the median shear modulus was estimated to be 4.14±0.63kPa.

Graphic Jump LocationF5 :

(a) Spatio-temporal representation of the normalized axial displacement field at a given depth of a 0.5%-agar phantom. The x-axis represents the time whereas the y-axis represents the lateral position. (b) Reconstructed shear modulus map (color scale) superimposed on the B-mode image (gray scale).

Heterogeneous Phantom

We performed experiments on a heterogeneous phantom to investigate the advantages of pulse compression in a more complex medium. Figure 6 shows a B-mode image of the phantom, made from a 0.5%-agar background and a 1%-agar inclusion. This structural image exhibits very low contrast between the inclusion and the background.

Graphic Jump LocationF6 :

B-mode image of a heterogeneous phantom made from 0.5%-agar background and 1%-agar inclusion. The white dash line delineates the inclusion.

We first performed measurements using a 3-kHz single-frequency excitation for the comparison with our previous work.26Figure 7(a) shows the displacement field detected at a given depth in the (t,x) domain. The distortions of the shear wavefronts can be observed on most cycles. They are due to interference between the incident shear wave and reflections of the shear wave at the boundaries of the inclusion. Multiwave superposition due to reflections can significantly affect the reconstruction, as shown in Fig. 7(c): the corresponding shear modulus map exhibits strong artifacts in the inclusion (unexpected soft region inside the inclusion and soft stripe near the right edge) as well as downstream of the inclusion (soft stripes on the right side of the region of interest). The median shear modulus of the expected inclusion region and the background estimated from this map are 6.86±4.26kPa and 4.86±1.53kPa, respectively.

Graphic Jump LocationF7 :

Results using a 3 kHz, 8-ms long, excitation on a 0.5%-agar phantom containing a 1%-agar inclusion. (a and b) Spatio-temporal representation of the normalized axial displacement field without (a) and with (b) time-gating. (c and d) Shear modulus map (color scale) superimposed on the B-mode image (gray scale) reconstructed without (c) and with (d) time-gating.

The effects of the reflected shear waves can be reduced using a time-gating method, as described in Sec. 2.4.2. The resulting time-gated displacement field is shown in Fig. 7(b). The corresponding shear modulus map is shown in Fig. 7(d). Only minor artifacts have been suppressed and the primary ones remain. This suggests that the first wavefront has too large a period (0.33 ms) to avoid interfering with the reflected waves. The median shear modulus of the expected inclusion region and the background estimated from this map are 7.20±5.11kPa and 4.51±1.30kPa, respectively.

Figure 8 summarizes the results obtained on the heterogeneous phantom using the pulse compression method. Figure 8(a) shows the displacements detected at a given depth resulting from the chirp excitation. Note that only the first 6 ms of the 8-ms long signal is displayed here. Similar to the single-frequency excitation, interferences between the incident waves and the reflected waves can be observed on the later cycles of the chirp. Figure 8(b) shows the displacement field after pulse compression. The initial chirp has been compressed to a main broadband pulse of about 0.21-ms width. Again, side lobes and distortions appear on the later cycles as a consequence of shear wave reflections at the boundaries of the inclusion. However, time-gating can be applied to isolate the earliest wavefront [Fig. 8(c)]. The corresponding shear modulus map is displayed in Fig. 8(d) with the exact same color scale as that of Figs. 7(c) and 7(d). Artifacts within the inclusion have been eliminated, allowing clear identification of the inclusion as a stiffer region than the background with high elastic contrast. The median shear modulus of the inclusion and the background are 14.81±5.12kPa and 5.00±0.79kPa, respectively.

Graphic Jump LocationF8 :

Heterogeneous phantom consisting of a 1%-agar background and a 2%-agar inclusion. (a) Spatio-temporal representation of the raw normalized displacement field. (b and c) Spatio-temporal representation of the displacement field after compression, without (b) and with (c) time-gating. (d) Shear modulus map reconstructed from the compressed, time-gated displacement field.

We have implemented digital pulse compression with an OCT-based shear wave elastography technique. We have performed the experiments on tissue-mimicking phantoms. A broadband ([1 to 5] kHz) frequency-modulated mechanical excitation was applied to phantoms through a piezoelectric actuator. The subsequent displacement field was detected and recorded using a PhS-OCT system. The inverse filtering was performed on the recorded displacements to spatio-temporally compress the initial propagating chirp into a broadband propagating pulse. The pulse propagation was then used to either reconstruct maps of the local shear modulus or compute dispersion curves (see below).

We have demonstrated that pulse compression significantly improves the SNR of detected displacements: a 12.8 dB increase has been obtained in a homogeneous agar phantom. PhS-OCT already has a great sensitivity (nanometer scale) to detect axial displacements compared to other imaging modalities such as ultrasound or MRI. Increasing the SNR will still be valuable for in vivo applications and particularly for ocular tissues. For instance, the cornea is both a stiff and fragile tissue that cannot undergo large displacements.

We have also demonstrated that pulse compression is an efficient method to simplify the reconstruction of shear modulus maps in heterogeneous media. The incident shear wave can indeed be reflected by the boundaries of elastic heterogeneities. For the case of nonlocalized (spatial and temporal) excitations, reflected waves have a high probability of interfering with the incident wave, leading to artifacts in straightforward time-of-flight reconstruction. In contrast, pulse compression enables time-gating on a localized (spatial and temporal) signal to isolate the incident waves from reflected waves. For instance, for the heterogeneous phantom presented in this study, the expected shear modulus values are, respectively, about 4 and 15 kPa for the 0.5%-agar background and the 1%-agar inclusion according to the values reported in the literature.30 The shear modulus of the inclusion estimated from the nonlocalized excitation (Fig. 7) is strongly underestimated (7.20 kPa) and has a high spatial relative variance (71% of the shear modulus), whereas that estimated from the pulse compression approach (Fig. 8) is much more accurate (14.81 kPa) and has a lower spatial relative variance (34% of the shear modulus). The shear modulus contrast between the inclusion and the background is much higher using the pulse compression method (50%) than using single-frequency excitation (23%).

Other refinements in postprocessing can be added to time-gating to improve the reconstruction. The shear modulus map shown in Fig. 8 (reconstructed after pulse compression and time-gating) still exhibits heterogeneities downstream of the inclusion. These apparent heterogeneities are likely caused by diffraction of the shear wave by the inclusion. In the previous work on elastographic reconstruction for shear waves, directional filters have been applied to displacement fields to reduce the effects of multidirectional propagation due to diffraction.31,32 Using a two-dimensional Fourier representation of displacement fields, plane waves propagating in different directions occupy different regions of Fourier space. Consequently, directional filters can be applied to eliminate propagation directions that are highly angled compared to the incident one. Even without these more sophisticated reconstruction approaches, we have proven the impact of pulse compression. Optimizing reconstructions with more sophisticated algorithms will be considered in future work.

An additional feature of pulse compression is the possibility to perform shear wave spectroscopy with single-shot acquisition. Frequency-dependent information such as the phase velocity (i.e., propagation speed of each frequency contained in the wave) and the attenuation can indeed be extracted from a broadband shear wave. The details can be found in Refs. 33 and 34. Although the samples used in this study did not exhibit any significant dispersion over the frequency range of the applied chirp, this feature will be useful for investigating biological tissue. For instance, viscosity can be extracted from the dispersion curves. Besides, thin-layered tissues (e.g., the cornea, vessels walls, skin layers) constrain shear wave propagation to guided modes that induce strong dispersive effects.5,34 Accessing the dispersion will thus be highly valuable, especially for OCT-based elastography, since OCT tends to examine small tissues due to its high spatial resolution.

Finally, we used a setup combining a mechanical actuator and PhS-OCT similar to our previous studies.26 The feasibility of performing elastography measurements in vivo with such a setup has been demonstrated in our previous studies on human skin.20 Nevertheless, pulse compression can be extended to a wide range of dynamic elastography methods using different shear sources and different detection modalities.

This work was supported in part by NIH RO1EB016034, R01CA170734, R01EB009682, R01HL093140, R01DC010201, Life Sciences Discovery Fund 3292512, the Coulter Translational Research Partnership Program, and the Department of Bioengineering at the University of Washington.

Tanter  M. et al., “Quantitative assessment of breast lesion viscoelasticity: initial clinical results using supersonic shear imaging,” Ultrasound Med. Biol.. 34, (9 ), 1373 –1386 (2008). 0301-5629 CrossRef
Palmeri  M. L. et al., “Quantifying hepatic shear modulus in vivo using acoustic radiation force,” Ultrasound Med. Biol.. 34, (4 ), 546 –558 (2008). 0301-5629 CrossRef
Bavu  E. et al., “Noninvasive in vivo liver fibrosis evaluation using supersonic shear imaging: a clinical study on 113 hepatitis C virus patients,” Ultrasound Med. Biol.. 37, (9 ), 1361 –1373 (2011). 0301-5629 CrossRef
de Korte  C. L. et al., “Characterization of plaque components with intravascular ultrasound elastography in human femoral and coronary arteries in vitro,” Circulation. 102, (6 ), 617 –623 (2000). 0009-7322 CrossRef
Couade  M. et al., “Quantitative assessment of arterial wall biomechanical properties using shear wave imaging,” Ultrasound Med. Biol.. 36, (10 ), 1662 –1676 (2010). 0301-5629 CrossRef
Muthupillai  R., Ehman  R. L., “Magnetic resonance elastography,” Nat. Med.. 2, (5 ), 601 –603 (1996). 1078-8956 CrossRef
Ophir  J. et al., “Elastography a quantitative method for imaging the elasticity of biological tissues,” Ultrason. Imaging. 13, , 111 –134 (1991). 0161-7346 CrossRef
O’Donnell  M. et al., “Internal displacement and strain imaging using ultrasonic speckle tracking,” IEEE Trans. Ultrason. Ferroelectr. Freq. Control. 41, (3 ), 314 –325 (1994). 0885-3010 CrossRef
Sun  C., Standish  B., Yang  V. X. D., “Optical coherence elastography: current status and future applications,” J. Biomed. Opt.. 16, (4 ), 043001  (2011). 1083-3668 CrossRef
Rogowska  J. et al., “Optical coherence tomographic elastography technique for measuring deformation and strain of atherosclerotic tissues,” Heart. 90, (5 ), 556 –562 (2004). 1355-6037 CrossRef
Ko  H.-J. et al., “Optical coherence elastography of engineered and developing tissue,” Tissue Eng.. 12, (1 ), 63 –73 (2006). 1937-3341 CrossRef
Wang  R. K., Ma  Z., Kirkpatrick  S. J., “Tissue Doppler optical coherence elastography for real time strain rate and strain mapping of soft tissue,” Appl. Phys. Lett.. 89, (14 ), 144103  (2006). 0003-6951 CrossRef
Wang  R. K., Kirkpatrick  S., Hinds  M., “Phase-sensitive optical coherence elastography for mapping tissue microstrains in real time,” Appl. Phys. Lett.. 90, (16 ), 164105  (2007). 0003-6951 CrossRef
Kennedy  B. F. et al., “Strain estimation in phase-sensitive optical coherence elastography,” Biomed. Opt. Express. 3, (8 ), 1865 –1879 (2012). 2156-7085 CrossRef
Wang  R. K., Nuttall  A. L., “Phase-sensitive optical coherence tomography imaging of the tissue motion within the organ of Corti at a subnanometer scale: a preliminary study,” J. Biomed. Opt.. 15, (5 ), 056005  (2010). 1083-3668 CrossRef
Adie  S. G. et al., “Audio frequency in vivo optical coherence elastography,” Phys. Med. Biol.. 54, (10 ), 3129 –3139 (2009). 0031-9155 CrossRef
Kennedy  B. F. et al., “In vivo dynamic optical coherence elastography using a ring actuator,” Opt. Express. 17, (24 ), 21762 –21772 (2009). 1094-4087 CrossRef
Liang  X. et al., “Dynamic spectral-domain optical coherence elastography for tissue characterization,” Opt. Express. 18, (13 ), 14183 –14190 (2010). 1094-4087 CrossRef
Manapuram  R. K. et al., “In vivo estimation of elastic wave parameters using phase-stabilized swept source optical coherence elastography,” J. Biomed. Opt.. 17, (10 ), 100501  (2012). 1083-3668 CrossRef
Li  C. et al., “Determining elastic properties of skin by measuring surface waves from an impulse mechanical stimulus using phase-sensitive optical coherence tomography,” J. R. Soc. Interface. 9, (70 ), 831 –841 (2012). 1742-5689 CrossRef
Li  C. et al., “Noncontact all-optical measurement of corneal elasticity,” Opt. Lett.. 37, (10 ), 1625 –1627 (2012). 0146-9592 CrossRef
Adie  S. G. et al., “Spectroscopic optical coherence elastography,” Opt. Express. 18, (25 ), 25519 –25534 (2010). 1094-4087 CrossRef
Cook  C. E., “Pulse compression-key to more efficient radar transmission,” Proc. IRE. 48, (3 ), 310 –316 (1959). 0096-8390 CrossRef
Kroszczynski  J. J., “Pulse compression by means of linear-period modulation,” Proc. IEEE. 57, (7 ), 1260 –1266 (1969). 0018-9219 CrossRef
O’Donnell  M., “Coded excitation system for improving the penetration of real-time phased-array imaging systems,” IEEE Trans. Ultrason. Ferroelectr. Freq. Control. 39, (3 ), 341 –351 (1992). 0885-3010 CrossRef
Song  S. et al., “Shear modulus imaging by direct visualization of propagating shear waves with phase-sensitive optical coherence tomography,” J. Biomed. Opt.. 18, (12 ), 121509  (2013). 1083-3668 CrossRef
Jia  Y. et al., “Doppler optical coherence tomography imaging of local fluid flow and shear stress within microporous scaffolds,” J. Biomed. Opt.. 14, (3 ), 034014  (2009). 1083-3668 CrossRef
Royer  D., Dieulesaint  E., Elastic Waves in Solids I: Free and Guided Propagation. ,  Springer-Verlag ,  Berlin, Heidelberg, New York  (2000).
McLaughlin  J., Renzi  D., “Shear wave speed recovery in transient elastography and supersonic imaging using propagating fronts,” Inverse Probl.. 22, (2 ), 681 –706 (2006). 0266-5611 CrossRef
Hamhaber  U. et al., “Comparison of quantitative shear wave MR-elastography with mechanical compression tests,” Magn. Reson. Med.. 49, (1 ), 71 –77 (2003). 0740-3194 CrossRef
Deffieux  T. et al., “On the effects of reflected waves in transient shear wave elastography,” IEEE Trans. Ultrason. Ferroelectr. Freq. Control. 58, (10 ), 2032 –2035 (2011). 0885-3010 CrossRef
Song  P. et al., “Fast shear compounding using directional filtering and two-dimensional shear wave speed calculation,” in  Proc. of the IEEE Joint UFFC, EFTF, PRF Symposium ,  Prague, Czech Republic , Vol. 1 (2013).
Deffieux  T. et al., “Shear wave spectroscopy for in vivo quantification of human soft tissues visco-elasticity,” IEEE Trans. Med. Imaging. 28, (3 ), 313 –322 (2009). 0278-0062 CrossRef
Nguyen  T.-M. et al., “Assessment of viscous and elastic properties of sub-wavelength layered soft tissues using shear wave spectroscopy: theoretical framework and in vitro experimental validation,” IEEE Trans. Ultrason. Ferroelectr. Freq. Control. 58, (11 ), 2305 –2315 (2011). 0885-3010 CrossRef

Biographies of the authors are not available.

© The Authors. Published by SPIE under a Creative Commons Attribution 3.0 Unported License. Distribution or reproduction of this work in whole or in part requires full attribution of the original publication, including its DOI.

Citation

Thu-Mai Nguyen ; Shaozhen Song ; Bastien Arnal ; Emily Y. Wong ; Zhihong Huang, et al.
"Shear wave pulse compression for dynamic elastography using phase-sensitive optical coherence tomography", J. Biomed. Opt. 19(1), 016013 (Jan 17, 2014). ; http://dx.doi.org/10.1117/1.JBO.19.1.016013


Figures

Graphic Jump LocationF1 :

Setup for dynamic elastography combining a phase-sensitive optical coherence tomography system for detection and a mechanical shear source. (x,z) are the lateral (perpendicular to the incident light beam) and axial (along the incident light beam) coordinates, respectively.

Graphic Jump LocationF2 :

Displacement field in a 0.5%-agar phantom resulting from a chirp excitation. (a) Axial displacements (color scale) superimposed on the B-mode image (gray scale) at different instants. (b) Normalized axial displacements profile (upper) and spectrum (lower) at one location.

Graphic Jump LocationF3 :

Pulse compression applied to a 0.5%-agar phantom. (a) Ideal output pulse (red thick line) derived from the autocorrelation function (black line) as defined in Eq. (3). (b) Inverse filter obtained from Eq. (4).

Graphic Jump LocationF4 :

Displacement field in a 0.5%-agar phantom after compression. (a) Normalized axial displacements (color scale) superimposed on the B-mode image (gray scale) at different instants. The white dash line delineates the region of interest where the compression has been applied. (b) Axial displacements profile (upper) and spectrum (lower) at one location.

Graphic Jump LocationF5 :

(a) Spatio-temporal representation of the normalized axial displacement field at a given depth of a 0.5%-agar phantom. The x-axis represents the time whereas the y-axis represents the lateral position. (b) Reconstructed shear modulus map (color scale) superimposed on the B-mode image (gray scale).

Graphic Jump LocationF6 :

B-mode image of a heterogeneous phantom made from 0.5%-agar background and 1%-agar inclusion. The white dash line delineates the inclusion.

Graphic Jump LocationF7 :

Results using a 3 kHz, 8-ms long, excitation on a 0.5%-agar phantom containing a 1%-agar inclusion. (a and b) Spatio-temporal representation of the normalized axial displacement field without (a) and with (b) time-gating. (c and d) Shear modulus map (color scale) superimposed on the B-mode image (gray scale) reconstructed without (c) and with (d) time-gating.

Graphic Jump LocationF8 :

Heterogeneous phantom consisting of a 1%-agar background and a 2%-agar inclusion. (a) Spatio-temporal representation of the raw normalized displacement field. (b and c) Spatio-temporal representation of the displacement field after compression, without (b) and with (c) time-gating. (d) Shear modulus map reconstructed from the compressed, time-gated displacement field.

Tables

References

Tanter  M. et al., “Quantitative assessment of breast lesion viscoelasticity: initial clinical results using supersonic shear imaging,” Ultrasound Med. Biol.. 34, (9 ), 1373 –1386 (2008). 0301-5629 CrossRef
Palmeri  M. L. et al., “Quantifying hepatic shear modulus in vivo using acoustic radiation force,” Ultrasound Med. Biol.. 34, (4 ), 546 –558 (2008). 0301-5629 CrossRef
Bavu  E. et al., “Noninvasive in vivo liver fibrosis evaluation using supersonic shear imaging: a clinical study on 113 hepatitis C virus patients,” Ultrasound Med. Biol.. 37, (9 ), 1361 –1373 (2011). 0301-5629 CrossRef
de Korte  C. L. et al., “Characterization of plaque components with intravascular ultrasound elastography in human femoral and coronary arteries in vitro,” Circulation. 102, (6 ), 617 –623 (2000). 0009-7322 CrossRef
Couade  M. et al., “Quantitative assessment of arterial wall biomechanical properties using shear wave imaging,” Ultrasound Med. Biol.. 36, (10 ), 1662 –1676 (2010). 0301-5629 CrossRef
Muthupillai  R., Ehman  R. L., “Magnetic resonance elastography,” Nat. Med.. 2, (5 ), 601 –603 (1996). 1078-8956 CrossRef
Ophir  J. et al., “Elastography a quantitative method for imaging the elasticity of biological tissues,” Ultrason. Imaging. 13, , 111 –134 (1991). 0161-7346 CrossRef
O’Donnell  M. et al., “Internal displacement and strain imaging using ultrasonic speckle tracking,” IEEE Trans. Ultrason. Ferroelectr. Freq. Control. 41, (3 ), 314 –325 (1994). 0885-3010 CrossRef
Sun  C., Standish  B., Yang  V. X. D., “Optical coherence elastography: current status and future applications,” J. Biomed. Opt.. 16, (4 ), 043001  (2011). 1083-3668 CrossRef
Rogowska  J. et al., “Optical coherence tomographic elastography technique for measuring deformation and strain of atherosclerotic tissues,” Heart. 90, (5 ), 556 –562 (2004). 1355-6037 CrossRef
Ko  H.-J. et al., “Optical coherence elastography of engineered and developing tissue,” Tissue Eng.. 12, (1 ), 63 –73 (2006). 1937-3341 CrossRef
Wang  R. K., Ma  Z., Kirkpatrick  S. J., “Tissue Doppler optical coherence elastography for real time strain rate and strain mapping of soft tissue,” Appl. Phys. Lett.. 89, (14 ), 144103  (2006). 0003-6951 CrossRef
Wang  R. K., Kirkpatrick  S., Hinds  M., “Phase-sensitive optical coherence elastography for mapping tissue microstrains in real time,” Appl. Phys. Lett.. 90, (16 ), 164105  (2007). 0003-6951 CrossRef
Kennedy  B. F. et al., “Strain estimation in phase-sensitive optical coherence elastography,” Biomed. Opt. Express. 3, (8 ), 1865 –1879 (2012). 2156-7085 CrossRef
Wang  R. K., Nuttall  A. L., “Phase-sensitive optical coherence tomography imaging of the tissue motion within the organ of Corti at a subnanometer scale: a preliminary study,” J. Biomed. Opt.. 15, (5 ), 056005  (2010). 1083-3668 CrossRef
Adie  S. G. et al., “Audio frequency in vivo optical coherence elastography,” Phys. Med. Biol.. 54, (10 ), 3129 –3139 (2009). 0031-9155 CrossRef
Kennedy  B. F. et al., “In vivo dynamic optical coherence elastography using a ring actuator,” Opt. Express. 17, (24 ), 21762 –21772 (2009). 1094-4087 CrossRef
Liang  X. et al., “Dynamic spectral-domain optical coherence elastography for tissue characterization,” Opt. Express. 18, (13 ), 14183 –14190 (2010). 1094-4087 CrossRef
Manapuram  R. K. et al., “In vivo estimation of elastic wave parameters using phase-stabilized swept source optical coherence elastography,” J. Biomed. Opt.. 17, (10 ), 100501  (2012). 1083-3668 CrossRef
Li  C. et al., “Determining elastic properties of skin by measuring surface waves from an impulse mechanical stimulus using phase-sensitive optical coherence tomography,” J. R. Soc. Interface. 9, (70 ), 831 –841 (2012). 1742-5689 CrossRef
Li  C. et al., “Noncontact all-optical measurement of corneal elasticity,” Opt. Lett.. 37, (10 ), 1625 –1627 (2012). 0146-9592 CrossRef
Adie  S. G. et al., “Spectroscopic optical coherence elastography,” Opt. Express. 18, (25 ), 25519 –25534 (2010). 1094-4087 CrossRef
Cook  C. E., “Pulse compression-key to more efficient radar transmission,” Proc. IRE. 48, (3 ), 310 –316 (1959). 0096-8390 CrossRef
Kroszczynski  J. J., “Pulse compression by means of linear-period modulation,” Proc. IEEE. 57, (7 ), 1260 –1266 (1969). 0018-9219 CrossRef
O’Donnell  M., “Coded excitation system for improving the penetration of real-time phased-array imaging systems,” IEEE Trans. Ultrason. Ferroelectr. Freq. Control. 39, (3 ), 341 –351 (1992). 0885-3010 CrossRef
Song  S. et al., “Shear modulus imaging by direct visualization of propagating shear waves with phase-sensitive optical coherence tomography,” J. Biomed. Opt.. 18, (12 ), 121509  (2013). 1083-3668 CrossRef
Jia  Y. et al., “Doppler optical coherence tomography imaging of local fluid flow and shear stress within microporous scaffolds,” J. Biomed. Opt.. 14, (3 ), 034014  (2009). 1083-3668 CrossRef
Royer  D., Dieulesaint  E., Elastic Waves in Solids I: Free and Guided Propagation. ,  Springer-Verlag ,  Berlin, Heidelberg, New York  (2000).
McLaughlin  J., Renzi  D., “Shear wave speed recovery in transient elastography and supersonic imaging using propagating fronts,” Inverse Probl.. 22, (2 ), 681 –706 (2006). 0266-5611 CrossRef
Hamhaber  U. et al., “Comparison of quantitative shear wave MR-elastography with mechanical compression tests,” Magn. Reson. Med.. 49, (1 ), 71 –77 (2003). 0740-3194 CrossRef
Deffieux  T. et al., “On the effects of reflected waves in transient shear wave elastography,” IEEE Trans. Ultrason. Ferroelectr. Freq. Control. 58, (10 ), 2032 –2035 (2011). 0885-3010 CrossRef
Song  P. et al., “Fast shear compounding using directional filtering and two-dimensional shear wave speed calculation,” in  Proc. of the IEEE Joint UFFC, EFTF, PRF Symposium ,  Prague, Czech Republic , Vol. 1 (2013).
Deffieux  T. et al., “Shear wave spectroscopy for in vivo quantification of human soft tissues visco-elasticity,” IEEE Trans. Med. Imaging. 28, (3 ), 313 –322 (2009). 0278-0062 CrossRef
Nguyen  T.-M. et al., “Assessment of viscous and elastic properties of sub-wavelength layered soft tissues using shear wave spectroscopy: theoretical framework and in vitro experimental validation,” IEEE Trans. Ultrason. Ferroelectr. Freq. Control. 58, (11 ), 2305 –2315 (2011). 0885-3010 CrossRef

Some tools below are only available to our subscribers or users with an online account.

Related Content

Customize your page view by dragging & repositioning the boxes below.

Related Book Chapters

Topic Collections

Advertisement
  • Don't have an account?
  • Subscribe to the SPIE Digital Library
  • Create a FREE account to sign up for Digital Library content alerts and gain access to institutional subscriptions remotely.
Access This Article
Sign in or Create a personal account to Buy this article ($20 for members, $25 for non-members).
Access This Proceeding
Sign in or Create a personal account to Buy this article ($15 for members, $18 for non-members).
Access This Chapter

Access to SPIE eBooks is limited to subscribing institutions and is not available as part of a personal subscription. Print or electronic versions of individual SPIE books may be purchased via SPIE.org.