Expand this Topic clickable element to expand a topic
Skip to content
Optica Publishing Group

OCT based corneal densitometry: the confounding effect of epithelial speckle

Open Access Open Access

Abstract

Corneal densitometry is a clinically validated method for objectively assessing the transparency of stroma. The technique is currently dominated by Scheimpflug technology. Still, optical coherence tomography (OCT), in which examination of the statistical properties of corneal speckle is undertaken, has also been considered to assess corneal densitometry. In-vivo, the stroma is observed via the epithelium. However, the effect of this external layer on stromal densitometry has not been considered as yet. This study aims to quantify the influence of epithelium integrity on corneal OCT densitometry. OCT images from eleven freshly enucleated porcine eyes before and after epithelial debridement were used. OCT densitometry was investigated at different stromal depths using four metrics of speckle statistics. Results indicate that there exist statistically significant differences in speckle statistics for a given stromal depth depending on the presence or absence of the epithelium. The estimation error in speckle statistics can reach over 20% depending on the stromal depth. The anterior stroma densitometry values are the ones most affected by epithelial integrity. In conclusion, if OCT densitometry stromal parameters are to be considered in absolute terms, it is essential to consider the confounding effect of the epithelial layer in the analysis.

© 2023 Optica Publishing Group under the terms of the Optica Open Access Publishing Agreement

1. Introduction

Objectively assessing corneal transparency is essential to evaluate ocular health [1,2]. In a clinical environment, there is a growing demand for objective and robust estimators of corneal transparency, cumulatively termed corneal densitometry. It is often used for assessing pathological corneas in cases such as keratoconus [2,3] and Fuchs endothelial dystrophy [4] but also corneas after certain surgical procedures such as, for example, cataract and refractive surgeries [5,6] and Bowman layer transplantation [7]. Any influence of the epithelial layer on densitometry measures has not been considered as yet, despite being known that the integrity of epithelium is often compromised in those pathological cases [8,9]. Corneal densitometry estimators are based on evaluating the light backscattered from the cornea (in short, corneal backscatter). Whenever light from the source reaches the tissue to be imaged (here, the cornea), part of it is backscattered toward the detector to form an image. The final image, consequently, depends on how light travels inside the cornea and how much of it is backscattered. As the cornea is not an isolated structure, besides corneal tissue integrity, corneal backscatter might also be influenced by confounding factors such as eye biometry [10] or eye tilt [11].

Currently, clinical assessment of corneal densitometry is dominated by Scheimpflug imaging technology, whereas Optical Coherence Tomography (OCT), despite having a better axial resolution, has been given much less attention for performing this task [1214]. This imbalance is mostly driven by the commercial success of Scheimpflug-based corneal densitometry.

OCT is based on interferometry, which produces speckle. Even though the speckle has been traditionally treated as noise that degrades the quality of OCT image, there exists experimental evidence of speckle being signal-carrying [15]. Therefore, speckle statistics can be applied to infer indirectly about the microstructure of the tissue. In corneal research, previous works have applied the statistical modeling of corneal speckle to investigate corneal integrity with pathological conditions [16,17] and corneal tissue alterations as a consequence of hydration deprivation [12,18]. Corneal speckle statistics (a term reserved for OCT images) and corneal densitometry (a general term commonly used clinically) refer to the same concept, i.e., corneal backscatter. To date, mimicking Scheimpflug-based corneal densitometry, only the mean speckle intensity has been considered for OCT-based corneal densitometry [13]. However, as shown earlier [1619], other statistics of OCT speckle carry clinically useful information that could also be utilized for OCT densitometry.

The cornea is a multi-layered structure. The epithelium is the most anterior corneal layer. Underneath the epithelium rests the stroma, which makes up approximately 90% of the corneal thickness [20]. Corneal epithelium and stroma have different reflective properties. In regular corneal tomography, regardless of whether it is captured with Scheimpflug or OCT technology, the epithelium appears as a distinct thin layer on the top of the stroma. Hence, image pixels corresponding to the epithelium contain different signal information than that from pixels representing the stroma. In practice, corneal densitometry is limited to evaluating information from corneal stroma [1,13]. However, there exist validated computational models [21,22] and experimental evidence [23] on how the scattering properties are altered in multi-layered tissue structures in OCT images, suggesting that epithelium might play a substantial role in corneal densitometry estimation.

