|
1.IntroductionModulation 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–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–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 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–13,15–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.Methodology2.1.MTF Review and Comparison of 1-D and 2-D MTFIn incoherent illumination systems, the optical transfer function (OTF) is the normalized autocorrelation of the system exit pupil, which can be expressed as where is the system’s exit pupil function, is wavelength, and is the distance from the exit pupil plane to the image plane. MTF is the absolute value of the system’s OTF, i.e., . 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 where is the spatial frequency, and . 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 where is the 1-D MTF value at spatial frequency .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 where and are the fundamental frequencies along the horizontal and vertical directions, and . 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.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 . 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 where is the 2-D MTF of the system at frequency . Based on Eq. (6), the output signal of the system is If we focus the distribution on the horizontal direction () at the center of a row, say , then Eq. (7) becomes Equations (7) and (8) give the two directions and the () 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 () direction of the 2-D output will influence the MTF result; this will not happen to the 1-D target when one analyzes the () 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 (). 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 and 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 ApproachAs 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 where output PSD and input PSD are the PSDs of the final images and the original images on the LCD. and 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, is determined by where is the PSD value measured with random images in the screen, while 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 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 where is a set of vectors, and means the ’th vector. is a set of matrices and represents the ’th matrix of . and are the numbers of vectors or matrices of and where abs is the function for obtaining an absolute value.We also discuss the relationship between the variances and peak to peak values of and , and the relationship between the variances and peak to peak values of and the , as are shown in Fig. 5. The variances are stable when increasing the number of vectors and matrices, which implies that vectors and 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 Results3.1.Experimental Setup and Experimental ProcedureThe 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., and are the pixel size of the LCD and CCD, respectively, and is the reduction ratio of the optical system. According to the sampling theorem, to avoid aliasing effects one should ensure since the highest spatial frequency of the LCD measured in the image plane is and the Nyquist frequency of the CCD is . In our experiment, monitors 1 and 2 are Samsung SyncMaster E1920W LCDs whose pixel size is 0.2835 mm with a total of . The CCD placed in the image plane is a ¼ in. Sony ICX618 monochrome CCD. Its total number of pixels is with each pixel size being . 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 . 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 () of the experimental setup. Then, random images are displayed on monitor 2 one by one. These random images are generated in MATLAB by using the functions named rand and imshow, and we set 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., , and images of them are used to calculate the 1-D PSD, i.e., . 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 DiscussionA 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 () 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 . Thus , which approximates . 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 is equal to that of the 1-D MTF at frequency . 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 . The 2-D MTF curve will be lower than the ideal curve when and will be higher than it when . 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.ConclusionsThe 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. AcknowledgmentsThis 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). ReferencesM. BornE. Wolf, Principles of Optics, 7th ed.CUP Archive, Cambridge
(1999). Google Scholar
X. Chenet al.,
“Sensor modulation transfer function measurement using band-limited laser speckle,”
Opt. Express, 16
(24), 20047
–20059
(2008). http://dx.doi.org/10.1364/OE.16.020047 OPEXFF 1094-4087 Google Scholar
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). http://dx.doi.org/10.1364/OL.33.000533 OPLEDP 0146-9592 Google Scholar
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). http://dx.doi.org/10.1080/01431161.2010.488254 IJSEDK 0143-1161 Google Scholar
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). http://dx.doi.org/10.1117/1.2430503 OPEGAR 0091-3286 Google Scholar
G. D. BoremanS. Yang,
“Modulation transfer function measurement using three-and four-bar targets,”
Appl. Opt., 34
(34), 8050
–8052
(1995). http://dx.doi.org/10.1364/AO.34.008050 APOPAI 0003-6935 Google Scholar
J. R. Rogers,
“Three-bar resolution versus MTF: how different can they be anyway?,”
Proc. SPIE, 7071 707109
(2008). http://dx.doi.org/10.1117/12.800403 PSISDG 0277-786X Google Scholar
X. Zhanget al.,
“Measuring the modulation transfer function of image capture devices: what do the numbers really mean?,”
Proc. SPIE, 8293 829307
(2012). http://dx.doi.org/10.1117/12.912989 PSISDG 0277-786X Google Scholar
E. MaromB. MilgromN. Konforti,
“Two-dimensional modulation transfer function: a new perspective,”
Appl. Opt., 49
(35), 6749
–6755
(2010). http://dx.doi.org/10.1364/AO.49.006749 APOPAI 0003-6935 Google Scholar
H. HaimN. KonfortiE. Marom,
“Performance of imaging systems analyzed with two-dimensional target,”
Appl. Opt., 51
(25), 5966
–5972
(2012). http://dx.doi.org/10.1364/AO.51.005966 APOPAI 0003-6935 Google Scholar
H. HaimN. KonfortiE. Marom,
“Optical imaging systems analyzed with a 2D template,”
Appl. Opt., 51
(14), 2739
–2746
(2012). http://dx.doi.org/10.1364/AO.51.002739 APOPAI 0003-6935 Google Scholar
A. Danielset al.,
“Random transparency targets for modulation transfer function measurement in the visible and infrared regions,”
Opt. Eng., 34
(3), 860
–868
(1995). http://dx.doi.org/10.1117/12.190433 OPEGAR 0091-3286 Google Scholar
S. M. Backmanet al.,
“Random target method for fast MTF inspection,”
Opt. Express, 12
(12), 2610
–2615
(2004). http://dx.doi.org/10.1364/OPEX.12.002610 OPEXFF 1094-4087 Google Scholar
A. H. LettingtonQ. H. Hong,
“Measurement of the discrete modulation transfer-function,”
J. Mod. Opt., 40
(2), 203
–212
(1993). http://dx.doi.org/10.1080/09500349314550221 JMOPEW 0950-0340 Google Scholar
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). http://dx.doi.org/10.1117/1.3476335 OPEGAR 0091-3286 Google Scholar
E. Levyet al.,
“Modulation transfer function of a lens measured with a random target method,”
Appl. Opt., 38
(4), 679
–683
(1999). http://dx.doi.org/10.1364/AO.38.000679 APOPAI 0003-6935 Google Scholar
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). http://dx.doi.org/10.2971/jeos.2013.13029 1990-2573 Google Scholar
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). http://dx.doi.org/10.1364/AO.39.000515 APOPAI 0003-6935 Google Scholar
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). http://dx.doi.org/10.1364/JOT.80.000294 JOTEE4 1070-9762 Google Scholar
J. Q. KangX. M. ChengQ. Hao,
“Wavefront aberration detection: a new perspective,”
Appl. Mech. Mater., 496 1164
–1168
(2014). http://dx.doi.org/10.4028/www.scientific.net/AMM.496-500 AMMPDO 1662-7482 Google Scholar
BiographyJiqiang 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. |