JBO Letters

Phase-resolved acoustic radiation force optical coherence elastography

[+] Author Affiliations
Wenjuan Qi

University of California, Irvine, Beckman Laser Institute, Irvine, 1002 Health Sciences Road East, California 92612

University of California, Irvine, Department of Chemical Engineering and Materials Science, California 92697

Ruimin Chen, Qifa Zhou

University of Southern California, Department of Biomedical Engineering, NIH Ultrasonic Transducer Resource Center, Los Angeles, California 90089

Lidek Chou, Gangjun Liu, Jun Zhang

University of California, Irvine, Beckman Laser Institute, Irvine, 1002 Health Sciences Road East, California 92612

Zhongping Chen

University of California, Irvine, Beckman Laser Institute, Irvine, 1002 Health Sciences Road East, California 92612

University of California, Irvine, Department of Chemical Engineering and Materials Science, California 92697

University of California, Irvine, Department of Biomedical Engineering, California 92697

J. Biomed. Opt. 17(11), 110505 (Nov 02, 2012). doi:10.1117/1.JBO.17.11.110505
History: Received June 21, 2012; Revised September 24, 2012; Accepted September 26, 2012
Text Size: A A A

Open Access Open Access

Abstract.  Many diseases involve changes in the biomechanical properties of tissue, and there is a close correlation between tissue elasticity and pathology. We report on the development of a phase-resolved acoustic radiation force optical coherence elastography method (ARF-OCE) to evaluate the elastic properties of tissue. This method utilizes chirped acoustic radiation force to produce excitation along the sample’s axial direction, and it uses phase-resolved optical coherence tomography (OCT) to measure the vibration of the sample. Under 500-Hz square wave modulated ARF signal excitation, phase change maps of tissue mimicking phantoms are generated by the ARF-OCE method, and the resulting Young’s modulus ratio is correlated with a standard compression test. The results verify that this technique could efficiently measure sample elastic properties accurately and quantitatively. Furthermore, a three-dimensional ARF-OCE image of the human atherosclerotic coronary artery is obtained. The result indicates that our dynamic phase-resolved ARF-OCE method can delineate tissues with different mechanical properties.

Figures in this Article

Many diseases involve change in the biomechanical properties of tissue. Cancerous tumors usually appear as hard nodules, and their stiffness differs from the surrounding tissue.1 A calcified plaque is markedly stiffer than the microenvironment around it. Elastography, quantitatively interrogating tissue mechanical properties by applying force to the tissue, has attracted increasing research interest. The beauty of elastography is that it is capable of delineating tissues with different mechanical properties even when their structures appear very similar. Magnetic resonance elastography and ultrasonic elastography are two mature imaging modalities being used for measuring mechanical properties at the tissue and organ level. Optical coherence elastography (OCE),26 like all forms of optical coherence tomography (OCT), has a superior micrometer scale resolution and is therefore suitable for imaging subtle mechanical changes in the early stages of disease.

Excitation and detection are generally two characteristics of an OCE system. OCE falls into two main categories with respect to the excitation method: static/quasi-static excitation OCE24 and dynamic excitation OCE.57 The former applies compression to the subject statically/quasi-statically and measures the relaxation dynamics, whereas the latter dynamically excites the subject using various waveforms and detects the induced displacement in a gated time window, which enables a much higher imaging speed. The recent reported dynamic OCE, featuring in direct-contact mechanical excitation method and a 5-kHz line rate, is applied for three-dimensional (3-D) skin imaging.7 In terms of the detection method, there are two widely used algorithms: speckle tracking and the phase-resolved method. In comparison with the speckle tracking algorithm, which suffers from low accuracy due to decorrelation, the phase-resolved method provides sensitive displacement measurement at nanometer resolution. This is especially important for intravascular imaging, because better sensitivity in displacement measurement minimizes the force that must be applied to a vessel wall to quantify tissue mechanical properties. Phase-resolved OCE had been investigated previously for shear wave/surface acoustic wave (SAW) propagation detection,8,9 but this method requires separation between excitation and detection points, which limits the speed and spatial resolution. Adaptation of this method is difficult for in vivo real-time two-dimensional (2-D) and 3-D vascular elastography.