In a recent experiment based on corneal OCT imaging, performed on freshly enucleated porcine eyes, a statistically significant difference in stromal corneal speckle parameter (the contrast ratio) was found in eyes before and after epithelium removal, without modifying any experimental condition other than the removal of the epithelial layer [19]. This experimental result suggests that epithelium, a layer through which light necessarily travels to access the stroma, plays a key role in corneal backscatter.

The current work aims to quantify the influence of epithelium integrity on corneal densitometry estimates and assess whether this should be taken into consideration in clinical applications.

2. Methods

2.1 Dataset

Retrospective corneal OCT data of 11 porcine eyeballs, collected earlier in a study on crosslinking, were exported from OCT device (SOCT REVO, Optopol, Zawiercie, Poland). Eyeball preparation and data collection process are described in detail in [19]. To summarize, intact eyeballs with no damages, scars nor corneal edema were obtained from a certified abattoir. Each eyeball, regularly moistened with phosphate-buffered saline, was imaged three times before and three times immediately after the epithelial debridement (denoted hereinafter epi-on and epi-off mode, respectively). Porcine eyes were imaged ex-vivo, so there were no eye movements. Also, all samples were imaged within the depth-of-focus of the instrument. Each three OCT acquisitions were taken consecutively without re-aligning the instrument. The 5 mm single B-scan scanning protocol was used. It composes of 12,032 raw A-scans, from which an OCT image is formed with the fixed width of 1,538 pixels and a varying depth of approximately 733 pixels, corresponding to 5 mm and about 2 mm, respectively. Hence, the approximate resolution was 308 pixels/mm horizontally and 366 pixels/mm vertically.

2.2 Image processing

The main purpose was to compare the corneal OCT densitometry parameters for the same stromal regions in epi-on and epi-off mode acquisitions and to assess how the corneal epithelium layer influences OCT densitometry in the stroma. Firstly, the image segmentation algorithm described earlier was employed [18]. In order to analyze corneal backscatter for different stromal depths, a moving region of interest (ROI) with a width of 1 mm, depth of 100 µm, and a downward step of 10 µm, corresponding to 90% overlap, was applied (see Fig. 1). The maximum downward shift of ROI was adapted to the smallest central corneal thickness from the examined corneas (i.e., 526 µm) and fixed for all eyeballs. Such a setting resulted in 43 ROI positions in the range [0, 420] µm (further denoted as stromal depths). Sample axial position within the OCT image was noted by segmenting the anterior corneal surface and finding its highest position (i.e., that of the corneal apex). The epithelium layer in epi-on images was delineated manually using digital calipers and the first position of the ROI in epi-on samples was just under the measured epithelium. The mean $\pm$ standard deviation of measured epithelial thickness was 24 $\pm$ 2 pixels, corresponding to approximately 66.7 $\pm$ 5.9 µm. For epi-off samples, as there was no epithelium, the anterior edge of the first ROI was included in the anterior corneal edge (see Fig. 1). This allowed for the analysis of the same stromal regions for epi-on and epi-off modes. In the preprocessing stage, each B-scan underwent a normalization process [24] in which, first, the pixel values were divided by the maximum value of 255. Subsequently, the inverse-log transformation was applied. The transformed pixel values were further subjected to normalization to lead speckle values between 0 and 1.

 figure: Fig. 1.

Fig. 1. Illustrative corneal OCT images for a porcine eye before (left, epi-on) and after (right, epi-off) epithelial debridement. White frames indicate the moving region of interest (ROI). Both ROIs for epi-on and epi-off mark out the same stromal regions set further for the OCT densitometry analysis. The horizontal dashed lines represent the position of the sample with respect to the imaging system.

Download Full Size | PDF

2.3 Statistical analyses

