30 October 2014 Measurement and comparison of one- and two-dimensional modulation transfer function of optical imaging systems based on the random target method
Author Affiliations +
Optical Engineering, 53(10), 104105 (2014). doi:10.1117/1.OE.53.10.104105
Abstract
One-dimensional modulation transfer function (1-D MTF) has been generally calculated to evaluate the image quality of optical imaging systems, such as the horizontal MTF and vertical MTF. These MTFs can be measured by the use of some mature ways. However, the information of 1-D MTF for performance evaluation may not enough for the systems handling two-dimensional (2-D) targets of high resolution, thus discussing 2-D MTF will be necessary. We investigate the measurement method for the 1-D and 2-D MTF of optical imaging systems based on the random target method, and the characteristics of 2-D MTF and 1-D MTF in terms of MTF values and cutoff frequency are also noted.
Kang, Hao, and Cheng: Measurement and comparison of one- and two-dimensional modulation transfer function of optical imaging systems based on the random target method

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.12.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,67.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/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.910.11.12.13,1516.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)

OTF(ξ,η)=P(x,y)P*(x+λziξ,y+λziη)dxdyP(x,y)dxdy,
where P(x,y) is the system’s exit pupil function, λ is wavelength, and zi is the distance from the exit pupil plane to the image plane. MTF is the absolute value of the system’s OTF, i.e., MTF=|OTF|. This is the rigorous definition of MTF. Although one rarely measures MTF through OTF because the OTF test itself is complicated, there is no denying that it is plausible.

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

(2)

Iinput1-D=12[1+sin(2πxf)],
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)

Ioutput1-D=12[1+sin(2πxf)·MTF1-D(f)]=12[1+sin(πxa)·MTF1-D(12a)],
where MTF1-D(f) is the 1-D MTF value at spatial frequency f.

Fig. 1

Typical sinusoidal targets: (a) one-dimensional (1-D) target, (b) two-dimensional (2-D) target.

OE_53_10_104105_f001.png

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)

Iinput2-D=12[1+sin(2πxfx)·sin(2πyfy)],
where fx and fy are the fundamental frequencies along the horizontal and vertical directions, and fx=fy=1/(2a). As for the 2-D target, we also need to obtain the output signals when input signals transmit through the optical system by analyzing its 2-D spectrum with 2-D Fourier transform, as in Eq. (5). Its spectrum is shown in Fig. 2.

(5)

H(ξ,η)=12{1+14[δ(ξfx)δ(ξ+fx)]·[δ(η+fy)δ(ηfy)]}=12{1+14[δ(ξ+fx,ηfy)+δ(ξfx,η+fy)δ(ξ+fx,η+fy)δ(ξfx,ηfy)]}.

Fig. 2

2-D frequency spectrum of a 2-D sinusoidal target.

OE_53_10_104105_f002.png

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=(fx2+fy2)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)

MTF2-D(fx,fy)=MTF1-D(fx2+fy2),
where MTF2-D(fx,fy) is the 2-D MTF of the system at frequency (fx,fy). Based on Eq. (6), the output signal of the system is

(7)

Ioutput2-D=12[1+sin(2πxfx)·sin(2πyfy)·MTF2-D(fx,fy)]=12[1+sin(πxa)·sin(πya)·MTF1-D(22a)].
If we focus the distribution on the horizontal direction (x^) at the center of a row, say y=a/2, then Eq. (7) becomes

(8)

Ioutput2-D=12[1+sin(πxa)·MTF1-D(22a)].
Equations (7) and (8) give the two directions and the (x^) direction’s expressions, respectively, of the 2-D output of the system. By comparing Eq. (3) with Eq. (8), one can note that although both of them have the same input expression, they are modulated differently which leads to different outputs. As the image brightness of a 2-D target changes in two orthogonal directions and a 1-D target changes in only one direction, the brightness change in the (y^) direction of the 2-D output will influence the MTF result; this will not happen to the 1-D target when one analyzes the (x^) direction’s MTF with the output images. Actually, this interaction of two orthogonal directions is more in line with the actual situation as most systems are used in a 2-D imaging scene which has abundant spatial frequencies distributed in 2-D frequency domain. Therefore, to evaluate the true resolving power of imaging systems, one should apply 2-D targets and adopt 2-D MTF as they can provide comprehensive modulation information about the imaging system.

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 (y^).

(9)

MTFx^2-D(f)=MTFx^1-D(2f).

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

(10)

fCF2-D=fCF1-D/2,
where fCF2-D and fCF1-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).

Fig. 3

A 256×256pixel random image.

OE_53_10_104105_f003.png

Fig. 4

PSD of a random image: (a) 1-D PSD and (b) 2-D PSD.

OE_53_10_104105_f004.png

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)

PSDoutput=[MTFsys×MTFtest]2×PSDinput,
where output PSD and input PSD are the PSDs of the final images and the original images on the LCD. MTFsys and MTFtest are the system inherent MTF of the system devices (including LCD and CCD) and the OSUT.

