## 1.

## Introduction

A number of papers are contributing to retinal imaging and image analysis.^{1}^{,}^{2} 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\mathrm{\Delta}=\mathrm{const}=0$ and their dispersion is $D\mathrm{\Delta}=\mathrm{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^{3}^{,}^{4} or wavelet theory^{5}6.7.^{–}^{8} 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:

## (1)

$$a(t)=A\text{\hspace{0.17em}}\mathrm{sin}(\phi +{\phi}_{0})=A(\mathrm{sin}\text{\hspace{0.17em}}\overline{\phi}+\delta \phi \text{\hspace{0.17em}}\mathrm{cos}\text{\hspace{0.17em}}\overline{\phi}-\frac{1}{2}{\delta}^{2}\phi \text{\hspace{0.17em}}\mathrm{sin}\text{\hspace{0.17em}}\overline{\phi}+\dots ),$$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 $\mathrm{\Delta}\omega $.

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 ${\omega}_{i}$ has expression

## (2)

$$D{a}_{i}=K({a}_{i},{a}_{i})=M{({a}_{i}-M{a}_{i})}^{2}=M({A}_{i}^{2}{\delta}^{2}{\phi}_{i}\text{\hspace{0.17em}}{\mathrm{cos}}^{2}\text{\hspace{0.17em}}{\overline{\phi}}_{i})\phantom{\rule{0ex}{0ex}}={A}_{i}^{2}\text{\hspace{0.17em}}{\mathrm{cos}}^{2}\text{\hspace{0.17em}}{\overline{\phi}}_{i}{\sigma}^{2}{\phi}_{i},$$The coefficient of correlation of the same frequency is equal to

The covariance of the two fixed different frequencies looks like

## (4)

$$K({a}_{i},{a}_{j})=M(\delta {a}_{i}\xb7\delta {a}_{j})=M({A}_{i}\xb7{A}_{j}\delta {\phi}_{i}\phantom{\rule{0ex}{0ex}}\times \delta {\phi}_{j}\text{\hspace{0.17em}}\mathrm{cos}\text{\hspace{0.17em}}{\overline{\phi}}_{i}\xb7\mathrm{cos}\text{\hspace{0.17em}}{\overline{\phi}}_{j})\phantom{\rule{0ex}{0ex}}={A}_{i}{A}_{j}\text{\hspace{0.17em}}\mathrm{cos}\text{\hspace{0.17em}}{\overline{\phi}}_{i}\xb7\mathrm{cos}\text{\hspace{0.17em}}{\overline{\phi}}_{j}M(\delta {\phi}_{i}\xb7\delta {\phi}_{j})=0,$$The coefficient of correlation of the two fixed different-frequency oscillations, applying Eq. (4), is equal to

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}_{ij}$ could be written as follows:^{9}

## (6)

$${a}_{ij}={a}_{i}+{a}_{j}+2\sqrt{{a}_{i}\xb7{a}_{j}}\gamma \text{\hspace{0.17em}}\mathrm{cos}\text{\hspace{0.17em}}\mathrm{\Delta}{\phi}_{ij},$$So we have such a linear expression of a sum of two different frequency interferences.

## (7)

$${a}_{ij}=\mathrm{sin}\text{\hspace{0.17em}}{\overline{\phi}}_{i}+\mathrm{sin}\text{\hspace{0.17em}}{\overline{\phi}}_{j}+2\sqrt{\mathrm{sin}\text{\hspace{0.17em}}{\overline{\phi}}_{i}\xb7\mathrm{sin}\text{\hspace{0.17em}}{\overline{\phi}}_{j}}\text{\hspace{0.17em}}\mathrm{cos}\text{\hspace{0.17em}}\mathrm{\Delta}{\overline{\phi}}_{ij}\phantom{\rule{0ex}{0ex}}+\mathrm{cos}\text{\hspace{0.17em}}{\overline{\phi}}_{i}\delta {\phi}_{i}+\mathrm{cos}\text{\hspace{0.17em}}{\overline{\phi}}_{j}\delta {\phi}_{j}+{(\mathrm{sin}\text{\hspace{0.17em}}{\overline{\phi}}_{i}\xb7\mathrm{sin}\text{\hspace{0.17em}}{\overline{\phi}}_{j})}^{-1/2}\mathrm{cos}\text{\hspace{0.17em}}\mathrm{\Delta}\overline{\phi}\phantom{\rule{0ex}{0ex}}\times (\mathrm{sin}\text{\hspace{0.17em}}{\overline{\phi}}_{j}\text{\hspace{0.17em}}\mathrm{cos}\text{\hspace{0.17em}}{\overline{\phi}}_{i}\delta {\phi}_{i}+\mathrm{sin}\text{\hspace{0.17em}}{\overline{\phi}}_{i}\text{\hspace{0.17em}}\mathrm{cos}\text{\hspace{0.17em}}{\overline{\phi}}_{j}\delta {\phi}_{j})\phantom{\rule{0ex}{0ex}}+2\sqrt{\mathrm{sin}\text{\hspace{0.17em}}{\overline{\phi}}_{i}\xb7\mathrm{sin}\text{\hspace{0.17em}}{\overline{\phi}}_{j}}(-\mathrm{sin}\text{\hspace{0.17em}}\mathrm{\Delta}{\overline{\phi}}_{ij}\delta {\phi}_{i}+\mathrm{sin}\text{\hspace{0.17em}}\mathrm{\Delta}{\overline{\phi}}_{ij}\delta {\phi}_{j}),$$An expression of the random error of the oscillations summing interference will be

