Next Article in Journal
Nearly Exact Discrepancy Principle for Low-Count Poisson Image Restoration
Previous Article in Journal
A Robust Tensor-Based Submodule Clustering for Imaging Data Using l12 Regularization and Simultaneous Noise Recovery via Sparse and Low Rank Decomposition Approach
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

On the Relationship between Corneal Biomechanics, Macrostructure, and Optical Properties

by
Francisco J. Ávila
*,
Maria Concepción Marcellán
and
Laura Remón
Departamento de Física Aplicada, Universidad de Zaragoza, 50009 Zaragoza, Spain
*
Author to whom correspondence should be addressed.
J. Imaging 2021, 7(12), 280; https://doi.org/10.3390/jimaging7120280
Submission received: 27 October 2021 / Revised: 15 December 2021 / Accepted: 16 December 2021 / Published: 18 December 2021
(This article belongs to the Special Issue Biomechanical Techniques for Biomedical Imaging)

Abstract

:
Optical properties of the cornea are responsible for correct vision; the ultrastructure allows optical transparency, and the biomechanical properties govern the shape, elasticity, or stiffness of the cornea, affecting ocular integrity and intraocular pressure. Therefore, the optical aberrations, corneal transparency, structure, and biomechanics play a fundamental role in the optical quality of human vision, ocular health, and refractive surgery outcomes. However, the inter-relationships of those properties are not yet reported at a macroscopic scale within the hierarchical structure of the cornea. This work explores the relationships between the biomechanics, structure, and optical properties (corneal aberrations and optical density) at a macro-structural level of the cornea through dual Placido–Scheimpflug imaging and air-puff tonometry systems in a healthy young adult population. Results showed correlation between optical transparency, corneal macrostructure, and biomechanics, whereas corneal aberrations and in particular spherical terms remained independent. A compensation mechanism for the spherical aberration is proposed through corneal shape and biomechanics.

1. Introduction

Corneal biomechanics is a branch of biomedical sciences that deals with the analysis of the stability of the tissue when an external load or pressure is applied [1,2] or when the intraocular pressure fluctuates. The biomechanical properties of the cornea are responsible for its shape and integrity, acting as a unique convergence point between balanced ductility to preserve aspherical geometry (and correct ocular refraction), stiffness to compensate the intraocular pressure, and an ultrastructure that allows optical transparency. Biomechanical properties of the cornea can be affected by systemic diseases such as diabetes [3,4] or sclerosis [5,6]. In particular, corneal keratoconus may compromise the biomechanical stability, modifying the microstructure [7], weakening mechanical strength [8] and leading to corneal protrusion [9,10], inducing optical aberrations [11,12] that reduce the quality of vision or lead to transplants in advanced stages. The clinical relevance of the study of corneal biomechanics reached special interest with the development of refractive surgery techniques to modify the optical power of the cornea by laser ablation [13] or lenticular extraction [14]. These techniques consist of modifying the lamellar structure of cornea, causing a redistribution of mechanical stress. The biomechanical response is expected to provide the correct corneal curvature [15], and together with optical transparency, normal vision.
On the other hand, corneal transparency has been explained from the hierarchical structure of the cornea. At the molecular scale, X-ray scattering revealed how the collagen ultrastructure within the stroma is responsible for the tridimensional microstructure and consequently for the macroscopic geometry and biomechanics [16].
To date, the maximum spatial resolution of structural hierarchy achieved in living human eyes has been the microscopic scale using two-photon scanning microscopy [17]. However, only confocal microscopes are currently clinically available, and although they allow visualization of the cellular matrix, they are invisible for the stromal architecture [18]. In this sense, Scheimpflug imaging provides excellent tomographic measurements of the macrostructure of the cornea as well as optical density (transparency) [19], which is widely reported in anterior segment analysis for the assessment of normal and keratoconus or ectatic corneas [20] or refractive surgery [21]. Corneal biomechanics is usually assessed employing air-puff tonometry [22]; also, the combination of Scheimpflug imaging and air-puff tonometry has been successfully integrated, bringing excellent results in dynamic assessment of corneal biomechanics [23].
As stated, the molecular organization of the corneal stroma controls the optical transparency, macroscopic shape, and structural stability (biomechanics) [16]. In this work, we will investigate if the relationship between corneal transparency and optical properties, geometry, and biomechanics is preserved at the macroscopic level of the hierarchical structure.
The biomechanical properties, corneal geometry, and optical properties and densitometry measurements were collected from 102 eyes of 51 young-adult healthy subjects using Scheimpflug imaging and air-puff tonometry.
This work focuses on the inter-relationships of corneal biomechanics, optical, and structural properties to bring a comprehensive macroscale characterization of the cornea that can also provide future predictive models of corneal biomechanics.

2. Materials and Methods

2.1. Participants

This research was reviewed by an independent Ethical Committee of Research of the Health Sciences Institute of Aragon (Spain) approved with reference: C.P.-C.I.PI20/377 (approval date: 14 July 2020). Measurements procedure and data collection were carried out according to the tenets of the Declaration of Helsinki. All participants were informed about the nature, risks, and possible adverse consequences of the study and signed an informed consent document. The ethnicity of the participants involved in this study was European Caucasian, all of them students from the School of Optics and Optometry of the University of Zaragoza (Spain). A total of 102 eyes from 51 healthy young subjects (mean age 24 ± 5) were analyzed using dual Placido–Scheimpflug imaging and air-puff tonometry systems. None of them presented ocular pathologies, corneal disorders, or abnormal intraocular pressure. Exclusion criteria included contact lens wearers due to the influence of contact lens wear in corneal optical density, thickness, and spherical aberration [24].

2.2. Experimental Measurements

Clinical measurements were carried out at the laboratory of Optometry of the Department of Applied Physics of the University of Zaragoza and conducted by an experienced clinical optometrist. Both eyes of all participants were analyzed in a sequential procedure: First, optical and geometrical properties (see Table 1) were acquired by a dual Scheimpflug–Placido disk imaging system, and next, corneal biomechanics was assessed using an air-puff tonometer device.

2.2.1. Dual Placido–Scheimpflug Imaging: Structural and Optical Parameters

The Galilei Dual Scheimpflug Analyzer (Ziemer Ophthalmic Systems AG, Port, Switzerland) is a clinical optical system that combines Placido Disk imaging and a revolving Scheimpflug camera providing simultaneous acquisition of corneal topography based on the internal and external surfaces including eccentricity, astigmatism, pachymetry (measures of central middle and peripheral cornea), three-dimensional analysis of the cornea, power measures, wavefront aberration, and optical densitometry [25]. Figure 1 shows an example of wavefront aberration mapping and anterior segment Scheimpflug image from a participant of our study. Table 1 summarizes those Galilei outputs selected for data analysis.

Optical Transparency Index Calculation

A Galilei device provides the degree of corneal and crystalline lens opacity by computing the relationship between the intensity of the illumination light and back-reflected rays expressed in standarized gray-scale units. It provides macroscopic visualization of the whole anterior segment of the eye by capturing full-angle imaging by rotating the Scheimpflug camera.
Scheimpflug imaging technology visualization is limited through strong scattering or opaque tissues such as the sclera, which is totally opaque due to irregular arrangement of collagen fibrils [16]. Nevertheless, toward the periphery of the iridocorneal angle in Scheimpflug images, a reference gray level from the sclera can be obtained from the Galilei outputs. In this work, the Optical Transparency Index (OTI) is defined and computed as:
OTI = 100 × [ 1 O D S l c e r a ( O D S c l e r a O D C o r n e a ) ] .
The OD of the cornea and sclera (ODCornea and ODSclera) was the average of OD values acquired at the same lateral location corresponding to 4 rotational images acquired at 0°, 45°, 90°, and 135° meridians. Regarding ODCornea, the lateral location was the optical axis reference, and the final OD value of each individual frame (i.e., corresponding to a given oriented meridian) was the average of the anterior, central, and posterior axial depths of the cornea. Figure 2 shows an example of an OD measurement at the posterior corneal location at a horizontal viewing angle.
OTI index ranges between 0 (total opacity) and 100 (absence of back-scattering light). The maximum value implies total corneal transmittance for the illumination wavelength underlying the understanding of corneal transparency.

2.2.2. Corneal Biomechanics Assessment

An Ocular Response Analyzer (ORA, Reichert Instruments, Depew, NY, USA) was employed to obtain measurements of corneal biomechanics. ORA is a non-invasive device based on air-puff applanation tonometry that measures intraocular pressure and corneal biomechanics, in particular corneal hysteresis (CH) and corneal resistance factor (CRF) parameters. Briefly, the corneal hysteresis can be defined as the energy dissipation when an external stress is applied, resulting in a time-dependent stain unlike purely elastic materials, which immediately recover their initial state once the stress stops. Whereas CH is a measurement of the viscoelasticity of the cornea, CRF is a measurement of the resistance (i.e., rigidity and/or elasticity) strongly associated to corneal thickness and then an indicator of the corneal pure elastic properties [26]. CRF and CH are related, but they do not describe the same biomechanical properties of the cornea.
CRF and CH are computed by quantifying the differential inward and outward corneal responses to an air pulse (see Figure 3) of approximately 24 milliseconds [27]. Once the first applanation is reached, the air pressure causes the cornea to move inward to a slight concavity and then back to a second applanation before recovering the natural shape.
Figure 3 shows an ORA measure from a participant of our study. From pressures P1 and P2 (corresponding to the first and second applanation pressures, respectively), CH and CRF are calculated as [28]:
CH = P1 − P2.
CRF = (P1 − 0.70 × P2) − 3.08.
From an elastodynamic point of view, the first applanation is given when the air pulse compensates for the intraocular and atmospheric pressures. The stress–strain response of the cornea undergoes elastic deformation, whereas the second applanation pressure is affected by energy dissipation; that is, P1 can be related to the pure elastic properties of the cornea, whereas P2 is used for viscoelasticity estimation. For instance, a pure elastic cornea will provide a null CH value and symmetrical pressure curve (see Figure 3) during an ORA measurement, which implies an applanation pressures equalization condition, P1 = P2. In that sense, a stiffness parameter has been reported as the ratio between the resultant pressure at the first applanation and the deflection amplitude [29]. In this study, Equations (1) and (2) were employed to calculate P1 and P2 from ORA outputs (CRF and CH).

2.2.3. Data Analysis

Dataset Clustering

