## 1.

## Introduction

The 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 three-dimensional (3-D) displays was first stated^{14} in the early 2000s. This is especially important for autostereoscopic 3-D 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 x-rays,^{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 3-D 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 3-D). 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 3-D displays based entirely on the moiré effect as the main physical principle, or to improve security.^{51}^{–}^{53} An image in a 3-D 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 self-check 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 3-D 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 two-dimensional (2-D) 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 Equation

The 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{1-2\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

1. Calculate the period and orientation of the moiré patterns in the parallel line gratings with their periods 1 and 1.1 mm.

2. The same for the identical gratings (period 2 mm) installed at the angles of 0 deg, 5 deg, and 10 deg.

3. Calculate the period of the moiré patterns in the square grids with the same periods and angles as in the problems 1 and 2.

4. The same for the line grating 1 mm and 30 deg-rhomboidal grid 1.1 mm at the angle of 25 deg between them.

5. Describe the shape of the moiré patterns in the circular (radial/concentric) grating + line grating.

## 3.

## Sinusoidal Coplanar Gratings

The 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 Gratings

An interaction between the gratings can be modeled mathematically by the multiplication of transparency functions of the gratings. For instance, a one-dimensional (1-D) 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 Grating

Superimposed 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 (kx-ky)+\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(x-y)+\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{1-2\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 short-period 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

1. Point out the terms of Eq. (15) responsible for the strongest functions given in Eqs. (19) and (20).

2. Write down the equation of strongest waves for the 60-deg rhomboidal grating combined with a line grating.

3. Find the period of residual moiré patterns at the optimal angle (you may use either of two branches).

4. Find the period of the moiré patterns in square grating 1 mm + line grating 1.1 mm at the angles of 0 deg, 5 deg, and 10 deg.

5. Determine the orientation of the moiré pattern in line gratings 1 mm + 1.2 mm installed at 10 deg.

## 4.

## Sinusoidal Noncoplanar Gratings

Practically, the layers may not lie in the same plane. For example, a typical structure of the autostereoscopic 3-D 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

1. Find the period of the moiré patterns in the identical parallel gratings (period 2 mm) installed with the gap of 1, 5, and 10 mm. The observer distance is 1 m.

2. Rewrite the equation for the moiré factor using one variable only. What could be its physical meaning?

3. Determine the gap based on the moiré period for:

4. Find the angle of the moiré image for the identical gratings in parallel planes (but not necessarily parallel orientation of gratings).

5. Obtain the homogeneous transformation matrices for the two cases shown in Fig. 12.

6. What is a condition to see the moiré patterns of the infinite period from the finite distance in the parallel players installed with a gap?

7. Is it physically possible for the moiré factor to be equal to one?

## 5.

## Cylindrical Moiré (Orthogonal Projection)

Consider a regular 3-D 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 re-expressed 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

1. The period of the moiré patterns for the grating period 1 mm and the diameter 10 and 20 mm; the observer distance is 0.5 m.

2. Draw a sketch of the patterns in wrapped skew gratings (neither vertical nor horizontal layout, but an intermediate inclination angle between 0 deg and 90 deg) based on Eqs. (35) and (36).

3. Compare the periods of the moiré patterns on the axis of the cylinder for the grating with the lines parallel to the axis of the cylinder and for the inclined grating with the correspondingly inclined parallel gratings (their gap = diameter of the cylinder).

## 6.

## Moiré Wave Vector as a Vector Sum

Now, 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 re-expressed 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

1. Rewrite the expression for the wave vector using the ratio of wave numbers.

2. Determine the minimum/maximum period from Eq. (49) depending on the angle.

3. Draw the moiré wave vector for three line gratings of arbitrary periods and angles.

4. Find the orientation of the moiré patterns at min/max period from exercise 2.

## 7.

## Basics of Two-Dimensional Fourier Transform

## 7.1.

### Spectra

A 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.

### Two-Dimensional Spectra

Figure 22 of the previous section shows the 1-D functions. A 1-D function of one variable (graphically, the “height” $y$ as a function of the abscissa $x$) can be generalized to a 2-D function of two variables (the “height” $z$ as a function of two independent coordinates $x$ and $y$).

The 2-D 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 1-D spectrum) and repeated (copied) in the orthogonal direction. A 2-D 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 2-D grid can be represented as a product of two 1-D 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

1. Describe the spectrum of the sinusoidal square grid (which is a superposition of two line gratings).

2. Describe the spectrum of the nonsinusoidal square grid.

3. Draw the map of peaks for three superposed sinusoidal gratings at the angle near 60 deg.

4. What is the 2-D Fourier transform of two orthogonal gratings with the periodic triangular transparency function?

5. How to find the phase of the wave from the Fourier coefficients?

6. Distance between the spectral peaks of the function $\mathrm{sin}(3x+4y)$.

## 8.

## Spectral Trajectories

Generally 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 1-D structure, the square grid is a 2-D 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

1. Draw a map of peaks in the case of three overlapped sinusoidal gratings installed at the angle near 60 deg.

2. The same for 30 deg.

3. Draw a sketch of the spectral trajectories of two line gratings for the running parameter $\rho $ and the angles $\alpha =15\text{\hspace{0.17em}\hspace{0.17em}}\mathrm{deg}$, and $\alpha =20\text{\hspace{0.17em}\hspace{0.17em}}\mathrm{deg}$.

4. Write equations of trajectories leaving the origin in Fig. 25(a).

5. Find the distance to the origin for the trajectories approaching the origin in Fig. 25(b).

## 9.

## Visual Effects on the Move and on the Rotation

In the case of a laterally moved 1-D 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 so-called 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 long-distance 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

1. Make a drawing of the square grid combined with the line grating near the angle of 45 deg. Describe the visual picture of the moiré patterns for a specific ratio of periods.

2. Describe the period of the moiré patterns in the gratings when the observer approaches the screen; consider the coplanar and noncoplanar gratings.

3. Describe the moiré picture for noncoplanar square grids.

4. Can two observers see identical moiré images, if they stand shoulder to shoulder? Behind each other?

## 10.

## Discussion and Conclusion

The 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 2-D 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 3-D, 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. Cruz-Diez, 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.

## Acknowledgments

This work was partially supported by the Cross-Ministry Giga KOREA Project through the Korean Government, Ministry of Science and ICT (MSIT); Grant No. GK17D0200, Development of Super Multi-View Display Providing Real-Time Interaction.

## References

**,” J. Opt. Soc. Am., 64 (10), 1287 –1294 (1974). https://doi.org/10.1364/JOSA.64.001287 JOSAAH 0030-3941 Google Scholar**

*Moiré: formation and interpretation***,” Opt. Commun., 11 (4), 378 –381 (1974). https://doi.org/10.1016/0030-4018(74)90240-5 OPCOB8 0030-4018 Google Scholar**

*Theoretical interpretation of the moiré pattern***,” Appl. Opt., 15 (9), 2223 –2227 (1976). https://doi.org/10.1364/AO.15.002223 APOPAI 0003-6935 Google Scholar**

*Geometric parameters of moiré fringes***,” Exp. Mech., 22 (11), 418 –433 (1982). https://doi.org/10.1007/BF02326823 EXMCAZ 0014-4851 Google Scholar**

*The moiré method—a review***,” Opt. Commun., 212 (4–6), 211 –216 (2002). https://doi.org/10.1016/S0030-4018(02)02017-5 OPCOB8 0030-4018 Google Scholar**

*Moiré-like patterns as a spatial beat frequency in photorefractive sinusoidal phase gratings superposition***,” Pure Appl. Opt., 3 133 –142 (1994). https://doi.org/10.1088/0963-9659/3/2/006 PAOAE3 0963-9659 Google Scholar**

*Moiré magnifier***,” Proc. SPIE, 4016 166 –170 (1999). https://doi.org/10.1117/12.373616 PSISDG 0277-786X Google Scholar**

*Moiré in a scanned image***,” Appl. Opt., 48 (11), 2178 –2187 (2009). https://doi.org/10.1364/AO.48.002178 APOPAI 0003-6935 Google Scholar**

*Color moiré pattern simulation and analysis in three-dimensional integral imaging for finding the moiré-reduced tilted angle of a lens array***,” in Proc. 10th Int. Display Research Workshops (IDW), 1401 –1404 (2003). Google Scholar**

*Moiré pattern reduction in full-parallax autostereoscopic imaging systems using two crossed lenticular plates as a viewing zone forming optics***,” J. Disp. Technol., 1 (1), 125 –140 (2005). https://doi.org/10.1109/JDT.2005.853354 IJDTAL 1551-319X Google Scholar**

*Three-dimensional imaging systems based on multiview images***,” Proc. IEEE, 94 (3), 591 –607 (2006). https://doi.org/10.1109/JPROC.2006.870696 IEEPAD 0018-9219 Google Scholar**

*Three-dimensional image sensing, visualization, and processing using integral imaging***,” Handbook of Optoelectronics, 1067 –1100 Taylor & Francis, London, United Kingdom (2006). Google Scholar**

*Three-dimensional display systems***,” J. Soc. Inf. Disp., 15 (11), 883 –888 (2007). https://doi.org/10.1889/1.2812988 JSIDE8 0734-1768 Google Scholar**

*Moiré interference in multilayered displays***,” Proc. IEEE, 97 (6), 1067 –1077 (2009). https://doi.org/10.1109/JPROC.2009.2016816 IEEPAD 0018-9219 Google Scholar**

*Progress in 3-D multiperspective display by integral imaging***,” J. Disp. Technol., 6 (10), 394 –403 (2010). https://doi.org/10.1109/JDT.2010.2045636 IJDTAL 1551-319X Google Scholar**

*Recent developments in 3-D imaging technologies***,” Appl. Opt., 50 (34), H87 –H115 (2011). https://doi.org/10.1364/AO.50.000H87 APOPAI 0003-6935 Google Scholar**

*Three-dimensional display technologies of recent interest: principles, status, and issues***,” Opt. Express, 21 (22), 26068 –26079 (2013). https://doi.org/10.1364/OE.21.026068 OPEXFF 1094-4087 Google Scholar**

*Analysis of moiré minimization in autostereoscopic parallax displays***,” Proc. R. Soc. London Ser. A, 246 (1246), 345 –368 (1958). https://doi.org/10.1098/rspa.1958.0144 Google Scholar**

*Moiré patterns on electron micrographs, and their application to the study of dislocations in metals***,” Microwave Opt. Technol. Lett., 49 (11), 2761 –2763 (2007). https://doi.org/10.1002/(ISSN)1098-2760 Google Scholar**

*Measurement on infrared radiation burning temperature field with moiré deviation tomography***,” Rev. Phys. Appl., 11 (1), 83 –87 (1976). https://doi.org/10.1051/rphysap:0197600110108300 RPHAAN 0035-1687 Google Scholar**

*X-ray interferometry and lattice parameter investigation***,” Phys. Rev. B, 81 125427 (2010). https://doi.org/10.1103/PhysRevB.81.125427 Google Scholar**

*Structural analysis of multilayer graphene via atomic moiré interferometry***,” Proc. Phys. Soc., 73 (1), 142 –144 (1959). https://doi.org/10.1088/0370-1328/73/1/129 Google Scholar**

*A simple method of calculating moiré patterns***,” Z. Angew. Math. Phys., 19 (2), 326 –333 (1968). https://doi.org/10.1007/BF01601476 Google Scholar**

*Gap effect on moiré patterns***,” J. Opt. Soc. Am., 66 (2), 87 –94 (1976). https://doi.org/10.1364/JOSA.66.000087 JOSAAH 0030-3941 Google Scholar**

*Characteristics of superposed patterns in optics***,” Jpn. J. Appl. Phys., 15 (3), 443 –456 (1976). https://doi.org/10.1143/JJAP.15.443 Google Scholar**

*Moiré profile prediction by using Fourier series formalism***,” Opt. Laser Eng., 3 (1), 15 –27 (1982). https://doi.org/10.1016/0143-8166(82)90012-4 Google Scholar**

*Moiré fringes***,” Exp. Mech., 5 (5), 154 –160 (1965). https://doi.org/10.1007/BF02324049 EXMCAZ 0014-4851 Google Scholar**

*Basic optical law in the interpretation of moiré patterns applied to the analysis of strains—Part 1***,” Exp. Mech., 7 (7), 289 –296 (1967). https://doi.org/10.1007/BF02327134 EXMCAZ 0014-4851 Google Scholar**

*Moiré topography of curved surfaces***,” J. Phys. E Sci. Instrum., 5 (3), 217 –219 (1972). https://doi.org/10.1088/0022-3735/5/3/010 Google Scholar**

*Vibration measurements by the moiré method***,” Optical Shop Testing, 653 –685 John Wiley & Sons, New York (1992). Google Scholar**

*Moiré and fringe projection techniques***,” Meas. Sci. Rev., 1 (3), 29 –32 (2001). Google Scholar**

*Inspection of surface by the moiré method***,” J. Phys. Conf. Ser., 274 (1), 012036 (2011). https://doi.org/10.1088/1742-6596/274/1/012036 Google Scholar**

*Dynamic moiré patterns for profilometry applications***,” in Proc. 33rd Int. Display Research Conf. (EuroDisplay), 54 –56 (2013). Google Scholar**

*A 3D moiré display***,” J. Soc. Inf. Disp., 22 (9), 482 –486 (2014). https://doi.org/10.1002/jsid.v22.9 JSIDE8 0734-1768 Google Scholar**

*Three-dimensional moiré display***,” Opt. Eng., 42 (9), 2646 –2652 (2003). https://doi.org/10.1117/1.1592803 Google Scholar**

*Sensitivity improvement in phase-shifted moiré interferometry using 1-D continuous wavelet transform image processing***,” Opt. Eng., 40 (9), 1978 –1983 (2001). https://doi.org/10.1117/1.1394738 Google Scholar**

*Influence of object roughness on specimen gratings for moiré interferometry***,” Opt. Eng., 31 (3), 434 –439 (1992). https://doi.org/10.1117/12.56108 Google Scholar**

*Polarization approach to high-sensitivity moiré interferometry***,” Opt. Eng., 24 (6), 944 –960 (1985). https://doi.org/10.1117/12.7973607 Google Scholar**

*Moiré deflectometry: a ray deflection approach to optical testing***,” Appl. Opt., 49 (4), 708 –713 (2010). https://doi.org/10.1364/AO.49.000708 APOPAI 0003-6935 Google Scholar**

*Extended dual-grating alignment method for optical projection lithography***,” Opt. Eng., 55 (4), 044101 (2016). https://doi.org/10.1117/1.OE.55.4.044101 Google Scholar**

*Method of alignment error control in free-form surface metrology with the tilted-wave-interferometer***,” ACM Trans. Graphics, 23 (3), 239 –247 (2004). https://doi.org/10.1145/1015706 ATGRDF 0730-0301 Google Scholar**

*Band moiré images***,” Opt. Commun., 282 (14), 2752 –2759 (2009). https://doi.org/10.1016/j.optcom.2009.04.007 OPCOB8 0030-4018 Google Scholar**

*Image hiding based on timeaveraging moiré***,” J. Disp. Technol., 4 (3), 332 –339 (2008). https://doi.org/10.1109/JDT.2008.921900 IJDTAL 1551-319X Google Scholar**

*About a moiré-less condition for non-square grids***,” Appl. Opt., 14 (1), 177 –179 (1975). https://doi.org/10.1364/AO.14.000177 APOPAI 0003-6935 Google Scholar**

*Moiré topography***,” Opt. Commun., 273 (2), 370 –378 (2007). https://doi.org/10.1016/j.optcom.2007.01.055 OPCOB8 0030-4018 Google Scholar**

*Improved algorithm for image encryption based on stochastic geometric moiré and its application***,” J. Phys. Conf. Ser., 77 012001 (2007). https://doi.org/10.1088/1742-6596/77/1/012001 JPCSDZ 1742-6588 Google Scholar**

*Moiré methods for the protection of documents and products***,” Proc. SPIE, 2951 126 –133 (2010). https://doi.org/10.1117/12.262419 PSISDG 0277-786X Google Scholar**

*Anti-counterfeiting features of artistic screening***,” J. Mod. Opt., 56 (9), 1103 –1118 (2009). https://doi.org/10.1080/09500340902994140 JMOPEW 0950-0340 Google Scholar**

*The role of Fourier theory and of modulation in the prediction of visible moiré effects***,” Opt. Express, 19 (27), 26065 –26078 (2011). https://doi.org/10.1364/OE.19.026065 OPEXFF 1094-4087 Google Scholar**

*Fourier domain interpretation of real and pseudo-moiré phenomena***,” Opt. Express, 21 (2), 1693 –1712 (2013). https://doi.org/10.1364/OE.21.001693 OPEXFF 1094-4087 Google Scholar**

*Theoretical estimation of moiré effect using spectral trajectories***,” J. Inf. Disp., (2018). https://doi.org/10.1080/15980316.2018.1452301 Google Scholar**

*Experimental amplitude, period and orientation of the moiré patterns in barrier 3D displays***,” J. Korean Phys. Soc., 57 (6), 1392 –1396 (2010). https://doi.org/10.3938/jkps.57.1392 KPSJAS 0374-4884 Google Scholar**

*Orientations and branches of moiré waves in three-dimensional displays***,” Opt. Express, 20 (3), 2163 –2177 (2012). https://doi.org/10.1364/OE.20.002163 OPEXFF 1094-4087 Google Scholar**

*Simulation and measurement of moiré patterns at finite distance***,” J. Korean Phys. Soc., 68 (9), 1075 –1082 (2016). https://doi.org/10.3938/jkps.68.1075 KPSJAS 0374-4884 Google Scholar**

*Moiré effect in cylindrical objects***,” in Proc. Optical Society of Korea Winter Annual Meeting, 5 (2017). Google Scholar**

*Moiré effect in single-layer cylindrical surface***,” in Proc. Spring Meeting of Korean Physical Society, (2017). Google Scholar**

*Towards physical simulation of the moiré effect in nanoparticles***,” Proc. SPIE, 3502 155 –162 (1998). https://doi.org/10.1117/12.317807 Google Scholar**

*Developing the spectral trajectories of major land cover change processes***,” Can. J. Remote Sens., 36 (S2), S243 –S256 (2010). https://doi.org/10.5589/m10-059 CJRSDP 0703-8992 Google Scholar**

*Estimating the age of cerrado regeneration using Landsat TM data***,” Comput. Technol. Appl., 3 (5), 353 –360 (2012). https://doi.org/10.17265/1934-7332/2012.05.004 Google Scholar**

*Estimation of moiré patterns using spectral trajectories in the complex plane***,” Opt. Express, 23 (20), 25597 –25607 (2015). https://doi.org/10.1364/OE.23.025597 OPEXFF 1094-4087 Google Scholar**

*Probability of the moiré effect in barrier and lenticular autostereoscopic 3D displays***,” Opt. Express, 22 (17), 20914 –20923 (2014). https://doi.org/10.1364/OE.22.020914 OPEXFF 1094-4087 Google Scholar**

*Formation mechanism and a universal period formula for the CCD moiré***,” Jpn. J. Appl. Phys., 47 (1), 193 –196 (2008). https://doi.org/10.1143/JJAP.47.193 Google Scholar**

*Investigation of pattern-induced brightness non-uniformity in active-matrix organic light-emitting diode displays***,” Opt. Express, 22 (3), 3128 –3136 (2014). https://doi.org/10.1364/OE.22.003128 OPEXFF 1094-4087 Google Scholar**

*An efficient simulation and analysis method of moiré patterns in display systems***,” Appl. Opt., 54 (6), 1521 –1527 (2015). https://doi.org/10.1364/AO.54.001521 APOPAI 0003-6935 Google Scholar**

*Quantitative measurement and control of optical moiré pattern in an autostereoscopic liquid crystal display system***,” in Proc. SICE Annual Conf., 1936 –1941 (2008). Google Scholar**

*A 3D display with variable depth moiré pattern***,” Appl. Opt., 50 (34), H153 –H158 (2011). https://doi.org/10.1364/AO.50.00H153 APOPAI 0003-6935 Google Scholar**

*Parameter design of a parallax barrier based on the color moiré patterns in autostereoscopic display***,” Opt. Express, 23 (11), 14114 –14125 (2015). https://doi.org/10.1364/OE.23.014114 OPEXFF 1094-4087 Google Scholar**

*Color moiré simulations in contact-type 3-D displays***,” Mol. Cryst. Liq. Cryst., 613 (1), 63 –68 (2015). https://doi.org/10.1080/15421406.2015.1032046 Google Scholar**

*Wall pattern design to reduce color moiré in autostereoscopic 3D display***,” Opt. Express, 24 (13), 14183 –14195 (2016). https://doi.org/10.1364/OE.24.014183 OPEXFF 1094-4087 Google Scholar**

*Fringe periods of color moirés in contact-type 3-D displays***,” Exp. Mech., 11 (4), 296 –302 (1971). https://doi.org/10.1007/BF02320582 EXMCAZ 0014-4851 Google Scholar**

*Deflection measurements using moiré gap effect***,” Exp. Mech., 13 (5), 209 –211 (1973). https://doi.org/10.1007/BF02322655 EXMCAZ 0014-4851 Google Scholar**

*A method to shift moiré fringes using gap effect***,” Opt. Lett., 13 (10), 934 –936 (1988). https://doi.org/10.1364/OL.13.000934 OPLEDP 0146-9592 Google Scholar**

*Study of the thickness of liquid layers by moiré deflectometry***,” J. Opt., 17 (3), 1 –7 (2015). https://doi.org/10.1088/2040-8978/17/3/035801 Google Scholar**

*Behaviors of moiré fringes induced by plate thickness***,” Jpn. J. Appl. Phys., 14 (7), 1065 –1066 (1975). https://doi.org/10.1143/JJAP.14.1065 Google Scholar**

*New image-rotation using moiré lenses***,” Sci. Am., 208 (5), 54 –63 (1963). https://doi.org/10.1038/scientificamerican0563-54 SCAMAC 0036-8733 Google Scholar**

*Moiré patterns***,” Appl. Opt., 4 (11), 1359 –1369 (1965). https://doi.org/10.1364/AO.4.001359 APOPAI 0003-6935 Google Scholar**

*Optical art***,” Perception, 7 (1), 21 –46 (1978). https://doi.org/10.1068/p070021 PCTNBA 0301-0066 Google Scholar**

*Op art and visual perception***,” (2018). andreaminini.myportfolio.com Google Scholar**

*Andrea Minini, illustrator & graphic designer***,” (2018) www.praz-delavallade.com/philippe-decrauzat March ). 2018). Google Scholar**

*Philippe Decrauzat, selected works***,” (2018) www.artsy.net/artist/carlos-cruz-diez/works March ). 2018). Google Scholar**

*Carlos Cruz-Diez, Venezuelan, b. 1923*## Biography

**Vladimir 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 3-D displays, image quality, moiré effect, wavelet transform, and nanoparticles.

**Sung-Kyu Kim** received his PhD from Quantum Optics Group of Physics, Korea University, in 2000. He spent two years as an invited research scientist at 3-D 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 3-D display systems, multiview image processing, digital holography, holographic optical elements, and multifocus 3-D 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.