Reduction of spherical-aberration impact in microscopy by wavefront coding

In modern high-NA optical scanning instruments, like scanning microscopes, the refractive-index mismatch between the sample and the immersion medium introduces a significant amount of spherical aberration when imaging deep inside the specimen, spreading out the impulse response. Since such aberration depends on the focalization depth, it is not possible to achieve a static global compensation for the whole 3D sample in scanning microscopy. Therefore a depth-variant impulse response is generated. Consequently, the design of pupil elements that increase the tolerance to this aberration is of great interest. In this paper we report a hybrid technique that provides a focal spot that remains almost invariant in the depth-scanning processing of thick samples. This invariance allows the application of 3D deconvolution techniques to that provide an improved recovery of the specimen structure when imaging thick samples. © 2009 Optical Society of America OCIS codes: (110.6880) Three-dimensional image acquisition; (100.2000) Digital image processing; (100.6890) Three-dimensional image processing; (180.6900) Three-dimensional microscopy. References and Links 1. J. Tsujiuchi, “Correction of optical images by compensation of aberrations and by spatial frequency filtering,” Prog. Opt. 2, 131–180 (1963). 2. R. Barakat, and A. Houston, “Transfer function of an annular aperture in the presence of spherical aberration,” J. Opt. Soc. Am. A 55(5), 538–541 (1965). 3. J. P. Mills, and B. J. Thompson, “Effect of aberrations and apodization on the performance of coherent optical systems. I. The amplitude impulse response,” J. Opt. Soc. Am. A 3(5), 694–703 (1986). 4. J. Ojeda-Castañeda, P. Andrés, and A. Díaz, “Annular apodizers for low sensitivity to defocus and to spherical aberration,” Opt. Lett. 11(8), 487–489 (1986). 5. J. Ojeda-Castañeda, P. Andrés, and A. Díaz, “Strehl ratio with low sensitivity to spherical aberration,” J. Opt. Soc. Am. A 5(8), 1233–1236 (1988). 6. S. Mezouari, and A. R. Harvey, “Phase pupil functions for reduction of defocus and spherical aberrations,” Opt. Lett. 28(10), 771–773 (2003). 7. C. J. R. Sheppard, and M. Gu, “Aberration compensation in confocal microscopy,” Appl. Opt. 30(25), 3563– 3568 (1991). 8. M. J. Booth, M. A. A. Neil, R. Juskaitis, and T. Wilson, “Adaptive aberration correction in a confocal microscope,” Proc. Natl. Acad. Sci. U.S.A. 99(9), 5788–5792 (2002). 9. M. J. Booth, M. Schwertner, T. Wilson, M. Nakano, Y. Kawata, M. Nakabayashi, and S. Miyata, “Predictive aberration correction for multilayer optical data storage,” Appl. Phys. Lett. 88(3), 031109 (2006). 10. C. Mauclair, A. Mermillod-Blondin, N. Huot, E. Audouard, and R. Stoian, “Ultrafast laser writing of homogeneous longitudinal waveguides in glasses using dynamic wavefront correction,” Opt. Express 16(8), 5481–5492 (2008). 11. C. Preza, and J. A. Conchello, “Depth-variant maximum likelihood restoration for three-dimensional fluorescence microscopy,” J. Opt. Soc. Am. A 21(9), 1593–1601 (2004). 12. J. Ojeda-Castañeda, R. Ramos, and A. Noyola-Isgleas, “High focal depth by apodization and digital restoration,” Appl. Opt. 27(12), 2583–2586 (1988). 13. E. R. Dowski, Jr., and W. T. Cathey, “Extended depth of field through wave-front coding,” Appl. Opt. 34(11), 1859–1865 (1995). 14. W. T. Cathey, and E. R. Dowski, “New paradigm for imaging systems,” Appl. Opt. 41(29), 6080–6092 (2002). #110948 $15.00 USD Received 4 May 2009; revised 23 Jun 2009; accepted 24 Jun 2009; published 24 Jul 2009 (C) 2009 OSA 3 August 2009 / Vol. 17, No. 16 / OPTICS EXPRESS 13810 15. M. Gu, Advanced optical imaging theory, Springer-Verlag (Berlin) 2000. 16. In a real microscopy experiments there is, still, a coverslip between the immersion liquid and the sample medium. Such slab may induce a constant amount of SA that does not depend on the focusing depth. This SA can be compensated statically either with the correction collar, or with a proper index-matching immersion medium. Therefore it is not necessary to consider it in the calculations. 17. The highly corrected objectives are equipped with a correction collar to allow adjustment of the central lens group position to coincide with fluctuations in cover-glass thickness. 18. C. J. R. Sheppard, and C. J. Cogswell, “Effects of aberrating layers and the tube length on confocal imaging properties,” Optik (Stuttg.) 87, 34–38 (1991). 19. P. Török, P. Varga, Z. Laczik, and G. R. Booker, “Electromagnetic diffraction of light focused through a planar interface between materials of mismatched refractive indices: an integral representation,” J. Opt. Soc. Am. A 12(2), 325–332 (1995). 20. V. N. Mahajan, “Gaussian apodization and beam propagation,” Prog. Opt. 49, 1–95 (2006). 21. I. Escobar, G. Saavedra, M. Martínez-Corral, and J. Lancis, “Reduction of the spherical aberration effect in highnumerical-aperture optical scanning instruments,” J. Opt. Soc. Am. A 23(12), 3150–3155 (2006). 22. J. Ojeda-Castaneda, J. E. A. Landgrave, and H. M. Escamilla, “Annular phase-only mask for high focal depth,” Opt. Lett. 30(13), 1647–1649 (2005). 23. D. D. Lowenthal, “Maréchal intensity criteria modified for Gaussian beams,” Appl. Opt. 13(9), 2126–2133 (1974). 24. M. Born, and E. Wolf, Principles of Optics. University Press (Cambrigde) 1999. Ch 13. 25. C. J. Cogswell, S. V. King, S. R. Prasanna Pavani, D. B. Conkey, and R. H. Cormack, “Microscope imaging capabilities improve using computational optics,” Proc. Focus on Microscopy 09 (http://focusonmicroscopy.org/2009/PDF/341_Cogswell.pdf) 26. S. V. King, A. R. Libertun, R. Piestun, C. J. Cogswell, and C. Preza, “Quantitative phase microscopy through differential interference imaging,” J. Biomed. Opt. 13(2), 024020 (2008).


