Optical coherence tomography-based topography determination of corneal grafts in eye bank cultivation

Abstract. Vision loss due to corneal injuries or diseases can be treated by transplantation of human corneal grafts (keratoplasty). However, quality assurance in retrieving and cultivating the tissue transplants is confined to visual and microbiological testing. To identify previous refractive surgery or morphological alterations, an automatic, noncontact, sterile screening procedure is required. Twenty-three corneal grafts have been measured in organ culture with a clinical spectral-domain optical coherence tomographer. Employing a biconic surface fit with 10 degrees of freedom, the radii of curvature and conic constants could be estimated for the anterior and posterior corneal surfaces. Thereupon, central corneal thickness, refractive values, and astigmatism have been calculated. Clinical investigations are required to elaborate specific donor–host matching in the future.


Introduction
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 ∼90%. 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. [9][10][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. [13][14][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 scanningbased 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. 24-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 R). 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. [31][32][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.

Experimental Setup
For the measurements, we used a clinical spectral-domain OCT system (SD-OCT) (Spectralis OCT, Heidelberg Engineering, Heidelberg, Germany) with a central wavelength of λ ¼ 870 nm. 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 25 cm 2 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 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 ∼23°C 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 3.893 μm across the corneal diameter of 8.3 mm [ Fig. 4(a)]. The horizontal resolution was 10.821 μm. In total, the measured volume consisted of a vertical sequence of 81 cross-sectional images. With a distance of 69 μm 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.

Data Analysis
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 (x-axis) and the central 5.037 mm from the vertical dimension (y-axis) were selected as ROI from the exported image files [ Fig. 4(b)]. The cropped images had a resolution of 568 × 496 pixels with a scale of 10.821 μm∕pixel laterally and 3.893 μm∕pixel 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.
For that reason, the background noise was eliminated with a (5 × 5) 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 R x and R y as well as the related conic constants κ x and κ y , we expanded the well-known equation for the height data zðx; yÞ 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 z-axis, the horizontal expanse of each cross-section as the x-axis, and the vertical sequence of cross-sections as the y-axis. Then the corneal volume data of scan 1 and scan 2 were combined in one coordinate system zðx; 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 x and y were transformed to x 0 and y 0 . The azimuthal rotation angle α of the orthogonal meridians around the z-axis was incorporated by multiplication of the coordinate vector X ¼ ½xy T with a rotation matrix R. Translational shifts x 0 , y 0 of the anticipatory apex z ¼ 0 from the origin of the coordinate system were addressed by subsequent multiplication with a translation matrix T and yielded the new coordinate vector X new ¼ XRT. 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 z 0 of the apex was subtracted from the resulting biconic equation. The final equation z model ðx 0 ; y 0 ; εÞ yields 10 degrees of freedom in terms of the corneal parameters R x , R y , κ x , and κ y and the spatial parameters x 0 , y 0 , z 0 , α, β, and γ, merged in the parameter vector ε ¼ ½R x ; R y ; κ x ; κ y x 0 ; y 0 ; z 0; α; β; γ: E Q -T A R G E T ; t e m p : i n t r a l i n k -; e 0 0 1 ; 3 2 6 ; 3 5 9 z model ðx 0 ; y 0 ; εÞ ¼ 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 (R x ¼ R y ¼ 7.8 mm for the anterior corneal surface and 6.4 mm for the posterior corneal surface) and conic constants (κ x ¼ κ y ¼ −0.2 anterior and −0.4 posterior) 36,37 to the measurement data. The initial spatial parameters of ε 0 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, z 0 was assessed for the average cornea using the RMSE from step A as an initial parameter for ε 1 . All the data points z i with residualðiÞ > RMSE · 1.5 were excluded. Afterward, the tilt angles β and γ were calculated for the average cornea with the z 0 value from step B. Data points z i with residualðiÞ > RMSE · 1.2 were discarded. The resulting data set z i;reduced ðx reduced ; y reduced Þ was used as the initial data set for steps D and E of the fitting procedure. The previously probed values for β, γ, x 0 , y 0 , and z 0 were set as initial parameters for step D, and the limits of their 95% confidence interval were set as lower and upper constraints. Then R x and R y as well as α were calculated iteratively until the fit converged, i.e., jε n iteration−1 −ε n iteration j ε n iteration−1 < 10 −6 for each parameter ε n of ε ¼ ðε 1 ; : : : ; ε m Þ. Following each iteration, data points with residualðiÞ > Q 0.0125 or residualðiÞ > Q 0.9875 with Q as the corresponding quantiles were ruled out. In the last fitting step, E, κ x and κ y were estimated while again using ε iteration−1 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 z i;reduced after each iteration.

Additional parameters
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 R f , the steeper one as R s , and the associated conic constants as κ f and κ s , 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 E Q -T A R G E T ; t e m p : i n t r a l i n k -; e 0 0 2 ; 3 2 6 ; 5 7 2 CCT ¼ jz anterior ðx anterior apex ; y anterior apex Þ − z posterior ðx anterior apex ; y anterior apex Þj: The refractive powers in the natural environment (interface air/cornea/aqueous humor) were calculated in both meridians K and K s for the anterior surface using Javal's keratometric approximation K ¼ n k −1 R anterior × 1000 with n k ¼ 1.3375. The posterior surface power was calculated by means of the surface power equation F ¼ n a −n c R posterior × 1000 with n a ¼ 1.336 (aqueous, inner medium) and n c ¼ 1.376 (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 TNP ¼ n c −1 R anterior × 1000 þ n a −n c R posterior × 1000 using the mean radius of curvature for each surface. Corneal astigmatism was calculated by the absolute difference of K or F.

Results
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  Table 1 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 AE standard deviation. Anterior keratometric powers are given with respect to Javal's keratometric index.  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
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 R f;anterior ¼ 7.890 mm, R s;anterior ¼ 7.705 mm, R f;posterior ¼ 6.538 mm, and R s;posterior ¼ 6.354 mm correspond to average horizontal and vertical corneal curvatures found in previous studies in vivo. 33,37,38 The mean CCT of 616 μm is slightly higher than the average value found in normal human corneas (534 μm). 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 (539.9 μm) 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 R f ¼ 9.214 mm and R s ¼ 6.099 μm 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 (κ f ¼ −1.074 and κ s ¼ −0.844). Although the posterior surface is within the average range, this cornea exhibits the highest anterior refractive power (D f ¼ 61.65 dpt).
With cornea 3, all corneal parameters except for the CCT (844.28 μm) were close to the median values. However, the ratios of the anterior to posterior curvature radius (0.794 and Table 2 Corneal parameters of the suspected post-LASIK cornea (cornea 1), a cornea with extreme curvature radii (cornea 2), and the thickest cornea (cornea 3).  Table 3 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. 0.784 for the flatter and the steeper meridian, respectively) are lower than the average ratio (0.82).

Fitting Procedure
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, 29 AE 20 iterations were necessary until the fit converged (minimum: 6, maximum: 95). The average RMSE was 15.13 μm AE 13.30 μm (minimum: 3.72 μm and maximum: 69.32 μm), and the highest values occurred when the number of iterations culminated.
Calculating the biconic surface with the average anterior radii of curvature (R f ¼ 8.030 mm and R s ¼ 7.624 mm) and conic constants (κ f ¼ −0.286 and κ s ¼ −0.252), a difference of Δκ ¼ −0.1 manifests in the peripheral regions of the measurement area (−2.83 mm and 2.83 mm) in an elevation difference of 1.7 μm. Considering the depth resolution of 3.893 μm and the mean RMSE of the fitting procedure, this calculation of R 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, circularshaped 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 surfacefitting 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 <3%.
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. [41][42][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. [48][49][50] 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 device 47,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.

Disclosures
The authors have no financial interest to declare in this manuscript.