MPE Mathematical Problems in Engineering 1563-5147 1024-123X Hindawi Publishing Corporation 10.1155/2015/863201 863201 Research Article Evaluation of Geometrical Modulation Transfer Function in Optical Lens System http://orcid.org/0000-0002-8919-5427 Tsai Cheng-Mu Li Mo Graduate Institute of Precision Engineering National Chung Hsing University 250 Kuo Kuang Road, Taichung 402 Taiwan nchu.edu.tw 2015 12102015 2015 20 09 2014 17 12 2014 12102015 2015 Copyright © 2015 Cheng-Mu Tsai. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

This paper presents ray tracing algorithms to evaluate the geometrical modulation transfer function (GMTF) of optical lens system. There are two kinds of ray tracings methods that can be applied to help simulate the point spread function (PSF) in the image plane, for example, paraxial optics and real ray tracings. The paraxial optics ray tracing is used to calculate the first-order properties such as the effective focal length (EFL) and the entrance pupil position through less cost of computation. However, the PSF could have a large tolerance by only using paraxial optics ray tracing for simulation. Some formulas for real ray tracing are applied in the sagittal and tangential line spread function (LSF). The algorithms are developed to demonstrate the simulation of LSF. Finally, the GMTF is evaluated after the fast Fourier transform (FFT) of the LSF.

1. Introduction

Optical products are very popular for people today. The camera function has also been widespread applied for various circumstances [1, 2]. Lens designers encounter more and more challenges for developing a lens with high image quality. Fortunately, advance computer technology can help the designers to simulate and analyze lens systems before implementing a lens. The major concern problem is how to optimize the lens parameters such as curvatures and thicknesses to get high image quality, for example, a high image resolution. Several optimizations have been proposed to improve the aberrations of lens systems . One simulation result, called modulation transfer function (MTF) , can be applied to determine the image resolution.

The MTF is the amplitude term of optical transfer function (OTF) that is similar to the transfer function of linear system. The transfer function is regarded as a major characteristic in the linear system. We cannot derive or evaluate the output signal without being given the transfer function. There are some proposals that present their methods to approach the transfer function . A simple one is to receive the impulse response at output as input being an impulse signal. This impulse response is related to the transfer function . Using the same procedure, a point source is respected as the impulse signal to help estimate the image response in a lens system. The image of point source shown in the image plane is called the point spread function (PSF) that is the inverse Fourier transform of the OTF. Therefore, the MTF applied to determine the image resolution can be derived from the amplitude term of the Fourier transform of the PSF.

The MTF can be evaluated from either the geometrical optics or diffraction calculation . Ray tracing methods are widely applied to simulate the image response of object point source in the lens system [12, 13]. Using the irradiance model, Lin and Liu  presented a MTF computation method without counting the number of rays traveling to each grid. Tseng et al.  proposed skew ray tracing method to simulate the geometrical PSF through a 4 × 4 homogeneous coordinate transformation matrix. In this paper, we describe the relation between the lens system and linear system. Finding the impulse response of the linear system is respected in the lens system to estimate the image response when the object is a point source. The developed algorithm of paraxial optic is used to calculate the first-order optical properties such as effective focal length (EFL) and the position of entrance pupil. An algorithm with real ray tracing method [14, 15] is developed to estimate the image of an object point (i.e., geometrical PSF).

Although the geometrical PSF can be used to determine the system performance, the MTF is always investigated to look into the image resolution of lens system. By the Fourier transform of the geometrical PSF, the MTF can be calculated to analyze the resolution of image in the lens systems. The organization of this paper is as follows. The relation between the lens system and linear system is explained in Section 2. Section 3 presents the ray tracing method for a lens system. The algorithms are developed to simulate a lens system in Section 4. Conclusions are drawn in Section 5.

2. Impulse Response Model of Optical Lens System

The transfer function of a linear system is to represent the relation between input and output signals, as shown in Figure 1. When the signal x ( t ) is input into a linear system with transfer function h ( t ) , the output signal y ( t ) is obtained based on the convolution theorem; that is, (1) y t = - x t h t - τ d τ . Equation (1) is used to evaluate the output signal in the linear system when given the transfer function and the input signal. The OTF of a lens system is similar to the transfer function of linear system. Lens designers usually determine the resolution of designed lens through the amplitude part of the OTF (e.g., MTF). It is important to evaluate the OTF for lens design. The issue is how to investigate the OTF from given lens parameters such as curvatures, thicknesses, and glass refractive indexes.

