
==== Front
Sci Rep
Sci Rep
Scientific Reports
2045-2322
Nature Publishing Group UK London

39223223
71480
10.1038/s41598-024-71480-6
Article
Fuchs’ Endothelial Corneal Dystrophy evaluation using a high-resolution wavefront sensor
Belda-Para Carolina 12
Velarde-Rodríguez Gonzalo 3
Marichal-Hernández José G. jmariher@ull.edu.es

2
Velasco-Ocaña Miriam 1
Trujillo-Sevilla Juan M. 1
Alejandre-Alba Nicolas 3
Rodríguez-Ramos José M. 12
1 Wooptix S.L., R &D Section, 38204 La Laguna, Spain
2 https://ror.org/01r9z8p25 grid.10041.34 0000 0001 2106 0879 ESIT, Industrial Engineering Department, Universidad de La Laguna, 38206 La Laguna, Spain
3 https://ror.org/049nvyb15 grid.419651.e 0000 0000 9538 1950 Ophthalmology Department, Hospital Universitario Fundación Jiménez Díaz, 28040 Madrid, Spain
2 9 2024
2 9 2024
2024
14 203692 7 2024
28 8 2024
© The Author(s) 2024
2024
https://creativecommons.org/licenses/by-nc-nd/4.0/ Open Access This article is licensed under a Creative Commons Attribution-NonCommercial-NoDerivatives 4.0 International License, which permits any non-commercial use, sharing, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if you modified the licensed material. You do not have permission under this licence to share adapted material derived from this article or parts of it. The images or other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by-nc-nd/4.0/.
This study aims to evaluate the applicability of the high-resolution WaveFront Phase Imaging Sensor (WFPI) in eyes with Fuchs’ Endothelial Corneal Dystrophy (FECD) through qualitative and quantitative analysis using a custom-designed Automatic Guttae Detection Method (AGDM). The ocular phase was measured using the t·eyede aberrometer and then was processed to obtain its High-Pass Filter Map (HPFM). The subjects were pathological and healthy patients from the Fundación Jiménez-Díaz Hospital (Madrid, Spain). The AGDM was developed and applied in pupils with 3 and 5 mm of diameter. A set of metrics were extracted and evaluated like the Root-Mean-Square error (RMS), Number of guttae, Guttae Area, and Area of Delaunay Triangulation (DT). Finally, a Support Vector Machine (SVM) model was trained to classify between pathological and healthy eyes. Quantitatively, the HPFM reveals a dark spots pattern according to the ophthalmologist’s description of the slit-lamp examination of guttae distribution. There were significant statistical differences in all the metrics when FECD and Healthy groups were compared using the same pupil size; but comparing both pupil sizes for the same group there were significant differences in most of the variables. This sensor is a value tool to objectively diagnose and monitor this pathology through wavefront phase changes.

Keywords

WaveFront sensor
Ocular aberrations
Fuchs' endothelial corneal dystrophy
Guttae
Image processing
Machine learning
Subject terms

Eye diseases
Corneal diseases
Biomedical engineering
Scientific data
Computer science
Imaging and sensing
Optical sensors
Characterization and analytical techniques
Imaging techniques
Wooptix S.L.issue-copyright-statement© Springer Nature Limited 2024
==== Body
pmcIntroduction

Fuchs’ Endothelial Corneal Dystrophy (FECD) is a bilateral, asymmetric, and low progression corneal disease characterized by morphologic changes in the endothelial cells (pleomorphism or polymegatism), the presence of guttae (excrescences from the Descemet Membrane) as its main hallmark, and edema in the severe stages1–3. The ophthalmologist Ernst Fuchs was the first to observe central corneal clouding in 1910th, who described this event as epithelial dystrophy4. Later studies based on this observation found that the endothelium had the predominant role because of the referred guttae1,5,6. This disease usually appears at the age of fifty years, and its progression will take 10-20 years, being more common in females than in males7–9. The literature states that it is associated with inherence and other ocular conditions such as glaucoma, keratoconus, or keratitis, while the prevalence varies among endemic populations, being higher in Caucasians than in the Asian population8,10. Although in mild stages the treatment is hypertonicity drops or pomade to reduce the appearance of edema, corneal transplantation is the procedure that restores the correct vision and reduces gradual complications7.

The usual exploration of this disease is performed using the slit-lamp to evaluate the status of the cornea (the presence of guttae, bullae, or corneal scarring). This device is employed to gradate the severity of the FECD by stages based on Krachmer’s classification regarding the number of guttae observed and their confluence, such as the presence of edema in the most severe stage11,12, as can be seen in Figure 2 of the Matthaei’s work13. Nevertheless, this classification depends on the subjective physician’s criteria and the illumination method. This fact decreases the diagnosis accuracy and viability, and highlights the importance of additional objective techniques to determine severity12,14,15. Several authors expressed the relevance of a quantitative observation of guttae distribution due to its benefit in the diagnosis and prognostic purposes, the effect of guttae on the quality of vision, and determining surgical indications for endothelial surgery in eyes with mild FECD16,17. Therefore, many studies have quantitatively assessed the distribution of guttae using slit-lamp images15,18, applying imaging treatment and morphological segmentation algorithms.

Additionally, another instrument commonly used in clinical practice for the evaluation of corneal status with this pathology is the specular microscope, which is also used to assess the endothelial cell morphology and guttae influence3,8,13, acquiring images of the endothelial cells mosaic of a concrete area of the cornea of 0.2 x 0.5 mm19. Many studies have quantitatively assessed the distribution of guttae using these sections16,17,20, but, when the disease is in its moderate stage, the guttae are so confluent that the specular microscope images are not valid.

Furthermore, some of these works enhance the relevance of guttae distribution on the quality of vision17,21, which is being explored through ocular aberrations22–24. Consequently, these studies reveal increased trefoil, coma, and high-order aberrations in FECD eyes using the estimated values from Pentacam elevation results (anterior and posterior corneal aberrations). Still, they report that the sources of the posterior corneal surface are unknown and are not able to distinguish whether guttae are the cause of this25. In the current literature, studies on aberrations in the FECD treat the phase only using Zernike polynomials (up to 6th or 8th order), and most are reconstructed from elevation data and not from direct phase measurements, leaving details hidden behind these orders that may be of utmost importance.

Recently, studies have been published using the WaveFront Phase Imaging (WFPI) sensor to compare keratoconic eyes with healthy eyes, which, apart from showing good reliability to obtain Zernike coefficients of up to 10th order26, is also able to discover in the ocular phase characteristic structures of these pathological eyes (banding patterns) that are not appreciated in healthy eyes, also analyzed quantitatively27. Hence, evaluating different ocular conditions is motivated because of its high lateral resolution (8.55 μm), and its reproducibility and repeatability already demonstrated against the gold-standard Hartmann-Shack28. Also, the technique implemented in this device has been validated in other sectors such as semiconductor metrology29, optical quality assessment30, and phase microscopy applied to biological samples31. All these fields require a high measurement precision that this sensor meets since it is capable of measuring Optical Path Differences (OPD) of a few nanometers caused by manufacturing defects or small refractive index changes on the same structure. Therefore, FECD is a clear candidate to be evaluated with this high-resolution sensor to unravel questions about how guttae affect the ocular wavefront phase.

