Abstract

Purpose. To compare the mean corneal power (Km) and total astigmatism (Ka) estimated by three methods: simulated keratometry (simK), true net power (TNP), and total corneal refractive power (TCRP) before and after femtosecond laser small incision lenticule extraction (SMILE) surgery. Methods. A retrospective, cross-sectional study. SimK, TNP, and TCRP from a Scheimpflug analyzer were obtained from 144 patients before and 6 months after SMILE surgery. Km and Ka were recorded as the mean of individual paracentral rings of 1.0 to 8.0 mm (R1 to R8). The surgically induced changes in Km (delta-simK, delta-TNP, and delta-TCRP) and Ka (delta-simKa, delta-TNPa, and delta-TCRPa) were compared to the changes in spherical equivalent of the cycloplegic refraction (delta-SE) and astigmatism (delta-RA). Results. Preoperatively, astigmatism values were greatest with simKa from R1 to R5 and greatest with TCRPa from R6 to R8. Astigmatism values were smallest with TNPa from R1 to R7. Postoperatively, astigmatism values were greatest with simKa from R1 to R5 and greatest with TCRPa from R6 to R8. Delta-TCRP3 and Delta-TCRP4 matched delta-SE most closely, and delta-TCRPa3 matched delta-RA most closely. Conclusions. TCRP proved to be the most accurate method in estimating corneal power and astigmatism both before and after SMILE surgery.

1. Introduction

With advancements in technology, the ability to measure corneal power and other corneal parameters is becoming more and more precise. With an increase in the number of cases of refractive surgery occurring each year, the need for accurate and reliable measures of corneal power will continue to become more apparent. Although several measures currently exist, their lack of congruity and agreement has led to the need to evaluate which measures are available to consistently provide the most accurate estimates of corneal power and astigmatism. The measurements of total corneal power and astigmatism by corneal tomography are influenced by several factors: inclusion of one or two corneal surfaces, refractive index, refractive effect (e.g., the spherical aberration), and the location of the principal planes [1, 2].

Conventional keratometry (simulated keratometry) measures the paracentral anterior surface of the cornea (usually 3 mm in diameter) and calculates corneal power using a standard keratometric index of 1.3375. The standard keratometric index is based on the assumption that the corneal thickness is constant and the ratio between the curvatures of the anterior and posterior surfaces is constant. True net power (TNP) is the corneal power calculated using Gaussian optics. TNP measures both the anterior and posterior surfaces of the cornea and uses the real refractive index values of air, cornea, and aqueous. It is assumed to be more accurate than simulated keratometry (simK). However, the Gaussian optics formula is a simplified formula that uses a paraxial approximation and assumes that the rays propagating through the posterior surface of the cornea are parallel. Following laser ablation of the cornea, the above assumptions are no longer valid due to a change in both the corneal thickness and the ratio of the anterior/posterior radius of curvature, leading to an error in the measurement of corneal power and astigmatism [15]. Moreover, neither simK nor TNP takes the refractive effect into consideration and could be inaccurate in measuring peripheral corneal power.

Total corneal refractive power (TCRP) is calculated using Snell’s law of refraction by exact ray tracing. All of the four factors mentioned above are taken into account. It is claimed to be more accurate than TNP and simK, especially in eyes following myopic laser ablation. Several studies comparing the accuracy of simK, TNP, and TCRP in the measurement of corneal power have been published [611]. However, no study has been conducted in the evaluation of corneal astigmatism. Moreover, past studies have primarily evaluated the central corneal power and the profile of the peripheral cornea has not yet been thoroughly investigated.

The purpose of this study was to investigate the corneal power and astigmatism distribution measured by three methods: simK, TNP, and TCRP, both before and after femtosecond laser small incision lenticule extraction (SMILE) surgery. To evaluate the accuracy of the three methods in the measurement of surgically induced changes in corneal power and astigmatism, a comparison was made between them and the changes in cycloplegic refraction following refractive surgery.

2. Methods

2.1. Patients

