*in vivo*resolution, which allows for the identification of morphological tissue structures. Optical coherence elastography (OCE) benefits from the superior spatial resolution of OCT and has promising applications, including cancer diagnosis and the detailed characterization of arterial wall biomechanics, both of which are based on the elastic properties of the tissue under investigation. We present OCE principles based on techniques associated with static and dynamic tissue excitation, and their corresponding elastogram image-reconstruction algorithms are reviewed. OCE techniques, including the development of intravascular- or catheter-based OCE, are in their early stages of development but show great promise for surgical oncology or intravascular cardiology applications.

## 1.

## Introduction

## 1.1.

### Elastography

Elastography is an experimental method and a noninvasive medical imaging technique that conveys new information about internal tissue structure and behavior under load and ultimately generates images called elastograms.^{1} Examples include elastograms of axial strains, lateral strains, and elastic moduli. Different degrees of stiffness are demonstrated as different shades of light and dark in an elastogram.^{2} Ideally, this elastography technique can extract the elastic modulus, which is a property of the constituent material describing the stiffness of an isotropic elastic material.^{3} There are several types of elastic moduli, where the main moduli used currently in elastography are the Young's modulus (the ratio of the uniaxial stress over the uniaxial strain) and shear modulus (the ratio of shear stress to the shear strain). Given any two of the six elastic moduli, the remaining modulus can be determined. Thus two elastic moduli can fully describe the elastic properties of a homogeneous and isotropic material.

The motivation underlying elastography stems from the fact that elastic properties of tissues are directly related to the underlying structure of the tissue and are strongly affected by pathological changes. For example, edema, ﬁbrosis, and calciﬁcation alter the elastic modulus of the extracellular tissue matrix;^{4} different tissues within atherosclerotic plaque have distinctive elastic properties;^{5} and cancerous tumors are often stiffer than benign and normal tissue.^{6} Elastography images contain information of local variations of the stiffness inside a region of interest and may provide additional clinical information to aid the identification of suspicious lesions and diagnosis of various disease states.^{7, 8, 9} Currently, there are several potential applications of elastography, which include assessing breast or brain tissue for malignancy^{10, 11} and atherosclerotic arteries for vulnerability to myocardial infarction.^{12}

## 1.2.

### Optical Coherence Elastography

Optical coherence tomography (OCT)–based elastography, which is termed optical coherence elastography (OCE),^{13} was ﬁrst introduced by Schmitt^{14} in 1998. OCE shows great potential for micron and submicron imaging applications because it benefits from the high resolution of OCT,^{15} while additionally providing the elastic properties of the sample. OCT provides structural images approaching histology, where the microstructures of biological tissues can be quantified based on the optical backscattering properties within the imaging region of interest. This beneficial high-resolution, noninvasive imaging modality allows OCE to evaluate the mechanics of intact tissue on a scale that cannot be offered by elastographies produced by the competing imaging modalities of ultrasound^{3, 16, 17, 18} or MRI (magnetic resonance imaging).^{19, 20} Recent advances of these technologies are pushing the resolution to higher values where high-frequency ultrasound^{21} can reach 30–70 μm and magnetic resonance microscopy^{22} has improved imaging resolution to ∼120 um. However, OCT has inherently high resolution, 1–10 μm *in vivo*; therefore, OCE has the potential to link microstructural data with biomechanical and functional information. In the last decade, several studies have demonstrated the ability of OCE to provide valuable information about the micromechanical properties of tissue. OCE has been used to investigate the microscopic deformation of gelatin tissue phantoms, porcine muscle, and human skin as a function of depth under compressive loads;^{23} define the properties of collagen organization in musculoskeletal tissue;^{24} provide the material properties of the cornea;^{25} develop strain maps of engineered tissue and a *Xenopus* tadpole at different stages of embryonic development;^{13} and define tumor and normal adipose human breast tissue borders.^{26} Finally, one of the current research directions in OCE imaging has been to develop intravascular OCE for rupture-prone (vulnerable) plaque detection.^{27}

## 1.3.

### Intravascular Elastography for Atherosclerotic Plaque Detection

The main components of atherosclerotic plaque include a large extracellular necrotic core with a thin ﬁbrous cap inﬁltrated by macrophages, all of which have varying biomechanical properties.^{28, 29} Rupture of the cap induces the formation of a thrombus, which may obstruct the coronary artery, causing an acute heart attack and frequently results in patient death.^{30, 31} Detailed characterization of the arterial wall biomechanics may provide complementary information to quantify lesion stability. Both *in vitro* and *in vivo* studies have revealed that strain is higher in fatty tissue components than in ﬁbrous plaques.^{32} The presence of a high-strain area that is surrounded by a low-strain region has a high predictive power to identify the vulnerable plaque with high sensitivity and speciﬁcity.^{33} Thus, changes in the mechanical properties of the tissue can be detected, quantiﬁed, and correlated with clinical symptoms and inﬂammation markers.

Intravascular elastography is particularly used to correlate elastograms with histological characteristics of the blood vessel wall. The ultrasound community has made signiﬁcant progress in developing intravascular ultrasound (IVUS) elastography,^{34} which has been successfully used to extract arterial radial strains with a spatial resolution of 200 μm. Schmitt
^{35} have used ultrasound elastography within carotid walls induced by the natural cardiac pulsation, and coupled to a simple inverse problem, to recover the wall elastic modulus at the blood pulsatility frequency (1 Hz). This was the ﬁrst *in vivo* axial strain and shear distributions of atherosclerotic carotid arteries measurement. In 2008, Schmitt
^{36} made further progress by adapting the dynamic microelastography method and formulating an inverse problem to study the radial viscoelasticity of the vessel wall. Although this has provided some success in plaque identiﬁcation, there exist several limitations because vulnerable atherosclerotic plaques have structural components (e.g., ﬁbrous caps) on the order of 50–200 μm, which lie below the resolvable limit of the IVUS imaging system. OCT is a potential technical solution to overcome these hurdles because it is not limited, in this regard, due to its inherent high resolution and provides noninvasive near-cellular-level imaging for plaque quantification.