Introduction
The acquisition of high-resolution images of 3D biological samples by optical microscopes requires the production of tightly focused spots, which are usually achieved by high numerical aperture (NA) immersion objectives.The refractive-index mismatch between the immersion liquid and the solution in which the sample is embedded introduces a significant amount of spherical aberration (SA) when imaging deep inside the specimen, spreading out the impulse response.Since such aberration depends on the focusing depth, the impulse response is degraded in a different manner for different slices inside the specimen.Thus, it is not possible to achieve a global compensation for the whole sample.Some methods were developed, in paraxial context, for the reduction of SA impact [1][2][3][4][5][6].These proposals are not directly applicable to optical microscopy since they do not take into account the tight focusing.More recent new proposals, based on non-paraxial equations, have been reported.The idea of altering the objective tubelength [7], or the use of adaptive optics [8][9][10], are very interesting.Another interesting approach is the use of a depth-variant stratum-based algorithm [11].This interesting approach allows the acquisition of aberration-free images taking advantage of the fact that it is possible to compute the deconvolution by using a small number of spatially invariant PSFs, each computed at different depth.
Here we focus on a technique that is very simple, but that is able to stabilize the SA impact over a large range of imaging depths.The technique is known as wavefront coding.This technique was pioneeringly suggested, in a different context, by Ojeda-Castañeda [12] and mainly developed by Cathey's group [13,14].Wavefront coding was developed for increasing the depth-of-field of imaging systems.A phase pupil mask is used to encode the wavefront emerging from the imaging system while the digital restoration produces the decoding.In this paper we take advantage of the fact that the axial values of the intensity spot obtained after focusing, through a plane interface, the light proceeding from a high-NA objective is proportional to the Fresnel transform of a properly mapped version of the objective aperture stop.This fact allows us to adapt the wavefront coding technique to the aim of reducing the impact of the SA induced during the axial scanning of the sample.The hybrid technique we report here provides a focal spot that remains almost invariant in the depth-scanning processing of thick samples.This invariance allows the application of 3D deconvolution techniques with the purpose of obtaining an improved recovery of the specimen structure when imaging thick samples.

Spherical aberration induced in the axial scanning
As is well known, the imaging features of an optical microscope can be accurately described through the 3D intensity point-spread-function (PSF), namely, by the 3D intensity distribution that is obtained in the neighborhood of the focal point of the microscope objective when it is illuminated by a monochromatic plane wave.Since, to obtain optical images with good resolution, the microscope objectives need to have high NA, the PSF cannot be calculated by means of paraxial diffraction formulas, but through the non-paraxial Debye's equation.According to the scalar Debye's formulation, and assuming that sine condition and axial symmetry hold, the intensity PSF is [15].
In this equation, ( ) P θ accounts for the amplitude transmittance at the objective exit pupil, and α is the maximum value for the aperture angle θ .Lateral and axial positions in the focal region are expressed through cylindrical coordinates ( ) , r z as measured from the Gaussian focus, F .The phase factor ( ) W θ accounts for potential phase distortions occurred during the focusing.
In many microscope observations a significant mismatch between refractive indices of the sample medium and the immersion liquid occurs.This mismatch induces a significant amount of SA.To analyze this effect, we consider the simple geometry in which the spherical wave emerging from the microscope objective is focused deeply through a planar interface between two media of different refractive indices (see Fig. 1) [16].To evaluate the phase distortions, we take into account that, according to Debye's formulation, the spherical wave can be seen as the superposition of plane waves whose directions fall inside the cone defined by the focus,  In Fig. 1 we have represented a plane-wave component of the field by a light-ray normal to the wavefront.We have named as the focusing depth, 0 d , the distance between the interface and the focus, F. On the contrary, we have named as the penetration depth, 0 ' d , the distance between the interface and the place where the wave effectively focuses.
In a typical microscopy procedure, the image of the different sections of the sample are obtained after axial scanning of the microscope.If at a given stage of the axial scanning the microscope is displaced a distance S z towards the sample, the focusing depth is S 0 S In The above equation has been obtained after expanding ( ) S ; W z θ into a power series and making a fourth-order approximation.Besides, we have introduced the well-known defocus and SA coefficients, which are given, respectively, by In Fig. 2 we show the results of a numerical simulation of the 3D intensity distribution in the focal region of a microscope objective when typical axial-scanning induced SA occurs.In our calculations we have considered a light microscope equipped with an immersion objective of NA = 1.4 (oil; 1 1.52 n = ), and a specimen whose refractive index is close to 2 1.33 n = , which is the index of the water solution where the specimen is embedded.We have assumed that the SA has been compensated, with a correction collar, for a penetration depth 0 '   In the second step of the simulation we have calculated the images of a spoke target set perpendicular to the optical axis of the microscope at the penetrations depths used in Fig. 2 3.A significant feature comes out from this figure: there is only a very small degradation of images of 2D samples.This is because, although the 3D PSF is strongly spread, the lateral PSF at the best image plane still has a sharp peak even in case of significant amount of SA.