144 patients who underwent surgery for the correction of myopia and myopic astigmatism by femtosecond laser small incision lenticule extraction surgery were recruited between March 2013 and December 2013 to participate in this retrospective study. The study was conducted at the Ophthalmology Department of the Eye and ENT Hospital, Shanghai, China. The inclusion criteria consisted of subjects who had spherical refractions between −1.00 and −10.00 diopters (D) and astigmatism between −0.25 and −5.00 D. Additionally, best-corrected distance visual acuity values of 20/25 or better were required along with a stable refraction over two years prior to surgery, an absence of other pathologic ocular conditions or relevant systemic diseases, and a minimum of 6 months’ follow-up. This study followed the tenets of the Declaration of Helsinki and was approved by the ethics committee of the EENT Hospital of Fudan University. Informed written consent was obtained from all subjects prior to participating in the study.

2.2. Surgical Technique

All surgeries were performed by the same surgeon (ZX). The femtosecond laser small incision lenticule extraction procedures were performed using a VisuMax femtosecond laser system (Carl Zeiss Meditec AG, Jena, Germany) following the surgical procedure outlined by Sekundo et al. [12]. The intended thickness of the upper tissue arcade (the cap) was 120 μm with an intended diameter of 7.5 mm. The diameter of the refractive lenticule was 6.2 mm to 6.5 mm. A single side cut of 90 degrees with a circumferential length of 2 mm was made in the superior position. No intraoperative or postoperative complications were observed.

2.3. Assessment of Refractive Changes and Corneal Power

The patients were followed 1 day, 1 week, 1 month, 3 months, and 6 months after surgery. The preoperative and 6 months’ postoperative measurement results were analyzed. Objective and subjective refraction assessments were performed before and after cycloplegia. The values were adjusted for the corneal plane using a vertex distance of 12 mm. Additionally, uncorrected and corrected distance visual acuities (UDVA and CDVA) were recorded.

Pentacam (Oculus GmbH, Wetzlar, Germany) imaging of the cornea was performed by the same experienced examiner, and three measurements were averaged for each result. Only the scans marked “OK” by the instrument were saved and analyzed. Corneal power was measured and shown in the power distribution display for the three methods: simK, TNP, and TCRP. simK was derived from the axial curvature map and was calculated with the keratometric index (). TNP is the corneal power calculated using the Gaussian formula as follows: where is the anterior corneal power, is the posterior corneal power, is the pachymetry result, and , which is the refractive index of the cornea. In the Pentacam program, TNP was calculated as the sum of and . TCRP was calculated as the corneal power using the ray-tracing method. Corneal thicknesses and curvatures of the anterior and posterior surfaces were obtained by Scheimpflug imaging. The calculation of TCRP is based on Snell’s law of refraction using real refractive index values (1 for air, 1.376 for the cornea, and 1.336 for aqueous). Corneal power is then determined by where and f is the focal length, referenced to the anterior corneal surface. Values (Ksteep, Kflat, and axis of Kflat) demonstrated by individual rings centered on the pupil center ranging from 1.0 mm to 8.0 mm (R1 to R8) and the mean values of the central 3 mm zone were recorded. The arithmetic mean of the steep and flat corneal axes and the difference between them (the astigmatism) were calculated and compared. The surgically induced changes in corneal power (delta-simK, delta-TNP, and delta-TCRP) were compared to the changes in cycloplegic refraction of spherical equivalent (delta-SE) which were calculated by subtracting the postoperative SE from the preoperative SE.

Vector analysis was performed to decompose the cylinder notation into two cross cylinder components, and [13, 14]. and are defined as where C is the magnitude of the cylinder power present and α is the axis in radians. represents the horizontal/vertical component of astigmatism and represents astigmatism with axes at 45° and 135°. Vector mean (poweraxis) is calculated as

The following two indices were used to compare corneal astigmatisms: the magnitude of the vector difference between two corneal astigmatisms and the absolute value of the angular difference between two corneal astigmatisms which was defined always to be an acute angle.

The surgically induced astigmatism (SIA) is the amount and axis of astigmatic change caused by surgery. It was determined as the vector difference between the pre- and postoperative refractive astigmatism (or corneal astigmatism) [13, 14]. The surgically induced changes in corneal astigmatism (delta-simKa, delta-TNPa, and delta-TCRPa) were compared to the changes in cycloplegic refraction of astigmatism using Bland-Altman analysis [15].

2.4. Statistical Analysis

