Purpose: To improve the effectivity of patient-specific finite element analysis (FEA) to predict refractive power change (RPC) in rigid Ortho-K contact lens fitting. Novel eyelid boundary detection is introduced to the FEA model to better model the effects of the lid on lens performance, and stress and strain outcomes are investigated to identify the most effective FEA components to use in modelling. Methods: The current study utilises fully anonymised records of 249 eyes, 132 right eyes, and 117 left eyes from subjects aged 14.1 ± 4.0 years on average (range 9 to 38 years), which were selected for secondary analysis processing. A set of custom-built MATLAB codes was built to automate the process from reading Medmont E300 height and distance files to processing and displaying FEA stress and strain outcomes. Measurements from before and after contact lens wear were handled to obtain the corneal surface change in shape and power. Tangential refractive power maps were constructed from which changes in refractive power pre- and post-Ortho-K wear were determined as the refractive power change (RPC). A total of 249 patient-specific FEA with innovative eyelid boundary detection and 3D construction analyses were automatically built and run for every anterior eye and lens combination while the lens was located in its clinically detected position. Maps of four stress components: contact pressure, Mises stress, pressure, and maximum principal stress were created in addition to maximum principal logarithmic strain maps. Stress and strain components were compared to the clinical RPC maps using the two-dimensional (2D) normalised cross-correlation and structural similarity (SSIM) index measure. Results: On the one hand, the maximum principal logarithmic strain recorded the highest moderate 2D cross-correlation area of 8.6 ± 10.3%, and contact pressure recorded the lowest area of 6.6 ± 9%. Mises stress recorded the second highest moderate 2D cross-correlation area with 8.3 ± 10.4%. On the other hand, when the SSIM index was used to compare the areas that were most similar to the clinical RPC, maximum principal stress was the most similar, with an average strong similarity percentage area of 26.5 ± 3.3%, and contact pressure was the least strong similarity area of 10.3 ± 7.3%. Regarding the moderate similarity areas, all components were recorded at around 34.4% similarity area except the contact pressure, which was down to 32.7 ± 5.8%. Conclusions: FEA is an increasingly effective tool in being able to predict the refractive outcome of Ortho-K treatment. Its accuracy depends on identifying which clinical and modelling metrics contribute to the most accurate prediction of RPC with minimal ocular complications. In terms of clinical metrics, age, Intra-ocular pressure (IOP), central corneal thickness (CCT), surface topography, lens decentration and the 3D eyelid effect are all important for effective modelling. In terms of FEA components, maximum principal stress was found to be the best FEA barometer that can be used to predict the performance of Ortho-K lenses. In contrast, contact pressure provided the worst stress performance. In terms of strain, the maximum principal logarithmic strain was an effective strain barometer.

Background: The aim is to investigate induced higher-order aberrations (HOA)s and astigmatism as a result of non-toric ortho-k lens decentration and utilise artificial intelligence (AI) to predict its magnitude and direction. Methods: Medmont E300 Video topographer was used to scan 249 corneas before and after ortho-k wear. Custom-built MATLAB codes extracted topography data and determined lens decentration from the boundary and midpoint of the central flattened treatment zone (TZ). An evaluation was carried out by conducting Zernike polynomial fittings via a computer-coded digital signal processing procedure. Finally, an AI-based machine learning neural network algorithm was developed to predict the direction and magnitude of TZ decentration. Results: Analysis of the first 21 Zernike polynomial coefficients indicate that the four low-order and four higher-order aberration terms were changed significantly by ortho-k wear. While baseline astigmatism was not correlated with lens decentration (R = 0.09), post-ortho-k astigmatism was moderately correlated with decentration (R = 0.38) and the difference in astigmatism (R = 0.3). Decentration was classified into three groups: ≤0.50 mm, reduced astigmatism by −0.9 ± 1 D; 0.5~1 mm, increased astigmatism by 0.8 ± 0.1 D; >1 mm, increased astigmatism by 2.7 ± 1.6 D and over 50% of lenses were decentred >0.5 mm. For lenses decentred >1 mm, 29.8% of right and 42.7% of left lenses decentred temporal-inferiorly and 13.7% of right and 9.4% of left lenses decentred temporal-superiorly. AI-based prediction successfully identified the decentration direction with accuracies of 70.2% for right and 71.8% for left lenses and predicted the magnitude of decentration with root-mean-square (RMS) of 0.31 mm and 0.25 mm for right and left eyes, respectively. Conclusions: Ortho-k lens decentration is common when fitting non-toric ortho-k lenses, resulting in induced HOAs and astigmatism, with the magnitude being related to the amount of decentration. AI-based algorithms can effectively predict decentration, potentially allowing for better control over ortho-k fitting and, thus, preferred clinical outcomes.

