
==== Front
J Ophthalmol
J Ophthalmol
JOPH
Journal of Ophthalmology
2090-004X
2090-0058
Hindawi

37483313
10.1155/2023/2261831
Research Article
A Novel Role for Corneal Pachymetry in Planning Cataract Surgery by Determining Changes in Spherical Equivalent Resulting from a Previous LASIK Treatment
https://orcid.org/0000-0003-3225-5035
Hammer Arthur arthur.hammer6@gmail.com
1
https://orcid.org/0000-0001-5297-2301
Heeren Tjebo F. C. 1
https://orcid.org/0000-0001-6635-5283
Angunawela Romesh 1 2
https://orcid.org/0000-0001-6536-0101
Marshall John 3
https://orcid.org/0009-0005-4129-3457
Saha Kamran 1
1Moorfields Eye Hospital NHS Foundation Trust, London, UK
2OCL Vision, London, UK
3University College London, Institute of Ophthalmology, London, UK
Academic Editor: Nilufer Yesilirmak

2023
31 5 2023
2023 22618318 3 2023
22 4 2023
27 4 2023
Copyright © 2023 Arthur Hammer et al.
2023
https://creativecommons.org/licenses/by/4.0/ This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.
Objectives

To provide a metric to differentiate between hyperopic and myopic ablation of a prior LASIK treatment based on the corneal pachymetry profile after laser vision correction (LVC).

Methods

Pachymetry data were retrospectively recovered from patients who had previous LASIK for refractive purposes between 2019 and 2020. Patients with any corneal disorder were excluded. Ablation spherical equivalent was predicted from the central to semiperipheral corneal thickness (CPT) ratio, both values were provided by using the Pentacam user interface software (UI), and values were computed from extracted raw pachymetry data.

Results

Data of 157 eyes of 81 patients were collected, of which data were analysed for 73 eyes of 73 patients to avoid concurrence of measurements in both eyes per subject (42% female; mean age 40.9; SD 12.8). The CPT ratio cutoff for distinction between myopic and hyperopic LASIK was 0.86 for Pentacam UI data. Sensitivity and specificity were 0.7 and 0.95, respectively. Accuracy increased with computation of the CPT ratio based on extracted raw data with sensitivity and specificity of 0.87 and 0.99, respectively. There was a marked linear correlation between the CPT ratio and the ablation spherical equivalent (R2 = 0.93).

Conclusions

CPT ratio cutoffs can correctly classify if a cornea previously had a hyperopic versus myopic LASIK surgery and estimate the ablation spherical equivalent of such treatment. This could prove useful for increased accuracy of intraocular lens (IOL) calculations for patients with no historical data of their prior LVC surgery at the time of cataract surgery planning.
==== Body
pmc1. Introduction

Intraocular lens (IOL) power calculations are a central element to today's cataract surgery and accuracy to predict postoperative refraction increases with technological advances. Today's most used formulas are vergence formulas using biometric measurements such as keratometry readings and axial length to estimate the effective lens position [1]. Newer formulas based on ray tracing [2] and artificial intelligence [3–7] are gaining in popularity but do not represent mainstream usage among cataract surgeons [8]. The refractive power of the cornea can be calculated using keratometry measurements, which are based on an estimated constant ratio between the anterior and posterior corneal curvature [9].

After laser vision correction (LVC), the relationship between the anterior and posterior corneal curvature is modified and the artificially assigned value of the keratometric index of the cornea changes accordingly [9]. IOL calculation formulas that do not account for this change result in higher refractive errors after cataract surgery, in particular, a hyperopic shift after myopic excimer laser treatments and to a lesser extent a myopic shift after hyperopic excimer laser treatments [10–12]. Other sources of refractive surprise can be related to anterior corneal curvature measurements outside of the treatment zone and wrong prediction of the effective lens position based on anterior corneal curvature measurements [9].

New formulas have been developed to address this issue. Some of those formulas, such as Haigis L and Barrett true K, can be used without LVC historical data such as the pre-LVC refraction. Average values and regression analysis based on LVC cohort's data can be used in such instances [9]. However, those formulas still require the surgeon to input the type of LVC correction received: hyperopic versus myopic.

