Multi-meridian corneal imaging of air-puff induced deformation for improved detection of biomechanical abnormalities

Improving early detection of keratoconus by Non Contact Tonometry. A computational study and new biomarkers proposal

Keratoconus is a progressive ocular disorder affecting the corneal tissue, leading to irregular astigmatism and decreased visual acuity. The architectural organization of corneal tissue is altered in keratoconus, however, data from ex vivo testing of biomechanical properties of keratoconic corneas are limited and it is unclear how their results relate to true mechanical properties in vivo. This study explores the mechanical properties of keratoconic corneas through numerical simulations of non-contact tonometry (NCT) reproducing the clinical test of the Corvis ST device. Three sensitivity analyses were conducted to assess the impact of corneal material properties, size, and location of the pathological area on NCT results. Additionally, novel asymmetry-based indices were proposed to better characterize corneal deformations and improve the diagnosis of keratoconus. Our results show that the weakening of corneal material properties leads to increased deformation amplitude and altered biomechanical response. Furthermore, asymmetry indices offer valuable information for locating the pathological tissue. These findings suggest that adjusting the Corvis ST operation, such as a camera rotation, could enhance keratoconus detection and provide insights into the relative position of the affected area. Future research could explore the application of these indices in detecting early-stage keratoconus and assessing the fellow eye's risk for developing the pathology.

Experimental evaluation of corneal stress-optic coefficients using a pair of force test

BACKGROUND

Topography and tomography are valuable techniques for measuring the corneal shape, but they cannot directly assess its internal mechanical stresses. And nonuniform corneal stress plays a crucial biomechanical role in the progression of diseases and postoperative changes. Given the cornea's inherent transparency, analyzing corneal stresses using the photoelasticity method is highly advantageous. However, quantification of photoelasticity faces challenges in obtaining the stress-optic coefficient due to wrinkles caused by the non-spherical geometry during tensional experiments.

OBJECTIVE

In this study, we propose an innovative experimental setup aimed at generating a gradient field of simple shear stress and achieving surface flatness during corneal stretching experiments, enabling the acquisition of the stress-optic coefficient through comparison with numerical results.

METHODS

Our designed setup applies fluid pressure and force couples on the cornea. The internal fluid pressure maintains the corneal shape, preventing wrinkles, while the force couples create a stress field leading to isochromatic fringes.

RESULTS

We successfully measured the stress-optic coefficients of the porcine anisotropic cornea in ex-vivo as 1.87 × 10-9 (horizontal) and 1.97 × 10-9 (vertical) (m2/N). Each isochromatic fringe order represents a shear stress range of 6.05 × 104 Pa under a low tension.

CONCLUSIONS

This study establishes a significant connection between corneal photoelastic patterns and the quantification of corneal stress by enabling direct measurement through advanced photoelastic visualization technology for clinical applications.

The extraction and application of antisymmetric characteristics of the cornea during air-puff perturbations

BACKGROUND

A non-contact tonometer is used to measure intraocular pressure, and studies have primarily relied on apex displacements to assess corneal properties. However, previous studies have overlooked the asymmetric characteristics of lateral corneal perturbations, leading to a gap in understanding of the lateral mechanical properties and its application.

METHOD

To investigate these lateral perturbations, we designed an experiment to sequentially record the corneal profiles when two consecutive air-puffs were applied at the center of the same cornea within a short period. Moreover, we used modal decomposition to decompose anterior surface profiles into symmetric and antisymmetric modes to comprehensively analyze the asymmetric characteristics. To extract mechanical properties, we utilized high-pass frequency analysis (>250 Hz) to filter out noise and errors.

RESULTS

Symmetric modes between the two consecutive air-puffs exhibited major similarities during vibration; however, antisymmetric modes exhibited minor differences in lateral perturbations of asymmetric vibration. The antisymmetric modes might be related to air-puff misalignment and mechanical properties. Through applying frequency analysis, the mechanical properties could be proven at high frequencies and misalignment shown at low frequencies. Furthermore, we compared the corneal vibration profiles of 259 healthy participants and 50 patients with keratoconus. Their properties showed that the antisymmetric modes of the keratoconus group exhibited a completely opposite direction of deformation compared to that in the healthy group.