In this paper, we propose a dynamic phase-resolved acoustic radiation force optical coherence elastography (ARF-OCE) method, which is capable of high-speed, high-spatial-resolution, point-by-point elastogram mapping in the axial direction. By applying a modulated acoustic radiation force to the sample and detecting the phase induced by the vibration of the sample using the phase-resolved OCT method, ARF-OCE is able to determine the relative Young’s modulus in various phantoms and human coronary artery tissues. The high-speed and high-spatial-resolution phase-resolved ARF-OCE has the potential to measure tissue elasticity (Young’s modulus) in vivo quantitatively.

Tissue Phantom and Tissue Preparation

The side-by-side agarose tissue phantom was fabricated. Two different concentrations of agarose (3.5% and 7% weight/volume) were used to obtain two different Young’s moduli in each side of the phantom. The 3-mm-thick sample was sectioned for OCE imaging.

A section of postmortem human atherosclerotic coronary artery was obtained and stored at 0°C. The tissue was cut open and flattened prior to the ARF-OCE imaging.

ARF-OCE System

An 890-nm spectral-domain OCT system with a 20-kHz acquisition rate, as shown in Fig. 1, was used for vibration detection. The axial and lateral resolutions of the system were measured to be 3.5 and 14.8 μm, respectively. The signal-noise-ratio (SNR) of the system was 100 dB with 650 μW of sample arm power and a 50-μs A-line rate. The minimum detectable phase for this system was measured to be 1.5 mrad, which corresponded to a velocity sensitivity of 2.13 μm/s. For the purposes of elastography, the sample was placed on the surface of a focused ultrasonic transducer, where the acoustic force was measured as roughly uniform by a highly accurate hydrophone. The transducer featured a resonant frequency of 4 MHz and was driven by a square wave modulated RF signal (50% duty cycle amplitude modulation) to generate periodical acoustic radiation pushes in the longitudinal direction on to the sample through a thin layer of ultrasound gel. The modulation frequency and amplitude were chosen for each sample in such a way that the phase difference induced between adjacent A-lines was large, to enhance the sensitivity, but less than |π|, to avoid phase wrapping. The phase difference was obtained from consequent A-lines while the system was scanning continuously. The oversampled data points were averaged by a factor of 4 and then used to calculate the phase shift for each A-line. From the phase shift information, the instantaneous axial velocity can be extracted as Display Formula

υ(x,z,t)=Δϕ(x,z,t)λ0/4πnτ,(1)
where the light source center wavelength λ0 is 890 nm, and the adjacent A-line time interval τ is 50 μs. x and z represent the lateral and axial location, Δϕ denotes the phase shift between two adjacent A-lines, and n is the tissue refraction index, which is assumed to be 1.4 in this paper. Within a certain customized time window, the displacement Δd and the axial strain ε of the sample are expressed asDisplay Formula
Δd=t1t2Δϕ(x,z,t)λ04πnτdt(2)
Display Formula
ε=t1t2Δϕ(x,z,t)λ04πnτz0dt,(3)
where z0 is the original thickness of the sample prior to ultrasonic stimulation. Young’s modulus is an important parameter that characterizes the elastic properties of the sample. It is defined by Display Formula
E=σε=F/Aε,(4)
where σ is the axial (or normal) stress (force per area) acting on the sample, F is the acoustic radiation force out of the ultrasonic transducer, and A is the sample surface area. Since the acoustic radiation force was uniform and the OCT probe beam size is much less than the ultrasound beam waist, the strain applied in the experiment is small enough (less than 0.1) that the material can be assumed as elastic, and Eq. (4) can be used as the first-order approximation to quantify the relative dynamic Young’s modulus in our dynamic OCE method. In this experiment, ARF-OCE of both the phantom and the human coronary artery samples were constructed while the samples were under vibration excitation.

Grahic Jump LocationF1 :

Schematic diagram of the spectral-domain ARF-OCE system.

