On identification of human eye retinas by the covariance analysis of their digital images

Jonas Skeivalas; Eimuntas K. Parseliunas

doi:10.1117/1.OE.52.7.073106

15 July 2013 On identification of human eye retinas by the covariance analysis of their digital images

Jonas Skeivalas, Eimuntas K. Parseliunas

Author Affiliations +

Optical Engineering, Vol. 52, Issue 7, 073106 (July 2013). https://doi.org/10.1117/1.OE.52.7.073106

Abstract

An identification of human eye retinas by applying the covariance function and wavelet theory is presented. The estimations of the autocovariance functions of the two digital images or single image are calculated according to random functions, based on the vectors created from the digital image pixels. The estimations of the pixel’s vectors are calculated by spreading the pixel arrays of the digital images into single column. During the changing of the scale of the digital image, the wave frequencies of the colors of the single pixels are prekept, and the influence of the change of a scale in the procedures of the calculations of the covariance functions does not occur. The Red, Green, Blue (RGB) color model of the colors spectrum for the encoding of the digital images was applied. The influence of the RGB spectrum components and the tensor of colors on the estimations of the covariance functions were analyzed. The identity of the digital images is estimated by analysis of the changes of the correlation coefficient values in the corresponding diapason.

1. Introduction

A number of papers are contributing to retinal imaging and image analysis.¹^,² Retinal researchers and practicians make more and more wide use of the digital images, which causes the need to improve the image processing algorithms. In this paper we stress on an identification of the digital images by applying the photogrammetric methods and random functions theory. The spatial positions of the digital image pixels are defined by the spatial region of the frequencies of color waves, i.e., by radiometric level, applying the Red, Green, Blue (RGB) coding format of the colors spectrum. The theoretical model is based on the stationary random function taking into account that the errors of the color wave frequencies are random and that they are of the same accuracy, i.e., a mean of errors is $M Δ = const = 0$ and their dispersion is $D Δ = const$ , and that the covariance function of the digital images depends only on the difference of the arguments, i.e., on the pixel quantize interval. The estimations of the covariance function of the two digital images or the autocovariance function of the single image are calculated according to random functions, based on the vectors, created from the digital image pixels. The Fourier transformation³^,⁴ or wavelet theory⁵^–⁸ is used for processing the digital signals. The main goal of this article is to provide the opportunity for a continuous improvement of the core algorithms, driven by performance of the covariance analysis approach. Such algorithms could be used in various areas of research and practice, including public and clinical health, biomedicine, security systems, etc.

2. Covariance Model of the Light Frequencies Spectrum

Let us analyze the autocovariance and covariance theoretical models of the white light spectrum, which is combination of the colors of different frequencies. The various color systems [RGB, hue, saturation and value (HSV), luma information, in-phase, quadrature (YIQ), hue, lightness and saturation (HLS), etc.] are used to define the digital images. The RGB system is used widely. The image of this system can be easily redone to the image of the other color systems, in which one of the component will be light signal and the other components define the color. It is used in the visualization, because we know that the human eye retina is more sensitive to the changes of the light strength than to the changes of the color itself.

Let us apply a linear expression of the harmonic oscillations equation to describe the light frequencies spectrum:

Eq. (1)

a (t) = A \sin (φ + φ_{0}) = A (\sin \bar{φ} + δ φ \cos \bar{φ} - \frac{1}{2} δ^{2} φ \sin \bar{φ} + \dots),

where

A

is an amplitude of the frequencies,

φ = ω t

is a phase,

φ_{0}

is an initial phase,

ω = 2 π f

is a cycle frequency,

f

is a frequency,

t

is a time moment,

M φ = \bar{φ}

is a mean of a phase,

δ φ = φ - \bar{φ}

is the random error of a phase.

In real world and in practice, the monochromatic (single fixed-frequency) oscillations do not exist. So, analyzing a single color from the RGB color system, we understand that it is a mix of single color frequencies in the very narrow interval of frequencies $Δ ω$ .

We will use the first-order members of Eq. (1) in the calculations of the covariance expressions, because an influence of the higher-order members is negligible. The dispersion $D a_{i}$ and, certainly, covariance $K (a_{i}, a_{i})$ of the monochromatic oscillations of frequency $ω_{i}$ has expression

