## 1.

## Introduction

Modulation transfer function (MTF) has been widely used to evaluate the image quality of imaging systems because it works in a direct way to reveal the performance of imaging systems.^{1}2.^{–}^{3} MTF describes the spatial frequency transfer characteristics of imaging systems from object space to image space by employing the contrast value of a series of sine waves with varying spatial frequencies.

There are a lot of mature ways to measure the MTF of optical imaging systems. According to their different templates, MTF test methods can be categorized into different types, such as the knife edge method,^{4} the slit method,^{5} the bar target method,^{6}7.^{–}^{8} etc. Generally, these methods’ templates have a common feature in that their image brightness changes in one direction and remains constant in the orthogonal direction, e.g., the U.S. Air Force resolution target. Naturally, the MTF measured with these templates will contain the modulation information in only one direction, such as the horizontal MTF and vertical MTF. In this paper, these MTFs are called one-dimensional (1-D) MTFs and their corresponding templates are called 1-D templates.

In 2010, Marom et al.^{9} proposed a new kind of two-dimensional (2-D) template which looks like a checkerboard measuring the “true characterization” of some systems, for example, the barcode scanning system. Marom et al. explained that a 2-D template is different from a 1-D template because its image brightness changes in two orthogonal directions, and the result of a 1-D contrast transfer function (CTF) for a 2-D template will be affected by the brightness changes in another direction, which will not occur with the 1-D template. Haim et al. proved theoretically and experimentally that the CTF measured by the 2-D template is lower than that obtained with 1-D templates at every spatial frequency and at its cutoff frequency while the former template is $1/\sqrt{2}$ times the latter one.^{10}^{,}^{11} To maintain consistency, in this paper we adopt Marom’s 2-D template definition in which the brightness of 2-D templates changes in two orthogonal directions. The MTF that simultaneously contains the modulation information in two orthogonal directions is called the 2-D MTF.

Recently, the random target method, a new MTF test method, has been applied for image quality evaluation. The method was first proposed by Daniels et al.^{12} in which a random image is used as the template, and then it is imaged on a charge coupled device (CCD) by the optical system under test (OSUT). By using the power spectral density (PSD) values of the image captured by a CCD and those of the original image, one can obtain the optical system’s MTF with Fourier spectral analysis. The test method is convenient and fast,^{13} because the random image has the characteristic of shift invariance,^{14} and one can easily achieve automatic measurement. The previous works focus mostly on the 1-D MTF test.^{9}10.11.12.^{–}^{13}^{,}^{15}16.17.^{–}^{18} However, the random image can also be used as a 2-D target since its image brightness changes in two orthogonal directions and is available for measuring 2-D MTF. In 2013, Evtikhiev et al.^{19} proposed a new random target to reduce the noise impact for the 2-D MTF measurement.

In this paper, we investigate the measurement for 2-D MTF with the random target method as well as 1-D MTF. A new approach is applied to reduce noise by the use of a significantly large number of random images. The results of the measurement for 1-D and 2-D MTF show that 2-D MTF can be more comprehensive than 1-D MTF when describing the true modulation features of optical imaging systems.

Section 2 analyzes the 1-D and 2-D MTF theoretically and elaborates on the mathematical principles for 1-D and 2-D MTF measurement by using the random target method. Section 3 presents the procedures of the experiment and the results. Section 4 contains the conclusions.

## 2.

## Methodology

## 2.1.

### MTF Review and Comparison of 1-D and 2-D MTF

In incoherent illumination systems, the optical transfer function (OTF) is the normalized autocorrelation of the system exit pupil, which can be expressed as

## (1)

