- Research
- Open access
- Published:
Corneal power evaluation after myopic corneal refractive surgery using artificial neural networks
BioMedical Engineering OnLine volume 15, Article number: 121 (2016)
Abstract
Background
Efficacy and high availability of surgery techniques for refractive defect correction increase the number of patients who undergo to this type of surgery. Regardless of that, with increasing age, more and more patients must undergo cataract surgery. Accurate evaluation of corneal power is an extremely important element affecting the precision of intraocular lens (IOL) power calculation and errors in this procedure could affect quality of life of patients and satisfaction with the service provided. The available device able to measure corneal power have been tested to be not reliable after myopic refractive surgery.
Methods
Artificial neural networks with error backpropagation and one hidden layer were proposed for corneal power prediction. The article analysed the features acquired from the Pentacam HR tomograph, which was necessary to measure the corneal power. Additionally, several billion iterations of artificial neural networks were conducted for several hundred simulations of different network configurations and different features derived from the Pentacam HR. The analysis was performed on a PC with Intel® Xeon® X5680 3.33 GHz CPU in Matlab® Version 7.11.0.584 (R2010b) with Signal Processing Toolbox Version 7.1 (R2010b), Neural Network Toolbox 7.0 (R2010b) and Statistics Toolbox (R2010b).
Results and conclusions
A total corneal power prediction error was obtained for 172 patients (113 patients forming the training set and 59 patients in the test set) with an average age of 32 ± 9.4 years, including 67% of men. The error was at an average level of 0.16 ± 0.14 diopters and its maximum value did not exceed 0.75 dioptres. The Pentacam parameters (measurement results) providing the above result are tangential anterial/posterior. The corneal net power and equivalent k-reading power. The analysis time for a single patient (a single eye) did not exceed 0.1 s, whereas the time of network training was about 3 s for 1000 iterations (the number of neurons in the hidden layer was 400).
Background
After the introduction of the excimer laser in the 1983, refractive surgery had a huge development and the number of patients who undergo this surgery is still growing [1]. Many of these patients had developed a significant cataract, therefore underwent phacoemulsification and intraocular lens (IOL) implantation [2]. The most important factor affecting accuracy of IOL calculating in these cases is the corneal power evaluation [3]. One of most diffuse devices for keratometry (corneal refraction, measurement of the radius of the anterior corneal surface curvature) is Pentacam HR. The Oculus Pentacam HR™ is a rotating Scheimpflug video camera that generates images from the anterior surface of the cornea to the posterior surface of the lens [4–6]. The number of images acquired from the Scheimpflug camera is closely dependent on the type of the Pentacam version. It can range from a few to a dozen of images, e.g. it can be performed at angles of 48°–228°, 55°–235°, 62°–242°, etc. [7, 8]. In addition to this, the following information is also obtained: corneal thickness, tangential curvature, axial/sagittal curvature, elevation, true next power, keratomic power deviation, anterior chamber depth, refractive power, refractive pachymetry, equiv K-reading power, total corneal refractive power [9, 10]. These information are typically obtained for patients before refractive surgery. Then, the tests are repeated after surgery. If the patient further develop a cataract, there is the needing of corneal power measurement to perform IOL power calculation. In naïve eyes, this type of measurement able to provide no, or very small, residual refractive defect after lens implantation so there is no need to wear eyeglasses [11–13]. However, the situation is quite different in the case of cataract patients who have been previously subjected to myopic refractive eye surgery: in these patients, the Oculus Pentacam and the other available devices do not allow for fully correct IOL power calculation [14, 15]. Therefore, many authors [16–18] are looking for an algebraic method using formulas (e.g. interpolation and approximation) [19, 20], which would help to reduce the residual defect. However, even if these methods improved a lot the accuracy of IOL power calculation, there still are improvement needing [9]. The relationship between the performed measurements (corneal thickness, tangential curvature, axial/sagittal curvature, elevation, true next power, keratomic power deviation etc.) and the lens power after corneal refractive surgery [21–24] is not precisely defined [25–36], so this article proposes the use of artificial intelligence methods. Of the various known methods of artificial intelligence, the authors have proposed artificial neural networks. In addition to other methods such as discriminant analysis, decision trees, k-menas [37, 38] and naive Bayes classifier [39–42], neural networks enable corneal power prediction on the basis of measurements performed with the Pentacam device. The proposed proprietary analysis algorithm and the results are presented later in this article.
Materials
Input data were obtained from 172 patients aged between 19 and 55 years (mean age of 32 years with the standard deviation of 9.4 years), including 67% of men. The study was conducted within routine tests performed in patients treated surgically in Centro Grandi Apparecchiature, Seconda Università di Napoli, Napels, Italy. All the patients were informed about the study, which was conducted in accordance with the Declaration of Helsinki. The measurements were made using the Pentacam HR keratometer with 1–19r11 software version (Fig. 1 shows examples of the results obtained) before refractive surgery, 1 month after the surgery and 3 months later. The following characteristics were measured: gender—feature w(1), age—feature w(2), left/right eye—feature w(3), cornea front/back-features w(4) and w(5) as the mean values of measurements for the angle of 90° and 180°, tangential anterior/posterior-features w(6) and w(7), analogously the mean values for the measurement in 4 points (similarly for the other features), net power of the cornea—feature w(8), total corneal refractive power—w(9), equiv K-reading power—feature w(10). They were acquired three times, namely prior to refractive surgery, 1 month after the surgery and 3 months later. Sample results obtained before the surgery for 5 patients are shown in Table 1. The mean, minimum and maximum values as well as standard deviation of the mean of the measured parameters for all patients are shown in Table 2. In addition to these data, the optical power of the cornea was calculated after the adjustment and correction of Pentacam errors—feature w(11) (Tables 1, 2). Feature w(11) will further constitute ground truth for the proposed structure of neural networks. The proposed new algorithm was implemented on a PC with Intel® Xeon® X5680 3.33 GHz CPU in Matlab® Version 7.11.0.584 (R2010b) with Signal Processing Toolbox Version 7.1 (R2010b), Neural Network Toolbox 7.0 (R2010b) and Statistics Toolbox (R2010b).
Methods
An approach using machine learning was proposed in the article. Initial tests confirmed the usefulness of neural networks with back-propagation of errors in predicting the value of feature w(11) on the basis of the features from w(1) to w(10). The proposed method was divided into two stages: training and testing. These two stages are directly related to the division of the entire database (the results obtained from 172 patients) into two sets: the training set containing 2/3 of the total number of patients, i.e. 113 patients, and the test set containing the remaining 1/3 of the patients, i.e. 59 patients. This type of division into the training and test set is typical of machine learning. In the present work, tests were conducted for neural networks with back-propagation of errors containing from 1 to 10 inputs and one output from the log-sigmoid transfer function. The networks provided one hidden layer. The tests were performed for different numbers of neurons in the hidden layer s from s = 10 to 500 every 10 neurons. The training process, unless stated otherwise, was performed in each case one million times. Every time after 1000 iterations, the order of individual patients in the training was random. In addition, every 1000 iterations, the error for the values obtained from the network \(w_{i,j}^{s} \left( { 1 1} \right)\) was compared to ground truth—w i (11), defined as:
where: \(w_{i,j}^{s} \left( { 1 1} \right)\)—the value obtained from the neural network output for i-th patient and j-th randomisation, w i (11)—the value of ground truth for i-th patient, s—the number of neurons in the hidden layer.
Since there is no evidence in the available literature [43] to what extent the features from w(1) to w(10) affect the outcome—tests were conducted (a million iterations) for each possible combination of the features w. In total, several billion iterations were performed for different combinations of features, different randomisations of cases in the training and test vector, and for different numbers of neurons in the hidden layer. A block diagram of the entire processing is shown in Fig. 2.
Examples of the results of network testing are shown in Fig. 3. Figure 3a shows the values of the error δ(j) for j randomisations (j∈(1,12000) and s = 150) and the features in the training vector w(3), w(4) and w(5). In the case shown in Fig. 3a, each randomisation and calculation of the error is performed after 1000 iterations of network training. Examples of the differences between the power of the cornea predicted by the neural network and the compared power (ground truth—see Table 1) for 59 test cases are shown in Fig. 3b. The graphs in Fig. 3 confirm the need for additional measurement of the maximum value δ max of the error δ(j) for j randomisations.
The minimum values are not significant from the point of view of analysis of results (since there is the mean value and standard deviation of the mean). This is due to the practical usefulness of the results of prediction and the attitude of patients who expect the same quality of vision (without glasses) as before cataract surgery.
Therefore, Table 3 shows the results of measurement of the mean value of the error δ(j) for j = 1000 randomisations, where for each randomisation there are 1000 iterations of network training and s = 150. The value of s equal to 150, which is the number of neurons in the hidden layer, was adopted on the basis of the measurements of the error δ(j) for all features. The relationships and participation of individual features were confirmed in the correlation matrix created for all features w which contributed to providing the best results (in accordance with Table 3).
As mentioned above, the tests were performed for different numbers of neurons in the hidden layer from s = 10 do 390 every 20 neurons. 500 iterations were performed every j-th randomisation. The graph in Fig. 5 shows changes in the error δ(j) as a function of different number of neurons in the hidden layer for all the features w. Detailed analysis of the results obtained and their interpretation are presented in the next section.
Results and discussion
The results obtained, especially Table 3 and Fig. 4, lead to some practical conclusions:
-
the most common feature in the top 10 results is w(10), which is equiv K-reading power;
-
the features w(1) and w(3) are not found in any of the 10 top results. In other words, the left/right eye and gender do not affect the obtained results and so they are not practically important for predicting the correct value of the optical power. A similar situation is with the patient’s age which is taken into consideration only in three cases (see Table 3);
-
the smallest number of features involved in correct corneal power calculation (with an error δ of less than 0.4 dioptres) is 3. These are the following combinations of features: w(6), w(7) and w(10) as well as w(7), w(8) and w(10). Medical interpretation of occurrence of these features results directly from the participation of tangential anterial/posterior, net power of the cornea, and equiv K-reading power in correct corneal power calculation [the feature w(11)];
-
the number of neurons in the hidden layer ranging from 10 to 590 provides the same accuracy of cornea power prediction (Fig. 5). There is only a change in the time of network training—the value of the variable j (for the same number of iterations).
The results can be compared to the results obtained by other authors. The publications on estimating the exact corneal power offer various formulas and corrections to be afterwards applied to the lens power calculation. Previous published papers about same topic showed corneal power evaluation after refractive surgery not to be reliable, presenting a wide range of errors, depending by the device tested and the study population [1–4, 9, 10, 27]. Compared to the results obtained in this study:
-
the presented methods are significantly different from each other in terms of mathematical formulas, which confirms that the artificial intelligence approach presented in this paper is correct;
-
the error obtained when using neural networks stems from the measurement noise which in turn results from the applied measurement method. Therefore, there is no possibility of analytical (assuming a random nature of methodology errors) reduction of the residual error in corneal power prediction to a value below 0.16 ± 0.14 dioptres. It is also impossible by using artificial intelligence—a histogram in Fig. 6 shows the exact value of the difference between prediction results obtained from neural networks with back-propagation of errors and the ground truth values;
-
the residual difference in corneal power evaluation obtained lead to no clinically significant variation of IOL to implant in eye facing cataract surgery after myopic refractive surgery.
-
analysis time for a single patient (single eye) for PC with Intel® Xeon® X5680 3.33 GHz CPU in Matlab® Version 7.11.0.584 (R2010b) with Signal Processing Toolbox Version 7.1 (R2010b), Neural Network Toolbox 7.0 (R2010b) and Statistics Toolbox (R2010b) does not exceed 0.1 s. At the same time, the time of network training is about 3 s for 1000 iterations (when the number of neurons in the hidden layer s = 400).
On this basis, it is possible to formulate the following conclusions presented in the next section.
Conclusions
The new approach to predicting the corneal power in patients previously undergone myopic refractive surgery can be characterized in the following points:
-
the proposed method of corneal power prediction is fully automatic and does not require operator intervention;
-
the presented method provides for 172 test patients (113 patients forming the training set and 59 patients in the test set) the error of corneal power prediction at the level of 0.16 ± 0.14 dioptres with a maximum value of 0.75 dioptres;
-
the presented method provides the above results for the following three input parameters (obtained with the Pentacam): tangential anterial/posterior, corneal net power and equiv K-reading power.
Corneal power prediction for IOL calculation in patients who had previously undergone surgical correction of myopic refractive defects should be made the same way as in relation to the optical measurement path. On the one hand, the applied approach and mathematical relationships implemented in the Penatacam tomograph require correction. Such correction is proposed by many authors. On the other hand, as shown in this article, artificial intelligence in the form of neural networks enables to reduce significantly the prediction error. Regardless of the results obtained, there is also the aforementioned random noise resulting from resolution and non-linearity errors (and also the measurement method itself). Thus, the presented method of artificial intelligence does not exhaust this very interesting topic which has many practical applications.
Abbreviations
- ILM:
-
internal limitans membrane
- RPE:
-
retinal pigment epithelium
- D:
-
dioptre
- std:
-
standard deviation of the mean
References
Aristeidou A, Taniguchi EV, Tsatsos M, Muller R, McAlinden C, Pineda R, Paschalis EI. The evolution of corneal and refractive surgery with the femtosecond laser. Eye Vis (Lond). 2015;14(2):12.
Klein BEK, Klein R, Lee KE, Gagnon RE. Incidence of age-related cataract over a 15-year interval; the Beaver Dam Eye Study. Ophthalmology. 2008;115:477–82.
Hoffer KJ. Intraocular lens power calculation after previous laser refractive surgery. J Cataract Refract Surg. 2009;35:759–65.
Hashemi H, Mehravaran S. Corneal changes after laser refractive surgery for myopia: com parison of Orbscan II and Pentacam findings. J Cataract Refract Surg. 2007;33(5):841–7.
Savini G, Barboni P, Profazio V, Zanini M, Hoffer KJ. Corneal power measurements with the Pentacam Scheimpflug camera after myopic excimer laser surgery. J Cataract Refract Surg. 2008;34(5):809–13.
Read SA, Collins MJ, Iskander DR, Davis BA. Corneal topography with Scheimpflug imaging and videokeratography: comparative study of normal eyes. J Cataract Refract Surg. 2009;35(6):1072–81.
Özyol P, Özyol E. Agreement between swept-source optical biometry and Scheimpflug-based topography measurements of anterior segment parameters. Am J Ophthalmol. 2016;169:73–8.
McAlinden C, Khadka J, Pesudovs K. A comprehensive evaluation of the precision (repeatability and reproducibility) of the Oculus Pentacam HR. Invest Ophthalmol Vis Sci. 2011;52(10):7731–7.
Lanza M, Paolillo E, Gironi Carnevale UA, Lanza A, Irregolare C, Mele L, Bifani M. Central corneal thickness evaluation in healthy eyes with three different optical devices. Cont Lens Anterior Eye. 2015;38(6):409–13.
Bialer OY, Kaiserman I, Bahar I. Accuracy of Scheimpflug Holladay equivalent keratometry readings after corneal refractive surgery in the absence of clinical history. Ophthalmic Res. 2014;52(4):217–23.
Olsen T. Improved accuracy of intraocular lens power calculation with Zeiss IOL Master. Acta Ophthalmol Scand. 2007;85:84–7.
Behnding A, Montan P, Stenevi U, Kugelberg M, Zetterström C, Lundström M. Aiming for emmetropia after cataract surgery: Swedish National Cataract Register study. J Cataract Refract Surg. 2012;38:1181–6.
Alio JL, Abdelghany AA, Buenaga RF. Management of residual refractive error after cataract surgery. Curr Opin Ophthalmol. 2014;25:291–7.
De Bernardo M, Capasso L, Caliendo L, Paolercio F, Rosa N. IOL power calculation after corneal refractive surgery. Biomed Res Int. 2014;2014:658350.
Hoffer KJ. Intraocular lens power calculation for eyes after refractive keratotomy. J Refract Surg. 1995;11(6):490–3.
Rosa N, De Bernardo M, Borrelli M, Lanza M. New factor to improve reliability of the clinical history method for intraocular lens power calculation after refractive surgery. J Cataract Refract Surg. 2010;12(36):2123–21128.
Shammas HJ, Shammas MC. No-history method of intraocular lens power calculation for cataract surgery aftermyopic laser in situ keratomileusis. J Cataract Refract Surg. 2007;33(1):31–6.
Wang L, Tang M, Huang D, Weikert MP, Koch DD. Comparison of newer intraocular lens power calculation methods for eyes after corneal refractive surgery. Ophthalmology. 2015;122(12):2443–9.
Potvin R, Hill W. New algorithm for intraocular lens power calculations after myopic laser in situ keratolileusis based on rotating Scheimpflug camera data. J Cataract Refract Surg. 2015;41(2):339–47.
Diehl JW, Yu F, Olson MD, Moral JN, Miller KM. Intraocular lens power adjustment nomogram after laser in situ keratomileusis. J Cataract Refract Surg. 2009;35(9):1587–90.
Latkany RA, Chokshi AR, Speaker MG, Abramson J, Soloway BD, Yu G. Intraocular lens calculations after refractive surgery. J Cataract Refract Surg. 2005;31(3):562–70.
Jarade EF, AbiNader FC, Tabbara KF. Intraocular lens power calculation following LASIK: determination of the new effective index of refraction. J Refract Surg. 2006;22(1):75–80.
Rosa N, Capasso L, Lanza M, Furgiuele D, Romano A. Reliability of the IOLMaster in measuring corneal power changes after photorefractive keratectomy. J Cataract Refract Surg. 2004;30(2):409–13.
Rosa N, Furgiuele D, Lanza M, Capasso L, Romano A. Correlation of changes in refraction and corneal topography after photorefractive keratectomy. J Refract Surg. 2004;20(5):478–83.
Savini G, Barboni P, Carbonelli M, Ducoli P, Hoffer KJ. Intraocular lens power calculation after myopic excimer laser surgery: selecting the best method using available clinical data. J Cataract Refract Surg. 2015;41(9):1880–8.
Lanza M, Iaccarino S, Cennamo M, Lanza A, Coen G. New Scheimpflug camera device in measuring corneal power changes after myopic laser refractive surgery. Cont Lens Anterior Eye. 2015;38(2):115–9.
Ianchulev T, Hoffer KJ, Yoo SH, Chang DF, Breen M, Padrick T, Tran DB. Intraoperative refractive biometry for predicting intraocular lens power calculation after prior myopic refractive surgery. Ophthalmology. 2014;121(1):56–60.
Savini G, Bedei A, Barboni P, Ducoli P, Hoffer KJ. Intraocular lens power calculation by ray-tracing after myopic excimer laser surgery. Am J Ophthalmol. 2014;157(1):150–3.
Sanctis U, Damiani F, Brusasco L, Grignolo F. Refractive error after cataract surgery combined with descemet stripping automated endothelial keratoplasty. Am J Ophthalmol. 2013;156(2):254–9.
Yang R, Yeh A, George MR, Rahman M, Boerman H, Wang M. Comparison of intraocular lens power calculation methods after myopic laser refractive surgery without previous refractive surgery data. J Cataract Refract Surg. 2013;39(9):1327–35.
Arriola-Villalobos P, DÃaz-Valle D, Garzòn N, Ruiz-Medrano J, Fernández-Perez C, Alejandre-Alba N, GarcÃa-Feijoò J. Effect of pharmacologic pupil dilation on OLCR optical biometry measurements for IOL predictions. Eur J Ophthalmol. 2014;24(1):53–7.
Awwad ST, Kilby A, Bowman RW, Verity SM, Cavanagh HD, Pessach Y, McCulley JP. The accuracy of the double-K adjustment for third-generation intraocular lens calculation formulas in previous keratorefractive surgery eyes. Eye Contact Lens. 2013;39(3):220–7.
Saiki M, Negishi K, Kato N, Ogino R, Arai H, Toda I, Dogru M, Tsubota K. Modified double-K method for intraocular lens power calculation after excimer laser corneal refractive surgery. J Cataract Refract Surg. 2013;39(4):556–62.
Potvin R, Hill W. New algorithm for post-radial keratotomy intraocular lens power calculations based on rotating Scheimpflug camera data. J Cataract Refract Surg. 2013;39(3):358–65.
Saiki M, Negishi K, Kato N, Arai H, Toda I, Torii H, Dogru M, Tsubota K. A new central-peripheral corneal curvature method for intraocular lens power calculation after excimer laser refractive surgery. Acta Ophthalmol. 2013;91(2):e133–9.
Fram NR, 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–101.
Porwik P, Doroz R, Orczyk T. The k-NN classifier and self-adaptive Hotelling data reduction technique in handwritten signatures recognition. Pattern Anal Appl. 2015;18(4):983–1001.
Porwik P, Wieclaw L. A new approach to reference point location in fingerprint recognition. IEICE Electron Express. 2004;1(18):575–81.
Koprowski R, Lyssek-Boron A, Nowinska A, et al. Selected parameters of the corneal deformation in the Corvis tonometer. Biomed Eng Online. 2014;13(55):1.
Koprowski R. Automatic method of analysis and measurement of additional parameters of corneal deformation in the Corvis tonometer. Biomed Eng Online. 2014;13(150):1.
Koprowski R, Teper S, Wrobel Z, et al. Automatic analysis of selected choroidal diseases in OCT images of the eye fundus. Biomed Eng Online. 2013;12(117):1.
Foster KR, Koprowski R, Skufca JD. Machine learning, medical diagnosis, and biomedical engineering research—commentary. Biomed Eng Online. 2014;13(94):1.
Nielsen M. Neural networks and deep learning. 2014 (free on line).
Authors’ contributions
RK suggested the algorithm for signal processing and implemented it and analysed the data. ML and CI were responsible for carrying out the measurements in patients and making the database available. All authors read and approved the final manuscript.
Acknowledgements
None.
Competing interests
The author declares that he has no competing interests.
Availability of data and materials
The data supporting the conclusions of this article are included within the article. Any queries regarding these data may be directed to the corresponding author.
Consent for publication
Any individual persons data were not applicable.
Ethics and consent
Not applicable. No tests, measurements or experiments were performed on humans as part of this work.
Author information
Authors and Affiliations
Corresponding author
Rights and permissions
Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.
About this article
Cite this article
Koprowski, R., Lanza, M. & Irregolare, C. Corneal power evaluation after myopic corneal refractive surgery using artificial neural networks. BioMed Eng OnLine 15, 121 (2016). https://doi.org/10.1186/s12938-016-0243-5
Received:
Accepted:
Published:
DOI: https://doi.org/10.1186/s12938-016-0243-5