Data segmentation is usually applied for dimensionality reduction. In this work, OTI calculations computed from OD measurements made by the Galilei analyzer were found to be discrete values ranging from 78.5 to 81.5 (0.5 step size) for all subjects involved in this study. Table 2 shows the number of subjects (N) that presented a representative OTI value within the range found experimentally, those data are normally distributed (Shapiro–Wilk test), which motivated data clustering for graphical representation and data analysis. Each cluster contains the information of the mean values of each corneal parameter and its standard deviation. In that sense, graphical representation and statistical analysis have been performed on clustered data.
Figure 4 shows the standard deviation of the computed OTI value as a function of the number of subject per discrete OTI value, all of the values were below 1% of error. Therefore, any significant variability in the variance of the clusters will be due to that corneal parameter to be correlated. In that sense, data clustering provides separate visualization of the entire dataset as a function of the meaningful features.
Collected data were stored into an Excel spreadsheet and then migrated to Origin Lab software (Origin Lab Corp., Northampton, MA, USA) for graphical representations and statistical analysis.
Statistical analysis consisted of Spearman Rank Order Correlation and linear regression analysis in order to establish or discard significant relationships between geometrical and optical parameters. The Shapiro–Wilk test was used to test the normality and make valid the inferences of linear regression analysis (power of performed tests = 0.995 with alpha = 0.05). The significance of the linear regression models was also indicated by the F-test. Limits of agreement were used to quantify the agreement between those parameters in graphical representations. Statistics was performed using the advanced statistical tool of Origin Lab software.

3. Results

This section shows results of our study that are structured in the analysis of the relationship between optical density of the cornea and three main factors: optical aberrations, biomechanics, and macroscopic geometry of the cornea.

3.1. Optical Density and Corneal Aberrations

Figure 1 showed an example of the total aberration wavefont map and the corresponding Scheimpflug image from a participant measured with a dual Galilei imaging system. Total and spherical aberration RMS and OTI average values (mean of 102 eyes measurements) are shown in Table 3. Since spherical aberration is the dominant high-order term at the cornea [30], it was evaluated apart from the total aberration RMS values. The statistical analysis (Spearman correlation) revealed no correlation between optical aberrations (neither total nor spherical term) and OTI values of the cornea. Whereas aberrations govern the correct focus of light and retinal image quality [31], it seems to be unrelated to the transparency of the cornea.

3.2. Optical Density and Corneal Macro-Structure

The cornea and sclera are composed mainly of type-I fibrillar collagen whose structural organization at the molecular scale is responsible for optical transparency. Due to the hierarchical structure of the cornea [32], the macroscopic structure is a consequence of the microscopic arrangement. This section investigates if the macroscopic structure of the cornea plays a role in the optical density or it disappears as a corneal transparency factor within the hierarachical organization. The total corneal astigmatism, anterior and posterior eccentricity, central, middle and peripheral thickness were measured and compared with the OTI index for all subjects.
Table 4 shows the mean values of thickness, corneal astigmatism, and eccentricity computed from the 102 measured eyes. The statistical analysis revealed that the corneal thickness and optical density are not related in young healthy subjects; however, both total corneal astigmatism and posterior eccentricity (i.e., the inner surface of the cornea) showed strong correlation with OTI parameters. These results imply that in the absence of pathological or physiological (aging) scattering contributions at the cornea, its optical transparency is related to the shape regardless of the thickness. Figure 5 shows the graphical representation of the cluster sampled OTI values as a function of total corneal astigmatism (Figure 5a) and posterior eccentricity (Figure 5b). The statistical analysis revealed positive correlation (Spearman, R2 = 0.87) between OTI and total corneal astigmatism (TCA), the regression analysis (overall F-test p = 0.003) confirmed that OTI and TCA are related.
Regarding the posterior eccentricity (PE) of the cornea, Figure 5b shows negative correlation (Spearman, R2 =0.94) between OTI and PE; the results of the linear regression analysis provided the significant relationship (overall F-test p =0.003). However, OTI and anterior eccentricity were found to be independent corneal parameters (see Table 3).

3.3. Optical Density and Corneal Biomechanics

This subsection investigates the impact of biomechanics on corneal transparency. Figure 6 shows a biomechanical image that represents the dynamic stress–strain response of the cornea during an air-puff applanation measurement at ORA device. The maximum deformation occurs at approximately 10 milliseconds. It can be observed how the distribution of the pressure curve is not symmetric but lopsided; this skewness is due to the viscoelastic nature of the cornea, while symmetric distributions correspond to corneas characterized by purely elastic biomechanical properties.
As described in the Methods section, the ORA device provides CH and CRF biomechanical parameters. From them, we derived the applanation pressures P1 and P2, which are related to the pure elastic (P1) and viscoelastic (P2) properties of the cornea. Figure 7 shows the graphical representation of the clustered sampling OTI values as a function of P1 (Figure 7a) and P2 (Figure 7b). The statistical analysis revealed significant (negative) correlations of OTI with P1 and P2 (R2 = 0.96 and R2 = 0.95, respectively). The error bars that correspond to the standard deviation of the cluster mean value are within the prediction bands of the regression analysis (overall F-tests p = 0.001 and p = 0.003, respectively). The statistical relationship between P1 and P2 with OTI proves that the optical density of the cornea is related to both elastic and viscosity properties from a macroscopic approach.