CONCLUSIONS

Our proposed algorithm not only extracts antisymmetric characteristics but also offers valuable insights into decompose misalignment and mechanical properties of healthy and keratoconus corneas, presenting a new perspective for corneal biomechanics.

Long-range frequency-domain optical delay line based on a spinning tilted mirror for low-cost ocular biometry

Optical biometers are routinely used to measure intraocular distances in ophthalmic applications such as cataract surgery planning or myopia monitoring. However, due to their high cost and reduced transportability, access to them for screening and surgical planning is still limited in low-resource and remote settings. To increase patients' access to optical biometry we propose a novel low-cost frequency-domain optical delay line (FD-ODL) based on an inexpensive stepper motor spinning a tilted mirror, for integration into a time-domain (TD)-biometer, amenable to a compact footprint. In the proposed FD-ODL, the axial scan range and the A-scan rate are decoupled from one another, as the former only depends on the spinning mirror tilt angle, while the A-scan rate only depends on the motor shaft rotational speed. We characterized the scanning performance and specifications for two spinning mirror tilt angles, and compared them to those of the standard, more expensive FD-ODL implementation, employing a galvanometric scanner for group delay generation. A prototype of the low-cost FD-ODL with a 1.5 deg tilt angle, resulting in an axial scan range of 6.61 mm and an A-scan rate of 10 Hz was experimentally implemented and integrated in a dual sample beam optical low-coherence reflectometry (OLCR) setup with a detour unit to replicate the measurement window around the anterior segment and the retina. The intraocular distances of a model eye were measured with the proposed low-cost biometer and found to be in good agreement with those acquired by a custom swept-source optical coherence tomography (SS-OCT) system and two commercial biometers, validating our novel design.

Optomechanical assessment of photorefractive corneal cross-linking via optical coherence elastography

Purpose: Corneal cross-linking (CXL) has recently been used with promising results to positively affect corneal refractive power in the treatment of hyperopia and mild myopia. However, understanding and predicting the optomechanical changes induced by this procedure are challenging. Methods: We applied ambient pressure modulation based optical coherence elastography (OCE) to quantify the refractive and mechanical effects of patterned CXL and their relationship to energy delivered during the treatment on porcine corneas. Three different patterned treatments were performed, designed according to Zernike polynomial functions (circle, astigmatism, coma). In addition, three different irradiation protocols were analyzed: standard Dresden CXL (fluence of 5.4 J/cm2), accelerated CXL (fluence of 5.4 J/cm2), and high-fluence CXL (fluence of 16.2 J/cm2). The axial strain distribution in the stroma induced by ocular inflation (Δp = 30 mmHg) was quantified, maps of the anterior sagittal curvature were constructed and cylindrical refraction was assessed. Results: Thirty minutes after CXL, there was a statistically significant increase in axial strain amplitude (p < 0.050) and a reduction in sagittal curvature (p < 0.050) in the regions treated with all irradiation patterns compared to the non-irradiated ones. Thirty-6 hours later, the non-irradiated regions showed compressive strains, while the axial strain in the CXL-treated regions was close to zero, and the reduction in sagittal curvature observed 30 minutes after the treatment was maintained. The Dresden CXL and accelerated CXL produced comparable amounts of stiffening and refractive changes (p = 0.856), while high-fluence CXL produced the strongest response in terms of axial strain (6.9‰ ± 1.9‰) and refractive correction (3.4 ± 0.9 D). Tripling the energy administered during CXL resulted in a 2.4-fold increase in the resulting refractive correction. Conclusion: OCE showed that refractive changes and alterations in corneal biomechanics are directly related. A patient-specific selection of both, the administered UV fluence and the irradiation pattern during CXL is promising to allow customized photorefractive corrections in the future.

Optical beam scanner with reconfigurable non-mechanical control of beam position, angle, and focus for low-cost whole-eye OCT imaging

Whole-eye optical coherence tomography (OCT) imaging is a promising tool in ocular biometry for cataract surgery planning, glaucoma diagnostics and myopia progression studies. However, conventional OCT systems are set up to perform either anterior or posterior eye segment scans and cannot easily switch between the two scan configurations without adding or exchanging optical components to account for the refraction of the eye's optics. Even in state-of-the-art whole-eye OCT systems, the scan configurations are pre-selected and cannot be dynamically reconfigured. In this work, we present the design, optimization and experimental validation of a reconfigurable and low-cost optical beam scanner based on three electro-tunable lenses, capable of non-mechanically controlling the beam position, angle and focus. We derive the analytical theory behind its control. We demonstrate its use in performing alternate anterior and posterior segment imaging by seamlessly switching between a telecentric focused beam scan to an angular collimated beam scan. We characterize the corresponding beam profiles and record whole-eye OCT images in a model eye and in an ex vivo rabbit eye, observing features comparable to those obtained with conventional anterior and posterior OCT scanners. The proposed beam scanner reduces the complexity and cost of other whole-eye scanners and is well suited for 2-D ocular biometry. Additionally, with the added versatility of seamless scan reconfiguration, its use can be easily expanded to other ophthalmic applications and beyond.

Estimation of Crystalline Lens Material Properties From Patient Accommodation Data and Finite Element Models

Purpose

The mechanical properties of the crystalline lens are related to its optical function of accommodation, and their changes with age are one of the potential causes for presbyopia. We estimated the mechanical parameters of the crystalline lens using quantitative optical coherence tomography (OCT) imaging and wavefront sensing data from accommodating participants and computer modeling.

Methods

Full-lens shape data (from quantitative swept-source OCT and eigenlens representation) and optical power data (from Hartmann-Shack wavefront sensor) were obtained from 11 participants (22-30 years old) for relaxed accommodation at infinity and -4.5 D accommodative demand. Finite element models of lens, capsular bag, zonulae, and ciliary body were constructed using measured lens geometry and literature data, assuming a 60-mN radial force. An inverse modeling scheme was used to determine the shear moduli of the nucleus and cortex of the lens, such that the simulated deformed (maximally stretched) lens matched the participant's lens at -4.5 D.

Results

The shear moduli of the nucleus and cortex were 1.62 ± 1.32 and 8.18 ± 5.63 kPa, on average, respectively. The shear modulus of the nucleus was lower than that of the cortex for all participants evaluated. The average of the two moduli per participant was statistically significantly correlated with age (R2 = 0.76, P = 0.0049).

Conclusions

In vivo imaging and mechanical modeling of the crystalline lens allow estimations of the crystalline lens' mechanical properties. Differences between nuclear and cortical moduli and their dependency with age appear to be critical in accommodative function and likely in its impairment in presbyopia.

Dynamic topography analysis of the cornea and its application to the diagnosis of keratoconus

PROPOSE

To establish a dynamic topography analysis method which simulates the dynamic biomechanical response of the cornea and reveals the variations of such response within the corneal surface, and thereafter to propose and clinically evaluate new parameters for the definite diagnosis of keratoconus.

METHODS

58 normal (Normal) and 56 keratoconus (KC) subjects were retrospectively included. Personalized corneal air-puff model was established using corneal topography data by Pentacam for each subject, and the dynamic deformation under air-puff loading was simulated using finite element method, which then enabled calculations of corneal biomechanical parameters of the entire corneal surface along any meridian. Variations in these parameters across different meridians and between different groups were explored by two-way repeated measurement analysis of variance. New dynamic topography parameters were proposed as the range of the calculated biomechanical parameters within the entire corneal surface, and the AUC of ROC curve was used to compare the diagnostic efficiency of newly proposed and existing parameters.

RESULTS

Corneal biomechanical parameters measured in different meridians varied significantly which were more pronounced in KC group due to its irregularity in corneal morphology. Considering such between-meridian variations thus led to improved diagnostic efficiency of KC as presented by the proposed dynamic topography parameter rIR with an AUC of 0.992 (sensitivity: 91.1%, specificity: 100%), significantly better than the current topography and biomechanical parameters.

CONCLUSIONS