The feasibility of the ARF-OCE system was first tested on a side-by-side agarose phantom. An OCE image was acquired on the sample with 2048 A-lines over a scanning length of 6 mm. A B-mode OCT image of the phantom, in which the two sides of different concentration agarose cannot be clearly distinguished, is shown in Fig. 2(a). The phase shift between two successive A-lines was calculated using the phase-resolved method, and the resultant phase maps of the OCE image of the phantom are presented in Fig. 2(b) and 2(c). The displacement of the sample over time is proportional to the adjacent A-line phase shift at each lateral location and is calculated according to Eq. (2). As shown in Fig. 2(b) and 2(c), the boundary (red arrow) between the different concentrations of agarose in the phantom is clear in the OCE image, which demonstrates the advantage of OCE over OCT intensity imaging. The displacements of the two sides were 2.291 and 0.748 μm, respectively, which corresponded to 0.076% and 0.025% ARF induced strain, given the thickness of each side of the phantom was 3 mm. Because the radiation force was distributed equally in the region of interest, a 1:3.05 Young’s modulus ratio between the two sides of the phantom was obtained according to Eq. (4). The Young’s modulus of the agarose phantom was also measured under a compression test (MTS Synergie 100). The Young’s moduli for 3.5% and 7% agarose samples were 83.6 and 265.7 kPa, respectively. The results of the two methods agree with each other very well; the Young’s modulus of the 7% agarose sample is also triple the one for the 3.5% sample, which indicates that our ARF-OCE method is capable of quantitatively characterizing the mechanical properties of sample materials. Figure 2(d)2(f) shows the three-dimensional ARF-OCE performed on the agarose tissue phantom with 7% agarose on the right and 3.5% on the left. In the 3-D OCT image, shown in Fig. 2(d), it is difficult to differentiate the two sides of the phantom with different stiffness. However, the 3-D ARF-OCE image, shown in Fig. 2(e), clearly delineates the two materials, which differ in their stiffness. The color difference indicates different vibration levels. The pink color indicates a stronger vibration level, and the green-black color represents stiffer material, which vibrates more weakly under the same excitation. This result confirms our method’s capacity to provide an additional and effective contrast in material stiffness.

Grahic Jump LocationF2 :

(a) OCT intensity image of side-by-side agarose phantom. (b) The OCE image clearly reveals the boundary of the phantom when 500-Hz, 200-mV square wave modulated excitation was applied to the phantom. (c) Phase amplitude image. Red arrows indicate the boundary between the two sides, which differ in agarose concentration. (d) 3-D OCT, (e) AF-OCE, and (f) fused OCT and AF-OCE images of an agarose tissue-mimicking phantom. The phantom is excited with a 500-Hz square wave modulated RF signal. Scale bars: 500 μm.

To test the consistency of this method, we determined the displacement as a function of applied voltage, which is proportional to the acoustic radiation force out of the ultrasonic transducer (shown in Fig. 3). Both the hard and the soft samples acted linearly when the applied voltage increased. The Young’s modulus ratio was calculated as the average value of the ratio between two samples under each voltage. A 1:1.74 ratio was given by the phase-resolved ARF-OCE measurement, which agrees with the 1:1.83 ratio from the standard compression test.

Grahic Jump LocationF3 :

Axial displacement as a function of preamplified voltage to the transducer.

To demonstrate the potential of this technology in intravascular imaging, we have performed imaging of a section of atherosclerotic human cadaver coronary artery with the phase-resolved ARF-OCE (Fig. 4). An OCT intensity image, shown in Fig. 4(a), provided a general morphological view of the tissue, but no obvious evidence of atherosclerosis was found. However, in Fig. 4(b), a strong vibration phase contrast could be clearly imaged using ARF-OCE. The red-colored region indicated by the blue arrow in the ARF-OCE image is characterized by smaller phase and less vibration and, therefore, represents less elastic, stiffer tissue such as plaques. The stronger vibration representing softer tissue (dark area) is indicated by yellow arrows. In this ARF-OCE image, the stiffer tissue, corresponding to the atherosclerotic lesion in the blue box in the histological image [Fig. 4(c)], is clearly distinguishable from the softer (normal) tissue area. The transition between plaque and normal regions appears orange in color (indicated by the red arrow) and is characterized by intermediate elasticity. This preliminary result clearly shows that our ARF-OCE method has great potential to characterize tissue mechanical properties quantitatively and thereby delineate diseased tissue from normal tissue.

Grahic Jump LocationF4 :

(a) OCT structural and (b) ARF-OCE phase images of a human cadaver coronary artery under 500-Hz, 350-mV AM modulated excitation; an atherosclerotic lesion was identified as the red region corresponding to the blue box in (c) the histological image [corresponding to across blue line location in (a)] and (d) the close-up view of the atherosclerotic lesion. Scale bars: 1 mm.

