Secondary Logo

Journal Logo

Strain Imaging of Corneal Tissue With an Ultrasound Elasticity Microscope

Hollman, Kyle W. Ph.D.; Emelianov, Stanislav Y. Ph.D.; Neiss, Jason H. B.S.; Jotyan, Gagik Ph.D.; Spooner, Gregory J.R. Ph.D.; Juhasz, Tibor Ph.D.; Kurtz, Ron M. M.D.; O'Donnell, Matthew Ph.D.

Basic Investigations
Free

Purpose. We hypothesize that high-resolution elasticity measurements can guide corrective refractive surgery of the cornea. Elasticity measurements would improve surgical outcomes by adding biomechanical information not used in existing clinical nomograms. As an initial investigation, we determined the usefulness and evaluated the ability of our ultrasound elasticity microscope by measuring strain ex vivo in an intact porcine eye globe.

Methods. Strain was predicted with a finite element model guided by direct mechanical measurements of corneal elasticity. Next, a porcine cornea was deformed with a slitted plate while being imaged with ultrasound. For high spatial resolution, the ultrasound elasticity microscope uses a 50 MHz transducer with a 1.4 f/number. It produces high-quality conventional ultrasonic B-scans over large thicknesses by confocal processing. Strain was calculated from tracking speckle in these images after deformation. This technique is compatible with in vivo measurements.

Results. Compressional and expansional deformations were the same order of magnitude from −3.5% to as great as +3.5%. Strain imaging indicated the stroma expanded into the slit of the deformation plate while Bowman's layer compressed. This bipolar variation within a specimen is unusual. Within the stroma, a variation of strain with depth was measured suggesting a distribution of elasticity. Results compared favorably with the finite element model.

Conclusion. An ultrasound elasticity microscope can produce high-resolution strain images throughout the corneal depth. Various layers with different elastic properties appeared as different strains in the images.

From the Departments of Biomedical Engineering (K.W.H., S.Y.E., J.H.N., G.J.R.S., T.J., M.O.) and Ophthalmology (R.M.K.) and the Center for Ultrafast Optical Science (G.J.R.S., T.J.), University of Michigan, Ann Arbor, Michigan, U.S.A.; and the Research Institute for Particle and Nuclear Physics (G.J.), Budapest, Hungary.

Submitted March 28, 2001.

Revision received August 16, 2001.

Accepted August 22, 2001.

This project was supported in part by NIH grants DK47324, HL47401, and 2r444EY12340-02.

Address correspondence and reprint requests to Dr. K.W. Hollman, Department of Biomedical Engineering, University of Michigan, Room 3304 G.G. Brown, 2350 Hayward, Ann Arbor, MI 48109-2125, U.S.A.

Conventional imaging with high-frequency ultrasound is now used in ophthalmology for disease detection and biometry. 1,2 We hypothesize that ultrasound can play a much greater role. In particular, an ultrasound elasticity microscope has the potential to produce spatial maps of corneal elasticity. Such images could provide valuable information to guide the development and eventual adaptation of new refractive surgical procedures 3,4 and to improve the reliability of current procedures.

As the most common reason for ophthalmic evaluation, refractive error accounts for 90% of patient visits. In the past decade, refractive surgical techniques have met with considerable success and popularity, resulting in tremendous growth. Just in the United States, laser-assisted in situ keratomileusis (LASIK) is expected to reach 1.5 million procedures annually. 5 Although the complication rates for LASIK are low, 6,7 steep increases in the number of procedures performed will expose greater numbers of patients to risk. Even in successful cases, nonphysiologic optical surfaces created by a laser can lead to optical aberrations, loss of contrast sensitivity, and other visual side effects. 8

In LASIK, a mechanical blade, a microkeratome, cuts a flap of anterior corneal tissue to expose deeper stromal layers for subsequent excimer laser ablation. Because of the lamellar structure of corneas, fibrils cut during tissue removal redistribute their load to uncut fibrils. Corneal shape relaxes to a new equilibrium state to reflect this change in loading. The biomechanical response, along with changes in mechanical properties of corneal tissue, partially determines final corneal curvature. Alterations of these properties may contribute to medical complications, such as overcorrections, undercorrections, and ectasia. 9

