Open-access operating algorithms for commercial videokeratographer and improvement of corneal sampling

We present an algorithm to process images of reflected Placido rings captured by a commercial videokeratoscope. Raw data are obtained with no Cartesian-to-polar-coordinate conversion, thus avoiding interpolation and associated numerical artifacts. The method provides a characteristic equation for the device and is able to process around 6 times more corneal data than the commercial software. Our proposal allows complete control over the whole process from the capture of corneal images until the computation of curvature radii.

In vivo characterization of corneal biomechanics

Interest in corneal biomechanics has increased with the development of new refractive surgery techniques aimed at modifying corneal properties and a variety of surgical options for corneal ectasia management. The human cornea behaves as soft biological material. It is a viscoelastic tissue and its response to a force applied to it depends not only on the magnitude of the force, but also on the velocity of the application. There are concerns about the limitations to measuring corneal biomechanical properties in vivo. To date, 2 systems are available for clinical use: the Ocular Response Analyzer, a dynamic bidirectional applanation device, and the Corvis ST, a dynamic Scheimpflug analyzer device. These devices are useful in clinical practice, especially for planning some surgical procedures and earlier detection of ectatic conditions, but further research is needed to connect the clinical measurements obtained with these devices to the standard mechanical properties.

Comparative analysis of the relationship between anterior and posterior corneal shape analyzed by Scheimpflug photography in normal and keratoconus eyes

Background To analyze and compare the relationship between anterior and posterior corneal shape evaluated by a tomographic system combining the Scheimpflug photography and Placido-disc in keratoconus and normal healthy eyes, as well as to evaluate its potential diagnostic value. Methods Comparative case series including a sample of 161 eyes of 161 subjects with ages ranging from 7 to 66 years and divided into two groups: normal group including 100 healthy eyes of 100 subjects, and keratoconus group including 61 keratoconus eyes of 61 patients. All eyes received a comprehensive ophthalmologic examination including an anterior segment analysis with the Sirius system (CSO). Antero-posterior ratios for corneal curvature (k ratio) and shape factor (p ratio) were calculated. Logistic regression analysis was used to evaluate if some antero–posterior ratios combined with other clinical parameters were predictors of the presence of keratoconus. Results No statistically significant differences between groups were found in the antero–posterior k ratios for 3-, 5- and 7-mm diameter corneal areas (p ≥ 0.09). The antero–posterior p ratio for 4.5- and 8-mm diameter corneal areas was significantly higher in the normal group than in the keratoconus group (p < 0.01). The k ratio for 3, 5, and 7 mm was significantly higher in the keratoconus grade IV subgroup than in the normal group (p < 0.01). Furthermore, significant differences were found in the p ratio between the normal group and the keratoconus grade II subgroup (p ≤ 0.01). Finally, the logistic regression analysis identified as significant independent predictors of the presence of keratoconus (p < 0.01) the 8-mm anterior shape factor, the anterior chamber depth, and the minimal corneal thickness. Conclusions The antero-posterior k and p ratios are parameters with poor prediction ability for keratoconus, in spite of the trend to the presence of more prolate posterior corneal surfaces compared to the anterior in keratoconus eyes.

Deep Anterior Lamellar Keratoplasty versus Penetrating Keratoplasty for Macular Corneal Dystrophy: A Randomized Trial

Purpose: To compare outcomes of big-bubble deep anterior lamellar keratoplasty (DALK) and penetrating keratoplasty (PK) for macular corneal dystrophy. Design: Prospective, randomized, interventional case series. Methods: Setting: Single hospital. Patients: Eighty-two eyes of 54 patients requiring keratoplasty for the treatment of macular corneal dystrophy without endothelial involvement were included. Main outcome measures: Operative complications, uncorrected visual acuity, best-corrected visual acuity, contrast sensitivity function, higher-order aberrations, and endothelial cell density were evaluated. Results: The DALK and PK group consisted of 35 and 41 eyes, respectively. Best-corrected visual acuity after surgery was 20/40 or better 68.5% and 70.7% of the eyes in the DALK and PK groups, respectively (P > .05). No statistically significant differences between groups were found in contrast sensitivity function with and without glare for any spatial frequency (P > .05). Significantly higher levels of higher-order aberrations were found in the DALK group (P < .01). In both groups, a progressive and statistically significant reduction in endothelial cell density was found (P < .01). At the last follow-up, the mean endothelial cell loss was 18.1% and 26.9% in DALK and PK groups, respectively (P = .03). Graft rejection episodes were seen in 5 eyes (12.1%) in the PK group, and regrafting was necessary in 3 eyes (7.3%). Recurrence of the disease was documented in 5.7% and 4.8% of the eyes in the DALK and PK groups, respectively. Conclusions: Deep anterior lamellar keratoplasty with the big-bubble technique provided comparable visual and optical results as PK and resulted in less endothelial damage, as well as eliminating endothelial rejection in macular corneal dystrophy. Deep anterior lamellar keratoplasty surgery is a viable option for macular corneal dystrophy without endothelial involvement.