## (8)

$$\delta {a}_{ij}=\mathrm{cos}\text{\hspace{0.17em}}{\overline{\phi}}_{i}\delta {\phi}_{i}+\mathrm{cos}\text{\hspace{0.17em}}{\overline{\phi}}_{j}\delta {\phi}_{j}+{(\mathrm{sin}\text{\hspace{0.17em}}{\overline{\phi}}_{i}\xb7\mathrm{sin}\text{\hspace{0.17em}}{\overline{\phi}}_{j})}^{-1/2}\text{\hspace{0.17em}}\mathrm{cos}\text{\hspace{0.17em}}\mathrm{\Delta}{\overline{\phi}}_{ij}\phantom{\rule{0ex}{0ex}}\times (\mathrm{sin}\text{\hspace{0.17em}}{\overline{\phi}}_{j}\text{\hspace{0.17em}}\mathrm{cos}\text{\hspace{0.17em}}{\phi}_{i}\delta {\phi}_{i}+\mathrm{sin}\text{\hspace{0.17em}}{\overline{\phi}}_{i}\text{\hspace{0.17em}}\mathrm{cos}\text{\hspace{0.17em}}{\overline{\phi}}_{j}\delta {\phi}_{j})\phantom{\rule{0ex}{0ex}}+2\sqrt{\mathrm{sin}\text{\hspace{0.17em}}{\overline{\phi}}_{i}\xb7\mathrm{sin}\text{\hspace{0.17em}}{\overline{\phi}}_{j}}(-\mathrm{sin}\text{\hspace{0.17em}}\mathrm{\Delta}{\overline{\phi}}_{ij}\delta {\phi}_{i}+\mathrm{sin}\text{\hspace{0.17em}}\mathrm{\Delta}{\overline{\phi}}_{ij}\delta {\phi}_{j}).$$We can ignore the last member of Eq. (8), because its influence is not big and its value $2\sqrt{\mathrm{sin}\text{\hspace{0.17em}}{\overline{\phi}}_{i}\xb7\mathrm{sin}\text{\hspace{0.17em}}{\overline{\phi}}_{j}}(-\mathrm{sin}\text{\hspace{0.17em}}\mathrm{\Delta}{\overline{\phi}}_{ij}\delta {\phi}_{i}+\mathrm{sin}\text{\hspace{0.17em}}\mathrm{\Delta}{\overline{\phi}}_{ij}\delta {\phi}_{j})\to 0$, when $\delta {\phi}_{i}\approx \delta {\phi}_{j}$.

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

## (9)