Eq. (2)

D a_{i} = K (a_{i}, a_{i}) = M {(a_{i} - M a_{i})}^{2} = M (A_{i}^{2} δ^{2} φ_{i} \cos^{2} {\bar{φ}}_{i}) = A_{i}^{2} \cos^{2} {\bar{φ}}_{i} σ^{2} φ_{i},

where

M a_{i} = \sin {\bar{φ}}_{i}

,

M φ_{i} = {\bar{φ}}_{i}

is the symbol of a mean,

σ_{φ_{i}}

is the standard deviation.

The coefficient of correlation of the same frequency is equal to

Eq. (3)

r (a_{i}, a_{i}) = \frac{K (a_{i}, a_{i})}{σ_{a_{i}} \cdot σ_{a_{i}}} = 1 .

The covariance of the two fixed different frequencies looks like

Eq. (4)

K (a_{i}, a_{j}) = M (δ a_{i} \cdot δ a_{j}) = M (A_{i} \cdot A_{j} δ φ_{i} \times δ φ_{j} \cos {\bar{φ}}_{i} \cdot \cos {\bar{φ}}_{j}) = A_{i} A_{j} \cos {\bar{φ}}_{i} \cdot \cos {\bar{φ}}_{j} M (δ φ_{i} \cdot δ φ_{j}) = 0,

where

δ a_{i} = a_{i} - M a_{i}

,

M (δ φ_{i} \cdot δ φ_{j}) = M δ φ_{i} \cdot M δ φ_{j} = 0

, because

M δ φ_{i} = 0

is a mean of the independent random errors.

The coefficient of correlation of the two fixed different-frequency oscillations, applying Eq. (4), is equal to

Eq. (5)

r (a_{i}, a_{j}) = \frac{K (a_{i}, a_{j})}{σ_{a_{i}} \cdot σ_{a_{j}}} = 0 .

Because in nature and in practice the fixed-frequency oscillations do not exist, so for the analysis of the covariance of the light spectrum colors, we will use the signal compositions, applying their interference. The main summing equation of the interference of two frequencies $a_{i j}$ could be written as follows:⁹

Eq. (6)

a_{i j} = a_{i} + a_{j} + 2 \sqrt{a_{i} \cdot a_{j}} γ \cos Δ φ_{i j},

where

γ

is a coefficient of the frequency coherence and

Δ φ_{i j} = φ_{i} - φ_{j}

is the difference of the oscillation phases. In further calculations we will adopt

A_{i} = 1

and

γ = 1

, because the influence of the permanent multipliers do not come into play during determination of the correlation coefficients.

So we have such a linear expression of a sum of two different frequency interferences.

Eq. (7)

a_{i j} = \sin {\bar{φ}}_{i} + \sin {\bar{φ}}_{j} + 2 \sqrt{\sin {\bar{φ}}_{i} \cdot \sin {\bar{φ}}_{j}} \cos Δ {\bar{φ}}_{i j} + \cos {\bar{φ}}_{i} δ φ_{i} + \cos {\bar{φ}}_{j} δ φ_{j} + {(\sin {\bar{φ}}_{i} \cdot \sin {\bar{φ}}_{j})}^{- 1 / 2} \cos Δ \bar{φ} \times (\sin {\bar{φ}}_{j} \cos {\bar{φ}}_{i} δ φ_{i} + \sin {\bar{φ}}_{i} \cos {\bar{φ}}_{j} δ φ_{j}) + 2 \sqrt{\sin {\bar{φ}}_{i} \cdot \sin {\bar{φ}}_{j}} (- \sin Δ {\bar{φ}}_{i j} δ φ_{i} + \sin Δ {\bar{φ}}_{i j} δ φ_{j}),

where

M a_{i j} = \sin {\bar{φ}}_{i} + \sin {\bar{φ}}_{j} + 2 \sqrt{\sin {\bar{φ}}_{i} \cdot \sin {\bar{φ}}_{j}} \cdot γ \cdot \cos Δ {\bar{φ}}_{i j}

.

An expression of the random error of the oscillations summing interference will be

Eq. (8)