The aim of this work is to objectively evaluate the applicability of the mentioned sensor using a custom-designed Automatic Guttae Detection Method (AGDM) capable of distinguishing between healthy and pathological eyes with FECD using phase information (concretely using the high-pass filter map).

Material and methods

Device: t·eyede aberrometer

In this study, the ocular phase is measured using the WFPI sensor (Wooptix S.L., La Laguna, Tenerife, Spain) implemented in the t·eyede aberrometer thoroughly described in previous works28,32,33. The sensor captures two intensity images shifted at the same distance (6.35 mm) on both sides of the conjugated pupil plane to recover the ocular phase, as shown in Fig. 1. The main advantage of this wavefront phase sensor is its 8.55 μm of lateral resolution obtained according to the number of pixels of a Charged-coupled device (CCD) sensor, which are converted into reliable measurement points. The system has also been optimized to prevent spurious light and speckle from reaching the sensor, and the numerous iterations performed by the algorithm make the phase recovery resistant to noise28,31.Fig. 1 WFPI sensor principle applied in ocular aberrometry.

This device uses a Superluminescent Diode (SLD) as a light source with a wavelength of 780 nm and an exposure time in the pupil plane of 30 ms with a power of 0.78 mW (value under the security limit of the ANSI rules). A Badal system is included to correct patient defocus over a range of − 10.00 to +10.00 D. Five measurements are taken with a few seconds of difference during the phase capture session, so the patient must look at the Maltese Cross used as the fixation target and may rest on the chin-rest provided on the set-up. Examples of an intensity image of both healthy and FECD eyes were shown in Fig. 2.Fig. 2 One of the intensity images of a healthy eye (a), and eye with FECD (b).

Once the measurement is completed, the intensity images are processed using the WFPI software, which recovers the ocular phase of the patient’s natural pupil with high resolution. This process has been detailed (and exemplified) in the previous study with healthy eyes and with keratoconus27. In summary, from this phase is extracted a phase map composed of the first 65th Zernike’s polynomials without the coefficients of piston, tip/tilt, and defocus (this map is considered as Low-Order Map). Moreover, the rest of the phase displayed (from the 66th coefficient and beyond) is what we name the Extremely High-Order Map because this phase map contains the high-frequencies. Then, a Gaussian high-pass filter with a cut-off frequency of 25 μm is applied to the whole raw phase (the same filter is used for all eyes to ensure comparability), obtaining hence the High-Pass Filter Map (HPFM). This study is based on the latter map, as it represents the OPD caused by small ocular structures not visible due to the main wavefront change due to low-order ocular aberrations.

Patients and measurement protocol

This study includes both healthy (control group) and FECD eyes (pathological group) according to five ophthalmologists. The pathological group is composed of: isolated FECD, FECD with cataracts, and FECD with Monofocal Intraocular Lens (IOL) cases. Dilated and undilated pupils were measured, in case of dilated pupils, it was on purpose for their previous ophthalmologic consultation and not for the t·eyede study. The measurements were performed at the Hospital Universitario Fundación Jiménez Díaz (Madrid, Spain). Its Institutional Review Board approved the study protocol that complied with the tenets of the Declaration of Helsinki. Informed consent was voluntarily obtained from the subjects after the explanation of the nature of the study previous t·eyede measurements. The exclusion criteria were presenting other ocular pathologies, being under the legal age, being currently pregnant or breastfeeding, or participating in another interventional study within 30 days before the start of this study.

To include patients within this study, ophthalmologists relied on slit lamp examination and specular microscope sections. With the slit lamp they determined whether an eye was pathological or not according to the appearance of guttae and their distribution along the cornea, and if so, some ophthalmologists graded its severity according to Krachmer’s classification. Also, they evaluated the state of the corneal endothelium with the specular microscope, where, through its algorithm, the number of endothelial cells that remain is counted (as an example of this algorithm is shown in Fig. 3b). Since the Hospital Fundación Jiménez Díaz is a clinic external to Wooptix, what it was able to provide to the company for performing this study were the slit lamp examination descriptions and some specular microscope images (examples of different cases with the stages assigned by the ophthalmologists presented in Fig. 3).Fig. 3 Specular Microscope sections of: (a) Healthy, (b) FECD Grade 2, and (c) FECD Grade 5 eyes according to the ophthalmologists’ criteria.

The control group is formed by 43 healthy eyes of 25 patients (range age: 27–56 years; Gender: 11 Males, 14 Females) and the pathological group has 204 eyes of 120 patients (range age: 40–93 years; Gender: 39 Males, 81 Females). We do not consider any other demographic data relevant to this study.

Automatic Guttae detection method

As all the eyes diagnosed with FECD by the ophthalmologist presented a dark spots pattern in the HPFM, in order to objectivate the evaluation of this disease using the WFPI technique, we have developed an Automatic Detection Method based on image processing by morphological segmentation. Due to the number of images per eye (up to five valid images in the same session) and the large number of guttae per image, this method was made up with three objectives in mind: that the number and location of guttae detected for an eye should coincide in the different measures, that the method should work for the several levels of evolution of the disease, and that it should not require excessive parameters whose fine-tuning would affect the repeatability of the measurements.Fig. 4 Diagram depicting the internal processes, I/O and parameters of the Automatic Guttae Detection Method.

The flow of operations of the method and the parameters that govern it are shown in the diagram in Fig. 4. It depends on several fixed internal parameters but only on two external modifiable ones, Minimum Brightness and CrossCorrelation Threshold. The entire experiment was performed on all images and for a wide range of combinations of both parameters. Those parameters that gave the best results were retained. In this manner, the method is self-adjusted to the depth absolute values that t·eyede provides without user intervention.

The first part of the method consists of cutting the pupil diameter to 5 and/or 3 mm, as well as applying masks to cover unwanted reflections that may affect the phase representation in the HPFM, called preprocessing because it is prior to the image processing performed by the algorithm. The next step consists of thresholding the input image to select the pixels with an intensity lower than the self-adjusted parameter Minimum Brightness. Those pixels are then segmented into neighboring zones, from which the area is calculated. Zones whose areas lie between the 10th and 40th percentiles are clearly considered guttae. A gutta profile model is formed for the eye in consideration by averaging the areas that have passed the filters described up to this point. Then, a normalized cross-correlation is performed between the gutta model and the eye image. The cross-correlation output is thresholded with the second self-adjusted parameter of the method: CrossCorrelation Treshold. Finally, a non-maximal suppression is applied to keep the central location of each detected gutta.Fig. 5 A depiction of the inner workings of the method advancing from top to bottom. See text for details.

Figure 5 shows the different stages of the method. Firstly, Fig. 5a is an input eye image with the region of interest surrounded by a red rectangle (subsequently analyzed in detail), whereas Fig. 5b ghts those zones whose intensity is superior to the Minimum Brightness. Figure 5c shows the clear guttae candidates. Figure 5d represents the gutta shape model. Figure 5e shows the output of the normalized cross-correlation of the eye image with the gutta model and highlights the zones that exceed the parameter CrossCorrelation Threshold. Notice that depending on that threshold, the number of final guttae can be greater or lower than the clear guttae candidates. As previously discussed, the method self-adjusted this parameter. Finally, Fig. 5f reveals each peak that succeeds in the non-maximal values suppression, that is, the final location of each gutta center. Hence, the surroundings of the guttae centers are analyzed to extract the metrics reported below.