Here, to reduce the influence of temporal and spatial noises of the experimental setup, PSDoutput is determined by

(12)

PSDoutput=PSDcapturedPSDsystem,
where PSDcaptured is the PSD value measured with random images in the screen, while PSDsystem is the PSD value measured without random images in the screen.

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)

PSDoutput1-D(f)=[MTFsys1-D(f)·MTFtest1-D(f)]2·PSDinput1-D(f).

(14)

PSDoutput2-D(ξ,η)=[MTFsys2-D(ξ,η)·MTFtest2D(ξ,η)]2·PSDinput2-D(ξ,η).

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)

PSD1-D=1M1i=1M1{abs[FFT(Ai)]2/N1},

(16)

PSD2-D=1M2j=1M2{abs[FFT2(Bj)]2/N22},
where A is a set of 1×N1 vectors, and Ai means the i’th vector. B is a set of N2×N2 matrices and Bj represents the j’th matrix of B. M1 and M2 are the numbers of vectors or matrices of A and B where abs is the function for obtaining an absolute value.

We also discuss the relationship between the variances and peak to peak values of PSD1-D and M1, and the relationship between the variances and peak to peak values of PSD2-D and the M2, as are shown in Fig. 5. The variances are stable when increasing the number of vectors and matrices, which implies that M1=1024 vectors and M2=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).

Fig. 5

Relationship between pp value (normalized at 1 vector or matrix) and variance with the number of vectors and matrices.

OE_53_10_104105_f005.png

Fig. 6

Average results: (a) 1-D PSD averaged on 1024 vectors and (b) 2-D PSD averaged on 1024 matrixes.

OE_53_10_104105_f006.png

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.

Fig. 7

Schematic diagram of the experimental setup.

OE_53_10_104105_f007.png

According to the sampling theorem, to avoid aliasing effects one should ensure fmaxs<fNc since the highest spatial frequency of the LCD measured in the image plane is fmaxs=1/(2×l×M) and the Nyquist frequency of the CCD is fNc=1/(2×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×900pixels. The CCD placed in the image plane is a ¼ in. Sony ICX618 monochrome CCD. Its total number of pixels is 656×492 with each pixel size being 5.6μm×5.6μ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×0.2835×0.02)88.2cycles/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 (PSDsys) of the experimental setup. Then, M2=1024 random images are displayed on monitor 2 one by one. These N×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., PSDmeasured2-D, and M1=(M2/N) images of them are used to calculate the 1-D PSD, i.e., PSDmeasured1-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).

Fig. 8

Block diagram of experimental procedures.

OE_53_10_104105_f008.png

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 (x^) direction are plotted in Fig. 12.

Fig. 9

(a) Random image on the screen. (b) Image captured by the CCD.

OE_53_10_104105_f009.png

Fig. 10

Raw data of 1-D MTF and curve fitted with polynomials.

OE_53_10_104105_f010.png

Fig. 11

2-D MTF.

OE_53_10_104105_f011.png

Fig. 12

The comparison of 1-D and 2-D MTF.

OE_53_10_104105_f012.png

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.9cycles/mm. Thus 88.2/65.9=1.34, which approximates 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 2f. 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.38fc33.47cycles/mm. The 2-D MTF curve will be lower than the ideal curve when f<0.38fc and will be higher than it when f>0.38fc. 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

1. 

M. BornE. Wolf, Principles of Optics, 7th ed., CUP Archive, Cambridge (1999).Google Scholar

2. 

X. Chenet al., “Sensor modulation transfer function measurement using band-limited laser speckle,” Opt. Express 16(24), 20047–20059 (2008).OPEXFF1094-4087http://dx.doi.org/10.1364/OE.16.020047Google Scholar

3. 

Y. Zhanget al., “Dynamic modulation transfer function: a method to characterize the temporal performance of liquid-crystal displays,” Opt. Lett. 33(6), 533–535 (2008).OPLEDP0146-9592http://dx.doi.org/10.1364/OL.33.000533Google Scholar

4. 

X. B. Liet al., “An analysis of the knife-edge method for on-orbit MTF estimation of optical sensors,” Int. J. Remote Sens. 31(17–18), 4995–5010 (2010).IJSEDK0143-1161http://dx.doi.org/10.1080/01431161.2010.488254Google Scholar

5. 

J. T. OlsonR. L. EspinolaE. L. Jacobs, “Comparison of tilted slit and tilted edge superresolution modulation transfer function techniques,” Opt. Eng. 46(1), 016403 (2007).OPEGAR0091-3286http://dx.doi.org/10.1117/1.2430503Google Scholar

6. 

G. D. BoremanS. Yang, “Modulation transfer function measurement using three-and four-bar targets,” Appl. Opt. 34(34), 8050–8052 (1995).APOPAI0003-6935http://dx.doi.org/10.1364/AO.34.008050Google Scholar

