|
|
||||||||
1 From the Ophthalmology Service, Bichat Hospital, Rothschild Foundation, Paris VII University, Paris, France; the 2 Institute of Statistics, Pierre et Marie Curie University, Paris, France; and the 3 Massachusetts Eye and Ear Infirmary and Schepens Eye Research Institute, Harvard Medical School, Boston, Massachusetts.
| Abstract |
|---|
|
|
|---|
METHODS. A mathematical model of aspheric myopic corneal laser surgery was generated. The initial corneal surface was modeled as a conic section of apical radius R1 and asphericity Q1. The final corneal surface was modeled as a conic section of apical R2 and asphericity Q2, where R2 was calculated from the paraxial optical formula for a given treatment magnitude (D), and Q2 was the intended final asphericity. The aspheric profile of ablation was defined as the difference between the initial and final corneal profiles for a given optical zone diameter (S), and the maximal depth of ablation was calculated from these equations. Using the Taylor series expansion, an equation was derived that allowed the approximation of the central depth of ablation (t0) for various magnitudes of treatment, optical zone diameters, and asphericity. In addition to the Munnerlyn term (M), incorporating Munnerlyns approximation (-D · S2/3), the equation included an asphericity term (A) and a change of asphericity term (
). This formula (t0 = M + A +
) was used to predict the maximal depth of ablation and the limits of customized asphericity treatments in several theoretical situations.
RESULTS. When the initial and final asphericities were identical (no intended change in asphericity; Q1 = Q2;
= 0), the maximal depth of ablation (t0 = M + A) increased linearly with the asphericity Q1. To achieve a more prolate final asphericity (Q2 < Q1; dQ < 0;
> 0), the maximal depth of ablation (M + A +
) was increased. For treatments in which Q2 was intended to be more oblate than Q1 (Q2 > Q1; dQ > 0;
< 0), the maximal depth of ablation was reduced. These effects sharply increased with increasing diameters of the optical zone(s). Similarly, in the case of PRK, the differential increase in epithelial thickness in the center of the cornea compared with the periphery resulted in increased oblateness.
CONCLUSIONS. Aspheric profiles of ablation result in varying central depths of ablation. Oblateness of the initial corneal surface, intentional increase in negative asphericity, and enlargement of the optical zone diameter result in deeper central ablations. This may be of clinical importance in planning aspheric profiles of ablation in LASIK procedures to correct spherical aberration without compromising the mechanical integrity of the cornea.
| Introduction |
|---|
|
|
|---|
Enlargement of the ablation diameter may be helpful in reducing the optical aberrations after excimer photoablation, but this results in increased depth of ablation. Deep ablations could weaken the mechanical integrity of the cornea,9 as suggested by several reports of keratectasia after LASIK in cases with high myopia corrections.10 11 12 Some investigators have proposed aspheric patterns of ablation to minimize spherical aberrations.13 14 15 16 However, the influence of the aspheric ablations on the depth of ablation and on mechanical stability of the cornea are not known. In this study, we used a mathematical analysis to predict the theoretical maximal depth of ablation for customized aspheric ablations that would allow correction of myopia as well as adjustment of the final corneal asphericity to desired values. We investigated the influence of the initial corneal apical radius of curvature, initial asphericity, intended diopteric correction, diameter of treatment, and intended change in corneal asphericity on the maximal depth of ablation.
| Materials and Methods |
|---|
|
|
|---|
![]() | (1) |
The postoperative corneal profile was modeled as a conic section of apical radius R2 and shape factor Q2:
![]() | (2) |
![]() | (3) |
The intended change in asphericity (dQ) was calculated as the difference between the preoperative and the postoperative asphericities: dQ = (Q2 - Q1). When dQ < 0 (Q1 > Q2), the final corneal surface is more prolate (or less oblate). When dQ > 0 (Q1 < Q2), the final surface is more oblate (or less prolate). When dQ = 0, the initial and final corneal surfaces have the same asphericity.
The pattern of ablation within an optical zone of diameter (S) was calculated as the difference in sagittal height between corresponding points of the initial and final surfaces, intersecting at the edge of the optical zone (y = S/2), which corresponds to the material removed between two aspheric surfaces whose curvature difference results in the targeted change in apical power and asphericity (Fig. 1) . The depth of ablation is zero at y = S/2. It increases as y approaches 0 (Fig. 1) .
|
![]() | (4) |
The Taylor Series Expansion for Depth Approximation
By using the Taylor series expansion up to the second order, t0 could be approximated by:
![]() | (5) |
![]() | (6) |
![]() | (7) |
![]() | (8) |
![]() | (9) |
and the second-order paraxial binomial expansion (equation 7)
. The second term is the initial asphericity term (A) and is a function of the initial shape factor (Q1) and of the Munnerlyn approximation (equation 8)
. The third term is the asphericity change term (
) which is a function of the intended change in corneal asphericity (dQ), the diameter of the optical zone to the fourth power (S4), and the initial apical corneal radius of curvature to the third power (R13; equation 9
).
Comparison of Munnerlyn, Analytical, and Approximation Methods
The theoretical maximal depths of ablations were compared for R1 of 7.8 mm, Q1 of -0.2, and S of 6 mm, using our approximation method (equation 5)
, and the analytical method (equation 4)
. Comparisons were performed for -2, -6, and -10 D corrections, while varying Q2 between -0.6 and +0.2. The contributions of the Munnerlyn term (M; equation 7
), the asphericity term (A; equation 8
), and the asphericity change term (
; equation 9
) were also calculated for R1 of 7.8 mm, Q1 of 0.2, and Q2 of -0.2 (dQ of -0.4). They were compared for magnitudes of treatments of -3, -6, -9, -12, and -15 D.
Ablation Depth Comparisons for Various Values of Q1 and dQ
We calculated the depths of ablation for treatment magnitudes of -1 to -15 D as a function of the initial asphericity Q1. For dQ = 0 (Q1 = Q2), we varied Q1 between -0.7 and +0.5 and compared the depths of ablation. For situations in which more prolate asphericity was intended (dQ < 0), we calculated the depths of ablation for R1 = 7.8 mm, Q1 = -0.2, and S = 6 mm. The magnitude of treatment varied between -1 and -15 D, and the depths of ablation were compared for dQ = 0, -0.2, -0.4, and -0.6. We also calculated and tabulated the incremental increase in ablation depth resulting from intentional increase in prolateness of the cornea (dQ < 0; equation 9
) for treatment diameters of 4 to 8 mm and for R1 of 7.5, 7.8, and 8.1 mm.
Influence of Epithelial Hyperplasia after PRK on Corneal Asphericity
In a previous study,18
we showed that after conventional excimer laser treatment for myopia conforming to the Munnerlyn paraxial formula,1
the postoperative theoretical corneal asphericity can be accurately approximated by a best-fit conic section. In initially prolate corneas, we noted a discrepancy between the clinically reported oblateness after excimer laser surgery for myopia and the theoretical prediction of increased prolateness. The discrepancy may be related to laser nomogram departures from the Munnerlyn formula, low accuracy of videotopographic measurements, and wound healing (epithelial hyperplasia, stromal remodeling).
Epithelial remodeling may modify the specific effect induced by the myopia ablation and could account for the observed clinical trend to oblateness. We used equation 9 to calculate the theoretical change in corneal asphericity induced by the difference of epithelial thickness in the center compared to the periphery of the treatment zone after a -5-D PRK myopia treatment conforming to the Munnerlyn equation. We varied the increment of central epithelial thickness between 0 and 30 µm for initial asphericity of -0.2 and R1 of 7.8 mm. The induced increase in asphericity (oblateness) was determined for treatment diameters of 5 to 8 mm.
| Results |
|---|
|
|
|---|
and the second-order paraxial binomial expansion, were compared to theoretical aspheric treatments, by using the finite analysis method (equation 4)
and our method of approximation (equation 5)
. The theoretical maximal depths of tissue ablation for different magnitudes of spherical myopic, treatments as calculated using equation 5
, are very similar to the maximal depths of ablation calculated using equation 4
. Figure 2
shows the theoretical maximal depth of laser ablation using our approximation method (equation 5)
compared with the finite analysis method (equation 4)
in corneas with initial radius of curvature of 7.8 mm and preoperative asphericity (Q1) of -0.2. Both methods of calculation show that the maximal depth of ablation increases if an increase in negative asphericity is intended.
|
|
), are shown for intended corrections of -3, -6, -9, -12, and -15 D and for initial R1, Q1, and Q2 of 7.8 mm, +0.2, and -0.2, respectively. The Munnerlyn approximation underestimates the maximal depth of ablation compared with M of equation 6
. The intended change of asphericity (dQ of -0.4; more prolate) requires additional tissue ablation of 8.5 µm.
|
Figure 3B
shows the influence of initial asphericity on maximal depth of ablation for intended similar corrections. Figure 3B
shows that the depths of ablation needed to maintain preoperative asphericity (dQ = 0) are greatest for oblate corneas (Q1 = +0.2) and lowest for prolate corneas (Q1 = -0.2). The depths of ablation are increased when an intentional change in the prolate direction is intended (dQ = -0.4). This increase in depth (
) is determined by R1 and S, but is unrelated to Q1 or to the magnitude of diopteric correction (D).
Ablation Depth Comparisons for Various Values of Q1and dQ
Equation 5
allows estimation of the maximal depth of ablation for a myopic spherical treatment: in this case, Q1 = 0, and dQ = 0. Figure 4
shows the linear variations of the theoretical maximal depth of ablation as a function of the initial asphericity in oblate (Q1 > 0) and prolate (Q1 < 0) corneas.
|
(equation 6)
is positive, and the maximal depth of ablation is increased. Conversely, when the final asphericity is more oblate (Q2 > Q1), the maximal depth of ablation is decreased. In both cases,
increases with the optical zone diameter to the fourth power and decreases with the apical radius of curvature to the third power. The value of this additional depth as predicted by equation 6
is independent of any variations of the magnitude of correction. Table 2
shows the values of this additional depth for an intended change of asphericity (dQ) of 0.1 and for different values for various treatment diameters (S) and for initial radii of curvature (R1) representing normal (R1 = 7.8 mm), steep (R1 = 7.5 mm), and flat (R1 = 8.1 mm) human corneas.
|
|
|
| Discussion |
|---|
|
|
|---|
Our mathematical model predicts that achieving an increase in corneal prolateness after excimer laser surgery requires greater depth of central photoablation, which is independent of the initial asphericity. Furthermore, in corneas that are initially prolate (Q1 < 0) the depth of ablation necessary to maintain initial asphericity (dQ = 0) is less than that required to preserve asphericity of initially oblate or spherical corneas. Accordingly, for patients with initially oblate corneas (Q1 > 0) in whom an aspheric ablation profile is intended to generate a prolate postoperative corneal shape (Q2 < 0; dQ < 0), the maximal depth of tissue ablation increases substantially, given the original oblateness (positive asphericity, A) and the intentional reduction in asphericity (positive
; equation 8
). This concept is illustrated in Figure 3B
.
Our theoretical analysis shows that further depth limitations may arise from attempting to increase the treatment diameter S. This effect can be predicted from the Munnerlyn equation,1
but our analysis shows that this effect is exaggerated if an increase in negative asphericity is attempted in initially oblate corneas. This can be deduced from equations 6
, 8, and 9, which indicate that asphericity (A) and asphericity change (
) are both proportional to the fourth power of the treatment diameter (S).
In conventional noncustomized excimer laser surgery for myopia, the goal is to correct the refractive error using arc-based mathematical calculus. Paraxial spherical models correspond to a particular case of our model, in which initial and final corneal asphericities are assumed to be identical and equal to 1. Munnerlyn et al.1 derived from their paraxial model a simplified approximation of the maximal depth of ablation for myopic spherical corrections, (depth of ablation = diopters of correction x ablation diameter2/3), which is incorporated into equation 7 . The Munnerlyn approximation was achieved by using binomial expansion. However, equation 5 shows that the predicted theoretical depth calculated from the Munnerlyn approximation underestimates the actual theoretical depth, because the binomial expansion was taken up only to the first order. In addition to the Munnerlyn approximation, equation 7 incorporates a second term that allows better estimation of the maximal theoretical depth of ablation induced by paraxial profiles of myopia ablation that do not take asphericity into consideration (Q1 = Q2 = 0; dQ = 0; Table 1 ). The value of this term is proportional to the magnitude of treatment and to the fourth power of the treatment diameter, thus assuming greater clinical relevance in patients with large pupils and for magnitudes of treatment greater than 7 D (Fig. 3) .
The normal human cornea is not spherical. Despite its shortcomings, modeling the corneal shape in cross section as a conic section is a better approximation and has been widely used21 22 23 24 since its introduction by Mandell and St. Helen in 1971.25 Most normal human corneas conform to a prolate ellipse and flatten from the center to the periphery (negative asphericity; Q1 < 0), but some corneas are oblate and steepen from the center to the periphery (positive asphericity; Q1 > 0). Figure 4 shows that the maximal theoretical depth of ablation when the surgeon seeks to maintain the initial corneal asphericity (Q1 = Q2) is slightly reduced for prolate corneas (Q1 < 0), compared with spherical and oblate corneas, when all other parameters are identical.
Determining the ideal postoperative asphericity for a given eye and a given myopia correction is beyond the scope of this article. Using a model eye featuring aspheric ocular interfaces and a gradient refractive index within the lens, Patel et al.19 have predicted that optimal optical imagery is produced when the corneal profile is represented by a flattening ellipse whose asphericity is between -0.35 and -0.15. Two recent studies using mathematical modeling and ray-tracing techniques to determine the ideal low spherical aberration ablation profile for the correction of myopia found it to be deeper and steeper, suggesting a lower intended postoperative asphericity.14 15 Conversely, using an optical design software to build a two conic surface model of the cornea, Munger16 determined that the optimal postoperative corneal asphericity that would maintain the preoperative aberrations increased nonlinearly (i.e., became more oblate) as a function of the magnitude of refractive correction. Further studies involving the use of ray-tracing techniques or the collection of wavefront sensing data may help in determining the best postoperative corneal profile in a given patient. However, it seems reasonable to postulate that customized ablations should retain the physiologic prolate corneal shape. In a recent theoretical study, we demonstrated that after conventional myopic excimer laser treatment conforming to the Munnerlyn paraxial formula, the postoperative theoretical corneal asphericity could be accurately approximated by a best-fit conic section. We also found that for initially oblate corneas (Q1 > 0), oblateness increased (Q2 > Q1 > 0), whereas for prolate corneas (Q1 < 0), prolateness increased (Q2 < Q1 < 0) within the treated zone after myopia treatment.18 The present study is in agreement with these results: the theoretical maximal depth of ablation induced by a paraxial treatment (spherical assumption) is deeper than needed for a prolate cornea to maintain its prolateness.
In practice, however, the cornea becomes oblate after conventional refractive excimer laser treatment for myopia.5 7 Holladay et al.5 have recently suggested that the loss of negative asphericity may be the predominant factor in the functional decrease in vision. Our clinical experience confirms the results of this study, showing a significant association between increased postoperative asphericity and greater myopia correction. Because the patterns of ablation of the existing laser devices are proprietary, we do not have access to them, and thus we cannot study separately the respective specific roles of the patterns of ablation and the biological healing, so as to explain the clinical observation of increased postoperative oblateness. The latter may be due to variations of the applied fluence on the corneal surface, to the incorporation of laser pretreatment protocols intended to reduce the incidence of postoperative central islands, or to stromal and epithelial remodeling after surgery. Another explanation is that the laser may become less efficient as we move peripherally, and the depth centrally would not be changed but less tissue than planned peripherally would be removed.
Epithelial hyperplasia after PRK may be a predominant factor in explaining the discrepancy between the clinical findings and the theoretical predictions. Topographical patterns have been shown to change with time,26 and variations of the epithelial thickness have been associated with refractive regression occurring after LASIK and PRK.27 28 29 30 31 Figure 6 illustrates that in addition to its effect on the apical power, an increase in central corneal thickness during wound healing could induce a modification in the corneal asphericity. The extent of epithelial and stromal thickening during wound healing after PRK are greater than those after LASIK.27 28 29 32 The in vivo clinical observations that epithelial hyperplasia is more common in eyes treated with small ablation zone diameters or with high magnitudes of treatment29 are consistent with the predictions of our model. To our knowledge, no clinical study has either compared the modification in asphericity after LASIK and PRK or investigated the possible correlation between the variation in corneal asphericity, apical power, and central corneal thickness.
Two studies have used corneal topography (Holladay Diagnostic Summary; EyeSys Laboratories, Houston, TX) to determine the corneal asphericity after excimer laser refractive surgery. In the study by Hersh et al.,7 mean asphericity for all patients 1 year after myopic PRK was Q2 = +1.05 (p2 = +2.05); preoperative asphericity was not reported. The mean preoperative corneal asphericity (Q1) measured under similar conditions, was reported to be -0.16 by Holladay et al.5 All corneas changed from a prolate to an oblate shape (mean Q2 - of +0.47), 6 months after LASIK for myopia. The shift toward oblateness was greater after PRK than after LASIK.
One limitation of our approach is the contribution of the crystalline lens to the reduction of optical aberrations, especially in that age-related lens changes may affect the determination of the ideal asphericity.33 34 35 The cornea would have to be progressively more prolate with age to compensate. Taking these clinical observations into consideration, certain features of our mathematical model may have to be modified to compensate for the postoperative trend toward increased oblateness. One possibility is to increase the reduction of postoperative asphericity by an amount similar to that reported in previous clinical studies.5 7 Based on Table 2 , an aspheric profile of ablation designed to preempt an oblate shift of +1.0 after LASIK would require an additional ablation depth of approximately 20 µm (optical zone diameter = 6 mm) compared with a Munnerlyn-based noncustomized ablation. Although this approach may improve the predictability of postoperative asphericity, it may not be sufficient, because the additional ablation, may exacerbate the biological healing and induce more regression after PRK, or may compromise corneal stability after LASIK, especially for large optical zone diameters and for high myopia corrections.
Another limitation of our theoretical analysis is that it is based on a static-shape subtraction model in which the postoperative corneal shape is determined only by the difference between the preoperative shape and the ablation profile. The biological effects of healing and the variations of the applied fluence at the cornea are not considered. Furthermore, our model neglects the influence of the transition zone. The increased curvature at the edges of the treated zone may introduce substantial optical aberrations under conditions of dim illumination. This increases the demand for larger treatment diameters, which would increase dramatically the depth of ablation (equations 6-9) .
In summary, our model provides a basis for predicting the variation in theoretical maximal depth of ablation induced by aspheric custom ablations to correct myopic refraction errors. Increasing negative asphericity without increasing the risk of ectasia for high magnitudes of treatment may be achieved by reducing the treatment diameter. The reduction of the optical zone diameter, however, may induce undesirable optical edge effects and may counterbalance the positive effect of restoring the prolate shape of the central cornea. Future studies of the relationships between optimal asphericity, based on the classic Q value, and wavefront aberration and further experimental work and clinical trials are necessary to compliment our theoretical calculations to refine the profiles of ablation and allow adequate control of postoperative corneal asphericity.
| Footnotes |
|---|
Submitted for publication July 30, 2001; revised October 26, 2001; accepted November 6, 2001.
Commercial relationships policy: N.
The publication costs of this article were defrayed in part by page charge payment. This article must therefore be marked "advertisement" in accordance with 18 U.S.C.
1734 solely to indicate this fact.
Corresponding author: Dimitri T. Azar, Director, Cornea and Refractive Surgery Services, Massachusetts Eye and Ear Infirmary, 243 Charles Street, Boston, MA 02114; dazar{at}meei.harvard.edu
| References |
|---|
|
|
|---|
This article has been cited by other articles:
![]() |
C. C. K. Chan and B. S. Boxer Wachler A Comparison of CustomCornea Myopia Algorithms for Wavefront-Guided Laser In Situ Keratomileusis Arch Ophthalmol, August 1, 2008; 126(8): 1067 - 1070. [Abstract] [Full Text] [PDF] |
||||
![]() |
C. Villa, R. Gutierrez, J. R. Jimenez, and J. M. Gonzalez-Meijome Night vision disturbances after successful LASIK surgery Br. J. Ophthalmol., August 1, 2007; 91(8): 1031 - 1037. [Abstract] [Full Text] [PDF] |
||||
![]() |
D. Gatinel, J. Malet, T. Hoang-Xuan, and D. T. Azar Corneal Asphericity Change after Excimer Laser Hyperopic Surgery: Theoretical Effects on Corneal Profiles and Corresponding Zernike Expansions Invest. Ophthalmol. Vis. Sci., May 1, 2004; 45(5): 1349 - 1359. [Abstract] [Full Text] [PDF] |
||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |