Non-Orthogonal Corneal Astigmatism among Normal and Keratoconic Brazilian and Chinese populations

Minimum Corneal Diameter and Anterior Steep Axis Curvature Share the Same Meridian: A Novel Finding

PURPOSE

To define the external scleral sulcus (ESS) on a Scheimpflug image and use it for a morphometric analysis of corneal diameter (CD).

DESIGN

Retrospective, cross-sectional study of pediatric Asian-Indian eyes.

METHODS

One random eye of 353 subjects between 5 and 18 years underwent 25-scan Pentacam HR imaging. For all scans, densitometry values along the anterior corneal edge were recorded and differentiated. The peaks on the differentiated curve were chosen as the ESS points, and this distance between them was called CD. Vertical (vCD), maximum (maxCD), minimum (minCD) CD and their meridians were defined. Multiple regression models (MRMs) with CD and other Pentacam parameters were built to predict astigmatism and its axis, mean keratometry (Kmean), and Belin/Ambrósio enhanced ectasia display deviation (BAD-D). MRMs were validated using intraclass correlation coefficient (ICC). Estimated horizontal CD (hCD) was validated against digital caliper measurement using ICC.

RESULTS

The ICC (95% CI) between caliper and hCD was 0.96 (0.93, 0.97). MRM predictions (P < .001) used CD parameters, anterior chamber depth, corneal volume and distance from the corneal thinnest location to apex. These predictions achieved an ICC of 0.34 (0.18, 0.46), 0.82 (0.78, 0.86), 0.87 (0.84, 0.89), and 0.81 (0.76, 0.84), respectively. The astigmatism axis prediction depended on the minCD and maxCD meridians. Its within-subject SD (4.97°) was less than 2 consecutive Pentacam scan angles (7.2°).

CONCLUSIONS

The CD metric strongly correlated with the astigmatism axis, keratometry, and BAD-D. Its spatial description may be significant in corneal treatment planning and disease diagnoses.

Colombian Ocular Diseases Epidemiology Study (CODES): incidence and sociodemographic characterisation of keratoconus between 2015 and 2020

OBJECTIVE

To estimate the incidence and describe the demographic characteristics of keratoconus (KC) in Colombia using national health registry data between January 1st 2015 and December 31th 2020.

METHODS AND ANALYSIS

We conducted a nationwide, population-based study using the Integrated Social Protection Information System from the Colombian Ministry of Health, the unique official database in the country. We used the International Classification of Diseases code for KC (H186) to identify the number of new cases of KC and estimate the incidence rates both overall and according to age and sex. We made a standard morbidity ratio map to graph Colombia's morbidity risk of KC onset.

RESULTS

Of 50 372 424 subjects, 21 710 had KC between 2015 and 2020. However, due to the COVID-19 pandemic, all the incidence rates of this study were based on the 18 419 reported until 2019. The incidence rate in the general population was 10.36 (95% CI 10.08 to 10.64) per 100 000 inhabitants. The incidence peak among males was in their early 20s and females in their late 20s. The overall male to female incidence rate ratio was 1.60. Regarding the distribution of the disease, most cases were reported in Bogotá (48.64%), Antioquia (14.04%) and Cundinamarca (10.38%).

CONCLUSION

We performed the first nationwide, population-based study of KC in Latin America, finding distribution patterns similar to those reported in the literature. This study provides valuable information on the epidemiology of KC in Colombia, which is helpful in the development of policies for the diagnosis, prevention and treatment of the disease.

Limitations of Reconstructing Pentacam Rabbit Corneal Tomography by Zernike Polynomials

The study aims to investigate the likelihood of Zernike polynomial being used for reconstructing rabbit corneal surfaces as scanned by the Pentacam segment tomographer, and hence evaluate the accuracy of corneal power maps calculated from such Zernike fitted surfaces. The study utilised a data set of both eyes of 21 rabbits using a reverse engineering approach for deductive reasoning. Pentacam raw elevation data were fitted to Zernike polynomials of orders 2 to 20. The surface fitting process to Zernike polynomials was carried out using randomly selected 80% of the corneal surface data points, and the root means squared fitting error (RMS) was determined for the other 20% of the surface data following the Pareto principle. The process was carried out for both the anterior and posterior surfaces of the corneal surfaces that were measured via Pentacam scans. Raw elevation data and the fitted corneal surfaces were then used to determine corneal axial and tangential curvature maps. For reconstructed surfaces calculated using the Zernike fitted surfaces, the mean and standard deviation of the error incurred by the fitting were calculated. For power maps computed using the raw elevation data, different levels of discrete cosine transform (DCT) smoothing were employed to infer the smoothing level utilised by the Pentacam device. The RMS error was not significantly improved for Zernike polynomial orders above 12 and 10 when fitting the anterior and posterior surfaces of the cornea, respectively. This was noted by the statistically non-significant increase in accuracy when the order was increased beyond these values. The corneal curvature calculations suggest that a smoothing process is employed in the corneal curvature maps outputted by the Pentacam device; however, the exact smoothing method is unknown. Additionally, the results suggest that fitting corneal surfaces to high-order Zernike polynomials will incur a clinical error in the calculation of axial and tangential corneal curvature of at least 0.16 ± 01 D and 0.36 ± 0.02 D, respectively. Rabbit corneal anterior and posterior surfaces scanned via the Pentacam were optimally fitted to orders 12 and 10 Zernike polynomials. This is essential to get stable values of high-order aberrations that are not affected by Zernike polynomial fittings, such as comas for Intracorneal Ring Segments (ICRS) adjustments or spherical aberration for pre-cataract operations. Smoothing was necessary to replicate the corneal curvature maps outputted by the Pentacam tomographer, and fitting corneal surfaces to Zernike polynomials introduces errors in the calculation of both the axial and tangential corneal curvatures.