Screening Subclinical Keratoconus With Placido-Based Corneal Indices

Purpose. To assess in a sample of normal, keratoconic, and keratoconus (KC) suspect eyes the performance of a set of new topographic indices computed directly from the digitized images of the Placido rings. Methods. This comparative study was composed of a total of 124 eyes of 106 patients from the ophthalmic clinics Vissum Alicante and Vissum Almería (Spain) divided into three groups: control group (50 eyes), KC group (50 eyes), and KC suspect group (24 eyes). In all cases, a comprehensive examination was performed, including the corneal topography with a Placidobased CSO topography system. Clinical outcomes were compared among groups, along with the discriminating performance of the proposed irregularity indices. Results. Significant differences at level 0.05 were found on the values of the indices among groups by means of Mann-Whitney-Wilcoxon nonparametric test and Fisher exact test. Additional statistical methods, such as receiver operating characteristic analysis and K-fold cross validation, confirmed the capability of the indices to discriminate between the three groups. Conclusions. Direct analysis of the digitized images of the Placido mires projected on the cornea is a valid and effective tool for detection of corneal irregularities. Although based only on the data from the anterior surface of the cornea, the new indices performed well even when applied to the KC suspect eyes. They have the advantage of simplicity of calculation combined with high sensitivity in corneal irregularity detection and thus can be used as supplementary criteria for diagnosing and grading KC that can be added to the current keratometric classifications.

Technologies for Anatomical and Geometric Characterization of the Corneal Structure and Anterior Segment: A Review

Corneal and anterior segment imaging techniques have become a crucial tool in the clinical practice of ophthalmology, with a great variety of applications, such as corneal curvature and pachymetric analysis, detection of ectatic corneal conditions, anatomical study of the anterior segment prior to phakic intraocular lens implantation, or densitometric analysis of the crystalline lens. From the Placido-based systems that allow only a characterization of the geometry of the anterior corneal surface to the Scheimpflug photography-based systems that provide a characterization of the cornea, anterior chamber, and crystalline lens, there is a great variety of devices with the capability of analyzing different anatomical parameters with very high precision. To date, Scheimpflug photography-based systems are the devices providing the more complete analysis of the anterior segment in a non-invasive way. More developments are required in anterior segment imaging technologies in order to improve the analysis of the crystalline lens structure as well as the ocular structures behind the iris in a non-invasive way when the pupil is not dilated.

Repeatability and reproducibility of corneal thickness using SOCT Copernicus HR

Purpose: The aim of this study is to determine the reliability of corneal thickness measurements derived from SOCT Copernicus HR (Fourier domain OCT). Methods: Thirty healthy eyes of 30 subjects were evaluated. One eye of each patient was chosen randomly. Images were obtained of the central (up to 2.0 mm from the corneal apex) and paracentral (2.0 to 4.0 mm) cornea. We assessed corneal thickness (central and paracentral) and epithelium thickness. The intra-observer repeatability data were analysed using the intra-class correlation coefficient (ICC) for a range of 95 per cent within-subject standard deviation (SW) and the within-subject coefficient of variation (CW). The level of agreement by Bland–Altman analysis was also represented for the study of the reproducibility between observers and agreement between methods of measurement (automatic versus manual). Results: The mean value of the central corneal thickness (CCT) was 542.4 ± 30.1 μm (SD). There was a high intra-observer agreement, finding the best result in the central sector with an intra-class correlation coefficient of 0.99, 95 per cent CI (0.989 to 0.997) and the worst, in the minimum corneal thickness, with an intra-class correlation coefficient of 0.672, 95 per cent CI (0.417 to 0.829). Reproducibility between observers was very high. The best result was found in the central sector thickness obtained both manually and automatically with an intra-class correlation coefficient of 0.990 in both cases and the worst result in the maximum corneal thickness with an intra-class correlation coefficient of 0.827. The agreement between measurement methods was also very high with intra-class correlation coefficient greater than 0.91. On the other hand the repeatability and reproducibility for epithelial measurements was poor. Conclusion: Pachymetric mapping with SOCT Copernicus HR was found to be highly repeatable and reproducible. We found that the device lacks an appropriate ergonomic design as proper focusing of the laser beam onto the cornea for anterior segment scanning required that patients were positioned slightly farther away from the machine head-rest than in the setup for retinal imaging.