Effect of the spherical aberration in 3D scanning microscopy
We have just seen that the SA does not degrade the images of 2D objects significantly, provided that one selects the penetration depth adequately.However we have seen that SA strongly degrades the 3D PSF and therefore we infer that images of 3D objects will be degraded as well.To analyze this fact in a simple way it is convenient to study the influence of the SA in the axial PSF.To this end, we perform some mathematical manipulations of Eq.
(1).Then by defining the normalized coordinate ( ) and after performing the nonlinear mapping ( ) ( ) we straightforwardly obtain the axial response of the system as [21] ( ) Note that Eq. ( 7) establishes a Fresnel-transform relationship between the axial response of the system and the mapped pupil ( ) q ζ .The same result can be encountered in the analysis of 1D paraxial focusing systems when considering the case of a cylindrical lens illuminated by a monochromatic plane wave.When the 1D aperture stop (of transmittance ( ) 0 p x ) is placed, as usually, at the front focal plane of the cylindrical lens of focal length x f , the intensity distribution on planes parallel to the back focal plane is given by ( ) ( ) ( ) where we have introduced the well-known 1D defocus coefficient, defined as ∆ being the width of the slit.As pointed out in [22] in a paraxial context, it is apparent that the 1D defocus coefficient 20 x w plays in Eq. ( 8) a role similar to the one played by coefficient 40 w in Eq. ( 7).Thus, one can conclude that the axial intensity distribution produced by a high-NA scanning optical microscope in which certain amount of SA is induced, is the same as the transverse intensity distribution obtained at a defocused plane experiment.Since very efficient solutions have been given for decreasing the effect of defocusing in the 1D case (i.e., for increasing the depth-of-field of such systems), we can now apply the analogous techniques to efficiently desensitize the axial response of the objective to the sample-induced SA.

Application of the wavefront coding technique
A well-known method for extending the depth-of-field of imaging systems is wavefront coding (WFC).In the WFC technique a pupil mask designed for increasing the depth-of-field is used in combination with a digital image-restoration process.Thus, as a first stage in the WFC method, the defocused patterns are uniformly blurred over a large axial range by effect of the mask.The almost-constant transverse impulse response in this case allows a sharp reconstruction inside this extended range by a single step deconvolution.A successful proposal in this sense is the 1D cubic phase filter proposed by Dowski and Cathey [13].To adapt this technique for the problem of reducing the SA impact, we propose the use of a properly scaled radial version of the cubic phase filter, namely ( ) ( ) Following the suggestions made in [13], we have set a value for A much higher than 5. Thus we set A = 50.Since SA is of particular impact only in 3D imaging experiments, we have performed a simulation in which we consider a confocal scanning microscope.To perform the calculations we assume, again, that the spherical aberration is compensated, by the correction collar, for a penetration depth o ' 10 d m µ = .We calculated, as the squared modulus of Eq. ( 7), the confocal PSFs corresponding to seven equally-spaced scanning steps (see Fig. 4).In the video we compare, for varying penetration depth, the focal spots produced by the circular aperture (CA) and by the cubic filter (CF).We see that the PSFs obtained with the cubic filter are very different from the ones obtained with the circular aperture.In fact they are much more spread.But, in contrast with the circular-aperture case, they do not change significantly as the penetration depth increases.So, we can affirm that the cubic phase filter produces a spread focus that is insensitive to the amount of SA.In Table 1 we compare the values of S ' d and S d , and show the corresponding value of the SA coefficient.Note that there is not any proportional relation between the penetration depth and the induced spherical aberration [23].Any 2D section of the 3D image, for example the one corresponding to a given scanning depth S ' z , is calculated as the incoherent superposition of eight 2D images.Each of these eight images is obtained as the 2D convolution between the corresponding plate and the defocused transverse section of the 3D PSF associated to the selected S ' z .In Fig. 6(a) we show five lateral sections of the 3D image obtained with the circular aperture as the objective aperture stop.We clearly see that the deeper the scanning position the bigger the induced SA and therefore the higher the influence of adjacent plates in the 2D sections of the 3D image.We have applied a typical inverse filtering procedure to the 3D confocal image obtained with the circular aperture.We have used for the 3D deconvolution algorithm the PSF corresponding to an intermediate depth, S ' 10.00 z m µ = .In Fig. 6(b) we show five sections of the deconvolved 3D image.These images are, or course much worse than the original.This is because it is in fact nonsense to apply deconvolution to systems whose impulse response varies very fast.
In Fig. 6(c) we show the sections of the 3D confocal image but obtained by placing the cubic phase filter in the aperture stop of the objective.As expected, the images are, mainly for low values of the SA, much worse than the ones obtained with the circular aperture.However, we can infer from the figure that the influence of adjacent plates keeps fairly constant.Note that due to the non-symmetrical form of the PSFs, each section of the 3D image is influenced mainly by less deep plates.
Since the PSF obtained with the cubic-phase filter remains practically invariant over a wide range of focusing depths, the microscope in which the filter is inserted can be considered, in quite good approximation, as 3D linear and shift invariant, and therefore a good candidate for the application of deconvolution techniques.Thus we have applied a typical inverse filtering procedure to the 3D confocal image obtained with the cubic filter.As the 3D PSF of the system we have selected the 3D PSF corresponding to an intermediate focusing depth, namely S ' 10.00 z m µ = .After performing the 3D deconvolution we have obtained a new 3D image of the synthetic object.In Fig. 6(d) we show five 2D sections of the said object.We conclude from the figure that a wavefront coding technique provides 3D images that are highly immune to the SA induced by axial scanning.

