CC BY
20
FOCAL-PLANE FIELD WHEN LIGHTING DOUBLE-RING PHASE ELEMENTS
A.V. Ustinov12
1 Image Processing Systems Institute of RAS - Branch of the FSRC "Crystallography and Photonics " RAS, Samara, Russia,
2 Samara National Research University, Samara, Russia
Abstract
The focal-plane field amplitude is calculated when lighting double-ring phase elements by flat and Gaussian beams. Emerging conditions in the minimum or maximum centers, including flat-top maxima, are given. For the field amplitude, we obtain equations that define the radius of the first zero-intensity ring based on the deduced expressions. The root values are listed for several parameters of optical elements and incident beams due to the lack of analytical solutions. Numerical simulation results are given for flat incident beams; they are fully consistent with the theoretical calculations.
Keywords: phase optical elements, double-ring phase elements, focal spot size.
Citation: Ustinov AV. Focal-plane field when lighting double-ring phase elements. Computer Optics 2017; 41(4): 515-520. DOI: 10.18287/2412-6179-2017-41-4-515-520.
Acknowledgements: This work was financially supported by the Russian Foundation for Basic Research (grant No. 16-07-00825).
Introduction
One of the simplest diffractive optical elements (DOE) is a circular plate with thickness step corresponding to phase difference of 180 degrees. Despite its design simplicity, this element happens to have several practical applications.
In particular, such plates are used to align intensity of focused Gaussian beams [1]. Moreover, we mean not to focus on a ring/square, but somewhat differently - to obtain a domain with dimensions close to diffraction limits with approximately uniform intensity (not necessarily with sharp edges). DOE proposed in paper [1] in the form of a ring or square step is the simplest element among all considered elements used to solve similar problems [2 - 4]. In paper [5], two stepped elements are used for one-dimensional dynamic beam correction.
Another application for phase-step plates is sharp focusing. Phase steps are used in accordance with the first Hermite-Gaussian mode to excite an electric-field longitudinal component at the focus center; this was considered in papers [6 - 8]. It is shown in papers [9, 10] that a focal pattern and, accordingly, contribution of the longitudinal component to the central part of focal domain is influenced by mutual arrangement of polarization axis and a phase step curve for the plate. Besides, in paper [11], the phase-step plate is used to study polarization sensitivity of a near-field microscope.
Radial-step phase plates are also used for sharp focusing and in polarization transformations [12]. In this case, they are capable to perform similar functions like either of multi-ring or multi-level phase plates [13 - 16]. Iterative or other optimization algorithms are often used to calculate them. The advantage of simple elements (binary one- or two-ring) lies in their ability to analytically evaluate their activity [17, 18].
Sharp focusing of radial-polarized beams is considered in paper [17]. Interrelation between width of a diaphragm ring with the introduced radial phase step and focal spot dimensions and intensity values is also studied there. It is shown that due to destructive interference generated by rings with different phases, it is possible to skip scalar limits related to the first zero of the zero-order Bessel function. The
minimum focal spot size (FWHM = 0.331) is achieved when diaphragm ring width is 20 per cent of radius of a full aperture. In this case, intensity of side lobes does not exceed 30 per cent of a central peak. It is also shown that by introducing phase steps and simultaneously widening a ring aperture, it is possible to generate the focal spot not exceeding the size limit corresponding to a narrow ring aperture while intensity increases almost 6 times. Furthermore, side lobes embrace 35 per cent of the central peak.
Paper [18] cites evidence that two annular zones are sufficient for many problems. Analytical estimates are obtained for parameters of a two-zone element in off-axis scalar approximation; the element provides maximum values of interference maxima on optical axis. Zone radii coincided with zone-plate radii. However, a zone plate with free-space limiting numerical aperture has its central zone r: = 1.121 in radius, and remaining zones are represented by rings less than 1/2 wide. In this case, decaying waves get mostly through a peripheral part of the optical element; they slightly affect distribution near optical axis. Therefore, it is possible to discuss influence of only one or two central zones. It is shown numerically and analytically that an optical microelement, consisting of only two axially aligned annular zones, can be used for sharp focusing of laser radiation. Moreover, the greatest focusing degree is achieved when the central zone is 1/2 in radius.
Here, as in paper [18], influence of the total element (instead of the ring element) is regarded in far-field diffraction. Though a more accurate off-axis model is used in paper [18], analytical calculation was carried out only for the optical axis field. The transverse in-plane field corresponding to the lens focal plane is analyzed in this paper.
1. Theoretical study of double-ring DOEs
Let's consider a phase DOE with radius r2; its structure is shown in Fig. 1. The phase is equal to zero in the internal circle of radius r1 and to 180 degrees in a ring between the radii r1 and r2. Let's consider the field being generated in the lens focal plane when lighting this element.
In this paper, we examine two types of incident beams: flat-top and Gaussian beams.
Fig. 1. Drawing of a double-ring phase DOE
1.1. Focal_ field when lighting DOEs by _ flat-top beams
The field amplitude in the lens focus plane is determined by the radially symmetrical Fourier transform as follows:
A(p) = j A0(r ) J0r dr.
(1)
where Ao(r) is the field amplitude in the input plane, Jo(x) is the Bessel function of the zero order, k=2p/1 is the wavenumber, l is the wavelength of an illuminating beam, f is the focal length of a lens, which is combined with the a double-ring phase DOE. In our case, we have the following:
Ao(r ) =
1,0 < r < r ,1 -1, r1 < r < r2.
(2)
Substituting Eq. (2) in Eq. (1), we obtain the following expression for the amplitude:
A(p) = f kp
2r1J1
kPr1
- r2 J1
kPr2
f J I f
The zero-point amplitude is equal to the following: r2
A(p = 0) = r> -
(3)
(4)
A(0) = 0, when r? = r22/2; provided r? > r22/2 A(0) > 0,
otherwise A(0) < 0. Since
A'(p) = - j j Ao(r) Ji ^ f j r 2dr , A'(p = 0) = 0 ,
we see that a critical point will be always in zero. The second derivative is equal to the following:
kpr
k ^
A(P) = -I j J
r dr
A"(p = 0) = -|f
2 '
4
r4 ^ 8
wherefrom we obtain A "(0) = 0, when r/ = r24 / 2 ; provided rl4 > r24/2 A"(0) < 0 (maximum), otherwise A"(0) > 0 (minimum). For r/ = r24/2 an extremum shall be flat.
Thus, using the Taylor series expansion, we can prove that this will be the maximum.
Let us introduce the following notation:
m = rj r2(< 1). (5)
Thus, the following three configurations are possible (Fig. 2): the negative minimum (provided 0 < | < 1/\/2 ); the positive minimum (provided 1/ V2 <m< 1/^2 ), and
the positive maximum (provided 1/^2 <|m< 1). The fourth configuration, i.e. the negative maximum, is impossible.
a) b) c)
Fig. 2. Layout view of the amplitude near optical axis: configurations I (a), II (b), III (c)
The zero amplitude position is obtained by setting the right-hand side of Eq. (3) equal to zero; when x = kpr2 / f this gives the following equation:
2|m/1(|jx) = J( x). (6)
The first nonzero root of this equation is considered. The equation is solved numerically by using Bessel function tables, and then, if necessary, the root curve x0 (|) can be approached to the analytic curve. If we know x0, we can determine the radius of the dark ring as follows:
x0 •f
Po =-
kr
(7)
It is directly shown in Eq. (6) that at m = 0 and m = 1, xo =ji,i is the first root of the function Ji(x). Thus solving Eq. (6), we obtain the following values.
Table 1. Zero amplitude position versus ring radius ratio
m 0 0.1 0.2 0.3 0.4 0.5 0.6
X0 3.832 3.742 3.513 3.212 2.866 2.462 1.912
m 0.7 -T - 0 ' + 0 V2 0.8 1 w 0.9 1
X0 0.553 0 6.116 5.539 5.241 4.734 3.832
1.2. Focal_ field when lighting DOEs by Gaussian beams
In this case, the field amplitude in the input plane is as follows:
Ao(r ) = e
= „-(,2/202) I1,0 < r < r1:
-1, r1 < r < r2.
(8)
where o is a beam-width determining parameter. It is impossible to make exact calculations here; however, we can obtain rather close approximation. Breaking down the element into two rings, we get the following:
A (p) =
= j e* Jo | f | r dr - j e202 Jo ^f ] r dr = = B+-B- = 2B+- (B+ + B_).
krp
(9)
The sum in parentheses can be determined with a reasonable degree of accuracy using the tables given below,
r2 r
if the pupil's radius meets the condition r2 >~ (2^2.5)o when an upper limit could be replaced by an infinite limit. The following formula is given in tables presented in paper [19] (p. 186).
J e p"z Jv ( cx ) x dx =
c4P
c 8p
8 PJ
(n-1)/2
d I - ,
8 P
f c2 ^
(v+1)/2
8 IP
(10)
When v = 0, the formula is simply as follows:
r - 2 1
e px J ( cx ) x dx = — e
0 J 2p
c
"4P
Using Eq. (10a), we obtain the following:
B++ B_» J e"2°2 J01 krp | r dr = o
f
_ (*pp)2 ,2e f
(10a)
(11)
We can calculate the value of B+ with the specified degree of accuracy replacing a Gaussian function bell with a parabolic cone, when ri < o (a knee of the Gaussian curve). The case when this condition isn't well fulfilled is uninteresting: an exterior ring has no meaningful effect. So, the result is limited to one circular hole to be lighted.
Suppose the condition ri < o has been satisfied. Then the following approximation should be made with a small
accuracy over the range of r < n:
_
e 2o2 » a1 _a2r2. (12)
The values of a1 and a2 depend on used approximation and they are given below. Using Eq. (12), we receive the following:
2B+ » 2j(a1 -a2r2) J0
krp f
r dr =
= 2 kP J1 ( f J <a'_a2 '' )+ ( ê J2 ^ 2 ( f ) a-
The overall amplitude is equal to the following: , f1 r( k9r1
(13)
A = 2 kp J1 [ ^M (a1 -a2 r12 )-
- (*pp)2
+4\ ^ I JA ^ Ia-o2e 2f2
Ä kp
kPr1 f
(14)
This shows that for large values of p decreasing is determined with the first summand and the rate of decay is identical to that one valid for lighting with flat-top beams. When r1 = 0 (no phase step), we get a trivial response:
_ (kpp)2 A = _o2e 2f2 , which is the Gaussian beam.
Before proceeding any further, it's worthwhile to get at least one technique for finding the values a1 and a2. The simplest method, though probably not the best one in terms of its accuracy, is a collocation method. We equal zero and edge values as follows:
1 = a1;
e-(r.2/2°2) = ,
• ^
a1 = 1,
a2 = (1 - e-»-/2s2)) • l.(15)
Substituting Eq. (15) in Eq. (14), we obtain the following:
A = 2fr1 J1 kp
+ 4 \f kp
J
kPr1 f
kPr1 f
1 - e
-O2e
(kpo)2 2 f2
(16)
The value of A (p = 0) can be calculated exactly without any approximations made. It equals to the following:
A(p = 0) = o2 - 2o2e 2°2 +G2e
(17)
Providing however that r2 is larger (see the line before Eq. (10)), the latter summand shall be neglected:
A(p = 0) = s2 (1 - 2e-(r'/2°2) ).
(17a)
Configurations will be determined according to Eq. (17a) being similar for derivatives. We write out accurate values, just as in Eq. (17), but they are less convenient for analysis due to their complexity.
A(0) = 0, when r1 =o>/2ln2 » 1.1774o. (18)
At less values of n, we have A (0) < 0, otherwise we have
A (0) > 0.
The first derivative A' (p = 0) equals to zero and the second one is as follows:
A'(p = 0) = -\f
+0,5o2 (r22 + 2o2)e
O4 -G2(r12 + 2G2)e
r2 A
(19)
With disregard for the latter summand, we get the following:
( r2 ^
A'(p = 0) = -\f
O
O2 - (r12 + 2o2)e
(19a)
A'(0) = 0, when r1 » 1.832o , (20)
where 1.832 is the root of the equation 1 = (x2+2)e-x2/2. This corresponds to the flat maximum at a smaller vale of A" (0) > 0 (minimum), otherwise we have A" (0) < 0 (maximum). Configuration change takes place similarly as in case with flat-top beams (Fig. 2)
Configuration I: n < 1.1774o;
Configuration II: 1.1774o < n < 1.832o;
Configuration III: n > 1.832o.
0
2
e
202 +
e
2
0
r
0
2
+
It is difficult to compare edges obtained from fully exact expressions because of the fact that there are two radii, not one, in exact expression. However, we can compare them with those ones that will be used in Eq. (16). Using two summands of the Taylor series (Ji(x) ~ (x/2) - (x3/16), J2(x) ~ (x2/8) - (x4/96)), the near-axis amplitude can be conceived of as follows:
A
r2 A
1 + e
-o
kp
f
'1
24
1 + 2e 20
-
2
(21)
Though Eq. (21) is approximate, the zero value is exact; it equals to the following:
2 ( r,2 ^ r —1—
A(0) = 1 + e 2o2 -o2, 2
V J
wherefrom
A(0) = 0 at r: » 1.15o, (22)
where 1.15 is the root of the equation x2(1+e~x2/2) = 2; it is somewhat less than the factor given in Eq. (18). The second derivative is equal to the following:
A V f *
wherefrom
- ZL (1 + 2e-(r-/2°2) ) + o4 12 ^ >
A "(0) = 0 at n » 1.687o, (23)
where 1.687 is the root of the equation x4(1+2e-x2/2) = 12; it is less than the factor given in Eq. (20). Overall, the area of configuration II is narrower than that one obtained by exact calculation. It should be pointed out that the factor given in Eq. (22) and especially in Eq. (23) is larger than permitted; this is necessary to make the approximation mentioned in Eq. (12), on the basis of which Eq. (16) has been obtained, be good.
Zero values of Eq. (16) shall be determined using the below tables. With the use of
x = (kprjf); r =m-o we have A (p) = 0 be written as follows:
J1( x)
2m2 +4m
• e
x
2 J2(x)
x2 '
-m2/2
(1 - e-m2/2 ) = e-x2/(V)
(24)
The equation is solved for x at given m. If we know the root of x0, the radius of the dark ring is determined as follows:
x0 f p0 = — •Tm k-
(25)
In contrast to Eq. (7), a real factor is not x0 itself, but x0/That's why we added a proper line in the table of roots of equation (25) (implication of one more line will be clarified below).
Table 2. Zero amplitude versus circle-and-beam radius ratio
m 0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 1.1 1.15
x0 0 0.3043 0.51 0.669 0.79 0.879 0.936 0.961 0.948 0.888 0.755 0.464 0
x0/m ¥ 3.043 2.55 2.23 1.975 1.757 1.56 1.372 1.185 0.986 0.755 0.422 0
2mJ- in m 0 0.3035 0.5075 0.658 0.766 0.833 0.858 0.836 0.756 0.584 0 - -
In contrast to the flat case, when m was within the range from 0 to 1 (though it had another definition there), here m can theoretically possess any positive value, however really meaningful values are limited: when m > 2 ^ 2.5, this is equivalent to near-zero values with the amplitude sign reversal. Besides, Eq. (16) and therefore Eq. (25) are obtained on the assumption that m <~ 1. Anyway, the values larger than 1.15 (see. Eq. (22)) are not particularly of any interest.
It is apparent from the above table that with reducing m the radius p0 goes up with no limit as it must be. The amplitude constraint tends to the Gaussian curve that possesses no zeros.
With small m we can also determine the root value without the benefit of the tables. The right-hand side of the equation (24) at m®0 and with a fixed value of x tends to zero much faster than the left-hand side; therefore x should be small, too. In this case, we can make the following approximations on the left-hand side of the equation:
-m2/2
a- 1 -e-m2/2 »ml- JM » 1
2 x 2
J 2( x)
Then the left-hand side of the equation will be aver-agely equal to the following:
.2
m
2m2 •11+4m2 •1 = m2+1
2 8 2 4
so, we get the following equation:
m2 = e-*2/(2m2) ® x0 = 2m V- ln m. (26)
The ratio x0 / m = 2^/- ln m goes up with no limit when m®0.
2. Numerical simulation results
Numerical simulation was implemented under the flat incident beam with the following parameter: the outer radius is r2 = 1001 (l is the wavelength of incident radiation equaled to 1 ^m); the lens focal length is f = 4001. Let S denote a ratio between the side lobe intensity and the centered intensity. The results are given in Fig. 3-6.
As we can see, when the size of the central light spot decreases, the energy goes out of it into a peripheral ring.
2
+
2
4
+
e
2
x
40 20 0 -20
A ___
r^- y
a) -4 -2 0 2 4 x/k b) -4 -2 0 2 4 xA
Fig. 3. Focusing results for the element with phase step radius ri = 0.6r2: the element's interior-part field (dashed line, FWHM= 3.431, S = 0.015), the element's front-end field (solid line, FWHM= 1.741, S = 0.123) and their superposition (dotted line,
FWHM = 1.21, S = 1.6). Amplitudes (a) and intensities (b) are shown
Fig. 5 shows that it is also possible to obtain very small sizes for the central spot through its very low
6000
4000
2000
0
/ \
/ / ! / \\
w \
/ t > \
// \\ v\
// S
-4
0
4 x/k
Fig. 4. Focusing results for: ri = 0 (FWHM = 2.051, dotted line), ri = 0.25r2 (FWHM = 1.911 dashed line), and ri = 0.5r2 (FWHM = 1.511, solid line)
brightness and growing side lobes.
If we compare numerical simulation results with theoretical calculations, they should nearly coincide. Sources of discrepancy are as follows: computational errors in solving the Eq. (6); computation errors of numerical simulation; visual accuracy of determining the zero position in graphs obtained by simulation.
700 500 300 100
a)
+
-V
ft
0 2 4 xA b) -0.4 -0.2 0 0.2 0.4 x/k
Fig. 5. Focusing results for: n = 0.6r2 (FWHM = 1.221, dashed line), r1 = 0.65r2 (FWHM = 0.961, solid line), and n = 0.7r2
(FWHM = 0.351, dotted line) (a); central part for n = 0.7r2 (b)
By parameter values from Eq. (7) used in simulation,
1200
600
400
0
.V
/ \
~—^
N \ / f V\\
/// V \ i t
y v.. ___'
-4 -2 0 2 4 x/k
Fig. 6. Focusing results for: r1 = 0.75r2 (dashed line), r1 = 0.8r2 (solid line), ri = 0.85r2 (FWHM = 3.91, dotted line)
we get the following: p0 = (0,637x0) -1.
(27)
Values of the radius of the first dark ring calculated by Eq. (27) and determined according to certain intensity graphs obtained from simulation are shown in Table 3. Values of xo are taken from Table 1 or obtained by interpolation.
Table 3. Comparison of zero points discovered theoretically and by simulation
m 0 0.25 0.5 0.6 0.65 0.7 0.75 0.8 0.85
po as in (27) 2.44 »2.148 1.568 1.218 »0.785 0.352 »3.726 3.528 »3.289
po as in graph =2.5 =2.2 =1.5 =1.2 =0.85 =0.32 =3.75 =3.5 =3.3
Values of the radius are given in units of wavelength. Designated values mean the use of linear interpolation. Note that for m = 0.25 and 0.85 interpolation produces adequate accuracy; for m = 0.75 - mean accuracy, and for m = 0.65 - low accuracy. This is due to the fact that interpolation accuracy increases when a derivative of the interpolated function approximates a constant. It may be concluded that compatibility between theoretical calculations and numerical simulation results is very high.
Conclusion
Analysis of operation of double-ring optical elements allowed us to come to the following conclusions: 1. Three significant configurations for amplitude (two -for intensity) are possible near optical axis. Edge configurations are supposed to be the centered zero value
and the largest flat-topped value similar to the configuration mentioned in paper [1]. A type of configuration is determined by the ratio between inner and outer radii at flat beam incidence or by the ratio between inner radius and width of the Gaussian beam.
2. Equations are derived for the first zero intensity position. Tables of their roots are given.
3. Numerical simulation implemented for the case of incident beams has fully confirmed theoretical predictions.
References
[1] Khonina SN, Kotlyar VV, Skidanov RV, Soifer VA. Leveling the focal spot intensity of the focused Gaussian beam. J Mod Opt 2000; 47(5): 883-904. DOI: 10.1080/09500340008235098.
[2] Veldkamp WB. Laser beam profile shaping with interlaced binary diffraction gratings. Appl Opt 1982; 21(17): 32093212. DOI: 10.1364/A0.21.003209.
[3] Han C-Y, Ishii Y, Murata K. Reshaping collimated laser beams with Gaussian profile to uniform profiles. Appl Opt 1983; 22(22): 3644-3647. DOI: I0.1364/A0.22.003644.
[4] Cordingley J. Application of a binary diffractive optics for beam shaping a semiconductor processing by lasers. Appl Opt 1993; 32(14): 2538-2542. DOI: 10.1364/A0.32.002538.
[5] Passilly N, Fromager M, Mechin L, Gunther C, Eimer S, Mohammed-Brahim T, Ait-Ameur K. 1-D laser beam shaping using an adjustable binary diffractive optical element. Opt Commun 2004; 241(4-6): 465-473. DOI: 10.1016/j.optcom.2004.07.036.
[6] Novotny L, Sánchez EJ, Xie XS. Near-field imaging using metal tips illuminated by higher-order Hermite-Gaussian beams. Ultramicroscopy 1998; 71(1-4): 21-29. DOI: 10.1016/S0304-3991(97)00077-6.
[7] Khonina SN, Golub I. Optimization of focusing of linearly polarized light. Optics Letters 2011; 36(3): 352-354. DOI: 10.1364/OL.36.000352.
[8] Khonina SN, Nesterenko DV, Morozov AA, Skidanov RV, Soifer VA. Narrowing of a light spot at diffraction of linearly-polarized beam on binary asymmetric axicons. Optical Memory and Neural Networks (Information Optics), Allerton Press 2012; 21(1): 17-26. DOI: 10.3103/S1060992X12010043
[9] Khonina SN. Simple phase optical elements for narrowing of a focal spot in high-numerical-aperture conditions. Optical Engineering 2013; 52(9): 091711. DOI: 10.1117/1 .OE.52.9.Ü91711.
[10] Khonina SN, Alferov SV, Karpeev SV. Strengthening the longitudinal component of the sharply focused electric field by means of higher-order laser beams. Optics Letters 2013; 38(17): 3223-3226. DOI: 10.1364/OL.38.003223.
[11] Alferov SV, Khonina SN, Karpeev SV. Study of polarization properties of fiber-optics probes with use of a binary phase plate. J Opt Soc Am A 2014; 31(4): 802-807. DOI: 10.1364/J0SAA.31.000802.
[12] Bokor N, Davidson N. Tight parabolic dark spot with high numerical aperture focusing with a circular p phase plate. Opt Commun 2007; 270(2): 145-150. DOI: 10.1016/j.opt-com.2006.09.022.
[13] Helseth LE. Mesoscopic orbitals in strongly focused light. Opt Commun 2003; 224(4-6): 255-261. DOI: 10.1016/j.opt-com.2003.07.017.
[14] Jabbour TG, Kuebler SM. Vector diffraction analysis of high numerical aperture focused beams modified by two-and three-zone annular multi-phase plates. Opt Express 2006; 14(3): 1033-1043. DOI: 10.1364/OE.14.001033.
[15] Gao X, Wang J, Gu H, Xu W. Focusing properties of concentric piecewise cylindrical vector beam. Optik 2007; 118(6): 257-265. DOI: 10.1016/j.ijleo.2006.10.006.
[16] Xu Q, Chen J. The creation of double tight focus by a concentric multi-belt pure phase filter. Opt Commun 2012; 285(7): 1642-1645. DOI: 10.1016/j.optcom.2011.11.116.
[17] Khonina SN, Ustinov AV. Sharper focal spot for a radially polarized beam using ring aperture with phase jump. Journal of Engineering 2013; 2013: 512971. DOI: 10.1155/2013/512971.
[18] Khonina SN, Savelyev DA, Ustinov AV. Diffraction of laser beam on a two-zone cylindrical microelement. Computer Optics 2013; 37(2): 160-169.
[19] Prudnikov AP, Brychkov YuA, Marichev OI. Integrals and Series. Volume 2: Special Functions. Amsterdam: Gordon and Breach Science Publishers; 1983. ISBN: 2-88124-097-6.
Author's information
Andrey Vladimirovich Ustinov, (b. 1968) graduated from Kuibyshev Aviation Institute named after academician S.P. Korolyov (KuAI) on a specialty "Applied Mathematics" in 1991. Candidate of Physical and Mathematical Sciences (2016), works as the researcher in the IPSI RAS - Branch of the FSRC "Crystallography and Photonics" RAS. Research interests: diffractive optics; software design for modeling of optical elements operating; images processing, particularly images of hydrodynamic processes and biomedical images. E-mail: andr@smr.ru .
Code of State Categories Scientific and Technical Information (in Russian - GRNTI): 29.31.15. Received July 18, 2017. The final version - August 22, 2017.