δ a_{i j} = \cos {\bar{φ}}_{i} δ φ_{i} + \cos {\bar{φ}}_{j} δ φ_{j} + {(\sin {\bar{φ}}_{i} \cdot \sin {\bar{φ}}_{j})}^{- 1 / 2} \cos Δ {\bar{φ}}_{i j} \times (\sin {\bar{φ}}_{j} \cos φ_{i} δ φ_{i} + \sin {\bar{φ}}_{i} \cos {\bar{φ}}_{j} δ φ_{j}) + 2 \sqrt{\sin {\bar{φ}}_{i} \cdot \sin {\bar{φ}}_{j}} (- \sin Δ {\bar{φ}}_{i j} δ φ_{i} + \sin Δ {\bar{φ}}_{i j} δ φ_{j}) .

We can ignore the last member of Eq. (8), because its influence is not big and its value $2 \sqrt{\sin {\bar{φ}}_{i} \cdot \sin {\bar{φ}}_{j}} (- \sin Δ {\bar{φ}}_{i j} δ φ_{i} + \sin Δ {\bar{φ}}_{i j} δ φ_{j}) \to 0$ , when $δ φ_{i} \approx δ φ_{j}$ .

The covariance expression of the oscillations summing interference could be written as follows:

Eq. (9)

K (a_{i j}, a_{i k}) = M (δ a_{i j} \cdot δ a_{i k}) = \cos^{2} {\bar{φ}}_{i} σ_{φ_{i}}^{2} + {(\sin {\bar{φ}}_{i} \cdot \sin {\bar{φ}}_{j})}^{- 1 / 2} \cos Δ {\bar{φ}}_{i j} \cdot \sin {\bar{φ}}_{j} \cos^{2} {\bar{φ}}_{i} σ_{φ_{i}}^{2} + {(\sin {\bar{φ}}_{i} \cdot \sin {\bar{φ}}_{k})}^{- 1 / 2} \sin {\bar{φ}}_{k} \cos Δ {\bar{φ}}_{i k} \cos^{2} {\bar{φ}}_{i} σ_{φ_{i}}^{2} + \sin^{- 1} {\bar{φ}}_{i} {(\sin {\bar{φ}}_{j} \cdot \sin {\bar{φ}}_{k})}^{- 1 / 2} \cos Δ {\bar{φ}}_{i j} \cos Δ {\bar{φ}}_{i k} \times \sin {\bar{φ}}_{j} \sin {\bar{φ}}_{k} \cos^{2} φ_{i} σ_{φ_{i}}^{2},

where

M δ φ_{i} = 0

,

M (δ φ_{i} \cdot δ φ_{j}) = M (δ φ_{i} \cdot δ φ_{k}) = M (δ φ_{j} \cdot δ φ_{k}) = 0

because of the independent random errors multiplication product.

Further, Eq. (9) has the following expression:

Eq. (10)

K (a_{i j}, a_{i k}) = [1 + {(\sin {\bar{φ}}_{i} \cdot \sin {\bar{φ}}_{j})}^{- 1 / 2} \sin φ_{j} \cos Δ {\bar{φ}}_{i j}] \times [1 + {(\sin {\bar{φ}}_{i} \cdot \sin {\bar{φ}}_{k})}^{- 1 / 2} \sin φ_{k} \cos Δ {\bar{φ}}_{i k}] \cos^{2} {\bar{φ}}_{i} σ_{φ_{i}}^{2} .

By using the similar calculations, we could write the formulae of the dispersion of the oscillations summing interference:

Eq. (11)

D a_{i j} = σ_{a_{i j}}^{2} = {[1 + {(\sin {\bar{φ}}_{i} \cdot \sin {\bar{φ}}_{j})}^{- 1 / 2} \sin {\bar{φ}}_{j} \cos Δ {\bar{φ}}_{i j}]}^{2} \times \cos^{2} {\bar{φ}}_{i} σ_{φ_{i}}^{2} + {[1 + {(\sin {\bar{φ}}_{i} \cdot \sin {\bar{φ}}_{j})}^{- 1 / 2} \sin {\bar{φ}}_{i} \cos Δ {\bar{φ}}_{i j}]}^{2} \times \cos^{2} {\bar{φ}}_{j} σ_{φ_{j}}^{2} .