With the ever-increasing popularity of LVC surgery, the number of patients seeking cataract surgery who have had prior treatments such as LASIK, PRK, or SMILE is rising. However, in practice, a significant proportion of those patients do not have access to their historical refractive data or to their surgical records.

For this reason, IOL calculation formulas not using refractive historical data, such as Haigis L and Barrett true K (no history) are of paramount interest, but it has been shown that providing patients' refractive historical data improves the accuracy of such formulas [13]. Therefore, in addition to being able to estimate historical data such as the type of LVC correction (myopic versus hyperopic), estimating also the ablation spherical equivalent, in dioptres, applied to a cornea could prove useful to increase the accuracy of those formulas.

It has been shown that the corneal thickness profile from the central to the peripheral cornea of myopic patients from low myopia (<6D) to high myopia (>12D) does not show significant differences [14]. Central corneal thickness (CCT) variation in the population has been extensively studied and yielded controversial results as to the presence of a correlation between the refractive error and CCT [15].

Laser ablation profiles for myopic and hyperopic correction theoretically change the thickness profile of the cornea in an opposite fashion, with myopic treatments thinning the central cornea and hyperopic treatments thinning the semiperipheral cornea [16]. As a consequence, the ratio between the central and semiperipheral corneal thickness should represent the amount of ablation spherical equivalent (SE) related to tissue ablation. Epithelial fluctuations will be discussed later.

This is a pilot study investigating if this central to semiperipheral corneal thickness ratio (CPT ratio) would be suitable as a measure to predict if a cornea previously received a hyperopic versus myopic LASIK treatment and if this CPT ratio could estimate the ablation spherical equivalent of this LASIK treatment.

2. Methods

Retrospective pachymetric data were extracted from a corneal tomographer (Pentacam, OCULUS Optikgeräte GmbH, Münchholzhäuser Straße 29, 35582 Wetzlar, Germany) for a pool of patients treated by LASIK for refractive purposes over a 1-year period from 2019 to 2020. No patients with pure astigmatism were included in this study. For most of the patients in this study, the ablation sphere had a higher weight on the ablation spherical equivalent than the ablation cylinder, except for a few hyperopic patients for whom the ablation spherical equivalent was 0 dioptre, in which case the weight of the sphere and cylinder on the ablation spherical equivalent was equal. Pentacam images were acquired by skilled technicians. All procedures have been performed in the same surgical centre (OCL Vision, London, UK) using the same excimer laser platform (Schwind Amaris 1050 RS, SCHWIND eye-tech-solutions GmbH, Mainparkstraße), with a 6 to 7 mm optical zone, and the same wavefront optimised profile. Only patients without preexisting ocular pathology other than refractive errors were included. Patients with corneal scarring, corneal ectasia, corneal dystrophy, or any disorder affecting corneal homeostasis were excluded. Only corneas with one excimer laser treatment were included. No patients were using contact lenses after LVC. Pentacam images acquired at least 6 months after LASIK were used to avoid wound healing bias effects on the measurements.

The main research question was as follows: can a cutoff value of the CPT ratio alone determine if a cornea had a prior hyperopic or myopic treatment, defined as the ablation spherical equivalent of in dioptres? Hyperopic treatment was defined as 0 or more dioptres (D) of ablation spherical equivalent. 0 dioptre of ablation spherical equivalent has been attributed to the hyperopic group because we used the negative cylinder astigmatism notation. Therefore, only hyperopic treatments can lead to an ablation spherical equivalent of 0 dioptre. Receiver-operating characteristics (ROC) curves, specificity, sensitivity, and accuracy were calculated to assess the predictability of the model.

The formula used to calculate the CPT ratio was(1) CPT ratio=central corneal thicknesssemiperipheral corneal thickness at a given radius.

All data handling and statistical analysis were performed using statistical programming language R [17].