Algorithm for Correcting the Keratometric Error in the Estimation of the Corneal Power in Eyes With Previous Myopic Laser Refractive Surgery

Purpose: To calculate theoretically the errors in the estimation of corneal power when using the keratometric index (nk) in eyes that underwent laser refractive surgery for the correction of myopia and to define and validate clinically an algorithm for minimizing such errors. Methods: Differences between corneal power estimation by using the classical nk and by using the Gaussian equation in eyes that underwent laser myopic refractive surgery were simulated and evaluated theoretically. Additionally, an adjusted keratometric index (nkadj) model dependent on r1c was developed for minimizing these differences. The model was validated clinically by retrospectively using the data from 32 myopic eyes [range, −1.00 to −6.00 diopters (D)] that had undergone laser in situ keratomileusis using a solid-state laser platform. The agreement between Gaussian (PGaussc) and adjusted keratometric (Pkadj) corneal powers in such eyes was evaluated. Results: It was found that overestimations of corneal power up to 3.5 D were possible for nk = 1.3375 according to our simulations. The nk value to avoid the keratometric error ranged between 1.2984 and 1.3297. The following nkadj models were obtained: nkadj= −0.0064286r1c + 1.37688 (Gullstrand eye model) and nkadj = −0.0063804r1c + 1.37806 (Le Grand). The mean difference between Pkadj and PGaussc was 0.00 D, with limits of agreement of −0.45 and +0.46 D. This difference correlated significantly with the posterior corneal radius (r = −0.94, P < 0.01). Conclusions: The use of a single nk for estimating the corneal power in eyes that underwent a laser myopic refractive surgery can lead to significant errors. These errors can be minimized by using a variable nk dependent on r1c.

Correlation of the Corneal Toricity Between Anterior and Posterior Corneal Surfaces in the Normal Human Eye

Purpose: To evaluate the correlation of the magnitude of corneal toricity and power vector components of both corneal surfaces measured with a Scheimpflug photography-based system. Methods: A total of 117 healthy normal eyes of 117 subjects selected randomly with ages ranging from 7 to 80 years were included. All eyes received an anterior segment and corneal analysis with the Sirius system (CSO) evaluating the anterior and posterior mean toricity for 3 and 7 mm (aAST and pAST). The vector components J0 and J45 as well as the overall strength blur (B) were calculated for each keratometric measurement using the procedure defined by Thibos and Horner. Results: The coefficient of correlation between aAST and pAST was 0.52 and 0.62 and the mean anteroposterior ratio for toricity was 0.46 ± 0.39 and 0.57 ± 0.75 for 3 and 7 mm, respectively. These ratios correlated significantly with aAST, anterior corneal J0, and manifest refraction J0 (r ≥ 0.39, P < 0.01). The coefficient of correlation was 0.69 and 0.81 between anterior and posterior J0 for 3 and 7 mm, respectively. For J45, the coefficients were 0.62 and 0.71, respectively. The linear regression analysis revealed that the pAST and power vectors could be predicted from the anterior corneal data (R2 ≥ 0.40, P < 0.01). Conclusions: The toricity and astigmatic power vector components of the posterior corneal surface in the human healthy eye are related to those of the anterior and therefore can be predicted consistently from the anterior toricity and astigmatic power vectors.

Correlation of Anterior and Posterior Corneal Shape in Keratoconus

