Introduction
Corneal power and astigmatism have traditionally been determined using instruments that measure only the anterior corneal curvature. The total corneal power and astigmatism are subsequently estimated by means of the keratometric index (usually 1.3375), which takes into account the negative power of the posterior corneal surface. More than 30% of eyes display > 1.0 dioptre (D) of keratometric astigmatism (KA), which can limit the patient’s postoperative uncorrected visual acuity if not treated at the time of cataract surgery 1, 2. Thus, selecting the best method to measure corneal astigmatism (CA) is very important for minimising postoperative residual refractive astigmatism (RA).
The introduction of Scheimpflug tomography made it possible to measure the posterior corneal curvature and directly obtain the total corneal power and astigmatism (without relying on fictitious assumptions related to the keratometric index) 3. In normal eyes, posterior corneal astigmatism (PCA) mean values have been reported to range from 0.30 to 0.54 D, but the magnitude can be higher than 1 D in 5.7% of eyes 4, 5. As KA already accounts for the PCA through its keratometric index, the mean difference between total corneal astigmatism (TCA) and KA is lower (0.20–0.25 D) than the mean magnitude of PCA 6, 7. In good agreement with these data, corneal astigmatism as determined by TCA measurements corresponded better with manifest refractive cylinder than corneal astigmatism based on KA in phakic eyes 8.
To determine Scheimpflug tomography measurement that most accurately predicts the RA, we designed a study on non-toric pseudophakic eyes, where the influence of lenticular astigmatism can be excluded. Under the assumption that the TCA must be equal to the RA in these eyes, we compared manifest RA and different corneal astigmatism measurements, with the goal of determining which value is most closely associated with the refractive cylinder and thus, may be considered as the best option to calculate the power of toric intraocular lenses (IOLs).
Materials and Methods
Patients and surgical procedures
This retrospective study included a consecutive series of pseudophakic patients enrolled between November 2015 and January 2021. Each patient signed an informed consent document prior to enrolment, and all study methods adhered to the tenets described in the Declaration of Helsinki. The study protocol was reviewed and approved by the Institutional Review Board of Yeouido St. Mary Hospital (Seoul, South Korea). Exclusion criteria were a history of any ocular disease, previous ocular surgery (other than cataract surgery), and occurrence of intraoperative and/or postoperative complications; corrected distance visual acuity had to be higher than 20/40. Eyes with any degree of corneal astigmatism were enrolled in the study.
Phacoemulsification was performed through a 2.2-mm temporal incision, and a non-toric intraocular lens (IOL, ZCB00, Johnson & Johnson Vision, USA) was implanted into the bag.
Postoperative measurements
Measurements of RA and CA were performed six months after the surgery. Manifest refraction was measured in a plus-cylinder format using an automated refractometer with a vertex distance setting of zero mm (Canon RK-5, Tokyo, Japan), which was subsequently refined by subjective refraction. CA was assessed by measuring each eye once with an automated keratometer (Canon RK-5, Tokyo, Japan) and Scheimpflug tomography (Pentacam HR, software version 1.17r91; Oculus, Wetzlar, Germany). This was used to analyse the cornea via a 25-picture scan; only scans whose quality specification was graded as “OK” by the instrument software were included.
In addition to keratometric astigmatism (KA) using an automated keratometer, the following corneal power measurements using Scheimpflug tomography were evaluated:
Keratometric astigmatism (KA): This value represents the difference between the power of the steep and flat meridians. It is equivalent to the simulated K-value used in traditional corneal topography and is calculated by entering the corneal curvature radius into a thin-lens formula for paraxial imagery, which considers the cornea as a single refractive sphere. The corneal radii were converted into dioptric power values using the keratometric index of refraction (1.3375).
Total corneal refractive power (TCRP): This value corresponds to TCA, and it is automatically measured by ray tracing and calculated using the values for anterior corneal radius, posterior corneal radius, and corneal thickness. Snell’s law and the specific refractive indices of air, cornea, and aqueous humour were used to calculate the corneal power. Measurements at a 2.0-mm zone/ring, 3.0-mm zone/ring, and 4.0-mm zone/ring centred on both the pupil and the apex were considered for the statistical analysis.
IOL decentration and tilt were evaluated using Scheimpflug tomography, and cases in which IOL decentration was > 0.4 mm or tilt was > 5° were excluded from the study. IOL decentration and tilt were measured according to methods previously described by de Castro et al. 9 and Kranitz et al. 10. IOL decentration was determined as the distance between the IOL centre and the pupillary axis. Positive and negative horizontal decentration indicate temporal and nasal decentration, respectively, while positive and negative vertical decentration indicate superior and inferior decentration, respectively. In this study, we calculated the mean absolute value of IOL decentration. Total decentration is the magnitude of the vector resulting from the horizontal and vertical decentration. A positive IOL tilt around the x-axis indicates that the superior edge of the IOL is located on the front, while a positive tilt around the y-axis indicates that the temporal edge of the IOL moves backward when compared with the nasal edge. We also calculated the mean absolute value of IOL tilt.
Astigmatism analysis
Vectorial astigmatism analysis was performed using the Næser polar value method 11, 12. Specific calculations were performed as previously reported by Bregnhøj et al. (2015).
CA and RA can be described as two-dimensional vectors. For this purpose, the net astigmatism (M @ α), where M (M ≥ 0) is the astigmatic magnitude in dioptres (D) and α is the steep astigmatic direction in degrees (°), was transformed into two polar values in units of dioptres in the following general format
12:
α-Φ))
where Φ is the reference meridian (°). In the present study, we used a variable reference meridian Φ = steeper CA α to describe RA as a function of CA (Bregnhøj et al. 2015). Therefore, the reference meridian Φ varied for each corneal measurement modality and for each corneal measurement.
KP(Φ) is positive for astigmatic meridians α along plane Φ and negative for meridians along Φ+90. KP(Φ+45) is positive for astigmatic meridians α rotated in a counter-clockwise direction and negative for a clockwise rotation relative to Φ.
The difference in astigmatism (DA), that is, the difference between refractive and corneal astigmatism, was separately calculated for both polar values:
ΔKP(Φ) is positive for CA underestimation of RA and negative for overestimation. ΔKP(Φ+45) is positive for RA meridian rotated in a counter-clockwise direction and negative for a clockwise rotation relative to Φ.
The DA was reconverted to the standard net cylinder format using the following equations:
In this context, Equation 6 provides the meridian relative to the steep anterior corneal meridian. The average net astigmatism was calculated using the average values for KP(Φ) and KP (Φ+45). The average DA magnitude (not considering the meridian) was calculated by inserting the individual values in Equation 5, followed by an averaging procedure.
A low mean value of DA indicates a high accuracy, with the corneal astigmatism measurement well aligned in both meridian and magnitude, and the corneal plane manifesting as a refractive cylinder. Aggregate analysis of astigmatism was also performed using the astigmatism double angle polar tool 13.
All eyes were divided into three groups according to the orientation of the steep meridian (as measured with an automated keratometer), and the data for each group were separately analysed. In the with-the-rule (WTR) astigmatism group, the orientation of the steep meridian ranged from 61° to 120°, whereas in the against-the-rule (ATR) astigmatism group, it ranged from 151° to 180° or from zero to 30°. The cases with astigmatism meridian not belonging to the above two groups were classified as the oblique astigmatism group. (31° to 60° or 121° to 150°).
We calculated polar values along the variable meridians to investigate each corneal measurement modality for all astigmatic directions. This approach allowed us to differentiate between the most accurate (with the smallest DA values) corneal measurements.
Finally, we performed optimization of the corneal measurements based on multiple linear regression analysis. As previously reported, this optimization followed the general format
14:
Statistical analysis
All statistical analyses were performed using the IBM SPSS Statistics for Windows (version 21.0. Armonk, NY: IBM Corp). The normality of the data distribution was assessed using the Kolmogorov-Smirnov test. As the data were not universally normally distributed, non-parametric methods were consistently employed. One- and two-sample Wilcoxon signed-rank tests were performed to determine differences between refractive cylinder and corneal astigmatism and to evaluate whether the mean values of KP(Φ) and KP(Φ+45) were significantly different from zero. The Friedman test was also performed to determine differences in cylinder magnitude and polar values among all astigmatism measurements. Multiple linear regression tests were performed to build equations predicting KP(Φ)RA. Statistical significance was set at p < 0.05.
Discussion
In the present study, we found that several TCA values were in good agreement with RA. Moreover, we confirmed that corneal astigmatism values based on the curvature of both corneal surfaces (i.e. TCA) resembled RA more closely than those based on the anterior corneal curvature only (i.e. KA).
Ideally, corneal astigmatism measurement should provide us with consistently good results in eyes with WTR, ATR, and oblique astigmatism. In this regard, the value that gave the best result was the TCRP at 4.0 mm (zone/pupil). However, the residual astigmatism is small but not negligible, and should be considered when planning toric IOL implantation. A possible explanation may be related to imperfect measurements of corneal astigmatism by Scheimpflug tomography or to other unknown sources of error.
We found that the arithmetic mean DA determined by automated KA was lower than that determined by KA based on Scheimpflug tomography (
Table 3–6). This result is similar to that of a previous study
15. The relatively lower accuracy of Scheimpflug tomography may be because of the prolonged time required for image acquisition
16, due to the different diameters of the analysed zone or the different technologies used to measure the corneal curvature.
In the analysis on the whole cohort using the variable steep anterior corneal meridian as reference, the Scheimpflug TCRP measurement at 4.0 mm pupil/zone was the most accurate. Therefore, the present study suggests that this modality should be used in all Pentacam measurements to predict refractive astigmatism. The accuracy of this measurement improved only moderately as a result of optimization (
Table 9), further strengthening the case for a physiological representation of refractive astigmatism. In contrast, the accuracy of automated KA improved significantly during optimization, thereby reaching a similar accuracy as the optimized TCRP at 4.0 mm pupil/zone. For this group of patients and with these specific measuring devices, the accuracies of optimized automated KA and TCRP at 4.0 mm (pupil/zone) are therefore identical. This may not be true for other keratometers.
Our findings relative to TCA and KA are in good agreement with previous studies showing that TCA measurements mirror RA better than KA. For example, the accuracy of TCA (3.0 mm), as measured with both a Scheimpflug rotating camera and a colour light-emitting diode corneal topographer, was higher than that of KA in pseudophakic eyes with non-toric IOLs 17. Similarly, the prediction error in RA was lower with TCA (3.0 mm) than with KA in a sample of eyes implanted with toric IOLs 18. It will be interesting to see if TCA at 4.0 mm (which was not analysed in previous studies) will further improve the outcomes of toric IOL implantation.
In the present study, the measured KA values were optimized using a multiple linear regression formula, thereby predicting
KP(Φ)RA using
KP(Φ)KA and cos(2α), where α is the meridian of the steep corneal meridian (
Table 7).
(when corneal refractive index = 1.3315)
This model assumes that the directions of RA and CA are identical, as evidenced by the average values for ΔKP(Φ+45)
RA close to zero (as shown in
Table 3, the average ΔKP(Φ+45)
RA amounted to only -0.04 D). This correlation allows for an individual estimation of KP(Φ)
RA for any combination of KA magnitude and meridian, that is, for both eyes with WTR, oblique, and ATR astigmatism. For the WTR KA of 1.0 D along 90°, KP(Φ)
RA is predicted as 0.49 D and, therefore, the net RA net astigmatism as 0.49 D along 90°. The RA net astigmatism is predicted as 1.27 D along zero degrees for the ATR net cylinder 1.0 D along zero degrees. This correlation was developed in a previous study, based on a pseudophakic eye model of 184 eyes with a postoperative KA of 0.91 D, using a keratometric index of 1.3375
14 as
KP(Φ)RA = -0.09 + 0.68* KP(Φ)KA+ 0.33*cos(2α). In that study, the regression constant for
KP(Φ)KA was considerably lower than the similar value in this study, 0.68 versus 0.92, which might be explained by the lower KA values in the present study (
Table 2). In the present study, the regression constants for cos(2α) were quite similar, 0.33 versus 0.39, in the present study. The same methodology was used to optimize KA for toric IOL calculation
18, yielding the regression equation 0.103 + 0.836*KA + 0.457*cos2α. This equation is quite similar to the present correlation for KA, although the toric IOL optimization is far more complex with additional variables, such as exact vergence calculation and toric IOL effective power. The measurement precision is of the most importance in all optimization procedures, and might be improved by using the average of multiple, rather than single, determinations of refractive and corneal astigmatism. The necessity for regression constants of the optimized Scheimpflug tomography measurements in
Table 7 is a sign of imperfections of the measured values. These may be due to internal calculation algorithms, erroneous measurement methods, or unknown factors.
The present study has several strengths and limitations. Among the former, we ruled out confounding factors induced by lenticular astigmatism as well as by tilt and decentration; therefore, only pseudophakic eyes with a decentration < 0.4 mm and a tilt < 5 degrees were enrolled (the critical amounts of decentration and tilt have been reported to be 0.4~0.8 mm and 5~10 degrees, respectively) 19, 20. Furthermore, we used Scheimpflug tomography, whose reproducibility when taking measurements of IOL decentration and tilt has been shown to be high 9, 21. Regarding limitations, we did not measure TCA with other devices or technologies, which may offer interesting comparisons.
In conclusion, TCA measurements using Scheimpflug tomography were closer to the mean RA with respect to corneal astigmatism measurements derived from the anterior corneal surface only (i.e. KA). Among the several available options, the TCRP value centred on the pupil and including a 4.0-mm diameter zone was the one that most closely mirrored RA in pseudophakic eyes. Automated KA measurements provided similar accuracies after optimization. Thus, both values may be used for accurate toric IOL calculation.