The relation between input and output for a linear system.

From the linear system concept, we can obtain the transfer function by replacing the input x ( t ) signal as impulse δ ( t ) signal. Even though the transfer function h ( t ) is unknown, the output signal y ( t ) becomes h ( t ) when substituting the signal x ( t ) to the signal δ ( t ) into (1). The same procedure allows us to evaluate the OTF for lens systems. In a lens system, the object point (or point source) plays the role like the impulse signal. The image shown in image plane is the impulse response of lens system when the object is a point source. Figure 2 shows the relation between object points and their image responses. In general, the image response depends on the location of object in the lens system. The object point sitting on the optical axis always receives the image response with the spherical aberration. The image response with the coma aberration is always shown in the image plane as the object point at off-axis. Therefore, the OTFs are related to the field of view (FOV). The object points o 1 and o 2 shown in Figure 2 are located at different FOV, which result in the light propagating to the images i 1 and i 2 having optical path difference (OPD). The OPD would cause the OTF variation at the image plane. Therefore, lens designers usually investigate three to five FOV OTFs to determine the overall resolution of image. In order to evaluate the OTFs, the ray tracing methods are applied to simulate the image response.

The relation between the object points and their image responses.

3. Ray Tracing in Optical Lens System

A camera lens is used to catch an object with multiple points into its image plane. Each object point is like a point source to diverge its power toward the outside. In geometrical optics, the direction of wavefront propagation is regarded as the rays. The ray tracing is to investigate how the rays travel from an object point to image plane. There are two processes used in ray tracing, that is, refraction process and transfer process.

3.1. Refraction Process for Paraxial Optics

An optical lens system is usually assembled by multiple lens elements. In order to collect the rays from an object point at image plane, the materials in both sides of a surface must have different refractive indices, which results in the rays deviating at the conjunctive surface. Figure 3 shows an object point at optical axis, A ; a ray emerges that propagates from one material with refractive index n impinging the spherical surface with center point, C , and then is deviated into the other material with refractive n . The deviated ray must satisfy Snell’s law; that is, (2) n sin i = n sin i , where i and i are incident and refractive angles, respectively. When the angle i is small enough (e.g., i < 5 °), (2) can be approximated as (3) n i = n i . Equation (3) is called the paraxial approximation.

A ray traces from one material to the other material.

From the triangles PCA and PCA shown in Figure 3, the geometry gives i = α + u and i = α - u . Following the sign convention , the incident and refractive angles can be expressed as (4a) i = u - α , (4b) i = u - α . Substituting (4a) and (4b) into (3), we obtain (5) n u - α = n u - α . Under the paraxial approximation, we approach α - h / R . Substituting this approach into (5), we can get (6) n u - n u = - h n - n R . Equation (6) is used in refraction process for the paraxial optics. The refractive ray would strike another surface to make the next refraction process. Between two refraction processes, however, the transfer process needs to be investigated to find the height that the ray hits on the next surface.

3.2. Transfer Process for Paraxial Optics

After refraction process, a ray travels from one surface to another surface. There is only one material between two surfaces. Figure 4 shows a ray coming from an object point A impinges on the point P 1 of the first surface to make the refraction process with the angle u 1 . The ray would travel from the point P 1 to the point P 2 based on the angle u 1 . We are interested in the height between the optical axis and the point P 2 . In the paraxial approximation, the refraction angle can be expressed as (7a) u 1 = - h 1 l 1 , (7b) u 1 = - h 2 l 2 . From (7b), the h 2 is - l 2 u 1 . In general, the lens designer represents a lens system by using the curvature and thickness. Figure 4 shows the thickness between two surfaces is d = l 1 - l 2 . After some arrangement, we can obtain that the h 2 is (8) h 2 = h 1 + d u 1 . Equations (6) and (8) are applied to calculate the EFL and pupil position. However, the angles of real rays are usually not small enough to make the paraxial approximation, which results in the fact that we cannot apply the paraxial ray tracing to simulate the spot diagram. Instead, the real ray tracing is used to evaluate the lens performance.

A ray travels between two surfaces.

3.3. Real Ray Tracing for Transfer