Purpose: To evaluate the correlation of the mean curvature and shape factors of both corneal surfaces for different corneal diameters measured with the Scheimpflug photography–based system in keratoconus eyes. Methods: A total of 61 keratoconus eyes of 61 subjects, aged 14 to 64 years, were included in this study. All eyes received a comprehensive ophthalmologic examination including anterior segment and corneal analysis with the Sirius system (CSO): anterior and posterior mean corneal radius for 3, 5, and 7 mm (aKM, pKM), anterior and posterior mean shape factor for 4.5 and 8 mm (ap, pp), central and minimal corneal thickness, and anterior chamber depth. Results: Mean aKM/pKM ratio around 1.20 (range, 0.95–1.48) was found for all corneal diameters (P = 0.24). Weak but significant correlations of this ratio with pachymetric parameters were found (r between −0.28 and −0.34, P < 0.04). The correlation coefficient between aKM and pKM was ≥0.92 for all corneal diameters. A strong and significant correlation was also found between ap and pp (r ≥ 0.86, P < 0.01). The multiple regression analysis revealed that central pKM was significantly correlated with aKM, central corneal thickness, anterior chamber depth, and spherical equivalent (R2 ≥ 0.88, P < 0.01) and that 8 mm pp was significantly correlated with 8 mm ap and age (R2 = 0.89, P < 0.01). Conclusions: Central posterior corneal curvature and shape factor in the keratoconus eye can be consistently predicted from the anterior corneal curvature and shape factor, respectively, in combination with other anatomical and ocular parameters.

Corneal biomechanics: a review

Biomechanics is often defined as ‘mechanics applied to biology’. Due to the variety and complexity of the behaviour of biological structures and materials, biomechanics is better defined as the development, extension and application of mechanics for a better understanding of physiology and physiopathology and consequently for a better diagnosis and treatment of disease and injury. Different methods for the characterisation of corneal biomechanics are reviewed in detail, including those that are currently commercially available (Ocular Response Analyzer and CorVis ST). The clinical applicability of the parameters provided by these devices are discussed, especially in the fields of glaucoma, detection of ectatic disorders and orthokeratology. Likewise, other methods are also reviewed, such as Brillouin microscopy or dynamic optical coherence tomography and others with potential application to clinical practice but not validated for in vivo measurements, such as ultrasonic elastography. Advantages and disadvantages of all these techniques are described. Finally, the concept of biomechanical modelling is revised as well as the requirements for developing biomechanical models, with special emphasis on finite element modelling.

Corneal Biomechanics, Retinal Nerve Fiber Layer, and Optic Disc in Children

Purpose: To evaluate the possible associations between corneal biomechanical parameters, optic disc morphology, and retinal nerve fiber layer (RNFL) thickness in healthy white Spanish children. Methods: This cross-sectional study included 100 myopic children and 99 emmetropic children as a control group, ranging in age from 6 to 17 years. The Ocular Response Analyzer was used to measure corneal hysteresis (CH) and corneal resistance factor. The optic disc morphology and RNFL thickness were assessed using posterior segment optical coherence tomography (Cirrus HD-OCT). The axial length was measured using an IOLMaster, whereas the central corneal thickness was measured by anterior segment optical coherence tomography (Visante OCT). Results: The mean (±SD) age and spherical equivalent were 12.11 (±2.76) years and −3.32 (±2.32) diopters for the myopic group and 11.88 (±2.97) years and +0.34 (±0.41) diopters for the emmetropic group. In a multivariable mixed-model analysis in myopic children, the average RNFL thickness and rim area correlated positively with CH (p = 0.007 and p = 0.001, respectively), whereas the average cup-to-disc area ratio correlated negatively with CH (p = 0.01). We did not observe correlation between RNFL thickness and axial length (p = 0.05). Corneal resistance factor was only positively correlated with the rim area (p = 0.001). The central corneal thickness did not correlate with the optic nerve parameters or with RNFL thickness. These associations were not found in the emmetropic group (p > 0.05 for all). Conclusions: The corneal biomechanics characterized with the Ocular Response Analyzer system are correlated with the optic disc profile and RNFL thickness in myopic children. Low CH values may indicate a reduction in the viscous dampening properties of the cornea and the sclera, especially in myopic children.

Relationship among Corneal Biomechanics, Refractive Error, and Axial Length

