Tailoring the axial shape of the point spread function using the Toraldo concept

A novel procedure for shaping the axial component of the point spread function of nonparaxial focusing systems by use of phase-only pupil filters is presented. The procedure is based on the Toraldo technique for tailoring focused fields. The resulting pupil filters consist of a number of concentric annular zones with constant real transmittance. The number of zones and their widths can be adapted according to the shape requirements. Our method is applied to design filters that produce axial superresolution in confocal scanning systems. 2002 Optical Society of America OCIS codes: (100.6640) Superresolution; (110.1220) Apertures; (140.3300) Laser beam shaping; (170.1790) Confocal microscopy References and links 1. P. Jacquinot and B. Rozien-Dossier, “Apodisation” in Progress in Optics, E. Wolf, ed., Vol III (NorthHolland, Amsterdam, 1964). 2. A. Boivin, Théorie el Calcul des Figures de Diffraction de Révolution (Université de Laval, Quebec, 1964). 3. C. J. R. Sheppard and Z. S. Hegedus, "Axial behavior of pupil-plane filters," J. Opt. Soc. Am. A 5, 643-647 (1988). 4. M. Martínez-Corral, P. Andrés, J. Ojeda-Castañeda and G. Saavedra, "Tunable axial superresolution by annular binary filters. Application to confocal microscopy," Opt. Commun. 119, 491-498 (1995). 5. J. Campos, J. C. Escalera, C. J. R. Sheppard and M. J. Yzuel, "Axially invariant pupil filters," J. Mod. Opt. 47, 57-68 (2000). 6. C. J. R. Sheppard, M. D. Sharma and A. Arbouet, "Axial apodizing filters for confocal imaging," Optik 111, 347-354 (2000). 7. C. J. R. Sheppard, "Leaky annular pupils for improved axial imaging," Optik 99, 32-34 (1995). 8. M. Martínez-Corral, P. Andrés, C. J. Zapata-Rodríguez and M. Kowalczyk, "Three-dimensional superresolution by annular binary filters," Opt. Commun. 165, 267-278 (1999). 9. S. Grill and E. H. K. Stelzer, "Method to calculate lateral and axial gain factors of optical setups with a large solid angle," J. Opt. Soc. Am. A 16, 2658-2665 (1999). 10. M. A. A. Neil, R. Juskaitis, T. Wilson, Z. J. Laczik and V. Sarafis, "Optimized pupil-plane filters for confocal microscope point-spread function engineering," Opt. Lett. 25, 245-247 (2000). 11. M. Martínez-Corral, L. Muñoz-Escrivá, M. Kowalczyk and T. Cichocki, "One-dimensional iterative algorithm for three-dimensional point-spread function engineering," Opt. Lett. 26, 1861-1863 (2001). 12. G. Toraldo di Francia, "Nuovo pupille superresolventi," Atti Fond. Giorgio Ronchi 7, 366-372 (1952). 13. Z. S. Hegedus and V. Sarafis, "Superresolving filters in confocally scanned imaging systems," J. Opt. Soc. Am. A 3, 1892-1896 (1986). 14. T. R. M. Sales and G. M. Morris, "Diffractive superresolution elements," J. Opt. Soc. Am. A 14, 1637-1646 (1997). 15. T. R. M. Sales and G. M. Morris, "Axial superresolution with phase-only pupil filters," Opt. Commun. 156, 227-230 (1998). (C) 2002 OSA 14 January 2002 / Vol. 10, No. 1 / OPTICS EXPRESS 98 #39469 $15.00 US Received December 4, 2001; Revised January 9, 2002 16. C. J. R. Sheppard, G. Calvert and M. Wheatland, "Focal distribution for superresolving toraldo filters," J. Opt. Soc. Am. A 15, 849-856 (1998). 17. H. Wang and F. Gan, "High focal depth with a pure-phase apodizer," Appl. Opt. 40, 5658-562 (2001). 18. M. Gu, Advanced Optical Imaging Theory (Springer-Verlag, Berlin, 2000).