Data were extracted in two ways for separate analysis. First, the CPT ratio was calculated from the corneal thickness data taken from the Pentacam software user interface (Pentacam UI) at a radius of 0 and 3 mm (6 mm diameter), as shown in Figure 1 (Pentacam UI CPT ratio). Second, raw pachymetry data were extracted from the Pentacam using the manufacturer's csv file export function (raw data CPT ratio). Pachymetry measurements are stored in a grid-like fashion, for one measurement every 0.1 mm, from −7 mm to +7 mm in both the x and y directions. The correct coordinates for extracting the corneal thicknesses at specific radii and angles were obtained by means of simple trigonometry. The average corneal thickness was calculated for every five degrees from 0 to 360 degrees for radii ranging from 1 to 5 mm, in steps of 0.5 mm. Larger radii were not considered useful because measurements regularly had large areas of missing data (pachymetric data rarely reached the 14 mm diameter generously given by the Pentacam software). To explore which radius was the most useful to predict the ablation spherical equivalent from the CPT ratio, we fitted simple linear regression models to the data. As there was concurrence of measurements in both eyes per subject, we decided to only use one eye per patient, which was selected by random sampling. This was performed with the sample function using a seed for reproducibility. The best model performance was evaluated based on R-squared values, as well as on their root mean squared error and the Bayesian information criterion (BIC).

The same one-eyed data were randomly split into a training set (2/3 of the data) and a validation set (1/3 of the data). Receiver-operating characteristics (ROC) were fitted to the training data using the ROCR package [18], whereby “hyperopia” (ablation spherical equivalent greater than and equal to 0 dioptre) was defined as the positive characteristic. Accuracy was estimated as (true positive + true negative)/(predicted positive + predicted negative). Sensitivity was estimated as true positive/positive, and specificity was estimated as true negative/negative. Predictions were then made on the test data. This was repeated 30 times to obtain error estimations.

In a secondary research question, we assessed the performance of the prediction of the ablation spherical equivalent from the CPT ratio based on simple linear models developed from the one-eyed data. Visual exploration showed a marked linear correlation. Therefore, linear regression seemed a viable approach. The Breusch–Pagan test on the data residuals did not reveal significant heteroscedasticity (p value =0.63).

As the optical zone or ablation zone might influence the significance of the CPT ratio as a predictive factor, we also tested its significance as a predictive variable by including it as an independent measure in our model.

3. Results

Data were collected for 157 eyes of 81 patients. After the exclusion of the eyes that had another laser refractive surgery than LASIK (8 eyes), or if the Pentacam examination was performed less than six months after the LVC treatment (8 eyes), or if there was more than one laser treatment performed for the same eye (1 eye), 140 eyes of 73 patients remained and were used for random sampling of one eye per patient; therefore, 73 eyes of 73 patients remained for analysis (42% female; mean age 40.9; SD 12.8). The mean time of the Pentacam examination after laser treatment was one year (SD 0.2 years).

For the Pentacam UI CPT ratio (Figure 1), the mean cutoff value for the best distinction of previous hyperopic versus myopic ablation was 0.86 (SD 0.01). The mean accuracy was 0.89 (SD 0.06), the mean sensitivity was 0.7 (SD 0.22), and the mean specificity was 0.95 (SD 0.06). Sensitivity and specificity curves for each iteration are shown in Supplemental Figure 1.

Supplemental Figure 2 shows the linear regression for the raw data CPT ratio at different radii of semiperipheral corneal thickness, and the best radius for prediction was 2.5 mm. The mean cutoff value for the best distinction was 0.92 (SD 0.01). The mean accuracy was 0.96 (SD 0.04), the mean sensitivity was 0.87 (SD 0.17), and the mean specificity was 0.99 (SD 0.02).

Figure 2 shows the area under the curve (AUC), a measure of the accuracy of the prediction (1 being ideal). AUC was 0.97 for the Pentacam UI CPT ratio (Figure 2(a)) and 0.999 for the raw data CPT ratio (Figure 2(b)).

Figure 3 shows the linear regression between CPT ratio and ablation spherical equivalent. The variance of the residuals is slightly higher toward hyperopic values when using the Pentacam UI CPT ratio, showing a tendency of outliers in the hyperopic range with a CPT ratio indicating a myopic ablation (Figure 3(a)). However, when using the raw data CPT ratio, the variance is more uniformly distributed and smaller throughout the entire range, and almost no outliers can be seen. As such, the magnitude of the error is smaller (Figure 3(b)). The formula obtained from the linear regression model for the Pentacam UI CPT ratio (Figures 1 and 3(a)) was(2) ablation spherical equivalent in dioptres=−49.14+57.62 x CPT ratio.