The common lens system is to employ lens elements with spherical surfaces. Let us discuss a real ray transfer process between spherical surfaces . Figure 5 shows a ray with direction cosine ( L , M , N ) propagates from the front spherical surface to the back spherical surface. The two spherical surfaces have their local coordinate. The point the ray sits on the front surface is P - 1 = ( x - 1 , y - 1 , z - 1 ) . After the ray refracted the point P = ( x , y , z ) , the ray would arrive at a point P 0 = ( x 0 , y 0 , z 0 ) of the x - y plane in the coordinate of the back surface. We can find that the point P 0 ( x 0 , y 0 , z 0 ) is (9a) x 0 = x - 1 + L d - z - 1 N , (9b) y 0 = y - 1 + M d - z - 1 N , where d is the thickness between two spherical surfaces. The transfer process is to find the point P = ( x , y , z ) that the ray hits the back surface; that is, (10a) x = x 0 + L Δ , (10b) y = y 0 + M Δ , (10c) z = N Δ , where Δ will be discussed later.

A ray propagates between two spherical surfaces .

Figure 6 shows a sphere with the radius R . A point p lies on the surface of this sphere, and we are interested in the distance from the y -axis to the point p . From the Pythagorean theorem, we can find the distance z as (11) z = 1 2 c x 2 + y 2 + z 2 , where c is the radius curvature; for example, c = 1 / R . Substituting (10a) and (10b) into (11), we can obtain (12) 1 2 c x 0 + L Δ 2 + y 0 + M Δ 2 + N 2 Δ 2 = N Δ or (13) c 2 Δ 2 - 2 c N - c L x 0 + M y 0 Δ + c 2 x 0 2 + y 0 2 = 0 , where we define F = c ( x 0 2 + y 0 2 ) and G = N - c ( L x 0 + M y 0 ) . Finally Δ can be found as (14) Δ = G ± G 2 - c F 1 / 2 c or (15) Δ = F G + G 2 - c F . Equations (9a) and (9b), (10a), (10b), and (10c), and (15) are used to calculate the point P that the ray coming from the front surface arrives at the back surface. Equation (15) is used to develop the algorithm instead of (14). This is due to the fact that the term G 2 shown in (14) is much larger than the term c F . The numerical precision must be satisfied in (14) when running the simulation. If the numerical precision is not large enough, the simulation results will receive much tolerance. Instead of (14), (15) is applied to develop the algorithm for enhancing numerical precision.

A point p sits on the spherical surface.

3.4. Real Ray Tracing for Refraction

The refraction of real ray tracing follows Snell’s law. Figure 7 shows a ray with incident unit vector r respecting the normal line n is refracted to the refractive unit vector r . Let us represent Snell’s law by the vector cross product form ; that is, (16) n r × n = n r × n , where r and r are the unit vectors for the incident and refractive rays, respectively, and n is a unit vector to normalize the interface at a point that the ray is refracted. When n is made the cross product operator with (16), we can find (17) n r - n r = n cos i - n cos i n , where i and i are the incident and refractive angles, respectively. When the ray is refracted by a spherical surface, the normal vector n can be defined as (18) n = - c x , - c y , 1 - c z . Using (17) and (18), we can find the refractive vector r . The real ray tracing process is applied to determine the line spread function.

The presentation of Snell’s law on vector form.

4. Algorithm Development and Simulation Results

A lens system with three lens elements shown in Figure 8 was employed to find its LSF. Table 1 lists the lens parameters. We follow the sign convention in optical design. The distance measured from the vertex of surface to right is positive and conversely it is negative. If the center is located at the right of the surface, the radius is positive; for example, the surface 1 is 22.014. Since the center of surface 2 is located at its left, the radius is negative (e.g., −435.76). So far, the commercial optical software is still expensive so we developed the algorithms that can be used to evaluate the image resolution. Algorithm 1 explains the algorithm of the transfer for real ray tracing method while the algorithm of the refraction is described in Algorithm 2. These two algorithms were based on (9a), (9b), (10a), (10b), and (10c), (13), (15), and (17)-(18). There are 100 thousand rays emerging from the object point to pass through all lens elements into the image plane. Two rays, called marginal rays, from the object point would travel the rim of the entrance pupil. The angle formed by these two marginal rays would be equally divided into 100 thousand parts for initializing the direction cosines of rays for the real ray tracing process. The 100 thousand rays would pass all the lens elements and finally arrive at the image plane.

The simulation lens parameters.

Radius ( r ) Thickness (thi) Refractive index (nd)
Object Infinity Infinity 1
1 22.014 3.259 1.62
2 −435.76 6.008 1
3 −22.213 1 1.62
4 20.292 4.75 1
5 79.684 2.952 1.62
6 −18.390 42.208 1