To investigate the relationship between Ortho-K contact lens design parameters and refractive power change of the eye through a parametric mathematical representation. The current study utilises fully anonymized records of 249 eyes, 132 right eyes, and 117 left eyes from subjects aged 14.1 ± 4.0 years on average (range 9–38 years) which were selected for secondary analysis processing. The data were split into 3 groups (G1 up to 35 days wear, from 10 to 35 days, G2 up to 99 days wear, more than 35–99 days & G3 more than 100 days wear) according to the length of time, in days, that the lenses were worn. Corneal shape was measured before and after contact lens wear using the Medmont E300 topographer, from which height and distance files were read by a custom-built MATLAB code to construct the corneal anterior surface independently. Changes in refractive power pre and post-Ortho-K wear were determined using constructed tangential refractive power maps from which both centrally flattened and annular steepened zones were automatically bounded, hence used to determine the refractive power change. On average, flat Sim-K and steep Sim-K were reduced after Ortho-K lens wear by 1.6 ± 1.3 D and 1.3 ± 1.4 D respectively. The radius of the base curve was correlated with the mean central flattened zone power change strongly in G1 (R = 0.7, p < 0.001) and moderately in G2 (R = 0.4) and G3 (R = 0.4, p < 0.001). Hence, a strong correlation with the base curve was recorded in group G1 and moderate in G2 and G3. The reverse curve was very strongly correlated to the mean central flattened zone power change in G1 (R = 0.8, p < 0.001) and strongly correlated with G2 (R = 0.6, p < 0.001) and G3 (R = 0.7, p < 0.001). The reverse curve was also strongly correlated with the mean annular steepened zone power change among all groups G1, G2, and G3 (R = 0.7, R = 0.6 and R = 0.6) respectively (p < 0.001). Although the central corneal refractive power change was strongly correlated to the Ortho-K lens base curve, it characterized only 50% of the target power change. However, the annular steepened zone refractive power change appears to be a clearer predictor of target power change, as there appears to be a one-to-one inverse relationship with the target refractive power correction. Differences between these results and the literature may be a result of the topography software smoothing effect.

Soft contact lenses are usually characterised at room temperature, yet they function on the eye at body temperature, where their mechanics and optical performance can change. This study investigated whether soft silicone hydrogel lenses become more compliant in a physiological environment. Two silicone hydrogel materials (Definitive 74 and Unisil) were tested at 24 °C and 35 °C using uniaxial tensile and compression methods, with Ogden hyperelastic models fitted and finite element analysis performed on a realistic eye model. Both materials became more compliant at 35 °C, with Definitive 74 showing a larger modulus decrease (0.40 to 0.32 MPa) than Unisil (0.73 to 0.70 MPa). Finite element simulations indicated that these temperature-driven changes in compliance significantly affected refractive power, especially when the lens base curve exceeded the corneal radius by more than 5%. These findings demonstrate that soft silicone hydrogel lenses are indeed more compliant in a warm, hydrated environment, highlighting the need for physiologically relevant testing to inform design, fitting strategies, comfort, and vision outcomes.

Purpose: To evaluate whether corneal topography map stitching can fully substitute the traditional single-shot capture methods in clinical settings. Methods: This record review study involved the measurement of corneal surfaces from 38 healthy subjects using two instruments: the Medmont Meridia, which employs a stitching composite topography method, and the Eye Surface Profiler (ESP), a single-shot measurement device. Data were processed separately for right and left eyes, with multiple gaze directions captured by the Medmont device. Surface registration and geometric transformation estimation, including neighbouring cubic interpolation, were applied to assess the accuracy of stitched maps compared to single-shot measurements. Results: The study evaluated eye rotation angles and surface alignment between Medmont topography across various gaze directions and ESP scans. Close eye rotations were found in the right-gaze, left-gaze and up-gaze directions, with rotation angles of around 8°; however, the down-gaze angle was around 15°, almost twice other gaze rotation angles. Root mean squared error (RMSE) analysis revealed notable discrepancies, particularly in the right-, left-, and up-gaze directions, with errors reaching up to 98 µm compared to ESP scans. Additionally, significance analyses showed that surface area ratios highlighted considerable differences, especially in the up-gaze direction, where discrepancies reached 70% for both right and left eyes. Conclusions: Despite potential benefits, the findings highlight a significant mismatch between stitched and single-shot measured surfaces due to digital processing artefacts. Findings suggest that stitching techniques, in their current form, are not yet ready to substitute single-shot topography measurements fully. Although stitching helps fit large-diameter contact lenses, care should be taken regarding the central area, especially if utilising the stitched data for optimising optics or wavefront analysis.

Purpose: To investigate regional changes in tear film quality associated with orthokeratology (Ortho-K) lens wear using high-resolution spatial mapping and to evaluate the potential of artificial intelligence (AI) models in anticipating these changes. Methods: This study analysed tear film quality in 92 Ortho-K wearers divided into three groups based on lens wear duration (10–29 days, 30–90 days, and ≥91 days). Placido-based topographer measurement was used to generate regional tear film maps before and after treatment. A custom MATLAB pipeline enabled regional comparisons and statistical mapping. A feedforward neural network was trained to forecast local tear film quality using spatial data. Results: Single-value global mean metrics showed minimal changes in tear film quality across groups. However, regional mean mapping revealed significant mid-peripheral and peripheral deterioration over time, particularly in nasal and temporal corneal zones. These changes were often overlooked by global averaging and remained invisible through tear film breakup time (TBUT) measurements. The AI model predicted spatial tear quality with high accuracy (R ≥ 0.9 in testing), capturing nuanced regional variations. Conclusions: The regional analysis uncovers subtle, clinically relevant tear film disruptions caused by Ortho-K lens wear, particularly in peripheral areas. These insights challenge the adequacy of traditional single-value global mean assessments. The AI model demonstrates the potential for non-invasive, predictive evaluation of tear stability, supporting more personalised and effective Ortho-K care.