The Efficiency of Using Mirror Imaged Topography in Fellow Eyes Analyses of Pentacam HR Data

Purpose: To investigate the effectiveness of flipping left corneas topography and analysethem quantitively along with fellow right corneas based on the assumption that they are mirror images of each other. Methods: The study involved scanning both eyes of 177 healthy participants (aged 35.3 ± 15.8) and 75 keratoconic participants (aged 33.9 ± 17.8). Clinical tomography data were collected for both eyes using the Pentacam HR and processed by a fully automated custom-built MATLAB code. For every case, the right eye was used as a datum fixed surface while the left cornea was flipped around in the superior–inferior direction. In this position, the root-mean-squared difference (RMS) between the flipped left cornea and the right cornea was initially determined for both the anterior and posterior corneal surfaces. Next, the iterative closest point transformation algorithm was applied on the three-dimensional flipped cornea to allow the flipped left corneal anterior surface to translate and rotate, minimising the difference between it and the right corneal anterior surface. Then, the RMS differences were recalculated and compared. Results: A comparison of the dioptric powers showed a significant difference between the RMS of both the flipped left eyes and the right eyes in the healthy and the KC groups (p < 0.001). The RMS of the surfaces of the flipped left corneas and the right corneas was 0.6 ± 0.4 D among the healthy group and 4.1 ± 2.3 among the KC group. After transforming the flipped left corneas, the RMS was recorded as 0.5 ± 0.3 D and 2.4 ± 2 D among the healthy and KC groups, respectively. Conclusions: Although fellow eyes are highly related in their clinical parameters, they should be treated with care when one eye topography is flipped and processed with the other eye topography in an optic-related research analysis where translation might be needed. In KC, an asymmetric disease, it was observed that a portion of the asymmetry was due to a corneal apex shift interfering with the image acquisition. Therefore, transforming the flipped left eyes by rotation and translation results in a fairer comparison between the fellow KC corneas.

Effects of intracorneal ring segments implementation technique and design on corneal biomechanics and keratometry in a personalized computational analysis

The implementation of intracorneal ring segments (ICRS) is one of the successfully applied refractive operations for the treatment of keratoconus (kc) progression. The different selection of ICRS types along with the surgical implementation techniques can significantly affect surgical outcomes. Thus, this study aimed to investigate the influence of ICRS implementation techniques and design on the postoperative biomechanical state and keratometry results. The clinical data of three patients with different stages and patterns of keratoconus were assessed to develop a three-dimensional (3D) patient-specific finite-element model (FEM) of the keratoconic cornea. For each patient, the exact surgery procedure definitions were interpreted in the step-by-step FEM. Then, seven surgical scenarios, including different ICRS designs (complete and incomplete segment), with two surgical implementation methods (tunnel incision and lamellar pocket cut), were simulated. The pre- and postoperative predicted results of FEM were validated with the corresponding clinical data. For the pre- and postoperative results, the average error of 0.4% and 3.7% for the mean keratometry value ([Formula: see text]) were predicted. Furthermore, the difference in induced flattening effects was negligible for three ICRS types (KeraRing segment with arc-length of 355, 320, and two separate 160) of equal thickness. In contrast, the single and double progressive thickness of KeraRing 160 caused a significantly lower flattening effect compared to the same type with constant thickness. The observations indicated that the greater the segment thickness and arc-length, the lower the induced mean keratometry values. While the application of the tunnel incision method resulted in a lower [Formula: see text] value for moderate and advanced KC, the induced maximum Von Mises stress on the postoperative cornea exceeded the induced maximum stress on the cornea more than two to five times compared to the pocket incision and the preoperative state of the cornea. In particular, an asymmetric regional Von Mises stress on the corneal surface was generated with a progressive ICRS thickness. These findings could be an early biomechanical sign for a later corneal instability and ICRS migration. The developed methodology provided a platform to personalize ICRS refractive surgery with regard to the patient's keratoconus stage in order to facilitate the efficiency and biomechanical stability of the surgery.

Fibril density reduction in keratoconic corneas