To ensure predictability and safety for refractive surgery, accurate biomechanical modeling of corneal structure is needed. Advanced models of corneas require spatial distributions of material parameters, particularly the elastic (Young's) modulus. 10 Currently, elastic modulus data are only available from ex vivo studies. 11 Ultrasound imaging can provide individual in vivo measurements before, during, and after surgery. Consequently, ultrasound elasticity can measure important biomechanical parameters driving mechanical models of corneas. Model outputs may then be used to adapt and refine individual surgical procedures in real time for optimal outcomes.

Ultrasound does not image tissue elasticity (shear or Young's modulus) directly because the contrast mechanism is primarily related to changes in bulk modulus, rather than shear elastic modulus. Ultrasound imaging systems can indirectly estimate and image elastic properties of tissue using ultrasonic measurements of externally or internally induced internal tissue motion. The ultrasound elasticity microscope was developed to image tissue stiffness at high resolution. Tissue mimicking phantoms 12 and tissue engineered specimens 13 have been imaged with this instrument. Its characteristics and function are described in more detail elsewhere. 14

To test whether ultrasound elasticity imaging can provide useful information to drive mechanical models of corneal structures, we have adapted an existing elasticity microscope to produce simultaneous B-scan and elasticity images of the cornea. Using B-scan and strain images as qualitative measures, this system has been adjusted until it produced images similar to those published previously. 12,13 We report the results of initial investigations with this system on corneal elasticity imaging of an ex vivo, intact, porcine eye globe. Hypotheses tested on this experimental system have direct application to an eventual real-time, clinical system. This experiment is a critical initial test in the development of the ultrasound elasticity microscope for aiding refractive correction procedures of corneas.

Back to Top | Article Outline

MATERIALS AND METHODS

Ultrasound Elasticity Microscope

A diagram of the experimental setup is shown in Figure 1. Details of the scanning and ultrasound system are similar to those described elsewhere. 12,14 An intact porcine eye globe was embedded in stiff gelatin up to the limbus to simulate the natural orbital environment. The cornea was deformed by a flat plate with a narrow slit machined in it. For controlled deformation, the plate was coupled by an arm to a precision sled and stepper motor. Poor acoustic impedance matching between water and the metal plate would cause most of the ultrasonic signal to be reflected from this boundary. To avoid this signal loss, the transducer images through a 2.3-mm wide slit. To avoid reflection of side lobes from either side of the slit, 1-mm wide images were acquired in the center of the slit. Perpendicular to the imaging plane, the slit is nearly 20 mm long. This geometry reduces out-of-plane motion. 12,14Figure 2 shows corneal deformation. Tissue in contact with the plate undergoes compressional deformation, whereas tissue in the middle of the slit bulges into it.

FIG. 1.

FIG. 1.

FIG. 2.

FIG. 2.

To image strain, confocal images of the cornea were taken before and after deformation. Before measurement, the cornea was in slight contact with the deformation plate. Two surface deformation steps of 50 μm were then applied, capturing ultrasound data at no deformation and at each step. Speckle was tracked between these images with a correlation-based algorithm. 15–17 Tissue motion was computed from speckle tracking, and strains were determined by numerically differentiating displacement components.

To maintain high resolution throughout the cornea, we imaged with a 50-MHz center frequency, single element transducer. Its depth of field was 260 μm, but typical corneal thicknesses are greater than 1 mm. Speckle tracking outside the focal zone of a fixed focus transducer is poor because of speckle decorrelation, so confocal image processing was used. In this process, the transducer is scanned at multiple axial positions, and resulting signals in the depth of field are spliced together guided by correlation. 18,19 Without confocal imaging, determining strain through the entire thickness of a cornea would be difficult.

With an f/number of 1.4 and a quality factor of 2.25, the transducer has a theoretical resolution of 46 μm laterally and 37 μm axially for conventional B-scans. In previous experiments 12 with the same system, axial resolution was measured to be 52 μm for displacement and 71 μm for strain. Resolution is slightly reduced at each step in processing, because displacement is calculated from several points in a B-scan and strain is calculated from multiple displacements.

Back to Top | Article Outline

Direct Mechanical Elasticity and Finite Element Model

To corroborate strain images obtained with this instrument, we created a finite element simulation based on parameters obtained from direct mechanical measurements. A mechanical technique was used to measure elasticity distributions in large tissue samples. 20 This technique produces nearly direct measures of elasticity and is widely accepted. 21,22 Using a femtosecond corneal surgical laser, 4,23 five cuts were made at increasing corneal depths to produce a series of circular disks approximately 200 μm thick and 7 mm in diameter. In contrast to microkeratome cuts in which the slices have complex curvatures and variable thickness, femtosecond laser cutting can create custom geometries. To produce disks for mechanical measurements, a plane was cut parallel to the applanated anterior corneal surface. If these disks assume natural curvatures in the intact globe, resulting anterior and posterior surfaces would no longer be planar when the applanating contact lens is removed after resection. In practice, an unloaded tissue block has insufficient rigidity to revert to its loaded curvature and assumes a planar shape through surface tension when placed on a flat surface. The resulting disk sections have uniform thickness with parallel, planar surfaces. 24 Each sample was positioned horizontally on the rigid platform of a scale. Deformations were applied with an 8-mm diameter piston covering the entire surface of the sample. The measured Young's elastic modulus as a function of tissue depth is shown in Figure 3. Elasticity variations with depth are shown for 10%, 20%, 30%, and 40% deformations. Because of the difficulty in cutting into thick porcine corneas, we were only able to obtain samples for two thirds of the thickness. No measurements were obtained for the endothelium or Descemet's membrane. In direct mechanical measurements, spatial resolution is determined by the section depth and diameter. Unlike direct mechanical measurements, ultrasound based elasticity has better axial and lateral resolution, and it can measure strain in vivo.

FIG. 3.

FIG. 3.

Mechanical measurements produce a direct measure of elasticity, a material property independent of geometry. Axial strain depends on geometry, 13,25 so a finite element model is needed to convert elasticity from mechanical measurements into axial strain to compare with ultrasonic measurements. Finite element simulations describing static deformations of a linear elastic mechanical body induced by small externally applied deformation were performed using a commercially available finite element package (ABAQUS, HKS Inc., Pawtucket, RI). Corneal thickness was set to 1.30 mm as determined from the ultrasonic B-scan. Because exact elasticity (Young's modulus) distribution in the entire cornea was not known a priori, the analysis was performed for several cases. In case 1, the cornea was modeled as an object with two layers, where the anterior 160 μm layer was 2.2 times stiffer than the background tissue. In case 2, an additional 130 μm thick layer was introduced in the posterior. This layer was 1.8 times stiffer than the corneal stroma. Mechanical properties of all layers were guided by the direct mechanical measurements. In case 3, anterior and posterior layers in the cornea were 4.4 and 3.6 times stiffer than the background tissue. Although case 3 doubles anterior and posterior stiffness from case 2, mechanical measurements suggest these elastic stiffness values are plausible. In all cases, materials were assumed to be nearly incompressible, that is, having a Poisson ratio of 0.5.

A plane strain approximation of the deformation experiment was used in simulations. Such approximations are reasonable for our experimental conditions and sufficient for qualitative analysis and comparison. The model cornea was deformed from the top by downward motion of a rigid plate with a 2.3-mm gap located in the center of symmetry. The bottom boundary of the cornea was constrained in the vertical direction, and horizontal slippage was allowed at the bottom and top boundaries. Side boundaries of 25-mm long vertical segments were unconstrained. Also, the 2.3-mm gap of the top boundary was unconstrained. A Cartesian spatial grid with nonuniform sampling was used in simulations. Underneath the slit, in the imaging region, a 25-μm grid was used, and elsewhere a larger grid was sufficient. Further decrease of the grid spacing did not produce significant change and was unnecessary.

Back to Top | Article Outline

RESULTS

Results of the finite element models are shown in Figure 4. Positive strain represents expansional deformation, and negative strain designates compressional deformation. For case 1, the model predicts that the stroma, Descemet's membrane, and endothelial layer will expansionally deform and that the epithelium and Bowman's layer will generally compressionally deform under this load. For case 2, the endothelial and Descemet's layers will have significantly reduced deformation compared with case 1. Case 3 behaves similar to case 2, except there is greater deformation in the stroma. In all three cases, the model predicts bipolar strain (i.e., expansional and compressional deformation in the same sample under the conditions of our experiment). Such results are unusual, but reasonable, considering the layered structure of corneas and the mechanical response of the layers.

FIG. 4.

FIG. 4.

Results of ultrasonic measurements on a porcine cornea are shown in Figure 5. In the B-scan (Fig. 5A), the transducer is positioned at the top of the image, pointing toward page bottom. The top bright echo indicates the boundary between the coupling medium (water) and the epithelial surface. Below this, a second bright echo indicates the boundary of Bowman's layer. Lower yet, stromal tissue appears as uniform speckle. At the bottom of the image, Descemet's membrane and the endothelium appear as bright echoes. Commonly, fluid in the anterior chamber does not scatter ultrasound in this frequency range, so it appears as a region with no signal. For accurate strain imaging, a high fidelity conventional B-scan image is needed because strain imaging requires spatial differentiation of the B-scan. Our image compares favorably with published images of standard ultrasound techniques in corneas. 1,2,26

FIG. 5.

FIG. 5.

Figure 5B shows a strain image in which positive and negative strains are mapped. Mid-gray to white indicates positive strain (expansional deformation). Mid-gray to black indicates negative strain (compressional deformation). As we follow strain from the middle epithelium to the anterior stroma, it changes rapidly from compressional to expansional deformation in a short distance. This result showing bipolar strain is unusual but within reason and is predicted by the finite element model.

To emphasize the behavior of corneal layers, the strain image was split in two. Strain images were overlaid on top of the B-scan to reference strain to corneal structure. Figure 6A emphasizes stromal behavior. In this image, strain ranges from + 0.2% to + 3.5%. High strains are presented in yellow and white and lower strains in red. This figure indicates general bulging in the stroma, but it is not uniform. The anterior stroma bulges more than the posterior stroma. In Figure 6B, the epithelium, Bowman's layer, and endothelium are emphasized by the color map. Strain values range from −3.5% to + 0.2%. Lower strain magnitudes are presented in magenta, and higher strain in light and dark blue. The epithelium compressionally deforms, and the endothelium and Bowman's layer bulge slightly, if at all.

FIG. 6.

FIG. 6.

Figure 7 shows average nominal axial strain through the thickness of the cornea. As with the strain images, this figures shows an expanding stroma with a compressing epithelium, but in this figure, the bipolar nature of the strain is even more apparent. The depth dependence of corneal strain measured with ultrasound most closely follows case 3 (Fig. 4). Strain processing averages over depth, and so anterior portions of epithelial strain average with zero strain from the coupling medium. This blending can be seen in Figures 5B and 6, but it does not have physical significance. To emphasize only physically significant data, depth in Figure 7 begins slightly inside the epithelium. Variations in strain in Figures 5, 6, and 7 include physical variations in the specimen and measurement uncertainty. Similar variations appear in previous studies. 12–14 Previous measurements of elasticity in corneal tissue, including our mechanical measurements, provide only an average bulk measurement, not the high resolution of the ultrasonic elasticity microscope.

FIG. 7.

FIG. 7.

Back to Top | Article Outline

DISCUSSION

The mechanical behavior of corneal layers indicated by strain imaging is reasonable given our experience with strain in layered tissue mimicking phantoms 12,13 and what is known about corneal structure. 11,27,28 Previous studies have shown that Bowman's layer and the combination of Descemet's membrane and the endothelium are stiffer than stroma. Layered tissue phantoms are created by pouring gelatin into a square mold and letting it set. The mold constrains the phantom on all sides except the top. From the top, it is deformed with a slitted plate. Because it is incompressible and restrained everywhere except the slit, the tissue phantom bulges into the slit when deformed.

In our experiment, stroma behaves similar to a single layer phantom. Descemet's membrane, the endothelial layers, and fluid pressure in the eye globe limit movement of the stroma into the anterior chamber. Like a single-layer phantom, the only place left for the stroma to move is into the slit. In contrast to a phantom, Bowman's layer prevents the stroma from fully expanding into the slit. As it acts against stromal movement, Bowman's layer is compressed. If we did not already know the general stiffnesses of corneal structures from ex vivo data, we could deduce them from the strain map.

Strain measured with ultrasound is predicted well with the finite element model considering the simplicity of the model and the potentially complex nature of corneal elasticity. The model assumes layers of three distinct elastic moduli. In real biologic tissue, we might expect a transition from one area of stiffness to the next. Also, ultrasonic strain data hints at inhomogeneity and a depth distribution of stromal elasticity. Anisotropy in the corneal plane is known to exist 29 but is not accounted for herein. Changes in hydration may explain some differences between the model guided by mechanical measurements and the ultrasonic measurement of strain, although exposure to water for both measurement systems was similar. Hydration might affect the specific distribution of elasticity, but would not affect basic mechanisms enabling ultrasonic measurement of corneal deformation. To reduce effects from hydration changes during the experiment, measurement time was kept as minimal as allowed by equipment. In clinical settings, the corneal pump mechanisms would maintain eye hydration, and measurements would take a matter of seconds. The necessary speed for clinical measurement is fundamentally possible, but would require additional engineering and more investment in advanced equipment.

With a finite element model, we obtained strain from elasticity data, but this process can be inverted. An elasticity map can be reconstructed from the strain image. Elasticity reconstruction algorithms have been developed 25,30,31 and applied to tissue mimicking phantoms. 13 Reconstruction is more complicated with a multilayered structure such as the cornea and is a subject of current investigation. Elasticity reconstruction will produce material properties, independent of geometry and other experimental parameters, enabling specimen to specimen studies of corneal elasticity.

Equipment and measurement geometry were arranged to emphasize the normal axial strain component. The slit in the deformation plate is long in the out-of-plane dimension and narrow in the lateral dimension. This geometry greatly reduces out-of-plane displacement and its resulting strain. 12,25,30 Assuming a two-dimensional model, lateral strain should be equal to axial strain in magnitude, but lateral displacement is much smaller than its axial counterpart. We calculated lateral displacement from the same data set, but, because of its small magnitude, lateral displacement and strain had a poor signal-to-noise ratio. Shear strains are also present in this experiment, but they occur mostly underneath edges of the slit. Because of small lateral displacement, shear strain also would suffer from poor signal to noise. Calculating shear strain from axial displacement should be possible using tissue incompressibility. 16

Another subject of ongoing study is the integration of ultrasound elasticity measurements into femtosecond laser refractive surgical instrumentation. Although we have shown that strain can be measured in corneas, additional work is needed before clinical application. In one goal toward this end, eliminating the slit in the deformation plate would allow simultaneous application for laser surgery and deformation for strain imaging, but the plate material must meet a tight set of requirements. It must be acoustically and optically transparent, yet stiff enough to apply the required deformation. The FDA has already approved femtosecond lasers for cutting corneal flaps. As part of that procedure, 3,4 a glass plate deforms the cornea until it is flat. Pressure increase from this procedure is under 10 mm Hg, well under four times less than FDA-approved limits. This deformation is larger than what was applied in our experiment.

An important motivation for the development of strain imaging is elasticity mapping as input for models of the biomechanical response of the cornea. An analytic model has been developed, comparing favorably with clinical data from refractive surgery. 4,10 In this model, values of Young's modulus were chosen to fit clinical data. The model does not match unless Young's modulus varies with depth. With reconstruction of an elasticity image, ultrasound can determine depth dependent variation of mechanical strength in individual eyes.

With input from ultrasonic elasticity imaging, a biomechanical model could guide surgeons in determining how much tissue to remove to produce the desired corneal shape. These models can aid the refinement of existing nomograms (i.e., LASIK) and the development of new procedures (femtosecond laser correction and femtosecond laser keratotomy). In this way, ultrasound elasticity may improve the predictability of clinical outcomes. After surgery, ultrasonic elasticity can assess changes in biomechanical properties of tissue caused by refractive procedures and the wound healing response. With the noninvasive nature of ultrasound and the minimally invasive nature of femtosecond laser surgery, a second fine-tuning refractive surgery may be possible if needed. Such a system may be capable of refraction correction of greater accuracy than current LASIK procedures.

Back to Top | Article Outline

CONCLUSIONS

We have shown that ultrasound can produce high-resolution strain images of a porcine cornea. The ultimate goal of this research is in vivo measurement of elasticity in human corneas. Use of the ultrasound elasticity microscope on a whole intact eye globe shows the usefulness of this technique for in vivo measurements. Within the cornea, variation of strain with depth may indicate a nonuniform distribution of elasticity. The overall trend of strain values from ultrasound compare well with finite element modeling guided by direct mechanical measurements of Young's modulus. This independent measure gives confidence to ultrasonic strain data.

Back to Top | Article Outline

Acknowledgments:

The authors thank Ramon Erkamp and Xunchang Chen in the Department of Biomedical Engineering at the University of Michigan for help with the direct mechanical measurements and strain image processing, respectively. We also thank Dr. Andrei Skovoroda of the Institute for Mathematical Problems in Biology of the Russian Academy of Sciences for useful discussions.

Back to Top | Article Outline

REFERENCES

1. Foster FS, Pavlin CJ, Harasiewicz KA, et al. Advances in ultrasound biomicroscopy. Ultrasound Med Biol 2000; 26: 1–27.
2. Pavlin CJ, Foster FS. Ultrasound Biomicroscopy of the Eye. New York: Springer-Verlag, 1995.
3. Juhasz T, Frieder H, Kurtz RM, et al. Corneal refractive surgery with femtosecond lasers. IEEE J Sel Top Quant 1999; 5: 902–10.
4. Juhasz T, Djotyan G, Loesel FH, et al. Applications of femtosecond lasers in corneal surgery. Laser Phys 2000; 10: 495–500.
5. Maller BS. Market trends in refractive surgery. Int Ophthalmol Clin 2000; 40: 11–9.
6. Gimbel HV, van Westenbrugger JA, Penno EEA, et al. Simultaneous bilateral laser in situ keratomileusis: safety and efficacy. Ophthalmology 1999; 106: 1461–8.
7. Stulting RD, Carr JD, Thompson KP, et al. Complications of laser in situ keratomileusis for the correction of myopia. Ophthalmology 1999; 106: 13–20.
8. Eshbaugh CG, Polchack JR, Krueger RR. Optics of refractive surgery. Int Ophthalmol Clin 2000; 40: 1–10.
9. Probst LE. Complex cases with LASIK: advanced techniques and complication management. Thorofare, NJ: Slack, 1999.
10. Djotyan GP, Kurtz RM, Fernandez DC, et al. An analytical solvable model for biomechanical response of the cornea to refractive surgery. J Biomech Eng-TASME 2001; 123: 440–5.
11. Woo SLY, Kobayashi AS, Schlegel WA, et al. Nonlinear material properties of intact cornea and sclera. Exp Eye Res 1972; 14: 29–39.
12. Cohn NA, Emelianov SY, O'Donnell M. An elasticity microscope. Part II: experimental results. IEEE Trans UFFC 1997; 44: 1320–31.
13. Cohn NA, Kim BS, Erkamp RQ, et al. High-resolution elasticity imaging for tissue engineering. IEEE Trans UFFC 2000; 47: 956–66.
14. Cohn NA, Emelianov SY, Lubinski MA, et al. An elasticity microscope. Part I: methods. IEEE Trans UFFC 1997; 44: 1304–19.
15. O'Donnell M, Skovoroda AR, Shapo BM, et al. Internal displacement and strain imaging using ultrasonic speckle tracking. IEEE Trans UFFC 1994; 41: 314–25.
16. Lubinski MA, Emelianov SY, Raghavan KR, et al. Lateral displacement estimation using tissue incompressibility. IEEE Trans UFFC 1996; 43: 247–56.
17. Lubinski MA, Emelianov SY, O'Donnell M. Speckle tracking methods for ultrasonic elasticity imaging using short time correlation. IEEE Trans UFFC 1999; 46: 82–96.
18. Passmann C, Ermert H. A 100-MHz ultrasound imaging system for dermatologic and ophthalmologic diagnosis. IEEE Trans UFFC 1996; 43: 545–52.
19. Knapik DA, Starkoski B, Pavlin CJ, et al. A 100–200 MHz ultrasound biomicroscope. IEEE Trans UFFC 2000; 47: 1540–9.
20. Erkamp RQ, Wiggins P, Skovoroda AR, et al. Measuring the elastic modulus of small tissue samples. Ultrason Imaging 1998; 20: 17–28.
21. Parker KJ, Huang SR, Musulin RA, et al. Tissue response to mechanical vibrations for sonoelasticity imaging. Ultrasound Med Biol 1990; 16: 241–6.
22. Chen EJ, Novakofski J, Jenkins WK, et al. Young's modulus measurements of soft tissues with application to elasticity imaging. IEEE Trans UFFC 1996; 43: 191–4.
23. Kurtz RM, Spooner GJR, Sletten KR, et al. Femtosecond laser corneal refractive surgery. In: Rol PO, Joos KM, Manns F, et al., eds. Ophthalmic Technologies IX. Proc SPIE, 1999;3591:209–19.
24. Ratkay-Traub I, Juhasz T, Horvath C, et al. Ultra-short (femtosecond) laser surgery: initial use in LASIK flap creation. Ophthalmology Clinics of North America 2001; 14: 347–55.
25. Skovoroda AR, Emelianov SY, O'Donnell M. Tissue elasticity reconstruction based on ultrasonic displacement and strain images. IEEE Trans UFFC 1995; 42: 747–65.
26. Reinstein DZ, Silverman RH, Sutton HFS, et al. Very high-frequency ultrasound corneal analysis identifies anatomic correlates of optical complications of lamellar refractive surgery: anatomic diagnosis in lamellar surgery. Ophthalmology 1999; 106: 474–82.
27. Bryant MR, McDonnell PJ. Constitutive laws for biomechanical modeling of refractive surgery. J Biomech Eng 1996; 188: 473–81.
28. Patel S, Reinstein DZ, Silverman RH, et al. The shape of Bowman's layer in the human cornea. J Refract Surg 1998; 14: 636–40.
29. Kampmeier J, Radt B, Birngruber R, et al. Thermal and biomechanical parameters of porcine cornea. Cornea 2000; 19: 355–63.
30. Skovoroda AR, Emelianov SY, Lubinski MA, et al. Theoretical analysis and verification of ultrasound displacement and strain imaging. IEEE Trans UFFC 1994; 41: 302–13.
31. Skovoroda AR, Lubinski MA, Emelianov SY, et al. Reconstructive elasticity imaging for large deformations. IEEE Trans UFFC 1999; 46: 523–35.
Keywords:

Cornea; Corneal properties; Elasticity; Refractive surgery; Strain imaging; Ultrasonic microscopy; Ultrasound

© 2002 Lippincott Williams & Wilkins, Inc.