Novel biomarker of sphericity and cylindricity indices in volume-rendering optical coherence tomography angiography in normal and diabetic eyes: a preliminary study.

PURPOSE: Preliminary to evaluate geometric indices (vessel sphericity and cylindricity) for volume-rendered optical coherence tomography angiography (OCTA) in healthy and diabetic eyes. METHODS: Twenty-six eyes of 13 healthy subjects and 12 eyes of patients with central ischemic, non-proliferative diabetic retinopathy were included. OCTA volume and surface area of the foveal vessels were measured and compared to determine OCTA sphericity and cylindricity indices and surface efficiency (SE). RESULTS: The overall average OCTA volume in healthy was 0.49 ± 0.09 mm3 (standard deviation [SD]), compared to 0.44 ± 0.07 mm3 (SD) in the diabetic eyes (difference in means 0.06 mm3, p = 0.054). The overall average OCTA surface area in the healthy eyes was 87.731 ± 9.51 mm2 (SD), compared to 76.65 ± 13.67 mm2 (SD) in the diabetic eyes (difference in means 11.08 mm2, p = 0.021). In relation to total foveolar tissue volume, the proportion of blood vessels was 22% in healthy individuals and only 20% in diabetics. The difference between the groups was more pronounced with respect to the total OCTA surface area, with a decrease of 13% in diabetics. A diabetic eye was most likely using geometric vessel indices analysis if the sphericity value was ≥ 0.190, with a cylindricity factor of ≥ 0.001. Reproducibility of the method was good. CONCLUSIONS: A method for OCTA surface area and volume measurements was developed. The application of the novel OCTA sphericity and cylindricity indices could be suitable as temporal biomarker to characterize stable disease or disease progression and may contribute to a better understanding in the evolution of diabetic retinopathy.

Biomechanical Modeling of Pterygium Radiation Surgery: A Retrospective Case Study.

Pterygium is a vascularized, invasive transformation on the anterior corneal surface that can be treated by Strontium-/Yttrium90 beta irradiation. Finite element modeling was used to analyze the biomechanical effects governing the treatment, and to help understand clinically observed changes in corneal astigmatism. Results suggested that irradiation-induced pulling forces on the anterior corneal surface can cause astigmatism, as well as central corneal flattening. Finite element modeling of corneal biomechanics closely predicted the postoperative corneal surface (astigmatism error -0.01D; central curvature error -0.16D), and can help in understanding beta irradiation treatment. Numerical simulations have the potential to preoperatively predict corneal shape and function changes, and help to improve corneal treatments.

Biomechanical Modeling of Femtosecond Laser Keyhole endokeratophakia Surgery.

PURPOSE: To apply a finite element model to endokeratophakia and evaluate anterior and posterior corneal surface changes. METHODS: Spatial elevation data (Pentacam HR; Oculus, Wetzlar, Germany) were obtained for the front and back corneal surfaces of an eye prior to undergoing an endokeratophakia procedure. These were used to warp a spherical template finite element model of the cornea to create a patient-specific finite element mesh and the initial stress distribution was computed with an iterative approach. The finite element model (Optimeyes; Integrated Scientific Services, Biel, Switzerland) included non-linear elastic characteristics of the stroma. The endokeratophakia procedure was recreated in the model: a donor lenticule (-10.50 diopters [D], 5.75-mm zone, 127-µm thick) was inserted into a lamellar pocket (180-µm deep, 6.25-mm diameter) and two 2-mm small incisions were made at 150° and 330°. Anterior and posterior surfaces, computed by the finite element model, were compared to clinical data to assess accuracy and reliability of finite element modeling. RESULTS: The postoperative axial curvature produced by the model closely resembled the patient data; average curvature was 48.01 D clinically and 48.23 D in the simulation, and corneal astigmatism was 3.01 D clinically and 2.88 D in the simulation. The posterior best-fit sphere elevation map also matched the patient data, replicating inward bulging of the posterior surface by approximately 40 µm. Stress distribution modeling predicted a stress increase by 159.94% ± 73% in the cap and a stress decrease by 32.41% ± 21% in the stromal bed. CONCLUSIONS: Finite element modeling of the cornea reproduced the clinically observed anterior and posterior corneal surface changes following an endokeratophakia procedure. This case sets the stage for further study to refine and yield predictive finite element modeling for the evaluation of corneal refractive surgical procedures.

Microstructure-based numerical simulation of the mechanical behaviour of ocular tissue.