$$\mathrm{OTF}(\xi ,\eta )=\frac{{\iint}_{-\infty}^{\infty}P(x,y){P}^{*}(x+\lambda {z}_{i}\xi ,y+\lambda {z}_{i}\eta )\mathrm{d}x\text{\hspace{0.17em}}\mathrm{d}y}{{\iint}_{-\infty}^{\infty}P(x,y)\mathrm{d}x\text{\hspace{0.17em}}\mathrm{d}y},$$MTF describes the spatial frequency transfer characteristics of imaging systems defined as the ratio of the output to the input sinusoidal targets’ contrast value. In this way, one can easily obtain the MTF through detecting the contrast value of input and output sinusoidal targets, but this method is rarely used because it is hard to exactly and efficiently fabricate sinusoidal targets. In our work, to reveal the internal relationship of the 1-D and 2-D MTF of the optical system, the sinusoidal targets are applied by use of Fourier transform.

At first, we analyze the 1-D MTF with a typical 1-D sinusoidal template which is shown in Fig. 1(a) and its expression is

where $f$ is the spatial frequency, and $f=1/(2a)$. If we treat the OSUT as a linear time-invariant system, then the input sinusoidal signal will be modulated by the system’s MTF when it transmits through the system, and the output signal will be## (3)

$${I}_{\text{output}}^{1\text{-}\mathrm{D}}=\frac{1}{2}[1+\mathrm{sin}(2\pi xf)\xb7{\mathrm{MTF}}_{1\text{-}\mathrm{D}}(f)]=\frac{1}{2}[1+\mathrm{sin}(\frac{\pi x}{a})\xb7{\mathrm{MTF}}_{1\text{-}\mathrm{D}}(\frac{1}{2a}\left)\right],$$The 2-D case is more complex than that of the 1-D. A typical 2-D sinusoidal target is shown in Fig. 1(b), which can be expressed as

## (4)

$${I}_{\text{input}}^{2\text{-}\mathrm{D}}=\frac{1}{2}[1+\mathrm{sin}(2\pi x{f}_{x})\xb7\mathrm{sin}(2\pi y{f}_{y})],$$## (5)

$$H(\xi ,\eta )=\frac{1}{2}\{1+\frac{1}{4}[\delta (\xi -{f}_{x})-\delta (\xi +{f}_{x})]\phantom{\rule{0ex}{0ex}}\xb7[\delta (\eta +{f}_{y})-\delta (\eta -{f}_{y})]\}\phantom{\rule{0ex}{0ex}}=\frac{1}{2}\{1+\frac{1}{4}[\delta (\xi +{f}_{x},\eta -{f}_{y})+\delta (\xi -{f}_{x},\eta +{f}_{y})\phantom{\rule{0ex}{0ex}}-\delta (\xi +{f}_{x},\eta +{f}_{y})-\delta (\xi -{f}_{x},\eta -{f}_{y})]\}.$$Equation (5) shows that 2-D target contains different frequency components in the 2-D frequency domain, and indicates the modulation features of the imaging system in the 2-D target case when adopting the 2-D MTF.

In Fig. 2, one can easily note that four different frequency components uniformly distribute on a circle whose radius is $f={({f}_{x}^{2}+{f}_{y}^{2})}^{1/2}$. In terms of an ideal axis symmetrical system (aberration-free), its 2-D MTF values will be equal to all the frequency points on the circle which can be calculated with equation

## (6)

$${\mathrm{MTF}}_{2\text{-}\mathrm{D}}({f}_{x},{f}_{y})={\mathrm{MTF}}_{1\text{-}\mathrm{D}}\left(\sqrt{{f}_{x}^{2}+{f}_{y}^{2}}\right),$$## (7)

$${I}_{\text{output}}^{2\text{-}\mathrm{D}}=\frac{1}{2}[1+\mathrm{sin}(2\pi x{f}_{x})\xb7\mathrm{sin}(2\pi y{f}_{y})\xb7{\mathrm{MTF}}_{2\text{-}\mathrm{D}}({f}_{x},{f}_{y})]=\frac{1}{2}[1+\mathrm{sin}(\frac{\pi x}{a})\xb7\mathrm{sin}(\frac{\pi y}{a})\xb7{\mathrm{MTF}}_{1\text{-}\mathrm{D}}(\frac{\sqrt{2}}{2a}\left)\right].$$## (8)

