Epithelial remodeling as basis for machine-based identification of keratoconus.

Purpose: To develop and evaluate automated computerized algorithms for differentiation of normal and keratoconus corneas based solely on epithelial and stromal thickness data.

Methods: Maps of the corneal epithelial and stromal thickness were generated from Artemis-1 very high-frequency ultrasound arc-scans of 130 normal and 74 keratoconic subjects diagnosed by combined topography and tomography examination. Keratoconus severity was graded based on anterior curvature, minimum corneal thickness, and refractive error. Computer analysis of maps produced 161 features for one randomly selected eye per subject. Stepwise linear discriminant analysis (LDA) and neural network (NN) analysis were then performed to develop multivariate models based on combinations of selected features to correctly classify cases. The sensitivity, specificity, and area under the receiver operating characteristic curve (AUC) were determined for each classifier.

Results: Stepwise LDA resulted in a six-variable model that provided an AUC of 100%, indicative of complete separation of keratoconic from normal corneas. Leave-one-out analysis resulted in 99.2% specificity and 94.6% sensitivity. Neural network analysis using the same six variables resulted in an AUC of 100% for the training set. Test set performance averaged over 10 trials gave a specificity of 99.5 ± 1.5% and sensitivity of 98.9 ± 1.9%. The LDA function values correlated with keratoconus severity grade.

Conclusions: The results demonstrate that epithelial remodeling in keratoconus represents an independent means for differentiation of normal from advanced keratoconus corneas.