Translator Disclaimer
3 December 2019 Corneal thickness measurement by secondary speckle tracking and image processing using machine-learning algorithms
Author Affiliations +

Corneal thickness (CoT) is an important tool in the evaluation process for several disorders and in the assessment of intraocular pressure. We present a method enabling high-precision measurement of CoT based on secondary speckle tracking and processing of the information by machine-learning (ML) algorithms. The proposed configuration includes capturing by fast camera the laser beam speckle patterns backscattered from the corneal–scleral border, followed by ML processing of the image. The technique was tested on a series of phantoms having different thicknesses as well as in clinical trials on human eyes. The results show high accuracy in determination of eye CoT, and implementation is speedy in comparison with other known measurement methods.



Corneal thickness (CoT) measurement has gained relevance in recent years. Evaluation of CoT is performed in a wide range of disorders and procedures, such as refractive surgery, assessment of ocular hypertension, corneal disease, and glaucoma.1,2 It is also an increasingly important measurement for precise evaluation of patients’ intraocular pressure (IOP).3 Ultrasound pachymetry (UP) is the most commonly used method to measure CoT. Although UP is accurate, easy to use, and relatively inexpensive, an invasive procedure necessitates a cooperative patient and use of anesthetic eye drops, and it may also be a source of infection.4 An alternative way of measuring CoT is by noncontact imaging techniques. One such method is optical coherence tomography (OCT), which is similar to ultrasonography imaging but uses light instead of acoustic waves. Two-dimensional cross-sectional microscopic images of tissue are created by OCT from multiple scans of backscatter light versus depth.59 Another approach to CoT measurement is possible by the use of a rotating Scheimpflug camera, producing a 3-D scan of the anterior eye segment.10,11 The use of a specular microscope to photograph and examine human endothelial cells is also possible in CoT evaluation. This method involves applanation of the cornea with a dipping cone mounted on the objective.1214 An additional noncontact way to measure CoT is by optical low-coherence reflectometry (OLCR). This allows measurement of the amplitude and relative phase of backscattered or reflected light, although the presence of high frequency noise affects the measurements.15 OLCR consists of a Michelson optical interferometer, a light source, a beam splitter, and a light detector.16 Noncontact methods do not require topical anesthesia and the procedures are relatively simple; however, the results differ substantially from those of UP1721 and the reduced precision of CoT measurement could also affect IOP evaluation.

To measure IOP, opthalmologists often use applanation tonometers. Various studies show that the precision of applanation tonometers depends on CoT.2227 Normally, thinner corneas lead to lower IOP readings and thicker corneas to higher ones.3 Most of the current IOP measuring devices, such as the Goldmann applanation tonometer (GAT), the gold standard in IOP examination, enable the reading of IOP but lack the ability to measure CoT simultaneously, and have the additional drawback in needing two different devices to measure CoT and IOP. Combined measuring devices, such as noncontact tonometers that measure both IOP and CoT, have a lower precision level than that obtained by GAT.28

The current work presents machine-learning (ML) networks able to learn and classify the laser-based secondary speckle patterns (SSP) reflected from the tissues. The method allows evaluation of peripheral CoT with high accuracy and faster than obtainable by traditional methods. Recently, we used reflected SSP technology for semicontinuous noninvasive IOP measurement.29 It was noted that the measurement of corneal or scleral thickness is required to introduce corrections for noncontact IOP measurements. In the current research, CoT was measured for the first time using a novel technique based on speckle tracking and processing of information by ML algorithms. In this research, we combined CoT and IOP measurements on the same device using speckle-based technology. It allows including CoT for correction and accurate semicontinuous, noninvasive measurement of IOP.




Background of the Proposed Technology

The proposed technology contains further development for IOP evaluation, which is based on the tracking of tissue back reflected speckle patterns. A laser beam illuminates the tissue surface and the back scattered light is captured by a fast camera. When wavefronts are back reflected from a rough surface and are summed up on the detector plane, a random intensity pattern appears due to the interference phenomenon. Those patterns are known as secondary or subjective speckles. In our experiment, the diverse SSPs reflected from the series of phantoms having different thicknesses were captured by a fast imaging camera.

The underlying principle of the proposed method is that the tilting movement can be detect using only a laser and a defocused fast camera in order to transform this movement into a transversal shift of the speckle patterns.30,31 By doing so, the tissue’s vibrations cause the speckle patterns to shift in time, with the shift being proportional to the extent of tilting occurring in the illuminated surface.

The camera’s requirement for focal length F is30

Eq. (1)

where Δx is the pixel size in the detector, λ is the laser wavelength, Z2 is the distance between the illuminated surface and the secondary speckle (due to defocusing), Z3 is the distance between the secondary speckle and the imaging module, D is the diameter of the laser spot, and K is the minimum number of pixels, in which it is assumed that every speckle in this plane will be seen.