$${I}_{\text{output}}^{2\text{-}\mathrm{D}}=\frac{1}{2}[1+\mathrm{sin}(\frac{\pi x}{a})\xb7{\mathrm{MTF}}_{1\text{-}\mathrm{D}}(\frac{\sqrt{2}}{2a}\left)\right].$$Moreover, in ideal systems (aberration-free), we build the relationship between the 1-D and the 2-D MTF theoretically, shown in Eq. (9), which is also true in the vertical direction ($\widehat{y}$).

## (9)

$${\mathrm{MTF}}_{\widehat{x}}^{2\text{-}\mathrm{D}}(f)={\mathrm{MTF}}_{\widehat{x}}^{1\text{-}\mathrm{D}}\left(\sqrt{2}f\right).$$Furthermore, we can deduce from Eq. (9) that the 2-D MTF values are lower than the 1-D MTF values in the corresponding direction and the relationship of their cutoff frequency is

where ${f}_{\mathrm{CF}}^{2\text{-}\mathrm{D}}$ and ${f}_{\mathrm{CF}}^{1\text{-}\mathrm{D}}$ are the cutoff frequencies of 2-D and 1-D MTF, respectively. This conclusion is made for the case of ideal systems. In the following sections, we will discuss the values of the experimental result.## 2.2.

### Random Target Method and Noise Reduction Approach

As stated above, a 2-D target can provide comprehensive information for the 2-D MTF measurement of optical imaging systems in 2-D imaging scenes. However, a selective 2-D target will improve the measurement efficiency of the 2-D frequency. It is appropriate to use random images because, in addition to the advantages stated in Sec. 1, one can easily obtain its 2-D frequency spectrum through 2-D fast Fourier transform (FFT). A pixels random image is shown in Fig. 3 and its 1-D and 2-D PSD are shown in Figs. 4(a) and 4(b).

Random images are generated by a random number generator. Although a random target method is performed, a series of random images are displayed in a liquid crystal display (LCD) one by one which are imaged by the OSUT on a CCD. By calculating the ratio of the output image captured by the CCD to the input random image’s PSD, one can get the system’s MTF. The mathematical principle of this process could be described as

## (11)

$${\mathrm{PSD}}_{\text{output}}={[{\mathrm{MTF}}_{\mathrm{sys}}\times {\mathrm{MTF}}_{\text{test}}]}^{2}\times {\mathrm{PSD}}_{\text{input}},$$Here, to reduce the influence of temporal and spatial noises of the experimental setup, ${\mathrm{PSD}}_{\text{output}}$ is determined by

## (12)

$${\mathrm{PSD}}_{\text{output}}={\mathrm{PSD}}_{\text{captured}}-{\mathrm{PSD}}_{\text{system}},$$In 1-D and 2-D MTF measurements, we should calculate the PSDs of the input and output images, respectively. Then, Eq. (11) becomes

## (13)

$${\mathrm{PSD}}_{\text{output}}^{1\text{-}\mathrm{D}}(f)={[{\mathrm{MTF}}_{\mathrm{sys}}^{1\text{-}\mathrm{D}}(f)\xb7{\mathrm{MTF}}_{\text{test}}^{1\text{-}\mathrm{D}}(f)]}^{2}\xb7{\mathrm{PSD}}_{\text{input}}^{1\text{-}\mathrm{D}}(f).$$## (14)

$${\mathrm{PSD}}_{\text{output}}^{2\text{-}\mathrm{D}}(\xi ,\eta )={[{\mathrm{MTF}}_{\mathrm{sys}}^{2\text{-}\mathrm{D}}(\xi ,\eta )\xb7{\mathrm{MTF}}_{\text{test}}^{2-\mathrm{D}}(\xi ,\eta )]}^{2}\phantom{\rule{0ex}{0ex}}\xb7{\mathrm{PSD}}_{\text{input}}^{2\text{-}\mathrm{D}}(\xi ,\eta ).$$1-D and 2-D PSD can be obtained by applying 1-D FFT and 2-D FFT to the images. At the same time, the noise reduction method is applied to obtain a smooth 1-D MTF curve and 2-D MTF surface by taking average of the FFT results for a large number of images. Equations (15) and (16) describe the calculation process in detail

## (15)

$${\mathrm{PSD}}_{1\text{-}\mathrm{D}}=\frac{1}{{M}_{1}}\sum _{i=1}^{{M}_{1}}\{\mathrm{abs}{[\mathrm{FFT}({A}_{i})]}^{2}/{N}_{1}\},$$## (16)

$${\mathrm{PSD}}_{2\text{-}\mathrm{D}}=\frac{1}{{M}_{2}}\sum _{j=1}^{{M}_{2}}\{\mathrm{abs}{[\mathrm{FFT}2({B}_{j})]}^{2}/{N}_{2}^{2}\},$$We also discuss the relationship between the variances and peak to peak values of ${\mathrm{PSD}}_{1\text{-}\mathrm{D}}$ and ${M}_{1}$, and the relationship between the variances and peak to peak values of ${\mathrm{PSD}}_{2\text{-}\mathrm{D}}$ and the ${M}_{2}$, as are shown in Fig. 5. The variances are stable when increasing the number of vectors and matrices, which implies that ${M}_{1}=1024$ vectors and ${M}_{2}=1024$ matrices are sufficient for 1-D and 2-D MTF analysis. The average results of 1-D and 2-D PSD are shown in Figs. 6(a) and 6(b).

## 3.

## Experiment and Results

## 3.1.

### Experimental Setup and Experimental Procedure

The schematic diagram of the experimental setup is shown in Fig. 7. A host with two LCDs is adopted. Monitor 1 is used to show the interactive interfaces of applications and monitor 2 is used to present the random images. Thus the parameters are defined, e.g., $l$ and $d$ are the pixel size of the LCD and CCD, respectively, and $M$ is the reduction ratio of the optical system.

According to the sampling theorem, to avoid aliasing effects one should ensure ${f}_{\mathrm{max}}^{s}<{f}_{N}^{c}$ since the highest spatial frequency of the LCD measured in the image plane is ${f}_{\mathrm{max}}^{s}=1/(2\times l\times M)$ and the Nyquist frequency of the CCD is ${f}_{N}^{c}=1/(2\times d)$.

In our experiment, monitors 1 and 2 are Samsung SyncMaster E1920W LCDs whose pixel size is 0.2835 mm with a total of $1440\times 900\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{pixels}$. The CCD placed in the image plane is a ¼ in. Sony ICX618 monochrome CCD. Its total number of pixels is $656\times 492$ with each pixel size being $5.6\text{\hspace{0.17em}}\text{\hspace{0.17em}}\mu \mathrm{m}\times 5.6\text{\hspace{0.17em}}\text{\hspace{0.17em}}\mu \mathrm{m}$. Therefore, according to the constraints, the optical system’s reduction ratio could not be smaller than 0.0197. In this paper, the reduction ratio of the OSUT is 0.02 and the highest frequency of the screen measured in the image plane is $1/(2\times 0.2835\times 0.02)\approx \phantom{\rule{0ex}{0ex}}88.2\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{cycles}/\mathrm{mm}$.

The block diagram of the experiment is shown in Fig. 8. First of all, the CCD captures a series of images without random images displayed in the screen. Also, these captured images will be used to calculate the temporal and spatial noises (${\mathrm{PSD}}_{\mathrm{sys}}$) of the experimental setup. Then, ${M}_{2}=1024$ random images are displayed on monitor 2 one by one. These $N\times N$ random images are generated in MATLAB by using the functions named rand and imshow, and we set $N=256$ here because it is precise enough for the measurement, see Ref. 13. The original random image and the image captured by the CCD are saved in the host. Finally, all the images are used to calculate the 2-D PSD, i.e., ${\mathrm{PSD}}_{\text{measured}}^{2\text{-}\mathrm{D}}$, and ${M}_{1}=({M}_{2}/N)$ images of them are used to calculate the 1-D PSD, i.e., ${\mathrm{PSD}}_{\text{measured}}^{1\text{-}\mathrm{D}}$. After acquiring these PSD values, we can calculate the values of the optical system’s 1-D and 2-D MTF based on Eqs. (13) and (14).