Metrics description and analysis

A set of metrics is extracted to quantitatively compare healthy and pathological eyes, as well as to train an automatic binary classification model between both groups using pupil sizes of 3 and 5 mm. The calculated metrics are: the Root Mean Square error (RMS) (mrad) considering the values inside the pupil of the HPFM, the Number of Guttae Detected (accounts), Guttae Density (calculated as the number of guttae per mm2), Mean Guttae Diameter (μm), Mean Guttae Area (μm2), OPD Height (OPDH) (mrad) considered as the difference between the highest and the lowest point of the guttae profile, and Phase Roughness (mrad) calculated by equation 1 based on the Verma’s work equation to obtaining the roughness in ablation in corneas with Laser Vision Correction34.1 PhaseRoughness=RMSOPDH

In addition, Delaunay Triangulation (DT) is obtained regarding the guttae position, extracting as quantitative metrics the mean triangle area and the mean triangle area calculated after applying the Convex Hull algorithm. This mathematical tool consists of constructing triangles from a set of points (in this case, the located position of guttae) where the circumcircle of each triangle that forms the mesh must not contain any vertex of another triangle35. Therefore, when guttae are distant and scattered, the triangle area would increase (like in a healthy eye) whereas, if guttae are agglomerated (like in a FECD eye), the triangle area would decrease due to the small distance between guttae. Moreover, the Convex Hull algorithm is implemented to not ponder those external points of the map considered as false negative values.

Once the values have been obtained, the results of both groups and both sizes of students are statistically compared. For this purpose, first the Shapiro-Wilk test is used to assess the normality of each sample, and then the Mann-Whitney U test is used in case any of the samples do not meet the normality distribution. Statistical analysis was performed using the package scipy.stats (v.1.7.1) of Python (v. 3.9).

Automatic binary classifier model

As a further consideration, this work proposes the Support Vector Machine (SVM) tool from Matlab (v.R2022b) to create an automatic classifier to distinguish between healthy and pathological eyes. This method is applied to both pupil diameter analyses.

Previously, in order to avoid redundant information from the nine previous variables that can correlate with each other (like Mean Guttae Diameter and Mean Guttae Area), the Principal Component Analysis (PCA) was applied to generate a new set of variables (Principal Components(PC)) as a result of linear and orthogonal combinations of the original variables based on finding the maximum variance among the data set36,37, so the principal components will be ordered according to the largest amount of variance among the initial metrics. The number of PC is equal to that of the original set, and each eye (input) would be characterized by the scores of each PC. In this study, the SVM model’s behavior is evaluated using three and five Principal Components.

With the new representation of the sample, to train the SVM model, it is necessary to enter a training data set (composed of the 50% of each group) where each input was labeled as healthy or pathological (ophthalmologist criterion). As at least five measures per eye were taken during the measurement session, only one measure is randomly chosen to avoid the same eye appearing in both sets. Once the model is trained, it is applied to a test data set (the remaining 50%) to be the model itself which classifies each new input as healthy or pathological38,39.

Then, to compare if there are statistical differences between both classifications for the same eye, McNemar’s test is used. That is a statistical test for paired nominal data (the same eye has two independent classifications with the options: Healthy or FECD) where its null hypothesis states that the performance of both classifiers is equal (accepted when the p-value is higher than 0.05)40; also it is recommended to use when the machine learning model is only once fitted, like in the model example presented41. In addition, taking doctors’ diagnoses as the reference classification, to quantify the model hits and misses, the subsequent evaluation metrics for binary classification are obtained according to the equations presented in Rainio’s work42: Accuracy, defined as the proportion of correctly classified instances by the SVM model in the set of all the instances; Sensitivity, defined as the proportion of the truely positive instances from those cases that are pathological according to the doctors’ criteria; and Precision, defined as the instances classified as positive from those cases that are predicted as pathological by the SVM model.

Results

Qualitative evaluation of the Fuchs’ endothelial corneal dystrophy

Since part of the pathological eyes has not been graded into stages by ophthalmologists, it has been decided to present examples of Healthy eyes and FECD with and without severity grade in Fig. 6 to demonstrate the agreement between what the ophthalmologists described in their slit-lamp evaluations with what is observed in the HPFM, expressing both comparisons below. Figure 6a is a healthy eye with no remarkable event according to the ophthalmologist and any characteristic pattern in the HPFM; Fig. 6b is an isolated FECD eye with grade 2 and inferior-paracentral guttae (2 mm area) as reported by the doctor, also in the HPFM slight agglomeration of dark spots in the 2 mm below the center appear; Fig. 6c corresponds to an isolated FECD eye with grade 5 (any other comment from the ophthalmologist), while in the HPFM there are dark spots concentrated in aproximately the whole pupil (less frequent at the lower edge); Fig. 6d is another isolated FECD case with grade 5 with guttae in central 4–5 mm according to the slit-lamp exploration, as well as in the HPFM there are dark spots extended over approximately the entire pupil (edges less affected); Fig. 6e is a FECD eye with undefined grade and IOL, but the doctors referred central guttae, as can be seen in the HPFM, where dark spots are extended in the centre of the pupil with a zone of slight roughness change (right side) surrounded by the dark spots; at last, Fig. 6f is a FECD eye with undefined grade and incipient cataract without any comment from the ophthalmologist, but in the HPFM there are dark spots extended in the centre of the pupil with a zone of moderate roughness change (right side) surrounded by the dark spots.Fig. 6 HPFM of eyes regarding the opthtalmologist criteria: (a) Healthy left eye (OS), (b) FECD grade 2 right eye (OD), (c) FECD grade 5 left eye, (d) FECD grade 5 right eye, (e) FECD left eye (undefined grade) with Monocular Intraocular Lens, (f) FECD left eye (undefined grade) with incipient nuclear cataract. Representation with the patient’s natural pupil.

Furthermore, these maps can be directly comparable to those slit-lamp captures shown in Figure 2 of Matthaei and collegues13, and in slit-lamp retroilluminated images from other studies15,18,43.

Metrics results and comparisons

Regarding the metrics descriptors, their results are presented in Table 1. The values are obtained for both groups (Healthy and FECD) for both pupil diameter analyses (3 and 5 mm). The study of 3 mm includes 247 eyes (43 Healthy and 204 FECD), while the 5 mm analysis has 149 eyes (39 Healthy and 110 FECD). In this last case, the sample is lower because the patient’s natural pupil is smaller than such diameter. Table 1 Mean and standard deviation of all the metrics for both Healthy and FECD groups, and for both pupil diameter analyses (3 and 5 mm).