Purpose: To evaluate the relationship between different ocular and corneal biomechanical parameters in emmetropic and ametropic healthy white children. Methods: This study included 293 eyes of 293 healthy Spanish children (135 boys and 158 girls), ranging in age from 6 to 17 years. Subjects were divided according to the refractive error: control (emmetropia, 99 children), myopia (100 children), and hyperopia (94 children) groups. In all cases, corneal hysteresis (CH) and corneal resistance factor (CRF) were evaluated with the Ocular Response Analyzer system. Axial length (AL) and mean corneal power were also measured by partial coherence interferometry (IOLMaster), and central corneal thickness (CCT) and anterior chamber depth were measured by anterior segment optical coherence tomography (Visante). Results: Mean (±SD) CH and CRF were 12.12 (±1.71) and 12.30 (±1.89) mm Hg, respectively. Mean (±SD) CCT was 542.68 (±37.20) μm and mean (±SD) spherical equivalent was +0.14 (±3.41) diopters. A positive correlation was found between CH and CRF (p < 0.001), and both correlated as well with CCT (p < 0.0001). Corneal resistance factor was found to decrease with increasing age (p = 0.01). Lower levels of CH were associated with longer AL and more myopia (p < 0.001 and p = 0.001, respectively). Higher values of CH were associated with increasing hyperopia. Significant differences in CH were found between emmetropic and myopic groups (p < 0.001) and between myopic and hyperopic groups (p = 0.011). There were also significant differences in CRF between emmetropic and myopic groups (p = 0.02). Multiple linear regression analysis showed that lower CH and CRF significantly associated with thinner CCT, longer AL, and flatter corneal curvature. Conclusions: The Ocular Response Analyzer corneal biomechanical properties seem to be compromised in myopia from an early age, especially in high myopia.

New Approach for Correction of Error Associated With Keratometric Estimation of Corneal Power in Keratoconus

The aim of this study was to obtain the exact value of the keratometric index (nkexact) and to clinically validate a variable keratometric index (nkadj) that minimizes this error. Methods: The nkexact value was determined by obtaining differences (DPc) between keratometric corneal power (Pk) and Gaussian corneal power (PGauss c ) equal to 0. The nkexact was defined as the value associated with an equivalent difference in the magnitude of DPc for extreme values of posterior corneal radius (r2c) for each anterior corneal radius value (r1c). This nkadj was considered for the calculation of the adjusted corneal power (Pkadj). Values of r1c ∈ (4.2, 8.5) mm and r2c ∈ (3.1, 8.2) mm were considered. Differences of True Net Power with PGauss c , Pkadj, and Pk(1.3375) were calculated in a clinical sample of 44 eyes with keratoconus. Results: nkexact ranged from 1.3153 to 1.3396 and nkadj from 1.3190 to 1.3339 depending on the eye model analyzed. All the nkadj values adjusted perfectly to 8 linear algorithms. Differences between Pkadj and PGauss c did not exceed 60.7 D (Diopter). Clinically, nk = 1.3375 was not valid in any case. Pkadj and True Net Power and Pk(1.3375) and Pkadj were statistically different (P , 0.01), whereas no differences were found between PGauss c and Pkadj (P . 0.01). Conclusions: The use of a single value of nk for the calculation of the total corneal power in keratoconus has been shown to be imprecise, leading to inaccuracies in the detection and classification of this corneal condition. Furthermore, our study shows the relevance of corneal thickness in corneal power calculations in keratoconus.

Estimation of the Central Corneal Power in Keratoconus: Theoretical and Clinical Assessment of the Error of the Keratometric Approach

Purpose: The aim of this study was to analyze theoretically the errors in the central corneal power calculation in eyes with keratoconus when a keratometric index (nk) is used and to clinically confirm the errors induced by this approach. Methods: Differences (DPc) between central corneal power estimation with the classical nk (Pk) and with the Gaussian equation (PGauss c ) in eyes with keratoconus were simulated and evaluated theoretically, considering the potential range of variation of the central radius of curvature of the anterior (r1c) and posterior (r2c) corneal surfaces. Further, these differences were also studied in a clinical sample including 44 keratoconic eyes (27 patients, age range: 14–73 years). The clinical agreement between Pk and PGauss c (true net power) obtained with a Scheimpflug photography–based topographer was evaluated in such eyes. Results: For nk = 1.3375, an overestimation was observed in most cases in the theoretical simulations, with DPc ranging from an underestimation of 20.1 diopters (D) (r1c = 7.9 mm and r2c = 8.2 mm) to an overestimation of 4.3 D (r1c = 4.7 mm and r2c = 3.1 mm). Clinically, Pk always overestimated the PGauss c given by the topography system in a range between 0.5 and 2.5 D (P , 0.01). The mean clinical DPc was 1.48 D, with limits of agreement of 0.71 and 2.25 D. A very strong statistically significant correlation was found between DPc and r2c (r = 20.93, P , 0.01). Conclusions: The use of a single value for nk for the calculation of corneal power is imprecise in keratoconus and can lead to significant clinical errors.