An important application of intravascular OCT^{37, 38} is the detection of thin-cap ﬁbroatheroma (TCFA), the most important type of vulnerable plaques. Xu
^{39} have addressed fundamental issues that underlie the tissue characterization of OCT images obtained from coronary arteries. With high spatial resolution and, thus, the ability to visualize structures on the size scale of thin ﬁbrous caps, intravascular OCE has the potential to provide high-resolution characterization of strains in tissue lying within 1.0–1.5 mm of the lumen interface, the region most susceptible to plaque disruption.^{27} Chau
^{32} have used OCT as a basis for ﬁnite element analysis and obtained results comparable to histology methods for stress and strain analysis of atherosclerosis. Their results suggest that OCT-based ﬁnite element analysis may be a powerful tool for investigating coronary atherosclerosis and detecting vulnerable plaque. The capability of OCE to resolve TCFA is advantageous when compared to competing technologies, such as IVUS. Chan
^{27} and Khalil
^{40, 41} have performed several studies highlighting the benefits of intravascular OCE, which may play an important role in the elasticity estimation of vulnerable plaques and provide patient risk stratification for appropriate treatment or vascular targeted therapies.

## 1.4.

### Technique Development of Optical Coherence Elastography

When considering OCE as a diagnostic tool, two important issues must be considered. The first is the mechanical design of a loading system because an external load is usually required to excite tissues for a certain deformation. Therefore, a loading system must be able to apply a known and reproducible stress to avoid artifacts from nonuniform stress distributions, be easy to use, and be compatible with the existing time-domain (TD) or frequency-domain (FD) OCT systems. The second issue relates to data analysis for reconstruction of tissue elastograms, which is a process of combining imaging, elastography, and computational modeling to build maps of soft-tissue mechanical properties. In general, soft tissues are viscoelastic, anisotropic, and incompressible. No biological soft tissue exhibits a truly elastic response, but under some conditions, the assumption of linear elasticity is both reasonable and useful.^{42} For example, when the deformation is small (<0.1%), the nonlinearity of the elastic constants of soft tissues can be neglected.^{25} Under the assumption that the tissue is linear elastic, isotropic, and subjected to a constant stress ﬁeld, the strain ﬁeld could be interpreted as a relative measure of elasticity distribution, where the strain is large in a compliant (soft) tissue, and small in a rigid (hard) tissue.^{43} However, for anisotropic materials, such as arteries with asymmetric geometries and the presence of plaques leading to a nonhomogeneous arterial wall stress ﬁeld, the material property identiﬁcation is much more complicated than for isotropic materials.

Inverse problem-solving^{43, 44} methods are usually developed to solve for unknown material properties from measured displacements in anisotropic materials. Prior to solving the inverse problem, constitutive relations need to be defined to describe the behavior of individual tissues depending on the conditions of interest, including the parameters of the loading proﬁle, material properties, geometric feature of the structure, boundary conditions, and/or a combination of these behaviors. Then, if a set of experimentally measured data (displacement, velocity, strain, natural frequency, etc., obtained by OCE) are available for the structure of interest, the material property can be characterized by minimizing the sum of the squares of the errors or discrepancy between the experimental and computed structure behavior data.^{45, 46} Elastograms of strain or modulus can then be reconstructed. Although hurdles to widespread clinical acceptance exist, short-term engineering solutions are being investigated, such that OCE may play a very important role as a diagnostic tool in the near future.

This review summarizes recent advances in OCE, including technical issues involved in the design of loading schemes associated with OCT platforms along with theoretical analysis for displacement, strain, and elastic moduli elastogram reconstruction. Clinical applications including intravascular OCE are presented, followed by a discussion of the future role of OCE in terms of technical development and clinical utility.

## 2.

## Optical Coherence Elastography Approaches

Karimi
^{47} demonstrated an overall ﬂowchart of the elastography process shown as Fig. 1, where the main steps include the following:

Identify the structural characteristics of tissue (i.e., isotropic or anisotropic) and establish an appropriate theoretical framework with suitable restrictions on the possible constitutive relations connecting applied stresses or forces to strains or deformations.

Design appropriate experimental methods, including the detection of the natural deformation of tissue due to factors such as blood ﬂow or external excitation based on the theoretical framework. Displacement or strain elastograms can then be obtained through medical imaging modalities, including MRI, ultrasound, and OCT.

Obtain a reliable stress-strain constitutive relation and elastic modulus elastograms. This process is based on an iterative scheme to minimize the difference between the strain obtained by the constitutive model and experimental results.

The displacement or strain elastograms are usually obtained through directly processing structural images obtained by biomedical imaging modalities, such as ultrasound or OCT, while elastic moduli elastograms are reconstructed by further calculations. Similar to untrasound elastography, external forces are usually applied to deform soft tissues although ultrasound, specifically high-frequency ultrasound, can use its interrogating acoustic radiation force to stimulate soft tissues. There are two particular types of stimuli, static/quasi-static and dynamic stimuli. An OCE system includes two main subsystems: an OCT system and a loading system used to deform biological tissues. The following sections review OCE systems and the corresponding algorithms for elastogram reconstruction.

## 3.1.

### Static or Quasi-Static OCE Method

OCE takes advantage of the short coherence lengths of broadband light sources to image micron-level displacement and strain distributions of individual tissue elements as they are stretched, compressed, or vibrated by an external force.^{48} When static or quasi-static compression is applied to tissue, the majority of OCE techniques reported to date rely on speckle tracking^{49} to estimate the relative motion of subsurface structures imaged under different applied forces. The following sections present and discuss these OCE techniques in detail.

## 3.1.1.

#### Traditional Cross-Correlation-Based Method

The cross-correlation-based OCE methodology can work for both the TD and FD OCT platforms. Typically, a sequence of B-scan OCT images are acquired before and after deformation of an object in conjunction with a loading system. The loading system uses a piezoelectric transducer attached to a glass plate, through which the OCT imaging beam is focused. A computer-controlled translation stage moves a load cell and applies necessary compression to the sample against the ﬁxed glass slide,^{13, 50} so that the force applied to the tissue in the axial direction (*z* direction) can be known from the load cell and the actual displacement of the sample can be acquired from the translation stage. Displacements are then calculated by correlating image sets obtained prior to and immediately following the application of the load using Eq. 1 or equivalent forms of 2-D cross-correlation equations,

