Optimization of Wavefront Coding Phase Mask Applied to 5X40X MicroObjectives Simultaneously
 Author: Liu Jiang, Miao Erlong, Sui Yongxin, Yang Jianghuai
 Publish: Journal of the Optical Society of Korea Volume 19, Issue5, p487~493, 25 Oct 2015

ABSTRACT
A wavefront coding (WFC) technique provides an extension of the depth of field for a microscopy imaging system with slight loss of image spatial resolution. Through the analysis of the relationship between the incidence angle of light at the phase mask and the system pupil function, a mixing symmetrical cubic phase mask (CPM) applied to 5X40X microobjectives is optimized simultaneously based on pointspread function (PSF) invariance and nonzero mean values of the modulation transfer function (MTF) near the spatial cutoff frequency. Optimization results of the CPM show that the depth of field of these microobjectives is extended 310 times respectively while keeping their resolution. Further imaging simulations also prove its ability in enhancing the defocus imaging.

KEYWORD
Wavefront coding , Microobjective , Depth of field , Imaging system

I. INTRODUCTION
A microobjective plays a decisive role in microscopy imaging systems, determining the resolution and magnification. In order to achieve better image quality and high resolution, a microobjective is designed with relatively large aperture, which means the depth of field is small.
Wavefront coding (WFC) is an opticaldigital hybrid imaging technique developed to extend the depth of field of conventional optical systems [13]. By introducing a phase mask into the pupil of the system to modulate the wavefront of the incident light, it is possible to make the pointspread function (PSF) nearly invariant over a wide range of defocus. Then sharp images can be restored from both focus and defocus images by the same digital deconvolution filter, extending the depth of field.
In the past several years, many different kinds of WFC phase mask have been designed [4, 5]. Because of relative simplicity in manufacture, the cubic phase mask is the most widely used form in the WFC technique [68]. Various optimized methods have been used to improve the performance of the cubic phase mask (CPM) [912]. Zhang and Chen optimized the CPM based on Strehl ratio. Carles presented an analytical approach to determine the optimal CPM strength, and Liu made the stationary phase analysis of CPM. However, most optimized CPMs in microscopy systems are only applied to one power microobjective [13]. That means microobjectives with different powers each need a specific CPM, which is not practical and economical.
In this paper, we optimize a mixing symmetrical CPM. Through the analysis of the relationship between the incidence angle of light at the phase mask and system pupil function, we build a merit function (MF) based on the PSF invariance and nonzero mean values of the modulation transfer function (MTF) near the spatial cutoff frequency to optimize the CPM. Then, we analyze the ability of optimized CPM in extending the depth of field for different power microobjectives and its effect on spatial resolution. Finally, we show its enhancement on imaging of 5X40X microobjectives with different defocus distances.
II. THE THEORETICAL ANALYSIS OF WAVEFRONT CODING PHASE MASK OF MICROOBJECTIVES WITH DIFFERENT POWERS
A schematic of a wavefront coding microscopy imaging system is shown in Fig. 1.
A WFC phase mask is located at the system pupil plane after the microobjective. CCD detects an encoded image with the same blur degree at different defocus distances due to the modulation of the phase mask. Then, a sharp image will be restored by the same deconvolution filter, which is built through the system PSF.
Numerical aperture (NA) of a microobjective is closely related to the maximum incidence angle of light in object space, which determines its resolution. That means microobjectives with different powers are not the same. When introducing a WFC phase mask in a microscopy system, the incidence angle of light at the phase mask also changes following the converting of microobjectives. It is shown in Fig. 2. In order to extend the depth of field of microobjectives with different powers using one phase mask, the influence of different incidence angle of light at the phase mask is analyzed as follows.
We select a mixing symmetrical CPM due to its manufacturing convenience. Its surface shape is represented by
Where
x andy are the normalized coordinates,α andβ are the coefficients of the CPM.For an aberrationfree optical system with only defocus, its pupil function modulated by CPM is represented by
Where,
ϕ is the phase of the pupil wavefront,W _{020} is the coefficient of defocus andλ is the wavelength. If 2π /λ is a constant, we analyze that the pupil function is a simple form by ignoring that constant and representing the pupil function byThe difference of incidence angle at CPM could be regarded as a rotation of coordinates, shown in Fig. 3, and transformation formulas between two coordinates can be represented by
Where,
x ,y ,z are the original coordinates,x ′,y ′,z ′ are the transformed coordinates,θ is the difference of two different incidence angles at CPM.Substituting Eq. (4) into Eq. (1) and simplifying, the pupil function with different incidence angles is represented by
The coefficients of CPM in pupil function vary with the incidence angle. It is shown in Fig. 4. Coefficient
α is getting larger with incidence angle increasing both at meridian and sagittal directions. It grows faster at the meridian direction. But coefficientβ remains unchanged at the sagittal direction and is getting smaller at the meridian direction. Two conclusions can be drawn from this, 1) selecting suitableα andβ can restrain the influence of incidence angle; 2) there is a difference between meridian and sagittal directions.III. THE OPTIMIZATION OF WAVEFRONT CODING PHASE MASK OF MICROOBJECTIVES WITH DIFFERENT POWERS
Based on the invariance of PSF, the merit function (MF) for WFC phase mask is represented by
Where,
Std (PSF_{diffdefocus} )_{i,j} is the standard deviation of PSF at the same spatial location of microobjectives with different powers, defined to represent the PSF invariance,x andy are sampling numbers of meridian and sagittal directions, N is the sum of sampling numbers.Variation of coefficients of phase mask will introduce more aberration in the microscopy imaging system, which results in decrease of value of MTF near spatial cutoff frequency. That creates a decline of resolution of the optical system.
Thus, optimization of our CPM adds a punishment function (PF) in order to eliminate or decrease this impact, and taking the difference between meridian and sagittal directions into consideration. The PF is represented by
Where,
Mean (MTFM_{freq} ) andMean (MTFM_{freq} ) are mean values of MTF near the spatial cutoff frequency in meridian and sagittal directions, N is the number of sampling frequencies.Then, the merit function for CPM is represented by
Finally, through an optimization search algorithm, such as a genetic algorithm (GA) or simulated annealing algorithm (SA), we could solve for the optimal coefficients of the CPM.
IV. OPTIMIZATION RESULT AND DISCUSSION
Shape of optimized CPM is shown in Fig. 5.
We combine this phase mask with 5X, 10X, 20X and 40X microobjectives, computing their MTF in optical design software (ZEMAX). These results are shown in Fig. 6. The MTFs of all microobjectives with the phase mask have an acceptable value near the cutoff frequency, by which we mean not losing any high spatial frequency, maintaining the resolution of microscopy. In addition, the higher power the microobjective is, the larger is the extension of the ranges of depth of field with the optimized CPM. For 5X objective, its depth of field is extended from 75 µm (without CPM) to about 200 µm (with CPM). It can be seen from Fig. 6(c), the value of MTF for the optical system with CPM at defocus distance 200 μm is already very small. That demonstrates the maximal ability of the CPM for extending the depth of field, because the value of MTF will drop to zero probably when defocus distance is over 200 μm for a practical optical system with aberrations. Similarly, the optimized CPM multiplies the depth of field of 10X, 20X, 40X objectives five to ten times.
Since the microobjectives with optimized CPM maintain their proper spatial cutoff frequency, the shape images can be restored from the encoded images. Below we present the comparison of microscopy imaging simulation results of microobjectives with different powers at different defocus locations. These results are exhibited in Tables 14. The target image is a spoke pattern with numbers and marks at corners of which the size is 512 × 512 pixels. The coordinates in all figures in Tables are in pixels.
V. CONCLUSION
Based on the theoretical analysis of the relationship between the incidence angle of light and pupil function, a mixing symmetrical cubic wavefront phase mask has been optimized. The depth of field of 5X40X microobjectives optical system have multiplied 3 to 10 times larger after introducing optimized CPM. In addition, the results of microscopy imaging simulation with different defocus distance show its ability in enhancement for defocus imaging.

[FIG. 1.] Schematic of a wavefront coding microscopy system.

[FIG. 2.] The incidence angle at the phase mask in the optical path of microobjectives with different powers.

[]

[]

[]

[FIG. 3.] The incidence angle at phase mask and its coordinate transform.

[]

[]

[FIG. 4.] The incidence angle at phase mask and its coordinates transform.

[]

[]

[]

[FIG. 5.] Shape of WFC phase mask.

[FIG. 6.] The MTFs of 5X, 10X, 20X, 40X microobjectives without (a), (d), (g), (j) and with (b), (e), (h), (k) the CPM. (c), (f), (i), (l) are the partial enlarged views of (b), (e), (h), (k) near the spatial cutoff frequency.

[TABLE 1.] Imaging simulation results of 5X microobjectives

[TABLE 2.] Imaging simulation results of 10X microobjectives

[TABLE 3.] Imaging simulation results of 20X microobjectives

[TABLE 4.] Imaging simulation results of 40X microobjectives