The diagnosis of keratoconus may be affected by the significant variations of corneal biomechanical parameters due to corneal morphology irregularity. By considering such variations, the current study established the dynamic topography analysis process which benefits from the high accuracy of (static) corneal topography measurement while improving its diagnosis capacity. The proposed dynamic topography parameters, especially the rIR parameter, showed comparable or better diagnostic efficiency for KC than existing topography and biomechanical parameters, which can be of great clinical significance for clinics without access to instrument for biomechanical evaluations.

Effect of fixational eye movements in corneal topography measurements with optical coherence tomography

There is an increasing interest in applying optical coherence tomography (OCT) to quantify the topography of ocular structures. However, in its most usual configuration, OCT data is acquired sequentially while a beam is scanned through the region of interest, and the presence of fixational eye movements can affect the accuracy of the technique. Several scan patterns and motion correction algorithms have been proposed to minimize this effect, but there is no consensus on the ideal parameters to obtain a correct topography. We have acquired corneal OCT images with raster and radial patterns, and modeled the data acquisition in the presence of eye movements. The simulations replicate the experimental variability in shape (radius of curvature and Zernike polynomials), corneal power, astigmatism, and calculated wavefront aberrations. The variability of the Zernike modes is highly dependent on the scan pattern, with higher variability in the direction of the slow scan axis. The model can be a useful tool to design motion correction algorithms and to determine the variability with different scan patterns.

In vivo corneal elastography: A topical review of challenges and opportunities

Clinical measurement of corneal biomechanics can aid in the early diagnosis, progression tracking, and treatment evaluation of ocular diseases. Over the past two decades, interdisciplinary collaborations between investigators in optical engineering, analytical biomechanical modeling, and clinical research has expanded our knowledge of corneal biomechanics. These advances have led to innovations in testing methods (ex vivo, and recently, in vivo) across multiple spatial and strain scales. However, in vivo measurement of corneal biomechanics remains a long-standing challenge and is currently an active area of research. Here, we review the existing and emerging approaches for in vivo corneal biomechanics evaluation, which include corneal applanation methods, such as ocular response analyzer (ORA) and corneal visualization Scheimpflug technology (Corvis ST), Brillouin microscopy, and elastography methods, and the emerging field of optical coherence elastography (OCE). We describe the fundamental concepts, analytical methods, and current clinical status for each of these methods. Finally, we discuss open questions for the current state of in vivo biomechanics assessment techniques and requirements for wider use that will further broaden our understanding of corneal biomechanics for the detection and management of ocular diseases, and improve the safety and efficacy of future clinical practice.

In-vivo high-speed biomechanical imaging of the cornea using Corvis ST and digital image correlation

In-vivo corneal biomechanical characterization has gained significant clinical relevance in ophthalmology, especially in the early diagnosis of eye disorders and diseases (e.g. keratoconus). In clinical medicine, the air-puff-based tonometers such as Ocular Response Analyzer (ORA) and Corvis ST have been used in the in-vivo biomechanical testing. In the test, the high-speed dynamic deformation of the cornea under air-puff excitation is analyzed to identify the abnormities in the morphological and biomechanical properties of the cornea. While most existing measurements reflect the overall corneal biomechanical properties, in-vivo high-speed strain and strain rate fields at the tissue level have not been assessed. In this study, 20 subjects were classified into two different groups: the normal (NORM, N = 10) group and the keratoconus (KC, N = 10) group. Image sequences of the horizontal cross-section of the human cornea under air puff were captured by the Corvis ST tonometer. The macroscale mechanical response of the cornea was determined through image analysis. The high-speed evolution of full-field corneal displacement, strain, velocity, and strain rate was reconstructed using the incremental digital image correlation (DIC) approach. Differences in the parameters between the NORM and KC groups were statistically analyzed and compared. Statistical results indicated that compared with the NORM group, the KC corneas absorbed more energy (KC: 8.98 ± 2.76 mN mm; NORM: 4.79 ± 0.62 mN mm; p-value <0.001) with smaller tangent stiffness (KC: 22.49 ± 2.62 mN/mm; NORM: 24.52 ± 3.20 mN/mm; p-value = 0.15) and larger maximum deflection (KC: 0.99 ± 0.07 mN/mm; NORM: 0.92 ± 0.06 mN/mm; p-value <0.05) on the macro scale. Further, we also observed that The maximum displacement (KC: 1.17 ± 0.06 mm; NORM: 1.06 ± 0.07 mm; p-value <0.005), velocity (KC: 236 ± 29 mm/s; NORM: 203 ± 17 mm/s; p-value <0.01), shear strain (KC: 24.43 ± 2.59%; NORM: 20.26 ± 1.54%; p-value <0.001), and shear strain rate (KC: 69.74 ± 11.99 s-1; NORM: 54.84 ± 3.03 s-1; p-value <0.005) in the KC group significantly increased at the tissue level. This is the first time that the incremental DIC method was applied to the in-vivo high-speed corneal deformation measurement in combination with the Corvis ST tonometer. Through the image registration using incremental DIC analysis, spatiotemporal dynamic strain/strain rate maps of the cornea can be estimated at the tissue level. This is constructive for the clinical recognition and diagnosis of keratoconus at a more underlying level.

Estimation of the full shape of the crystalline lens in-vivo from OCT images using eigenlenses

Quantifying the full 3-D shape of the human crystalline lens is important for improving intraocular lens power or sizing calculations in treatments of cataract and presbyopia. In a previous work we described a novel method for the representation of the full shape of the ex vivo crystalline lens called eigenlenses, which proved more compact and accurate than compared state-of-the art methods of crystalline lens shape quantification. Here we demonstrate the use of eigenlenses to estimate the full shape of the crystalline lens in vivo from optical coherence tomography images, where only the information visible through the pupil is available. We compare the performance of eigenlenses with previous methods of full crystalline lens shape estimation, and demonstrate an improvement in repeatability, robustness and use of computational resources. We found that eigenlenses can be used to describe efficiently the crystalline lens full shape changes with accommodation and refractive error.

Co-axial acoustic-based optical coherence vibrometry probe for the quantification of resonance frequency modes in ocular tissue

We present a co-axial acoustic-based optical coherence vibrometry probe (CoA-OCV) for vibro-acoustic resonance quantification in biological tissues. Sample vibrations were stimulated via a loudspeaker, and pre-compensation was used to calibrate the acoustic spectrum. Sample vibrations were measured via phase-sensitive swept-source optical coherence tomography (OCT). Resonance frequencies of corneal phantoms were measured at varying intraocular pressures (IOP), and dependencies on Young´s Modulus (E), phantom thickness and IOP were observed. Cycling IOP revealed hysteresis. For E = 0.3 MPa, resonance frequencies increased with IOP at a rate of 3.9, 3.7 and 3.5 Hz/mmHg for varied thicknesses and 1.7, 2.5 and 2.8 Hz/mmHg for E = 0.16 MPa. Resonance frequencies increased with thickness at a rate of 0.25 Hz/µm for E = 0.3 MPa, and 0.40 Hz/µm for E = 0.16 MPa. E showed the most predominant impact in the shift of the resonance frequencies. Full width at half maximum (FWHM) of the resonance modes increased with increasing thickness and decreased with increasing E. Only thickness and E contributed to the variance of FWHM. In rabbit corneas, resonance frequencies of 360-460 Hz were observed. The results of the current study demonstrate the feasibility of CoA-OCV for use in future OCT-V studies.

In situ measurement of the stiffness increase in the posterior sclera after UV-riboflavin crosslinking by optical coherence elastography