Statistical analysis was performed using a commercially available statistical software package (SPSS, ver. 13.0; SPSS, Chicago, IL). Comparisons between the three methods were performed using analysis of variance (ANOVA) and least significant difference (LSD) analyses. Comparisons between refractive changes and corneal power changes were performed using a paired t-test. A value of less than 0.05 was considered statistically significant. The bias between delta-refractive astigmatism (delta-RA) and delta-corneal astigmatism (delta-Ka) was evaluated using 95% limits of agreement. 95% limits of agreement were the mean difference ± 1.96 × standard deviation of the difference [15].

3. Results

The right eyes of 144 patients (88 females (61.1%)) were included in the analyses. The mean age of patients was 26.3 ± 6.6 years, with a range of 18 to 41 years. The mean preoperative and postoperative cycloplegic spherical errors at the corneal plane were −5.73 ± 1.33 D (range: −1.71 to −8.12 D) and 0.12 ± 0.51 D (range: −1.95 to +1.53 D), respectively. The mean preoperative and postoperative cycloplegic astigmatisms at the corneal plane were −0.92 ± 0.61 D (range: −3.70 to 0.00 D) and −0.10 ± 0.21 D (range: −1.22 to 0.00 D), respectively. The mean preoperative CDVA (logarithm of minimal angle of resolution (logMAR)) was −0.02 ± 0.06 (range: −0.18 to 0.30). The mean postoperative UDVA and CDVA were −0.02 ± 0.08 (range: −0.18 to 0.40) and −0.05 ± 0.06 (range: −0.20 to 0.10), respectively. Delta-spherical equivalent (delta-SE) was 6.61 ± 1.37 D, and delta-refractive astigmatism (delta-RA) determined by vector analysis was (0.58 ± 0.61 D) × 90.8°, respectively.

3.1. Preoperatively

Significant differences were found between the three corneal powers from R1 to R8 (ANOVA: , Figure 1). TNP produced the smallest estimates among the three corneal powers from R1 to R8 (LSD: ). simKs were greater than TCRPs from R1 to R3, and TCRPs were greater than simKs from R6 to R8 (LSD: ). When calculating the mean of the central 3 mm zone, simK estimates (43.31 ± 1.46 D) were the greatest and TNP estimates (41.92 ± 1.40 D) were the smallest among the three corneal powers (LSD: ).

The vector means of three corneal astigmatisms and the comparisons between them are listed in Table 1. Astigmatism values were greatest with simKa from R1 to R5 and greatest with TCRPa from R6 to R8. Astigmatism values were smallest with TNPa from R1 to R7. The absolute values of angular difference between TCRPa and simKa were greatest with R1 and smallest with R4; the absolute values of angular difference between TCRPa and TNPa were greatest with R8 and smallest with R3 to R5; the absolute values of angular difference between simKa and TNPa were greatest with R1 and smallest with R4 and R5.

3.2. Postoperatively

Postoperatively, Kms of simK, TNP, and TCRP were all significantly decreased from the preoperative Kms for R1 to R8 (paired t-test, ). simKs were greater than TNPs for R1 to R8 (LSD: , Figure 2). simKs were greater than TCRPs for R1 to R6 and smaller than TCRP for R8 (LSD: ). TCRPs were greater than TNPs for R4 to R8 (LSD: ). The mean postoperative corneal power of the central 3 mm zone was greatest with simK (37.90 ± 1.96 D) and smallest with TNP (35.95 ± 2.00 D, LSD: ).

The vector means of postoperative corneal astigmatisms are listed in Table 2. When compared to the preoperative Ka, astigmatism values were decreased for R1 to R5 and increased for R6 to R8, as measured by all the three methods. Astigmatism values were greatest with simKa from R1 to R5 and greatest with TCRPa from R6 to R8. TNPa and TCRPa matched each other well from R1 to R5. The absolute values of angular difference in astigmatism between TCRPa and simKa were greatest with R2 and smallest with R6; the absolute values of angular difference in astigmatism between TCRPa and TNPa were greatest with R8 and smallest with R3 to R5; the absolute values of angular difference between simKa and TNPa were greatest with R2 and smallest with R6.

3.3. Surgically Induced Changes in Corneal Power (Delta-Km) and Astigmatism (Delta-Ka)