Introduction
Controlling the structure of the three-dimensional (3D) region that surrounds the focus of an optical system, i.e. the 3D point spread function (PSF), has long been the aim of many research efforts.In particular the design of techniques to overcome the limits in resolution imposed by diffraction in imaging systems is of great interest.Most of the reported apodization techniques have been aimed at improving the resolving capability of conventional twodimensional (2D) imaging system; and therefore improving the transverse resolution [1][2][3].However, when dealing with imaging systems designed to obtain the image of 3D samples, not only the transverse resolution but also the axial resolution is of great interest.In this context the design of pupil filters for obtaining axial superresolution was the purpose of many scientific contributions [3][4][5][6].Another purpose was tailoring the whole 3D PSF [7][8][9][10][11].Among all the techniques for designing pupil filters, one of the most innovative is the one reported by Toraldo di Francia [12].Toraldo showed that by subdividing the pupil into the proper number of concentric annular zones with constant complex transmittance, a band-limited transverse diffraction pattern of any shape may be obtained.The Toraldo method has been further applied for shaping the transverse PSF [13][14], and also the axial PSF [15][16][17].It should be noted however that in the later case the Toraldo method was not adapted to shape the axial PSF.On the contrary, the axial-shaping filters were obtained as a result of a comprehensive empirical study.
In this paper we present the axial form of Toraldo's pupil synthesis method.We show that by applying the Toraldo concept to a 1D properly nonlinearly mapped version of the pupil transmittance, one can control at will the positions of a finite number of zeros of the axial PSF.Moreover, since the axial PSF is a matter of interest in 3D imaging systems in which the lenses have high numerical aperture, such as confocal microscopy, we have extended the Toraldo method to the case of nonparaxially focused fields.

The axial response of apodized focusing systems
We start by considering the axial amplitude PSF of an aberration-free axially-symmetric focusing system in the nonparaxial scalar Debye approximation [18], that is where ( ) θ A represents the apodization function, i.e. the amplitude transmittance of the exit pupil, θ being the angular coordinate at the exit pupil plane as seeing from the focal point.z is the axial coordinate as measured from the focus of the system, α the semiaperture angle, λ the wavelength, and n the refractive index of the medium.
After performing the nonlinear mapping Eq. ( 1) can be rewritten as From Eq. ( 3) it follows that the axial amplitude PSF is governed by the 1D Fourier transform of a nonlinearly mapped version of the pupil amplitude transmittance.Note that this formula is valid to describe the axial behavior of both paraxial (low aperture) and nonparaxial (high aperture) focusing systems.To obtain a desired axial distribution it is necessary to solve an inverse problem and design a mapped function ( ) such that its 1D Fourier transform approximates the desired form.The optical implementation of a filter whose transmittance maps to such ( ) is in general a non-trivial task.An alternative approach is based on the Toraldo concept.As we show below it is possible to design ring-shaped pupil filters to control the position of the zeros of the axial-amplitude PSF.In the shaping procedure the influence of the angular aperture is two-fold.On one hand, for a given function ( ) the actual form of the apodization function explicitly depends on the value of α through the mapping of Eq. ( 2).On the other hand, the axial extent of the focal spot is determined by the value of α through the scale factor of Eq. ( 4).This means that, in general, the optimal pupil transmittance for a given purpose does not scale linearly with α , and then must be calculated explicitly for the NA value of interest.

The Toraldo concept
Toraldo di Francia [12] showed that the radii of the zero-intensity rings in the focal plane of a paraxial focusing system can be selected at will by using a pupil filters subdivided into concentric zones with constant transmittance.A thorough description of Toraldo's procedure can be found elsewhere [13,16].As remarked in Section 1, the aim of this work is to adapt the Toraldo concept to a new situation: the control of the shape of the axial PSF of nonparaxial focusing systems.We therefore outline the similarities and differences between the original Toraldo algorithm and our procedure.In the original algorithm the amplitude transmittance of the filter, ( ) r t , is subdivided into m concentric annular zones to control the radii of 1 − m rings of zero intensity.In our approach the function ( ) is also subdivided into constant-transmittance zones.Here it is important to take into account that if one wants to tailor the axial PSF with the constraint that the transverse PSF should remain almost invariant, the resulting function ( ) ζ q must be centrosymmetric [3,4].This implies that to control the positions of 1 − m axial zeros, the interval [-0.5, 0.5] should be divided into 1 2 − m subintervals such that in each subinterval the function ( ) ζ q is constant.On the basis of the above reasoning, the procedure for shaping the axial intensity of nonparaxially focused scalar fields is as follows: First, the function ( ) In Eq. ( 5) , and i k is the transmittance of the i th zone.According to Eq. ( 3) the axial amplitude PSF of this filter is given, apart from irrelevant constant factors, by the 1D Fourier transform of ( ) In the simplest case, corresponding to 2 = m , which allows one to control the position of one zero, the axial amplitude PSF is given by  (2).Note that the form of the pupil filter strongly depends on the value of α .