The formulae of the correlation coefficient of the two oscillations summing interference:

Eq. (12)

r (a_{i j}, a_{i k}) = \frac{K (a_{i j}, a_{i k})}{σ_{a_{i j}} \cdot σ_{a_{i k}}} .

In further calculations we will adopt $σ_{φ_{i}} = σ_{φ_{j}} = σ_{φ_{k}} = σ_{φ}$ .

The accuracy of the phase measurements of the white light colors is equable, i.e., the standard deviations of the oscillation phases of the discrete frequencies are equal, $σ_{φ_{i}} = \dots = σ_{φ_{k}} = σ_{φ}$ . It comes from the assumption that the distributions of the random errors of the spectrum components are asymptotically similar.

Equations (3), (5), and (12) show that correlation between oscillations of different frequencies does not exist, whereas it exists between oscillations of similar frequencies. The correlation between the compound oscillations of mixed frequencies exists in the case when these oscillations have the components of similar frequencies.

3. Model of the Covariance Functions of the Digital Images

To develop the theoretical model we will apply the presumption that the errors of the digital image pixel parameters are random. The random function is constructed by spreading the arrow of the digital image pixels according to columns into one-dimensional space along the same coordinate axis. In each column of the pixel arrays the trend of the corresponding column is eliminated. The parameters are the indexes of the pixel color intensities in the RGB format color spectrum. We will accept the random function constructed in such a way as a stationary function (in wide understanding), i.e., its mean $M [φ (t)] \to const$ and a covariance function depend on the argument difference only, $K_{φ} (τ)$ . The pixel arrows of the two segments of the single digital image or segments of the two images $h_{l} (u)$ ir $h_{j} (u + τ)$ , accepted as the realizations of the random functions, continuous covariance function $K_{h} (τ)$ could be written as follows:¹⁰^,¹¹

Eq. (13)

K_{h} (τ) = \frac{1}{T - τ} \int_{0}^{T - τ} δ h_{l} (u) δ h_{j} (u + τ) d u,

where

δ h_{l} (u)

,

δ h_{j} (u + τ)

is the centered pixel parameters segments,

u

is a parameter of the segment pixel,

T

is the length of the segment in conventional units,

τ = k \cdot Δ

is a changeable quantise interval,

Δ

is a value of the pixel parameter,

k

is the number of pixels in the quantise interval.