Conclusions
We have presented a hybrid technique for the reduction of impact of spherical aberration induced in the axial scanning inherent to microscopy observations.Note that, in general our design model should have taken into account that the radiation emitted by distant object points might be scattered or reflected by nearer object points.However, we have assumed that the first-order Born approximation holds, so that the multiple scattering and depletion of the incident beam are negligible [24].The main feature of the proposed method is that with a minimum modification of the objective architecture one can mitigate the SA induced during the axial scanning inherent to any 3D imaging procedure.We have illustrated the utility of our method by means of a numerically simulated 3D microscopy experiment.A question that remains open is how to reduce the noise artifacts inherent to the use of wavefront coding technique.In this sense, in further research we will consider the use of non-linear deconvolution procedures, as suggested in [25,26].
F, and the exit pupil.Each plane-wave component of the field emerging from the objective obeys Snell's law, at the interface.The resulting field is reconstructed as the superposition of refracted plane waves.

Fig. 1 .
Fig. 1.Scheme for the evaluation of the phase distortions occurred in the tight focusing through a refractive-index interface.Definition of nomenclature: focusing depth, 0 d , and penetration

Fig. 2 .
Fig. 2. Intensity distribution in the focal spot calculated for the axial displacements: S ' 0 z = , . The images were obtained by convolving the target with the lateral response at the best image plane (marked with a yellow line in the figure).Note that in contrast with what is commonly assumed [20], we have not set the best focal plane at 20 40 w w = − , but at the axial position of the PSF main peak.The images are shown in Fig.

Fig. 4 .
Fig. 4. Axial intensity distribution in the focal spot calculated for S ' 0 z = and other six equally-spaced scanning positions.In the video (Media 1) we compare, for varying penetration depth, the focal spots produced by the circular aperture (CA) and by the cubic filter (CF).
To illustrate the utility of the proposed wavefront coding method we have performed a numerical simulation of a confocal imaging experiment.As the 3D object we have used a synthetic object composed by 8 plane plates set at S ' are shown in Fig.5.

Fig. 5 .
Fig. 5.The synthetic 3D object used for the numerical experiment is composed by 8 equidistant plates.

Fig. 6 .
Fig. 6.Five sections of the 3D image obtained with: (a) The clear circular pupil as the aperture stop; (b) The clear circular pupil and 3D deconvolution; (c) The cubic phase filter as the aperture stop; (d) The proposed hybrid method in which the cubic filter is inserted and a 3D deconvolution algorithm is applied.
such case the penetration depth is S