Table 1 shows the ablation spherical equivalent predicted from the Pentacam UI CPT ratio values comprised between 0.67 and 0.95 as per the previous equation.

The Pentacam UI CPT ratio is a good measure to predict the ablation spherical equivalent (R2 = 0.83; p value <0.001) (Figure 3(a)). The prediction performance increased considerably when using the manually computed raw data CPT ratio at 2.5 mm radius (R2 = 0.93; p value <0.001) (Figure 3(b)). This means that the raw data CPT ratio at 2.5 mm radius explained 93% of the variance of the ablation spherical equivalent.

The optical zone diameter ranged from 6.0 to 7.0 mm (mean 6.4 mm; SD 0.2 mm), and the ablation zone diameter ranged from 6.8 to 7.8 mm (mean 7.6 mm; SD 0.4 mm). As shown in Supplemental Figure 3, both optical and ablation zones had no clear linear relationship with the ablation spherical equivalent; thus, we decided to not include them as independent variables in our linear model. As the performance of our simple linear model was already adequate, we decided not to fit an exponential model to the data.

4. Discussion

Intraocular lens selection is a critical aspect of cataract surgery planning and the most challenging for LVC patients for whom the current formulas still leave a likelihood of refractive surprise. Knowing the pre-LVC refractive status helps the prediction of cataract surgery refractive outcomes for LVC patients [14]. When using IOL calculation formulas such as Haigis-L and Barrett true K (no history), it is required to specify the type of prior LVC treatment (myopic versus hyperopic). In fact, one can observe the significant variations in IOL power suggestions, for the same patient's data, depending on the type of prior LVC treatment chosen at the time of IOL power calculation. To date, without records of the prior refractive status, the treating clinician is limited to subjective and arbitrary evaluation of sagittal and tangential curvature maps which could give clues if there might have been a prior myopic or hyperopic treatment [19]. However, this method of evaluation lacks standardisation and validation.

This pilot study shows that the post-LVC central to semiperipheral corneal thickness ratio (CPT ratio) can accurately predict the type of LASIK ablation profile (hyperopic vs myopic) that was used, as well as provide an estimate of the ablation spherical equivalent in dioptres. This could in fact be predicted from the ablation profiles used by most laser platforms with mostly central ablations for myopic treatments, therefore decreasing the CPT ratio, and inversely for hyperopic treatments consisting of more peripheral ablations, therefore increasing this ratio.

Our study shows that the CPT ratio can serve as a valid cutoff value to distinguish between prior myopic versus hyperopic LVC. This information is required even when using formulas not requiring pre-LVC refractive data such as Haigis-L and Barrett true K (no history). In addition, studies have shown that additional knowledge of pre-LVC refractive status increases the accuracy of IOL calculations [10, 14]. The pre-LVC spherical equivalent can be estimated from the manifest refraction at the time of cataract surgery planning and the estimated ablation spherical equivalent given by the CPT ratio as follows: (3) Estimated pre−LVC SED=Manifest refraction SED+Ablation SED estimated from the CPT ratio.

The ablation spherical equivalent could be reasonably estimated from the CPT ratio alone. This ratio can be calculated by cataract surgeons using corneal thicknesses displayed on the Pentacam user interface software or the corneal tomographer of their choice. Alternatively, it can be calculated from the raw pachymetry data for improved accuracy. Thus, using the CPT ratio can help clinicians make a decision affecting the refractive outcome of cataract surgery in LVC patients without historical refractive data by increasing the accuracy of biometric predictions.

Hashemi et al. [20] observed in an Iranian population that the corneal thickness decreased with age by 0.65 μm per year peripherally and 0.25 μm per year centrally. As a result, adjusting the CPT ratio for age of the population might help improve the accuracy of the prediction of the ablation spherical equivalent. However, we do not yet fully understand whether similar changes would be expected in ageing eyes after LVC.