<bold>Algorithm 1: </bold>The transfer algorithm for real ray tracing.

typeXYZ TProcess(typeXYZ p_1,typeLMN LMN,double d,double C,double cosI)

{

double x0,y0,z0;

double delta,F,G;

double tmp;

typeXYZ p;

x0=p_1.x+LMN.L/LMN.N (d-p_1.z); // x 0 = x - 1 + L ( d - z - 1 ) / N

y0=p_1.y+LMN.M/LMN.N (d-p_1.z); // y 0 = y - 1 + M ( d - z - 1 ) / N

z0=0;

F=C (x0 x0+y0 y0); // F = c ( x 0 2    + y 0 2 )

G=LMN.N-C (LMN.L x0+LMN.M y0); // G = N - c L x 0 + M y 0

tmp=G G-C F;

cosI=sqrt(tmp); // cos I = G 2 - C F 1 / 2

delta=F/(G+ cosI); // Δ = F / ( G + G 2 - c F 1 / 2    )

p.x=x0+LMN.L delta;     // x = x 0 + L Δ

p.y=y0+LMN.M delta;   // y = y 0 + M Δ

p.z=LMN.N delta;   // z = N Δ

return p;

}

<bold>Algorithm 2: </bold>The refraction algorithm for real ray tracing.

typeLMN RProcess(typeLMN LMN,typeXYZ p,double C,double n,double np,double cosI){

double cosIp;

double tmp;

double K;

typeLMN LMNp;

tmp=n/np;

tmp=1-tmp tmp (1-cosI cosI);

cosIp=sqrt(tmp); / / n cos I = + n 2 - n 2 1 - cos 2 I 1 / 2

K=(np cosIp-n cosI) C; / / K = c ( n cos I - n cos I )

tmp=n LMN.L-K p.x;

LMNp.L=tmp/np; / / n L = n L - K x

tmp=n LMN.M-K p.y;

LMNp.M=tmp/np; // n M = n M - K y

tmp=n LMN.N-K p.z+np cosIp-n cosI;

LMNp.N=tmp/np; n N = n N - K z + n cos I - n cos I

return LMNp;

}

The structure of the simulation lens system.

We collect these tracing rays on image plane to investigate the LSF. Sampling the rays on the image plane can count the number of the rays to imitate the power distribution. The LSP is sampled based on taking Δ S units apart to consist of M samples so that the total collected space is required as (19) S = M Δ S . And its corresponding discrete frequency domain is (20) Δ u = 1 M Δ S = 1 S . The fast Fourier transform (FFT) will be used to calculate the lens transform function. Therefore, we are concerned about the entire frequency range spanned by the M components of the FFT; that is, (21) Ω = M Δ u = 1 Δ S . The MTF in the proposed simulation lens is less than 500 cycles/mm. When S is set to 1 mm and M is 10000, we can evaluate the MTF with 10000 cycles/mm which is large enough to usual cases. Therefore, the S and M are set to 1 mm and 10000, respectively. Figure 9 shows the LSF as object point on the optical axis. The tangential and sagittal LSFs with fields 14° and 20° are shown in Figures 10 and 11, respectively. As we can see, the object point placed on the optical axis results in the image spot being symmetrical, making the spherical aberration. The sagittal and tangential images are different when the object point locates at the off-axis, which induces the coma aberration, as shown in Figures 10 and 11. Although the LSFs look very similar in Figures 9, 10(b), and 11(a), the resolutions are very different from the MTFs. Figure 12 shows the MTFs after running the FFT of the LSFs. Figure 13 shows the MTF results from the commercial software Zemax. The MTF shown in Figure 13 is acceptable around 25 cycles/mm, which is the same as the proposed simulation shown in Figure 12. The results shown in Figure 12 have x - and y -axis. As a result, from Figures 9 and 10, the LSF is symmetrical for x -axis while the LSF is asymmetrical for y -axis. Therefore, the MTFs are different between x -axis and y -axis. Besides, the transform function is related to the field angle. These two things make the MTF look very different. The MTF of the field 0° is the worst in Figure 12 due to the fact that the image of point source at field 0° has less power than others at the center of image, which is lower resolution than that of other fields. Figure 10(a) shows the power of the LSF is almost at the image center, resulting in the fact that the greatest MTF is at the x -axis of the field 14° to receive the best resolution. As we mentioned early, the OTF is related to the position of object point so that the image responses would be different at the image plane.

The simulated LSF at field 0°.