Scleral crosslinking may provide a way to prevent or treat myopia by stiffening scleral tissues. The ability to measure the stiffness of scleral tissues in situ pre and post scleral crosslinking would be useful but has not been established. Here, we tested the feasibility of optical coherence elastography (OCE) to measure shear modulus of scleral tissues and evaluate the impact of crosslinking on different posterior scleral regions using ex vivo porcine eyes as a model. From measured elastic wave speeds at 6 - 16 kHz, we obtained out-of-plane shear modulus value of 0.71 ± 0.12 MPa (n = 20) for normal porcine scleral tissues. After riboflavin-assisted UV crosslinking, the shear modulus increased to 1.50 ± 0.39 MPa (n = 20). This 2-fold change was consistent with the increase of static Young's modulus from 5.5 ± 1.1 MPa to 9.3 ± 1.9 MPa after crosslinking, which we measured using conventional uniaxial extensometry on tissue stripes. OCE revealed regional stiffness differences across the temporal, nasal, and deeper posterior sclera. Our results show the potential of OCE as a noninvasive tool to evaluate the effect of scleral crosslinking.

Micron-scale hysteresis measurement using dynamic optical coherence elastography

We present a novel optical coherence elastography (OCE) method to characterize mechanical hysteresis of soft tissues based on transient (milliseconds), low-pressure (<20 Pa) non-contact microliter air-pulse stimulation and micrometer-scale sample displacements. The energy dissipation rate (sample hysteresis) was quantified for soft-tissue phantoms (0.8% to 2.0% agar) and beef shank samples under different loading forces and displacement amplitudes. Sample hysteresis was defined as the loss ratio (hysteresis loop area divided by the total loading energy). The loss ratio was primarily driven by the sample unloading response which decreased as loading energy increased. Samples were distinguishable based on their loss ratio responses as a function loading energy or displacement amplitude. Finite element analysis and mechanical testing methods were used to validate these observations. We further performed the OCE measurements on a beef shank tissue sample to distinguish the muscle and connective tissue components based on the displacement and hysteresis features. This novel, noninvasive OCE approach has the potential to differentiate soft tissues by quantifying their viscoelasticity using micron-scale transient tissue displacement dynamics. Focal tissue hysteresis measurements could provide additional clinically useful metrics for guiding disease diagnosis and tissue treatment responses.

Estimation of scleral mechanical properties from air-puff optical coherence tomography

We introduce a method to estimate the biomechanical properties of the porcine sclera in intact eye globes ex vivo, using optical coherence tomography that is coupled with an air-puff excitation source, and inverse optimization techniques based on finite element modeling. Air-puff induced tissue deformation was determined at seven different locations on the ocular globe, and the maximum apex deformation, the deformation velocity, and the arc-length during deformation were quantified. In the sclera, the experimental maximum deformation amplitude and the corresponding arc length were dependent on the location of air-puff excitation. The normalized temporal deformation profile of the sclera was distinct from that in the cornea, but similar in all tested scleral locations, suggesting that this profile is independent of variations in scleral thickness. Inverse optimization techniques showed that the estimated scleral elastic modulus ranged from 1.84 ± 0.30 MPa (equatorial inferior) to 6.04 ± 2.11 MPa (equatorial temporal). The use of scleral air-puff imaging holds promise for non-invasively investigating the structural changes in the sclera associated with myopia and glaucoma, and for monitoring potential modulation of scleral stiffness in disease or treatment.

Heartbeat optical coherence elastography: corneal biomechanics in vivo

SIGNIFICANCE

Mechanical assessment of the cornea can provide important structural and functional information regarding its health. Current clinically available tools are limited in their efficacy at measuring corneal mechanical properties. Elastography allows for the direct estimation of mechanical properties of tissues in vivo but is generally performed using external excitation force.

AIM

To show that heartbeat optical coherence elastography (Hb-OCE) can be used to assess the mechanical properties of the cornea in vivo.

APPROACH

Hb-OCE was utilized to detect Hb-induced deformations in the rabbit cornea in vivo without the need for external excitation. Furthermore, we demonstrate how this technique can distinguish corneal stiffness between untreated (UT) and crosslinked (CXL) tissue.

RESULTS

Our results demonstrate that stiffness changes in the cornea can be detected using only the Hb-induced deformations in the cornea. Additionally, we demonstrate a statistically significant difference in strain between the UT and CXL corneas.

CONCLUSIONS

Hb-OCE may be an effective tool for assessing the mechanical properties of the cornea in vivo without the need for external excitation. This tool may be effective for clinical assessment of corneal mechanical properties because it only requires optical coherence tomography imaging and data processing.