Delta-Kms were greatest with delta-TCRP and smallest with delta-simK from R1 to R8 (LSD: , Figure 3). Delta-TCRP for R3 (6.63 ± 1.20 D) and R4 (6.64 ± 1.16 D) matched delta-SE (6.61 ± 1.37 D) most closely. The SIAs calculated by vector analyses, the means of differences between the vector values of delta-RAs and delta-Kas, the standard deviations, and the 95% limits of agreement are listed in Table 3. Delta-RA (0.58 × 90.8°) was smaller than delta-Ka of R1 and R2 measured by all the three methods. Delta-TCRPa of R3 (0.57 × 94.8°) matched delta-RA most closely. Delta-TNPa of R3 (0.52 × 94.4°) and delta-simKa of R3 (0.46 × 94.7) underestimated delta-RA.

4. Discussion

The need for an accurate and reliable measure of corneal power and astigmatism is becoming more apparent, especially in the constantly advancing setting of corneal refractive surgery. While several methods already exist to measure corneal power and astigmatism, many of them lack great consistency and accuracy between them. Some of the commonly used methods also make assumptions or use parameters in their estimations that are not valid following refractive surgery. The aim of this study was to evaluate the accuracy of three commonly used measures of corneal power and astigmatism before and after myopic corneal refractive surgery. The accuracy of these three measures was determined through comparisons between their estimated surgically induced changes in corneal power and astigmatism with the changes in cycloplegic refraction obtained clinically following refractive surgery.

In the current study, total corneal power was calculated by three methods: simK, TNP, and TCRP. In the central cornea prior to corneal refractive surgery, simK estimated the greatest values of corneal power while TNP estimated the smallest of the three measures. This has been confirmed in several studies in which TCRP was found to be most accurate, while simK overestimated and TNP underestimated the central corneal power [24]. This could be a result of the keratometric index error involved in the simK calculation and the paraxial approximation used in the TNP calculation. Moreover, the different reference planes used by the different methods may also contribute to the error. The ray-tracing method uses the anterior corneal surface as the reference plane, whereas the Gaussian optics formula uses the second principal plane immediately in the front of the cornea.

When analyzing the peripheral cornea, we found that TCRPs increased from the central to peripheral cornea, while both simK and TNP decreased from the central to peripheral cornea. The distribution of TCRPs reflected the positive spherical aberration of the cornea [16]. Consequently, TCRP should be adopted when an accurate measurement of peripheral corneal power is needed, as it most closely parallels the changes in the higher-order aberration across the corneal profile from the center to the periphery. Although TCRP has been claimed to be the most accurate measure of corneal power, there is currently a very limited number of studies investigating its application. Therefore, this approach needs to be validated and optimized before it can be adopted for routine use in the measurement of corneal power. Moreover, since the current intraocular lens (IOL) formulas are based on keratometry values, TCRP might not be utilized in these formulas until a calculation based on the ray-tracing method has been established.

The results of the current study obtained following refractive surgery were consistent with those of previous studies [311], demonstrating that simK significantly overestimated postoperative central corneal power while TNP significantly underestimated postoperative central corneal power. Although TCRPs slightly overestimated [6, 9] or underestimated [8] delta-SE, they matched delta-SE more closely than delta-simK or delta-TNP. The errors induced by the assumptions used in the calculations of simK and TNP were exaggerated by corneal ablation, which changed both the corneal thickness and the ratio of the anterior/posterior radius of curvature. We also found that delta-TCRP of the 3.0 mm and 4.0 mm rings matched delta-SE most closely. Savini et al. [6] found that the delta-TCRP of the 2.0 mm ring and the delta-TCRP of the central 3.0 mm zone provided the closest approximation to the delta-SE after myopic refractive surgery. Savini et al. [7] also demonstrated that TCRP calculated over the entrance pupil ensured the most accurate estimation of delta-SE. Oh et al. [8] established that the delta-TCRP of the 4 mm zone exhibited the least difference with the noncycloplegic delta-SE following myopic refractive surgery. The different results obtained in the studies mentioned above could be related to the variable sizes of the optical zone, the variable ablation profiles, and the different methods of calculation (pupil or apex center, ring or zone) used in the individual investigations. Therefore, the conclusion of an individual investigation should not be generalized across all procedures and should be applied cautiously to other laser surgical platforms.