To describe corneal OCT densitometry within each considered ROI, four measures of speckle statistics, further denoted as densitometry parameters, were calculated: Gamma distribution parameters (shape and scale), in a parametric approach, and mean ($\mu$) and contrast ratio (CR, the ratio of the standard deviation to the mean value), in a non-parametric approach [24]. The probability density function of the Gamma distribution can be described as:

$$f_X(x;\alpha,\beta) = \frac{1}{{\beta^\alpha}{\Gamma(\alpha)}}{x}^{\alpha-1}{e}^\frac{-x}{\beta}\quad\text{for}\quad x\ge 0,$$
where $\alpha >0$ is the shape parameter, $\beta >0$ is the scale parameter, and $x$ denotes the normalized image pixel intensity values in a specified ROI. For further analyses, for every eyeball, the median of the parameter estimates from the three epi-on and a median from three epi-off acquisitions were gathered. Additionally, the coefficient of variation (CV) was evaluated for every three measurements. All OCT densitometry parameters were tested for the null hypothesis of Gaussian distribution using the Shapiro-Wilk test. The level of significance was set to 0.05 for all tests employed in the study. The rejection of null hypotheses of normality was included in further statistical analysis. To ascertain whether there are differences in $\alpha >0$, $\beta >0$, $\mu$, and CR between epi-on and epi-off modes, Wilcoxon signed-rank test was used. Further, Spearman correlations were evaluated for epi-on versus epi-off parameters as well as partial correlations with controlling for stromal depths, subject number (eyeballs), difference in sample axial position between epi-on and epi-off modes, and the mean epithelial intensity in the epi-on mode.

To additionally assess the impact of the confounding effect of epithelial speckle on stroma densitometry, for each stromal depth the median error (expressed in percentages) was calculated as (e.g., for $\alpha$):

$$\left(1-m\left(\frac{\hat\alpha_{\text{epi-off}}}{\hat\alpha_{\text{epi-on}}}\right)\right)\times 100,$$
where $m$ denotes the sample median, whereas as $\hat \alpha _{\text {epi-off}}$ and $\hat \alpha _{\text {epi-on}}$ are the estimators of the shape parameter $\alpha$ for the epi-off and epi-on modes at the particular stromal depth, respectively. The median errors for the scale parameter $\beta$, the sample mean $\mu$, and CR were calculated in a similar fashion.

3. Results

The data for the analyses consisted of four densitometry parameters ($\alpha$, $\beta$, $\mu$, and CR) estimated for 11 eyeballs, for epi-on and epi-off modes, for 43 stromal depths. The mean coefficient of variation (CV) for three subsequent measurements, evaluated across eyeballs, modes and depths was less 5% (CV = 2.1%, 1.9%, 0.9% and 1.0% for $\alpha$, $\beta$, $\mu$ and CR, respectively). Figure 2 shows the boxplots for the considered parameters for the epi-on mode juxtaposed with those for the epi-off mode. Asterisks indicate stromal depths where the Wilcoxon sign rank test indicated statistically significant differences between epi-on and epi-off modes. It is evident that stromal OCT densitometry parameters change after the removal of the epithelium. These changes are not straightforward and depend on the sample depth. For the shape parameter $\alpha$, there are statistically significant differences at 270 and 280 µm and then at 340 µm and beyond. For the scale parameter $\beta$, there are statistically significant differences only at small sample depths (0–130 µm), where the stroma is close to the epithelial layer. For $\mu$, statistically significant differences are observed at both small (0–120 µm) and larger stromal depths (190–420 µm). For CR, the differences seem to be stronger for small sample depths (10–40 µm) and then for 360–420 µm.

 figure: Fig. 2.

Fig. 2. Boxplots of the four considered OCT densitometry parameters calculated for pixel intensities within the ROI: $\alpha$, $\beta$ (shape and scale of Gamma distribution), $\mu$ (mean) and CR (contrast ratio) calculated for different stromal depths. The value for each eyeball is the median of three repeated measurements. Each box is built of 11 values (for 11 porcine eyeballs). Circles denote outliers (Tukey’s fences with $k = 1.5$). Asterisks indicate stromal depths at which there are statistically significant differences between epi-on and epi-off mode. Epi-on and epi-off mode means the measurement of a porcine eyeball before and after epithelial debridement, respectively.

Download Full Size | PDF

Figure 3 shows the scatter plots for the median values of $\alpha$, $\beta$, $\mu$, and CR with the stromal depth indicated by the saturation of blue marks. The correlations of the four considered parameters between the epi-on and epi-off modes are apparent. At the same time, the scatter plots indicate that the densitometry parameters of the corneal stroma are confounded when viewed by the intact layer of epithelium. Table 1 shows the corresponding results of Spearman’s correlations and partial correlations of the four parameters in epi-on against epi-off modes. Despite incorporating stromal depths or the subject number as a controlling variable, the correlation is still maintained at the level of significance. Similarly, the difference in sample axial position that ranged between 0 and 17 pixels as well as the average intensity value in the epithelial region (mean $\pm$ standard deviation: $0.183\pm 0.029$) did not influence those correlations. For each of the parameters, the differences in correlation coefficients evaluated without a control variable (the first column in Table 1) and those with a given control variable (columns 2 to 5 of Table 1) were not substantial, indicating that all control variables are not confounding factors.

 figure: Fig. 3.

Fig. 3. Scatter plots of four considered OCT densitometry parameters calculated for pixel intensities within the ROI: $\alpha$, $\beta$ (shape and scale of Gamma distribution), $\mu$ (mean) and CR (contrast ratio) for epi-on and epi-off modes, where the stromal depth is indicated by the saturation of blue marks. Epi-on and epi-off mode means the measurement of a porcine eyeball before and after epithelial debridement, respectively. Each value is the median of three measurements taken for every eyeball. The dashed line indicates the 1:1 line.

Download Full Size | PDF

Tables Icon

Table 1. The results of Spearman’s correlations and partial correlations (Part. corr.) between the epi-on and epi-off modes, calculated for the four considered OCT densitometry parameters: $\alpha$, $\beta$, $\mu$ and CR. Partial correlations were performed with stromal depth (SD), the subject number (SN), difference in sample axial position between epi-on and epi-off mode (AP) and mean epithelium intensity (MEI) as the controlling variable (Ctrl). $r$ – correlation coefficient, $p$ – the p-value.a

Figure 4 shows the median error obtained for the four considered parameters for each stromal depth. An area of $\pm$5% error, assumed here as clinically negligible, has been highlighted in light blue. Of particular note is the existence of over 20% median error for the parameters $\beta$ and $\mu$.

 figure: Fig. 4.

Fig. 4. Median errors between epi-on and epi-off mode of four considered OCT densitometry parameters calculated for pixel intensities within the ROI: $\alpha$, $\beta$ (shape and scale of Gamma distribution), $\mu$ (mean) and CR (contrast ratio), showing the confounding effect of epithelial speckle when evaluating the densitometry on different stromal depths. Each value is the median of three measurements taken for every eyeball. An area of $\pm$5% error, assumed here as clinically negligible, has been highlighted in light blue.

Download Full Size | PDF

4. Discussion

OCT corneal densitometry is in its infancy, despite the technology having a substantial potential for this task when compared to currently dominated field of Scheimpflug-based densitometry. This study investigated the confounding effect of epithelial speckle on stromal densitometry. The structure of epithelium greatly differs from that of the stroma. Consequently, light backscatter behaves differently in these two layers. This is meaningful, because the stromal densitometry is seen through the epithelium in standard clinical conditions. The results from the current work show that the corneal densitometry results are affected by epithelium integrity and that the epithelial layer tends, depending on the parameter considered, to over- or underestimate the stromal readings (see Fig. 3).

Corneal OCT densitometry was investigated at different stromal depths using four speckle measures, and parametric and non-parametric approaches were considered. The parametric approach, in agreement with previous literature on statistical modeling of corneal OCT speckle [18,25], consisted of Gamma distribution with shape ($\alpha$) and scale ($\beta$) parameters. In the non-parametric approach, two parameters commonly used to describe OCT images were applied [19,24]: they are the mean pixel intensity ($\mu$) and the contrast ratio (CR). Independently of the metric used, statistically significant differences in corneal stroma densitometry were found depending on the presence (epi-on) or absence (epi-off) of the epithelium (see Fig. 2) and those differences were depth dependent.

If the epithelium were not affecting corneal densitometry readings, one would expect the lack of statistical significance in Fig. 2 and the points in the scatter plots of Fig. 3 to fall on the 1:1 line. However, this is not observed. For example, when it comes to parameters $\beta$ and $\mu$, the points corresponding to small stromal depths (i.e., close to the epithelium, darker points in Fig. 3) fall over the line 1:1. In other words, the values of $\beta$ and $\mu$ are overestimated in the epi-on mode. Contrarily, for $\alpha$, stromal readings are underestimated in deeper stromal layers (i.e., far from the epithelium, brighter points in Fig. 3).

Even though $\beta$ and $\mu$ represent parametric and non-parametric approaches, respectively, they have a similar mathematical origin because the mean value of the Gamma distribution is the product of $\alpha$ and $\beta$. Consequently, the information they carry is analogous (see Fig. 3 and 4). Similarly, $\text {CR}=1/\sqrt {\alpha }$, so $\alpha$ tends to be underestimated whereas CR is overestimated (Fig. 3). Hence, the parametric and non-parametric approaches carry different, yet complementary information.

There is a clinical significance of the obtained results as the examination of corneal speckle is of interest in studies involving the debridement of the epithelium [26] or the healing processes occurring in cornea [27,28]. In fact, any change in the optical properties of the epithelium (e.g., caused by epithelial damage, hypoxia, or changes in permeability) may also change the stromal speckle parameters. If the presence of the epithelium is not considered, then anterior corneal stromal readings could be over- or underestimated. The more significant bias occurs in the anterior stroma, likely because of the attenuation of the OCT signal at deeper layers. However, there are minor but statistically significant differences also observed in posterior stroma layers.

Beyond clinical practice, there is an emerging promising field of study based on depth-resolved corneal images with full-field OCT (FF-OCT), which allows histology-like analysis of stromal features [2931]. The results from the current study might also be of use for FF-OCT researchers, who are in the position to validate whether by increasing the axial resolution of OCT, similar confounding factor of epithelial layer on stromal densitometry parameters is evident.

The current study sheds a different light on the results of the experimental work performed earlier on crosslinking (CXL) using freshly enucleated porcine eyes [19]. In particular, this concerns the result of a CXL procedure corresponding to standard epi-off Dresden protocol, but without the exposure to UVA light. The statistically significant changes observed there in CR before the epithelium debridement and after the CXL (no UVA) procedure can be interpreted as changes between epi-on and epi-off modes, irrespective of CXL.

The study does not present strong limitations. Freshly enucleated porcine eyes with epithelial debridement were used in this research work. Nevertheless, it should be emphasized that the obtained results were achieved for the given set of porcine eyes and the particular type of imaging instrument with a given scanning protocol. The presented results highlight the problem of epithelium confounding densitometry results of the stroma. A similar effect of one tissue layer confounding the results of other tissue layer is expected to happen in the human eyes, but the extent of that effect is currently unknown. Another limitation of the study is that there were no standard Scheimpflug densitometry measures for comparison with those based on OCT images.

Concluding, in the case of estimating densitometry parameters of stroma, to avoid over- or underestimation, it is essential to consider the confounding effect of the epithelial layer in the analysis.

Acknowledgments

A part of this study was presented during ARVO Annual Meeting 2023 in April.

Disclosures

Authors declare no conflicts of interest.

Data availability

Data underlying the results presented in this paper are not publicly available at this time but may be obtained from the authors upon reasonable request.

References

1. A. M. Otri, U. Fares, M. A. Al-Aqaba, and H. S. Dua, “Corneal densitometry as an indicator of corneal health,” Ophthalmology 119(3), 501–508 (2012). [CrossRef]  

2. M. Koc, K. Tekin, M. I. Tekin, M. M. Uzel, P. Kosekahya, K. Ozulken, and P. Yilmazbas, “An early finding of keratoconus: increase in corneal densitometry,” Cornea 37(5), 580–586 (2018). [CrossRef]  

3. B. Lopes, I. Ramos, and R. Ambrósio Jr, “Corneal densitometry in keratoconus,” Cornea 33(12), 1282–1286 (2014). [CrossRef]  

4. M. Alnawaiseh, L. Zumhagen, G. Wirths, M. Eveslage, N. Eter, and A. Rosentreter, “Corneal densitometry, central corneal thickness, and corneal central-to-peripheral thickness ratio in patients with Fuchs endothelial dystrophy,” Cornea 35(3), 358–362 (2016). [CrossRef]  

5. F. Poyales, N. Garzón, J. Mendicute, I. Illarramendi, P. Caro, O. Já nez, F. Argüeso, and A. López, “Corneal densitometry after photorefractive keratectomy, laser-assisted in situ keratomileusis, and small-incision lenticule extraction,” Eye 31(12), 1647–1654 (2017). [CrossRef]  

6. S. Ishikawa, N. Kato, and M. Takeuchi, “Quantitative evaluation of corneal epithelial edema after cataract surgery using corneal densitometry: a prospective study,” BMC Ophthalmol. 18(1), 334 (2018). [CrossRef]  

7. S. Luceri, J. Parker, I. Dapena, L. Baydoun, S. Oellerich, K. van Dijk, and G. R. Melles, “Corneal densitometry and higher order aberrations after Bowman layer transplantation: 1-year results,” Cornea 35(7), 959–966 (2016). [CrossRef]  

8. D. Z. Reinstein, T. J. Archer, and M. Gobbe, “Corneal epithelial thickness profile in the diagnosis of keratoconus,” J. Refract. Surg. 25(7), 604–610 (2009). [CrossRef]  

9. M. Q. Salom ao, A. L. Hofling-Lima, B. T. Lopes, A. L. C. Canedo, D. G. Dawson, R. Carneiro-Freitas, and R. Ambrósio Jr, “Role of the corneal epithelium measurements in keratorefractive surgery,” Curr. Opin. Ophthalmol. 28(4), 326–336 (2017). [CrossRef]  

10. M. Miażdżyk, A. Consejo, and D. R. Iskander, “Assessing and compensating for the confounding factors in Scheimpflug-based corneal densitometry,” Biomed. Opt. Express 13(12), 6258–6272 (2022). [CrossRef]  

11. A. Consejo, M. Jiménez-García, J. J. Rozema, and A. Abass, “Influence of eye tilt on corneal densitometry,” Ophthalmic Physiol. Opt. 42(5), 1032–1037 (2022). [CrossRef]  

12. J. Wang, T. L. Simpson, and D. Fonn, “Objective measurements of corneal light-backscatter during corneal swelling, by optical coherence tomography,” Invest. Ophthalmol. Visual Sci. 45(10), 3493–3498 (2004). [CrossRef]  

13. X. Y. Wang, T. Q. Zhang, A. R. Rachwani, I. Blanco-Domínguez, C. Rocha de Lossada, A. M. Adán-Civiera, and J. Peraza-Nieves, “New algorithm for corneal densitometry assessment based on anterior segment optical coherence tomography,” Eye 36(8), 1675–1680 (2022). [CrossRef]  

14. M. Vilbert, R. Bocheux, H. Lama, C. Georgeon, V. Borderie, P. Pernot, K. Irsch, and K. Plamann, “Objective assessment of corneal transparency in the clinical setting: correction of acquisition artifacts in SD-OCT images and application to normal corneas,” in European Conference on Biomedical Optics, (Optica Publishing Group, 2021), pp. EW1C–5.

15. J. M. Schmitt, S. Xiang, and K. M. Yung, “Speckle in optical coherence tomography,” J. Biomed. Opt. 4(1), 95–105 (1999). [CrossRef]  

16. M. E. Danielewska, A. Antończyk, D. A. De Jesus, M. M. Rogala, A. Błońska, M. Ćwirko, Z. Kiełbowicz, and D. R. Iskander, “Corneal optical coherence tomography speckle in crosslinked and untreated rabbit eyes in response to elevated intraocular pressure,” Trans. Vis. Sci. Tech. 10(5), 2 (2021). [CrossRef]  

17. D. R. Iskander, M. A. Kostyszak, D. A. Jesus, M. Majewska, M. E. Danielewska, and P. Krzyzanowska-Berkowska, “Assessing corneal speckle in optical coherence tomography: a new look at glaucomatous eyes,” Optom. Vis. Sci. 97(2), 62–67 (2020). [CrossRef]  

18. D. A. Jesus and D. R. Iskander, “Assessment of corneal properties based on statistical modeling of OCT speckle,” Biomed. Opt. Express 8(1), 162–176 (2017). [CrossRef]  

19. M. E. Danielewska, M. A. Kostyszak, P. Sareło, M. Gąsior-Głogowska, M. Niemczyk, P. Prządka, A. Antończyk, Z. Kiełbowicz, and D. R. Iskander, “Indirectly assessing changes in corneal properties with OCT speckle after crosslinking in porcine eyes,” Exp. Eye Res. 219, 109051 (2022). [CrossRef]  

20. E. M. Espana and D. E. Birk, “Composition, structure and function of the corneal stroma,” Exp. Eye Res. 198, 108137 (2020). [CrossRef]  

21. L. Wang, S. L. Jacques, and L. Zheng, “MCML—Monte Carlo modeling of light transport in multi-layered tissues,” Comput. Methods Programs Biomed. 47(2), 131–146 (1995). [CrossRef]  

22. L. Thrane, M. H. Frosz, T. M. Jørgensen, A. Tycho, H. T. Yura, and P. E. Andersen, “Extraction of optical scattering parameters and attenuation compensation in optical coherence tomography images of multilayered tissue structures,” Opt. Lett. 29(14), 1641–1643 (2004). [CrossRef]  

23. L. Thrane, M. H. Frosz, D. Levitz, T. M. Jorgensen, C. B. Andersen, P. R. Hansen, J. Valanciunaite, J. Swartling, S. Andersson-Engels, A. Tycho, H. T. Yura, and P. E. Andersen, “Extraction of tissue optical properties from optical coherence tomography images for diagnostic purposes,” in Saratov Fall Meeting 2004: Optical Technologies in Biophysics and Medicine VI, vol. 5771 (SPIE, 2005), pp. 139–150.

24. M. Niemczyk and D. R. Iskander, “Statistical analysis of corneal OCT speckle: a non-parametric approach,” Biomed. Opt. Express 12(10), 6407–6421 (2021). [CrossRef]  

25. M. Niemczyk, M. E. Danielewska, M. A. Kostyszak, D. Lewandowski, and D. R. Iskander, “The effect of intraocular pressure elevation and related ocular biometry changes on corneal OCT speckle distribution in porcine eyes,” PLoS One 16(3), e0249213 (2021). [CrossRef]  

26. C. Wirbelauer, C. Scholz, H. Häberle, H. Laqua, and D. T. Pham, “Corneal optical coherence tomography before and after phototherapeutic keratectomy for recurrent epithelial erosions,” J. Cataract. & Refract. Surg. 28(9), 1629–1635 (2002). [CrossRef]  

27. A. Hixson, S. Blanc, and J. Sowka, “Monitoring keratitis resolution with optical coherence tomography,” Optom. Vis. Sci. 91(4), S40–S45 (2014). [CrossRef]  

28. B. J. Kaluzny, M. Szkulmowski, D. M. Bukowska, and M. Wojtkowski, “Spectral oct with speckle contrast reduction for evaluation of the healing process after PRK and transepithelial PRK,” Biomed. Opt. Express 5(4), 1089–1098 (2014). [CrossRef]  

29. E. Auksorius, D. Borycki, P. Stremplewski, K. Liżewski, S. Tomczewski, P. Niedźwiedziuk, B. L. Sikorski, and M. Wojtkowski, “In vivo imaging of the human cornea with high-speed and high-resolution Fourier-domain full-field optical coherence tomography,” Biomed. Opt. Express 11(5), 2849–2865 (2020). [CrossRef]  

30. K. Irsch, K. Grieve, M. Borderie, M. Vilbert, K. Plamann, D. Ghoubay, C. Georgeon, and V. Borderie, “Full-field optical coherence microscopy for histology-like analysis of stromal features in corneal grafts,” Journal of Visualized Experiments p. e57104 (2022).

31. R. Bocheux, P. Pernot, V. Borderie, K. Plamann, and K. Irsch, “Quantitative measures of corneal transparency, derived from objective analysis of depth-resolved corneal images, demonstrated with full-field optical coherence tomographic microscopy,” PLoS One 14(8), e0221707 (2019). [CrossRef]  

Data availability

Data underlying the results presented in this paper are not publicly available at this time but may be obtained from the authors upon reasonable request.

Cited By

Optica participates in Crossref's Cited-By Linking service. Citing articles from Optica Publishing Group journals and other participating publishers are listed here.

Alert me when this article is cited.


Figures (4)

Fig. 1.
Fig. 1. Illustrative corneal OCT images for a porcine eye before (left, epi-on) and after (right, epi-off) epithelial debridement. White frames indicate the moving region of interest (ROI). Both ROIs for epi-on and epi-off mark out the same stromal regions set further for the OCT densitometry analysis. The horizontal dashed lines represent the position of the sample with respect to the imaging system.
Fig. 2.
Fig. 2. Boxplots of the four considered OCT densitometry parameters calculated for pixel intensities within the ROI: $\alpha$, $\beta$ (shape and scale of Gamma distribution), $\mu$ (mean) and CR (contrast ratio) calculated for different stromal depths. The value for each eyeball is the median of three repeated measurements. Each box is built of 11 values (for 11 porcine eyeballs). Circles denote outliers (Tukey’s fences with $k = 1.5$). Asterisks indicate stromal depths at which there are statistically significant differences between epi-on and epi-off mode. Epi-on and epi-off mode means the measurement of a porcine eyeball before and after epithelial debridement, respectively.
Fig. 3.
Fig. 3. Scatter plots of four considered OCT densitometry parameters calculated for pixel intensities within the ROI: $\alpha$, $\beta$ (shape and scale of Gamma distribution), $\mu$ (mean) and CR (contrast ratio) for epi-on and epi-off modes, where the stromal depth is indicated by the saturation of blue marks. Epi-on and epi-off mode means the measurement of a porcine eyeball before and after epithelial debridement, respectively. Each value is the median of three measurements taken for every eyeball. The dashed line indicates the 1:1 line.
Fig. 4.
Fig. 4. Median errors between epi-on and epi-off mode of four considered OCT densitometry parameters calculated for pixel intensities within the ROI: $\alpha$, $\beta$ (shape and scale of Gamma distribution), $\mu$ (mean) and CR (contrast ratio), showing the confounding effect of epithelial speckle when evaluating the densitometry on different stromal depths. Each value is the median of three measurements taken for every eyeball. An area of $\pm$5% error, assumed here as clinically negligible, has been highlighted in light blue.

Tables (1)

Tables Icon

Table 1. The results of Spearman’s correlations and partial correlations (Part. corr.) between the epi-on and epi-off modes, calculated for the four considered OCT densitometry parameters: α , β , μ and CR. Partial correlations were performed with stromal depth (SD), the subject number (SN), difference in sample axial position between epi-on and epi-off mode (AP) and mean epithelium intensity (MEI) as the controlling variable (Ctrl). r – correlation coefficient, p – the p-value.a

Equations (2)

Equations on this page are rendered with MathJax. Learn more.

f X ( x ; α , β ) = 1 β α Γ ( α ) x α 1 e x β for x 0 ,
( 1 m ( α ^ epi-off α ^ epi-on ) ) × 100 ,
Select as filters


Select Topics Cancel
© Copyright 2024 | Optica Publishing Group. All rights reserved, including rights for text and data mining and training of artificial technologies or similar technologies.