UDC 535.42:517.443:51.74
DOI: 10.18698/0236-3933-2022-4-70-79
THE INFLUENCE OF THE APERTURE STOP FRACTAL SHAPE OF AN OPTICAL SYSTEM ON THE ILLUMINANCE DISTRIBUTION
V.I. Zavarzin S.B. Kaledin S.V. Yakubovskiy
zavarzin@bmstu.ru
sbkaledin@mail.ru
yakubovskiystas@mail.ru
Bauman Moscow State Technical University, Moscow, Russian Federation
Abstract
In the article the options for the application of aperture shapes with fractal properties in the design of optical systems are considered. Calculations of mathematical models of point spread functions of a diffraction-limited optical system are performed. The diffraction patterns of the light distribution in these systems are presented, and the point spread functions are considered for various shapes of the aperture stop. Analytical expressions are obtained for the light distribution depending on the pupil shape, which can be used to control the process of image formation. The pupil shape, which has the shape of an equilateral triangle, is chosen as the basic one, and the shape of the pupil as a "Koch snowflake" curve is also considered. Using the Fraunhofer integral, the dependences of the distribution of the spectral density of the complex amplitude on the aperture located on an opaque screen are derived in the Fraunhofer approximation and under the condition of illumination by a plane monochromatic wave. Using the relationship with the complex amplitude, the sought-for intensity distribution in the plane of the diffraction pattern is obtained. Taking into account the simplifications adopted in this article, the solution of the Fraunhofer integral is found, by setting the integration limits, depending on: the selected aperture profile, the coordinate system chosen for it, and the position of nodal points in this system
Keywords
Aperture stop, aperture, point spread function, Fraunhofer diffraction, fractal
Received 14.04.2022 Accepted 27.05.2022 © Author(s), 2022
Introduction. There is a range of practical tasks that require optical systems with the particular intensity distribution of point source image, when it is necessary to provide a "smooth" distribution of energy rather than minimum spot size. When it comes to practice, optical systems with circular apertures are used
in most cases. Nevertheless, there are optical systems with apertures of more complex shapes [1], which require evaluating their impact on image quality.
Formulation of the problem. In this case, a possible option could lie in using an aperture of a particular shape, on account of which an appropriate point spread function (PSF), which characterizes system resolution [2-12], will be formed in an optical system. Particularly, in our opinion, a fractal structure may be a suitable shape, as it is demonstrated in studies [13-17]. The possibility of using suchlike shapes when designing an optical system allows to laying properties that affect system resolution in advance before the synthesis of its components. This way the necessity of complicating the optical system is eradicated. However, the main difficulty of taking the approach lies in getting the analytical dependencies needed for the calculation of light distribution passing through such an aperture. The process of forming a picture could be regulated with their help.
This study aims is to examine the impact of the aperture's shape on the formation of the point spread function of a diffraction-limited optical system.
Applicability condition, assumptions, problem-solving method. The solution to the problem is viewed in the initial approximation. We regard the original optical system as a diffraction-limited one which forms an aberrationless image of a point source in the Fraunhofer diffraction approximation. A flat wavefront from a distant monochromatic source is incident on the aperture stop. The aperture stop plane is considered an opaque screen with a hole, and the plane of image formation is taken as the observation plane or the plane of the diffraction pattern. Spectral density distribution of complex amplitude versus aperture on an opaque screen during Fraunhofer diffraction and illumination by a plane monochromatic wave (Fig. 1), according to [18, 19], is defined by the formula:
U (vx, vy ) ~ {j^U(x, y)exp[-j2n(vxX +vyy )]dxdy, (1)
where vx x / X = ^ / (Xz0), vy = yy / ^ = ^ / (Xz0) are spatial frequencies; U (x, y) is complex wave amplitude at a point P (x, y); X is wavelength; yx, yy are diffraction angles; z0 is distance between planes.
In the calculations, we make the following assumptions: the size of the hole on the screen and the size of the observation area are much smaller than the distance from the screen to the observation plane; the size of the hole on the screen is large compared to the wavelength of the incident light, and the light source is monochromatic and located at such a distance from the screen that the light incident on the screen has an almost flat wavefront and a constant amplitude. Therefore, for a monochromatic plane wave we take U (x, y) = E = const.
Fig. 1. Diffraction of a monochromatic wave in the general case
The advantage of the Fraunhofer diffraction approximation lies in the correlation of spatial frequencies with diffraction angles, which allows one to make the transition from spectral density to the most complex amplitude by the corresponding change of coordinates, rather than calculating the inverse Fourier transform instead. Thus, speaking of the distribution of the spectral density of the complex amplitude from the aperture, we can also judge the distribution of the complex field amplitude in the plane of the diffraction pattern. Similar reasoning is used for the intensity distribution.
Accordingly, the intensity distribution in the diffraction pattern plane is associated with the distribution of the complex amplitude by the known expression [18, 19]:
l(vx, v y ) = U (vx, v y )U* (vx, V y ), (2)
where U (vx, vy ) is complex conjugate spectrum of complex amplitude.
Thus, the solution to the diffraction problem is reduced to finding the integral (1) and substituting the resulting expression in the formula (2). Then, the corresponding coordinate replacement is carried out considering the diffraction approximation, and the intensity distribution function in the image of a point source, i.e., point spread function of a diffraction-limited system, is obtained. The solution of integral (1), taking into account the correctly specified integration limits, is determined by the aperture profile and the reference coordinate system selected for it.
Calculation of Fraunhofer diffraction by triangle aperture. At first, we consider the shape of an equilateral triangle as the basic shape of the aperture. Within the normal incidence of a plane wave on an equilateral triangular hole (Fig. 2), we have the coordinates of points 1 and 2.
The Euler formula has to be used to calculate the integral (1). We choose a counterclockwise direction of travel (positively oriented contour), and then apply Green's theorem [20]. The closed-loop integral is found as the sum of three integrals taken on each side of the triangle. After calculating each of the three amplitudes and their addition, we obtain the intensity distribution according to the formula (2):
Fig. 2. Triangle aperture
2 „2
J(v*, Vy ) = \ -3-
v 7 у 16^2 V2
sin С I — (vx - yßvy )
X Sin
л/3|
^ V y + yßv x ) - sin i — (v y + yßv x)
V
sin С I x + V3Vy j
V y
cos
2
an
(vy + ylbvx )
V
X Sin С
r a%
v T
(vx - y/3vy j I + sinc^an^vx + V3'vy jj - 2sinc(a2^vx)
(3)
In formula (3), we take the natural logarithm of the expression for the so-called "amplitude amplification effect". In addition, if a graph of level lines (top view in the sagittal plane) of the obtained dependence is made, then it is called the diffraction pattern of the light distribution in the image plane.
The graph of the tangential section was built using the expression (3) with the replacement of spatial frequencies by spatial coordinates denoted as x'(^) and y ' respectively, where y' = 0, is called the point spread function (PSF). The same graph characterizes the parameters of the scattering spot. All dependencies are normalized in amplitude.
Calculation of the Fraunhofer diffraction from the fractal aperture — "Koch snowflake". The mechanism of formation of the fractal is shown in Fig. 3 and presented in [21].
We assume that the initial triangle with sides a (see Fig. 2), the distribution of the field from which is already known, is the zero iteration count.
Fig. 3. The "Koch snowflake" formation (4 iteration steps)
У
Therefore, to find the intensity distribution formula in the case of the first iteration, it is necessary to consider the formation of a six-pointed star from the initial triangle (Fig. 4) and find the field distribution from it, that is, the field distribution from three more triangles in addition to the original one.
The distributions of light intensities in the image plane from diffraction by triangles 3-4-5, 6-7-8, and 9-10-11, calculated by the formula (1), are identical to each other.
In the subsequent formation of new triangles (the following iterations), the self-similarity of the field distribution from the fractal can be verified. Therefore, we formulate the final light distribution intensity during diffraction from the hole in the shape of the "Koch snowflake" curve of the n-th iteration, considering the number of newly formed triangles at each step n:
2 л
, V /\ у
\ \ /
\ / \ /
\ / \ /
\ / \ /
7 \ .
Fig. 4. Formation of a six-pointed star from the initial triangle
is (v x, V y ) =
3a2 E 2
16tc2vX
sin c
na
+ sin c
^(vx - V3vy)
- 2 sin c
x sin c
-^a (v x - yßv y )
(vx + V3Vy )
-y- (vx + л/э-Vy ) cos(navx)>+
+ X
n = 1
22(n - 3) a2E2 9n -1 n2v2x
sin c
na
2 • 3
- (v x + V3v y )
+ sin c
na
2 • 3
-(vx - V3vy)
- 2 sin c
na
2 • 3
- (v x + V3v y )
x sin c
Ш
2 • 3n
(vx - V3vy )
cos
Ш
2 • 3n
"V x
(4)
The diffraction pattern is shown in Fig. 5. As you can see, the distribution structure is uniform in all directions. The given distributions in all figures are calculated at a wavelength of X = 700 nm, side a = 2 ^m and amplitude E = 106.
Discussion of the obtained results. A comparison of the point spread functions for each of the cases considered is carried out with a classic round aperture without shielding and a rectangular aperture. We take the radius of the circumference of the round hole equal
to R = a / V3 and assume that both the triangular shape of the pupil and the fractal shape of the "Koch snowflake" curve are inscribed
in the same circle. The intensity distribution from the round pupil in the Fraunhofer diffraction is written as [18]:
s2
Fig. 5. Diffraction pattern of light
distribution from the aperture in the shape of a "Koch snowflake" of the 4th iteration
i(v x, v y ) = | E nR
2 /i(2rcvrR) 2%vrR
(5)
where vr
vx + V2
vy is radial spatial frequency.
Accordingly, the intensity distribution from the rectangular pupil in the Fraunhofer diffraction is
J(v x, v y ) = (Eab )2 sin c2 (aKV x )sin c2 (bnv y ),
(6)
where the ratio of length a to width b in the rectangle is 5 : 4.
The point spread functions of the diffraction-limited system for each shape of the aperture are presented, according to the definition given above, in Fig. 6. As can be seen in Fig. 6, the first PSF minimum for a triangular hole coincides
CIRCLE
RECTANGLE
TRIANGLE
SNOWFLAKE
Fig. 6. Point spread functions for cases when the aperture is round, rectangular, triangular and in the shape of the "Koch snowflake"
with the second PSF minimum for the round one, the scattering spot from the triangular hole is larger, its distribution is more uniform, and it reduces the effect of secondary maxima relative to the PSF for a circular aperture, which can be used in application packages as an image noise filtering function [2, 22, 23].
In the case of the "Koch snowflake", we have an even larger scattering spot than the PSF for the triangle, while the PSF for the snowflake goes around the maxima of the PSF for the circle and is even more uniform than for the triangle. PSF from a rectangular hole has specific large secondary maxima, which are larger than in the rest of the cases considered.
Thus, the shape of the aperture entails a specific change in the scattering spot, which allows one to use PSF to analyze the properties and image quality of the optical system [11, 12, 24, 25]. When giving the aperture shape fractal properties ("fractalization" of the shape) we have a smoother intensity distribution in the image space. Such intensity distribution can be applied, for example, in portrait lenses or laser systems, where high-resolution structures can play
a negative role for the user and, therefore, it may be necessary to "blur" the image (increase the scattering spot). To do so, it is unnecessary to bring in defocusing of the system, change its design and dimensions, or achieve an exponential decrease in the amplitude transmittance function, which also leads to a more uniform intensity distribution and an increase in the scattering spot. It is enough to replace the shape of the aperture stop from round to fractal, which speeds up the system design process to some extent.
Conclusion. A mathematical description of the point spread function formation process in the Fraunhofer diffraction approximation is proposed. A comparative analysis of the intensity distribution in the image of a point for a round, triangular and fractal ("Koch snowflake") apertures is carried out and it is proved that the fractal aperture shape can be used to increase the uniformity of the image illumination in the optical system.
Based on the developed methodology, it is possible to carry out a preliminary analysis of image quality by the described point spread functions even before designing the optical system itself.
REFERENCES
[1] Grebel H., Jggard D.L., Kim Y. Diffraction by fractally serrated apertures. J. Opt. Soc. Am. A, 1991, vol. 8, no. 1, pp. 20-26. DOI: https://doi.org/10.1364/JOSAA.8.000020
[2] Scepuro N.G., Krasin G.K., Kovalev M.S., et al. Determination of the point spread function of a computer-synthesized lens formed by a phase light modulator. Opt. Spec-trosc., 2020, vol. 128, no. 7, pp. 1030-1034.
DOI: https://doi.org/10.1134/S0030400X20070231
[3] Kompanets I.N., Andreev A.L. Microdisplays in spatial light modulators. Quantum Electron, 2017, vol. 47, no. 4, pp. 294-302. DOI: https://doi.org/10.1070/QEL16293
[4] King S.V., Doblas A., Patwary N., et al. Spatial light modulator phase mask implementation of wavefront encoded 3D computational-optical microscopy. Appl. Opt., 2015, vol. 54, no. 29, pp. 8587-8595. DOI: https://doi.org/10.1364/ao.54.008587
[5] Kovalev M.S., Krasin G.K., Odinokov S.B., et al. Measurement of wavefront curvature using computer-generated holograms. Opt. Express., 2019, vol. 27, no. 2, pp. 1563-1568. DOI: https://doi.org/10.1364/oe.27.001563
[6] Zurauskas M., Dobbie I.M., Parton R.M., et al. IsoSense: frequency enhanced sensor-less adaptive optics through structured illumination. Optica, 2019, vol. 6, no. 3, pp. 370379. DOI: https://doi.org/10.1364/optica.6.000370
[7] Ruchka P.A., Verenikina N.M., Gritsenko I.V., et al. Hardware/software support for correlation detection in holographic wavefront sensors. Opt. Spectrosc., 2019, vol. 127, no. 4, pp. 618-624. DOI: https://doi.org/10.1134/S0030400X19100230
[8] Goncharov D.S., Zlokazov E.Yu., Petrova E.K., et al. Features of the implementation of holographic invariant correlation filters based on a phase liquid-crystal space-time light modulator. Bull. Lebedev. Phys. Inst., 2019, vol. 46, no. 4, pp. 126-129.
DOI: https://doi.org/10.3103/S1068335619040055
[9] Fuentes J.L.M., Fernandez E.J., Prieto P.M., et al. Interferometric method for phase calibration in liquid crystal spatial light modulators using a self-generated diffraction-grating. Opt. Express., 2016, vol. 24, no. 13, pp. 14159-14171.
DOI: https://doi.org/10.1364/oe.24.014159
[10] Naprienko S.A., Medvedev P.N., Raevskikh A.N., et al. Diffraction-based investigation methods in analysing plastic strain zone underneath fracture surface. Herald of the Bauman Moscow State Technical University, Series Mechanical Engineering, 2019, no. 4 (127), pp. 97-110 (in Russ.).
DOI: http://dx.doi.org/10.18698/0236-3941-2019-4-97-110
[11] Kulakova N.N., Kaledin S.B., Sazonov V.N. Error analysis of IR lens focal length measured by a goniometric method. Herald of the Bauman Moscow State Technical University, Series Instrument Engineering, 2017, no. 4 (115), pp. 17-26 (in Russ.).
DOI: http://dx.doi.org/10.18698/0236-3933-2017-4-17-26
[12] Zavarzin V.I., Osipovich I.R. Measurement of modulation transfer function for large-sized lenses. Herald of the Bauman Moscow State Technical University, Series Mechanical Engineering, 2003, no. 2 (51), pp. 76-84 (in Russ.).
[13] Allain C., Cloitre M. Optical diffraction on fractals. Phys. Rev. B, 1986, vol. 33, no. 5, pp. 3566-3569. DOI: https://doi.org/10.1103/PhysRevB.33.3566
[14] Hou B., Xu G., Wen W., et al. Diffraction by an optical fractal grating. Appl. Phys. Lett, 2004, vol. 85, no. 35, pp. 6125-6127. DOI: https://doi.org/10.1063/L1840112
[15] Uno K., Uozumi J., Asakura T. Statistical properties of the Fraunhofer diffraction field produced by random fractals. Appl. Opt., 1993, vol. 32, no. 15, pp. 2722-2729. DOI: https://doi.org/10.1364/AO.32.002722
[16] Sakurada Y., Uozumi J., Asakura T. Diffraction fields of fractally bounded apertures. Opt. Rev., 1994, vol. 1, no. 1, pp. 3-7. DOI: https://doi.org/10.1007/s10043-994-0003-9
[17] Chabassier G., Angéli B., Heliodore F., et al. Optical wave diffraction on fractal objects. Pure Appl. Optics, 1992, vol. 1, no. 1, pp. 41-54.
DOI: https://doi.org/10.1088/0963-9659/1/1/005
[18] Born M., Wolf E. Principles of optics. Cambridge, Cambridge University Press, 1999.
[19] Zavarzin V.I., Kalyuzhny A.I. Numerical calculation of Kirchhoff integral in approximation of Fresnel diffraction. Herald of the Bauman Moscow State Technical University, Series Instrument Engineering, 2000, no. 3 (40), pp. 31-42 (in Russ.).
[20] Algazin O.D., Kopaev A.V. The solution of the mixed boundary value problem of Dirichlet — Neumann for the Poisson equation in a multidimensional. Herald of the Bauman Moscow State Technical University, Series Natural Sciences, 2016, no. 3 (66), pp. 42-56 (in Russ.). DOI: http://dx.doi.org/10.18698/1812-3368-2016-3-42-56
[21] Uozumi J., Kumura H., Asakura T. Fraunhofer diffraction by Koch fractals: the dimensionality. J. Mod. Opt., 1990, vol. 38, no. 7, pp. 1335-1347.
DOI: https://doi.org/10.1080/09500349114551501
[22] Dronnikova S.A., Gurov I.P. Image quality enhancement by processing of video frames with different exposure time. Nauchno-tekhnicheskiy vestnik informatsionnykh tekhnologiy, mekhaniki i optiki [Scientific and Technical Journal of Information Technologies, Mechanics and Optics], 2017, vol. 17, no. 13, pp. 424-430 (in Russ.). DOI: https://doi.org/10.17586/2226-1494-2017-17-3-424-430
[23] Vintaev V.N., Zhilenev M.Yu., Matorin S.I., et al. The technique of formation and correction of high-resolution space images. Nauchnyy rezultat. Informatsionnye tekhnologii [Research Result. Information Technologies], 2016, no. 1 (in Russ.).
DOI: https://doi.org/10.18413/2518-1092-2016-1-1-72-80
[24] Richardson W.H. Bayesian-based iterative method of image restoration. J. Opt. Soc. Am., 1972, vol. 62, no. 1, pp. 55-59. DOI: https://doi.org/10.1364/JOSA.62.000055
[25] Levin A. Blind motion deblurring using image statistics. In: Advances in neural information processing systems. Vol. 16. Cambridge, MIT Press, 2006, p. 841.
Zavarzin V.I. — Dr. Sc. (Eng.), Professor, Department of Laser and Optoelectronic Systems, Dean of the Department of Optoelectronic Instrumental Engineering, Bau-man Moscow State Technical University (2-ya Baumanskaya ul. 5, str. 1, Moscow, 105005 Russian Federation).
Kaledin S.B. — Cand. Sc. (Eng.), Assoc. Professor, Department of Laser and Optoelectronic Systems, Bauman Moscow State Technical University (2-ya Baumanskaya ul. 5, str. 1, Moscow, 105005 Russian Federation).
Yakubovskiy S.V. — Post-Graduate Student, Department of Laser and Optoelectronic Systems, Bauman Moscow State Technical University (2-ya Baumanskaya ul. 5, str. 1, Moscow, 105005 Russian Federation).
Please cite this article as:
Zavarzin V.I., Kaledin S.B., Yakubovskiy S.V. The influence of the aperture stop fractal shape of an optical system on the illuminance distribution. Herald of the Bauman Moscow State Technical University, Series Instrument Engineering, 2022, no. 4 (141), pp. 70-79. DOI: https://doi.org/10.18698/0236-3933-2022-4-70-79