The covariance function $K_{h}^{'} (τ)$ based on the measurement results could be estimated according to the following formulae:

Eq. (14)

K_{h}^{'} (τ) = K_{h}^{'} (k) = \frac{1}{n - k} \sum_{i = 1}^{n - k} δ h_{l} (u_{i}) δ h_{j} (u_{i + k}),

where

n

is the total number of discrete intervals.

Equation (14) could be applied in the form of the autocovariance and intercovariance function. In the case of the autocovariance function, the segments $h_{l} (u)$ and $h_{j} (u + τ)$ are parts of the single digital image, and in the case of intercovariance function, these segments are parts of the two different images.

The estimation of the normalized covariance function is

Eq. (15)

R_{h}^{'} (k) = \frac{K_{h}^{'} (k)}{K_{h}^{'} (0)} = \frac{K_{h}^{'} (k)}{σ_{h}^{' 2}},

where

σ_{h}^{'}

is an estimation of the standard deviation of the random function.

To eliminate the trend in the $i$ column of the image pixels array, we can use the formula

Eq. (16)

δ H_{i} = H_{i} - e \cdot {\bar{h}}_{i}^{T} = (δ h_{i 1}, δ h_{i 2}, \dots, δ h_{i m}),

where

δ H_{i}

is an array of reduced pixels of the

i

digital image, in which column trend was eliminated;

H_{i}

is an array of the pixels of the

i

image,

e

is a unit vector, dimension of which is (

n \times 1

);

n

is the number of rows of the

i

array,

{\bar{h}}_{i}

is a vector of the means of the

i

pixels array column,

δ h_{i j}

is a

j

column (vector) of the reduced pixels of the

i

array.

The vector of the means of the columns of the $i$ pixels array could be calculated according to

Eq. (17)

{\bar{h}}_{i}^{T} = \frac{1}{n} e^{T} \cdot H_{i}

or

Eq. (18)

{\bar{h}}_{i} = \frac{1}{n} H_{i}^{T} \cdot e .

The realization of the random function of the $i$ pixels array of the digital image in the form of vectors has an expression

Eq. (19)

δ h_{i} = (\begin{matrix} δ h_{i 1} \\ δ h_{i 2} \\ \dots \\ δ h_{i m} \end{matrix}) = {(δ h_{i 1}^{T} δ h_{i 2}^{T} \dots δ h_{i m}^{T})}^{T} .

An estimation of the autocovariance matrix of the $i$ pixels array of the digital image looks like

Eq. (20)

K^{'} (δ H_{i}) = \frac{1}{n - 1} δ H_{i}^{T} δ H_{i} .

An estimation of the covariance matrix of the two digital images or two pixels arrays of the single digital image could be written as follows:

Eq. (21)

K^{'} (δ H_{i}, δ H_{j}) = \frac{1}{n - 1} δ H_{i}^{T} δ H_{j},

where the dimensions of the arrays

δ H_{i}

,

δ H_{j}

should be equal.

Applying the theory of covariance functions, the influence of the RGB format color spectrum components on the expressions of the covariance functions of the digital images was analyzed. Also, the expressions of the covariance functions of the digital images were estimated using the RGB colors continuous spectrum in the sense of the color tensor.

The estimations of the covariance matrixes $K^{'} (δ H_{i})$ and $K^{'} (δ H_{i}, δ H_{j})$ are reduced to the estimations of the correlation coefficient matrixes $R^{'} (δ H_{i})$ and $R^{'} (δ H_{i}, δ H_{j})$ .¹⁰^,¹¹

Eq. (22)

R^{'} (δ H_{i}) = D_{i}^{- 1 / 2} K^{'} (δ H_{i}) D_{i}^{- 1 / 2},

Eq. (23)

R^{'} (δ H_{i}, δ H_{j}) = D_{i j}^{- 1 / 2} K^{'} (δ H_{i}, δ H_{j}) D_{i j}^{- 1 / 2},

where

D_{i}

,

D_{i j}

are the diagonal matrixes of the main diagonal members of the corresponding covariance matrix estimations

K^{'} (δ H_{i})

and

K^{'} (δ H_{i}, δ H_{j})

.

4. Results of the Experiment and Analysis

The digital images of the human eye [right oculus dexter (OD) and left oculus sinister (OS)] retinas were used in the analysis. The digital images were taken by ordinary photo camera and coded in JPEG format. Example of the OD retina’s digital image is presented in Fig. 1. The calculations were executed by computer procedures written in MATLAB.

Fig. 1

A digital image of an OD retina.

The results of the calculations are presented in Table 1 and in Figs. 2 Fig. 3 Fig. 4 Fig. 5 Fig. 6 Fig. 7 Fig. 8 through 9. In the calculations of the covariance functions, the values of quantise interval were changed from 1 pixel till $n / 2 pixels$ (in our case, $n = 120000$ , average number of pixels in the segment of a digital image). An analysis was done applying all RGB color tensors and its components—red R, green G, and blue B colors.

Table 1

The mean values of the correlation coefficients between the pixels arrays of the digital images.

Combinations of the arrays of the digital images of the right OD and left OS eyes	Mean value of the correlation coefficient
	Red R	Green G	Blue B	RGB tensor
Between OD pixel columns	0.54	$- 0.15$	$- 0.17$	0.98
Between OS pixel columns	0.40	0.13	$- 0.11$	0.98
Between OD and OS pixel arrays	$- 0.42$	$- 0.097$	0.11	0.95

Fig. 2

The changes of the correlation coefficients between the five columns of pixel arrays of the digital image of the OD retina.

Fig. 3

The changes of the correlation coefficients between the five columns of pixel arrays of the digital image of the OS retina.

Fig. 4

The changes of the correlation coefficients between the five columns of pixel arrays of the digital images of the OD and OS retinas.

Fig. 5

A normalized autocovariance function of the digital image of the OD retina.

Fig. 6

A normalized autocovariance function of the digital image of the OS retina.

Fig. 7

A normalized covariance function of the digital images of the OD and OS retinas.

Fig. 8

The changes in percents of the matrix of the covariance coefficients of OD and OS retinas.

Fig. 9

Graphical view of the spatial matrix of the correlation coefficients of the digital images of the OD and OS retinas.

In Table 1 the summarized mean values of the correlation coefficients between the pixels arrays of the digital images of the right OD and left OS eye retinas are shown. The calculations were executed along the single RGB color spectrum components and according to RGB spectrum tensor.

In the calculations of the normalized covariance functions, the results are presented only for red color of RGB spectrum, because the changes of the results for other spectrum colors are negligible.

From the data in Table 1 we can see that the values of the correlation coefficients between OD and OS retina arrays are near zero along all RGB spectrum colors, except red color, where correlation coefficients between OD and OS pixel columns are $r_{OD} = 0.54$ and $r_{OS} = 0.40$ . This means that retinas absorb (accept) all the color oscillations at the same level, with an exception R color. The values of the correlation coefficients along all the RGB color tensors are near one and this means there is a very good relation between both retinas along all the RGB frequencies complex.

The change ( $- 0.3 < r_{2} < 1.0$ ) of the correlation coefficients between the five columns of pixel arrays of the digital image of the OD retina is shown in Fig. 2, the change ( $- 0.4 < r_{2} < 1.0$ ) of the correlation coefficients between the five columns of pixel arrays of the digital image of the OS retina is shown in Fig. 3, and the change ( $- 0.7 < r_{4} < 0.3$ ) of the correlation coefficients between the five columns of pixel arrays of the digital images of the OD and OS retinas is presented in Fig. 4.

A normalized autocovariance function of the digital image of the retina OD is presented in Fig. 5. It describes a change of the correlation due to changes of the quantise interval $k$ between pixels. The values of the correlation coefficient are decreasing from $r = 1.0$ at $k = 0$ to $r < 0$ , 1 when $k = 13000$ . A normalized autocovariance function of the digital image of the retina OS is presented in Fig. 6. The values of the correlation coefficient are changing from $r = 1.0$ at $k = 0$ to $r < 0.1$ when $k = 10000$ . Last results show the decreasing of the correlation coefficient values to $r \to 0$ at a large value of $k = 10000$ .

A normalized covariance function of the digital images of the OD and OS retinas is shown in Fig. 7.

The values of covariance function at $k = 0$ are not big ( $- 0.25 < r < 0.25$ ); however, when quantise interval is increasing, the values of covariance function are slightly increasing too, and when $k = 20,000$ they are changing in the interval $- 0.3 < r < 0.45$ .

A scatter in percentages of the correlation coefficient matrix values of the digital images of the OD and OS retinas is shown in Fig. 8.

The negative correlation occupies about 50% of both array areas. Graphical view of the spatial matrix of the correlation coefficients of the digital images of the OD and OS retinas is presented in Fig. 9.

5. Conclusions

1. The mean correlation between the pixel columns of the arrays of the digital images of the human eye OD and OS retinas along the single RGB component is not big or even negative. However, the values of the correlation coefficients between pixel columns of the OD and OS arrays applying total RGB color tensors are near one. This means that OD and OS retina correlation is very strong.
2. A normalized autocovariance function of the digital image allows us to determine the change of the correlation depending on the pixel quantise interval. The estimations of the normalized covariance functions do not differ much when applying the different RGB spectrum components. The values of a normalized covariance function of the digital images of the OD and OS retinas are decreasing very slowly and approach zero $r \to 0$ at a large value of the quantise interval $k = 10000$ . This shows a very big correlation between the pixels, which are close to each other. The values of a normalized covariance function of the digital images of the OD and OS retinas are not big ( $- 0.2 < r < 0.45$ ) in all the quantise intervals. This shows that correlation between OD and OS retinas is weak.
3. It was detected that negative correlation between the pixels occupies $\sim 50 %$ of both OD and OS retina array areas and that shows the graphical view of the spatial matrix of the correlation coefficients.

References

1.

N. Pattonaet al., “Retinal image analysis: concepts, applications, and potential,” Prog. Retin. Eye Res., 25 (1), 99 –127 (2006). http://dx.doi.org/10.1016/j.preteyeres.2005.07.001 PRTRES 1350-9462 Google Scholar

2.

M. D. AbràmoffM. K. GarvinM. Sonka, “Retinal imaging, and image analysis,” IEEE Rev. Biomed. Eng., 3 169 –208 (2010). http://dx.doi.org/10.1109/RBME.2010.2084567 1937-3333 Google Scholar

3.

N. KardoulasA. C. BirdA. I. Lawan, “Geometric correction of SPOT and Landsat imagery: a comparison of map and GPS derived control points,” Photogramm. Eng. Rem. Sens., 62 (10), 1173 –1177 (1996). PERSDV 0099-1112 Google Scholar

4.

M. EkstromA. McEwen, “Adaptive box filters for removal of random noise from digital images,” Photogramm. Eng. Rem. Sens., 56 (4), 453 –458 (1990). PERSDV 0099-1112 Google Scholar

5.

G. Horgan, “Wavelets for SAR image smoothing,” Photogramm. Eng. Rem. Sens., 64 (12), 1171 –1177 (1998). PERSDV 0099-1112 Google Scholar

6.

B. HuntT. W. RyanF. A. Gifford, “Hough transform extraction of cartographic calibration marks from aerial photography,” Photogramm. Eng. Rem. Sens., 59 (7), 1161 –1167 (1993). PERSDV 0099-1112 Google Scholar

7.

J. P. Antoine, “Wavelet analysis of signals and images. A grand tour,” Revista Ciencias Matematicas (La Habana), 18 113 –143 (2000). Google Scholar

8.

D. E. DutkayP. E. T. Jorgensen, “Wavelets on fractals,” Rev. Mat. Iberoamericana, 22 131 –180 (2006). 0213-2230 Google Scholar

9.

J. Skeivalas, Theory and Practice of GPS Networks, Technika, Vilnius (2008). Google Scholar

10.

J. Skeivalas, “An accuracy of determination of the covariation of random values,” Geodesy and Cartography, 25 (4), 156 –158 (1999). 1392-1541 Google Scholar

11.

J. SkeivalasR. Kizlaitis, “The application of photogrammetric numerical methods to the analysis of magnetic resonance images,” Geodesy and Cartography, 35 (2), 50 –54 (2009). http://dx.doi.org/10.3846/1392-1541.2009.35.50-54 1392-1541 Google Scholar

Biography

Jonas Skeivalas is a author of 3 monographs and more than 150 research papers. His research interests include processing of measurements with respect to tolerances, adjustment of geodetic networks.

Eimuntas Parseliunas is a author of two teaching books and more than 50 scientific papers. He participated in many intern conferences. His research interests include graphs theory in geodesy, adjustment of geodetic networks, and geoinformation systems.

CC BY: © 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.

Citation Download Citation

Jonas Skeivalas and Eimuntas K. Parseliunas "On identification of human eye retinas by the covariance analysis of their digital images," Optical Engineering 52(7), 073106 (15 July 2013). https://doi.org/10.1117/1.OE.52.7.073106

Published: 15 July 2013

Access the abstract

JOURNAL ARTICLE
7 PAGES

DOWNLOAD PAPER SAVE TO MY LIBRARY

GET CITATION

CITATIONS

Cited by 15 scholarly publications.

Explore citations on Lens.org

KEYWORDS

Retina

RGB color model

Eye

Image segmentation

Image analysis

Optical engineering

Wavelets

1.

Introduction

2.

Covariance Model of the Light Frequencies Spectrum

Eq. (1)

Eq. (2)

Eq. (3)

Eq. (4)

Eq. (5)

Eq. (6)

Eq. (7)

Eq. (8)

Eq. (9)

Eq. (10)

Eq. (11)

Eq. (12)

3.

Model of the Covariance Functions of the Digital Images

Eq. (13)

Eq. (14)

Eq. (15)

Eq. (16)

Eq. (17)

Eq. (18)

Eq. (19)

Eq. (20)

Eq. (21)

Eq. (22)

Eq. (23)

4.

Results of the Experiment and Analysis

Fig. 1

Table 1

Fig. 2

Fig. 3

Fig. 4

Fig. 5

Fig. 6

Fig. 7

Fig. 8

Fig. 9

5.

Conclusions

References

Biography

Show All Keywords

Keywords/Phrases

Search In:

Publication Years