3.4. Experimental Expression to Predict the Human Corneal Transparency through Macroscopic Parameters

The results above described deal with the separated investigation of the relationship between optical density and corneal parameters that describe the optical, biomechanical, and macro-structural properties of the cornea. In this last subsection, the presence (or lack of) relationships between those characterization parameters are investigated.
Regarding the macro-structural approach, Figure 8a shows a regression plot of PE versus TCA cluster mean values. The statistical results (R2 = 0.88, p = 0.0028) revealed that as the total corneal astigmatism decreases, the posterior eccentricity significantly increases.
In this sense, it is obvious that the more direct the influence of the macro-structural shape on cornel aberrometry measurements, the less evident the relationship found between biomechanics and corneal aberrations. However, our findings did not show any statistical relationship in either total or spherical aberration (which is the dominant high-order term in the cornea) with respect to corneal biomechanics in young healthy subjects.
In that sense, corneal astigmatism is one of the low-order aberration terms that plays a fundamental role in ocular refraction; however, in terms of wave-front distortion, it does not increase the total RMS enough to show the dependence of total corneal aberrations with OTI or corneal biomechanics. However, understanding total corneal astigmatism as a measure of corneal asymmetry (structural parameter), Figure 8b shows the linear regression plot of applanation pressures (i.e., P1 and P2) versus total corneal astigmatism. The regression analysis revealed a strong dependence of the required pressures at ORA for first (R2 = 0.94, p = 0.001) and second ((R2 = 0.89, p = 0.002) applanations as the corneal astigmatism increases; that is, the higher the corneal asymmetry (measured by astigmatism), the lower the required air pressure to flatten the cornea. It is worth highlighting that the higher correlation value corresponds to the first applanation pressure, which is related to pure elastic properties unlike P2, which is related to viscous properties.
Finally, multiple regression analysis revealed that corneal transparency can be modeled from those macroscopic parameters that showed significant correlation. OTI parameter can be predicted from a linear combination (significant dependence) of the following variables: P1(p = 0.13), P2 (p = 0.22), TCA (p = 0.15), and PE (p = 0.021) (analysis of the variance of the multiple linear regression, p = 0.018) as:
OTI = 236,895 − [(7.817 × P1) + (0.561 × P2) + (4.458 × TCA) + (18.497 × PE)].

4. Discussion and Conclusions

Whereas the role of the corneal ultrastructure in the three-dimensional architecture at microscale, corneal transparency, and mechanical stability is well understood [16], the inter-relationships of those properties at the macroscale within the hierarchical corneal structure has not been investigated.
Considering that most of the medical devices for corneal analysis offer macro-structural resolution and information, in this study, we explore the relationship between the biomechanics, shape, and optical properties of the cornea in healthy young subjects. We analyzed wavefront aberration, optical density, geometrical parameters (thickness, eccentricity, and total corneal astigmatism) and corneal biomechanics using dual Placido–Scheimpflug analyzer and air-puff applanation tonometry systems, respectively.
Corneal aberrations and, in particular the spherical term, seem to have a compensation mechanism to keep stable the retinal image quality in the presence of intraocular scattering [31]; also, in traumatic injuries, corneal opacity appears together with high-order aberrations [33]. In particular, the decrease in optical density and corneal spherical aberration is correlated in contact lens wearers due to corneal swelling [24].
Although corneal astigmatism is one of the most representative low-order aberrations, in our work, we employed the total corneal astigmatism as a measure of the asymmetry in corneal shape (macro-structure parameter) instead of a wavefront distortion, since the contribution of astigmatism aberration term to the total aberration RMS was not enough to establish a relationship between corneal aberrations and OTI.
Regarding structural aspects, the corneal thickness and optical density change in dry eye syndrome, diabetes, and glaucoma [34] are correlated in edema processes such as corneal swelling associated with contact lens wear [24]. However, our findings showed that corneal thickness does not affect the optical transparency in healthy young subjects.
In addition, corneal surface irregularities are associated not only to optical aberrations but also to light scattering [35]. In that sense, our study also included as macro-structural parameters total corneal astigmatism and anterior and posterior eccentricities. As shown in the Results section, the higher the OTI value, the higher the corneal astigmatism and the lower the posterior eccentricity of the cornea. It is worth mentioning that only posterior corneal eccentricity depends on corneal astigmatism [36], and as shown in the Results section, the total corneal astigmatism and posterior corneal eccentricity were negatively correlated (R2 = 0.88, p = 0.028), which implies that the optical transparency and corneal macrostructure are related by means of the relative shape of the cornea independently of the thickness.
Focusing at the molecular scale, the structure of Type-I corneal collagen is responsible for the optical transparency and three-dimensional arrangement of the stroma [16]. The three-dimensional arrangement of the cornea determines its shape (geometry) and biomechanics [37]. Changes in corneal structure associated to aging or disease factors alter corneal biomechanics [38]. Thus, in our work, we also explored if biomechanics and corneal transparency are actually related at the macroscopic level of study. Results showed (see Figure 6) that both elastic- (P1) and viscoelastic-related parameters (P2) are strongly correlated with OTI; the viscoelastic parameter is more weakly correlated, which implies the dominance of the elastic property in corneal transparency.
Finally, a statistical relationship between biomechanics and corneal geometry (measured by total astigmatism) was found. That is, the corneal astigmatism, elastic, and viscoelastic properties of the cornea are related. It is worth mentioning that the dynamic corneal response depends not only on intrinsic biomechanical properties but also on intraocular pressure and corneal geometry [39]. In that sense, the ORA device provides corrected biomechanical measurements from intraocular pressure and central corneal thickness [23] but does not consider corneal astigmatism as a geometrical parameter affecting biomechanical assessment.
The relationship between corneal aberrations (in particular spherical term) and biomechanics has been previously stated in keratoconic eyes [40]. However, in young healthy subjects, the spherical aberration and corneal biomechanics are not related.
The independence of spherical aberration of corneal biomechanics in young healthy subjects could be explained by a feedback cycle compensation mechanism that occurs in the posterior cornea considering that biomechanical measurements are carried out at the anterior surface of the cornea: our results showed on the one hand that posterior eccentricity significantly increases as total astigmatism decreases, and on the other hand, the lower the total astigmatism of the cornea, the higher the applanation pressures at ORA (i.e., P1 and P2).
These results are consistent with the study reported by Li et al. [41], since their work concluded that the posterior corneal surface plays an important role in compensating for spherical aberration of the anterior corneal surface.
To conclude, within the hierarchical structure of the cornea, the nanoscopic scale leads to optical transparency, whereas the microscopic architecture models corneal biomechanics that impact the macroscale [16]. Previously to our work, Garzón et al. [42] reported a study on corneal densitometry and its correlation with aging, corneal thickness and curvature, and refractive error using Scheimpflug imaging; they found no correlation between optical transparency and refractive parameters. In our study, we expanded the analysis of those factors that could have an impact on corneal transparency.
In conclusion, if the macroscopic structure of the cornea is connected through optical and geometrical properties, its optical transparency can be modelled. Optical transparency measured through macroscale approaches is related to corneal biomechanics; in particular, the elastic property seems to be the dominant contribution. In addition, corneal astigmatism affects the biomechanical measurements in the sense that less applanation pressure is required as the total astigmatism increases, so future corrections must be calibrated into air-puff measurements to establish a compensation for asymmetric corneas and in short, avoid underestimation of the corneal biomechanics assessment.
Future research will expand this dataset to be incorporated into a predictive convolutional neural network including aging and pathological conditions to help us to develop predictive models in terms of the inter-relationships of the optical properties, structure, and biomechanics of the cornea.

Author Contributions

Conceptualization, F.J.Á.; methodology, F.J.Á.; validation, F.J.Á., L.R. and M.C.M.; formal analysis, F.J.Á.; investigation, F.J.Á.; data curation, F.J.Á.; writing—original draft preparation, F.J.Á.; writing—review and editing, F.J.Á., L.R. and M.C.M.; supervision, F.J.Á. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Institutional Review Board Statement

The study was conducted according to the guidelines of the Declaration of Helsinki and approved by the Ethics Committee of Research of the Health Sciences Institute of Aragon (Spain) with reference: C.P.-C.I. PI20/377. Approval date: 14 July 2020.

Informed Consent Statement

Informed consent was obtained from all subjects involved in the study.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Piñero, D.P.; Alcón, N. Corneal biomechanics: A review. Clin. Exp. Optom. 2015, 98, 107–116. [Google Scholar] [CrossRef]
  2. King, S.; Hafezi, F. Corneal biomechanics—A review. Opthalmic. Physiol. Opt. 2017, 37, 240–252. [Google Scholar] [CrossRef] [Green Version]
  3. Del Buey, M.A.; Casas, P.; Caramello, C.; López, N.; de la Rica, M.; Subirón, A.B.; Lanchares, E.; Huerva, V.; Grzybowski, A.; Ascaso, F.J. An update on corneal biomechanical cornal differences between type-2 diabetic and nondiabetic patients. J. Ophthalmol. 2019, 2019, 7645352. [Google Scholar] [PubMed]
  4. Beato, J.N.; Esteves-Leandro, J.; Reis, D.; Falcão, M.; Rosas, V.; Ângela, C.; Falcao-Reis, F. Structural and Biomechanical Corneal Differences between Type 2 Diabetic and Nondiabetic Patients. J. Ophthalmol. 2019, 2019, 3764878. [Google Scholar] [CrossRef]
  5. Ayhan, Z.; Kaya, M.; Ozturk, T.; Arikan, G.; Birlik, M. Association between Skin Thickness Measurements with Corneal Biomechanical Properties and Dry Eye Tests in Systemic Sclerosis. Ocul. Immunol. Inflamm. 2019, 27, 1138–1143. [Google Scholar] [CrossRef] [PubMed]
  6. Zimmermann, N.; Brandt, S.; Brünner, J.; Erb, C. Evaluation changes in corneal biomechanics in patients with systemic sclerosis. Klin. Monbl. Augenheilkd. 2019, 236, 806–815. [Google Scholar]
  7. Parker, J.S.; van Dijk, K.; Melles, G.R. Treatment options for advanced keratoconus: A review. Surv. Ophthalmol. 2018, 60, 459–480. [Google Scholar] [CrossRef]
  8. Roy, A.S.; Shetty, R.; Kummelil, M.K. Keratoconus: A biomechanical perspective on loss of corneal stiffness. Indian J. Ophthalmol. 2013, 61, 392–393. [Google Scholar] [CrossRef] [PubMed]
  9. Binder, P.S.; May, C.H.; Grant, S.C. An Evaluation of Orthokeratology. Ophthalmology 1980, 87, 729–744. [Google Scholar] [CrossRef]
  10. Ambrósio, R.; Correira, F.F.; Lopes, B.; Salomão, M.Q.; Luz, A.; Dawson, D.G.; Elsheikh, A.; Vinciguerra, R.; Vinciguerra, P.; Roberts, C.J. Corneal biomechanics in ectatic diseases: Refractive surgery implications. Open Ophthalmol. J. 2017, 11, 176–193. [Google Scholar] [CrossRef] [Green Version]
  11. Aksoy, S.; Akkaya, S.; Özkurt, Y.; Kurna, S.; Açıkalın, B.; Şengör, T. Topography and higher order corneal aberrations of the fellow eye in unilateral kerato-conus. Turk J. Ophthalmol. 2017, 47, 249–254. [Google Scholar] [CrossRef]
  12. Naderan, M.; Jahanrad, A.; Farjadnia, M. Ocular, corneal, and internal aberrations in eyes with keratoconus, forme fruste keratoconus, and healthy eyes. Int. Ophthalmol. 2018, 38, 1565–1573. [Google Scholar] [CrossRef]
  13. Wilkinson, J.M.; Cozine, E.W.; Kahn, A.R. Refractive Eye Surgery: Helping Patients Make Informed Decisions About LASIK. Am. Fam. Physician 2017, 95, 637–644. [Google Scholar] [PubMed]
  14. Zhang, Y.; Shen, Q.; Jia, Y.; Zhou, D.; Zhou, J. Clinical outcomes of SMILE and FS-LASIK used to treat myopia: A meta análisis. J. Refract. Surg. 2016, 32, 256–265. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  15. Sultan, M.B.; Mansberger, S.L.; Lee, P.P. Understanding the Importance of IOP Variables in Glaucoma: A Systematic Review. Surv. Ophthalmol. 2009, 54, 643–662. [Google Scholar] [CrossRef]
  16. Meek, K.M.; Knupp, C. Corneal structure and transparency. Prog. Retin. Eye Res. 2015, 49, 1–16. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  17. Ávila, F.J.; Gambín, A.; Artal, P.; Bueno, J.M. In vivo two-photon microscopy of the human eye. Sci. Rep. 2019, 9, 10121. [Google Scholar] [CrossRef] [Green Version]
  18. Petroll, W.M.; Robertson, D.M. In Vivo Confocal Microscopy of the Cornea: New Developments in Image Acquisition, Reconstruction, and Analysis Using the HRT-Rostock Corneal Module. Ocul. Surf. 2015, 13, 187–203. [Google Scholar] [CrossRef] [Green Version]
  19. Otri, A.M.; Fares, U.; Al-Aqaba, M.A.; Dua, H.S. Corneal Densitometry as an Indicator of Corneal Health. Ophthalmology 2012, 119, 501–508. [Google Scholar] [CrossRef] [PubMed]
  20. Belin, M.W.; Ambrósio, R. Scheimpflug imaging for keratoconus and ectatic disease. Indian J. Ophthalmol. 2013, 61, 401–406. [Google Scholar] [CrossRef] [PubMed]
  21. Boulze-Pankert, M.; Dariel, R.; Hoffart, L. Corneal Scheimpflug Densitometry Following Photorefractive Keratectomy in Myopic Eyes. J. Refract. Surg. 2016, 32, 788–791. [Google Scholar] [CrossRef] [PubMed]
  22. Qin, X.; Tian, L.; Zhang, H.; Chen, X.; Li, L. Evaluation of corneal elastic modulus based on Corneal Visualization Scheimpflug Technology. Biomed. Eng. Online 2019, 18, 42. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  23. Luce, D.A. Determining in vivo biomechanical properties of the cornea with an ocular response analyzer. J. Cataract. Refract. Surg. 2005, 31, 156–162. [Google Scholar] [CrossRef] [PubMed]
  24. Motevasseli, T.; Baradaran-Rafii, A.; Yazdizadeh, F.; Karimian, F.; Fekri, S.; Baradaran-Rafii, A. Comparison between two scheimpflug anterior segment analyzers. J. Ophthalmic Vis. Res. 2017, 12, 23–29. [Google Scholar] [CrossRef] [Green Version]
  25. Piñero, D.; Alcón, N. In vivo characterization of corneal biomechanics. J. Cataract. Refract. Surg. 2014, 40, 870–877. [Google Scholar] [CrossRef]
  26. Kaushik, S.; Pandav, S.S. Ocular Response Analyzer. J. Curr. Glaucoma Pract. 2012, 6, 17–19. [Google Scholar] [CrossRef]
  27. Lau, W.; Pye, D. A Clinical Description of Ocular Response Analyzer Measurements. Investig. Opthalmology Vis. Sci. 2011, 52, 2911–2916. [Google Scholar] [CrossRef] [Green Version]
  28. Roberts, C.J.; Mahmoud, A.M.; Bons, J.P.; Hossain, A.; Elsheikh, A.; Vinciguerra, R.; Vinciguerra, P.; Ambrósio, R., Jr. Introduction of two novel stiffness parameters and interpretation of air puff bio-mechanical deformation parameters with a Dynamic Scheimpflug analyzer. J. Refract. Surg. 2017, 33, 266–273. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  29. Llorente, L.; Barbero, S.; Merayo, J.; Marcos, S. Total and corneal optical aberrations induced by laser in situ keratomileusis for hypero-pia. J. Refract. Surg. 2004, 20, 203–216. [Google Scholar] [CrossRef]
  30. Pérez, G.M.; Manzanera, S.; Artal, P. Impact of scattering and spherical aberration in contrast sensitivity. J. Vis. 2009, 9, 19. [Google Scholar] [CrossRef]
  31. Ávila, F.J.; Ares, J.; Collados, M.V.; Marcellán, M.C.; Remón, L. An Update on Corneal Imaging Techniques: From Macroscale to Nanostructure. SN Compr. Clin. Med. 2019, 2, 1–10. [Google Scholar] [CrossRef] [Green Version]
  32. Shimizu, E.; Yamaguchi, T.; Tsubota, K.; Shimazaki, J. Corneal high-order aberrations in eyes with corneal scar after traumatic perfora-tion. Eye Contact Lens 2019, 45, 124–131. [Google Scholar] [CrossRef] [PubMed]
  33. Lombardo, M.; Rosati, M.; Pileri, M.; Schiano-Lomoriello, D.; Serrao, S. evaluation of corneal optical properties in subjects wearing hydrogel etafilcon A contact lenses and the effect of administering mannitol-enriched sodium hyaluronate ophthalmic solution. Clin. Ophthalmol. 2014, 8, 2347–2354. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  34. Wegener, A.R.; Meyer, L.; Hauch, C.; Sacchetti, M.; Plateroti, R.; Lambiase, A. Changes of Corneal Transparency and Thickness in Dry–Eye Syndrome Single and in Combination With Other Ocular Diseases. Investig. Opthalmol. Vis. Sci. 2005, 46, 4473. [Google Scholar]
  35. Spadea, L.; Maraone, G.; Verboschi, F.; Vingolo, E.M.; Tognetto, D. Effect of corneal light scatter on vision: A review of the literature. Int. J. Ophthalmol. 2016, 9, 459–464. [Google Scholar] [CrossRef]
  36. Mohammadi, S.-F.; Khorrami-Nejad, M.; Hamidirad, M. Posterior corneal astigmatism: A review article. Clin. Optom. 2019, 11, 85–96. [Google Scholar] [CrossRef] [Green Version]
  37. Whitford, C.; Studer, H.; Boote, C.; Meek, K.M.; Elsheikh, A. Biomechanical model of the human cornea: Considering shear stiffness and regional variation of collagen anisotropy and density. J. Mech. Behav. Biomed. Mater. 2015, 42, 76–87. [Google Scholar] [CrossRef]
  38. Blackburn, B.J.; Jenkins, M.W.; Rollins, A.M.; Dupps, W.J. A Review of Structural and Biomechanical Changes in the Cornea in Aging, Disease, and Photochemical Crosslinking. Front. Bioeng. Biotechnol. 2019, 7, 66. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  39. Ariza-Gracia, M.A.; Zurita, J.F.; Piñero, D.P.; Rodriguez-Matas, J.F.; Calvo, B. Coupled biomechanical response of the cornea assessed by non-contact tonom-etry. A simulation study. PLoS ONE 2015, 10, e0121486. [Google Scholar] [CrossRef] [Green Version]
  40. Piñero, D.P.; Alio, J.L.; Barraquer, R.I.; Michael, R.; Jiménez, R. Corneal biomechanics, refraction and corneal aberrometry in Keratoconus: An inte-grated study. Investig. Opthalmol. Vis. Sci. 2010, 51, 1948–1955. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  41. Li, X.; Wang, Y.; Dou, R. Aberration compensation between anterior and posterior corneal surfaces after Small incision lenticule extraction and Femtosecond laser-assisted laser in-situ keratomileusis. Ophthalmic Physiol. Opt. 2015, 35, 540–551. [Google Scholar] [CrossRef]
  42. Garzón, N.; Poyales, F.; Illarramendi, I.; Mendicute, J.; Jáñez, Óscar; Caro, P.; López-Muñoz, A.; Argüeso, F. Corneal densitometry and its correlation with age, pachymetry, corneal curvature, and refraction. Int. Ophthalmol. 2017, 37, 1263–1268. [Google Scholar] [CrossRef]
Figure 1. Total corneal wavefront aberration map (left) and anterior segment Scheimpflug image (right) from a volunteer of the study.
Figure 1. Total corneal wavefront aberration map (left) and anterior segment Scheimpflug image (right) from a volunteer of the study.
Jimaging 07 00280 g001
Figure 2. Optical density measurement at the posterior corneal location and horizontal viewing of the Scheimpflug camera.
Figure 2. Optical density measurement at the posterior corneal location and horizontal viewing of the Scheimpflug camera.
Jimaging 07 00280 g002
Figure 3. ORA measurement from a participant in our study. P1, P2, Max P, and CH correspond to the first and second applanation, maximum pressure, and corneal hysteresis, respectively.
Figure 3. ORA measurement from a participant in our study. P1, P2, Max P, and CH correspond to the first and second applanation, maximum pressure, and corneal hysteresis, respectively.
Jimaging 07 00280 g003
Figure 4. Standard deviation of OTI values as a function of the number of subjects per cluster.
Figure 4. Standard deviation of OTI values as a function of the number of subjects per cluster.
Jimaging 07 00280 g004
Figure 5. Mean clustered OTI values as function of total corneal astigmatism (TCA) (a) and posterior eccentricity (PE) (b) for all subjects. Standard deviation of the clustered data, equations, confident, and prediction bands of the regression analysis are included.
Figure 5. Mean clustered OTI values as function of total corneal astigmatism (TCA) (a) and posterior eccentricity (PE) (b) for all subjects. Standard deviation of the clustered data, equations, confident, and prediction bands of the regression analysis are included.
Jimaging 07 00280 g005
Figure 6. Dynamical representation of corneal applanation as a function of time from a volunteer of the study. Air pulse pressure is scaled in arbitrary units and shown in the bottom right corner legend.
Figure 6. Dynamical representation of corneal applanation as a function of time from a volunteer of the study. Air pulse pressure is scaled in arbitrary units and shown in the bottom right corner legend.
Jimaging 07 00280 g006
Figure 7. Mean clustered OTI values as function of first (a) and second applanation pressures (b) at the ORA device for all subjects. Standard deviation of the clustered data, equations, confidence, and prediction bands of the regression analysis are included.
Figure 7. Mean clustered OTI values as function of first (a) and second applanation pressures (b) at the ORA device for all subjects. Standard deviation of the clustered data, equations, confidence, and prediction bands of the regression analysis are included.
Jimaging 07 00280 g007
Figure 8. Mean clustered posterior eccentricity versus total corneal astigmatism values (a) and mean clustered applanation pressures at ORA versus total corneal astigmatism values (b). Linear regression fits are included.
Figure 8. Mean clustered posterior eccentricity versus total corneal astigmatism values (a) and mean clustered applanation pressures at ORA versus total corneal astigmatism values (b). Linear regression fits are included.
Jimaging 07 00280 g008
Table 1. Structural and optical parameters from the Galilei system considered of interest in our study. The total and spherical aberration wavefront were numerically evaluated by its root mean square (RMS) value.
Table 1. Structural and optical parameters from the Galilei system considered of interest in our study. The total and spherical aberration wavefront were numerically evaluated by its root mean square (RMS) value.
Structural ParametersOptical Parameters
Central ThicknessTotal Aberration RMS
Middle ThicknessSpherical Aberration RMS
Peripheral Thickness
Anterior Eccentricity
Posterior EccentricityOptical Density
Total Corneal Astigmatism
Table 2. Clustering of the OTI values as a function of the number of subjects (N).
Table 2. Clustering of the OTI values as a function of the number of subjects (N).
OTINCluster
78.5141
79.092
79.5223
80.0164
80.5155
81.0126
81.587
Table 3. Optical transparency index (OTI); Total aberration RMS (TA RMS); Spherical aberration RMS (SA RMS); and Spearman correlation analysis for the data collected from all participants of the study (* p < 0.05); ** p < 0.005).
Table 3. Optical transparency index (OTI); Total aberration RMS (TA RMS); Spherical aberration RMS (SA RMS); and Spearman correlation analysis for the data collected from all participants of the study (* p < 0.05); ** p < 0.005).
OTITA RMS * (µm)SA RMS ** (µm)Spearman’s *Spearman’s **
80.0 ± 0.91.45 ± 0.28−0.15 ± 0.05Failed, p = 0.38.Failed, p = 0.26.
Table 4. Mean values of geometrical parameters and correlation results with OTI index.
Table 4. Mean values of geometrical parameters and correlation results with OTI index.
Structural ParameterLocationMean ValueSpearman’s
Corneal ThicknessCentral554 ± 30 µmFailed, p = 0.58.
Middle601 ± 30 µmFailed, p = 0.42.
Peripheral674 ± 69 µmFailed, p = 0.29.
Corneal AstigmatismGlobal0.87 ± 0.34 Dpt.R2 = 0.87
Anterior EccentricityGlobal0.24 ± 0.16Failed, p = 0.39
Posterior EccentricityGlobal0.42 ± 0.12R2 = 0.94
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Ávila, F.J.; Marcellán, M.C.; Remón, L. On the Relationship between Corneal Biomechanics, Macrostructure, and Optical Properties. J. Imaging 2021, 7, 280. https://doi.org/10.3390/jimaging7120280

AMA Style

Ávila FJ, Marcellán MC, Remón L. On the Relationship between Corneal Biomechanics, Macrostructure, and Optical Properties. Journal of Imaging. 2021; 7(12):280. https://doi.org/10.3390/jimaging7120280

Chicago/Turabian Style

Ávila, Francisco J., Maria Concepción Marcellán, and Laura Remón. 2021. "On the Relationship between Corneal Biomechanics, Macrostructure, and Optical Properties" Journal of Imaging 7, no. 12: 280. https://doi.org/10.3390/jimaging7120280

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop