A pupil filter is a useful technique for modifying the light intensity distribution near the focus of an optical system to realize depth of field (DOF) extension and superresolution. In this paper, we proposed a new design of the phase only pupil filter by using Zernike polynomials. The effect of design parameters of the new filters on DOF extension and superresolution are discussed, such as defocus Strehl ratio (S.R.), superresolution factor (G) and relative first side lobe intensity (M). In comparison with the other two types of pupil filters, the proposed filter presents its advantages on controlling both the axial and radial light intensity distribution. Finally, defocused imaging simulations are carried out to further demonstrate the effectiveness and superiority of the proposed pupil filter on DOF extension and superresolution in an optical imaging system. In numerous applications it is necessary to extend the depth of field (DOF) and improve the resolution for an optical system, such as microlithography [1] , optical data storage [2] , and microscopy imaging [3 , 4] . One of the most popular techniques to achieve DOF extension and superresolution includes introducing a pupil filter at the pupil plane of an optical system to modify the radial and axial spot size of the point spread function (PSF). Since Toraldo made a seminal study on a radial-symmetric filter [5] , a lot of work has been done on the design of a pupil filter. Although an amplitude only filter can provide effective superresolution and DOF extension, it still has two main difficulties. The first one is the reduction of image brightness and the other refers to the fabrication, especially for continuously varying amplitude functions [6 , 7] . Therefore, in recent years phase only filter design has achieved a rapid development based on the parabolic approximation of the PSF for the general filters introduced by Sheppard [8 , 9] . Typically, the design strategies of pupil filters include circular obstructions [10 - 12] , continuously varying functions, such as Bessel functions [13 , 14] and Gaussian functions [15] , and special forms, such as a spoke wheel filter [16] . In this paper, a phase only filter is proposed by using Zernike polynomials. Radial and axial spot sizes of an optical system with and without the designed filter are analyzed relatively based on scalar diffraction theory. Then, compared with the other two typical phase filters (three-zones and Gaussian functions), we show its advantages on DOF extension and superresolution. Finally, we further demonstrate the availability and superiority of our proposed phase only filter through defocus imaging simulation. According to the Born and Wolf’s theory [17] , at the focal region, the amplitude distribution of a converging monochromatic spherical wavefront passing through the pupil plane can be written as: Where ρ is the normalized radial coordinate of the pupil plane, P represents the pupil function, J₀ represents Bessel function of the first kind of order zero, u and v are the radial and axial dimensionless optical coordinates, given by: Where NA is the numerical aperture of the optical system, λ is the wavelength of the incident light, r and z are the usual radial and axial distances. Since most practical imaging systems use incoherent illumination, the light intensity distribution near the focal region can be calculated by its amplitude distribution, written as: Where * denotes the complex conjugate. After introducing the phase pupil filter into the optical system, the light intensity distribution will be modulated by the phase changing of the pupil function to realize the DOF extension and superresolution. These can be evaluated by the following three parameters: These parameters can be written as: Where z_defocus is the defocus distance, I_max (0, z_defocus ) is the intensity at the focal point with defocus distance z_defocus , I_max(0, 0) is the intensity at the focal point of without defocus, R_{central lobe} is the radius of central lobe of intensity, v_{first-side lobe} is the radial optical coordinate of maximum intensity at the first side lobe. The pupil functions of an optical system with phase only pupil filter can be written as: Where Ψ(ρ) is phase function of the pupil filter. Compared with conventional power functions or Gaussian functions, there are some advantages for pupil filter design by using Zernike functions since they are closely related with optical aberrations. The phase functions with Zernike functions now can be written as: Where K is the diffraction order, here set to 1, N is the number of Zernike coefficients, Z_i is the i^th Zernike standard polynomial and A_i is the corresponding coefficient, ρ is the normalized radial coordinate of the pupil plane and φ is the angular coordinate. The first 16 polynomials of the Zernike functions are shown in Table 1 [18] . Here 1^st, 4^th, 5^th, 9^th, 12^th, 16^th polynomials of the Zernike functions are chosen to build the phase functions of the pupil filter by considering both availability for spot size modulation and the ease of fabrication. According to the relationships between Zernike functions and Seidel aberrations, roles that these polynomials play in the phase function can be analyzed reasonably. The optimization process is described as follows: A series of coefficients of chosen Zernike polynomials is used to calculate the PSF of the optical system with the phase pupil filter, and the M can be obtained from the calculated PSF. If the value of M is larger than the threshold (here the threshold is 0.2), select the next series of coefficients to re-calculate the PSF since the optical system is no longer suitable for imaging with a bigger M . If the value of M is not larger than the threshold, continue calculating defocus Strehl ratio S.R . and superresolution gain G of the optical system. Then, choose proper weight coefficients κ₁ and κ₂ to calculate the comprehensive evaluation result of this series of coefficients of Zernike polynomials. The coefficients κ₁ and κ₂ are used to balance super-resolution and DOF extension. In this paper, we are more concerned about DOF extension, so the κ₁ equals 0.7 and κ₂ equals 0.3. Repeat the above steps, selecting the best comprehensive evaluation result. The corresponding coefficients are the optimization results. It can be presented as follows by C-type codes: Iteration: for{( A_i ) = Coef( i ), i ++} end ( A_i ) _opt = Coef ( Maxima ( Value )) In the optimization, Coef is a matrix containing all the combinations of coefficients within a certain range, Value_sr and Value_g are used to keep the defocus Strehl ratio S.R . and superresolution gain G with different coefficients under M constraint, κ₁ and κ₂ are the weight coefficients. Thus, optimized coefficients of the proposed Zernike phase only pupil filter can be obtained by a global optimization algorithm, as shown below: A ₁ = 0.250, A ₄ = −0.233, A ₅ = −0.126, A ₉ = 0.317, A ₁₂ = −0.071, A ₁₆ = 0.042. By introducing the designed pupil filter into a conventional optical system with clear pupil, the improvement for imaging quality can be analyzed numerically through calculating the intensity distribution near the focal plane. The radial intensity distributions at the focal point are shown in Fig. 1 . The coordinates of the figure are in wavelengths. As Fig. 1 shows, the size of the central lobe is compressed by the designed pupil filter compared with the optical system of a clear pupil, which means there is resolution enhancement. The ability of extending DOF of the Zernike phase only pupil filter can be reflected on the axial intensity distribution across the focal region, as shown in Fig. 2 . Similarly, the coordinates are in wavelength. As can be seen from Fig. 2 , the length of the bright spot is much longer than that of a conventional optical system, which characterizes the DOF extension. Then, through comparison with the other two pupil filters, advantages of the Zernike phase only pupil filter are shown as follows. Here, forms of the two pupil filters are three zones [10] and Gaussian function [15] , their phase functions can be respectively written as: Three zones: Where ρ₁ and ρ₂ are the radiuses of the first and second rings, here, ρ₁ = 0.64, ρ₂ = 0.82 [10] . Gaussian functions: Where w is the width of half height, and n is order of Gaussian functions, here, w = 0.12, n = 8 [15] . The comparison of radial intensity distributions at the focal point are shown in Fig. 3 . It is clearly seen that the intensity distribution with the Zernike phase only pupil filter has the smallest radius of central lobe, which indicates its advantages on superresolution for a conventional optical system. Besides, the maximum intensity at the first side lobe of the proposed pupil filter is constrained quite low, which benefits from the 5^th, 12^th polynomials of the Zernike functions. The comparison of axial intensity distributions across the focal plane are shown in Fig. 4 . The value of the transverse dotted line equals 0.8, which is used to indicate effective defocus Strehl ratio. That means the imaging quality is acceptable at a certain defocus distance if the corresponding normalized intensity is greater than this value. It can be seen in Fig. 4 , the optical system with Zernike phase only pupil filter can keep its normalized intensity greater than 0.8 over the longest range of defocus distance. The summary of the ability of DOF extension and superresolution for these three pupil filters is shown in the table below. In Table 1 , effective depth of focal D_S.R. ≥ 0.8 of the clear pupil is normalized and superresolution factor G is inverted to show the resolution enhancement more intuitively. Compared with conventional optical system, DOF is enlarged 4.15 times and resolution is enhanced 1.32 times when introducing the Zernike phase only pupil filter. In order to demonstrate the ability of DOF extension and superresolution of our designed pupil filter further, imaging simulations are carried out. The virtual object used in the simulations is a resolution testing panel of line pairs with different spatial frequencies. Meanwhile, the defocus distance added in the simulations is from 0λ to 3λ. Results are shown in Table 2 . From the first line of these figures in Table 2 , the resolution of the optical system is clearly enhanced with all the three pupil filters compared to that with the clear pupil. In performance on DOF extension, the Zernike phase only pupil filter show its advantages. Its imaging results keep quite good quality over the entire defocus distance, which is very useful in numerous applications. In this paper, we have proposed an optical phase only pupil filter using Zernike functions, used to realize DOF extension and superresolution. Based on analyzing the relationship between Seidel aberrations and Zernike polynomials, design parameters of the proposed pupil filter are optimized. Compared with other two typical pupil filters, the Zernike pupil filter shows its advantages on both DOF extension and superresolution, enlarging DOF 4.15 times and enhancing resolution 1.32 times. Finally, imaging simulations with different defocus distances further demonstrate its effectiveness and superiority on improving the performance of the conventional optical system.