
1.IntroductionThe moiré effect is a physical phenomenon of linear optics. The moiré patterns appear as a result of an interaction between transparent layers of a repeated structure^{1} when superposed layers are viewed through. Several examples of the moiré effect are shown in Fig. 1 (a bridge, a facade of a building, an air conditioner grid, a textile curtain, etc.). The visual appearance of the moiré patterns depends on characteristics of gratings and on the location of the observer.^{2} In the case of two gratings, the patterns look like an added series of repeated stripes, which could look bright, vivid, and sometimes even unpredictable. The moiré effect is not unknown in the literature. A general description of the moiré effect can be found in several books (Refs. 1–3), which also include many useful examples. A perfectly illustrated book, Ref. 4 is full of excellent moiré images. There are also research papers, specifically Refs. 56.7.8.9.–10, which describe various aspects of the moiré effect. Examples of the moiré effect in digital devices are shown in Fig. 2. In gratings of similar layouts, the moiré patterns reproduce the structure of the gratings as in the moiré magnifier,^{11} see Fig. 3. This figure also shows that with all other conditions equal, the period of the moiré patterns is the same in gratings of different geometric layouts. Generally speaking, visual displays can be either dynamic or static. The dynamic displays are CRT TV, LCD, and OLED. The examples of static displays are a printed picture or photograph, poster, and postcard. In visual displays, the moiré effect may create an unwanted and sometimes unexpected image of bands (as in Fig. 1) or patterns (as in Fig. 3). Such additional visual background decreases the contrast on the display screen and consequently reduces the image quality; therefore in imaging and displays, the moiré effect is an undesirable adverse visual effect.^{12}^{,}^{13} The elimination or the reduction of the moiré patterns in displays is an important issue of improvement in the visual quality. The necessity of the minimization of the moiré patterns in threedimensional (3D) displays was first stated^{14} in the early 2000s. This is especially important for autostereoscopic 3D displays,^{15}^{–}^{23} where the moiré effect may often occur due to their typical design. The moiré effect was observed not only under visible light but also under rays and beams of different nature, for instance, in electron beams,^{23} infrared light,^{24} and xrays,^{25} as well as at the nanoscale in graphene layers.^{26} Finding characteristics of patterns based on the layouts of the layers is a direct problem.^{27}^{–}^{31} An inverse problem is to find the position (shape) of another grating based on the observed moiré patterns when positions of an observer and of one grating are known together with the parameters of the grating; this is a topic of the 3D shape measurement.^{2}^{,}^{32}^{–}^{37} In contrast, it appears to be possible to arrange the moiré patterns so the resulting pattern would carry useful informational content (including 3D). The first solution of this previously unknown moiré problem was proposed in Ref. 38 and implemented in Ref. 39. Also, the moiré effect can be used for other purposes, such as the moiré interferometry,^{40}^{–}^{42} the moiré deflectometry,^{43} the optical alignment,^{44}^{,}^{45} the visual security (cryptography),^{46}^{,}^{47} and many other applications. That is to say, this is a twofold moiré problem that on the one hand, deals with minimization of the moiré effect to improve the image quality;^{48} and on the other hand, with its maximization to measure distances,^{49}^{,}^{50} as well as to display meaningful images^{38} in 3D displays based entirely on the moiré effect as the main physical principle, or to improve security.^{51}^{–}^{53} An image in a 3D moiré display is shown in Fig. 4. In each application example mentioned in this paragraph, the understanding of the behavior of the moiré patterns is needed. There are many approaches to investigate the moiré effect. In some cases, solutions can be found from direct analytical considerations. The indicial method can be used to find locations of the characteristic points (minima or maxima) of the patterns. The indicial and direct approaches require simplification, such as a simplified structure, small angles, smooth functions (nearly sinusoidal), etc. Dealing with spectra represents a generalized approach.^{54}^{–}^{56} Using spectral trajectories,^{57} the behavior of the patterns can be estimated geometrically. As the equations of the trajectories are derived from the geometric characteristics of the gratings, the estimation can be made without calculations of spectra. This makes the spectral trajectories suitable for an interactive computer simulation. The tutorial is based on the authors’ knowledge and comprehension of the moiré phenomenon and experience in investigating it for a decade. It includes several methods to find characteristics of the moiré patterns in various situations that are presented. After sections, we provide exercises. We hope that the exercises can make readers to think deeply, and this way to improve their understanding of the considered topics. Implied is a selfcheck as, for example, the same result obtained by an independent method. This is probably one of the best methods of verification. Nevertheless, whenever readers would send their solutions to the authors, the readers’ notes and comments will be considered and replied to with the greatest pleasure. The tutorial is arranged as follows: Secs. 2–5 deal with the moiré effect directly, in the spatial domain, whereas Secs. 6–8 deal with spectra (the spectral domain). Namely, in Sec. 2, we explain the indicial equation method. In Secs. 3 and 4, we describe the plain coplanar and noncoplanar sinusoidal gratings. Section 5 gives an example of the moiré effect in a 3D object, a cylinder. Then, Secs. 6 and 7 present such fundamental issues as the wave vector of the moiré patterns in terms of the vector sum and the basics of twodimensional (2D) Fourier transform. Section 8 describes the spectral trajectories. After discussion about the visual effects in displaced or rotated plane gratings in Sec. 9, a conclusion finalizes the tutorial. 2.Indicial EquationThe indicial equation is an analytical method to calculate the characteristic locations of fringes, implying that the locations of the minima/maxima of gratings are known. In this approach, a line grating is modeled by a series of thin lines (a sketch of a family of lines), i.e., some kind of a wireframe of the maxima (or minima) only, with the intensity profile ignored. The lines of the moiré bands connect the intersections of these families of lines of gratings. As soon as the equations of these intersections can be calculated based on the given equations of the families, technically the equations of the moiré bands can be found as analytical expressions. As an easy example, consider two families of the equidistant parallel lines in the $xy$plane: the horizontal lines with period ${T}_{1}$ and the slanted lines with period ${T}_{2}$ (rotated at the angle $\alpha $). The sketch lines of the gratings are shown by solid lines in Fig. 5; the moiré patterns comprise the third family of lines with the period ${T}_{m}$ at the angle $\theta $ (shown by dotted lines connecting the intersections in a shortest way). These spots or relatively low visual densities are visually connected as brighter spaces between darker areas, and the moiré patterns appear. The periods of bright and dark lines are identical. Without loss of generality, consider the maxima. The equations of two families are Eq. (2)$$x\text{\hspace{0.17em}}\mathrm{sin}\text{\hspace{0.17em}}\alpha +y\text{\hspace{0.17em}}\mathrm{cos}\text{\hspace{0.17em}}\alpha =n{T}_{2},$$Excluding $m$ and $n$ from Eqs. (1)–(3), we obtain the equation of the $q$’th intersection, i.e., the $q$’th moiré line Eq. (4)$$x{T}_{1}\text{\hspace{0.17em}}\mathrm{sin}\text{\hspace{0.17em}}\alpha +y({T}_{2}{T}_{1}\text{\hspace{0.17em}}\mathrm{cos}\text{\hspace{0.17em}}\alpha )={T}_{1}{T}_{2}q.$$This equation represents the family of the parallel lines of lower spatial frequency (whose period is longer than the period of either grating). From the equation of the straight line, Eq. (4) can be rewritten in the normal form Eq. (5)$$y\text{\hspace{0.17em}}\mathrm{sin}\text{\hspace{0.17em}}\theta +x\text{\hspace{0.17em}}\mathrm{cos}\text{\hspace{0.17em}}\theta p=0.$$Eq. (6)$$\mathrm{tan}\text{\hspace{0.17em}}\theta =\frac{{T}_{1}\text{\hspace{0.17em}}\mathrm{sin}\text{\hspace{0.17em}}\alpha}{{T}_{2}{T}_{1}\text{\hspace{0.17em}}\mathrm{cos}\text{\hspace{0.17em}}\alpha}=\frac{\mathrm{sin}\text{\hspace{0.17em}}\alpha}{\frac{{T}_{2}}{{T}_{1}}\mathrm{cos}\text{\hspace{0.17em}}\alpha},$$Eq. (7)$${T}_{m}=\frac{{T}_{1}{T}_{2}}{\sqrt{{T}_{1}^{2}2{T}_{1}{T}_{2}\text{\hspace{0.17em}}\mathrm{cos}\text{\hspace{0.17em}}\alpha +{T}_{2}^{2}}}=\frac{{T}_{2}}{\sqrt{12\frac{{T}_{2}}{{T}_{1}}\text{\hspace{0.17em}}\mathrm{cos}\text{\hspace{0.17em}}\alpha +{\left(\frac{{T}_{2}}{{T}_{1}}\right)}^{2}}}.$$In the case of the identical gratings (${T}_{1}={T}_{2}$), Eqs. (6) and (7) are simplified. The angle becomes Eq. (8)$$\mathrm{tan}\text{\hspace{0.17em}}\theta =\frac{\mathrm{sin}\text{\hspace{0.17em}}\alpha}{1\mathrm{cos}\text{\hspace{0.17em}}\alpha}=\mathrm{cot}\text{\hspace{0.17em}}\frac{\alpha}{2},$$Eq. (9)$${T}_{m}=\frac{{T}_{1}}{2\text{\hspace{0.17em}}\mathrm{sin}\text{\hspace{0.17em}}\frac{\alpha}{2}}.$$Equations (6) and (7) correspond to Eq. (2.9) from Ref. 1 with ${\theta}_{1}=0$, whereas Eqs. (8) and (9) correspond to Eq. (2.10) there. A visual illustration of the moiré effect in the identical gratings is shown in Fig. 6 ($\alpha =5\text{\hspace{0.17em}\hspace{0.17em}}\mathrm{deg}$). 2.1.Exercises
3.Sinusoidal Coplanar GratingsThe profile of the gratings also can be taken into account. For example, the analytical expressions of the moiré bands can be obtained for the sinusoidal gratings. The moiré patterns are the patterns of a longer period. The estimation of their visual appearance is based on the wave numbers. 3.1.Two Line GratingsAn interaction between the gratings can be modeled mathematically by the multiplication of transparency functions of the gratings. For instance, a onedimensional (1D) sinusoidal line grating (whose intensity profile along certain direction is a sinusoidal function) can be described by its transparency function $t=(1+\mathrm{cos}\text{\hspace{0.17em}}\mathbf{k}\xb7\mathbf{x})/2$, where $\mathbf{k}$ is the wave vector (inversely proportional to the period), and $\mathbf{x}$ is the coordinate axis. A sinusoidal grating is shown in Fig. 7(a). A result of the interaction between two gratings Fig. 7(a) can be written as Eq. (10)$${t}_{12}={t}_{1}{t}_{2}=\frac{1}{2}[1+\mathrm{cos}({\mathbf{k}}_{1}\xb7\mathbf{x})]\xb7\frac{1}{2}[1+\mathrm{cos}({\mathbf{k}}_{2}\xb7\mathbf{x})]\phantom{\rule{0ex}{0ex}}=\frac{1}{4}[1+\mathrm{cos}({\mathbf{k}}_{1}\xb7\mathbf{x})+\mathrm{cos}({\mathbf{k}}_{2}\xb7\mathbf{x})+\mathrm{cos}({\mathbf{k}}_{1}\xb7\mathbf{x})\mathrm{cos}({\mathbf{k}}_{2}\xb7\mathbf{x})]\phantom{\rule{0ex}{0ex}}=\frac{1}{4}\{1+\mathrm{cos}({\mathbf{k}}_{1}\xb7\mathbf{x})+\mathrm{cos}({\mathbf{k}}_{2}\xb7\mathbf{x})+\frac{1}{2}\text{\hspace{0.17em}}\mathrm{cos}[({\mathbf{k}}_{1}+{\mathbf{k}}_{2})\xb7\mathbf{x}]+\frac{1}{2}\text{\hspace{0.17em}}\mathrm{cos}[({\mathbf{k}}_{1}{\mathbf{k}}_{2})\xb7\mathbf{x}]\},$$For the moiré effect, we need the smallest wave number. In Eq. (10), the first term is the constant bias (a “DC” term), the second and third terms are the gratings themselves, whereas the fourth and fifth terms represent combinational spatial frequencies (sum and difference). In Eq. (10), the term with the smallest wave number (and with the lowest spatial frequency) is the fourth term proportional to $\mathrm{cos}[({\mathbf{k}}_{1}{\mathbf{k}}_{2})\xb7\mathbf{x}]$. 3.2.Square Grid + Line GratingSuperimposed gratings model a typical structure of many displays. A square grid can be thought as a pair of overlapped orthogonal gratings, whereas the barrier (or lenticular) plate can be modeled by a line grating. Accordingly, for two superposed rectangular gratings, we need two such pairs, four gratings total. For the square grid superposed with the line grating, we need three gratings. Examples of the grating and the grid are shown in Fig. 7. The transparency function of overlapped gratings obtained by the multiplication is Eq. (11)$${t}_{123}={t}_{1}{t}_{2}{t}_{3}=\frac{1}{16}[1+\mathrm{cos}({\mathbf{k}}_{1x}\xb7{\mathbf{x}}_{1})][1+\mathrm{cos}({\mathbf{k}}_{1y}\xb7{\mathbf{y}}_{1})][1+\mathrm{cos}({\mathbf{k}}_{2x}\xb7{\mathbf{x}}_{2})].$$The vector relationship Eq. (11) can be rewritten in a scalar form. For that purpose, the dot products in Eq. (11) can be rewritten as follows: Eq. (12)$${\mathbf{k}}_{1x}\xb7{\mathbf{x}}_{1}=\rho k(\mathrm{1,0})\xb7(x,y)=\rho kx\phantom{\rule{0ex}{0ex}}{\mathbf{k}}_{1y}\xb7{\mathbf{y}}_{1}=\rho k(\mathrm{0,1})\xb7(x,y)=\rho ky\phantom{\rule{0ex}{0ex}}{\mathbf{k}}_{2x}\xb7{\mathbf{x}}_{2}=k(\mathrm{cos}\text{\hspace{0.17em}}\alpha ,\mathrm{sin}\text{\hspace{0.17em}}\alpha )\xb7(x,y)=k(x\text{\hspace{0.17em}}\mathrm{cos}\text{\hspace{0.17em}}\alpha +y\text{\hspace{0.17em}}\mathrm{sin}\text{\hspace{0.17em}}\alpha ),$$Eq. (14)$${t}_{123}={t}_{1}{t}_{2}{t}_{3}=\frac{1}{16}(1+\mathrm{cos}\text{\hspace{0.17em}}\rho kx)(1+\mathrm{cos}\text{\hspace{0.17em}}\rho ky)\{1+\mathrm{cos}[k(x\text{\hspace{0.17em}}\mathrm{cos}\text{\hspace{0.17em}}\alpha +y\text{\hspace{0.17em}}\mathrm{sin}\text{\hspace{0.17em}}\alpha )]\}\phantom{\rule{0ex}{0ex}}=\frac{1}{16}(1+\mathrm{cos}\text{\hspace{0.17em}}\rho kx+\mathrm{cos}\text{\hspace{0.17em}}\rho ky+\mathrm{cos}\text{\hspace{0.17em}}\rho kx\text{\hspace{0.17em}}\mathrm{cos}\text{\hspace{0.17em}}\rho ky)\{1+\mathrm{cos}[k(x\text{\hspace{0.17em}}\mathrm{cos}\text{\hspace{0.17em}}\alpha +y\text{\hspace{0.17em}}\mathrm{sin}\text{\hspace{0.17em}}\alpha )]\}\phantom{\rule{0ex}{0ex}}=\frac{1}{16}[1+\mathrm{cos}\text{\hspace{0.17em}}\rho kx+\mathrm{cos}\text{\hspace{0.17em}}\rho ky+\frac{1}{2}\text{\hspace{0.17em}}\mathrm{cos}\text{\hspace{0.17em}}\rho (kxky)+\frac{1}{2}\text{\hspace{0.17em}}\mathrm{cos}\text{\hspace{0.17em}}\rho (kx+ky)]\{1+\mathrm{cos}[k(x\text{\hspace{0.17em}}\mathrm{cos}\text{\hspace{0.17em}}\alpha +y\text{\hspace{0.17em}}\mathrm{sin}\text{\hspace{0.17em}}\alpha )]\}.$$Similarly, applying the equation for the product of cosines several times, we have Eq. (15)$${t}_{123}=\frac{1}{16}[1+\mathrm{cos}\text{\hspace{0.17em}}\rho kx+\mathrm{cos}\text{\hspace{0.17em}}\rho ky+\mathrm{cos}(\mathrm{cos}\text{\hspace{0.17em}}\alpha kx+\mathrm{sin}\text{\hspace{0.17em}}\alpha ky)]\phantom{\rule{0ex}{0ex}}+\frac{1}{32}\left\{\begin{array}{l}\mathrm{cos}\text{\hspace{0.17em}}\rho k(xy)+\mathrm{cos}\text{\hspace{0.17em}}\rho k(x+y)\\ +\mathrm{cos}[(\rho +\mathrm{cos}\text{\hspace{0.17em}}\alpha )kx+\mathrm{sin}\text{\hspace{0.17em}}\alpha ky]+\mathrm{cos}[(\rho \mathrm{cos}\text{\hspace{0.17em}}\alpha )kx\mathrm{sin}\text{\hspace{0.17em}}\alpha ky]\\ +\mathrm{cos}[\mathrm{cos}\text{\hspace{0.17em}}\alpha kx+(\rho +\mathrm{sin}\text{\hspace{0.17em}}\alpha )ky]+\mathrm{cos}[\mathrm{cos}\text{\hspace{0.17em}}\alpha kx+(\rho \mathrm{sin}\text{\hspace{0.17em}}\alpha )ky]\end{array}\right\}\phantom{\rule{0ex}{0ex}}+\frac{1}{64}\left\{\begin{array}{l}\mathrm{cos}[(\rho +\mathrm{cos}\text{\hspace{0.17em}}\alpha )kx+(\rho +\mathrm{sin}\text{\hspace{0.17em}}\alpha )ky]+\mathrm{cos}[(\rho \mathrm{cos}\text{\hspace{0.17em}}\alpha )kx+(\rho +\mathrm{sin}\text{\hspace{0.17em}}\alpha )ky]\\ +\mathrm{cos}[(\rho +\mathrm{cos}\text{\hspace{0.17em}}\alpha )kx+(\rho +\mathrm{sin}\text{\hspace{0.17em}}\alpha )ky]+\mathrm{cos}[(\rho \mathrm{cos}\text{\hspace{0.17em}}\alpha )kx+(\rho \mathrm{sin}\text{\hspace{0.17em}}\alpha )ky]\end{array}\right\}.$$The following identity means the rotation of coordinates by the angle $\theta =\mathrm{arctan}\text{\hspace{0.17em}}b/a$ Eq. (16)$$\mathrm{cos}(ax+by)=\mathrm{cos}\left(\sqrt{{a}^{2}+{b}^{2}}{x}^{\prime}\right)=\mathrm{cos}(k{x}^{\prime}),$$Applying Eq. (17) to Eq. (15), we can obtain the wave numbers of all terms of Eq. (15). These wave numbers are listed below in the same order as in Eq. (15) Eq. (18)$${k}_{i}=0,\rho ,\rho ,1,\sqrt{2}\xb7\rho ,\sqrt{2}\xb7\rho ,\phantom{\rule[0.0ex]{1em}{0.0ex}}\sqrt{{\rho}^{2}+2\rho \text{\hspace{0.17em}}\mathrm{cos}\text{\hspace{0.17em}}\alpha +1},\phantom{\rule{0ex}{0ex}}\sqrt{{\rho}^{2}2\rho \text{\hspace{0.17em}}\mathrm{cos}\text{\hspace{0.17em}}\alpha +1},\phantom{\rule[0.0ex]{1em}{0.0ex}}\sqrt{{\rho}^{2}+2\rho \text{\hspace{0.17em}}\mathrm{sin}\text{\hspace{0.17em}}\alpha +1},\phantom{\rule[0.0ex]{1em}{0.0ex}}\sqrt{{\rho}^{2}2\rho \text{\hspace{0.17em}}\mathrm{sin}\text{\hspace{0.17em}}\alpha +1},\sqrt{2{\rho}^{2}+2\rho \text{\hspace{0.17em}}\mathrm{cos}\text{\hspace{0.17em}}\alpha 2\rho \text{\hspace{0.17em}}\mathrm{sin}\text{\hspace{0.17em}}\alpha +1},\phantom{\rule[0.0ex]{1em}{0.0ex}}\phantom{\rule{0ex}{0ex}}\sqrt{2{\rho}^{2}2\rho \text{\hspace{0.17em}}\mathrm{cos}\text{\hspace{0.17em}}\alpha +2\rho \text{\hspace{0.17em}}\mathrm{sin}\text{\hspace{0.17em}}\alpha +1},\sqrt{2{\rho}^{2}+2\rho \text{\hspace{0.17em}}\mathrm{cos}\text{\hspace{0.17em}}\alpha +2\rho \text{\hspace{0.17em}}\mathrm{sin}\text{\hspace{0.17em}}\alpha +1},\phantom{\rule[0.0ex]{1em}{0.0ex}}\sqrt{2{\rho}^{2}2\rho \text{\hspace{0.17em}}\mathrm{cos}\text{\hspace{0.17em}}\alpha 2\rho \text{\hspace{0.17em}}\mathrm{sin}\text{\hspace{0.17em}}\alpha +1},\phantom{\rule[0.0ex]{1em}{0.0ex}}i=1,\dots ,14.$$Equation (18) shows that there are eight different nonconstants (i.e., depending on the angle) wave numbers among 14 terms of Eq. (15). The symmetry of the problem suggests that the angular range (domain) to be considered in the case of the square grid is [0 deg, 45 deg]. All calculations can be made within this domain only; the angles outside it can be reduced into the domain by means of finding the remainder of the division by 45. None of eight nonconstant functions equals 0 within the domain. Only two functions fall down to 0 at the edges; the function $\sqrt{{\rho}^{2}2\rho \text{\hspace{0.17em}}\mathrm{cos}\text{\hspace{0.17em}}\alpha +1}=0$ at the left edge ($\alpha =0$) with ${\rho}_{01}=0$, and the function $\sqrt{2{\rho}^{2}2\rho \text{\hspace{0.17em}}\mathrm{cos}\text{\hspace{0.17em}}\alpha 2\rho \text{\hspace{0.17em}}\mathrm{sin}\text{\hspace{0.17em}}\alpha +1}$ at the right edge ($\alpha =45\text{\hspace{0.17em}\hspace{0.17em}}\mathrm{deg}$) with ${\rho}_{11}=1/\surd 2$. Both cases describe infinitely long waves, which represent the strongest moiré waves because according to our visibility assumption^{48} that the waves with longer periods are better visible (this assumption is confirmed later experimentally by the direct amplitude measurements^{58}). The corresponding wave numbers are Eq. (19)$${k}_{01}(\rho ,\alpha )=\sqrt{{\rho}^{2}2\rho \text{\hspace{0.17em}}\mathrm{cos}\text{\hspace{0.17em}}\alpha +1},$$Eq. (20)$${k}_{11}(\rho ,\alpha )=\sqrt{2{\rho}^{2}2\rho \text{\hspace{0.17em}}\mathrm{cos}\text{\hspace{0.17em}}\alpha 2\rho \text{\hspace{0.17em}}\mathrm{sin}\text{\hspace{0.17em}}\alpha +1}.$$The functions defined in Eqs. (19) and (20) are continuous and monotonous; one of them is rising, another falling. Therefore, they must intersect within the domain at certain ${\alpha}_{0}$, and the period at that point is definitely shorter than at the edges, and the moiré patterns are less visible. Thus, the moiré patterns can be minimized. In the left “half” of the domain ($0<\alpha <{\alpha}_{0}$), the function defined in Eq. (19) prevails, i.e., it has lower values which correspond to the longer wavelengths. In the right “half” of the domain (${\alpha}_{0}<\alpha <45\text{\hspace{0.17em}\hspace{0.17em}}\mathrm{deg}$), the function defined in Eq. (20) prevails. The best visible patterns at the left edge of the domain are represented by the first function. At larger angles, the visibility falls down (because the wavelength shortens) until the intersection point. After that point, the best visible patterns are represented by the second function; the visibility increases and reaches another maximum at the right edge of the domain. Therefore in general, the visibility of the moiré patterns is lowest at the intersection point ${\alpha}_{0}$. Examples of the moiré patterns are shown in Fig. 8; the corresponding gratings are shown in Fig. 7. Figure 8 shows the moiré patterns at the angles near the edges of the domain and near its middle for the size ratios of 1 and 0.707 (namely, the angles slightly deviated from 0 deg, 45 deg, and 27 deg). Note that the period of the patterns in Figs. 8(b)–8(e) is not much longer than the period of the gratings. Let us find the intersection of the functions defined by Eqs. (19) and (20). Substitute ${\rho}_{1}$ mentioned above into the first function, and ${\rho}_{2}$ into the second one. At the intersection, the two functions are equal, i.e., Eq. (21)$$\sqrt{12\text{\hspace{0.17em}}\mathrm{cos}\text{\hspace{0.17em}}\alpha +1}=\sqrt{1\frac{2}{\sqrt{2}}\text{\hspace{0.17em}}\mathrm{cos}\text{\hspace{0.17em}}\alpha \frac{2}{\sqrt{2}}\text{\hspace{0.17em}}\mathrm{sin}\text{\hspace{0.17em}}\alpha +1},$$Eq. (22)$$\mathrm{cos}\text{\hspace{0.17em}}\alpha (1\mathrm{sin}\text{\hspace{0.17em}}\alpha )=\frac{1}{2}.$$For the new variable $S=\mathrm{sin}\text{\hspace{0.17em}}\alpha $, Eq. (22) can be rewritten as a quartic equation This equation has the only real root $S=0.442$ within the domain; the corresponding angle is Note that this angle is very close to Eq. (25)$${\alpha}_{1/2}=\mathrm{arctan}\text{\hspace{0.17em}}\frac{1}{2}=26.565\text{\hspace{0.17em}\hspace{0.17em}}\mathrm{deg}.$$From a proper (long enough) distance, the shortperiod patterns at the angle ${\alpha}_{0}$ (or ${\alpha}_{1/2}$) are unrecognizable and thus are effectively eliminated, see Figs. 8(b) and 8(e). For the equations in the polar coordinates refer to Ref. 59. 3.2.1.Exercises
4.Sinusoidal Noncoplanar GratingsPractically, the layers may not lie in the same plane. For example, a typical structure of the autostereoscopic 3D displays contains two parallel layers separated by an air gap. Such a nonplanar layout makes the moiré patterns alive and vivid. The visual picture may look dissimilar from different directions, and the movement of the patterns sometimes looks unexpected; moreover, their behavior may seem unpredictable. This section is based on Refs. 28 and 60. Consider two layers observed from a finite distance $l$ as shown in Fig. 9. According to Ref. 60, the period of the moiré patterns in this case (two line gratings) is where $\rho $ is defined in Eq. (13) and $s$ is the geometric factor where $l$ is the coordinate of the observer and $\mathrm{\Delta}l$ is the gap between the layers.This expression for the period can be rewritten in terms of the moiré magnification factor; by its physical meaning, the moiré factor is an amplification coefficient (or a gain factor) for the period of the grating. The moiré factor corresponding to Eq. (26) is Note that when $\rho =s$, the moiré factor Eq. (28) formally reaches infinity, which means the infinitely long period. However practically, we cannot measure an infinite period. The maximum value of the period is always limited by the size of a screen, where the patterns are observed; this is because the patterns with the period longer than that screen cannot be recognized as periodic waves at all. Equation (28) is graphically shown in Fig. 10. The particular case of Eq. (28) for the identical gratings with $\rho =1$ is as follows: Examples of moiré patterns are shown in Fig. 11 for the same distance, but different gaps between gratings in Figs. 11(a) and 11(b), as well as for the same gap but different distances in Figs. 11(c) and 11(d). A lateral displacement of the gratings (or the camera) does not affect the period; however in this case, the patterns are displaced laterally by the following value: where ${x}_{1}$ and ${x}_{2}$ are the displacements of the gratings. The physical meaning of this equation is that the displacement of a grating by its period results in the displacement of the patterns by their period.Note that our equations were obtained for the camera axis orthogonal to the plain gratings. The picture may look completely different in an nonorthogonal case. Compare two photographs of the same plane pavement in Fig. 12 taken by the same camera, but with the different angles between the camera axis and the surface. In the first case (the camera axis perpendicular to the surface), the shape and size of all bricks (rectangles) in the photograph are identical. In the second case (the camera with the axis deviated from the perpendicular), all quadrilaterals (projections of the rectangular bricks) are different; it means a strong dependence of their shape and size on coordinates. 4.1.Exercises
5.Cylindrical Moiré (Orthogonal Projection)Consider a regular 3D object: a cylinder made of a wrapped periodic mesh, see Fig. 13. This section is based on Refs. 6162.–63. The cylinder can be logically split into two halves, front and rear, which are made of the same mesh. When we look through meshes with a gap between them, the moiré patterns may appear, as described in Sec. 4. To calculate the period of the moiré patterns, we use the orthogonal projection of the points lying on the surface of the cylinder onto a screen which is parallel to the $xy$plane. Correspondingly, the projection lines are always parallel to the $z$axis; the coordinate $x$ can be treated as an impact parameter. Let $R$ be the radius of the cylinder, $L$ the distance from the screen $S$ to the center of the circle. An example projection line in Fig. 14 connects three circular dots, two on the circle, and one on the screen. The point, where this line crosses the $x$axis, is the orthogonal projection of two points from the cylindrical surface onto the screen $S$. In this section, all distances are measured along the lines parallel to the $z$axis. The distance $\mathrm{\Delta}l$ between the gratings along the $z$coordinate displaced by $x$ is equal to the chord of the circle. From the equation of the cylinder, the length of the chord is Then, the distance $l$ from the observer to the first grating is $L$ minus one half of the chord With using the polar angle $\phi =\mathrm{arcsin}(x/R)$ of polar coordinates, Eqs. (31) and (32) can be reexpressed as follows: The variables $l$ and $\mathrm{\Delta}l$ for Eq. (29) describing the moiré effect in the parallel noncoplanar gratings are found. However in the case of a spatial object as a cylinder, the gratings are not parallel, and therefore the orientation of the wave vector in space has to be additionally taken into account. Correspondingly, two types of gratings on the surface of the cylinder should be analyzed: the grating with the vertical wave vector orthogonal to the base of the cylinder shown in Fig. 15(a) and the grating with the horizontal wave vector parallel to the base shown in Fig. 15(b). In the other words, the former grating consists of the identical circles uniformly displaced along the $y$axis; the latter consists of the vertical lines uniformly distributed along the circle, the base of the cylinder. These two cases are well separated (this fact is confirmed experimentally in Ref. 61) and therefore can be considered in isolation from one another. In the orthogonal projection of the first grating onto the $yz$plane, the projected period is equal to the period of the grating, see Fig. 16(a). However, the projected period of the second grating is not constant along the $x$axis, as shown in Fig. 16(b), and depends on the local inclination of the cylindrical surface, which is proportional to the cosine of the polar angle. Therefore in this case, the projected period should be multiplied by $\mathrm{cos}\text{\hspace{0.17em}}\phi $. Based on Eqs. (33) and (34), the moiré periods of these cases are calculated as follows: Eq. (35)$${T}_{v}=\frac{l}{\mathrm{\Delta}l}=\frac{T}{2}(\frac{L}{R\text{\hspace{0.17em}}\mathrm{cos}\text{\hspace{0.17em}}\phi}1),$$Eq. (36)$${T}_{h}=\frac{l}{\mathrm{\Delta}l}\text{\hspace{0.17em}}\mathrm{cos}\text{\hspace{0.17em}}\phi =\frac{T}{2}(\frac{L}{R}\mathrm{cos}\text{\hspace{0.17em}}\phi ).$$The photographs of the moiré patterns in cylinders with the horizontal and vertical wave vectors are shown in Fig. 17. The corresponding moiré magnification factors are given by the following equations: Eq. (37)$${\mu}_{v}=\frac{1}{2}(\frac{L}{R\text{\hspace{0.17em}}\mathrm{cos}\text{\hspace{0.17em}}\phi}1),$$A graphical illustration of the theoretical dependences Eqs. (37) and (38) is shown in Fig. 18. Note that $\mathrm{cos}\text{\hspace{0.17em}}\phi =x/R$. Across the radius, the horizontal moiré factor deviates $<1\%$, whereas the vertical one raises more than twice. 5.1.Exercises
6.Moiré Wave Vector as a Vector SumNow, consider the moiré effect in the spectral domain, where we deal with the reciprocal distances or, equivalently, with the wave vectors. Generally speaking, the wave vector of the moiré patterns is that linear combination of the wave vectors of gratings (with coefficients $+1$ and $1$), which is closest to the origin. It means either the sum or the difference; whichever is shorter, see Fig. 19. Typically (but not always), it is the difference between the wave vectors of two gratings. Let ${\mathbf{k}}_{1}$ and ${\mathbf{k}}_{2}$ be the wave vectors of the gratings and $\alpha $ the angle between them. Based on these vectors, we will find the moiré wave vector ${\mathbf{k}}_{\mathrm{m}}$ and the moiré orientation angle $\theta $. Figure 19 is an illustration for the wave vectors with close wave numbers but different orientations. It can be seen that for close wave numbers and $\alpha >90\text{\hspace{0.17em}\hspace{0.17em}}\mathrm{deg}$, it is the summation; otherwise (i.e., when $\alpha <90\text{\hspace{0.17em}\hspace{0.17em}}\mathrm{deg}$), it is the subtraction. Figure 19(b) actually shows the layout of the wave vectors for Eq. (10) in Sec. 3. Figure 20 shows the map of the spectrum of two overlapped sinusoidal gratings [Fig. 19(b)]. Among all combinational components, the linear combination with the smallest wave number is the moiré wave vector. The extremely important concept of the visibility circle^{1} models the human visual system in the spectral domain. The visible vectors lie within the visibility circle, the vectors outside the visibility circle are invisible, see Fig. 21. Consider the triangle with the sides ${\mathbf{k}}_{1}$ and ${\mathbf{k}}_{m}$ and the angle $2\pi \theta $ between them in the case of the subtraction shown in Fig. 19(b). From the law of cosines for the side ${\mathbf{k}}_{2}$ Eq. (39)$${k}_{2}^{2}={k}_{1}^{2}+{k}_{m}^{2}2{k}_{1}{k}_{m}\text{\hspace{0.17em}}\mathrm{cos}(2\pi \theta )$$Eq. (40)$$\mathrm{cos}(2\pi \theta )=\frac{{k}_{1}^{2}+{k}_{m}^{2}{k}_{2}^{2}}{2{k}_{1}{k}_{m}}.$$From the law of cosines for the side ${\mathbf{k}}_{m}$, we have the modulus of the moiré wave vector Eq. (41)$${k}_{m}^{2}={k}_{1}^{2}+{k}_{2}^{2}2{k}_{1}{k}_{2}\text{\hspace{0.17em}}\mathrm{cos}\text{\hspace{0.17em}}\alpha $$Eq. (42)$$\mathrm{cos}(2\pi \theta )=\frac{{k}_{1}{k}_{2}\text{\hspace{0.17em}}\mathrm{cos}\text{\hspace{0.17em}}\alpha}{\sqrt{{k}_{1}^{2}+{k}_{2}^{2}2{k}_{1}{k}_{2}\text{\hspace{0.17em}}\mathrm{cos}\text{\hspace{0.17em}}\alpha}}.$$Therefore, Eq. (43)$$\mathrm{sin}(2\pi \theta )=\sqrt{1{\mathrm{cos}}^{2}(2\pi \theta )}=\frac{{k}_{2}\text{\hspace{0.17em}}\mathrm{sin}\text{\hspace{0.17em}}\alpha}{\sqrt{{k}_{1}^{2}+{k}_{2}^{2}2{k}_{1}{k}_{2}\text{\hspace{0.17em}}\mathrm{cos}\text{\hspace{0.17em}}\alpha}}.$$Recall that Thus, the orientation of the wave vector of the moiré patterns is Eq. (46)$$\mathrm{tan}\text{\hspace{0.17em}}\theta =\frac{\mathrm{sin}\text{\hspace{0.17em}}\theta}{\mathrm{cos}\text{\hspace{0.17em}}\theta}=\frac{\mathrm{sin}(2\pi \theta )}{\mathrm{cos}(2\pi \theta )}=\frac{{k}_{2}\text{\hspace{0.17em}}\mathrm{sin}\text{\hspace{0.17em}}\alpha}{{k}_{1}{k}_{2}\text{\hspace{0.17em}}\mathrm{cos}\text{\hspace{0.17em}}\alpha}.$$The two Eqs. (41) and (46) comprise the full solution (the wave number and the orientation). The relation for the periods is an inverse wave vector Eq. (41) Eq. (47)$${T}_{m}=\frac{{T}_{1}{T}_{2}}{\sqrt{{T}_{1}^{2}+{T}_{2}^{2}2{T}_{1}{T}_{2}\text{\hspace{0.17em}}\mathrm{cos}\text{\hspace{0.17em}}\alpha}}.$$The orientation Eq. (46) can be also rewritten in terms of the periods Eq. (48)$$\mathrm{tan}\text{\hspace{0.17em}}\theta =\frac{\frac{1}{{T}_{2}}\text{\hspace{0.17em}}\mathrm{sin}\text{\hspace{0.17em}}\alpha}{\frac{1}{{T}_{1}}\frac{1}{{T}_{2}}\text{\hspace{0.17em}}\mathrm{cos}\text{\hspace{0.17em}}\alpha}=\frac{{T}_{1}\text{\hspace{0.17em}}\mathrm{sin}\text{\hspace{0.17em}}\alpha}{{T}_{2}{T}_{1}\text{\hspace{0.17em}}\mathrm{cos}\text{\hspace{0.17em}}\alpha}.$$This is actually the same expression as Eq. (6) in Sec. 2. Moreover, Eqs. (47) and (48) can be reexpressed in terms of the ratio of periods $\rho ={T}_{2}/{T}_{1}$ [compare with the definition Eq. (13) in Sec. 3.2] as follows: Eq. (49)$${T}_{m}={T}_{2}\frac{1}{\sqrt{1+{\rho}^{2}2\rho \text{\hspace{0.17em}}\mathrm{cos}\text{\hspace{0.17em}}\alpha}},$$Eq. (50)$$\mathrm{tan}\text{\hspace{0.17em}}\theta =\frac{\mathrm{sin}\text{\hspace{0.17em}}\alpha}{\rho \mathrm{cos}\text{\hspace{0.17em}}\alpha}.$$Furthermore, Eq. (49) can be expressed in terms of the moiré factor Eq. (51)$$\mu =\frac{1}{\sqrt{1+{\rho}^{2}2\rho \text{\hspace{0.17em}}\mathrm{cos}\text{\hspace{0.17em}}\alpha}}.$$6.1.Exercises
7.Basics of TwoDimensional Fourier Transform7.1.SpectraA relation between variables can be described mathematically in two ways: either by the functional dependence or by the spectrum; corresponding examples are given in Fig. 22. Both ways characterize the same relation from different perspectives. For example, Figs. 22(a) and 22(b) describe the sinusoidal wave, whereas Figs. 22(c) and 22(d) the square wave. The Fourier spectra of real symmetric functions are real and symmetric, although in general, the spectra of arbitrary functions (including the real but nonsymmetric functions) are complex. In the tutorial we consider the power spectra, i.e., the modules of the Fourier coefficients, which are always real. For the periodic sinusoidal ${f}_{1}$ and rectangular ${f}_{2}$ functions shown in Fig. 22, we have Eq. (53)$$\{\begin{array}{l}{f}_{2}(x)=\{\begin{array}{ll}1,& x<1/2\\ 0,& 1/2<x<1\end{array}\\ {f}_{2}(x+T)={f}_{2}(x)\end{array},$$Eq. (55)$${F}_{2}(k)=\frac{1}{2}+\frac{2}{\pi}\sum _{n=\mathrm{1,3},5,\dots}^{\infty}\frac{1}{n}\text{\hspace{0.17em}}\mathrm{sin}(2\pi nk).$$The previous expression can be rewritten in terms of sinc function defined as follows: (Sometimes you can meet an alternative definition $\mathrm{sinc}(x)=\mathrm{sin}(\pi x)/\pi x$ which differs in the normalization of the coordinate.) The rewritten Eq. (55) is as follows:Eq. (57)$${F}_{2}(k)=\frac{1}{2}+4k\sum _{n=\mathrm{1,3},5,\dots}^{\infty}\text{\hspace{0.17em}}\mathrm{sinc}(2\pi nk).$$The sinusoidal grating has three spectral components (one of them is a constant term, while two others represent a sinusoidal wave itself). The rectangular grating has many spectral components; theoretically, an infinite number of the decayed components. The decay rate of the Fourier coefficients depends on the smoothness of the function^{64} and particularly, for a piecewise continuous function is $1/n$. The spectrum of a symmetric square wave contains only odd harmonic frequencies, see Fig. 22(d), where all even harmonics are equal to zero. Figure 22(b) shows that the power spectrum of a sinusoidal grating can be thought as a limited (central) part of the spectrum of a rectangular grating. 7.2.TwoDimensional SpectraFigure 22 of the previous section shows the 1D functions. A 1D function of one variable (graphically, the “height” $y$ as a function of the abscissa $x$) can be generalized to a 2D function of two variables (the “height” $z$ as a function of two independent coordinates $x$ and $y$). The 2D case looks somewhat complicated, but not very sophisticated and still understandable. In two dimensions, the spectrum of a line grating is spread along a slant straight line (an abscissa of a 1D spectrum) and repeated (copied) in the orthogonal direction. A 2D function represents a surface and can be displayed, for instance, as a “map” with colors for the height; the white color may mean the lowest value (say, zero), while the black color means the highest value (say, one). Such maps of the plane waves (sinusoidal and square profile) are drawn in Figs. 23(a) and 23(c) together with their profiles; the spectra are shown in Figs. 23(b) and 23(d), respectively. In many cases, a 2D grid can be represented as a product of two 1D gratings. For illustration, refer to Fig. 7 and to Eqs. (10) and (11). Many useful details about discrete transforms can be found in Ref. 56. 7.2.1.Exercises
8.Spectral TrajectoriesGenerally speaking, parameters of the gratings do not always remain constant and may change. This change causes the change of the spectrum. In the case of an incremental change of a parameter, a set of several spectra represents a richer picture of the behavior of the patterns. For example, when a sinusoidal line grating whose initial spectrum is shown in Fig. 24(a) turned around, its spectrum is also rotated, see Fig. 24(b). Several overlapped spectra of an incrementally rotated grating are shown in Fig. 24(c) for the angles 0 deg to 40 deg with increment 5 deg; these look like a set of discrete points. Schematically, such sets of spectra can be drawn by continuous trajectories in the spectral domain, see Fig. 24(d). Furthermore, one may consider the spectral domain with the axes $u$, $v$ as the complex plane with $u$ and $v$ as the real and the imaginary parts of a complex number $z$. Then, all locations and shapes in the complex plane and spectral domain are identical, but the calculations become convenient and simple because the complex numbers allow applying many powerful mathematical theorems. This section is based on Ref. 57. Consider a superposition of two rectangular grids, which can be often met in practice. The equation of the spectral peaks of four gratings arranged in two layers (each layer consists of two orthogonal gratings) is Eq. (58)$${T}_{2\times 2}=k({p}_{1}{\sigma}_{1}+i{p}_{2})+k\rho ({p}_{3}{\sigma}_{3}+i{p}_{4}){e}^{i\alpha},$$A picture of trajectories Fig. 24(d) shows where the spectral peaks can be located when the angle varies. The spectral trajectories are known in many areas, refer, for instance, to Refs. 65 and 66. For the moiré effect, the spectral trajectories were first time proposed in Ref. 57. In the case of two square gratings and the running angle, the trajectory derived from Eq. (58) is as follows: It can be proven that the trajectories [Eq. (59)] are either circular arcs or segments of straight lines. Different dimensions of gratings (implying that a line grating is a 1D structure, the square grid is a 2D structure) yield two particular cases The trajectories for two line gratings and two square grids [Eqs. (60) and (59)] are shown in Fig. 25. These trajectories were observed in experiments.^{57} Many examples of trajectories of the sinusoidal gratings for other running parameters can be found in Ref. 67. For the moiré effect, it is important that there can be trajectories leaving the visibility circle, approaching, entering, or crossing it, as well as the trajectories always outside or always inside the visibility circle. The visibility circle is shown in Fig. 25 by thin dashed line. 8.1.Exercises
9.Visual Effects on the Move and on the RotationIn the case of a laterally moved 1D grating, the corresponding displacement of the visible moiré patterns is also lateral and given by Eq. (30) in Sec. 4. Similarly, the equation for displacement of the visible moiré patterns in the case of the moved observer can be derived from Ref. 60 as follows: In both cases (the moved grating or the moved observer), the displacement of the moved object by one of its periods causes the displacement of the visual picture by one period of the moiré patterns. This visual picture repeats periodically. Correspondingly, the laterally moved observer will repeatedly see the same visual picture at each period of the patterns. In the case of identical gratings ($\rho =1$), we have from Eq. (62) Equation (63) represents socalled moiré mirror effect in the identical gratings, which results in the displacement of the moiré patterns equal to the displacement of an observer; i.e., the patterns literally follow the observer’s movement, as his/her reflection in a plain mirror, see Fig. 26. When the gratings are turned around, i.e., the angle between the gratings changes, the period of the moiré patterns also changes. In the case of the sinusoidal square gratings, almost certainly there could be two maxima at 0 and at 45 deg. For nonsinusoidal gratings, there could be several maxima at the intermediate rational angles (whose tangents are rational numbers). At the maxima, the axis of the moiré patterns is parallel to the axis of the rotated grating and their period is maximal.^{68} This can be explained in the following way. Although the spectral trajectory passes the neighborhood of the origin, the wave number reaches a minimum at that point, where the trajectory crosses the line connecting the origin and the center of the trajectory. According to Eq. (26), the period remains finite for any relation between the parameters $\rho $ and $s$, except for the case of $\rho =s$. In other words, the maximum period characterizes the relation between $\rho $ and $s$. From this perspective, the value of the maximum magnification factor gives an estimate of this relation. Theoretically, the maximum period depends on the particular ratio of periods of the gratings, and for the coplanar gratings with the integer ratios, the period of the moiré patterns is infinite. There can exist several local maxima. For illustration, refer to the computer simulation.^{69} The moiré factor can only be infinite in the gratings with parallel wave vectors, when $\rho =s$, i.e., orThis means that to obtain the infinite moiré factor, the size should correspond to the distance and gap. Before and after the “critical” distance, the moiré factor is finite. When a longdistance observer approaches to the gratings, the moiré factor increases, then reaches infinity at the critical distance, and finally decreases. In the identical gratings at zero angle ($\alpha =0$), the patterns are parallel to the gratings and the moiré factor is equal to $z/d$. In identical gratings ($\rho =1$) with an arbitrary angle, both moiré factor and orientation vary as follows: Eq. (66)$$\mu =\frac{1}{\sqrt{1+{s}^{2}2s\text{\hspace{0.17em}}\mathrm{cos}\text{\hspace{0.17em}}\alpha}},$$Eq. (67)$$\mathrm{tan}\text{\hspace{0.17em}}\phi =\frac{s\text{\hspace{0.17em}}\mathrm{sin}\text{\hspace{0.17em}}\alpha}{s\text{\hspace{0.17em}}\mathrm{cos}\text{\hspace{0.17em}}\alpha 1}.$$These equations mean that at a long distance, the moiré factor can be large, depending on the angle. When the observer approaches, the moiré factor drops down practically to zero, always remaining finite. On that move (approach), the moiré patterns rotate by $\sim 90\text{\hspace{0.17em}\hspace{0.17em}}\mathrm{deg}$ from almost orthogonal orientation to the parallel one. These situations can be graphically illustrated as follows. The moiré factor in the identical gratings as a function of the distance under two conditions: $d=\mathrm{const}$ with $\alpha $ parameter is shown in Fig. 27; the graphs for $\alpha =\mathrm{const}$ with $d$ parameter look similar. In nonidentical gratings, the moiré factor theoretically can reach the infinity. The cases of nonidentical gratings (with $\alpha $ as parameter) and the parallel gratings (with $\rho $ as parameter) installed at the gap $d=1$ are shown in Fig. 28. In the latter case, the moiré orientation is unchanged because practically we cannot distinguish between the directions 0 deg and 180 deg. The moiré factor and the orientation of the moiré patterns in the nonidentical gratings installed at a small angle are shown in Fig. 29. 9.1.Exercises
10.Discussion and ConclusionThe tutorial covers several approaches to understand the moiré effect, to study it, and, particularly, to obtain the characteristics of the visible moiré patterns in the spatial and spectral domains. The physical meaning of equations is explained. The cross references between sections are made. The experimental evidences are provided in figures. The following topics are covered: the indicial equation, the moiré wave vector, the sinusoidal coplanar and noncoplanar gratings, the moiré effect in the cylinder, as well as the 2D Fourier transform, the moiré spectra, and the spectral trajectories. The visual effects in the displaced or rotated gratings are discussed. These topics collected in one article describe the moiré effect from various perspectives in a variety of scenarios. This gives readers a flexible opportunity to find solutions of practical problems using this or that approach. For further reading, we would like to provide some additional references to the moiré effect in displays.^{70}^{–}^{72} Amazing, but the moiré effect, traditionally considered as a negative visual effect in displays, can be used to generate images including 3D, refer to Refs. 38, 39, and 73. The color moiré effect is not considered in the tutorial. However, somebody interested can continue reading (Refs. 7475.76.77.–78). As far as the gap effect is concerned, the paper^{28} is already mentioned. Some more papers on the gap effect must be mentioned, too, Refs. 7980.–81. The related moiré rotation effect is considered in Ref. 82. Many sources on the moiré art are available.^{83}^{–}^{85} Among the books related to the moiré art, we would highly recommend the book,^{4} as well as a good interactive illustration.^{86} There are several modern painters and artists such as A. Minini, P. Dickens, P. Decrauzat, and C. CruzDiez, who use the moiré effect in their works; many of them are presented on the websites.^{87}^{–}^{90} The tutorial is intended for a wide audience, from beginners to specialists. Someone discovers the beauty of the moiré effect; someone finds details of the patterns, somebody else reveals a new approach. The authors believe that the tutorial can be useful for everybody who would read it. AcknowledgmentsThis work was partially supported by the CrossMinistry Giga KOREA Project through the Korean Government, Ministry of Science and ICT (MSIT); Grant No. GK17D0200, Development of Super MultiView Display Providing RealTime Interaction. ReferencesI. Amidror, The Theory of the Moiré Phenomenon, Volume I: Periodic Layers, 2nd ed.SpringerVerlag, London, United Kingdom
(2009). Google Scholar
K. Patorsky and M. Kujawinska, Handbook of the Moiré Fringe Technique, Elsevier, Amsterdam, The Netherland
(1993). Google Scholar
C. A. Walker, Handbook of Moiré Measurement, IOP Publishing, Bristol, Pennsylvania
(2004). Google Scholar
N. Carsten, Moiré Index, Gestalten Verlag, Berlin
(2010). Google Scholar
O. Bryngdahl,
“Moiré: formation and interpretation,”
J. Opt. Soc. Am., 64
(10), 1287
–1294
(1974). https://doi.org/10.1364/JOSA.64.001287 JOSAAH 00303941 Google Scholar
S. Yokozeki,
“Theoretical interpretation of the moiré pattern,”
Opt. Commun., 11
(4), 378
–381
(1974). https://doi.org/10.1016/00304018(74)902405 OPCOB8 00304018 Google Scholar
S. Yokozeki, Y. Kusaka and K. Patorski,
“Geometric parameters of moiré fringes,”
Appl. Opt., 15
(9), 2223
–2227
(1976). https://doi.org/10.1364/AO.15.002223 APOPAI 00036935 Google Scholar
C. A. Sciammarella,
“The moiré method—a review,”
Exp. Mech., 22
(11), 418
–433
(1982). https://doi.org/10.1007/BF02326823 EXMCAZ 00144851 Google Scholar
P. A. M. dos Santos,
“Moirélike patterns as a spatial beat frequency in photorefractive sinusoidal phase gratings superposition,”
Opt. Commun., 212
(4–6), 211
–216
(2002). https://doi.org/10.1016/S00304018(02)020175 OPCOB8 00304018 Google Scholar
O. Kafri and I. Glatt, The Physics of Moiré Metrology, Wiley, New York
(1990). Google Scholar
M. C. Hutley et al.,
“Moiré magnifier,”
Pure Appl. Opt., 3 133
–142
(1994). https://doi.org/10.1088/09639659/3/2/006 PAOAE3 09639659 Google Scholar
M. Dohnal,
“Moiré in a scanned image,”
Proc. SPIE, 4016 166
–170
(1999). https://doi.org/10.1117/12.373616 PSISDG 0277786X Google Scholar
Y. Kim et al.,
“Color moiré pattern simulation and analysis in threedimensional integral imaging for finding the moiréreduced tilted angle of a lens array,”
Appl. Opt., 48
(11), 2178
–2187
(2009). https://doi.org/10.1364/AO.48.002178 APOPAI 00036935 Google Scholar
J.Y. Son et al.,
“Moiré pattern reduction in fullparallax autostereoscopic imaging systems using two crossed lenticular plates as a viewing zone forming optics,”
in Proc. 10th Int. Display Research Workshops (IDW),
1401
–1404
(2003). Google Scholar
J.Y. Son and B. Javidi,
“Threedimensional imaging systems based on multiview images,”
J. Disp. Technol., 1
(1), 125
–140
(2005). https://doi.org/10.1109/JDT.2005.853354 IJDTAL 1551319X Google Scholar
A. Stern and B. Javidi,
“Threedimensional image sensing, visualization, and processing using integral imaging,”
Proc. IEEE, 94
(3), 591
–607
(2006). https://doi.org/10.1109/JPROC.2006.870696 IEEPAD 00189219 Google Scholar
N. S. Holliman,
“Threedimensional display systems,”
Handbook of Optoelectronics, 1067
–1100 Taylor & Francis, London, United Kingdom
(2006). Google Scholar
G. Bell, R. Craig and T. Simmiss,
“Moiré interference in multilayered displays,”
J. Soc. Inf. Disp., 15
(11), 883
–888
(2007). https://doi.org/10.1889/1.2812988 JSIDE8 07341768 Google Scholar
R. MartinezCuenca et al.,
“Progress in 3D multiperspective display by integral imaging,”
Proc. IEEE, 97
(6), 1067
–1077
(2009). https://doi.org/10.1109/JPROC.2009.2016816 IEEPAD 00189219 Google Scholar
J.Y. Son et al.,
“Recent developments in 3D imaging technologies,”
J. Disp. Technol., 6
(10), 394
–403
(2010). https://doi.org/10.1109/JDT.2010.2045636 IJDTAL 1551319X Google Scholar
J. Hong et al.,
“Threedimensional display technologies of recent interest: principles, status, and issues,”
Appl. Opt., 50
(34), H87
–H115
(2011). https://doi.org/10.1364/AO.50.000H87 APOPAI 00036935 Google Scholar
L. Kong, G. Jin and T. Wang,
“Analysis of moiré minimization in autostereoscopic parallax displays,”
Opt. Express, 21
(22), 26068
–26079
(2013). https://doi.org/10.1364/OE.21.026068 OPEXFF 10944087 Google Scholar
G. A. Bassett, J. W. Menter and D. W. Pashley,
“Moiré patterns on electron micrographs, and their application to the study of dislocations in metals,”
Proc. R. Soc. London Ser. A, 246
(1246), 345
–368
(1958). https://doi.org/10.1098/rspa.1958.0144 Google Scholar
Y. Hongbing et al.,
“Measurement on infrared radiation burning temperature field with moiré deviation tomography,”
Microwave Opt. Technol. Lett., 49
(11), 2761
–2763
(2007). https://doi.org/10.1002/(ISSN)10982760 Google Scholar
U. Bonse, W. Grae and G. Materlik,
“Xray interferometry and lattice parameter investigation,”
Rev. Phys. Appl., 11
(1), 83
–87
(1976). https://doi.org/10.1051/rphysap:0197600110108300 RPHAAN 00351687 Google Scholar
D. L. Miller et al.,
“Structural analysis of multilayer graphene via atomic moiré interferometry,”
Phys. Rev. B, 81 125427
(2010). https://doi.org/10.1103/PhysRevB.81.125427 Google Scholar
G. L. Rogers,
“A simple method of calculating moiré patterns,”
Proc. Phys. Soc., 73
(1), 142
–144
(1959). https://doi.org/10.1088/03701328/73/1/129 Google Scholar
C. A. Sciammarella and F.P. Chiang,
“Gap effect on moiré patterns,”
Z. Angew. Math. Phys., 19
(2), 326
–333
(1968). https://doi.org/10.1007/BF01601476 Google Scholar
O. Bryngdahl,
“Characteristics of superposed patterns in optics,”
J. Opt. Soc. Am., 66
(2), 87
–94
(1976). https://doi.org/10.1364/JOSA.66.000087 JOSAAH 00303941 Google Scholar
K. Patorski, S. Yokozeki and T. Suzuki,
“Moiré profile prediction by using Fourier series formalism,”
Jpn. J. Appl. Phys., 15
(3), 443
–456
(1976). https://doi.org/10.1143/JJAP.15.443 Google Scholar
S. Yokozeki,
“Moiré fringes,”
Opt. Laser Eng., 3
(1), 15
–27
(1982). https://doi.org/10.1016/01438166(82)900124 Google Scholar
C. A. Sciammarella,
“Basic optical law in the interpretation of moiré patterns applied to the analysis of strains—Part 1,”
Exp. Mech., 5
(5), 154
–160
(1965). https://doi.org/10.1007/BF02324049 EXMCAZ 00144851 Google Scholar
P. S. Theocaris,
“Moiré topography of curved surfaces,”
Exp. Mech., 7
(7), 289
–296
(1967). https://doi.org/10.1007/BF02327134 EXMCAZ 00144851 Google Scholar
P. S. Theocaris and S. A. Paipetis,
“Vibration measurements by the moiré method,”
J. Phys. E Sci. Instrum., 5
(3), 217
–219
(1972). https://doi.org/10.1088/00223735/5/3/010 Google Scholar
K. Creath, J. C. Wyant,
“Moiré and fringe projection techniques,”
Optical Shop Testing, 653
–685 John Wiley & Sons, New York
(1992). Google Scholar
J. Bartl, R. Fira and M. Hain,
“Inspection of surface by the moiré method,”
Meas. Sci. Rev., 1
(3), 29
–32
(2001). Google Scholar
G. N. de Oliveira, M. E. de Oliveira and P. A. M. dos Santos,
“Dynamic moiré patterns for profilometry applications,”
J. Phys. Conf. Ser., 274
(1), 012036
(2011). https://doi.org/10.1088/17426596/274/1/012036 Google Scholar
V. Saveljev and S.K. Kim,
“A 3D moiré display,”
in Proc. 33rd Int. Display Research Conf. (EuroDisplay),
54
–56
(2013). Google Scholar
V. Saveljev and S.K. Kim,
“Threedimensional moiré display,”
J. Soc. Inf. Disp., 22
(9), 482
–486
(2014). https://doi.org/10.1002/jsid.v22.9 JSIDE8 07341768 Google Scholar
H. Liu, A. N. Cartwright and C. Basaran,
“Sensitivity improvement in phaseshifted moiré interferometry using 1D continuous wavelet transform image processing,”
Opt. Eng., 42
(9), 2646
–2652
(2003). https://doi.org/10.1117/1.1592803 Google Scholar
A. Martínez, R. RodríguezVera and J. A. Rayas,
“Influence of object roughness on specimen gratings for moiré interferometry,”
Opt. Eng., 40
(9), 1978
–1983
(2001). https://doi.org/10.1117/1.1394738 Google Scholar
L. Sałbut, K. Patorski and M. Kujawiñska,
“Polarization approach to highsensitivity moiré interferometry,”
Opt. Eng., 31
(3), 434
–439
(1992). https://doi.org/10.1117/12.56108 Google Scholar
O. Kafri and I. Glatt,
“Moiré deflectometry: a ray deflection approach to optical testing,”
Opt. Eng., 24
(6), 944
–960
(1985). https://doi.org/10.1117/12.7973607 Google Scholar
W. Chen et al.,
“Extended dualgrating alignment method for optical projection lithography,”
Appl. Opt., 49
(4), 708
–713
(2010). https://doi.org/10.1364/AO.49.000708 APOPAI 00036935 Google Scholar
J. Li, H. Shen and R. Zhu,
“Method of alignment error control in freeform surface metrology with the tiltedwaveinterferometer,”
Opt. Eng., 55
(4), 044101
(2016). https://doi.org/10.1117/1.OE.55.4.044101 Google Scholar
R. D. Hersch and S. Chosson,
“Band moiré images,”
ACM Trans. Graphics, 23
(3), 239
–247
(2004). https://doi.org/10.1145/1015706 ATGRDF 07300301 Google Scholar
M. Ragulskis and A. Aleksa,
“Image hiding based on timeaveraging moiré,”
Opt. Commun., 282
(14), 2752
–2759
(2009). https://doi.org/10.1016/j.optcom.2009.04.007 OPCOB8 00304018 Google Scholar
V. Saveljev et al.,
“About a moiréless condition for nonsquare grids,”
J. Disp. Technol., 4
(3), 332
–339
(2008). https://doi.org/10.1109/JDT.2008.921900 IJDTAL 1551319X Google Scholar
C. Chiang,
“Moiré topography,”
Appl. Opt., 14
(1), 177
–179
(1975). https://doi.org/10.1364/AO.14.000177 APOPAI 00036935 Google Scholar
D. Post, B. Han and P. Ifju, High Sensitivity Moiré: Experimental Analysis for Mechanics and Materials, Springer, New York
(1994). Google Scholar
M. Ragulskis, A. Aleksa and L. Saunoriene,
“Improved algorithm for image encryption based on stochastic geometric moiré and its application,”
Opt. Commun., 273
(2), 370
–378
(2007). https://doi.org/10.1016/j.optcom.2007.01.055 OPCOB8 00304018 Google Scholar
I. Amidror, S. Chosson and R. D. Hersch,
“Moiré methods for the protection of documents and products,”
J. Phys. Conf. Ser., 77 012001
(2007). https://doi.org/10.1088/17426596/77/1/012001 JPCSDZ 17426588 Google Scholar
V. Ostromoukhov et al.,
“Anticounterfeiting features of artistic screening,”
Proc. SPIE, 2951 126
–133
(2010). https://doi.org/10.1117/12.262419 PSISDG 0277786X Google Scholar
I. Amidror and R. D. Hersch,
“The role of Fourier theory and of modulation in the prediction of visible moiré effects,”
J. Mod. Opt., 56
(9), 1103
–1118
(2009). https://doi.org/10.1080/09500340902994140 JMOPEW 09500340 Google Scholar
K. Patorski, K. Pokorski and M. Trusiak,
“Fourier domain interpretation of real and pseudomoiré phenomena,”
Opt. Express, 19
(27), 26065
–26078
(2011). https://doi.org/10.1364/OE.19.026065 OPEXFF 10944087 Google Scholar
I. Amidror, Mastering the Discrete Fourier Transform in One, Two or Several Dimensions, Springer, London, United Kingdom
(2013). Google Scholar
V. Saveljev and S.K. Kim,
“Theoretical estimation of moiré effect using spectral trajectories,”
Opt. Express, 21
(2), 1693
–1712
(2013). https://doi.org/10.1364/OE.21.001693 OPEXFF 10944087 Google Scholar
V. Saveljev and S.K. Kim,
“Experimental amplitude, period and orientation of the moiré patterns in barrier 3D displays,”
J. Inf. Disp.,
(2018). https://doi.org/10.1080/15980316.2018.1452301 Google Scholar
V. Saveljev,
“Orientations and branches of moiré waves in threedimensional displays,”
J. Korean Phys. Soc., 57
(6), 1392
–1396
(2010). https://doi.org/10.3938/jkps.57.1392 KPSJAS 03744884 Google Scholar
V. Saveljev and S.K. Kim,
“Simulation and measurement of moiré patterns at finite distance,”
Opt. Express, 20
(3), 2163
–2177
(2012). https://doi.org/10.1364/OE.20.002163 OPEXFF 10944087 Google Scholar
V. Saveljev,
“Moiré effect in cylindrical objects,”
J. Korean Phys. Soc., 68
(9), 1075
–1082
(2016). https://doi.org/10.3938/jkps.68.1075 KPSJAS 03744884 Google Scholar
V. Saveljev and J. Kim,
“Moiré effect in singlelayer cylindrical surface,”
in Proc. Optical Society of Korea Winter Annual Meeting,
5
(2017). Google Scholar
V. Saveljev and J. Kim,
“Towards physical simulation of the moiré effect in nanoparticles,”
in Proc. Spring Meeting of Korean Physical Society,
(2017). Google Scholar
G. Strang, Computational Science and Engineering, WellesleyCambridge University, Wellesley, Massachusetts
(2007). Google Scholar
C. Huang et al.,
“Developing the spectral trajectories of major land cover change processes,”
Proc. SPIE, 3502 155
–162
(1998). https://doi.org/10.1117/12.317807 Google Scholar
P. S. CostaPereira and P. Maillard,
“Estimating the age of cerrado regeneration using Landsat TM data,”
Can. J. Remote Sens., 36
(S2), S243
–S256
(2010). https://doi.org/10.5589/m10059 CJRSDP 07038992 Google Scholar
V. Saveljev and S.K. Kim,
“Estimation of moiré patterns using spectral trajectories in the complex plane,”
Comput. Technol. Appl., 3
(5), 353
–360
(2012). https://doi.org/10.17265/19347332/2012.05.004 Google Scholar
V. Saveljev and S.K. Kim,
“Probability of the moiré effect in barrier and lenticular autostereoscopic 3D displays,”
Opt. Express, 23
(20), 25597
–25607
(2015). https://doi.org/10.1364/OE.23.025597 OPEXFF 10944087 Google Scholar
J. Li et al.,
“Formation mechanism and a universal period formula for the CCD moiré,”
Opt. Express, 22
(17), 20914
–20923
(2014). https://doi.org/10.1364/OE.22.020914 OPEXFF 10944087 Google Scholar
K.B. Kim et al.,
“Investigation of patterninduced brightness nonuniformity in activematrix organic lightemitting diode displays,”
Jpn. J. Appl. Phys., 47
(1), 193
–196
(2008). https://doi.org/10.1143/JJAP.47.193 Google Scholar
S.J. Byun et al.,
“An efficient simulation and analysis method of moiré patterns in display systems,”
Opt. Express, 22
(3), 3128
–3136
(2014). https://doi.org/10.1364/OE.22.003128 OPEXFF 10944087 Google Scholar
Y. Zhou et al.,
“Quantitative measurement and control of optical moiré pattern in an autostereoscopic liquid crystal display system,”
Appl. Opt., 54
(6), 1521
–1527
(2015). https://doi.org/10.1364/AO.54.001521 APOPAI 00036935 Google Scholar
K. Nagasaki and Y. Bao,
“A 3D display with variable depth moiré pattern,”
in Proc. SICE Annual Conf.,
1936
–1941
(2008). Google Scholar
L. Kong et al.,
“Parameter design of a parallax barrier based on the color moiré patterns in autostereoscopic display,”
Appl. Opt., 50
(34), H153
–H158
(2011). https://doi.org/10.1364/AO.50.00H153 APOPAI 00036935 Google Scholar
B.R. Lee et al.,
“Color moiré simulations in contacttype 3D displays,”
Opt. Express, 23
(11), 14114
–14125
(2015). https://doi.org/10.1364/OE.23.014114 OPEXFF 10944087 Google Scholar
C.W. Jeon et al.,
“Wall pattern design to reduce color moiré in autostereoscopic 3D display,”
Mol. Cryst. Liq. Cryst., 613
(1), 63
–68
(2015). https://doi.org/10.1080/15421406.2015.1032046 Google Scholar
H. Lee et al.,
“Fringe periods of color moirés in contacttype 3D displays,”
Opt. Express, 24
(13), 14183
–14195
(2016). https://doi.org/10.1364/OE.24.014183 OPEXFF 10944087 Google Scholar
F. P. Chiang and B. Ranganayakamma,
“Deflection measurements using moiré gap effect,”
Exp. Mech., 11
(4), 296
–302
(1971). https://doi.org/10.1007/BF02320582 EXMCAZ 00144851 Google Scholar
F. P. Chiang and R. M. Juan,
“A method to shift moiré fringes using gap effect,”
Exp. Mech., 13
(5), 209
–211
(1973). https://doi.org/10.1007/BF02322655 EXMCAZ 00144851 Google Scholar
D. Yogev, S. Efrima and O. Kafri,
“Study of the thickness of liquid layers by moiré deflectometry,”
Opt. Lett., 13
(10), 934
–936
(1988). https://doi.org/10.1364/OL.13.000934 OPLEDP 01469592 Google Scholar
J. Kim et al.,
“Behaviors of moiré fringes induced by plate thickness,”
J. Opt., 17
(3), 1
–7
(2015). https://doi.org/10.1088/20408978/17/3/035801 Google Scholar
O. Mikami,
“New imagerotation using moiré lenses,”
Jpn. J. Appl. Phys., 14
(7), 1065
–1066
(1975). https://doi.org/10.1143/JJAP.14.1065 Google Scholar
G. Oster and Y. Nishijima,
“Moiré patterns,”
Sci. Am., 208
(5), 54
–63
(1963). https://doi.org/10.1038/scientificamerican056354 SCAMAC 00368733 Google Scholar
G. Oster,
“Optical art,”
Appl. Opt., 4
(11), 1359
–1369
(1965). https://doi.org/10.1364/AO.4.001359 APOPAI 00036935 Google Scholar
N. J. Wade,
“Op art and visual perception,”
Perception, 7
(1), 21
–46
(1978). https://doi.org/10.1068/p070021 PCTNBA 03010066 Google Scholar
A. Minini,
“Andrea Minini, illustrator & graphic designer,”
(2018). andreaminini.myportfolio.com Google Scholar
P. Decrauzat,
“Philippe Decrauzat, selected works,”
(2018) www.prazdelavallade.com/philippedecrauzat March ). 2018). Google Scholar
C. CruzDiez,
“Carlos CruzDiez, Venezuelan, b. 1923,”
(2018) www.artsy.net/artist/carloscruzdiez/works March ). 2018). Google Scholar
BiographyVladimir Saveljev is a research professor at the Myongji University, Yongin, South Korea. He received his PhD from Institute of Automation and Electrometry (IA&E), Siberian Branch of the Russian Academy of Sciences in 2014. After IA&E, he was with Korea Institute of Science and Technology and Hanyang University. He is the author of 150 papers and three book chapters. His current research interests include autostereoscopic 3D displays, image quality, moiré effect, wavelet transform, and nanoparticles. SungKyu Kim received his PhD from Quantum Optics Group of Physics, Korea University, in 2000. He spent two years as an invited research scientist at 3D TV Group of Telecommunications Advancement Organization, Japan. In 2001, he joined at Korea Institute of Science and Technology, where he is currently a principal research scientist. His research interests include optical design of autostereoscopic 3D display systems, multiview image processing, digital holography, holographic optical elements, and multifocus 3D display. Jaisoon Kim received his undergraduate degree in physics education from Seoul National University in 1980, his master’s degree in physics education from Seoul National University in 1987, and his doctorate in physics from Korea University in 1999. He has been a research scientist at OSC for 2000–2002 and a professor at Seoul National University for 2002–2009. Currently, he is a professor at the Myongji University, Korea. 