This paper aims to present a novel full-eye biomechanical material model that incorporates the characteristics of ocular tissues at microstructural level, and use the model to analyse the age-related stiffening in tissue behaviour. The collagen content in ocular tissues, as obtained using X-ray scattering measurements, was represented by sets of Zernike polynomials that covered both the cornea and sclera, then used to reconstruct maps of collagen fibril magnitude and orientation on the three-dimensional geometry of the eye globe. Fine-mesh finite-element (FE) models with eye-specific geometry were built and supported by a user-defined material model (UMAT), which considered the regional variation of fibril density and orientation. The models were then used in an iterative inverse modelling study to derive the material parameters that represent the experimental behaviour of ocular tissues from donors aged between 50 and 90 years obtained in earlier ex vivo studies. Sensitivity analysis showed that reducing the number of directions that represented the anisotropy of collagen fibril orientation at each X-ray scattering measurement point from 180 to 16 would have limited and insignificant effect on the FE solution (0.08%). Inverse analysis resulted in material parameters that provided a close match with experimental intraocular pressure-deformation behaviour with a root mean square of error between 3.6% and 4.3%. The results also demonstrated a steady increase in mechanical stiffness in all ocular regions with age. A constitutive material model based on distributions of collagen fibril density and orientation has been developed to enable the accurate representation of the biomechanical behaviour of ocular tissues. The model offers a high level of control of stiffness and anisotropy across ocular globe, and therefore has the potential for use in planning surgical and medical procedures.

Interdevice variability of central corneal thickness measurement.

PURPOSE: To evaluate variability of central corneal thickness measurement (CCT) devices using a hitherto unprecedented number of CCT devices. METHODS: CCT was measured consecutively in 122 normal corneas of 61 subjects with seven different devices using three distinct measurement technologies: Scheimpflug, Ultrasound, and Optical Coherence Tomography (OCT). Per device deviation from the mean CCT value per eye was used to determine which of the devices performed best, compared to the mean value. RESULTS: Cirrus OCT yielded the lowest deviation. Deviations of the individual devices from the mean CCT of each eye were (OS/OD) 12.8±5.0/14.9±9.4 µm for Topcon noncontact specular microscopy (NCSM), 11.3±5.9/10.6±7.3 µm for Pentacam, 10.7±5.2/10.4±4.8 µm for Spectralis OCT, 6.0±3.9/6.2±4.9 µm for Topcon DRI OCT, 5.1±3.4/5.9±10.3 µm for AngioVue OCT, 4.8±4.1/5.7±4.6 µm for US pachymetry, and 4.2±3.2/5.7±4.6 µm for Cirrus OCT. The maximum differences between US pachymetry and the other devices were very high (up to 120 µm). CONCLUSION: Central corneal thickness may be under- or overestimated due to high interdevice variations. Measuring CCT with one device only may lead to inappropriate clinical and surgical recommendations. OCT showed superior results.

Patient-specific finite-element simulation of the human cornea: a clinical validation study on cataract surgery.

The planning of refractive surgical interventions is a challenging task. Numerical modeling has been proposed as a solution to support surgical intervention and predict the visual acuity, but validation on patient specific intervention is missing. The purpose of this study was to validate the numerical predictions of the post-operative corneal topography induced by the incisions required for cataract surgery. The corneal topography of 13 patients was assessed preoperatively and postoperatively (1-day and 30-day follow-up) with a Pentacam tomography device. The preoperatively acquired geometric corneal topography - anterior, posterior and pachymetry data - was used to build patient-specific finite element models. For each patient, the effects of the cataract incisions were simulated numerically and the resulting corneal surfaces were compared to the clinical postoperative measurements at one day and at 30-days follow up. Results showed that the model was able to reproduce experimental measurements with an error on the surgically induced sphere of 0.38D one day postoperatively and 0.19D 30 days postoperatively. The standard deviation of the surgically induced cylinder was 0.54D at the first postoperative day and 0.38D 30 days postoperatively. The prediction errors in surface elevation and curvature were below the topography measurement device accuracy of ±5µm and ±0.25D after the 30-day follow-up. The results showed that finite element simulations of corneal biomechanics are able to predict post cataract surgery within topography measurement device accuracy. We can conclude that the numerical simulation can become a valuable tool to plan corneal incisions in cataract surgery and other ophthalmosurgical procedures in order to optimize patients' refractive outcome and visual function.