The far-field approximation needs to be kept by Z2:

Eq. (2)


The assumption behind our technique is that increase of tissue thickness, contributes the light scattering from an object, and consequently causes reduction of the average speckle size throw the image (see Fig. 1 for illustration). Furthermore, the overall speckle image undergoes several changes. The ML algorithms can find those features.

Fig. 1

Illustration of the proposed technique.


The reason that the size of the defocused speckles, being collected by the camera, depends on the amount of scattering of the inspected medium, is related to the fact that larger scattering coefficients increase the average path of light (diffusion of light) and thus generate on the surface of the scattering medium a larger spot of diffused light. In the far-field plane (defocused speckles), the spot of light on the surface of the medium sets the size of the collected speckles (larger spot yields smaller defocused speckles):

Eq. (3)

Thus the size of the defocused speckles is inversely dependent on the size of the diffused spot of light generated on the surface of the scattering medium.

This assumption is also supported in the scientific literature,32,33 for instance by the model presented in Ref. 32 where the scattering coefficient of the medium sets an effective location of a master and slave light sources that properly fulfill the light diffusion equation where the position of the master source from the surface of the scattering medium is set by

Eq. (4)

and of the position of the slave source is set by

Eq. (5)


Eq. (6)

g is the scattering anisotropy of the inspected medium, μs is the scattering coefficient, and n is an open parameter depending on the modeling as presented in Ref. 32. So the distance of the effective light source depends on the scattering properties of the inspected medium. This means that the size of the speckles varies since their dimensions depend on the effective distance as presented in Eq. (3).

This assumption of the dependence between the speckle size and the amount of scattering is also experimentally supported by Ref. 33 (relevant measurements are presented in Fig. S6 of the reference).


Speckle Image Transformation Tracking

The variation in speckle pattern images captured by fast camera, being related to the illuminated surface structure or thickness, can be analyzed by several methods, including cross-correlation or ML algorithms, as described below.


Cross correlation

The cross-correlation method was used to calculate correlation across the different axes of the image, given that two speckle pattern images differ only by an unknown shift. The cross-correlation function can be used to determine the image shift along the same axis. When the images match, the value of the cross-correlation function is maximized.

Autocorrelation is the cross correlation of a signal with itself. In order to calculate the average feature size, a normalized autocorrelation is defined34

Eq. (7)

where I represent the speckle intensity and u is the displacement.

The width of the autocorrelation peak under A=0.5 was found to be related to the average speckle size (in pixels).34

We tried cross correlation along the X axis for CoT partition and identified a correlation of 4 to 6 pixels only, which was found to be insufficient for successful CoT separation.


Machine-learning algorithms

ML35 refers to computer systems used to effectively perform a given task without using specific instructions while relying only on captured data (which in our case are the speckle image). ML algorithms create a mathematical model of sample data, known as “training data,” in order to make predictions or decisions on other types of data without being explicitly programmed to perform the assignment. We used two ML network architectures to recognize different CoT variations.

Convolution neural network with classification/regression layer architecture

For the first time, a convolution neural network (CNN)36 with classification/regression architecture was used to analyze nonlinear image features and establish an accurate predictor based on CNN training using SSP images.

The functions described below were applied in the network.

  • Rectified linear units.

Rectified linear units (ReLU) were applied as nonlinear functions for network training. These grant flexibility during training and show a good fit to image data. Activation is simply defined as threshold at zero.

  • Fully connected layer.

The fully connected layer function is applied at the end of the network before classification to reduce the output size to the number of classifications while preserving all the neurons data.37 It is the most accepted layer type for common neural networks: neurons between two adjacent layers are fully connected pairwise, yet no connection exists in a single layer.

  • Softmax layer.

The Softmax function takes an un-normalized vector and normalizes it into a probability distribution. In neural networks, it is employed to map un-normalized output to a probability distribution over a predicted output class.

ML methods were used to process and classify speckle pattern images with a goal of determining the object thickness. The CNN architecture used can be seen in Fig. 2.

Fig. 2

Schematic CNN architecture scheme.



Materials and Methods


Experimental Setup

A diagram of the experimental setup for phantoms having tissue-like properties is presented in Fig. 3.

Fig. 3

Experimental setup for remote testing of tissue-like phantoms.


The system contains a laser and a fast camera for monitoring SSP reflected back from the object. The camera’s focal length was 55 mm, with an F number of 2.8, and the illuminating beam was 3-mm in diameter. In the tissue-like phantoms a HJ 532-nm green laser was used for illumination and the speckle patterns created were recorded by a PixeLINK high-sampling-rate digital camera. The camera was positioned 35 cm away from the tested object (phantoms/human eye) and has no position constraints due to the wide angle of speckle diffraction. The illumination beam was directed at the center of the liquid phantom, placed in a Petri dish. For testing human eyes, a JDS Uniphase CW 1550 nm WDM DFB laser was positioned directly opposite the eyeball at distance of 32 cm. The laser beam was fixed at a selected location on the cornea. The power of the infrared laser was limited to the safety range (750  μW) approved by European Standard EN 60825-1. The speckle patterns reflected back from the cornea were captured using an EHD Imaging InGaAs IK1112 digital camera.

The tested objects were illuminated for 3 s and the camera captured the speckle images with a frequency of 200 frames per second. A total of 600 random speckle images were captured for each CoT that was tested in order to have a representative database of pictures for postprocessing analysis. However, the duration of illumination and the number of images could be reduced without affecting the accuracy of measurement. Increasing the camera frame rate will reduce the measurement time dramatically.


Liquid Phantom Preparation

The current work assumes that tissue-like phantoms (with the same scattering coefficient) having different thicknesses could be identified by analyzing SSP reflected from an object illuminated by a laser beam.

At the first stage, we tried to prove that our ML network has a stable network architecture that can detect tissue-like phantom with different thicknesses.

It is well established that the phantoms and human cornea can differ in their optical properties, but if we have high accuracy throughout our ML test, we may conclude that this tool is suitable for estimation of CoT in vivo, which was proved later experimentally.

The technique was examined with different coherent light wavelengths, exposure times, and resolutions to determine which parameters influence the measurement precision. Liquid phantoms with different thicknesses were prepared and placed in a Petri dish (35-mm diameter) to simulate a tissue. The phantoms were prepared using varying concentrations [8% for 532 nm (Ref. 38) and 1.5% for 650 nm (Ref. 39)] of Intralipid (IL) (Lipofundin MCT/LCT 20%, B. Braun Melsungen AG, Germany), used as a scattering component. The IL concentration was calculated according to Cubeddu et al.40 Seventeen phantoms having thicknesses from 0.4 to 0.8 mm, with 0.025-mm increments, were created and the Petri dish was measured as a reference. The selected thickness range covers the entire CoT variation of the human eye. Three repetitions were conducted for each wavelength, exposure time, and image resolution. Six hundred frames were taken as a database for our ML algorithm. Finally, the results of the experiment were analyzed and compared for different resolutions, cases, wavelengths, and network architectures.

The experiments showed that the wavelength and pixel resolution have a diminutive effect on the phantom thickness (PhT) measurement, compared with intensity and speckle pattern size and variation having significant influence on the results.

The experimental setup structure is presented in Table 1. Different wavelengths, exposure times, and resolutions were examined on different PhTs, and the mean absolute error was evaluated.

Table 1

Tissue like phantoms experimental setup and measurement error.

Test no.Illumination wavelength (nm)PhT range (μm)Thickness increment (μm)Exposure time (ms)Resolution (pixel)Mean abs. measurement error (μm)
1532450 to 800500.35128×1280
2400 to 775250.3532×327.55
3400 to 800250.664×649.98
4400 to 800250.664×6410.64
5400 to 800250.232×3212.75
6400 to 800250.232×3213.9
7400 to 800250.232×3215.07
8450 to 800500.3532×3218.75
9650400 to 800250.664×6418.28
10400 to 800250.664×6420.07
11400 to 800250.232×3220.69
12400 to 800250.232×3225.19

In order to show that our network is reliable and has good generalization to predict CoT, we let the network forecast several similar tests combined together. The predicted results are shown in Table 2.

Table 2

Combined experiments setup summary.

Test numberCombined experiments (from Table 1)Mean abs. measurement error (μm)
13 to 413.52
25 to 719.09
39 to 1025.96
411 to 1226.46


Human Eye Measurements

Measurements of CoT were performed in the Bar-Ilan Optometry Department. The experimental protocol was approved by the internal board of ethics committee of Bar-Ilan University. The applied laser system has a laser safety certificate for biomedical measurements on human eyes.


Reference eye cornea measurement by pachymetry

Pachymetry measurements of CoT were performed by a qualified technician with a view to obtaining reference data. Our ML methods were later compared with the pachymetry results in order to evaluate the accuracy of the speckle-based CoT measurement under eye safe laser illumination. The device used for the measurements was TOMEY TMS-5. The device accuracy41 according to the specification is 20  μm by the spherical accuracy factor.

10 eyes on subjects aged 23 to 72 were tested for CoT assessment. The left and right eye in each individual was tested three times and the median result was taken into account.

The CoT in the center and nasal side of the eye was analyzed using standard pachymetry equipment, as presented in Fig. 4.

Fig. 4

Result of eye CoT measurement by pachymetry.



In vivo human eye measurements

The 10 human eyes tested for CoT evaluation by pachymetry were also tested by our speckle-based method. The monitoring device was positioned at a distance of 32 cm from the eye. A diagram of the experimental setup for the measurement is presented in Fig. 5.