## 1

[TeX:] \documentclass[12pt]{minimal}\begin{document}\begin{eqnarray} \rho (x',z') = \frac{{\sum_{x,z} {[I_1 (x,z) - \overline {I_1 }][I_2 (x - x',z - z') - \overline {I_2 }]} }}{{\sqrt {\sum_{x,z} {[I_1 (x,z) - \overline {I_1 }]^2 [I_2 (x - x',z - z) - \overline {I_2 }]^2 } } }},\nonumber\\ \end{eqnarray}\end{document} $$\begin{array}{c}\hfill \rho ({x}^{\prime},{z}^{\prime})=\frac{{\sum}_{x,z}[{I}_{1}(x,z)-\overline{{I}_{1}}][{I}_{2}(x-{x}^{\prime},z-{z}^{\prime})-\overline{{I}_{2}}]}{\sqrt{{\sum}_{x,z}{[{I}_{1}(x,z)-\overline{{I}_{1}}]}^{2}{[{I}_{2}(x-{x}^{\prime},z-z)-\overline{{I}_{2}}]}^{2}}},\end{array}$$*x*

^{′},

*z*

^{′}) is the correlation coefficient, which is a function of displacements

*x*

^{′}and

*z*

^{′}within a predefined subset;

*I*

_{1}(

*x*,

*z*) is the pixel intensity of a point (

*x*,

*z*) and

*I*

_{2}(

*x*–

*x*

^{′},

*z*–

*z*) is the pixel intensity of the point on the deformed image;

*x*is the transverse location, and

*z*is the axial (depth) location within the images. [TeX:] $\overline {I_1 } $ $\overline{{I}_{1}}$ and [TeX:] $\overline {I_2 } $ $\overline{{I}_{2}}$ are the average subset pixel values. From the resultant deformation pattern, strain images are obtained by numerically differentiating the axial displacement component assuming the applied stress was uniform or had a known distribution. In this case, the 2-D cross-correlation is applied to the

*xz*plane(B-scan of OCT) instead of the

*yz*or

*xy*plane because B-scans are recorded in real time, whereas images in other directions can only be obtained when the 3-D OCT scan is finished. Because of the well-developed cross-correlation algorithms, researchers have shown interest in cross-correlation-based OCE methods (Schmitt,

^{14}Duncan and Kirkpatrick,

^{49}Ko,

^{51}Rogowska

^{52}).

Cross-correlation-based OCE has shown promise as a technique to provide additional contrast for disease detection. However, several hurdles exist that may limit its usefulness and adoption into daily imaging procedures for diagnosis, disease screening, or progression monitoring. First, correlation-based algorithms are based on the maximization of a correlation coefficient that is determined by examining pixel-intensity array subsets on two or more corresponding images. During the correlation calculation, the size of the correlation window must be chosen carefully. If the window is too small, then the displacement results may not be accurate. However, large windows tend to average out the difference in displacements of small features in images and decrease the ability of speckle tracking to make microstructural assessments. This is a time-consuming high-complexity computational procedure, making real-time measurements and characterization difficult. Second, it is required that the resultant deformation of tissue must be large enough, where the peak of correlation coefficient shift can be quantified, but small enough such that the correlation of the speckle patterns does not lie within the potential sources of measurement error: (*i*) excessive surface displacement (either in-plane, out-of-plane, or rotational) due to vibration, (*ii*) excessive surface tilt, or (*iii*) plastic deformation, which may induce speckle-pattern decorrelation.^{53} Thus far, strain is calculated from the gradient of assessed displacement ﬁelds and the differentiation operation has a high sensitivity to local errors. Thus, cross-correlation approaches are not applicable when the deformation is either very large or very small.

Kirkpatrick
^{50} used a maximum-likelihood estimation method to deal with small deformations, by which the local velocity of every pixel point in the image was assessed and the strain elastograms were then obtained. Instead of using cross-correlations, they also proposed a phase-sensitive OCE to process large deformations.

## 3.1.2.

#### Phase-Sensitive Optical Coherence Elastography

Phase-sensitive OCE (PSOCE) is based on the measurement of phase changes induced by tissue motions in complex OCT images. Kirkpatrick
^{50} and Wang
^{54} developed a tissue Doppler OCE (tDOCE) approach to obtain large movement of speckles from Doppler frequency shift induced by instantaneous tissue motion. Theoretical basis of tDOCE is similar with Doppler OCT.^{55, 56} In the spectral domain, the spectral density of combined reference and sample beams is^{50}

## 2

[TeX:] \documentclass[12pt]{minimal}\begin{document}\begin{eqnarray} I_d (\upsilon) &=& S(\upsilon)\left[ 1 + \int R(\tau)d\tau + 2R_{{\rm self}}\right.\nonumber\\ &&\left. +\, 2\int {\sqrt {R(\tau)} \cos (2\pi \upsilon \tau + \phi _0)d\tau } \right], \end{eqnarray}\end{document} $$\begin{array}{ccc}\hfill {I}_{d}\left(\upsilon \right)& =& S\left(\upsilon \right)\left[1+\int R\left(\tau \right)d\tau +2{R}_{\mathrm{self}}\right.\hfill \\ & & \left.+\phantom{\rule{0.16em}{0ex}}2\int \sqrt{R\left(\tau \right)}\mathrm{cos}(2\pi \upsilon \tau +{\phi}_{0})d\tau \right],\hfill \end{array}$$*R*(τ) is the normalized intensity backscattered from the scatter at time τ, ϕ

_{0}is the initial phase value at τ = 0, and

*S*(υ) is the spectral density of the light source. The fourth term in the brackets represents the spectral interference between the sample and the reference and is the term that contributes to the OCT signal. Discard the first three terms and perform Fourier transform; the depth information regarding the local scatters in the sample is obtained:

## 3

[TeX:] \documentclass[12pt]{minimal}\begin{document}\begin{equation} I(z) = A(z)\exp [ - i\Phi (z)]. \end{equation}\end{document} $$I\left(z\right)=A\left(z\right)\mathrm{exp}[-i\Phi (z\left)\right].$$*z*is translated by a distance Δ

*d*(

*z*) during the time interval Δ

*t*between two successive A-scans, then the induced phase change will be

## 4

[TeX:] \documentclass[12pt]{minimal}\begin{document}\begin{equation} \Delta \Phi (z,t) = 2nk\Delta d(z), \end{equation}\end{document} $$\Delta \Phi (z,t)=2nk\Delta d\left(z\right),$$*k*is the wavenumber of the light source and

*n*is the refractive index of the sample. The depth-resolved tissue velocity

*v*(

*z,t*) in the beam direction could then be obtained from the phase change, ΔΦ(

*z,t*) by

## 5

[TeX:] \documentclass[12pt]{minimal}\begin{document}\begin{equation} v(z,t) = \frac{{\Delta d}}{{\Delta t}} = \frac{{\Delta \Phi (z,t)\lambda }}{{4\pi n\Delta t}},\end{equation}\end{document} $$v(z,t)=\frac{\Delta d}{\Delta t}=\frac{\Delta \Phi (z,t)\lambda}{4\pi n\Delta t},$$## 6

[TeX:] \documentclass[12pt]{minimal}\begin{document}\begin{equation} \dot \varepsilon (z,t) = \frac{{v(z,t)}}{{z_0 }} = \frac{{\Delta \Phi (z,t)\lambda }}{{4\pi nz_0 \Delta t}},\end{equation}\end{document} $$\stackrel{\u0307}{\varepsilon}(z,t)=\frac{v(z,t)}{{z}_{0}}=\frac{\Delta \Phi (z,t)\lambda}{4\pi n{z}_{0}\Delta t},$$*z*

_{0}is the initial depth of the sample before the tissue movement. The instantaneous tissue deformation and strain rate could then be imaged in real time, and the total displacement and strain over a time period can be obtained by integration of the instantaneous displacement and strain over the elapsed time. This approach was expanded to detecting phase change between successive OCT-B scans so that small deformation can also be measured.

^{57}tDOCE is faster than cross-correlation-based OCE methods, but the sensitivity is fundamentally limited by the phase sensitivity of the OCT system. Example images of velocity, displacement, strain rate, and cumulative strain maps of a phantom are illustrated in Fig. 2,

^{57}where the phantom was compressed repeatedly by a triangular wave.

The techniques mentioned above can obtain deformation or strain elastograms, which may be enough to show the stiffness difference between healthy and diseased tissues in scenarios such as cancer diagnosis. However, for the case of atherosclerotic plaque characterization, the absolute value of modulus of elasticity needs to be known. Moreover, as intrinsic mechanical properties of materials, the elastic moduli are thought to be more advantageous for wider clinical applications.^{1} Inverse problem approaches can be applied to aid the reconstruction of distributed elastic moduli elastograms.

## 3.1.3.

#### Tissue Elasticity Reconstruction and Inverse Problem

Traditionally, elastography imaging protocols typically assume a “constant stress ﬁeld,” implying that an elastogram is equivalent to a modulus image. If the applied force is known and the strain is obtained from OCE, then the Young's modulus can be estimated, assuming a tissue is uniform, incompressible, and an isotropic linear elastic material. With the application of uniformly applied stress in the axial direction, the local stress is equal to the applied stress (constant throughout the sample), which can be calculated to yield the Young's modulus:^{52}

## 7

[TeX:] \documentclass[12pt]{minimal}\begin{document}\begin{equation} E = \frac{\sigma }{\varepsilon }, \end{equation}\end{document} $$E=\frac{\sigma}{\varepsilon},$$*E*is the Young's modulus, σ is stress, and ɛ is strain. Researchers use this simple, direct linear elastic, isotropic model for two main reasons. First, elasticity estimation is a discrete numerical effort, where elasticity values are assigned to individual elements based on displacement maps generated via image processing. Second, unlike a nonlinear constitutive model, elasticity estimations provide a comprehensive map of material properties.

^{40}However, as previously mentioned, biological tissues are anisotropic and their response to specific external excitation has a viscous as well as an elastic component. Therefore, the determination of both elastic and viscoelastic properties of soft tissues is important for full characterization of tissue. In this case, a more elaborate solution strategy is required to solve the inverse problem of elastic modulus imaging.

Karimi
^{47} have performed sophisticated attempts on using inverse problem approaches to generate elastic moduli elastograms considering tissue's nonlinear behavior. Following the flowchart shown in Fig. 1, they implemented the Mooney-Rivlin model,^{58} a hyperelastic material model, to estimate the mechanical behavior of the corresponding regions in their finite element method (FEM) model of an aortic specimen, including normal vessel wall, ﬁbrous plaque, and the presence of lipid pools. An iteration scheme, multipoint scheme, and a genetic algorithm are used for material parameters estimation by minimizing the difference between the strain map obtained with elastography and the strain map calculated from the FEM model. The inverse problem technique^{59} has experienced initial preclinical success in several areas, where the biomechanical properties of tissue may be used to image and quantify disease formation. For example, Sridhar
^{60} applied inverse problem solutions to investigate how disease-related changes to molecular bonding within stromal tissues affect the broad spectrum of viscoelastic responses^{61} and proposed a new approach to minimize the mean-squared error between computed and measured displacements.

The greatest potential for the inverse problem technique aided OCE lies in patient risk stratification for management of vulnerable plaques. Floch
^{62} summarized different methods to constrain the optimization algorithm of solving inverse problems and studied the preconditioning of the algorithm based on the best estimation of the plaque component contours. However, technical investigations into the meaning of associated parameters, such as the role of the magnitude of the regularization parameter and extensive *in vivo* testing, are required prior to the acceptance of involving inverse problem approaches in standard clinical imaging scenarios.

## 3.2.

### Dynamic Optical Coherence Elastography Method

In dynamic excitation, samples are excited by various mechanical wave forms, where the distributions of dynamic displacements (time or space) within a small region of the samples are imaged to capture dynamic OCE signals. These OCE systems have high resolution in both axial and transverse directions with the major advantage of not requiring a cross-correlation algorithm.

Liang
^{26, 63, 65} and Xing and Boppart^{64} proposed a dynamic optical coherence elastography method in which tissue was excited by optomechanical waves, such as a step or sinusoidal wave form or acoustic pressure, coaxially with the optical beam, with a piezoelectric rod actuator operating in a frequency range from a 100 Hz to kilohertz. The Young's moduli were then extracted from the detected phase-resolved signals. In detail, the sample and mechanical wave driver were modeled as a Voigt body, via the following equation:^{63}

## 8

[TeX:] \documentclass[12pt]{minimal}\begin{document}\begin{equation} m\ddot x + (\gamma _0 + \gamma)\dot x + (k_0 + k)x = F - mg, \end{equation}\end{document} $$m\stackrel{\u0308}{x}+({\gamma}_{0}+\gamma )\stackrel{\u0307}{x}+({k}_{0}+k)x=F-mg,$$*x*is the displacement of samples;

*k*

_{0}and γ

_{0}are the spring constant and the damping parameter for the mechanical wave driver, respectively;

*k*and γ are the spring constant and the coefficient of viscosity of the samples, and

*m*denotes the mass of sample stage and the samples. As described in the previous section, the phase variation of OCT signals is related to displacement of samples being measured based on Eq. 4, so the phases obtained from OCT can be fitted to the solution of Eq. 8. During the fitting process, parameters in Eq. 8 are solved. Then the elastic modulus of the sample can be calculated. Using this technique, an elasticity map of tissue containing both tumor and normal adipose breast tissue were obtained and displayed in Fig. 3.

^{26}The measured elasticity map is shown in Fig. 3c, where elasticity differences on the tumor side of the OCE map are observed.

Acoustic radiation force are currently popularly used to generate an internal load noninvasively.^{66} Liang
^{67} presented an internal mechanical excitation method using acoustic radiation forces and developed an acoustomotive OCE based on the similar mechanical model introduced above. The acoustomotive OCE has the potential for imaging biomechanical properties *in vivo*. However, using these dynamic methodologies, every pixel in the OCE image or biomechanical map requires wave-equation ﬁtting and, as a result, the processing time is too long for real-time applications. In addition, the feasibility of accurately detecting vibration amplitude in the presence of lateral scatterer motion, such as when shear waves are excited, needs to be investigated further for applicability in this model.

Adie
^{68} proposed an OCE technique based on analysis of the dynamic interferometry caused by sample excitation at audio frequencies. This approach is suitable for quantitative measurement of tissue elastic properties in the 100 Hz to kilohertz frequency range, which matches the typical frequency response range of many biological tissues. The theoretical basis of this method is to expand the detected photocurrent signal of interest as a series of Bessel functions resulting in

## 9

[TeX:] \documentclass[12pt]{minimal}\begin{document}\begin{eqnarray} i_1 (z,t) &=& 2\rho \sqrt {I_{\rm R} I_{\rm S} } \left(\cos \phi _{{\rm DC}} \left[ J_0 \left({\phi _S } \right)\vphantom{\sum\limits_{n = 1}^\infty}\right.\right.\nonumber\\ &&\left.\left. +\, 2\sum\limits_{n = 1}^\infty {J_{2n} \left({\phi _S } \right)\cos \left({4\pi n\Omega t} \right)} \right] \right.\nonumber\\ &&\left. { -\, \sin \phi _{{\rm DC}} \left\{ {2\sum\limits_{n = 0}^\infty {J_{2n + 1} \left({\phi _{\rm S} } \right)\sin \left[ {2\pi \left({2n + 1} \right)\Omega t} \right]} } \right\}} \right),\hskip-6pt\nonumber\\ \end{eqnarray}\end{document} $$\begin{array}{ccc}\hfill {i}_{1}(z,t)& =& 2\rho \sqrt{{I}_{\mathrm{R}}{I}_{\mathrm{S}}}\left(\mathrm{cos}{\phi}_{\mathrm{DC}}\left[{J}_{0}\left({\phi}_{S}\right)\phantom{\sum _{n=1}^{\infty}}\right.\right.\hfill \\ & & \left.+\phantom{\rule{0.16em}{0ex}}2\sum _{n=1}^{\infty}{J}_{2n}\left({\phi}_{S}\right)\mathrm{cos}\left(4\pi n\Omega t\right)\right]\hfill \\ & & \left.-\phantom{\rule{0.16em}{0ex}}\mathrm{sin}{\phi}_{\mathrm{DC}}\left\{2\sum _{n=0}^{\infty}{J}_{2n+1}\left({\phi}_{\mathrm{S}}\right)\mathrm{sin}\left[2\pi \left(2n+1\right)\Omega t\right]\right\}\right),\hfill \end{array}$$## 10

[TeX:] \documentclass[12pt]{minimal}\begin{document}\begin{eqnarray} i_2 (z,t) &=& 2\rho \sqrt {I_{\rm R} I_{\rm S} } \left(\sin \phi _{{\rm DC}} \left[ J_0 \left({\phi _S } \right)\vphantom{\sum\limits_{n = 1}^\infty}\right.\right.\nonumber\\ &&\left.\left. +\, 2\sum\limits_{n = 1}^\infty {J_{2n} \left({\phi _S } \right)\cos \left({4\pi n\Omega t} \right)} \right] \right.\nonumber\\ &&\left. { +\, \cos \phi _{{\rm DC}} \left\{ {2\sum\limits_{n = 0}^\infty {J_{2n + 1} \left({\phi _{\rm S} } \right)\sin \left[ {2\pi \left({2n + 1} \right)\Omega t} \right]} } \right\}} \right),\hskip-6pt\nonumber\\ \end{eqnarray}\end{document} $$\begin{array}{ccc}\hfill {i}_{2}(z,t)& =& 2\rho \sqrt{{I}_{\mathrm{R}}{I}_{\mathrm{S}}}\left(\mathrm{sin}{\phi}_{\mathrm{DC}}\left[{J}_{0}\left({\phi}_{S}\right)\phantom{\sum _{n=1}^{\infty}}\right.\right.\hfill \\ & & \left.+\phantom{\rule{0.16em}{0ex}}2\sum _{n=1}^{\infty}{J}_{2n}\left({\phi}_{S}\right)\mathrm{cos}\left(4\pi n\Omega t\right)\right]\hfill \\ & & \left.+\phantom{\rule{0.16em}{0ex}}\mathrm{cos}{\phi}_{\mathrm{DC}}\left\{2\sum _{n=0}^{\infty}{J}_{2n+1}\left({\phi}_{\mathrm{S}}\right)\mathrm{sin}\left[2\pi \left(2n+1\right)\Omega t\right]\right\}\right),\hfill \end{array}$$*J*

_{n}is the

*n*’th-order Bessel function, Ω is the frequency of the harmonic excitation,

*I*

_{S}and

*I*

_{R}are the sample and reference optical intensities, respectively, ϕ

_{DC}is the quasi-static interferometric phase governed by the mean axial position of the scatterer, ρ is the detector responsivity, and ϕ

_{S}is expressed as

## 11

[TeX:] \documentclass[12pt]{minimal}\begin{document}\begin{equation} \phi _{\rm S} = \frac{{4\pi }}{\lambda }d\left({z,\Omega } \right), \end{equation}\end{document} $${\phi}_{\mathrm{S}}=\frac{4\pi}{\lambda}d\left(z,\Omega \right),$$*d*is the local, generally frequency-dependent, axial displacement amplitude so the mapping of

*d*represents the OCE displacement elastogram. It can be seen that

*d*can be determined from an interferometric measurement of ϕ

_{S}, which can in turn be achieved through the piecewise mapping of the measured ratios of Bessel harmonic terms to the vibration amplitude based on Eqs. 9, 10. This method however also involves extensive calculations to extract the vibration harmonics and is prone to noise. Moreover, elastograms of elastic moduli must be calculated from the obtained displacement based on the assumption of uniformly applied stress or using inverse problem-solving methods.

Because dynamic OCE is an active area of research, several additional techniques are in the beginning stages of feasibility demonstration. Huang
^{69} used a pressure-controlled air jet as an indenter to compress the soft tissue from the same side of imaging, in a noncontact manner and utilize the OCT signals to extract the deformation. Kennedy
^{70} have designed an appropriate sample arm based on a ring actuator and also enabled excitation and imaging to be achieved from the same side of the sample. Oldenburg and Boppart^{71} have proposed a method using embedded magnetomotive nanotransducers to excite soft tissues. However, experimental and computational techniques are still needed to further improve all of these techniques. Because of acquisition speed and the need for coaxial optical beam and vibration excitation, the *in vivo* application of dynamic OCE are limited to exposed skin, such as human hands.

As previously mentioned, OCE shows great advantage over ultrasound elastography on vascular disease diagnosis. However, for applications of OCE to vascular tissue characterization, special processing methodologies are usually required.

## 3.

## Intravascular Optical Coherence Elastography

## 4.1.

### Variational Method

Chan
^{72} proposed a variational method to improve deformation and strain estimation in OCE. The variational method takes a data-driven approach in which velocity ﬁltering takes place during the correlation maximization process itself. In this manner, all the available correlation data are taken into consideration in the robust velocity estimation. Speciﬁcally, the variational method incorporates prior knowledge about velocity ﬁelds in the pulsating arterial wall, where the overall variational energy functional is

## 12

[TeX:] \documentclass[12pt]{minimal}\begin{document}\begin{eqnarray} E\left[ {V(x,y)} \right] = aE_{\rm D} \left[ {V(x,y)} \right] + bE_{\rm S} \left[ {V(x,y)} \right] + cE_{\rm I} \left[ {V(x,y)} \right]\hskip-6pt\nonumber\\ \end{eqnarray}\end{document} $$\begin{array}{c}\hfill E\left[V(x,y)\right]=a{E}_{\mathrm{D}}\left[V(x,y)\right]+b{E}_{\mathrm{S}}\left[V(x,y)\right]+c{E}_{\mathrm{I}}\left[V(x,y)\right]\end{array}$$*V*(

*x*,

*y*) is the velocity ﬁeld,

*E*

_{D}[

*V*(

*x*,

*y*)] is data ﬁdelity,

*E*

_{S}[

*V*(

*x*,

*y*)] is strain ﬁeld smoothness, and

*E*

_{I}[

*V*(

*x*,

*y*)] is arterial wall incompressibility. These three terms are combined based on a certain weight:

*a*,

*b*, and

*c*. When the overall variational energy is minimized, the desired velocity ﬁeld is obtained. Following the velocity estimation, the deformation matrix can then be calculated at every point within the OCE imaging volume.

The velocity estimation of an intact aortic segment was obtained experimentally by mounting the specimen horizontally and applying mechanical load (lateral direction) with a one-dimensional translation micromanipulator to approximate circumferential stretching. This method produced improved results when compared to the traditional cross-correlation method, but *a priori* knowledge about arterial wall behavior was required for robust estimation of tissue motion and strain. These constraints may limit the overall usefulness of this technique in a clinical scenario based on the inhomogeneity of human tissue.

In 2005, Khalil
^{40} improved their OCE technique with the introduction of the nonlinear least-squares method, which minimizes the difference between computed and measured mechanical responses of soft tissue to estimate elasticity ﬁelds of soft tissues. This was an important step toward gaining quantitative elasticity information from *in vivo* OCT imaging of coronary arteries.

In 2006, the inverse elasticity problem approach^{73} was used comprising of FEM-derived mechanics directly within the image registration, which is the process of transforming the different sets of data into one coordinate system. Using a gradient-based search algorithm, tissue elasticity images were reconstructed given the internal tissue displacements based on subsequent imaging as the tissue is deformed. Results of *in vitro* OCE measurement of a porcine aorta injected with lipid are shown in Fig. 4. A fatty region was identiﬁed at the center of the arterial wall, showing a higher contrast in both the elastogram and the modulus image, further demonstrating the potential of OCE in vascular disease diagnosis by characterizing atherosclerotic plaque.

## 4.2.

### Optical Coherence Tomography Alternating-Line Elastography

To measure the radial displacement of vascular tissue, van Soest
^{74, 75} developed a catheter-based OCT alternating-line elastography (OCTALE) technique. It employs the correlation between the overlapping point-spread functions of adjacent A-scans, instead of subsequent B-mode imaging. The B-mode data are then separated into the odd and even A-scans, where acoustic forces are applied synchronously with the A-scan rate of the OCT system. This alternated sequence of A-scans, correlating to a low- or high-pressure scenario, yields odd and even displacement fields. The OCE image can be calculated by optimizing a one-dimensional cross-correlation function^{75}

## 13

[TeX:] \documentclass[12pt]{minimal}\begin{document}\begin{eqnarray} &&\hskip-6pt\rho _i \left(\delta \right)\nonumber\\ &&\hskip3pt = \frac{{\sum {\left[ {I_b \left({r_i } \right) - \left\langle {I_b \left({r_i } \right)} \right\rangle } \right] \cdot \left[ {I_{\rm s} \left({r_i + \delta } \right) - \left\langle {I_{\rm s} \left({r_i + \delta } \right)} \right\rangle } \right]} }}{{\sqrt {\sum {\left[ {I_b \left({r_i } \right) - \left\langle {I_b \left({r_i } \right)} \right\rangle } \right]^2 \cdot \sum {\left[ {I_{\rm s} \left({r_i + \delta } \right) - \left\langle {I_{\rm s} \left({r_i + \delta } \right)} \right\rangle } \right]^2 } } } }},\nonumber\\ \end{eqnarray}\end{document} $$\begin{array}{ccc}& & {\rho}_{i}\left(\delta \right)\hfill \\ & & \phantom{\rule{3.0pt}{0ex}}=\frac{\sum \left[{I}_{b}\left({r}_{i}\right)-\u2329{I}_{b}\left({r}_{i}\right)\u232a\right]\xb7\left[{I}_{\mathrm{s}}\left({r}_{i}+\delta \right)-\u2329{I}_{\mathrm{s}}\left({r}_{i}+\delta \right)\u232a\right]}{\sqrt{\sum {\left[{I}_{b}\left({r}_{i}\right)-\u2329{I}_{b}\left({r}_{i}\right)\u232a\right]}^{2}\xb7\sum {\left[{I}_{\mathrm{s}}\left({r}_{i}+\delta \right)-\u2329{I}_{\mathrm{s}}\left({r}_{i}+\delta \right)\u232a\right]}^{2}}},\hfill \end{array}$$*m*-lines by

*n*pixels. The tissue displacement is obtained when the cross-correlation coefficient ρ attains its maximum. Radial strain can then be computed from the displacement via the following equation:

## 14

[TeX:] \documentclass[12pt]{minimal}\begin{document}\begin{eqnarray} \varepsilon _{\rm r} \left({r,\theta } \right) &=& \frac{d}{{dr}}\Delta \left({r,\theta } \right) \Rightarrow \varepsilon _{\rm r},ij\nonumber\\ & =& \frac{1}{{r_{{\rm step}} }}[ {\Delta ({r_i,\theta _j }) - \Delta({r_{i + j},\theta _j })}], \end{eqnarray}\end{document} $$\begin{array}{ccc}\hfill {\varepsilon}_{\mathrm{r}}\left(r,\theta \right)& =& \frac{d}{dr}\Delta \left(r,\theta \right)\Rightarrow {\varepsilon}_{\mathrm{r}},ij\hfill \\ & =& \frac{1}{{r}_{\mathrm{step}}}\left[\Delta \left({r}_{i},{\theta}_{j}\right)-\Delta \left({r}_{i+j},{\theta}_{j}\right)\right],\hfill \end{array}$$*r*

_{step}is the radial distance in pixels between the centers of subsequent windows, Δ is the displacement. The performance of the ALE algorithm applied to a simulated displacement/strain field is shown in Fig. 5. The OCT image of a vessel phantom used for simulation is shown in Fig. 5f. The displacement field as estimated by ALE is shown in Fig. 5a, which can be compared to the simulated input, displayed in Fig. 5d. Strain field calculated based on ALE is shown in Fig. 5b, which can be compared to the simulated strain field in Fig. 5e. The peak correlation value is shown in Fig. 5c. Comparing the pairs [Figs. 5a and 5d to 5b and 5e], one can see that the displacement and strain field are well reconstructed using the ALE algorithm.

The OCTALE method is less sensitive to motion artifacts and nonuniform catheter rotation than the traditional cross-correlation-based method, but it is difficult to apply in practice. The applied localized force to the vascular wall at a rate of exactly half the OCT line-scan frequency must be applied, and the frequency of the applied pressure oscillation must be locked to the instrument's line rate, which brings many requirements for the load and control systems design. However, ALE expands the potential of OCT as a powerful tool for intracoronary vulnerable plaque detection.

There are still challenges faced in deployment of OCE for *in vivo* catheter-based imaging. The cross-correlation between frames is not suitable for cardiovascular applications because the displacement and strain patterns induced by a change in ﬂuid pressure in a tubular geometry differ from those resulting from a mechanical force applied in a planar geometry. Additionally, when one operates an OCE probe inside a living patient, while imaging through a rotating catheter, the images and resulting stress/strain maps are more susceptible to artifacts due to motion and optical disturbances. Arteries are typically characterized by asymmetric geometries, which together with the presence of plaques with varying compositions lead to inhomogeneous stress ﬁelds in the arterial wall. Furthermore, arterial tissues exhibit a nonlinear behavior and experience high cyclic strain levels in each single cycle of a heart beat. These strain fields can vary up to 10% at the common carotid artery or carotid bifurcation.^{76, 77} Viscoelastic models might be needed for dynamic analysis, especially for the lipid pool component of the plaques to fully characterize tissue in an effort to predict which plaques are truly vulnerable to stratify appropriate treatments or targeted cardiovascular therapies.

## 4.

## Discussion

The purpose of imaging elastic tissue properties is to support the detection, localization, and differential diagnosis of lesions. The quality of elastograms is determined in terms of resolution, dynamic range, and signal-to-noise ratios. These attributes are largely dependant on the ability to obtain quality estimations of tissue strain. For all elastography techniques, the resolution is determined by the underlying imaging modality. OCT has superior spatial resolution (1−10 μm, in tissue), when compared to other imaging modalities, to a depth of a few millimeters resulting in a high-resolution OCE technique. The ability to provide high-resolution imaging *in vivo* is one of the most signiﬁcant advances attributable to OCT. With even higher resolution and faster data aquisition techniques on the horizon and simpliﬁcation of the procedural requirements, OCT will become more accessible to a greater number of clinical centers worldwide.^{78}

Elastograms provide additional information of stiffness of soft tissues to OCT images, which can be especially important in medicine for early diagnosis of diseases before any morphological symptoms occur. Different schemes need to be designed for different measurement purposes. For example, breast elastography typically deals with a large single nodule that is usually <10 times stiffer than the normal breast tissue and static compression forces are often used.^{10} However, arterial elastography is more challenging than breast-tissue elastography due to the smaller length scales involved, time-varying blood ﬂow, and presence of multiple plaques with different sizes, geometries, and compositions. Modulus of lipid pools in atherosclerotic plaques can alone span four orders of magnitude.^{79} Moreover, a time-varying pressure is most commonly exerted on the arterial plaques.^{47}

Currently, there are still the following hurdles to overcome before the widespread adoption of OCE for clinically relevant *in vivo* applications: (*i*) The basic assumption for correlation-based approaches is that the speckle in the OCT image produced from the coherence volume is stable, except for a lateral translation between successive B-scans. This assumption limits these techniques to detecting motions between the successive B-scans less than the spatial coherence length of the light source used. For biomedical diagnostics, there may be potential sources of error and speckle decorrelation associated with substantial speckle pattern shifts during deformation of the soft tissues. These characteristics contribute to high noise components and must be addressed prior to the development of a robust elastography technique. (*ii*) The application of controlled transmural pressures is required to make accurate OCE measurements *in vivo*. This will require the development of new OCT probes and can potentially be achieved with an intravascular pressure balloon or calibrated saline ﬂushes. (*iii*) Kinematic and mechanical models suitable for target tissues are required. The approximation of certain tissue types being elastic and incompressible is reasonable. However, diverse soft tissues do not exhibit the linear stress-strain response. Therefore, nonlinear theory of material behavior must be implemented in order to quantify soft tissues, such as atherosclerotic or severely hypertensive arteries, where the tensile properties will not necessarily be completely uniform. (*iv*) Well-designed and precisely controlled experimental methods are required to investigate the anisotropic behavior of biomaterials in an effort to reduce imaging artifacts. (*v*) Sophisticated numerical methods, including ﬁnite element method and inverse problem method to estimate material parameters, must be continually improved to provide a model that can accurately describe specific *in vivo* environments. (*vi*) The packaging of such devices for intraluminal or internal use remains to be a critical engineering problem that must be solved before multicenter clinical testing can be performed as a first step to widespread clinical adoption of OCE. It is important to note that although these technical hurdles exist, engineering solutions are possible in the near term.

A large number of studies have demonstrated the potential of using OCT to investigate and facilitate the diagnosis of disease, direct therapeutic intervention, and image *in vivo* arterial pathology in humans.^{56, 80, 81, 82} Beneﬁting from high-speed *in vivo* OCT imaging, OCE system also has great potential to measure 3-D biomechanical properties of tissues and conduct *in vivo* studies. With cellular-level resolution and several millimeters of imaging penetration into highly scattering tissues, OCE has the unique ability to noninvasively measure *in vivo* tissue biomechanical properties. Possible clinical applications of OCT elastography include differentiation of hard and soft masses during tissue biopsy, imaging of arterial plaque composition, and evaluation of wound healing. OCT elastography may also aid efforts to solve fundamental problems in biomechanics, such as the role of microscopic deformations in the development of the embryo.^{48}

OCE may play a major role in the biomedical imaging field specifically in intravascular and atherosclerotic tissue imaging to help the characterization of tissue biomechanical properties. Mechanical properties of vascular tissue could provide diagnostic information about a range of vascular diseases, from pulmonary hypertension to coronary atherosclerosis. The *in vivo* intracoronary OCE would augment the high-resolution imaging capability within formation on mechanical properties of arterial walls. Monitoring of atherosclerosis is vital for determining the effectiveness of pharmaceutical treatments that try to reduce the rupture proneness of an artery (e.g., by reducing the amount of lipid or by stiffening the lipid). A large stiffness variation near the lumen is generally assumed mechanically unstable.^{83} Therefore, the creation of an arterial Young's modulus image through OCE might be possible to establish a diagnostic tool for assessment of plaque vulnerability.

Because of the advantage of high resolution of OCE, additional applications are expected to be explored in areas such as tissue engineering to interpret microstructural reorganization of connective tissues during cancer progression from the macroscopic deformation patterns in viscoelastic images, molecular, and cellular biology,^{32} and shear stress over smooth muscle cells in a deformable arterial wall.^{84} Another application may be the material property inspection of aneurysms.^{85} The majority of aneurysms form in the cerebral arteries and the abdominal aorta. A cerebral aneurysm,^{86} as a blood-ﬁlled dilation (balloonlike bulge) in the wall of the cerebral artery, can cause devastating clinical consequences. It is widely thought that aneurysms rupture when hemodynamically induced wall stress exceeds wall strength. Hence, OCE may be very helpful for biomechanical analysis of aneurysms.^{87} Identiﬁcation of the mechanism responsible for aneurysm rupture may be used as clinical predictors of the rupture. This area of research remains to be the focus of many studies as stroke continues to be a leading cause of morbidity and mortality in humans. As the OCE technique continues to evolve, it may elucidate a better understanding of both disease and injury and in so doing contribute to the advancement of health care.

## References

*in vivo*,” Eur. J. Ultrasound 3, 49–70 (1996). 10.1016/0929-8266(95)00134-4 Google Scholar

*in vivo*: demonstration of feasibility,” Phys. Med. Biol. 54, 1191–1207 (2009). 10.1088/0031-9155/54/5/006 Google Scholar

*In vivo*Doppler optical coherence tomography of mucocutaneous telangiectases in hereditary hemorrhagic telangiectasia,” GI Endosc. 58, 591–598 (2003). 10.1067/S0016-5107(03)01978-3 Google Scholar

*In vivo*dynamic optical coherence elastography using a ring actuator,” Opt. Express 17, 21762–21772 (2009). 10.1364/OE.17.021762 Google Scholar