7. 

J. R. Rogers, “Three-bar resolution versus MTF: how different can they be anyway?,” Proc. SPIE 7071, 707109 (2008).PSISDG0277-786Xhttp://dx.doi.org/10.1117/12.800403Google Scholar

8. 

X. Zhanget al., “Measuring the modulation transfer function of image capture devices: what do the numbers really mean?,” Proc. SPIE 8293, 829307 (2012).PSISDG0277-786Xhttp://dx.doi.org/10.1117/12.912989Google Scholar

9. 

E. MaromB. MilgromN. Konforti, “Two-dimensional modulation transfer function: a new perspective,” Appl. Opt. 49(35), 6749–6755 (2010).APOPAI0003-6935http://dx.doi.org/10.1364/AO.49.006749Google Scholar

10. 

H. HaimN. KonfortiE. Marom, “Performance of imaging systems analyzed with two-dimensional target,” Appl. Opt. 51(25), 5966–5972 (2012).APOPAI0003-6935http://dx.doi.org/10.1364/AO.51.005966Google Scholar

11. 

H. HaimN. KonfortiE. Marom, “Optical imaging systems analyzed with a 2D template,” Appl. Opt. 51(14), 2739–2746 (2012).APOPAI0003-6935http://dx.doi.org/10.1364/AO.51.002739Google Scholar

12. 

A. Danielset al., “Random transparency targets for modulation transfer function measurement in the visible and infrared regions,” Opt. Eng. 34(3), 860–868 (1995).OPEGAR0091-3286http://dx.doi.org/10.1117/12.190433Google Scholar

13. 

S. M. Backmanet al., “Random target method for fast MTF inspection,” Opt. Express 12(12), 2610–2615 (2004).OPEXFF1094-4087http://dx.doi.org/10.1364/OPEX.12.002610Google Scholar

14. 

A. H. LettingtonQ. H. Hong, “Measurement of the discrete modulation transfer-function,” J. Mod. Opt. 40(2), 203–212 (1993).JMOPEW0950-0340http://dx.doi.org/10.1080/09500349314550221Google Scholar

15. 

A. Fernandez-OliverasA. M. PozoM. Rubino, “Comparison of spectacle-lens optical quality by modulation transfer function measurements based on random-dot patterns,” Opt. Eng. 49(8), 083603 (2010).OPEGAR0091-3286http://dx.doi.org/10.1117/1.3476335Google Scholar

16. 

E. Levyet al., “Modulation transfer function of a lens measured with a random target method,” Appl. Opt. 38(4), 679–683 (1999).APOPAI0003-6935http://dx.doi.org/10.1364/AO.38.000679Google Scholar

17. 

F. A. Navas-Moyaet al., “Measurement of the optical transfer function using a white-dot pattern presented on a liquid-crystal display,” J. Eur. Opt. Soc. 8, 13029 (2013).1990-2573http://dx.doi.org/10.2971/jeos.2013.13029Google Scholar

18. 

B. T. TeipenD. L. MacFarlane, “Liquid-crystal-display projector-based modulation transfer function measurements of charge-coupled-device video camera systems,” Appl. Opt. 39(4), 515–525 (2000).APOPAI0003-6935http://dx.doi.org/10.1364/AO.39.000515Google Scholar

19. 

N. N. EvtikhievV. V. KrasnovS. N. Starikov, “A method of generating amplitude masks with a constant power spectra and using them to measure the two-dimensional modulation-transfer functions of optical systems,” J. Opt. Technol. 80(5), 294–300 (2013).JOTEE41070-9762http://dx.doi.org/10.1364/JOT.80.000294Google Scholar

20. 

J. Q. KangX. M. ChengQ. Hao, “Wavefront aberration detection: a new perspective,” Appl. Mech. Mater. 496, 1164–1168 (2014).AMMPDO1662-7482http://dx.doi.org/10.4028/www.scientific.net/AMM.496-500Google Scholar

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.

Jiqiang Kang, Qun Hao, Xuemin Cheng, "Measurement and comparison of one- and two-dimensional modulation transfer function of optical imaging systems based on the random target method," Optical Engineering 53(10), 104105 (30 October 2014). http://dx.doi.org/10.1117/1.OE.53.10.104105
Submission: Received ; Accepted
JOURNAL ARTICLE
7 PAGES


SHARE
KEYWORDS
Modulation transfer functions

Imaging systems

Charge-coupled devices

Optical imaging

Spatial frequencies

LCDs

Modulation

RELATED CONTENT

MTF-measurement-of-LCDs-by-a-linear-CCD-imager...
Proceedings of SPIE (November 01 1997)
The-AWMC-metric-for-zoom-lens-assessment
Proceedings of SPIE (October 19 2012)
Evaluation-of-color-display-by-polychromatic-MTF
Proceedings of SPIE (April 26 2002)

Back to Top