The accuracy of TCRP in the measurement of astigmatism in virginal cornea has been proved in several studies where TCRP has been used as the gold standard [1719]. We found that simK produced the largest and TNP produced the smallest estimates of the three corneal astigmatism measurements in the central cornea. We also found that TCRP and TNP matched each other much better than simK, demonstrating that the error involved in the power calculation of simK also exists in the astigmatic calculation. Roh et al. [20] compared astigmatism measurements using automated keratometry versus TCRP on cataract patients and found that TCRPs were more accurate in astigmatism measurements in highly irregular corneas. Calossi [16] found that the difference in astigmatism between simK and TCRP depends on the axis orientation. simK overestimates TCRP (mean 0.11 ± 0.22 D) in the eyes with with-the-rule astigmatism, and simK underestimates TCRP (mean 0.26 ± 0.31 D) in the eyes with against-the-rule astigmatism. Consequently, the eyes of with-the-rule and against-the-rule (as well as oblique) astigmatisms should be evaluated separately. Unfortunately, since the majority of the eyes in the current study had with-the-rule astigmatism and the sample sizes were very small for the other two types of astigmatism, any comparison between them is difficult. Further study with more eyes with against-the-rule and oblique astigmatisms will be done to clarify this question.

To our knowledge, this is the first study evaluating the profile of total corneal astigmatism as a function of corneal diameters. We found that astigmatism measured by simK increased from 1 mm to 4 mm and then decreased. No significant difference was found for TNP from 1 mm (1.18 D) to 8 mm (1.34 D). Astigmatism measured by TCRP increased from 1 mm (1.04 D) to 8 mm (1.50 D). An increase in corneal astigmatism from the center to the periphery has also been reported by Mathur et al. [21]. Consequently, TCRP should be used in the measurement of peripheral corneal astigmatism, for example, in the evaluation of contact lens fittings, since the peripheral corneal contour could influence the fitting success of contact lenses [22, 23].

As for the postoperative corneal astigmatism, we found that delta-Ka of R1 and R2 measured by all the three methods overestimated delta-RA. Delta-TCRPa of R3 matched delta-RA most closely while both delta-TNPa3 and delta-simKa3 underestimated delta-RA. Thus, TCRP proved to be accurate in the estimates of both corneal power and corneal astigmatism after refractive surgery. The current results may provide significant contributions to contact lens fitting following refractive surgery, as they can provide practitioners with a tool to evaluate changes in the central and peripheral corneal astigmatisms as a function of postoperative manifest refractive changes [24, 25].

We also found that all postoperative Kms (from R1 to R8) were decreased from the preoperative TCRPs. This is beyond the size of the lenticule, as well as the size of the cap. In contrast, the postoperative astigmatism values were decreased from the preoperative Ka from R1 to R5 and increased from R6 to R8. The creation of the cap and the postoperative changes in peripheral corneal contour after the removal of the lenticule may contribute to the peripheral changes in corneal power [26]. Additional studies applying other imaging systems, such as anterior segment optical coherence tomography (AS-OCT), could possibly contribute to solving this question. The changes in corneal power beyond the size of the lenticule can create a transitional zone-like area, which may help improve the visual quality.

A limit of this study was that the higher order aberrations (HOAs) were not included in the analyses. Several studies have reported significant increases in HOAs following SMILE surgery, especially increases in the coma. The change in the corneal asphericity, as well as the change in pupil diameter, may have an effect on the postoperative corneal astigmatism [2729]. Further investigations, including pre- and postoperative HOAs, would help to clarify this problem.

In conclusion, simK overestimated the central corneal power (and astigmatism) and TNP underestimated the central corneal power (and astigmatism) both before and after refractive surgery. TCRP proved to be the most accurate method among them, especially after refractive surgery. The distribution of TCRPs reflected the profile of positive spherical aberration across the cornea. Delta-TCRP of R3 and R4 matched the cycloplegic delta-SE most closely. Delta-TCRPa of R3 matched delta-RA most closely. Although the determination of TCRP has been claimed to be the most accurate measure of corneal power, this approach needs to be further validated and optimized before it can be adopted for routine use in the measurement of corneal power and astigmatism in a clinical setting.

Conflicts of Interest

The authors declare that there is no conflict of interest regarding the publication of this paper.

Authors’ Contributions

Yishan Qian and Yan Liu contributed equally to this article and thus are considered co-first authors.

Acknowledgments

This study was supported by the Natural Science Foundation of China (61108050).