This study aims to estimate the reduction in collagen fibril density within the central 6 mm radius of keratoconic corneas through the processing of microstructure and videokeratography data. Collagen fibril distribution maps and topography maps were obtained for seven keratoconic and six healthy corneas, and topographic features were assessed to detect and calculate the area of the cone in each keratoconic eye. The reduction in collagen fibril density within the cone area was estimated with reference to the same region in the characteristic collagen fibril maps of healthy corneas. Together with minimum thickness and mean central corneal refractive power, the cone area was correlated with the reduction in the cone collagen fibrils. For the corneas considered, the mean area of keratoconic cones was 3.30 ± 1.90 mm². Compared with healthy corneas, fibril density in the cones of keratoconic corneas was lower by as much as 35%, and the mean reduction was 17 ± 10%. A linear approximation was developed to relate the magnitude of reduction to the refractive power, minimum corneal thickness and cone area (R² = 0.95, p < 0.001). Outside the cone area, there was no significant difference between fibril arrangement in healthy and keratoconic corneas. The presented method can predict the mean fibril density in the keratoconic eye's cone area. The technique can be applied in microstructure-based finite-element models of the eye to regulate its stiffness level and the stiffness distribution within the areas affected by keratoconus.

Characterization of cone size and centre in keratoconic corneas

A novel method to locate the centre of keratoconus (KC) and the transition zone between the pathological area and the rest of the corneal tissue is proposed in this study. A spherical coordinate system was used to generate a spherical height map measured relative to the centre of the optimal sphere fit, and normal to the surface. The cone centre was defined as the point with the maximum height. Second derivatives of spherical height were then used to estimate the area of pathology in an iterative process. There was mirror symmetry between cone centre locations in both eyes. The mean distance between cone centre and corneal apex was 1.45 ± 0.25 mm (0.07-2.00), the mean cone height normal to the surface was 37 ± 23 µm (2-129) and 75 ± 45 µm (5-243) in the anterior and posterior surfaces, respectively. There was a significant negative correlation between the cone height and the radius of the sphere of optimal fit (p < 0.05 for both anterior and posterior surfaces). On average, posterior cone height was larger than the corresponding anterior cone height by 37 ± 24 µm (0-158). The novel method proposed can be used to estimate the cone centre and area, and explore the changes in anterior and posterior corneal surfaces that take place with KC progression. It can help improve understanding of keratoconic corneal morphology and assist in developing customized treatments.

Keratoconic eyes with stable corneal tomography could benefit more from custom intraocular lens design than normal eyes

We investigated whether eyes with keratoconic corneal tomography pattern could benefit more from aberration correction with custom intraocular lenses (IOLs) than normal cataractous eyes despite the effect of misalignment on the correction of aberrations. Custom IOLs (cIOLs) were calculated for twelve normal and twelve keratoconic eyes using personalized numerical ray tracing models. The Stiles-Crawford weighted root-mean-square spot-size (wRMS) at the virtual fovea was evaluated for cIOLs and aberration-neutral IOLs (nIOLs) in a simulated clinical study with 500 virtual IOL implantations per eye and per IOL. IOL misalignment (decentration, tilt, rotation) and pupillary ectopia (4.5 mm iris aperture) were varied upon each virtual implantation. The nIOLs achieved average wRMS of 16.4 ± 4.3 μm for normal, and 92.7 ± 34.4 μm for keratoconic eyes (mean ± standard deviation). The cIOLs reduced the average wRMS to 10.3 ± 5.8 μm for normal, and 28.5 ± 18.6 μm for keratoconic eyes. The cIOLs produced smaller wRMS than nIOLs in most virtual implantations (86.7% for normal and 99.4% for keratoconic eyes). IOL misalignment resulted in larger wRMS variations in the keratoconus group than in the normal group. Custom freeform IOL-optics-design may become a promising option for the correction of advanced aberrations in eyes with non-progressive keratoconic corneal tomography pattern.

Non-Orthogonal Refractive Lenses for Non-Orthogonal Astigmatic Eyes

Purpose: To present a novel design method for non-orthogonal lenses to reduce the problem of residual astigmatism in non-orthogonal, astigmatic eyes Methods: A method to create spectacle trial lenses with non-orthogonal power axes was developed based on a novel optimised light ray-tracing algorithm rather than conventional lens design methods which could not fully eliminate spherical aberration. Using this method, three sets of refraction trial lenses were made with the angles between power axes of each set controlled at 80°, 70° and 60°, respectively. Within each set, the cylindrical power varied from -1.00 D to -6.00 D in 1.00 D steps in addition to a -0.50 D lens. Computer-based numerical simulation of the lenses optical performance was carried out to apply orthogonal and non-orthogonal lenses on simulated astigmatic eyes. Subsequently, three clinical trial cases were investigated. Results: Computer-simulated optical performance of non-orthogonal lenses showed the ability to achieve high performance in correcting non-orthogonal astigmatism. Subsequently, three patients with irregular astigmatism were refracted with the non-orthogonal lens sets, and clinically observed improvement at least two lines in the LogMAR chart was achieved in all three cases, compared with correction with orthogonal lenses, along with subjective improvement in image quality. Conclusions: Non-orthogonal astigmatism, which is commonly ignored by current eye prescription systems, is taken into account in this study in the design of spectacle and soft contact lenses. The new approach considers the possible non-orthogonal positions of the eye's two optical power meridians and appears to be better able to correct the vision of irregular astigmatic eyes and significantly reduce residual astigmatism.