Skip to main content

Corneal power evaluation after myopic corneal refractive surgery using artificial neural networks

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).

Fig. 1
figure 1

Sample result of measuring the corneal power in the Pentacam device: a the results of measurement using the rotating Scheimpflug video camera, b cornea front and back, lens front and back, c anonymised patient data, d corneal thickness

Table 1 The values of measured parameters for 5 patients (all units, unless otherwise stated, are in dioptres)
Table 2 The mean, minimum and maximum values as well as standard deviations of the mean of the measured parameters for all 172 patients (the values of the feature w(1) are expressed in years, the values of the other features are expressed in dioptres)

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:

$$\delta \left( j \right) = \frac{1}{59}\mathop \sum \limits_{i = 1}^{59} \left| {w_{i,j}^{s} (11) - w_{ i} (11)} \right|$$
(1)

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.

Fig. 2
figure 2

Block diagram of the proposed algorithm for data analysis. This diagram is divided into three areas: data acquisition including the data of patients before refractive surgery, 1 month after surgery and 3 months later; data analysis including training and modifications of the structure of neural networks; the results obtained when the best features were selected

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.

Fig. 3
figure 3

Sample results of network testing: a 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). Each randomisation and calculation of the error are performed after 1000 iterations of network training; b the absolute difference between the power of the cornea predicted by the neural network and the reference power (ground truth) for 59 test cases

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).

Table 3 The first 10 smallest mean values of the error δ for 59 test patients for different configurations of features, j = 1000 randomisations and 1000 iterations for each randomisation and the number of neurons in the hidden layer s = 150

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:

Fig. 4
figure 4

Correlation matrix created for all the features w participating in obtaining the best results in accordance with Table 3. The artificial colour palette has been added to facilitate reading of the correlation values

  • 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).

    Fig. 5
    figure 5

    Graph of changes in the error δ(j) as a function of different numbers of neurons s in the hidden layer for all features w. The measurements were performed every 20 neurons for 500 iterations (j = 1)

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;

    Fig. 6
    figure 6

    Histogram obtained for 59 test patients for network configurations providing the best results. The values on the x-axis are the differences between the prediction results obtained using neural networks with back-propagation of errors and the ground truth values, whereas the values on the y-axis refer to the number of cases (number of patients)

  • 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

  1. 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.

    Article  Google Scholar 

  2. 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.

    Article  Google Scholar 

  3. Hoffer KJ. Intraocular lens power calculation after previous laser refractive surgery. J Cataract Refract Surg. 2009;35:759–65.

    Article  Google Scholar 

  4. 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.

    Article  Google Scholar 

  5. 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.

    Article  Google Scholar 

  6. 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.

    Article  Google Scholar 

  7. Ö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.

    Article  Google Scholar 

  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.

    Article  Google Scholar 

  9. 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.

    Article  Google Scholar 

  10. 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.

    Google Scholar 

  11. Olsen T. Improved accuracy of intraocular lens power calculation with Zeiss IOL Master. Acta Ophthalmol Scand. 2007;85:84–7.

    Article  Google Scholar 

  12. 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.

    Article  Google Scholar 

  13. Alio JL, Abdelghany AA, Buenaga RF. Management of residual refractive error after cataract surgery. Curr Opin Ophthalmol. 2014;25:291–7.

    Article  Google Scholar 

  14. De Bernardo M, Capasso L, Caliendo L, Paolercio F, Rosa N. IOL power calculation after corneal refractive surgery. Biomed Res Int. 2014;2014:658350.

    Article  Google Scholar 

  15. Hoffer KJ. Intraocular lens power calculation for eyes after refractive keratotomy. J Refract Surg. 1995;11(6):490–3.

    Google Scholar 

  16. 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.

    Article  Google Scholar 

  17. 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.

    Article  Google Scholar 

  18. 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.

    Article  Google Scholar 

  19. 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.

    Article  Google Scholar 

  20. 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.

    Article  Google Scholar 

  21. 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.

    Article  Google Scholar 

  22. 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.

    Google Scholar 

  23. 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.

    Article  Google Scholar 

  24. 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.

    Google Scholar 

  25. 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.

    Article  Google Scholar 

  26. 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.

    Article  Google Scholar 

  27. 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.

    Article  Google Scholar 

  28. 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.

    Article  Google Scholar 

  29. 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.

    Article  Google Scholar 

  30. 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.

    Article  Google Scholar 

  31. 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.

    Article  Google Scholar 

  32. 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.

    Article  Google Scholar 

  33. 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.

    Article  Google Scholar 

  34. 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.

    Article  Google Scholar 

  35. 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.

    Article  Google Scholar 

  36. 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.

    Article  Google Scholar 

  37. 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.

    Article  MathSciNet  Google Scholar 

  38. Porwik P, Wieclaw L. A new approach to reference point location in fingerprint recognition. IEICE Electron Express. 2004;1(18):575–81.

    Article  Google Scholar 

  39. 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.

    Google Scholar 

  40. 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.

    Google Scholar 

  41. 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.

    Google Scholar 

  42. Foster KR, Koprowski R, Skufca JD. Machine learning, medical diagnosis, and biomedical engineering research—commentary. Biomed Eng Online. 2014;13(94):1.

    Google Scholar 

  43. Nielsen M. Neural networks and deep learning. 2014 (free on line).

Download references

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

Authors

Corresponding author

Correspondence to Robert Koprowski.

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.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

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

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12938-016-0243-5

Keywords