Alternatively, IOL calculation methods such as ray tracing do not require historical data and do not depend on keratometry indices. Ray tracing allows the calculation of the real anterior and posterior corneal curvatures and is not associated with measurement radius issues as calculations can be done on any diameter. Also, ray tracing allows corneal aberrations to be included in the IOL calculations [2, 21]. Unfortunately, so far, ray tracing has shown controversial results in the literature and is not yet widely available [22, 23].

Intraoperative aberrometry has been used to improve IOL choice during cataract surgery [9, 24]. It has been shown that for LVC patients without historical data, intraoperative aberrometry using the Optiwave Refractive Analysis (ORA) System (WaveTec Vision Systems Inc, Aliso Viejo, CA) wavefront aberrometer improved the postoperative refractive prediction error. However, this system is also not widely available, is costly and has shown controversial results [25].

One possible constraint of the present study is that all data analysed were recovered from one surgical centre with one excimer laser using the same wavefront-optimised profile for LASIK treatments. Variations are to be expected between our dataset and other surgical centres using different lasers, surgical procedures (PRK, trans-PRK, LASEK, SMILE, etc...) and ablation profiles (conventional, wavefront-optimised, wavefront-guided, topography-guided, Q-based, etc...), and those variables will change with technological advances. Further studies with data from multiple centres using different lasers, treatment modalities and ablation profiles are needed to validate our results in other settings.

In addition, it has been shown previously that the excimer laser stromal ablation and the flap side cut in the case of LASIK both induce biomechanical changes in the cornea inducing refractive changes which are, to date, difficult to predict with a current lack of a comprehensive biomechanical model. When collagen fibres and lamellae are cut, stress redistribution among the uncut fibres induces a hyperopic shift. As such, the stroma thins where the lamellae are cut and the surrounding stroma thickens [26].

It must be noted that, as shown previously, the epithelial thickness distribution changes after LVC. In particular, the epithelium overlying ablation areas thickens proportionally to the ablation depth and the anterior corneal curvature flattening rate. Conversely, the epithelium thins overlying areas of anterior corneal curvature steepening. Thus, epithelial remodelling participates in the regression toward the prior ametropia. [27–29]. Therefore, using a device (e.g. anterior segment OCT, very high-frequency ultrasound) capable of subtracting the epithelial thickness from the total corneal thickness and thus allowing the stromal thickness analysis could improve the CPT ratio's accuracy. As a result, this could improve the accuracy of the distinction between prior hyperopic and myopic LVC and their estimated ablation spherical equivalent.

Furthermore, tomography measurements show interdevice variability which in turn can affect the predictability of the CPT ratio [30]. More studies are needed to establish a reliable CPT ratio using other corneal tomographer devices than the Pentacam. In addition, corneal tomographers also show, to some extent, intradevice measurements variability [31].

As shown in our analysis of the optical/ablation zones, we could not find a significant correlation with the CPT ratio, nor with the ablation spherical equivalent. Also, as previously described in the literature, most treatment planning systems of excimer lasers use optical zones ranging between 6 and 6.5 mm with a general trend toward larger optical zones. Those are comparable to the optical zones used in this study [32–35]. Therefore, the principle of the CPT ratio may be applied to excimer lasers from other manufacturers than the one used in this study.

In conclusion, we have shown that the CPT ratio can differentiate, with reasonable accuracy, between a prior myopic versus hyperopic LASIK treatment and even predict the ablation spherical equivalent of such treatments. Whilst the present study did not include any presbyopic treatment, similar considerations would apply. Incorporation of this metric into current corneal tomographers, anterior segment OCTs, biometers, and other devices allowing analysis of the corneal thickness profile could help inform clinicians regarding prior LVC and be utilised to improve cataract surgery outcomes in LVC patients without historical refractive data or surgical records.

Acknowledgments