The simulated LSF at field 14°: (a) x -axis, (b) y -axis.

The simulated LSF at field 20°: (a) x -axis, (b) y -axis.

Modulation transfer function by proposed simulation.

Modulation transfer function by using Zemax.

5. Conclusions

This paper presents that the transfer function of a lens system is similar to that of a linear system. A transfer function of the linear system can be evaluated through inputting an impulse signal, making the fact that the output would be the transfer function. The same concept of linear system is used to respect a point source at object space as the impulse signal to simulate its response at image space. The paraxial ray tracing is applied to calculate the essential lens properties like the EFL and pupil position. However, the image of a point source would reveal a large tolerance when only using paraxial optics ray tracing for simulation. The real ray tracing method is applied to evaluate the LPF. The algorithms for ray tracing were also developed to simulate the LPF of a lens system. After the FFT process, we can evaluate the MTF to determine the resolution of a lens system.

Conflict of Interests

There is no conflict of interests regarding the publication of this manuscript.

Liang M. Huang X. Chen C.-H. Zheng G. Tokuta A. Robust calibration of cameras with telephoto lens using regularized least squares Mathematical Problems in Engineering 2014 2014 9 689429 10.1155/2014/689429 2-s2.0-84896298371 Wang L. Duan F. Lv K. Fisheye-lens-based visual sun compass for perception of spatial orientation Mathematical Problems in Engineering 2012 2012 15 460430 10.1155/2012/460430 2-s2.0-84859736477 Fang Y.-C. Tsai C.-M. MacDonald J. Pai Y.-C. Eliminating chromatic aberration in Gauss-type lens design using a novel genetic algorithm Applied Optics 2007 46 13 2401 2410 10.1364/AO.46.002401 2-s2.0-34250320515 Fang Y. C. Tsai C. M. Miniature lens design and optimization with liquid lens element via genetic algorithm Journal of Optics A: Pure and Applied Optics 2008 10 7 9 075304 10.1088/1464-4258/10/7/075304 2-s2.0-47249102925 Chen C. C. Tsai C. M. Fang Y. C. Optical design of LCOS optical engine and optimization with genetic algorithm Journal of Display Technology 2009 5 8 293 305 10.1109/jdt.2009.2021540 2-s2.0-68849083770 Zhang W. Aljasem K. Zappe H. Seifert A. Highly flexible MTF measurement system for tunable micro lenses Optics Express 2010 18 12 12458 12469 2-s2.0-77953606025 Levy E. Peles D. Opher-Lipson M. Lipson S. G. Modulation transfer function of a lens measured with a random target method Applied Optics 1999 38 4 679 683 10.1364/ao.38.000679 2-s2.0-0033072592 de la Barrière F. Druart G. Guérineau N. Taboury J. Primot J. Deschamps J. Modulation transfer function measurement of a multichannel optical system Applied Optics 2010 49 15 2879 2890 10.1364/ao.49.002879 2-s2.0-77956005938 Baddour N. Multidimensional wave field signal theory: transfer function relationships Mathematical Problems in Engineering 2012 2012 27 478295 10.1155/2012/478295 MR2853325 2-s2.0-84855594181 Chen K. F. Li Y. F. On the integration schemes of retrieving impulse response functions from transfer functions Mathematical Problems in Engineering 2010 2010 9 143582 10.1155/2010/143582 MR2764978 Li M. Lim S. C. Chen S. Exact solution of impulse response to a class of fractional oscillators and its stability Mathematical Problems in Engineering 2011 2011 9 657839 10.1155/2011/657839 MR2740336 2-s2.0-79251547368 Lin P. D. Liu C. S. Geometrical MTF computation method based on the irradiance model Applied Physics B: Lasers and Optics 2011 102 1 243 249 10.1007/s00340-010-4349-3 2-s2.0-79851509389 Tseng K.-H. Kung C. Liao T.-T. Chang H.-P. Calculation of modulation transfer function of an optical system by using skew ray tracing Transactions of the Canadian Society for Mechanical Engineering 2009 33 3 429 442 2-s2.0-77955784593 Wang Y. Hopkins H. H. Ray-tracing and aberration formulate for a general optical system Journal of Modern Optics 1992 39 9 1897 1938 10.1080/09500349214551981 Welford W. T. Aberrations of Optical Systems 1986 CRC Press Mouroulis P. Macdonald J. Geometrical Optics and Optical Design 1997 Oxford, UK Oxford University Press