Fig. 5

Setup for measurement of corneal thick-CoT in the human eye.


As shown in the previous phantom experiments (Sec. 3.2), the method has very low wavelength dependence and therefore the tested eyes were illuminated by a JDS Uniphase CW 1550 nm WDM DFB eye safe laser. An EHD Imaging InGaAs IK1112 digital camera in defocused mode was used to capture reflections from the cornea. Focusing was performed on a focal plane that fulfilled the far-field conditions of diffraction with respect to the reflecting surface, which in the current case was the nasal side of the eye at the cornea-sclera limit. The pupil of the eye absorbs most of the energy so we could not test the entire cornea surface due to safety regulations limiting laser power dispassion to be below 1 mW.

Each frame of the camera’s output represents an SSP. The speckle patterns were analyzed with our ML algorithm using MATLAB (The MathWorks Inc., Massachusetts) to evaluate CoT.

The results were compared with those of the pachymetry measurements.




Corneal Thickness Evaluation by CNN with Classification Layer

A classification network containing speckle patterns recording was developed as a first step toward evaluation of the CoT by ML.

Splitting videos captured from the illuminated phantoms were subdivided into 18 categories in the range of 400 to 800  μm (covering the human CoT range) with increments of 25  μm. A green laser with a wavelength of 532 nm and 32×32  pixels resolution camera were used for the experiment. All combinations of the three described types of network architecture and two training options were tested after being found suitable for the task. We trained and evaluated CNNs using the different architectures described above. CNN classifiers were tested to generate probability scores. The network contained the following sections: imageInputLayer, convolution2dLayer, batchNormalizationLayer, ReLU layer, and maxPooling2dLayer, found to be the best for training (Fig. 6).

Fig. 6

Accuracy of network versus iteration while training.


The blue curve represents the training data fit, and the black curve represents the test data fit. The result shows 70.83% accuracy on test data (20% of data was used in the test). The maximum classification error of the network on test data was 75  μm, and the mean absolute error along the entire prediction vector was 8.125  μm.

We then evaluated the confusion matrix of test labels in comparison with the predictions on test data of our network. The results show a strong true label prediction match over different classes (strong confusion matrix diagonal), indicating good classification results of the test data.

As a result of our high-quality developed network, a very low absolute mean error was obtained on the test data—below 10  μm, matching the current standard equipment error.


Corneal Thickness Evaluation by CNN with Regression Layer

To analyze the nonlinear image features and to establish the most accurate predictor based on CNN network architecture, the regression network was applied. The regression layer could find the best fit to the nonlinear data while training. The experimental data were collected in a tensor and vector containing matching labels that describe the origin of the data. The tensor was constructed from the different videos recorded for different CoTs with one color channel (gray level image). The data were shuffled randomly so that the network would not “desire” obtaining one thickness from the others. The experimental data were split: 80% for training and 20% for test/validation purposes to examine if the CNN model with a total of 2160 random speckle pattern images were predicted on the test, which makes the prediction very reliable.

The images ware normalized to the range (0, 1) for fast and accurate training. Normalization helps the network to train better on different features with a different scale of values because the steps are constant across the training on normalized values. The data were uniformly distributed across the training data after shuffling randomly. The image input was constructed to have 32×32  pixels with one color channel (the speckle is gray level-based because the laser is in one wavelength only).

The following network architecture was used for an input gray image of 128×128  pixels.

  • 1. Five convolution layers were used for the network with filters of size 3×3 (pixels) convolving over the image to capture small features.

  • 2. After each convolution layer, a batch normalization layer was used to speed up the training process (as explained regarding the gradient).

  • 3. The ReLU nonlinear function was used to learn the nonlinear features of the image and obtain the best fit of data.

  • 4. The network size was reduced by max pooling, which takes a 2×2  pixels window from the image and averages the pixels, decreasing the network by a factor of 2 in each dimension.

  • 5. A fully connected layer was used to connect all output neurons from the previous layer to one output and attach it to the regression layer, which is optimized to reduce the root-mean-square error (RMSE) of the network by optimizing the network weights across the training data.

The optimizer used for the regression layer was stochastic gradient descent and the network was trained for 30 epochs, with an initial learning rate of 1×103 and a learning rate drop factor of 0.1 after 20 epochs.

The validation data were tested along the process and showed a small training-test gap, reflecting a low generalization error in the model. The training results with regression are presented in Fig. 7.

Fig. 7

Training results of phantoms regression network. Blue: training data fit. Black: test data fit.


The training results of the regression network show a fit error as low as 0.04, with a low generalization gap between training and test results, indicating a high-quality network. The mean absolute prediction error obtained from the test results (20% of data, 2160 images) after 30 epochs was 24.740  μm, with an RMSE of 0.04.

The test data accuracy for different threshold values is presented in Fig. 8.

Fig. 8