We have developed a dynamic phase-resolved ARF-OCE method for accessing tissue mechanical properties along the axial direction. The feasibility of the method was first tested on different tissue-mimicking phantoms. The ratio of Young’s moduli of 3.5% and 7% agarose concentration phantoms obtained from the aforementioned method corresponded very closely with the standard compression test result. Measurements using this method were also confirmed to be consistent by analyzing displacement under varying radiation forces. We have performed in vitro imaging of an atherosclerotic human cadaver coronary artery specimen with ARF-OCE. The atherosclerosis was found in the ARF-OCE image, while no obvious evidence appeared in the OCT image. This result shows that our ARF-OCE method has great potential to delineate diseased tissue from normal tissue and quantitatively characterize mechanical properties of tissues.

The authors thank Claire Robertson for the phantom compression test, Jiawen Li for preparing the tissue, and the individuals who donated their bodies and tissues for the advancement of education and research. This work was supported by the National Institutes of Health (EB-10090, EY-021519, HL-105215, HL-103764, and EB-015890), the Air Force Office of Scientific Research (FA9550-04-0101), and the Arnold and Mabel Beckman Foundation. Dr. Chen has a financial interest in OCT Medical Imaging Inc., which, however, did not support this work.

Garra  B. S. et al., “Elastography of breast lesions: initial clinical results,” Radiology. 202, (1 ), 79 –86 (1997). 0033-8419 
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
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
Adie  S. G. et al., “Spectroscopic optical coherence elastography,” Opt. Express. 18, (25 ), 25519 –25534 (2010). 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
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
Kennedy  B. F. et al., “In vivo three-dimensional optical coherence elastography,” Opt. Express. 19, (7 ), 6623 –6634 (2011). 1094-4087 CrossRef
Razani  M. et al., “Feasibility of optical coherence elastography measurements of shear wave propagation in homogeneous tissue equivalent phantoms,” Biomed. Opt. Express. 3, (5 ), 972 –980 (2012). 2156-7085 CrossRef
Li  C. et al., “Quantitative elastography provided by surface acoustic waves measured by phase-sensitive optical coherence tomography,” Opt. Lett.. 37, (4 ), 722 –724 (2012). 0146-9592 CrossRef
© 2012 Society of Photo-Optical Instrumentation Engineers

Citation

Wenjuan Qi ; Ruimin Chen ; Lidek Chou ; Gangjun Liu ; Jun Zhang, et al.
"Phase-resolved acoustic radiation force optical coherence elastography", J. Biomed. Opt. 17(11), 110505 (Nov 02, 2012). ; http://dx.doi.org/10.1117/1.JBO.17.11.110505


Figures

Grahic Jump LocationF1 :

Schematic diagram of the spectral-domain ARF-OCE system.

Grahic Jump LocationF2 :

(a) OCT intensity image of side-by-side agarose phantom. (b) The OCE image clearly reveals the boundary of the phantom when 500-Hz, 200-mV square wave modulated excitation was applied to the phantom. (c) Phase amplitude image. Red arrows indicate the boundary between the two sides, which differ in agarose concentration. (d) 3-D OCT, (e) AF-OCE, and (f) fused OCT and AF-OCE images of an agarose tissue-mimicking phantom. The phantom is excited with a 500-Hz square wave modulated RF signal. Scale bars: 500 μm.

Grahic Jump LocationF3 :

Axial displacement as a function of preamplified voltage to the transducer.

Grahic Jump LocationF4 :

(a) OCT structural and (b) ARF-OCE phase images of a human cadaver coronary artery under 500-Hz, 350-mV AM modulated excitation; an atherosclerotic lesion was identified as the red region corresponding to the blue box in (c) the histological image [corresponding to across blue line location in (a)] and (d) the close-up view of the atherosclerotic lesion. Scale bars: 1 mm.

Tables

References

Garra  B. S. et al., “Elastography of breast lesions: initial clinical results,” Radiology. 202, (1 ), 79 –86 (1997). 0033-8419 
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
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
Adie  S. G. et al., “Spectroscopic optical coherence elastography,” Opt. Express. 18, (25 ), 25519 –25534 (2010). 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
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
Kennedy  B. F. et al., “In vivo three-dimensional optical coherence elastography,” Opt. Express. 19, (7 ), 6623 –6634 (2011). 1094-4087 CrossRef
Razani  M. et al., “Feasibility of optical coherence elastography measurements of shear wave propagation in homogeneous tissue equivalent phantoms,” Biomed. Opt. Express. 3, (5 ), 972 –980 (2012). 2156-7085 CrossRef
Li  C. et al., “Quantitative elastography provided by surface acoustic waves measured by phase-sensitive optical coherence tomography,” Opt. Lett.. 37, (4 ), 722 –724 (2012). 0146-9592 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.