The transplantation of human corneal grafts (keratoplasty) has been performed for more than 100 years.1 In 2012, more than 184,500 corneal transplantations were carried out worldwide, about 90% of them were conducted as penetrating keratoplasty.2 About 3,628 keratoplasties were performed in Germany in 2012,3 while 2 years later, the number increased by more than 57% (5,725).4 Thus, keratoplasty is recognized as the oldest, most frequently performed, and most successful transplantation of human tissue with a success rate of .5,6 Retrieval, cultivation, and storage of human ocular tissue are performed by eye banks, which are also in charge of the quality assurance and provision of the corneal grafts for keratoplasty. The donor cornea is usually harvested as a corneoscleral button. After excision, it is transferred to the organ culture, where it is preserved in a supplemented tissue culture medium that contains growth factors, antibiotics, antimycotics, and buffer components. This preservation maintains the corneal metabolism and causes folding of the posterior surface [Fig. 1(a)], stromal swelling, and epithelial shedding.7,8 A few days before transplantation, the swelling is reversed by re-immersing the cornea into a second culture medium additionally enriched with dextran as a dehydrating agent [Fig. 1(a)].
Being part of the human tissue and organ transplant system, the corneal graft is subject to a sophisticated inspection procedure to assure transplantation eligibility.910.11.–12 The frequency of discarded grafts due to contamination, insufficient endothelium, and other reasons is about 30% as stated by Reinshagen et al.3,4 Donor cornea examination is usually performed with a slit lamp biomicroscope and a specular microscope. Assessments include visual inspection for alterations in tissue homogeneity as well as the determination of the endothelial cell density. In addition, donor liquor and swabs undergo microbiological examination to exclude contamination of the corneal button. Pathological changes, such as opacification, scars, or inflammations, can be identified by these techniques. However, manifestations of previous intraocular or anterior segment surgery, such as refractive corneal procedures (arcuate cuts or intrastromal refractive surgery), laser photoablation surgery [e.g., laser in-situ keratomileusis (LASIK) and photorefractive keratectomy (PRK)], or morphological alterations such as pterygia, superficial scars, (subclinical) keratoconus, or corneal dystrophies, may pass visual inspection unnoticed.1314.–15 Such pathological grafts may cause complications during surgery and moreover impair the achievable visual outcome.16,17 Supplementary screening methods are required which are able to visualize the structure of the corneal tissue and to calculate the surface topography or tomography. An objective, semiautomated screening method may assist standardization. An increasing number of potential donors have had refractive surgery and, therefore, the corneas may be unsuitable for penetrating keratoplasty, but not for endothelial lamellar keratoplasty.18 This requires an increase in the sensitivity and specificity of screening methods.
While identifying donor corneas potentially being unsuitable for transplantation as part of the eye banks’ quality assurance procedure, the corneal grafts should stay within the sterile culture flask [Fig. 2(a)] filled with a culture medium for maintaining a sterile environment throughout any quality assessment procedure. Ultrasound biomicroscopy requires contact of the ultrasound probe with the measurement object, either with direct or with immersion fluid, and, therefore, appears unsuitable. Optical measurement techniques, however, are contact-free and can offer a better spatial resolution. The surfaces of the culture flask are not equipped with antireflective coating, therefore, they reflect a perceptible portion of light stronger than the signal received from the corneal button. Techniques based on the reflection of placido disc mires are, therefore, unsuitable. In addition, particles from the corneal metabolism accumulate within the culture medium and lead to increased scattering. Reflections and scattering affect the image contrast. Furthermore, the optical resolution must be sufficient to measure density gradients within the stromal tissue and to derive the topography of the corneal anterior and posterior surfaces in detail while also providing data on the corneal stroma. Tomographic imaging based on slit projection principle and optical coherence tomography (OCT) provide a high spatial resolution.19,20 In addition, OCT-based imaging outperforms slit-projection when dealing with scattering particles and more intrinsic noise in the image data.21,22
Initial attempts for the detection of ablation patterns induced by previous keratorefractive surgery (e.g., LASIK, PRK, and radial keratotomy) were performed using the Orbscan (Bausch + Lomb, Rochester), a Placido and slit scanning-based topographer, on enucleated complete eye globes in air.23 Other investigations employed OCT-based methods because they enable the measurement of corneoscleral buttons under sterile conditions.2425.26.–27 LASIK interfaces could be detected within the reflection pattern of the stromal tissue. Beyond that, Lin et al. estimated anterior and posterior curvatures by means of a parabolic fit. Janunts et al. showed that keratoconus and extreme topographic irregularities could be identified by analyzing the thickness distribution of the cornea.28
To compute geometric parameters for classifying donor corneas, an appropriate surface model has to be chosen. The first approach was limited to spherical models with restriction to 1 degree of freedom (one radius of curvature ). With the use of conic sections, the model is extended by another parameter, the conic constant . It describes the contour of a conic section and thus considers the typical prolate ellipsoidal shape of the corneal surfaces.29 However, the normal human cornea is not rotationally symmetric.30 Therefore, the model needs to be expanded to a more general representation treating toroidal surfaces and two radii of curvature in orthogonal meridians. The biconic model is based on an explicit equation that allows four degrees of freedom, namely two radii of curvature and two conic constants. In this manner, both the cornea’s asphericity and toricity are being addressed.3132.–33
The purpose of this study was to develop an OCT-based measurement procedure for automatic classification of corneal grafts that can be carried out during cultivation under sterile conditions.
Material and Methods
For the measurements, we used a clinical spectral-domain OCT system (SD-OCT) (Spectralis OCT, Heidelberg Engineering, Heidelberg, Germany) with a central wavelength of . To allow imaging of the cornea, it was equipped with an add-on objective lens (Anterior Segment Module, Heidelberg Engineering, Heidelberg, Germany). The rectangular polystyrene tissue cultivation flasks (Primaria Canted-Neck Cell Culture Flask, Corning Inc., New York) containing the corneoscleral grafts [Fig. 2(a)] were mounted in front of the OCT system [Fig. 2(b)].
The corneal grafts were gathered from the LIONS cornea bank Saar-Lor-Lux, Trier/Westpfalz in Homburg/Saar and measured inside the sterile culture flask. The graft is fixated onto a polyethylene cornea holder with four brackets [Bausch + Lomb, Heidelberg, Germany; Fig. 2(a)]. The measurements were conducted at room temperature of and within 1 day after reimmersement because corneal metabolism particles obfuscate the culture medium, and the folding distorts the posterior surface topography. After measurement, the culture flasks were returned to the storing environment at 36°C. The storage and transportation were done according to the German regulations for cornea banking. According to these the allowed temperature range in the second medium is between 10°C and 40°C.10
The data were collected on the basis of a volume reconstruction of cross-sectional images through the posterior surface of the cornea (Fig. 3). The OCT was manually aligned to capture the central part of the cornea. The average penetration depth was 1.927 mm and the axial resolution across the corneal diameter of 8.3 mm [Fig. 4(a)]. The horizontal resolution was . In total, the measured volume consisted of a vertical sequence of 81 cross-sectional images. With a distance of between each cross-section, the area had a vertical expanse of 5.52 mm. After acquisition of the vertical array (scan 1), the scan window of the OCT was rotated by 90 deg with the sample in place, and the acquisition was repeated (scan 2). The two volume scans were exported as video files (image stacks) and combined in one coordinate system.
The data analysis was implemented in MATLAB (The MathWorks Inc., Massachusetts). The distinct steps of the procedure are shown in Fig. 5. At first, the region of interest (ROI) was applied to cutoff mirror artifacts in the periphery arising from the Fourier transform of the OCT fringe signal.34 The central 5.66 mm from the horizontal dimension of the corneal diameter (-axis) and the central 5.037 mm from the vertical dimension (-axis) were selected as ROI from the exported image files [Fig. 4(b)]. The cropped images had a resolution of with a scale of laterally and axially. Afterward, the cropped images were denoised, and the corneal volume was reconstructed as described in the following section. The two orthogonal measurements were placed in a common coordinate system [Fig. 4(a)] following an edge detection of the corneal surfaces. With the aid of a biconic surface fit, the two orthogonal radii of curvature and conic constants were determined.
Denoising and edge detection
Prior to the edge detection, the image data had to be refined from the remaining artifacts of the flask wall, reflections, and particles while maintaining the anterior and posterior corneal surfaces. For that reason, the background noise was eliminated with a () median filter. The number of particularly bright spots, such as reflections from the central corneal surfaces or mirrored artifacts from the back wall of the flask varied along the corneal diameter. The contrast range was adjusted by remapping the gray values of each image linearly to the interval [0, 1]. Clustered background noise components and mirror artifacts were selectively removed by a morphological opening. The size of the structuring element was chosen to be large enough to remove salty noise while being as small as possible avoiding corneal edge distortion.
Variations in image brightness toward the periphery were corrected by adjusting the contrast in relation to the distance from the center. Remaining particles were removed by deleting small objects [Fig. 6(c)]. Missing parts of the corneal tissue due to the shadowing were filled up by a top-hat-transformation. This morphological operation performs another opening with a disk-shaped structuring element with a radius of 45 pixels. It was applied to the inverted image. As a result, the shadowings were opened as bright spots. The resulting image was subtracted from the original image. In this manner, the shadowed areas were closed. Finally, the edges were detected using the Canny method,35 starting with a low gradient threshold of 0.25 and including connected edge pixels up to a high gradient threshold of 0. The resulting edges were assembled as surface data.
Biconic surface fit
To determine two orthogonal radii of curvature and as well as the related conic constants and , we expanded the well-known equation for the height data of a biconic surface.
At first, the detected edges of the anterior and posterior surfaces from each scan were transformed into spatial coordinates. The penetration depth was defined as the -axis, the horizontal expanse of each cross-section as the -axis, and the vertical sequence of cross-sections as the -axis. Then the corneal volume data of scan 1 and scan 2 were combined in one coordinate system with and .
To extract the required parameters, the biconic surface fit had to be independent from the spatial alignment of the corneoscleral button within the culture flask in relation to the measuring system.
For that reason, the spatial coordinates and were transformed to and . The azimuthal rotation angle of the orthogonal meridians around the -axis was incorporated by multiplication of the coordinate vector with a rotation matrix . Translational shifts , of the anticipatory apex from the origin of the coordinate system were addressed by subsequent multiplication with a translation matrix T and yielded the new coordinate vector . The slant of the surface plane due to the varying fixation of the culture flask in relation to the OCT device, and the floating of the corneal graft was included by adding two tilted planes with the tilt angles and . Furthermore, the depth offset of the apex was subtracted from the resulting biconic equation. The final equation yields 10 degrees of freedom in terms of the corneal parameters , , , and and the spatial parameters , , , , , and , merged in the parameter vector :
The values of the corneal parameters for a given spatial setting were derived with an iterative nonlinear least squares fit with trust region algorithm. In the first step A of the fitting procedure (Fig. 7), the spatial parameters were estimated by fitting a rotationally symmetrical anterior and posterior surfaces with average radii of curvature ( for the anterior corneal surface and 6.4 mm for the posterior corneal surface) and conic constants ( anterior and posterior)36,37 to the measurement data. The initial spatial parameters of were set to 0. The root-mean-square error (RMSE) was calculated, and outliers were excluded from the data set.
In the second step B, was assessed for the average cornea using the RMSE from step A as an initial parameter for . All the data points with were excluded. Afterward, the tilt angles and were calculated for the average cornea with the value from step B. Data points with were discarded. The resulting data set was used as the initial data set for steps D and E of the fitting procedure.
The previously probed values for , , , , and were set as initial parameters for step D, and the limits of their 95% confidence interval were set as lower and upper constraints. Then and as well as were calculated iteratively until the fit converged, i.e., for each parameter of . Following each iteration, data points with or with as the corresponding quantiles were ruled out. In the last fitting step, E, and were estimated while again using as initial parameters and the limits of their 95% confidence interval as constraints until the fit converged. Once more, the data points of the upper and lower 1.25% quantiles of the residual were excluded from the initial data set after each iteration.
The orientation of the corneoscleral button on the corneal holder after excision and re-immersion is unknown. The flatter radius of curvature has been declared as , the steeper one as , and the associated conic constants as and , respectively.
From the topographic data, the central corneal thickness (CCT) was determined by calculating the absolute distance between the anterior apex and the posterior surface
The refractive powers in the natural environment (interface air/cornea/aqueous humor) were calculated in both meridians and for the anterior surface using Javal’s keratometric approximation with . The posterior surface power was calculated by means of the surface power equation with (aqueous, inner medium) and (cornea, outer medium) to compare the outcome with clinical assessment of a physiological cornea using a Scheimpflug tomographer. The total refractive power or “true net power” (TNP) was computed by using the mean radius of curvature for each surface. Corneal astigmatism was calculated by the absolute difference of or .
To validate the measurement principle and fitting procedure, a rotationally symmetric model cornea made from polymethylmethacrylate (PMMA) was measured at first. The resulting parameter matched the theoretical specifications of the model lens with a relative error of , , , , and . Then 23 corneal grafts from 18 donors (11 left and 12 right eyes) were screened. One graft was anamnestically suspected to have undergone LASIK surgery. The measurements were conducted 1 day after the transfer to the deswelling culture medium, and the transplantation was scheduled within the next 4 days. The corneal radii of curvature, values, and the estimated keratometric power values and astigmatism are shown in Table 1. The CCT was (455 to ).
Radii of curvature R, conic constants κ, and keratometric power for the flat and steep meridians of the anterior and posterior corneal surfaces (n=23 transplants). Data are given as mean±standard deviation. Anterior keratometric powers are given with respect to Javal’s keratometric index.
Three extreme cases are presented in Table 2. Cornea 1 represents a historically suspected post-LASIK cornea. Cornea 2 represents the highest amount of astigmatism, whereas cornea 3 represents the thickest graft. The corresponding refractive powers and CCT are displayed in Table 3.
Corneal parameters of the suspected post-LASIK cornea (cornea 1), a cornea with extreme curvature radii (cornea 2), and the thickest cornea (cornea 3).
|Rf (mm)||Rs (mm)||κf||κs|
Refractive powers and CCT of the suspected post-LASIK cornea (cornea 1), a cornea with extreme curvature radii (cornea 2), and the thickest cornea (cornea 3). Anterior keratometric powers are given with respect to Javal’s keratometric index.
|Parameters||Cornea 1||Cornea 2||Cornea 3|
In this study, the geometrical properties of corneal grafts were assessed by means of an OCT-based screening method. The donor corneas were maintained in their sterile cultivation flasks within the eye bank’s storage procedure. Using a biconic surface fit, two orthogonal radii of curvature and the corresponding conic constants were elicited from the anterior and posterior corneal surfaces, respectively. Refractive power and astigmatism were calculated in a way that allows direct comparison with in vivo measured corneas. The corneal thickness, particularly at the anterior apex (CCT), was derived from the determined corneal volume model. The median radii of curvature of , , , and correspond to average horizontal and vertical corneal curvatures found in previous studies in vivo.33,37,38 The mean CCT of is slightly higher than the average value found in normal human corneas ().39 This effect can be attributed to the cultivation process of the graft as the culture medium increases stromal hydration.40 We conducted our measurements within 1 day after the grafts had been reimmersed in the deswelling medium, and in some cases, this time span might have been insufficient. The deswelling of organ-cultured corneal tissue was investigated by Schnitzler et al. who found that the condition of the endothelium might have a significant effect on the deswelling aside from the osmotic effect of the culture medium.27 Thus, measuring corneal thickness before reimmersement and at a standardized time after reimmersement might be one indicative factor for endothelium activity.
Schnitzler et al. reported a minimal corneal thickness after 24 h of deswelling.27 We chose a time span of 1 day as a good compromise between the time for deswelling and residual OCT image contrast.
One cornea was anamnestically reported to have undergone refractive surgery (cornea 1). The type of refractive surgery, however, was unknown. However, the corneal parameters, refractive power, and the CCT () did not differ significantly from the median values. The visual inspection of cornea 1 did not reveal any signs for refractive surgery. Previous studies found that LASIK interfaces can be seen in the short term in the stromal reflectivity data as bright spots. For longer periods of time, the interface intensity fades and manifests in a lower internal reflectivity.26 Lin et al. found the LASIK corneas to have flatter anterior curvature radii, which could not be detected in cornea 1 of this study.
Cornea 2 with the flattest anterior curvature radius showed a large discrepancy between the two curvature radii and and is, therefore, causing the extreme outlier within the astigmatism values (20.84 D). The corresponding conic constants of the anterior surface also show extreme values ( and ). Although the posterior surface is within the average range, this cornea exhibits the highest anterior refractive power ().
With cornea 3, all corneal parameters except for the CCT () were close to the median values. However, the ratios of the anterior to posterior curvature radius (0.794 and 0.784 for the flatter and the steeper meridian, respectively) are lower than the average ratio (0.82).
Due to the manual mounting of the cultivation flask in front of the objective of the OCT, the exact position of the corneal graft was unknown and not reproducible. By involving rotation, tilt angles, and translation parameters, the automatic fitting procedure was made independent of the alignment of the donor cornea to the corneal holder and its positioning within the culture flask.
Because of the rotationally nonsymmetrical toroidal surface of the cornea, a general biconic representation was elected as a mathematical model for the nonlinear least squares fitting procedure. In contrast to the ellipsoidal model, the individual parameters are independent from each other.41 The denoising routine yielded a smooth edge that eliminates particles and reflections and preserves the corneal surface. During the fit, we used two different refinement criteria to exclude outliers in the residual data. The first restricted the data to 1.5·RMSE, the second one to 1.2·RMSE. These criteria were arbitrarily chosen as a compromise between quality of the fit and convergence of the fit while still considering a sufficient amount of data. On average, iterations were necessary until the fit converged (minimum: 6, maximum: 95). The average RMSE was (minimum: and maximum: ), and the highest values occurred when the number of iterations culminated.
Calculating the biconic surface with the average anterior radii of curvature ( and ) and conic constants ( and ), a difference of manifests in the peripheral regions of the measurement area ( and 2.83 mm) in an elevation difference of . Considering the depth resolution of and the mean RMSE of the fitting procedure, this calculation of and is rather inaccurate compared with clinical measurement methods. The average human corneal diameter (horizontal limbus-to-limbus distance) is 11 mm.36 Screening a larger, circular-shaped area would significantly improve the calculation of the conic constant. Preliminarily estimating the unknown location of the corneal graft in relation to the OCT device by fitting an average cornea reduces the number of fitting parameters. Especially in the peripheral regions, a difference in the elevation data can almost equally be assigned to either different corneal parameters or varying tilt and rotation angles. This may lead to an inconclusive fitting result. To eliminate this drawback, the spatial parameters have been estimated beforehand in steps A–D of the fitting procedure (Fig. 6). The accuracy of this anticipatory approach has been verified with the PMMA model lens. Tilt, translation (especially relative to the optical axis of the OCT’s objective lens), and asymmetry of the cornea will affect the quality of the measurement data and the surface-fitting procedure, as they might induce significant refraction or displacement of the OCT beam. In the fitting procedure tilt, translation and corneal asymmetry will lead to errors in the biconic fit due to fitting the data to different data points. Furthermore, keratometric values will be compromised as well. Therefore, some care should be taken for a properly placed specimen during the image acquisition procedure. This is, in fact, one of the limitations in the study, the upright measurement procedure using a clinical OCT and the fact that the cornea might most probably not be in its natural geometric state, will affect our results (compare cornea 2). Alternative fixation techniques for corneal grafts are available that might help to reduce this problem. In addition, a marker on the corneal graft may assist proper alignment to improve the results. However, regulatory decrees do not allow an easy change of the fixation procedure. Another limitation arises from the OCT technique. One should take care of correcting the refraction of the rays at the first refractive surface of the cornea, which is especially important when measuring the cornea in air due to the high difference in the refractive indices between cornea and air. Exact measurement of geometrical parameters also requires correction of refraction at the interface between culture medium and cornea. However, we assumed that this effect would be rather small due to the small difference between the refractive indices of these media. The resulting error in thickness measurement due to neglecting the optical path differences between the surrounding medium and the corneal tissue was estimated to be .
The use of Zernike polynomials for the reconstruction of corneal surfaces either as a global approach or as an additional decomposition of the residual error from a basic surface fit has been widely discussed in the literature.4142.–43 While optical aberrations can be sufficiently modeled by Zernike decomposition and might be beneficial in the detection of keratoconus,44,45 the coefficients yield no direct information about geometrical properties such as curvature or asphericity. Moreover, the number of coefficients necessary to minimize the RMSE when representing abnormal corneas increases rapidly.
Optical Coherence Tomography Technique
In comparison with other measurement methods like x-ray or magnetic resonance techniques, OCT offers the best spatial resolution and does not harm the examined tissue.20 Unlike confocal microscopy, it is independent from the aperture of the system and provides a higher penetration depth. For that reason, OCT has become a widespread measurement method in ophthalmology,46 and it has previously been applied for screening corneal grafts aside from imaging the cornea in vivo.22,24,26,28,40,47 Considering the necessity of a sterile environment for the graft, the use of a contact-free measurement technique is obligatory. Because of the reflective behavior of the flask walls, other noncontact, high-resolution measurement methods, such as slit-projection or Placido topography would require the exposure of the tissue to air. Ultrasound techniques necessitate the application of an immersion medium, which is less comfortable.
The imaging modality of the OCT device exerts an influence on the signal-to-noise ratio and duration of the measurement. SD-OCT imaging offers several advantages over time-domain OCT (TD-OCT) devices. The Spectralis OCT is a Fourier-domain OCT (FD-OCT) and has a broadband light source, which simultaneously collects the measurement signal of multiple wavelengths. Therefore, it works much faster than the sequential image acquisition of TD-OCT devices. Although FD-OCT devices show difficulties in eliminating speckle and electron noise, several studies have shown that they offer a higher sensitivity and better spatial resolution, especially for the transparent tissue of the eye, but not in highly scattering tissue.4849.–50
Using SD-OCT, eye banks’ donor corneas could be observed under sterile conditions in their culture flasks to detect major curvature abnormalities. Using a biconic surface fit with 10 degrees of freedom, the curvature radii and conic constants could be estimated for the anterior and posterior corneal surfaces.
To further elaborate the screening method, the accuracy can be refined by expanding the scan range. The development of a custom-built, portable screening device47,49,51,52 would facilitate the integration of the screening procedure into the eye banks’ cultivation procedure. Moreover, the donor–host matching could be improved by labeling the meridians of the corneoscleral button and extraction of dioptric values before excision. The influence of the excision, cultivation, and transplantation procedure on the imaging quality of the corneal graft could be assessed and compared with the clinical outcome of the transplantation. In this way, the transplants could be selected specifically and aligned during surgery to compensate preexisting host astigmatism. The intensity gradient of the stromal tissue and the thickness profile should be taken into account as well to identify previous LASIK and (preclinical) keratoconus on the corneal grafts.
We gratefully acknowledge the financial support provided by the Dr. Rolf M. Schwiete Foundation. We would like to thank Mr. Herbert Junger from Heidelberg Engineering for the support with the anterior segment module and the Lions Cornea Bank Saar/Trier/Westpfalz for providing the corneal grafts.
Angela Damian is a scientist at the Institute of Experimental Ophthalmology, Saarland University in Homburg/Saar. She received her bachelor’s degree in electrical engineering (energy and automation) from Frankfurt University of Applied Sciences in 2013 and her master’s degree in biomedical engineering from the University of Applied Sciences in Saarbruecken in 2015. Her research interests include the development of automatic optical measurement technology, diagnostic systems for the human eye, electrical systems, and neural engineering.
Berthold Seitz is a full professor and the chairman of the Department of Ophthalmology, Saarland University Medical Center in Homburg/Saar. He completed his residency in ophthalmology in 1993 and he is a renowned expert in corneal diseases, corneal surgery, eye banking, and corneal transplantation. He was appointed as the head of the University Eye Hospital at the Saarland Medical Center in 2006. His current research interests include keratoconus and new therapeutic options for corneal pathologies.
Achim Langenbucher is a full professor and the head of the Institute of Experimental Ophthalmology, Saarland University in Homburg/Saar. He received his diploma degree in electrical engineering and his doctoral degree from the University of Erlangen-Nuremberg in 1995. He was appointed as an associate professor for medical optics at the Institute of Medical Physics in Erlangen in 2005. He is the chairman of Experimental Ophthalmology, Saarland University since 2009. His research interests include laser therapy in ophthalmology and individualized ophthalmology.
Timo Eppig is an adjunct professor at the Institute of Experimental Ophthalmology, Saarland University in Homburg/Saar. He received his diploma degree in optical technology and image processing from the Darmstadt University of Applied Sciences in 2005 and his doctoral degree from the University Erlangen-Nuremberg in 2009. He works as a research associate at the Institute of Experimental Ophthalmology at Saarland University since 2009. His research interests include individualized ophthalmology, automated diagnostics, intraocular lenses, and optical design.