Test data accuracy as a function of corneal thick-CoT threshold.


In 99.21% of the cases, the network predicts the correct value 100  μm absolute error from the real thickness of the tested cornea. This shows very high accuracy with the 2160 images tested after network training. The mean predictor error is only 24.74  μm and this can be improved by higher division of the CoT range in the experiment.

To make sure that the network is reliable and has good generalization of CoT, we predicted the data from an additional experiment, in which the network had not been trained. This means that the network was trained in the first data set and tested in the additional data set obtained separately. The results are presented in Table 3.

Table 3

Results of training in additional data set.

ExperimentAbsolute mean error (μm)RMSETest accuracy (%)

Table 3 shows high-accuracy close to the original test data, low RMSE fit error, and good mean prediction error on test results when using the network on new experimental data. This demonstrates good generalization of the network, with reliability for future experiments and predictions.


Measurement of Human Eye Corneal Thickness


Pachymetry eye thickness measurements

We used pachymetry as reference measurements in order to compare the results with the speckle-based technique. The eye thickness measurement was repeated three times and the mean value, variance, and median were extracted. The test results for the 10 eyes (nasal side) are presented in Table 4.

Table 4

CoT pachymetry measurements on the 10 human eyes, μm. Each eye measured three times and the mean value, variance, and median were calculated.

Eye no.Measurement 1Measurement 2Measurement 3Mean valueVarianceMedian

We used pachymetry as a reference to our measurements while comparing the results with our speckle-based technique. The center of the cornea has higher absorption than the nasal side of the eye and we could not use a laser safe tool to test the center of the cornea. Thus we used only the nasal side of the subject’s eye. The eye thickness measurement was repeated three times and the mean value, variance, and median were extracted. The test results for 10 eyes (nasal side) are presented in Table 4.

Table 4 shows that the pachymetry device has a variance of 117.98  μm, with a low error rate of 15.5  μm (mean error of measurement between the minimum and maximum values from the three measurements conducted). This error can occur due to the human factor in measuring the CoT or due to the equipment itself, both being critical factors in the measurement, with no way to distinguish between the errors. The mean and median are equally good in estimating the experimental error.


In vivo speckle-based corneal thickness testing with machine-learning data processing

All subjects tested by the pachymetry method were tested with the speckle-based method and ML algorithm. The secondary speckle images were recorded using the setup as explained. In order to neglect eye motion factor, training was conducted on the delta between consecutive images. The delta also helps to focus only on the speckle pattern inside the image and eliminate any other constant objects that can affect the training process. Training only on the speckle pattern without the background prevents the network from converging into the object background instead of the speckle pattern in which we are interested.

We took five eyes as training data with an equal number of images (600 frames for each subject). An equal distribution of training data along the different CoTs ensures equal training for each CoT without preferring one thickness over others.

The training network results are presented in Fig. 9, showing the RMSE fit of the network while training versus the number of iterations after 10 epochs of training on the entire dataset. A low RMSE gives a better fit for data.

Fig. 9

Training results of human eyes regression network. Blue line: training data. Black line: validation data.


The optimized RMSE value of 0.031 was found after 10 epochs of training, and the results of training and test data were good.

In order to test the quality of our model, a second test was performed, in which a prediction was made on two new eyes that had not been used for training (which the network was not familiar with). Those eyes were tested with the same network architecture trained in the previous test. For reference purposes, the CoT was evaluated by pachymetry, giving values of 716 and 719 nm. The testing part contained 1200 images with different random speckle patterns. Considering the generalization gap and the quality of prediction, the results show high accuracy of the method. The result showed the accuracy of the network along all test images while using different permissible thresholds of the allowed error. For example, if the threshold is set to 50  μm and the regression network predicts a CoT with a ±50-μm error, it is considered a good prediction and anything above is considered a false prediction. If 50% of the test data is below this threshold, the prediction accuracy is 50%. For a threshold of 75  μm, the accuracy is already above 95% of data inside this boundary. For a 100-μm threshold, we achieved 100% accuracy for all our data. A summary of the regression network on human eyes is presented in Table 5.

Table 5

Validation results of network regression on human eyes.

Test numberAbsolute mean error (μm)RMSE (root mean squared error)Accuracy for 50 nm threshold error (%)

A comparison with known measurement tools such as OLCR, UP, specular microscopy (SM), and Pentacam (Scheimpflug method) can be seen in Table 6.

Table 6

Comparison between our technique and known methods.

CoT measurement methodsOLCRUPSMPentacamOur method SSP-based ML analysis
Repeatability (%)1.513.463.144.23N.A
Measurement time (s)18.55.613.545.7<1  s can be reduced dramatically with higher frame rate and GPU unit
Bland–Altman plot fitGood fit for UP, Pentacam, and OLCRGood fit for UP, Pentacam, and OLCRPoor fit to other methodsGood fit for UP, Pentacam, and OLCRGood fit for TOMEY TMS-5 26  μm±20  μm

The time of our speckle-based measurement with ML was 1/200  s per frame (it can be reduced dramatically) and shows a good fit for TOMEY TMS-5. We found a mean prediction absolute difference of 26  μm with a small ML eye database as described, and the device accuracy of TOMEY was found to be 20  μm.

Thus it can be assumed that the total true network error can be somewhere around 46  μm (measurement error of 26  μm at most, depending on the 20-μm optical lab accuracy as described).

To conclude, we developed a very accurate measurement tool for CoT estimation based on the ML network.



CoT measurement is an important factor for accurate evaluation of IOP. A technique using a recording of SSPs and evaluation by a regression network based on CNN architecture was introduced in order to estimate CoT on the scleral border. The technique was also shown to apply successfully on human eye as an example of biological tissue. It was found that CoT evaluation could be treated as a classification problem. The different classes have different features as showed by our ML network and by the used scanning filters. The problem of CoT estimation can also be treated as a regression problem, giving accurate estimation. The ML-based method was demonstrated to be reliable in terms of accuracy (26  μm of mean fit error) being on the level of pachymetry precision. The measurement is simple and fast, with the test taking only a few seconds. The testing period can be shortened by selecting a camera with higher frequency. It was shown that speckle-based measurement of IOP and CoT could be done simultaneously on the same device, whereas the CoT value could be used as a correction factor for IOP evaluation.


The authors declare that there are no conflicts of interest related to the research presented in this article.


The authors are grateful to Inbar Yariv and Hamootal Duadi, Engineering Department, Bar-Ilan University, for their assistance in the phantom preparation process.



M. Böhnke et al., “Continuous non-contact corneal pachymetry with a high speed reflectometer,” J. Refractive Surg., 14 (2), 140 –146 (1998). Google Scholar


I. A. Chaudhry, “Measurement of central corneal thickness in health and disease,” Saudi J. Ophthalmol., 23 (3–4), 179 –180 (2009). Google Scholar


R. P. Copt, R. Thomas and A. Mermoud, “Corneal thickness in ocular hypertension, primary open-angle glaucoma, and normal tension glaucoma,” Arch. Ophthalmol., 117 14 –16 (1999). AROPAW 0003-9950 Google Scholar


G. W. Belovay and I. Goldberg, “The thick and thin of the central corneal thickness in glaucoma,” Eye, 32 915 –923 (2018). 12ZYAS 0950-222X Google Scholar


C. J. Pavlin, M. D. Sherar and F. S. Foster, “Subsurface ultrasound microscopic imaging of the intact eye,” Ophthalmology, 97 244 –250 (1990). OPANEW 0743-751X Google Scholar


R. C. Youngquist, S. Carr and D. E. N. Davies, “Optical coherence-domain reflectometry: a new optical evaluation technique,” Opt. Lett., 12 158 –160 (1987). OPLEDP 0146-9592 Google Scholar


D. Huang et al., “Micron-resolution ranging of cornea anterior chamber by optical reflectometry,” Lasers Surg. Med., 11 419 –425 (1991). LSMEDI 0196-8092 Google Scholar


K. Takada et al., “New measurement system for fault location in optical waveguide devices based on an interferometric technique,” Appl. Opt., 26 (9), 1603 –1606 (1987). APOPAI 0003-6935 Google Scholar


E. A. Swanson et al., “High-speed optical coherence domain reflectometry,” Opt. Lett., 17 (2), 151 –153 (1992). OPLEDP 0146-9592 Google Scholar


B. Lackner et al., “Repeatability and reproducibility of central corneal thickness measurement with Pentacam, Orbscan, and ultrasound,” Optom. Vision Sci., 82 892 –899 (2005). OVSCET 1040-5488 Google Scholar


T. M. Rabsilber, R. Khoramnia and G. U. Auffarth, “Anterior chamber measurements using Pentacam rotating Scheimpflug camera,” J. Cataract Refract. Surg., 32 456 –459 (2006). Google Scholar


P. Maurice, “Self-closing valve device for implantation in the human body,” U.S. Patent 3,402,710, Hydra Power Corp. (1968).


W. M. Bourne and H. E. Kaufman, “Specular microscopy of human corneal endothelium in vivo,” Am. J. Ophthalmol., 81 319 –323 (1976). AJOPAA 0002-9394 Google Scholar


R. A. Laing, M. M. Sandstrom and H. M. Leibowitz, “In vivo photomicrography of the corneal endothelium,” Arch. Ophthalmol., 93 143 –145 (1975). AROPAW 0003-9950 Google Scholar


D. P. Davé and T. E. Milner, “Optical low-coherence reflectometer for differential phase measurement,” Opt. Lett., 25 227 –229 (2000). OPLEDP 0146-9592 Google Scholar


B. R. Masters, “Early development of optical low-coherence reflectometry and some recent biomedical applications,” J. Biomed. Opt., 4 (2), 236 –248 (1999). JBOPFO 1083-3668 Google Scholar


S. M. Nam et al., “Accuracy of RTVue optical coherence tomography, Pentacam, and ultrasonic pachymetry for the measurement of central corneal thickness,” Ophthalmology, 117 (11), 2096 –2103 (2010). OPANEW 0743-751X Google Scholar


F. Bao et al., “Comparison and evaluation of central corneal thickness using 2 new noncontact specular microscopes and conventional pachymetry devices,” Cornea, 33 576 –581 (2014). CORNDB 0277-3740 Google Scholar


L. Borrego-Sanz et al., “Agreement between central corneal thickness measured using Pentacam, ultrasound pachymetry, specular microscopy and optic biometer Lenstar LS 900 and the influence of intraocular pressure,” Ophthalmologica, 231 (4), 226 –235 (2014). OPHTAD 0030-3755 Google Scholar


J. J. Garcia-Medina et al., “Comparative study of central corneal thickness using Fourier-domain optical coherence tomography versus ultrasound pachymetry in primary open-angle glaucoma,” Cornea, 32 9 –13 (2013). CORNDB 0277-3740 Google Scholar


M. Ang et al., “Comparison of anterior segment optical tomography parameters measured using a semi-automatic software to standard clinical instruments,” PLoS One, 8 e65559 (2013). POLNCL 1932-6203 Google Scholar


N. Ehlers and F. K. Hansen, “Central corneal thickness in low-tension glaucoma,” Acta Ophthalmol., 52 740 –746 (1974). Google Scholar


N. Ehlers, F. K. Hansen and H. Aasved, “Biometric correlations of corneal thickness,” Acta Ophthalmol., 53 652 –659 (1975). Google Scholar


N. Ehlers, T. Bramsen and S. Sperling, “Applanation tonometry and central corneal thickness,” Acta Ophthalmol., 53 (1), 34 –43 (1975). Google Scholar


M. M. Whitacre, R. A. Stein and K. Hassanein, “The effect of corneal thickness on applanation tonometry,” Am. J. Ophthalmol., 115 592 –596 (1993). AJOPAA 0002-9394 Google Scholar


A. Yildiz and T. Yasar, “Comparison of Goldmann applanation, non-contact, dynamic contour and tonopen tonometry measurements in healthy and glaucomatous eyes, and effect of central corneal thickness on the measurement results,” Med. Glas., 15 (2), 152 –157 (2018). Google Scholar


A. Khatri et al., “Influence of central corneal thickness (CCT) on the intraocular pressure (IOP) measurements taken from Goldmann Applanation tonometer, tonopen, and airpuff tonometer,” Birat J. Health Sci., 3 532 –536 (2019). Google Scholar


S. Y. Hsu et al., “Comparisons of intraocular pressure measurements: Goldmann applanation tonometry, noncontact tonometry, Tono-Pen tonometry, and dynamic contour tonometry,” Eye, 23 1582 –1588 (2009). 12ZYAS 0950-222X Google Scholar


A. Bennett et al., “Intraocular pressure remote photonic biomonitoring based on temporally encoded external sound wave stimulation,” J. Biomed. Opt., 23 117001 (2018). JBOPFO 1083-3668 Google Scholar


Z. Zalevsky et al., “Simultaneous remote extraction of multiple speech sources and heart beats from secondary speckles pattern,” Opt. Express, 17 (24), 21566 –21580 (2009). OPEXFF 1094-4087 Google Scholar


Z. Zalevsky and J. Garcia, “Motion detection system and method,” U.S. Patent 8,638,991, Universtitat de Valencia and Bar Ilan University (2014).


D. Piao and S. Patel, “Simple empirical master-slave dual-source configuration within the diffusion approximation enhances modeling of spatially resolved diffuse reflectance at short-path and with low scattering from a semi-infinite homogeneous medium,” Appl. Opt., 56 (5), 1447 –1452 (2017). APOPAI 0003-6935 Google Scholar


A. Schwarz et al., “Coherent elastographic tomography via time multiplexing,” Quantum Electron. J., 49 (1), 35 –42 (2019). Google Scholar


H. Schreier, J. J. Orteu and M. A. Sutton, Image Correlation for Shape, Motion and Deformation Measurements: Basic Concepts, Theory and Applications, Springer-Verlag, Boston, Massachusetts (2009). Google Scholar


C. M. Bishop, Pattern Recognition and Machine Learning (Information Science and Statistics), Springer-Verlag, Berlin, Heidelberg (2006). Google Scholar


“Convolutional neural networks for visual recognition,” (2019) April 2019). Google Scholar


J. Long, E. Shelhamer and T. Darrell, “Fully convolutional networks for semantic segmentation,” in Proc. IEEE Conf. Comput. Vision and Pattern Recognit., 3431 –3440 (2015). Google Scholar


A. J. Costantino et al., “Determining the light scattering and absorption parameters from forward-directed flux measurements in cardiac tissue,” J. Biomed. Opt., 22 076009 (2017). JBOPFO 1083-3668 Google Scholar


S. L. Jacques, “Optical properties of biological tissues: a review,” Phys. Med. Biol., 58 R37 –R61 (2013). PHMBA7 0031-9155 Google Scholar


R. Cubeddu et al., “A solid tissue phantom for photon migration studies,” Phys. Med. Biol., 42 1971 –1979 (1997). PHMBA7 0031-9155 Google Scholar


Aviya Bennett received her BSc and MSc degrees in electrical engineering at Bar Ilan University (BIU) in 2005 and 2016, respectively. In the interim she worked at Intel Electronics Corporation. Currently, she is working toward her PhD studies under Professor Zeev Zalevsky’s supervision. She received the “Intel prize for Academic Excellence” in 2015. Her research interests include optical remote sensing, biomedical optics, and nanophotonics. She has published four refereed journal papers and two patents.

Elnatan Davidovitch received his BSc and MSc degrees in electrical engineering from Bar-Ilan University (BIU) at 2019 under the supervision of Professor Zeev Zalevsky. He worked at Intel Corporation for three years. Currently, he works at Orbotech Corporation, since 2017 as algorithm engineer in the field of machine learning and classical vision algorithms. His research studies and expertise include machine learning, deep learning, biomedical optics, and fluorescence microscopy.

Yafim Beiderman received his MSc degree in mechanics from the Polytechnic Institute, Odessa, Ukraine, his PhD from Agricultural Academy, Latvia. Since 2015, he has been with Bar-Ilan University, Israel. His research involves the development, testing, and evaluation of biomedical sensors, noninvasive laser speckle-based measurement of blood and intraocular pressure, breast cancer early detection, fiber based noncontact estimation of vital bio-signs, and alternative ways for transmitting information to the brains of blind people.

Sergey Agdarov received his MSc degree in mechanics and mathematics from Tajiki State University of Dushanbe, Tajikistan, in 1987. Since 2011, he has worked as a researcher in the faculty of engineering at Bar-Ilan University, Ramat-Gan, Israel. His research involves the development, testing, and evaluation of biomedical optical sensors, noninvasive measurement of IOP, breast cancer early detection, silicon-based mechanic-photonic wavelength conversion, and fiber sensors for noncontact estimation of vital bio-signs.

Yevgeny Beiderman received his BSc degree in mechanical engineering from Tel Aviv University in 1999, his MSc degree in intermediate studies in engineering, Tel Aviv University in 2002, and his PhD in mathematics in collaboration with of electro-optical track at the Engineering Faculty of Bar-Ilan University in 2010. He has published more than 20 scientific publications and registered more than 10 patents. He has received numerous prizes due to his research developed during his stay at Bar-Ilan University.

Avital Moshkovitz is a licensed optometrist and currently a PhD student of vision science in the lab of Uri Polat at Bar Ilan University. Her BSc and MSc degrees were in optometry and vision science, from Bar Ilan University.

Uri Polat neuroscientist specializing in brain plasticity, visual processing and optometry. He received his PhD in brain research from the Weizmann Institute of Science, Rehovot, Israel, and fellowship at the Smith-Kettlewell Eye Research Institute, San-Francisco, USA. He is serving as a reviewer in many scientific journals and is an editorial board member for Nature’s Scientific Reports. Currently he is the head of the School of Optometry and Vision Science at Bar-Ilan University.

Zeev Zalevsky received his BSc degree and direct PhD in electrical engineering from Tel-Aviv University in 1993 and 1996, respectively. Currently, he is a full professor in the faculty of engineering at Bar-Ilan University, Israel. He is the dean of engineering and a board member of the nanophotonic center there. His major fields of research are optical super-resolution, biomedical optics, nanophotonics, electro-optical devices, and microwave photonics. He has published more than 800 papers, 9 books, and has more than 100 patents.

© The Authors. Published by SPIE under a Creative Commons Attribution 4.0 Unported License. Distribution or reproduction of this work in whole or in part requires full attribution of the original publication, including its DOI.
Aviya Bennett, Elnatan Davidovitch, Yafim Beiderman, Sergey Agdarov, Yevgeny Beiderman, Avital Moshkovitz, Uri Polat, and Zeev Zalevsky "Corneal thickness measurement by secondary speckle tracking and image processing using machine-learning algorithms," Journal of Biomedical Optics 24(12), 126001 (3 December 2019).
Received: 26 June 2019; Accepted: 4 November 2019; Published: 3 December 2019

Back to Top