Metrics	3 mm	5 mm	
Healthy	FECD	Healthy	FECD	
RMS (mrad)	0.347 ± 0.066	0.81 ± 0.24*	0.334 ± 0.060*	0.77 ± 0.23*	
Num. Guttae (accounts)	136.1 ± 56.6*	341.1 ± 90.9*	367.7 ± 148.3*	908 ± 268*	
Guttae Density (accounts/mm2)	19.25 ± 8.00	48.3 ± 12.9*	18.73 ± 7.55	46.2 ± 13.7*	
Guttae Diameter (μm)	22.72 ± 2.12	26.34 ± 2.52	22.64 ± 1.57	26.42 ± 2.22	
Guttae Area (μm2)	409.1 ± 80.1	550.1 ± 104.7	404.4 ± 56.5	552.0 ± 91.7	
OPDH (mrad)	0.96 ± 0.33	1.32 ± 0.37*	0.90 ± 0.30	1.19 ± 0.35*	
Phase Roughness (mrad)	11.65 ± 3.25	23.33 ± 9.09	11.65 ± 3.23	23.24 ± 9.18	
Area DT (px2)	1,520 ± 930	555 ± 218*	1,625 ± 887	602 ± 309*	
Area DT with Convex Hull (px2)	1,559 ± 838*	614 ± 216	1,729 ± 925*	642 ± 276	
The * mark indicates that the p-value obtained in the U Mann–Whitney test is less than 0.05 for the same group comparing both pupil diameters. FECD, RMS, OPDH, and DT are the abbreviations of Fuchs’ Endothelial Corneal Dystrophy, Root Mean Square error, Optical Path Difference Height, and Delaunay Triangulation, respectively.

The samples have been evaluated statistically using the Shapiro-Wilk test to check their normality distribution. As most cases presented a p-value lower than 0.05, the hypothesis that the samples tend to the normal distribution is rejected, and the U Mann-Whitney test is applied to compare them between groups and pupil diameters. When both groups are compared using the same pupil diameter, the p-value is lower than 0.05 for all the metrics, indicating significant statistical differences between healthy and pathological eyes in both pupil diameters. Otherwise, when different pupil diameters are compared for the same group, the statistical differences are represented in Table 1 with the * mark.

In addition, to show the stability of the phase in the same capture session (five measures with a few seconds of difference), we add an example of a healthy eye, two FECD eyes at different stages of the disease, an one ungraded FECD eye according to the ophthalmologist criteria. Therefore, three HPFM of each eye are presented in Fig. 7, and Table 2 quantitatively demonstrates the variation between the different measurements of the same eye in the same session.Fig. 7 HPFM of different measures from the same eye in a single session: (a) Healthy eye, (b) FECD grade 2, (c) FECD grade 5, (d) FECD undefined grade with Monofocal Intraocular Lenses. Pupil size of 5 mm.

Table 2 Mean and Standard Deviation of a set of metrics calculated from the three images of each eye (different severity examples) in Fig. 7 to show the stability between measures using a pupil diameter of 5 mm. FECD, RMS, and DT are the abbreviations of Fuchs’ Endothelial Corneal Dystrophy, Root Mean Square error, and Delaunay Triangulation, respectively.

Description	RMS (mrad)	Num. Guttae (accounts)	Phase Roughness (mrad)	Area DT (px2)	
Healthy	0.382 ± 0.014	258.2 ± 53.6	11.60 ± 1.06	1,946 ± 323	
FECD Grade 2	0.490 ± 0.056	746.8 ± 40.42	16.71 ± 3.88	670.7 ±35.5	
FECD Grade 5	0.571 ±0.037	1,283.6 ± 71.6	14.60 ± 1.85	392.55 ± 19.86	
FECD Undefined grade	0.662 ± 0.045	1,039.5 ± 125.1	18.02 ± 2.38	490.8 ± 56.2	

Automatic binary classifier model evaluation

An example of the PCA results comparing the first three dimensions using the pupil diameter of 3 mm is represented in Fig. 8. In this case, the first principal component is mostly influenced by the results obtained with the metrics of Delaunay Triangulation (associated with the guttae distribution), the second principal component with the metric of Mean Guttae Area, and the third component is mostly influenced by the metric Number of Guttae Detected.Fig. 8 (a) 3D Plot of the first PCA dimensions; (b) Scores for the 1st against the 2nd PCA dimensions; (c) Scores for the 1st against the 3rd PCA dimensions.

Regarding the statistical comparison between both classifications applying McNemar’s test, the p-values obtained for the 3 mm pupil size were 0.67 and 0.51 using 3 and 5 PCs respectively, meanwhile for the 5 mm pupil size using 3 PCs was 0.56 and 0.66 with the 5 PCs. Then, the null hypothesis is accepted for all the cases, indicating that there are no statistical differences between both classifications.

On the other hand, the evaluation metrics of the SVM model have also been calculated for these 4 analyses considering the physicians’ labels as true. The accuracy, sensitivity, and precision results for the 3 mm pupil analysis with 3 PCs were: 92.8%, 96.9%, and 94.8% respectively, meanwhile using 5 PCs the results were: 92.8%, 96.1%, and 95.2%. Also, for the 5 mm pupil analysis with 3 PCs the values were: 90.2%, 94.8%, and 92.1% respectively, and 90.7%, 94.8%, and 92.7% using the 5 PCs.

Discussions

This study illustrates the applicability of the WFPI Sensor to evaluate Fuchs’ Endothelial Dystrophy by means of a qualitative and quantitative analysis using metrics extracted from the ocular wavefront phase. The usefulness of the current sensor demonstrates the relevance of considering the phase beyond extreme orders to uncover minor details such as guttae. Reviewing the literature, no previous aberrometer has revealed a pattern like these from the wavefront phase, being one of their main disadvantages representing the phase up to the 10th order of Zernike’s polynomials as well as recalculating the corneal wavefront from elevation data.

First of all, examples of a healthy and FECD eyes HPFM with different corneal status according to the ophthalmologist’s description are shown in Fig. 6. In these cases, there is concordance between what is described by the ophthalmologist in his examination using the slit lamp (distribution of guttae) with the pattern of dots observed in the HPFM. Likewise, these patterns are similar to the distribution of guttae shown by other authors with slit-lamp15,18,43, and as clear examples are the comparisons between our maps (d) and (f) in Fig. 6 with the slit lamp images (c) and (e) in Figure 2 of Matthaei’s work respectively. In the first example, the guttae are agglomerated in the center of the pupil but do not exist in the periphery; while in the second case there is an area of changing roughness surrounded by guttae that Matthaei describes as a bullae. On the other hand, given the novelty of the sensor and its applicability in ophthalmology, but previously validated in artificial and human eyes27,28, it has been decided to show examples of stability of the wavefront phase measurement between 5 captures (taken consecutively a few seconds apart with blinks in between) to demonstrate that it is not random noise generated by the sensor or type or fluctuating structures such as tear deposits or vacuoles due to their variable position (Fig. 7). We have also calculated the mean and standard deviation of some metrics between the measurements for each eye, showing that the standard deviation is not critical (Table 2).

Regarding the statistical differences obtained in the metrics extracted from the Automatic Guttae Detection Method, we can hypothesize that: the metric Number of Guttae Detected increments when the pupil diameter increases, but the metric Guttae Density does not change in the Healthy group, indicating that in this group the regions detected are not dependent of the pupil diameter while in the FECD group the value was slightly decreased due to the highest concentration of guttae in the central region that tends to spread toward the periphery, as stated in all reviews1,3,7–9. The metrics Mean Guttae Diameter and Mean Guttae Area are statistically different between healthy and pathological eyes independently of pupil diameter, but, for the same group, the value remains regardless of pupil size, indicating that each group has its own characteristics. In this case, although no literature characterizes these parameters with wavefront phase information, these results can be compared with other studies that use different methods to calculate it: Soh and colleagues obtained a mean value of guttae size of 407 ± 155 μm218, being closer to ours obtained in the healthy group rather than in the pathological group, and Dong and colleagues reported guttae diameters in a range of 20 to 60 μm44 in agreement with our values obtained in both groups. Concerning the OPDH metric, there were statistical differences in the FECD group expressing that it decreases when the external area is evaluated, launching the hypothesis that this could be associated with the dystrophy severity, while this metric remains constant in healthy eyes between different pupil sizes. The Phase Roughness does not change in any of the groups despite pupil size but differs from group to group, as was also illustrated in Vaitinadapoulé’s work, demonstrating that the roughness increases in the FECD eyes against the healthy ones measuring corneal nanotopography by chromatic confocal microscopy45. The results with DT metric match with those of the Guttae Density, indicating that the guttae distribution changes tend to spread along the pupil (explaining the triangle area increment) in the FECD group meanwhile the distribution is uniform in healthy eyes. Otherwise, when Convex Hull method is applied, the DT mean area illustrates that removing external points affects the Healthy group (increasing the triangle area) but not the FECD eyes due to the high amount of guttae also in the periphery because these spots are not considered as outsiders. Nevertheless, due to the different methodologies of acquisition, the quantitative comparisons between works are not completely reliable, only tentative.

After this detailed analysis, the usability of the automatic method is supported by the number of metrics that are easy to interpret and are distinguishable between healthy and pathological eyes. Moreover, this method is able to directly acquire more features than conventional images, such as phase roughness. This parameter should be evaluated in more detail in future studies, as it has been seen in some pathological eyes as there is a change of roughness in certain regions, so it should be evaluated individually as it could be key in the early detection of the appearance of edema in mild cases that could be missed in slit lamp examination. There are already studies focused on this aspect using the Anterior Segment Spectral Domain Optical Coherence Tomography (OCT) technique46, where the presence of edema at different stages of the pathology is mainly evaluated by means of a specific number of corneal transversal sections. But in these cases, the distribution of the guttae is not relevant, which should be taken into account since it has been shown the relation between the guttae distribution and corneal edema47. Moreover, Soh and colleagues in their Figure 418 represent the area with greater accumulation of guttae regarding the laterality of the eye, being clearly superior in the inferotemporal pupil zone, which agrees with the regions where the HPFM changes the roughness in right eyes (on the left side of the map), and left eyes (on the right side of the map), like are exemplified in Fig. 6e, f. With this in mind, the device presented in this study could quantitatively monitor pathology over time based on the distribution of guttae and can assess when this change in roughness occurs.

Notice that throughout the paper we have been talking about the Optical Path Difference Height because we want to make it clear that we are not measuring the morphological height or depth of these per se, but the phase change they produce. That justifies the scale with which the HPFM is represented, together with the fact that when recovering the phase by measuring the intensity propagation, the details that are below the axis resolution of the system (not the lateral resolution) are mitigated but not disappeared, as vanDam refers in his work48.

Concerning the PCA and the SVM model, the example ploted in the Fig. 8 with the first 3 PCs for 3 mm pupil diameter shows that both groups are in two separate regions, and the principal components that distinguish between both groups are obtained from the metrics that would be easily interpreted by ophthalmologists: guttae distribution, guttae area, and regions detected. The pathologic group tends to be more crowded than the healthy group. Given the clear distribution of values of both groups, we have chosen to perform the automatic binary classification using the SVM method, thinking that it would be the best fit since it generates a plane that tries to separate both groups with the greatest possible difference. This model has achieved a sensitivity of approximately 95%, and there is no difference between its classification and that of the ophthalmologists, being more sensitive than other works with 91% using only images to train the classifier model49. However, the remaining percentage would be influenced by the fact that each eye has been evaluated by a single ophthalmologist, having participated five in total. So, the examination and diagnosis depend on the experience of each one, where it may be the case that incipient cases of FECD can be classified as healthy (False Negatives) or vice versa (False Positives). Also, we should consider that the recovered phase is from all the ocular media and not just from the cornea. Hence, particles, vacuoles, or pigment could be detected by the automatic method as guttae in a healthy eye, labeling it as FECD (False Positives). Moreover, the pathological sample includes patients with cataracts that could hide guttae where both media overlap and the method would classify an eye with FECD as a healthy eye (False Negatives).

On the other hand, this automatic classifier is trained only using metrics, but we are aware that these accuracy results can be improved by applying algorithms that allow training neural networks combining the results of the metrics together with the HPFM, and that is also able to distinguish between different pathologies by taking advantage of the results of previous studies, as other authors have already achieved using only images50. Despite this, the purpose of this study is to present the automatic gutta detection method, its development, and its ability to obtain objective results, leaving in the background the initial proposal of a binary classification model.

In summary, the objective quantification of guttae is an important prerequisite for the Fuchs’ Dystrophy assessment and treatment, and it would allow a better understanding of the origin of visual impairment, as the authors of the references advocate. In the same way, the literature raises the hypothesis that guttae are the source of the aberrations of the posterior corneal surface, although the cause is so far unknown. Because of the presented results, using this high-resolution sensor, it has been possible to process the phase so guttae could be extracted, demonstrating the relevance of analyzing the wavefront phase further low orders, as well as the impact of these protrusions in the ocular aberrations. Alongside this, a reliable objective analysis of the guttae could be performed with an automatic detection method, being a valid tool for the diagnosis and monitoring of the pathology over time.

Acknowledgements

J.G.M.-H. was partially funded by “Research agreement on ophthalmology Wooptix-ULL, 2023”.

Author contributions

C.B.P. contributed to all stages of the project, J.G.M.H. contributed to all stages of the project, M.V.O. contributed to all stages of the project except data acquisition, J.M.T.S. contributed to the software used in the work and the manuscript review, J.M.R.R. contributed to all stages of the project except data acquisition, G.V.R. contributed to the acquisition of the data and the manuscript review, N.A.A. contributed to the acquisition of the data and the manuscript review.

Data availability

The datasets generated during and/or analyzed during the current study are not publicly available due to private company reasons but are available from the corresponding author upon reasonable request.

Competing interests

The authors declare the following competing interests: C.B.P., employment. M.V.O., employment. J.M.T.S., employment. J.M.R.R., funding. J.G.M.H., non-financial competing interests. G.V.R., non-financial competing interests. N.A.A., non-financial competing interests.

Publisher's note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
==== Refs
References