The authors would like to thank surgeons at OCL Vision (London, UK. https://www.oclvision.com/), who kindly provided all refractive and pachymetry data necessary for this study. The authors would also like to thank Dr Samuel Arba Mosquera (Schwind Eyetech Solutions., GmbH) for critical information regarding excimer laser ablation profiles.

Data Availability

The data used to support the findings of this study are included within the article.

Conflicts of Interest

The authors declare that they have no conflicts of interest.

Supplementary Materials

Supplementary Materials Supplemental Figure 1: postoperative central to semiperipheral corneal thickness (CPT) ratio, based on precomputed values provided by Pentacam UI software (A) and on extracted pachymetry raw data (B). Training was repeated on 2/3 of 30 independent random samples of the one-eyed data (one random eye for each participant). Supplemental Figure 2: ablation spherical equivalent in dioptres as a function of the postoperative central to semiperipheral corneal thickness (CPT) ratio, shown for different radii from 1 mm to 5 mm (one random eye per patient). CPT ratio was computed based on extracted pachymetry raw data. The best prediction performance was found at 2.5 mm. Supplemental Figure 3: ablation zone and optical zone relation to the ablation spherical equivalent. The ablation zone seemed to show a rather complex relation to the previous treatment (where y is not a function of x), and the optical zone vs the ablation spherical equivalent showed a rather sinusoidal relation. Thus, in both cases, there was no simple way of including them as variables in our linear model.

Figure 1 Pachymetry page of the Pentacam tomographer user interface (UI) software. The table at the bottom of the page shows the average pachymetry at the center of the cornea and at different diameters in the periphery. The values shown in red are used to calculate the CPT ratio in this study.

Figure 2 Receiver-operating characteristics (ROC) curves, based on precomputed values provided by Pentacam UI software output (a) and on extracted raw pachymetry data (b). Training was repeated on 2/3 of 30 independent random samples of the one-eyed data (one random eye for each participant). AUC: area under the curve.

Figure 3 Spherical equivalent of LASIK treatments as a function of the postoperative central to semiperipheral corneal thickness (CPT) ratio. (a) The CPT ratio calculated from the precomputed values of central and semiperipheral (at 3 mm radius) corneal thicknesses provided by the Pentacam UI software and (b) the CPT ratio calculated from the extracted raw data values of central and semiperipheral (at 2.5 mm radius) corneal thicknesses provided by the .csv export function of the Pentacam software. The computed mean cutoff is shown for both panels. Using raw data increased accuracy substantially (b) with dots in the extremes (lower left and upper right quadrants) showing less spread from the ideal linear regression.

Table 1 Ablation spherical equivalent (SE) predicted from the central to semiperipheral corneal thickness ratio (CPT ratio) calculated from the precomputed values obtained from the pentacam user interface software.

CPT ratio	SE (D)	
0.67	−10.5	
0.69	−9.4	
0.71	−8.2	
0.73	−7.1	
0.75	−5.9	
0.77	−4.8	
0.79	−3.6	
0.81	−2.5	
0.83	−1.3	
0.85	−0.2	
0.87	1.0	
0.89	2.1	
0.91	3.3	
0.93	4.5	
0.95	5.6	
D = dioptre; CPT = central to semiperipheral; SE = spherical equivalent.
==== Refs
1 Xia T. Martinez C. E. Tsai L. M. Update on intraocular lens formulas and calculations Asia-Pacific journal of ophthalmology (Philadelphia, Pa) 2020 9 3 186 193 10.1097/apo.0000000000000293 32501896
2 Preussner P.-R. Wahl J. Lahdo H. Dick B. Findl O. Ray tracing for intraocular lens calculation Journal of Cataract & Refractive Surgery 2002 28 8 1412 1419 10.1016/s0886-3350(01)01346-3 2-s2.0-0036688211 12160812
3 Carmona González D. Palomino Bautista C. Accuracy of a new intraocular lens power calculation method based on artificial intelligence Eye 2021 35 2 517 522 10.1038/s41433-020-0883-3 32346109
4 Li T. Stein J. D. Nallasamy N. AI-powered effective lens position prediction improves the accuracy of existing lens formulas medRxiv: The Preprint Server for Health Sciences 2020 106
5 Yamauchi T. Tabuchi H. Takase K. Masumoto H. Use of a machine learning method in predicting refraction after cataract surgery Journal of Clinical Medicine 2021 10 5 p. 1103 10.3390/jcm10051103
6 Darcy K. Gunn D. Tavassoli S. Sparrow J. Kane J. X. Assessment of the accuracy of new and updated intraocular lens power calculation formulas in 10 930 eyes from the UK National Health Service Journal of Cataract & Refractive Surgery 2020 46 1 2 7 10.1016/j.jcrs.2019.08.014 32050225
7 Connell B. J. Kane J. X. Comparison of the Kane formula with existing formulas for intraocular lens power selection BMJ open ophthalmology 2019 4 1 e000251 10.1136/bmjophth-2018-000251 2-s2.0-85066730805
8 Kane J. X. van Heerden A. Atik A. Petsoglou C. Accuracy of 3 new methods for intraocular lens power selection Journal of Cataract & Refractive Surgery 2017 43 3 333 339 10.1016/j.jcrs.2016.12.021 2-s2.0-85018349906 28410714
9 Savini G. Hoffer K. J. Intraocular lens power calculation in eyes with previous corneal refractive surgery Eye and vision (London, England) 2018 5 1 p. 18 10.1186/s40662-018-0110-5
10 Feiz V. Mannis M. J. Garcia-Ferrer F. Intraocular lens power calculation after laser in situ keratomileusis for myopia and hyperopia: a standardized approach Cornea 2001 20 8 792 797 10.1097/00003226-200111000-00003 2-s2.0-0034764626 11685053
11 Seitz B. Langenbucher A. Nguyen N. X. Kus M. M. Küchle M. Underestimation of intraocular lens power for cataract surgery after myopic photorefractive keratectomy Ophthalmology 1999 106 4 693 702 10.1016/s0161-6420(99)90153-7 2-s2.0-0033498941 10201589
12 Odenthal M. T. P. Eggink C. A. Melles G. Pameyer J. H. Geerards A. J. M. Beekhuis W. H. Clinical and theoretical results of intraocular lens power calculation for cataract surgery after photorefractive keratectomy for myopia Archives of Ophthalmology 1960 120 4 431 438 10.1001/archopht.120.4.431 2-s2.0-0036303315
13 Savini G. Hoffer K. J. Barrett G. D. Results of the Barrett True-K formula for IOL power calculation based on Scheimpflug camera measurements in eyes with previous myopic excimer laser surgery Journal of Cataract & Refractive Surgery 2020 46 7 1016 1019 10.1097/j.jcrs.0000000000000205 32271267
14 Ortiz S. Mena L. Rio-San Cristobal A. Martin R. Relationships between central and peripheral corneal thickness in different degrees of myopia Journal of optometry 2014 7 1 44 50 10.1016/j.optom.2013.03.005 2-s2.0-84894905491 24646900
15 AlMahmoud T. Priest D. Munger R. Jackson W. B. Correlation between refractive error, corneal power, and thickness in a large population with a wide range of ametropia Investigative Opthalmology & Visual Science 2011 52 3 1235 1242 10.1167/iovs.10-5449 2-s2.0-79955940089
16 Azar D. Refractive surgery 3rd edition - kniha 2021 https://www.elsevier.com/books/refractive-surgery/azar/978-0-323-54769-7
17 Project R. The R project for statistical computing 2021 https://www.r-project.org/
18 Sing T. Sander O. Beerenwinkel N. Lengauer T. ROCR: visualizing classifier performance in R Bioinformatics 2005 21 20 3940 3941 10.1093/bioinformatics/bti623 2-s2.0-27544491192 16096348
19 Shah R. S. Khandelwal S. S. Goshe J. M. Haberman I. D. Randleman J. B. Comparative postoperative topography pattern recognition analysis using axial vs tangential curvature maps Journal of Cataract & Refractive Surgery 2020 46 10 1368 1373 10.1097/j.jcrs.0000000000000264 32483077
20 Hashemi H. Asgari S. Mehravaran S. Emamian M. H. Shariati M. Fotouhi A. The distribution of corneal thickness in a 40- to 64-year-old population of Shahroud, Iran Cornea 2011 30 12 1409 1413 10.1097/ico.0b013e31822018dd 2-s2.0-81355160425 21975441
21 Gjerdrum B. Gundersen K. G. Lundmark P. O. Aakre B. M. Refractive precision of ray tracing IOL calculations based on OCT data versus traditional IOL calculation formulas based on reflectometry in patients with a history of laser vision correction for myopia Clinical Ophthalmology 2021 15 845 857 10.2147/opth.s298007 33664562
22 Lwowski C. van Keer K. Adas M. Ray-tracing calculation using scheimpflug tomography of diffractive extended depth of focus IOLs following myopic LASIK Journal of Refractive Surgery 2021 37 4 Thorofare, NJ, USA 231 239 10.3928/1081597x-20210114-02 34038658
23 Zhang J. Shao J. Zheng L. Shen Y. Zhao X. Comparative clinical accuracy analysis of the newly developed ZZ IOL and four existing IOL formulas for post-corneal refractive surgery eyes BMC Ophthalmology 2021 21 1 p. 231 10.1186/s12886-021-01991-7 34034681
24 Fram N. R. Masket S. Wang L. Comparison of intraoperative aberrometry, OCT-based IOL formula, haigis-L, and masket formulae for IOL power calculation after laser vision correction Ophthalmology 2015 122 6 1096 1101 10.1016/j.ophtha.2015.01.027 2-s2.0-84930043870 25766733
25 Christopher K. L. Patnaik J. L. Ifantides C. Time utilization and refractive prediction enhancement associated with intraoperative aberrometry use during cataract surgery Clinical Ophthalmology 2021 15 531 539 10.2147/opth.s287573 33603331
26 Dupps W. J. J. Wilson S. E. Biomechanics and wound healing in the cornea Experimental Eye Research 2006 83 4 709 720 10.1016/j.exer.2006.03.015 2-s2.0-33747884984 16720023
27 Reinstein D. Z. Archer T. J. Gobbe M. Silverman R. H. Coleman D. J. Epithelial thickness after hyperopic LASIK: three-dimensional display with Artemis very high-frequency digital ultrasound Journal of Refractive Surgery 2010 26 8 555 564 10.3928/1081597x-20091105-02 2-s2.0-78449242405 19928697
28 Reinstein D. Z. Archer T. J. Gobbe M. Change in epithelial thickness profile 24 hours and longitudinally for 1 year after myopic LASIK: three-dimensional display with Artemis very high-frequency digital ultrasound Journal of Refractive Surgery 2012 28 3 195 201 10.3928/1081597x-20120127-02 2-s2.0-84857890654 22301100
29 Kanellopoulos A. J. Asimellis G. Longitudinal postoperative lasik epithelial thickness profile changes in correlation with degree of myopia correction Journal of Refractive Surgery 2014 30 3 166 171 10.3928/1081597X-20140219-01 2-s2.0-84896947189 24576651
30 Turner J. M. Purslow C. Murphy P. J. Comparison of Javal-Schiøtz keratometer, Orbscan IIz and Pentacam topographers in evaluating anterior corneal topography Clinical and Experimental Optometry 2022 1–8 1 8 10.1080/08164622.2022.2067470
31 Gharieb H. M. Shalaby H. S. Othman I. S. p&gt;Repeatability and interchangeability of topometric, anterior chamber and corneal wavefront data between two scheimpflug camera devices</p&gt Clinical Ophthalmology 2020 14 3801 3810 10.2147/opth.s274303 33177806
32 Haw W. W. Manche E. E. Large optical ablation zone using the VISX S2 smoothscan excimer laser Journal of Cataract & Refractive Surgery 2000 26 12 1742 1747 10.1016/s0886-3350(00)00707-0 2-s2.0-0034536889 11134873
33 Shin D. H. Lee Y. W. Song J. E. Choi C. Y. Comparison of refractive outcomes after photorefractive keratectomy with different optical zones using Mel 90 excimer laser BMC Ophthalmology Jul. 2020 20 1 p. 270 10.1186/s12886-020-01537-3
34 Moshirfar M. Huynh R. Bundogji N. Comparison of 6.0 mm versus 6.5 mm optical zone on visual outcomes after LASIK Journal of Clinical Medicine 2021 10 17 p. 3776 10.3390/jcm10173776
35 Camellin M. Arba Mosquera S. Aspheric optical zones: the effective optical zone with the SCHWIND AMARIS Journal of Refractive Surgery 2011 27 2 135 146 10.3928/1081597x-20100428-03 2-s2.0-79951864433 20481411