$$K({a}_{ij},{a}_{ik})=M(\delta {a}_{ij}\xb7\delta {a}_{ik})={\mathrm{cos}}^{2}\text{\hspace{0.17em}}{\overline{\phi}}_{i}{\sigma}_{{\phi}_{i}}^{2}\phantom{\rule{0ex}{0ex}}+{(\mathrm{sin}\text{\hspace{0.17em}}{\overline{\phi}}_{i}\xb7\mathrm{sin}\text{\hspace{0.17em}}{\overline{\phi}}_{j})}^{-1/2}\text{\hspace{0.17em}}\mathrm{cos}\text{\hspace{0.17em}}\mathrm{\Delta}{\overline{\phi}}_{ij}\xb7\mathrm{sin}\text{\hspace{0.17em}}{\overline{\phi}}_{j}\text{\hspace{0.17em}}{\mathrm{cos}}^{2}\text{\hspace{0.17em}}{\overline{\phi}}_{i}{\sigma}_{{\phi}_{i}}^{2}\phantom{\rule{0ex}{0ex}}+{(\mathrm{sin}\text{\hspace{0.17em}}{\overline{\phi}}_{i}\xb7\mathrm{sin}\text{\hspace{0.17em}}{\overline{\phi}}_{k})}^{-1/2}\text{\hspace{0.17em}}\mathrm{sin}\text{\hspace{0.17em}}{\overline{\phi}}_{k}\text{\hspace{0.17em}}\mathrm{cos}\text{\hspace{0.17em}}\mathrm{\Delta}{\overline{\phi}}_{ik}\text{\hspace{0.17em}}{\mathrm{cos}}^{2}\text{\hspace{0.17em}}{\overline{\phi}}_{i}{\sigma}_{{\phi}_{i}}^{2}\phantom{\rule{0ex}{0ex}}+{\mathrm{sin}}^{-1}\text{\hspace{0.17em}}{\overline{\phi}}_{i}{(\mathrm{sin}\text{\hspace{0.17em}}{\overline{\phi}}_{j}\xb7\mathrm{sin}\text{\hspace{0.17em}}{\overline{\phi}}_{k})}^{-1/2}\text{\hspace{0.17em}}\mathrm{cos}\text{\hspace{0.17em}}\mathrm{\Delta}{\overline{\phi}}_{ij}\text{\hspace{0.17em}}\mathrm{cos}\text{\hspace{0.17em}}\mathrm{\Delta}{\overline{\phi}}_{ik}\phantom{\rule{0ex}{0ex}}\times \mathrm{sin}\text{\hspace{0.17em}}{\overline{\phi}}_{j}\text{\hspace{0.17em}}\mathrm{sin}\text{\hspace{0.17em}}{\overline{\phi}}_{k}\text{\hspace{0.17em}}{\mathrm{cos}}^{2}\text{\hspace{0.17em}}{\phi}_{i}{\sigma}_{{\phi}_{i}}^{2},$$Further, Eq. (9) has the following expression:

## (10)

$$K({a}_{ij},{a}_{ik})=[1+{(\mathrm{sin}\text{\hspace{0.17em}}{\overline{\phi}}_{i}\xb7\mathrm{sin}\text{\hspace{0.17em}}{\overline{\phi}}_{j})}^{-1/2}\mathrm{sin}\text{\hspace{0.17em}}{\phi}_{j}\text{\hspace{0.17em}}\mathrm{cos}\text{\hspace{0.17em}}\mathrm{\Delta}{\overline{\phi}}_{ij}]\phantom{\rule{0ex}{0ex}}\times [1+{(\mathrm{sin}\text{\hspace{0.17em}}{\overline{\phi}}_{i}\xb7\mathrm{sin}\text{\hspace{0.17em}}{\overline{\phi}}_{k})}^{-1/2}\text{\hspace{0.17em}}\mathrm{sin}\text{\hspace{0.17em}}{\phi}_{k}\text{\hspace{0.17em}}\mathrm{cos}\text{\hspace{0.17em}}\mathrm{\Delta}{\overline{\phi}}_{ik}]{\mathrm{cos}}^{2}\text{\hspace{0.17em}}{\overline{\phi}}_{i}{\sigma}_{{\phi}_{i}}^{2}.$$By using the similar calculations, we could write the formulae of the dispersion of the oscillations summing interference:

## (11)

$$D{a}_{ij}={\sigma}_{{a}_{ij}}^{2}={[1+{(\mathrm{sin}\text{\hspace{0.17em}}{\overline{\phi}}_{i}\xb7\mathrm{sin}\text{\hspace{0.17em}}{\overline{\phi}}_{j})}^{-1/2}\mathrm{sin}\text{\hspace{0.17em}}{\overline{\phi}}_{j}\text{\hspace{0.17em}}\mathrm{cos}\text{\hspace{0.17em}}\mathrm{\Delta}{\overline{\phi}}_{ij}]}^{2}\phantom{\rule{0ex}{0ex}}\times {\mathrm{cos}}^{2}\text{\hspace{0.17em}}{\overline{\phi}}_{i}{\sigma}_{{\phi}_{i}}^{2}+{[1+{(\mathrm{sin}\text{\hspace{0.17em}}{\overline{\phi}}_{i}\xb7\mathrm{sin}\text{\hspace{0.17em}}{\overline{\phi}}_{j})}^{-1/2}\mathrm{sin}\text{\hspace{0.17em}}{\overline{\phi}}_{i}\text{\hspace{0.17em}}\mathrm{cos}\text{\hspace{0.17em}}\mathrm{\Delta}{\overline{\phi}}_{ij}]}^{2}\phantom{\rule{0ex}{0ex}}\times {\mathrm{cos}}^{2}\text{\hspace{0.17em}}{\overline{\phi}}_{j}{\sigma}_{{\phi}_{j}}^{2}.$$The formulae of the correlation coefficient of the two oscillations summing interference:

## (12)

$$r({a}_{ij},{a}_{ik})=\frac{K({a}_{ij},{a}_{ik})}{{\sigma}_{{a}_{ij}}\xb7{\sigma}_{{a}_{ik}}}.$$In further calculations we will adopt ${\sigma}_{{\phi}_{i}}={\sigma}_{{\phi}_{j}}\phantom{\rule{0ex}{0ex}}={\sigma}_{{\phi}_{k}}={\sigma}_{\phi}$.

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, ${\sigma}_{{\phi}_{i}}=\dots ={\sigma}_{{\phi}_{k}}={\sigma}_{\phi}$. 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[\phi (t)]\to \mathrm{const}$ and a covariance function depend on the argument difference only, ${K}_{\phi}(\tau )$. 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+\tau )$, accepted as the realizations of the random functions, continuous covariance function ${K}_{h}(\tau )$ could be written as follows:^{10}^{,}^{11}

## (13)

$${K}_{h}(\tau )=\frac{1}{T-\tau}{\int}_{0}^{T-\tau}\delta {h}_{l}(u)\delta {h}_{j}(u+\tau )\mathrm{d}u,$$The covariance function ${K}_{h}^{\prime}(\tau )$ based on the measurement results could be estimated according to the following formulae:

## (14)

$${K}_{h}^{\prime}(\tau )={K}_{h}^{\prime}(k)=\frac{1}{n-k}\sum _{i=1}^{n-k}\delta {h}_{l}({u}_{i})\delta {h}_{j}({u}_{i+k}),$$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+\tau )$ 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

## (15)

$${R}_{h}^{\prime}(k)=\frac{{K}_{h}^{\prime}(k)}{{K}_{h}^{\prime}(0)}=\frac{{K}_{h}^{\prime}(k)}{{\sigma}_{h}^{\prime 2}},$$To eliminate the trend in the $i$ column of the image pixels array, we can use the formula

## (16)

$$\delta {H}_{i}={H}_{i}-e\xb7{\overline{h}}_{i}^{T}=(\delta {h}_{i1},\delta {h}_{i2},\dots ,\delta {h}_{im}),$$The vector of the means of the columns of the $i$ pixels array could be calculated according to

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

## (19)

$$\delta {h}_{i}=\left(\begin{array}{c}\delta {h}_{i1}\\ \delta {h}_{i2}\\ \dots \\ \delta {h}_{im}\end{array}\right)={(\delta {h}_{i1}^{T}\delta {h}_{i2}^{T}\dots \delta {h}_{im}^{T})}^{T}.$$An estimation of the autocovariance matrix of the $i$ pixels array of the digital image looks like

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:

where the dimensions of the arrays $\delta {H}_{i}$, $\delta {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}^{\prime}(\delta {H}_{i})$ and ${K}^{\prime}(\delta {H}_{i},\delta {H}_{j})$ are reduced to the estimations of the correlation coefficient matrixes ${R}^{\prime}(\delta {H}_{i})$ and ${R}^{\prime}(\delta {H}_{i},\delta {H}_{j})$.^{10}^{,}^{11}

## (23)

$${R}^{\prime}(\delta {H}_{i},\delta {H}_{j})={D}_{ij}^{-1/2}{K}^{\prime}(\delta {H}_{i},\delta {H}_{j}){D}_{ij}^{-1/2},$$## 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.

The results of the calculations are presented in Table 1 and in Figs. 2Fig. 3Fig. 4Fig. 5Fig. 6Fig. 7Fig. 8 through 9. In the calculations of the covariance functions, the values of quantise interval were changed from 1 pixel till $n/2\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{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 |

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}_{\mathrm{OD}}=0.54$ and ${r}_{\mathrm{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=\mathrm{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

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

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

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).PERSDV0099-1112Google Scholar

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

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

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).PERSDV0099-1112Google Scholar

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

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

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

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

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