1. Adamis A Filatov V Tripathi B Tripathi R Fuchs endothelial dystrophy of the cornea Surv. Ophthalmol. 1993 38 149 168 10.1016/0039-6257(93)90099-s 8235998
Adamis, A., Filatov, V., Tripathi, B. & Tripathi, R. Fuchs endothelial dystrophy of the cornea. Surv. Ophthalmol. 38, 149–168. 10.1016/0039-6257(93)90099-s (1993).8235998 10.1016/0039-6257(93)90099-s
2. Mustonen R In vivo confocal microscopy of fuchs’ endothelial dystrophy Cornea 1998 17 493 503 10.1097/00003226-199809000-00006 9756443
Mustonen, R. et al. In vivo confocal microscopy of fuchs’ endothelial dystrophy. Cornea 17, 493–503. 10.1097/00003226-199809000-00006 (1998).9756443 10.1097/00003226-199809000-00006
3. Elhalis H Azizi B Jurkunas U Fuchs endothelial corneal dystrophy Ocul. Surf. 2010 8 173 184 10.1016/s1542-0124(12)70232-x 20964980
Elhalis, H., Azizi, B. & Jurkunas, U. Fuchs endothelial corneal dystrophy. Ocul. Surf. 8, 173–184. 10.1016/s1542-0124(12)70232-x (2010).20964980 10.1016/s1542-0124(12)70232-x
4. Fuchs E Dystrophia epithelialis corneae Albrecht von Graefes Archiv fur klinische und experimentelle ophthalmologie 1910 76 478 508 10.1007/BF01986362
Fuchs, E. Dystrophia epithelialis corneae. Albrecht von Graefes Archiv fur klinische und experimentelle ophthalmologie 76, 478–508. 10.1007/BF01986362 (1910).10.1007/BF01986362
5. Kirby D Excrescences of the central area of descemet’s membrane Arch. Ophthalmol. 1925 54 588 91
Kirby, D. Excrescences of the central area of descemet’s membrane. Arch. Ophthalmol. 54, 588–91 (1925).
6. Gifford S Epithelial dystrophy of the cornea and its relation to endothelial dystrophy Am. J. Ophthalmol. 1926 9 81 85 10.1016/S0002-9394(26)90951-3
Gifford, S. Epithelial dystrophy of the cornea and its relation to endothelial dystrophy. Am. J. Ophthalmol. 9, 81–85. 10.1016/S0002-9394(26)90951-3 (1926).10.1016/S0002-9394(26)90951-3
7. Ali M Cho K Srikumaran D Fuchs dystrophy and cataract: Diagnosis, evaluation and treatment Ophthalmol. Ther. 2023 12 691 704 10.1007/s40123-022-00637-1 36637659
Ali, M., Cho, K. & Srikumaran, D. Fuchs dystrophy and cataract: Diagnosis, evaluation and treatment. Ophthalmol. Ther. 12, 691–704. 10.1007/s40123-022-00637-1 (2023).36637659 10.1007/s40123-022-00637-1
8. Vedana G Villarreal J Jun A Fuchs endothelial corneal dystrophy: Current perspectives Clin. Ophthalmol. 2016 10 321 330 10.2147/OPTH.S83467 26937169
Vedana, G., Villarreal, J. & Jun, A. Fuchs endothelial corneal dystrophy: Current perspectives. Clin. Ophthalmol. 10, 321–330. 10.2147/OPTH.S83467 (2016).26937169 10.2147/OPTH.S83467
9. Zhang J Patel D The pathophysiology of fuchs’ endothelial dystrophy - a review of molecular and cellular insights Exp. Eye Res. 2015 130 97 105 10.1016/j.exer.2014.10.023 25446318
Zhang, J. & Patel, D. The pathophysiology of fuchs’ endothelial dystrophy - a review of molecular and cellular insights. Exp. Eye Res. 130, 97–105. 10.1016/j.exer.2014.10.023 (2015).25446318 10.1016/j.exer.2014.10.023
10. Kobashi H Kamiya K Shimizu K Factors influencing visual acuity in fuchs’ endothelial corneal dystrophy OPTOMETRY VISION SCI. 2018 95 21 26 10.1097/OPX.0000000000001157
Kobashi, H., Kamiya, K. & Shimizu, K. Factors influencing visual acuity in fuchs’ endothelial corneal dystrophy. OPTOMETRY VISION SCI. 95, 21–26. 10.1097/OPX.0000000000001157 (2018).10.1097/OPX.0000000000001157
11. Krachmer J Corneal endothelial dystrophy Arch. Ophthalmol. 1978 96 2036 10.1001/archopht.1978.03910060424004 309758
Krachmer, J. Corneal endothelial dystrophy. Arch. Ophthalmol. 96, 2036. 10.1001/archopht.1978.03910060424004 (1978).309758 10.1001/archopht.1978.03910060424004
12. Repp D Hodge D Baratz K McLaren J Patel S Fuchs’ endothelial corneal dystrophy subjective grading versus objective grading based on the central-to-peripheral thickness ratio Ophthalmology 2013 120 687 694 10.1016/j.ophtha.2012.09.022 23369486
Repp, D., Hodge, D., Baratz, K., McLaren, J. & Patel, S. Fuchs’ endothelial corneal dystrophy subjective grading versus objective grading based on the central-to-peripheral thickness ratio. Ophthalmology 120, 687–694. 10.1016/j.ophtha.2012.09.022 (2013).23369486 10.1016/j.ophtha.2012.09.022
13. Matthaei, M. et al. Fuchs endothelial corneal dystrophy: Clinical, genetic, pathophysiologic, and therapeutic aspects. In Movshon, J. & Wandell, B. (eds.) Annu. Rev. Vis. Sci., vol. 5 of Annual Review of Vision Science, 151–175, 10.1146/annurev-vision-091718-014852 (Annual review of vision science, 2019).
14. Zoega G Prevalence and risk factors for cornea guttata in the reykjavik eye study Ophthalmology 2006 113 565 569 10.1016/j.ophtha.2005.12.014 16581419
Zoega, G. et al. Prevalence and risk factors for cornea guttata in the reykjavik eye study. Ophthalmology 113, 565–569. 10.1016/j.ophtha.2005.12.014 (2006).16581419 10.1016/j.ophtha.2005.12.014
15. Eghrari A Retroillumination photography analysis enhances clinical definition of severe fuchs corneal dystrophy Cornea 2015 34 1623 10.1097/ICO.0000000000000656 26488628
Eghrari, A. et al. Retroillumination photography analysis enhances clinical definition of severe fuchs corneal dystrophy. Cornea 34, 1623. 10.1097/ICO.0000000000000656 (2015).26488628 10.1097/ICO.0000000000000656
16. Giasson CJ Morphometry of cells and guttae in subjects with normal or guttate endothelium with a contour detection algorithm Eye Contact Lens 2005 31 158 165 10.1097/01.icl.0000165286.05080.23 16021003
Giasson, C. J. et al. Morphometry of cells and guttae in subjects with normal or guttate endothelium with a contour detection algorithm. Eye Contact Lens 31, 158–165. 10.1097/01.icl.0000165286.05080.23 (2005).16021003 10.1097/01.icl.0000165286.05080.23
17. Watanabe S Relationship between corneal guttae and quality of vision in patients with mild fuchs’ endothelial corneal dystrophy Ophthalmology 2015 122 2103 2109 10.1016/j.ophtha.2015.06.019 26189189
Watanabe, S. et al. Relationship between corneal guttae and quality of vision in patients with mild fuchs’ endothelial corneal dystrophy. Ophthalmology 122, 2103–2109. 10.1016/j.ophtha.2015.06.019 (2015).26189189 10.1016/j.ophtha.2015.06.019
18. Soh Y Peh G Naso S Kocaba V Mehta J Automated clinical assessment of corneal guttae in fuchs endothelial corneal dystrophy Am. J. Ophthalmol. 2021 221 260 272 10.1016/j.ajo.2020.07.029 32730910
Soh, Y., Peh, G., Naso, S., Kocaba, V. & Mehta, J. Automated clinical assessment of corneal guttae in fuchs endothelial corneal dystrophy. Am. J. Ophthalmol. 221, 260–272. 10.1016/j.ajo.2020.07.029 (2021).32730910 10.1016/j.ajo.2020.07.029
19. Hara M Morishige N Chikama T Nishida T Comparison of confocal biomicroscopy and noncontact specular microscopy for evaluation of the corneal endothelium Cornea 2003 22 512 515 10.1097/00003226-200308000-00005 12883342
Hara, M., Morishige, N., Chikama, T. & Nishida, T. Comparison of confocal biomicroscopy and noncontact specular microscopy for evaluation of the corneal endothelium. Cornea 22, 512–515. 10.1097/00003226-200308000-00005 (2003).12883342 10.1097/00003226-200308000-00005
20. Fujimoto H Quantitative regional differences in corneal endothelial abnormalities in the central and peripheral zones in fuchs’ endothelial corneal dystrophy Investig. Ophthalmol. Vis. Sci. 2014 55 5090 5098 10.1167/iovs.14-14249 25061116
Fujimoto, H. et al. Quantitative regional differences in corneal endothelial abnormalities in the central and peripheral zones in fuchs’ endothelial corneal dystrophy. Investig. Ophthalmol. Vis. Sci. 55, 5090–5098. 10.1167/iovs.14-14249 (2014).25061116 10.1167/iovs.14-14249
21. Patel S Amin S Baratz K Kane K McLaren J Anterior corneal aberrations in relation to severity of fuchs endothelial dystrophy Investig. Ophthalmol. Vis. Sci. 2014 55 3580 3580
Patel, S., Amin, S., Baratz, K., Kane, K. & McLaren, J. Anterior corneal aberrations in relation to severity of fuchs endothelial dystrophy. Investig. Ophthalmol. Vis. Sci. 55, 3580–3580 (2014).
22. Arora R Analysis of corneal scheimpflug densitometry and ocular wavefront aberrations post descemet stripping automated endothelial keratoplasty Eye Contact Lens 2022 48 242 249 10.1097/ICL.0000000000000868 34860724
Arora, R. et al. Analysis of corneal scheimpflug densitometry and ocular wavefront aberrations post descemet stripping automated endothelial keratoplasty. Eye Contact Lens 48, 242–249. 10.1097/ICL.0000000000000868 (2022).34860724 10.1097/ICL.0000000000000868
23. van Dijk K Optical quality of the cornea after descemet membrane endothelial keratoplasty Am. J. Ophthalmol. 2014 158 71 79 10.1016/j.ajo.2014.04.008 24784873
van Dijk, K. et al. Optical quality of the cornea after descemet membrane endothelial keratoplasty. Am. J. Ophthalmol. 158, 71–79. 10.1016/j.ajo.2014.04.008 (2014).24784873 10.1016/j.ajo.2014.04.008
24. Waldrop W Gillings M Robertson D Petroll W Mootha V Lower corneal haze and aberrations in descemet membrane endothelial keratoplasty versus descemet stripping automated endothelial keratoplasty in fellow eyes for fuchs endothelial corneal dystrophy Cornea 2020 39 1227 1234 10.1097/ICO.0000000000002416 32639312
Waldrop, W., Gillings, M., Robertson, D., Petroll, W. & Mootha, V. Lower corneal haze and aberrations in descemet membrane endothelial keratoplasty versus descemet stripping automated endothelial keratoplasty in fellow eyes for fuchs endothelial corneal dystrophy. Cornea 39, 1227–1234. 10.1097/ICO.0000000000002416 (2020).32639312 10.1097/ICO.0000000000002416
25. Wacker K McLaren J Amin S Baratz K Patel S Corneal high-order aberrations and backscatter in fuchs’ endothelial corneal dystrophy Ophthalmology 2015 122 1645 1652 10.1016/j.ophtha.2015.05.005 26050543
Wacker, K., McLaren, J., Amin, S., Baratz, K. & Patel, S. Corneal high-order aberrations and backscatter in fuchs’ endothelial corneal dystrophy. Ophthalmology 122, 1645–1652. 10.1016/j.ophtha.2015.05.005 (2015).26050543 10.1016/j.ophtha.2015.05.005
26. Velarde-Rodriguez G Ultra-high resolution optical aberrometry in patients with keratoconus: A cross-sectional study Ophthalmol. Ther. 2023 12 1569 1582 10.1007/s40123-023-00684-2 36856979
Velarde-Rodriguez, G. et al. Ultra-high resolution optical aberrometry in patients with keratoconus: A cross-sectional study. Ophthalmol. Ther. 12, 1569–1582. 10.1007/s40123-023-00684-2 (2023).36856979 10.1007/s40123-023-00684-2
27. Belda-Para C Comparing the clinical applicability of wavefront phase imaging in keratoconus versus normal eyes Sci. Rep. 2024 14 9984 10.1038/s41598-024-60842-9 38693352
Belda-Para, C. et al. Comparing the clinical applicability of wavefront phase imaging in keratoconus versus normal eyes. Sci. Rep. 14, 9984. 10.1038/s41598-024-60842-9 (2024).38693352 10.1038/s41598-024-60842-9
28. Bonaque-González S The optics of the human eye at 8.6 μm resolution Sci. Rep. 2021 11 23334 10.1038/s41598-021-02653-w 34857820
Bonaque-González, S. et al. The optics of the human eye at 8.6 m resolution. Sci. Rep. 11, 23334. 10.1038/s41598-021-02653-w (2021).34857820 10.1038/s41598-021-02653-w
29. Trujillo-Sevilla, J. M., Casanova Gonzalez, O., Bonaque-González, S., Gaudestad, J. & Rodríguez Ramos, J. M. High-resolution wave front phase sensor for silicon wafer metrology. In Photonic Instrumentation Engineering VI, vol. 10925, 109250I, 10.1117/12.2505764. International Society for Optics and Photonics (SPIE, 2019).
30. Trujillo-Sevilla J Velasco-Ocaña M Bonaque-González S Belda-Para C Rodríguez-Ramos J Wavefront phase measurement of striae in optical glass Appl. Opt. 2022 61 3912 3918 10.1364/AO.450219 36256060
Trujillo-Sevilla, J., Velasco-Ocaña, M., Bonaque-González, S., Belda-Para, C. & Rodríguez-Ramos, J. Wavefront phase measurement of striae in optical glass. Appl. Opt. 61, 3912–3918. 10.1364/AO.450219 (2022).36256060 10.1364/AO.450219
31. Cairós C Refractive index estimation in biological tissues by quantitative phase imaging Opt. Mater. 2023 142 114087 10.1016/j.optmat.2023.114087
Cairós, C. et al. Refractive index estimation in biological tissues by quantitative phase imaging. Opt. Mater. 142, 114087. 10.1016/j.optmat.2023.114087 (2023).10.1016/j.optmat.2023.114087
32. Bonaque-González S New high-resolution wavefront sensing ophthalmic technique for the characterization of ocular optics Ophthalmic Technol. XXX 2020 11218 117 123 10.1117/12.2548595
Bonaque-González, S. et al. New high-resolution wavefront sensing ophthalmic technique for the characterization of ocular optics. Ophthalmic Technol. XXX 11218, 117–123. 10.1117/12.2548595 (2020).10.1117/12.2548595
33. Bonaque-González S Extremely high resolution ocular aberrometry up to 2.4 million points Investig. Ophthalmol. Vis. Sci. 2019 60 603 603
Bonaque-González, S. et al. Extremely high resolution ocular aberrometry up to 2.4 million points. Investig. Ophthalmol. Vis. Sci. 60, 603–603 (2019).
34. Verma S Hesser J Arba-Mosquera S Optimum laser beam characteristics for achieving smoother ablations in laser vision correction Investig. Ophthalmol. Vis. Sci. 2017 58 2021 2037 10.1167/iovs.16-21025 28384723
Verma, S., Hesser, J. & Arba-Mosquera, S. Optimum laser beam characteristics for achieving smoother ablations in laser vision correction. Investig. Ophthalmol. Vis. Sci. 58, 2021–2037. 10.1167/iovs.16-21025 (2017).28384723 10.1167/iovs.16-21025
35. Yang X Cui W A novel spatial clustering algorithm based on delaunay triangulation J. Softw. Eng. Appl. 2010 03 141 149 10.4236/jsea.2010.32018
Yang, X. & Cui, W. A novel spatial clustering algorithm based on delaunay triangulation. J. Softw. Eng. Appl. 03, 141–149. 10.4236/jsea.2010.32018 (2010).10.4236/jsea.2010.32018
36. Greenacre M Principal component analysis Nat. Rev. Methods Primers 2022 2 100 10.1038/s43586-022-00184-w
Greenacre, M. et al. Principal component analysis. Nat. Rev. Methods Primers 2, 100. 10.1038/s43586-022-00184-w (2022).10.1038/s43586-022-00184-w
37. MathWorks. Principal component analysis of raw data. https://es.mathworks.com/help/stats/pca.html.
38. MathWorks. Support vector machines for binary classification. https://mathworks.com/help/stats/support-vector-machines-for-binary-classification.html.
39. Cristianini N Shawe-Taylor J An Introduction to Support Vector Machines and Other Kernel-based Learning Methods 2000 Cambridge Cambridge University Press
Cristianini, N. & Shawe-Taylor, J. An Introduction to Support Vector Machines and Other Kernel-based Learning Methods (Cambridge University Press, Cambridge, 2000).
40. McNemar Q Note on the sampling error of the difference between correlated proportions or percentages Psychometrika 1947 12 153 157 10.1007/bf02295996 20254758
McNemar, Q. Note on the sampling error of the difference between correlated proportions or percentages. Psychometrika 12, 153–157. 10.1007/bf02295996 (1947).20254758 10.1007/bf02295996
41. Dietterich T Approximate statistical tests for comparing supervised classification learning algorithms Neural Comput. 1998 10 1895 1923 10.1162/089976698300017197 9744903
Dietterich, T. Approximate statistical tests for comparing supervised classification learning algorithms. Neural Comput. 10, 1895–1923. 10.1162/089976698300017197 (1998).9744903 10.1162/089976698300017197
42. Rainio O Teuho J Klén R Evaluation metrics and statistical tests for machine learning Sci. Rep. 2024 10.1038/s41598-024-56706-x 38977765
Rainio, O., Teuho, J. & Klén, R. Evaluation metrics and statistical tests for machine learning. Sci. Rep.[SPACE]10.1038/s41598-024-56706-x (2024).38977765 10.1038/s41598-024-56706-x
43. Rosenblum P Stark W Maumenee I Hirst L Maumenee A Hereditary fuchs’ dystrophy Am. J. Ophthalmol. 1980 90 455 462 10.1016/s0002-9394(14)75011-1 6968504
Rosenblum, P., Stark, W., Maumenee, I., Hirst, L. & Maumenee, A. Hereditary fuchs’ dystrophy. Am. J. Ophthalmol. 90, 455–462. 10.1016/s0002-9394(14)75011-1 (1980).6968504 10.1016/s0002-9394(14)75011-1
44. Dong W Zou L Pan Z Wang L Morphologic characteristics of cornea in fuchs endothelial dystrophy observed by confocal microscopy [Zhonghua yan ke za zhi] Chin. J. Ophthalmol. 2004 40 465 470
Dong, W., Zou, L., Pan, Z. & Wang, L. Morphologic characteristics of cornea in fuchs endothelial dystrophy observed by confocal microscopy. [Zhonghua yan ke za zhi] Chin. J. Ophthalmol. 40, 465–470 (2004).
45. Vaitinadapoulé H Nanotopography by chromatic confocal microscopy of the endothelium in fuchs endothelial corneal dystrophy, pseudophakic bullous keratopathy and healthy corneas Br. J. Ophthalmol. 2023 10.1136/bjo-2023-323297
Vaitinadapoulé, H. et al. Nanotopography by chromatic confocal microscopy of the endothelium in fuchs endothelial corneal dystrophy, pseudophakic bullous keratopathy and healthy corneas. Br. J. Ophthalmol.[SPACE]10.1136/bjo-2023-323297 (2023).10.1136/bjo-2023-323297
46. Iovino C Corneal endothelium features in fuchs’ endothelial corneal dystrophy: A preliminary 3d anterior segment optical coherence tomography study PLoS ONE 2018 13 1 12 10.1371/journal.pone.0207891
Iovino, C. et al. Corneal endothelium features in fuchs’ endothelial corneal dystrophy: A preliminary 3d anterior segment optical coherence tomography study. PLoS ONE 13, 1–12. 10.1371/journal.pone.0207891 (2018).10.1371/journal.pone.0207891
47. Shah K Eghrari A Vanner E O’Brien T Koo E Scheimpflug corneal densitometry values and severity of guttae in relation to visual acuity in fuchs endothelial corneal dystrophy Cornea 2022 41 692 698 10.1097/ico.0000000000002762 35175018
Shah, K., Eghrari, A., Vanner, E., O’Brien, T. & Koo, E. Scheimpflug corneal densitometry values and severity of guttae in relation to visual acuity in fuchs endothelial corneal dystrophy. Cornea 41, 692–698. 10.1097/ico.0000000000002762 (2022).35175018 10.1097/ico.0000000000002762
48. van Dam M Lane R Extended analysis of curvature sensing JOSA A 2002 19 1390 1397 10.1364/JOSAA.19.001390 12095207
van Dam, M. & Lane, R. Extended analysis of curvature sensing. JOSA A 19, 1390–1397. 10.1364/JOSAA.19.001390 (2002).12095207 10.1364/JOSAA.19.001390
49. Eleiwa T Elsawy A Ozcan E Abou Shousha M Automated diagnosis and staging of fuchs’ endothelial cell corneal dystrophy using deep learning Eye and Vision 2020 7 44 10.1186/s40662-020-00209-z 32884962
Eleiwa, T., Elsawy, A., Ozcan, E. & Abou Shousha, M. Automated diagnosis and staging of fuchs’ endothelial cell corneal dystrophy using deep learning. Eye and Vision 7, 44. 10.1186/s40662-020-00209-z (2020).32884962 10.1186/s40662-020-00209-z
50. Elsawy A Multidisease deep learning neural network for the diagnosis of corneal diseases Am. J. Ophthalmol. 2021 226 252 261 10.1016/j.ajo.2021.01.018 33529589
Elsawy, A. et al. Multidisease deep learning neural network for the diagnosis of corneal diseases. Am. J. Ophthalmol. 226, 252–261. 10.1016/j.ajo.2021.01.018 (2021).33529589 10.1016/j.ajo.2021.01.018
