|
The biomechanical properties of the cornea can be significantly altered by several ocular diseases, such as keratoconus.1 Hence, assessing corneal biomechanical properties are required for detecting the onset and monitoring progression of the ocular diseases as well as assessing the effects of therapeutic interventions. Previous studies have shown the feasibility of commercially available instruments such as the ocular response analyzer and the CorVis,2 and emerging techniques such as Brillouin microscopy3 for estimating corneal biomechanical properties. Although these methods can provide important measurements that can reflect some corneal biomechanical features, quantitative assessment of the corneal viscoelasticity is still a challenge. Optical coherence elastography (OCE) is an emerging technique with a great potential for noninvasive measurements of the local biomechanical properties of tissues with high spatial and temporal resolutions.4 Similar to other elastographic techniques, such as magnetic resonance elastography5 and ultrasound elastography,6 OCE combines the corresponding imaging modality with mechanical loading. We previously used OCE in different applications, such as soft tissue tumor detection,7 cornea,8 and cardiac muscle9 elasticity estimation. However, quantification of tissue viscoelastic properties from OCE measurements requires the selection of a proper model that can accurately map the parameters of measured elastic waves to, e.g., Young’s modulus. In this study, we demonstrate, for the first time to the best of our knowledge, application of a modified Rayleigh–Lamb frequency equation (RLFE) to quantify the viscoelastic properties of porcine corneas from OCE measurements. Based on the temporal displacement profiles of a focused air-pulse-induced elastic wave, phase velocities over a range of angular frequencies were calculated and fitted using the RLFE to extract the corneal elasticity and viscosity. A validation study was performed on agar phantoms before the RLFE method was applied to the porcine corneas. The home-built OCE system is composed of a focused air-pulse delivery system10 and a phase-stabilized swept source optical coherence tomography (PhS-SSOCT) system.11 The system utilized a broadband swept laser source (HSL2000, Santec, Inc., California) with a central wavelength of , bandwidth of , scan rate of 30 kHz, and experimentally measured phase stability of during cornea experiments. The focused air-pulse delivery system was comprised of a controller with a signal input for synchronization, a solenoid-controlled air gate, and an air-pulse port with a flat edge and diameter of . The system is capable of delivering a short duration focused air-pulse () with a Gaussian profile to induce a small amplitude deformation (order of ) within the cornea tissue. The air-pulse was positioned away from the apex of the cornea to ensure that the elastic wave propagated across the apex and away from the surface of the cornea. Imaging of the elastic wave propagation was performed by synchronizing 501 M-mode images at successive positions across tissue surface.12 The phase velocity of the elastic wave was calculated by performing a fast Fourier transform on the elastic wave temporal displacement profiles at each measurement position and imaged depth, producing a phase shift, , at a specific angular frequency, , and depth. For each depth, the distance, , of the elastic wave propagation was calculated, taking the curvature of the cornea into account.13 The phase velocity of the elastic wave at the angular frequency, , is expressed as: , which was obtained by least-squares linear fitting. The resulting phase velocities at each frequency were then averaged over the imaged depths in the sample. Validation experiments were first performed on 2% agar phantoms (, ) with a cylindrical shape: diameter and height . The phantoms were assumed to be homogeneous and isotropic with a negligible viscosity.14 The agar phantoms were placed on a ring support such that the central part of the bottom surface was free of the contact. The excitation position was at the center of the top surface of the phantom and all OCE measurements were taken in the central region. Assuming free boundary conditions, the Rayleigh-Lamb characteristic equation for an anti-symmetric wave has the form15 In Eq. (1), is half the thickness of the phantom (), and is the wave number at angular frequency . For Eq. (1), the following relationships can be expressed with where is the material density, and and are the Lamé constants with the relationships and , respectively. is the Young’s modulus and is a real number as the phantom viscosity is assumed to be negligible and is the Poisson’s ratio of the material. The compressional wave velocity and the shear wave velocity are and , respectively.By solving Eq. (1), we obtained the numerical relationship between and , which was then compared with the OCE experimental data to quantify the Young’s modulus. Elasticity values obtained from OCE measurements were then validated on agar phantoms (same concentration 2% , ) using uniaxial mechanical compressional tests (Model 5943, Instron Corp., Massachusetts). Figure 1 shows the agar phantom elasticity assessment by RLFE and validation by mechanical testing. Figure 1(a) presents a typical temporal displacement profile from the agar phantoms as measured by OCE. The duration of displacement [ as marked in Fig. 1(a)] was assumed as a half wavelength for determining the lower limit of frequencies to be used in the RLFE. For the agar phantoms, angular frequencies lower than were due to low frequency noise and were not used for calculations. Figure 1(b) depicts the phase velocities calculated by the RLFE and as measured by OCE. From the RLFE and phase velocity data, the Young’s modulus of the 2% agar phantoms was estimated at 160 kPa. Deviations of the OCE data from the RLFE may be due to the fact that the Lamb wave model is strictly based on a thin layer geometry, to which the phantoms do not strictly adhere. Figure 1(c) shows the Young’s moduli of the agar phantoms as estimated by the RLFE and measured by uniaxial mechanical compression testing. The RLFE results agreed well with the mechanical testing measurement, which demonstrates the feasibility of this method for quantitative elasticity assessment. Fig. 1(a) A typical temporal displacement profile of a single point on a 2% agar phantom. (b) Optical coherence elastography (OCE) measurements of the elastic wave dispersion for 2% agar phantom fitted with Rayleigh–Lamb frequency Eq. (1). (c) Quantitative results of the Young’s modulus of 2% agar phantom assessed by Rayleigh–Lamb frequency equation (RLFE) and mechanical testing. ![]() Next, experiments were performed on porcine corneas () in whole eyeball configuration with artificially controlled intraocular pressure, . The average thickness of the corneas was (de-epithelialized). The RLFE used for agar phantoms was modified because the cornea has two distinct features. First, the viscosity of the cornea must be taken into account. We assumed a Kelvin–Voigt viscoelastic model, which led to a complex dynamic shear modulus of the form: , where is the shear viscosity and is the imaginary unit. Second, the fluid effect caused by the aqueous humor must be considered. Although the cornea has a convex shape, it was treated as a flat disk for simplification (the effect of the corneal curvature on the estimated values of and is currently the topic of our investigation). In addition, the boundary conditions were assumed as zero stress on the corneal anterior surface and equal stresses and vertical displacements between the corneal tissue and aqueous humor on the posterior surface. For this case, the Lamb wave cannot be considered as symmetric or anti-symmetric. Based on these assumptions, the RLFE for the cornea was then modified as with andHere, is half of the porcine corneal thickness, is the fluid density, is the cornea density,16 and is the speed of sound in water. Figure 2 presents the corneal viscoelasticity assessment by the modified RLFE from OCE measurements. A typical temporal displacement profile as measured by OCE from a porcine cornea is shown in Fig. 2(a). Similar to the agar phantoms, the displacement duration of marked in Fig. 2(a) was used to determine the lower frequency limit of to be adopted in the modified RLFE. Figure 2(b) shows the estimated Young’s modulus () and shear viscosity () of three porcine corneas at an artificially controlled IOP of 20 mmHg, which were estimated by fitting the phase velocities from the OCE measurements with the modified RLFE numerical results. The estimated Young’s modulus is of the same order as measured by atomic force microscopy,17 which suggests that the proposed method could provide accurate quantitative viscoelastic characterization of the cornea under the proper assumptions. Fig. 2(a) A typical temporal displacement profile from a porcine cornea. (b) Rayleigh–Lamb wave dispersion for porcine corneas () fitted with modified Rayleigh–Lamb frequency Eq. (4), with estimated Young’s modulus () and shear viscosity (). ![]() To summarize, in this study, the viscoelasticity of porcine corneas was quantitatively assessed by fitting the phase velocities of an elastic wave measured by PhS-SSOCE to a modified Rayleigh-Lamb wave model, which accounted for the fluid effect on the corneal posterior surface. RLFE validation was conducted on 2% agar phantoms by uniaxial mechanical compression testing, which verified the accuracy of this method. The Young’s modulus of porcine corneas at 20 mmHg IOP was estimated to be and the shear viscosity as . The results indicate that the combination of PhS-SSOCE and the modified Rayleigh-Lamb characteristic equation may be potentially useful for assessing corneal viscoelasticity in vivo. AcknowledgmentsThis work was supported, in part, by Grant Nos. 1R01EY022362, 1R01HL120140, and U54HG006348 from the NIH and PRJ71TN from DOD/NAVSEA. ReferencesZ. L. Hanet al.,
“Biomechanical and refractive behaviors of keratoconic cornea based on three-dimensional anisotropic hyperelastic models,”
J. Refract. Surg., 29
(4), 282
–290
(2013). http://dx.doi.org/10.3928/1081597X-20130318-08 JRSUEY 0883-0444 Google Scholar
Z. L. Hanet al.,
“Air puff induced corneal vibrations: theoretical simulations and clinical observations,”
J. Refract. Surg., 30
(3), 208
–213
(2014). http://dx.doi.org/10.3928/1081597X-20140212-02 JRSUEY 0883-0444 Google Scholar
G. Scarcelliet al.,
“Brillouin microscopy of collagen crosslinking: noncontact depth-dependent analysis of corneal elastic modulus,”
Invest. Ophthalmol. Vis. Sci., 54
(2), 1418
–1425
(2013). http://dx.doi.org/10.1167/iovs.12-11387 IOVSDA 0146-0404 Google Scholar
S. WangK. V. Larin,
“Optical coherence elastography for tissue characterization: a review,”
J. Biophoton.,
(2014). http://dx.doi.org/10.1002/jbio.201400108 JBOIBX 1864-063X Google Scholar
R. Muthupillaiet al.,
“Magnetic resonance elastography by direct visualization of propagating acoustic strain waves,”
Science, 269
(5232), 1854
–1857
(1995). http://dx.doi.org/10.1126/science.7569924 SCIEAS 0036-8075 Google Scholar
J. Ophiret al.,
“Elastography: a quantitative method for imaging the elasticity of biological tissues,”
Ultrason. Imaging, 13
(2), 111
–134
(1991). http://dx.doi.org/10.1177/016173469101300201 ULIMD4 0161-7346 Google Scholar
S. Wanget al.,
“Computational analysis of optical coherence tomography images for the detection of soft tissue sarcomas,”
Proc. SPIE, 8580 85800T
(2013). http://dx.doi.org/10.1117/12.2006638 PSISDG 0277-786X Google Scholar
J. Liet al.,
“Air-pulse OCE for assessment of age-related changes in mouse cornea in vivo,”
Laser Phys. Lett., 11 065601
(2014). http://dx.doi.org/10.1088/1612-2011/11/6/065601 1612-2011 Google Scholar
S. Wanget al.,
“Noncontact quantitative biomechanical characterization of cardiac muscle using shear wave imaging optical coherence tomography,”
Biomed. Opt. Express, 5
(7), 1980
–1992
(2014). http://dx.doi.org/10.1364/BOE.5.001980 BOEICL 2156-7085 Google Scholar
S. Wanget al.,
“A focused air-pulse system for optical-coherence-tomography-based measurements of tissue elasticity,”
Laser Phys. Lett., 10
(7), 075605
(2013). http://dx.doi.org/10.1088/1612-2011/10/7/075605 1612-2011 Google Scholar
R. ManapuramV. ManneK. Larin,
“Development of phase-stabilized swept-source OCT for the ultrasensitive quantification of microbubbles,”
Laser Phys., 18
(9), 1080
–1086
(2008). http://dx.doi.org/10.1134/S1054660X08090144 LAPHEJ 1054-660X Google Scholar
S. WangK. V. Larin,
“Shear wave imaging optical coherence tomography (SWI-OCT) for ocular tissue biomechanics,”
Opt. Lett., 39
(1), 41
–44
(2014). http://dx.doi.org/10.1364/OL.39.000041 OPLEDP 0146-9592 Google Scholar
S. WangK. V. Larin,
“Noncontact depth-resolved micro-scale optical coherence elastography of the cornea,”
Biomed. Opt. Express, 5
(11), 3807
–3821
(2014). http://dx.doi.org/10.1364/BOE.5.003807 BOEICL 2156-7085 Google Scholar
A. Ahmadet al.,
“Magnetomotive optical coherence elastography using magnetic particles to induce mechanical waves,”
Biomed. Opt. Express, 5
(7), 2349
–2361
(2014). http://dx.doi.org/10.1364/BOE.5.002349 BOEICL 2156-7085 Google Scholar
K. F. Graff, Wave Motion in Elastic Solids, Courier Dover Publications, Mineola, New York
(1975). Google Scholar
J. Kampmeieret al.,
“Thermal and biomechanical parameters of porcine cornea,”
Cornea, 19
(3), 355
–363
(2000). http://dx.doi.org/10.1097/00003226-200005000-00020 CORNDB 0277-3740 Google Scholar
J. Seifertet al.,
“Distribution of Young’s modulus in porcine corneas after riboflavin/UVA-induced collagen cross-linking as measured by atomic force microscopy,”
Plos One, 9
(1), e88186
(2014). http://dx.doi.org/10.1371/journal.pone.0088186 1932-6203 Google Scholar
|