Application to Confocal Microscopy
The use of the three-zone Toraldo filter gives rise to an axial distribution in which the narrowness of the central lobe is accompanied by the enlargement of outer lobes.When dealing with some 3D imaging systems in which the axial resolution is of great importance, e.g.confocal fluorescence microscopes, this problem can be overcome.This is because the intensity PSF of the confocal system results from the product of the PSFs of the illuminating and the collecting systems.In a confocal fluorescence system the illumination of the sample and the collection of light are two statistically-independent processes.If one modifies the transmittance of the pupil of one of the arms, the corresponding PSF is shaped while the PSF of the other arm remains unaffected.What we propose here is to place the Toraldo filter in the illumination pupil and a circular aperture in the collection (to avoid light losses in the collection process).In this case, the resulting confocal-system PSF is narrowed in the axial direction, compared to the corresponding to the nonapodized confocal system, and the high axial sidelobes disappear.However, in fluorescence applications we should consider that the presence of high sidelobes in the illumination PSF could result in significant photo-bleaching of the 3D specimen at the positions of the strong secondary peaks.This problem can be overcome by designing filters with higher number of zones.Specifically we found that with a seven-zone Toraldo filter one can control the effect of photo-bleaching by sending the huge axial sidelobe further away.Note that in this case the design procedure is adapted to the thickness of the fluorescent sample, to ensure that the axial sidelobes remain outside the sample during imaging.In Fig. 3 we have plotted the mapped transmittance of the seven-zone filter, together with its normalized axial intensity PSF.In this case we see a significant reduction in the width of the central lobe, and samples up to Here also the PSF of the system results from the product between the PSFs of the illuminating and collecting processes.However, in this case the height of the axial sidelobes does not constitute a problem.The requirements for the design are: (a) to narrow the central lobe of the illumination axial PSF, in order to reduce the width of the central lobe of the confocal PSF; (b) to position the zeros of the illumination axial PSF so that they coincide with the maxima of the collection axial-PSF sidelobes.This minimizes the effect of the sidelobes in the confocal PSF, which was not the case in other approaches for obtaining axial superresolution [4,15].Taking into account these requirements we designed the seven-zone pupil filter shown in Fig. 4(a).In this figure we also show the contour plots of the 3D PSF for a bright-field confocal instrument with two circular pupils (Fig. 4(b)) and for the confocal system but with the Toraldo filter in the illumination set (Fig. 4(c)).
Note that the use of the seven-zone Toraldo filter narrows the intensity PSF, alog the optical axis, and therefore provides an important improvement of the optical sectioning capacity of the microscope.

Conclusions
The axial form of the Toraldo pupil synthesis method has been presented.The reported method allows the design of filters with the ability to control the positions of the zeros of the axial PSF of a lens, while preserving the structure of the transverse PSF.We show that, given a desired shape for the axial PSF, the actual form of the Toraldo filter explicitly depends on the aperture of the lens.The resulting filters, whose main feature is their simplicity, can be applied to improve the performance of 3D image formation systems, such as fluorescence and brightfield reflection confocal microscopes.
imaged without bleaching by the sidelobes.Of course, if one wants to send the (C) 2002 OSA 14 January 2002 / Vol. 10, No. 1 / OPTICS EXPRESS 102 huge axial sidelobe further away it is necessary to design a Toraldo filter composed by a higher number of zones.A practical rule for the design procedure is that the position of the huge sidelobe is approximately 1 + = m z N .Another application of Toraldo filters is in reflection bright-field confocal microscopy.

Fig. 4 .
Fig. 4. (a) Seven-zone Toraldo filter designed to obtain axial superresolution in reflection confocal microscopy.The filter was calculated for the case of °= α 5 67.; (b) 3D intensity PSF of a reflection confocal microscope with two circular pupils; (c) as b) but with the seven-zone Toraldo filter in the illumination system.