efficiency of an algorithm derived from corneal tomography parameters to distinguish highly susceptible corneas to ectasia from healthy

Conditions: Keratoconus, Artificial Intelligence, Support Vector Machine

Summary: the objective of this study was to identify and build an algorithm through an imaging process using a support vector machine (svm) with the tomography variables of cases with, kc, highly susceptible corneas to ectasia (hsce), and healthy corneas and to compare this algorithm to bad-d (belin_ambrosio display) and prfi (pentacam random forest index). the study included 148 eyes with kc, 351 with healthy corneas, and 88 eyes with hsce.

patients were considered to be very asymmetric (vae-nt) if the diagnosis of ectasia was confirmed in one eye based on the previously described criteria and the fellow eye had a normal front surface curvature (topometric) map. objective criteria for considering normal topography was applied for defining the cases of vae-nt, including objective front surface curvature metrics derived from pentacam. normal topography was rigorously considered based on objective criteria of a maximum curvature kmax (steepest front keratometry) \<47.2 diopters, a paracentral inferior-superior (i-s value) asymmetry value at 6 mm (3-mm radii) \< 1.45, and keratoconus percentage index (kisa%) score \< 60. the cutoff point used to discriminate normal corneas and vae-nt from kc corneas was the maximal posterior elevation (\< 29 micron. this cutoff point had been determined in a previous study, using the same instrument and the same setting (30). the posterior elevation map was displayed with a 5-mm color-coded scale, and maximal posterior elevation was measured manually using the cursor in the central 5 mm. the study included 148 eyes with kc, 351 with healthy corneas, and 88 eyes with suspected kc. the patients were divided into three groups: * control group - normal eyes (cg): 351 eyes without kc of 351 patients who underwent lasik or photorefractive keratectomy (prk), stable after at least 18 months of follow-up, without any changes in the posterior elevation at the 18-month pentacam in relation to the preoperative exam (surgeries performed in 2012-2018). the inclusion criteria for being a normal case was to have normal corneas on the general eye examination in both eyes, including normal slit lamp biomicroscopy, corrected distance visual acuity of 20/20 or better, overall subjective normal topography and tomography examinations.. our objective topographic criteria were: both eyes with a kisa% index of less than 60%, kmax of 47.2 d or less, and i-s difference of less than 1.45 d. because no truly established tomographic parameter(s)/cut-off(s) for differentiating normal from keratoconus suspect eyes exist, we adapted our classification for normal eyes to the recent publication by ambrósio et al.(31) by adding the criterion of "overall subjective normal topography and tomography examinations" based on the evaluation of experienced refractive surgeon (gcaj). only one eye was randomly selected for further statistical analysis. the cg included one eye randomly selected from 323 patients with normal cornea; one eye was randomly included per patient to avoid selection bias related to the use of both eyes from the same patient * very assimetric ectasia with normal topography group (vae-nt g): 88 eyes of 88 patients with very asymmetric ectasia with normal topography (vae-nt) in one eye and frank ectasia (vae-e) in the fellow eye. the inclusion criteria followed previous studies (28, 32, 33) eyes in this group with insufficient topographic findings to meet diagnostic criteria for keratoconus, and following features normal-appearing cornea on slit-lamp biomicroscopy, keratometry, retinoscopy. these cases were the less affected eye (fellow eye) of a keratoconic patient was included if the following criteria were met: kisa% index of less than 60%, i-s difference of less than 1.45 d, and kmax of 47.2 d or less (ie, same topographic criteria as in normal eyes, except than in normal eyes, both eyes of the patient met the criteria). these patients can be considered with corneas highly susceptible to ectasia. * very asymmetric eyes with ectasia (vae-e): the fellow eyes of the vae-nt displaying a kisa% index of greater than 100% and at least one of the following biomicroscopic signs: vogt striae, fleischer ring, or focal stromal thinning. oculus implemented their own staging system into the pentacam software, which should mimic the amsler/krumeich systems: the topographic keratoconus classification (tkc) (34). tkc classifies kc into four stages (plus four intermediate stages) and identifies other corneal pathologies, such as corneal refractive surgery or pellucid marginal degenerative (pmd). tkc classification of pentacam tomographic index showing algum grau de kc variando de 1-4. * keratoconus group (kcg): 148 patients (one eye each) with bilateral clinical kc. the kcg included one eye randomly selected from 148 patients with keratoconus; one eye was randomly included per patient to avoid selection bias related to the use of both eyes from the same patient. the inclusion criteria were the same as for vae-e, except that both eyes of the patient met the ectasia criteria. all subjects underwent complete eye examination as well as refraction assessment, biomicroscopy, retinoscopy, fundoscopy, topography, and tomography assessment. all patients were assessed at the visum eye center between january 2012 and january 2018. this study adhered to the tenets of the declaration of helsinki and was approved by the research ethics committee of the sao jose do rio preto faculty of medicine. all patients were informed about the objectives of the study, and they signed written informed consent forms before being enrolled. external validation was conducted with 140 patients, whose data were not included in building the algorithm. they met the same inclusion criteria as the others, with a total of 82 eyes of 82 patients with healthy corneas, 19 eyes of 19 patients with vae-nt, and 39 eyes of 39 patients with kc. pentacam tomography: all eyes were examined by rotating scheimpflug corneal and anterior segment tomography (pentacam hr; oculus gmbh, wetzlar, germany). image quality was checked so that only cases with acceptable-quality images were included in the study. an experienced fellowship-trained corneal specialist (gcaj) reviewed all the cases so that they were correctly classified in the kc and vae-nt groups. the raw data (u12 files) were obtained from all cases; therefore, the same customized software (version 1.20r118) was used to process all the export files, and all scheimpflug variables were directly downloaded from the pentacam software using the "call-all" function. mathematical algorithm: to build the equation extracted from svm, 58 variables were used, some of them were extracted from the spreadsheet.. after the construction of these 58 feature vectors (fv), an svm-derived index was created, which was called the corneal tomography multivariate index derived from a support vector machine (ctmvi). considering that each patient represents a point on a cartesian plane with 58 dimensions (each coordinate representing one of the 58 fv), the role of svm is to find the hyperplane that best separates the cg, kcg, and vae-nt g subjects. a hyperplane is algebraically described by a linear equation; in this case, there are 59 coefficients, 58 of which are related to the fv and one independent coefficient representing the bias (which is a possible parallel dislocation of a given hyperplane). the analyzed fv were: arc (3 mm zone): anterior radius of curvature in the 3.0 mm zone centered on the thinnest location of the cornea; artmax: ambrosio relational thickness maximum; artmin: ambrosio relational thickness minimum; bad d: belin/ambrosio enhanced ectasia total deviation value ;bad daa: deviation of the art average; bad dam: deviation of the art max; bad db: deviation of back elevation difference map; bad de: deviation from the posterior elevation at the thinnest considering bfs 8 mm; bad df: deviation of front elevation difference map; bad df: deviation of minimum thickness; bad dk: deviation from kmax; bad dp: deviation of average pachymetric progression;bad dr: deviation from the more negative value on the relative thickness map; bad dy: deviation from the vertical displacement of the thinnest point from the apex; c.vol d 3mm: corneal volume of 3 mm diameter area; c.vol d 5mm: corneal volume of 5 mm diameter area; c.vol d 7mm: corneal volume of 7 mm diameter area; c.vol d 10mm: corneal volume of 10 mm diameter area; d2 mm / pachy min: the quotient of d2 mm / pachy min; d2 mm: average corneal thickness of 2 mm circle centered on the thinnest location; d4 mm / pachy min: the quotient of d4 mm / pachy min; d4 mm: average corneal thickness of 4 mm circle centered on the thinnest location; d6 mm / pachy min: the quotient of d6mm / pachy min; d6 mm: average corneal thickness of 6 mm circle centered on the thinnest location; d8 mm / pachy min: the quotient of d8 mm / pachy min; d8 mm: average corneal thickness of 8 mm circle centered on the thinnest location; ele b bfs 8 mm max. 4 mm zone: elevation parameter derived from the back surface centered at the point with highest value within the 4 mm (diameter) using the 8 mm best-fit sphere; ele b bfs 8mm apex: elevation parameter derived from the back surface centered at the apex calculated using the 8 mm best-fit sphere; ele b bfs 8mm thinnest: elevation parameter derived from the back surface centered at the thinnest point using the 8 mm best-fit sphere; ele b bfte 8 mm max. 4 mm zone: elevation parameter derived from the back surface centered at the point with highest value within the 4 mm (diameter) using the 8 mm best-fit toric ellipsoid; ele b bfte 8mm apex: elevation parameter derived from the back surface centered at the apex calculated using the 8 mm best-fit toric ellipsoid; ele b bfte 8mm thinnest: elevation parameter derived from the back surface centered at the thinnest point using the 8 mm best-fit toric ellipsoid; ele f bfs 8 mm max. 4 mm zone: elevation parameter derived from the front surface centered at the point with highest value within the 4 mm (diameter) using the 8 mm best-fit sphere; ele f bfs 8mm apex: elevation parameter derived from the front surface centered at the apex calculated using the 8 mm best-fit sphere; ele f bfs 8mm thinnest: elevation parameter derived from the front surface centered at the thinnest point using the 8 mm best-fit sphere; ele f bfte 8 mm max. 4 mm zone: elevation parameter derived from the front surface centered at the point with highest value within the 4 mm (diameter) using the 8 mm best-fit toric ellipsoid; ele f bfte 8mm apex: elevation parameter derived from the front surface centered at the apex calculated using the 8 mm best-fit toric ellipsoid; ele f bfte 8mm thinnest: elevation parameter derived from the front surface centered at the thinnest point using the 8 mm best-fit toric ellipsoid; iha: index highest asymmetry; ihd: index highest decentration; isv: index of surface variance; iva: index of vertical asymmetry; ki: keratoconus index; pachy min: corneal thickness at the thinnest point; pachy min y: position of minimum corneal thickness in relation of y axis centered on cornea apex; prc (3mm zone): posterior radius of curvature in the 3.0 mm zone centered on the thinnest location of the cornea; rel pachy min: relative corneal thickness at the thinnest point; rms hoa (cb): root mean square of high order aberration of cornea back; rms hoa (cf): root mean square of high order aberration of cornea front; rms hoa (cornea): root mean square of high order aberration of total cornea; rpiavg: average pachymetric progression index; rpimax: maximum pachymetric progression index; rpimin: minimum pachymetric progression index; z 3 -1 (cb): 3rd order vertical coma aberration cornea back; z 3 -1 (cf): 3rd order vertical coma aberration of cornea front; z 3 -1 (cornea): 3rd order vertical coma aberration total cornea; z 5 -1 (cb): 5th order vertical coma aberration of cornea back; z 5 -1 (cf): 5th order vertical coma aberration of cornea front. all zernike measurements were made for a corneal diameter of 6 mm.

Eligibility criteria:
inclusion criteria: patients were considered to be very asymmetric (vae-nt) if the diagnosis of ectasia was confirmed in one eye based on the previously described criteria and the fellow eye had a normal front surface curvature (topometric) map. objective criteria for considering normal topography was applied for defining the cases of vae-nt, including objective front surface curvature metrics derived from pentacam. normal topography was rigorously considered based on objective criteria (27, 28) of a maximum curvature kmax (steepest front keratometry) \<47.2 diopters, a paracentral inferior-superior (i-s value) asymmetry value at 6 mm (3-mm radii) \< 1.45, and keratoconus percentage index (kisa%) score \< 60 and (29). the cutoff point used to discriminate normal corneas and vae-nt from kc corneas was the maximal posterior elevation (\< 29 µm). this cutoff point had been determined in a previous study, using the same instrument and the same setting. the posterior elevation map was displayed with a 5-mm color-coded scale, and maximal posterior elevation was measured manually using the cursor in the central 5 mm. exclusion criteria: the following exclusion criteria were adopted: history of ocular trauma; chronic use of eye medication, glaucoma; dry eye syndrome; corneal scarring; neurotrophic keratopathy; severe meibomian gland dysfunction; vulnerable state owing to physical or mental illness and with language-related difficulties; pregnant or breastfeeding.

Keywords: keratoconus, ectasia, support vector machine, artificial intelligence