## 3.2.

### Results and Discussion

A random image on the screen and its image captured by the CCD are shown in Fig. 9. The 1-D MTF in the horizontal direction is shown in Fig. 10 and the 2-D MTF is shown in Fig. 11. The 1-D MTF and the curve of 2-D MTF in the ($\widehat{x}$) direction are plotted in Fig. 12.

In Fig. 12, one can easily note that for a practical system, the 2-D MTF values are lower than the 1-D MTF values at every spatial frequency point. The 1-D MTF’s value is 0.023 at the cutoff frequency point. If we set 0.023 as the cutoff value, then the cutoff frequency of the 2-D MTF is $65.9\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{cycles}/\mathrm{mm}$. Thus $88.2/65.9=1.34$, which approximates $\sqrt{2}$. These results are consistent with the conclusions for an ideal system.

According to Eq. (9), another curve, the ideal curve, is also illustrated in pink in Fig. 12. The value for the ideal curve at frequency $f$ is equal to that of the 1-D MTF at frequency $\sqrt{2}f$. Also, when comparing with the 2-D MTF curve (in blue), these two curves do not strictly obey the relationship in Eq. (9) except at one point, which is located at about $0.38{f}_{c}\approx 33.47\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{cycles}/\mathrm{mm}$. The 2-D MTF curve will be lower than the ideal curve when $f<0.38{f}_{c}$ and will be higher than it when $f>0.38{f}_{c}$. These results are not completely consistent with the conclusions for an ideal system when considering various aberrations in the practical systems since these aberrations will be reflected in the 2-D MTF measurement results.^{20}

## 4.

## Conclusions

The relationship between 1-D and 2-D MTF has been discussed in theory and investigated in an experiment based on the random target method. The theoretical and experimental results show that the 2-D MTF’s values and cutoff frequency are lower than their counterparts in 1-D MTF, and that 2-D MTF can provide more comprehensive information than the 1-D MTF does since it covers the 2-D frequency domain.

The major contributions of this work are that, to the best of our knowledge, we take the average of the 2-D FFT results of a sufficient number of random images; this is used as a noise reduction method in the 2-D MTF measurement of optical imaging systems for the first time. Also, we discuss the characteristics of 2-D MTF and 1-D MTF in terms of MTF values and cutoff frequency through comparing 2-D MTF with 1-D MTF in the same direction.

## Acknowledgments

This research was supported by the grant from the National Natural Science Foundation of China (Nos. 61275003 and 51327005), the Research Fund for Shenzhen Key Laboratory of LED Packaging (ZDSY20120619141243215), and Guangdong Project (No. 2012B091100014).

## References

## Biography

**Jiqiang Kang** received his master’s degree in precise instrumentation from Tsinghua University, Beijing, China, in 2014. He is currently a PhD candidate at Hong Kong University.

**Qun Hao** received her PhD degree in precise instrumentation from Tsinghua University, Beijing, China, in 1998. She is currently a professor at Beijing Institute of Technology. Her research interests are freeform surface measurement, target pointing, and tracking, on which she has published at least 100 papers.

**Xuemin Cheng** received her PhD degree in photoelectric engineering from Beijing Institute of Technology, Beijing, China, in 2004. She is currently an associate professor at Tsinghua University, an OMNERC member, and also a committee member of the Opto-electronic-mechanical Committee in the Chinese Optical Society. Her research interests are optimization algorithms, and intelligent design of optical systems, on which she has published at least 41 papers.