Advances in Optical TechnologiesVolume 2013 (2013), Article ID 360450, 6 pageshttp://dx.doi.org/10.1155/2013/360450

Research Article

Far-field Diffraction Properties of Annular Walsh Filters

Pubali Mukherjee¹ and Lakshminarayan Hazra²

¹MCKV Institute of Engineering, 243 G.T. Road, Liluah, Howrah 711204, India²Department of Applied Optics and Photonics, University of Calcutta, 92 A.P.C. Road, Kolkata 700 009, India

Received 22 March 2013; Revised 3 June 2013; Accepted 2 July 2013

Academic Editor: Augusto Belendez

Copyright © 2013 Pubali Mukherjee and Lakshminarayan Hazra. 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.

Abstract.

Annular Walsh filters are derived from the rotationally symmetric annular Walsh functions which form a complete set of orthogonal functions that take on values either +1 or −1 over the domain specified by the inner and outer radii of the annulus. The value of any annular Walsh function is taken as zero from the centre of the circular aperture to the inner radius of the annulus. The three values 0, +1, and −1 in an annular Walsh function can be realized in a corresponding annular Walsh filter by using transmission values of zero amplitude (i.e., an obscuration), unity amplitude and zero phase, and unity amplitude and phase, respectively. Not only the order of the Walsh filter but also the size of the inner radius of the annulus provides an additional degree of freedom in tailoring of point spread function by using these filters for pupil plane filtering in imaging systems. In this report, we present the far-field amplitude characteristics of some of these filters to underscore their potential for effective use in several demanding applications like high-resolution microscopy, optical data storage, microlithography, optical encryption, and optical micromanipulation.

1. Introduction

Annular apertures and different types of ring-shaped apertures continue to be investigated for catering to the growing exigencies in diverse applications, for example, high resolution microscopy, optical data storage, microlithography, optical encryption, and optical micromanipulation [1–5]. Not only for obvious energy considerations but also for their higher inherent potential in delivering complex far-field amplitude distributions, annular phase filters are being investigated in different contexts [6–9].

A systematic study on the use of phase filters on annular pupils can be conveniently carried out with the help of annular Walsh filters derived from the annular Walsh functions. Walsh functions form a closed set of normal orthogonal functions over a given finite interval and take on values +1 or −1, except at a finite number of points of discontinuity, where they take the value zero [10, 11]. The order of a Walsh function is directly related to the number of its zero crossings or phase transitions within the specified domain, and they constitute a closed set of normal orthogonal functions over the specified interval. They have the interesting property that an approximation of a continuous function over a finite interval by a finite set of Walsh functions leads to a piecewise constant approximation to the function. Walsh filters of various orders may be obtained from corresponding Walsh functions, by realizing transmission values of +1 and −1 by 0 and phase, respectively. Incidentally it may be noted that binary phase filters are being explored for many interesting light distributions [12].

Walsh functions have been used in the field of signal coding and transmission and in allied problems of information processing [13]. Two-dimensional Walsh functions in the usual rectangular coordinates have been used in digital image processing applications [14]. For treatment of problems of optical imaging, Walsh functions in polar co-ordinates have been utilised [15]. For systems with rotational symmetry about the axis, radial Walsh functions [16] have been developed as a special case of Walsh functions in polar co-ordinates and they were proved useful in the treatment of apodization problems [17, 18]. It has also been observed that not only the transverse amplitude distribution on the far-field plane but also the axial distribution of amplitude/intensity in the far-field is significantly modified in presence of radial Walsh filters on unobscured apertures [19].

Annular Walsh functions are a generalization of radial Walsh functions. For a specific central obscuration ratio, they are a complete set of orthogonal functions over the annulus. Annular Walsh filters may be considered to have ternary transmission values. The three values of transmission are zero amplitude over the central obscuration, unity amplitude and zero phase for value +1, and unity amplitude and phase for value −1 in corresponding annular Walsh functions.

In this paper, we present some results of our investigations on the far-field amplitude characteristics of annular Walsh filters of orders 0, 1, 2, and 3 with different central obscuration ratios. After a brief description of annular Walsh functions in the next section, Section 3 presents the mathematical expression for the amplitude distribution in the far-field of an exit pupil with annular Walsh filters on it. Some interesting numerical results and our observations on the same are put forward in the last two sections.

2. Annular Walsh Functions

To define annular Walsh function of index and argument over an annular region with and 1 as inner and outer radii, respectively, it is necessary to express the integer in the form where are the bits, 0 or 1 of the binary numeral for , and is the power of 2 that just exceeds . For all in , is defined as where The orthogonality condition implies that where is the Kronecker delta defined as Figure 1 shows the first four annular Walsh functions for central obscuration in two dimensions. Figure 2 presents values of the functions , , along the radius in an azimuthal direction. It should be noted that the order of the functions is equal to the number of zero crossings, or sign changes of the function in the interval , and locations of the points of zero crossings for members of the set of functions , are given by In general, for the first , where and is a positive integer, and annular Walsh functions , , the zero crossings are located at The inner and outer radii of the annulus is and . Note that the set of zero crossing locations, , , consists of all zero crossing locations required for specifying members of this particular set of Walsh functions. An individual member of this set of Walsh functions will have the same number of zero crossings as its order.

Figure 1: Annular Walsh functions in two dimensions of order for central obscuration .

Figure 2: Annular Walsh functions of order along radius for central obscuration .

For computational purposes it is often convenient to express an annular Walsh function as where . For a particular value of , the integer is taken such that just exceeds . are zero-one functions, also known as Walsh block functions, defined as

For all annular Walsh functions of order , the same values of , are to be used. They are unique for the particular value of , and their values are given in (7). are the elements of a Hadamard matrix whose elements are +1 or −1 [13]. Table 1 gives the Hadamard matrix for . The Hadamard matrices for other values of can be obtained from the defining relation (2) of annular Walsh functions.

Table 1: The Hadamard matrix for .

+1
+1
+1
+1

+1
+1
−1
−1

+1
−1
−1
+1

+1
−1
+1
−1

3. Far-Field Amplitude Distribution

Figure 3 shows the image space of an axially symmetric imaging system. It is well known that the complex amplitude distribution on the image plane corresponding to an axial object point is given by the far-field diffraction pattern of the pupil function over the exit pupil. The complex amplitude at a point in the far-field due to a point object on the axis is, apart from the multiplicative constant, given by the Hankel transform of order zero of the pupil function [20, 21]: where is the circularly symmetric pupil function. The variable is the fractional coordinate for a point on the pupil sphere of center and radius , where is the axial position of the exit pupil. is equal to the ratio , where is the perpendicular distance of the point from the axis and is the maximum value of , so that . is the reduced diffraction variable given by where is the semiangular aperture of the imaging system, is the refractive index of the image space, is the propagation constant, and is the actual geometrical distance of point on the image plane from the center of the diffraction pattern.

Figure 3: Schematic representation of the image space parameters.

For a specific central obscuration ratio , annular Walsh filters of various orders can be obtained from the corresponding annular Walsh functions by realizing transmission values +1 and −1 by zero and phase filters, respectively. In presence of an annular Walsh filter on a uniform pupil, the pupil function is given by It may be noted that is binary (value either 0 or +1) only in the case of zero order annular Walsh function ; for all other orders, is ternary with value either 0, +1, or −1.

The far-field amplitude pattern of an annular Walsh filter is given by where

4. Numerical Results and Discussion

Figures 4–7 give the far-field amplitude distributions for annular Walsh filters , with central obscuration , and 0.9, respectively. In each figure, the amplitude distribution is shown for the annular Walsh filters of orders 0, 1, 2, and 3. Note that the amplitude values are normalized by the central amplitude for an unobstructed uniform aperture. For the sake of underscoring amplitude variations, the scales along the ordinate and abscissa of different curves could not be maintained the same; they were tailored to bring forth the specific features.

(a)

(b)

(c)

(d)

Figure 4: Far-field amplitude distributions for annular Walsh filters on an annular aperture with central obscuration for orders (a) , (b) , (c) , and (d) . The red graph in (a) represents amplitude distribution for an unobscured uniform aperture.

(a)

(b)

(c)

(d)

Figure 5: Far-field amplitude distributions for annular Walsh filters on an annular aperture with central obscuration for orders (a) , (b) , (c) , and (d) . The red graph in (a) represents amplitude distribution for an unobscured uniform aperture.

(a)

(b)

(c)

(d)

Figure 6: Far-field amplitude distributions for annular Walsh filters on an annular aperture with central obscuration for orders (a) , (b) , (c) , and (d) . The red graph in (a) represents amplitude distribution for an unobscured uniform aperture.

(a)

(b)

(c)

(d)

Figure 7: Far-field amplitude distributions for annular Walsh filters on an annular aperture with central obscuration for orders (a) , (b) , (c) , and (d) . The red graph in (a) represents amplitude distribution for an unobscured uniform aperture.

It is seen in Figures 4(a), 5(a), 6(a), and 7(a) that for the zero order annular Walsh filters, there is a central maximum of amplitude in the diffraction pattern, and the magnitude of the maximum is obviously dependent on the central obscuration . For the sake of comparison, amplitude distribution for an unobscured uniform aperture is shown alongside the same for zero order annular Walsh filters. It is well known that with higher values of obscuration, the central lobe narrows down significantly [20].

Figures 4–7 show that for Walsh filters of all orders other than zero, the central amplitude is zero. This is a consequence of orthogonality property of the Walsh functions. The dark center is surrounded by rings of oscillating amplitudes. For lower orders, this oscillation gradually decays after the first few rings, but for higher orders, the oscillations continue for many more rings. As expected from energy considerations, the amplitude of oscillation is significantly less in case of higher order filters compared to the lower order ones in each case of obscuration. From lower to higher obscuration, the ring with peak oscillatory amplitude gradually shifts away from the centre in case of higher order annular Walsh filters.

Use of ternary values of transmission 0, +1, and −1 in annular Walsh filters has opened up new possibilities for generation of interesting far-field patterns. Amplitude patterns shown in Figures 6(b)–6(d) and 7(b)–7(d) corresponding to annular Walsh filters of orders 1, 2, and 3 on a pupil with large obscuration have distinct characteristics that remain to be explored further for many practical applications, for example, optical encryption. Incidentally, it is also observed that for orders higher than 3, the typical nature of amplitude distribution is more pronounced, with the frequency of oscillations in amplitude becoming increasingly higher with concomitant decrease in peak amplitude.

5. Concluding Remarks

Optimal combinations of annular Walsh filters can help in realizing desired amplitude distributions on the far-field plane for applications using coherent illumination. An important application is to obtain “superresolution” effect by narrowing down of the central lobe and suppression of the side ripples. Solution of the related optimization problem can make effective use of orthogonality of annular Walsh functions. Incidentally, inverse problems on pupil synthesis for prespecified intensity characteristics in the far-field plane were tackled earlier for the case of incoherent illumination by the use of radial Walsh filters for unobscured apertures [18]. It is obvious that the availability of an additional degree of freedom, namely, the obscuration ratio , is likely to expand the scope of this inverse problem.

Ready availability of high efficiency spatial light modulators has facilitated the practical realization of phase filters. [22–24]. Nevertheless, it is obvious that, in general, practical implementation of filters with binary or ternary phase values is relatively easier than filters with continuously varying phase. Use of Walsh function-based analysis and synthesis of pupil plane filters has become useful, for it can circumvent the tricky problem of synthesizing continuously varying phase by providing alternative analytical treatments that deal with a finite number of discrete phase levels.

Finally, it remains to investigate the axial resolution characteristics of annular Walsh filters and also their usefulness in the synthesis of three-dimensional light fields [25].

Acknowledgment

Authors are indebted to the anonymous reviewers whose comments have significantly improved the quality of our presentation.

References

G. T. di Francia, “Nuovo pupille superrisolvente,” Atti della Fondazione Giorgio Ronchi, vol. 7, pp. 366–372, 1952.

E. H. Linfoot and E. Wolf, “Diffraction images in systems with an annular aperture,” Proceedings of the Physical Society B, vol. 66, no. 2, pp. 145–149, 1953.

A. Boivin, Théorie et Calcul des Figures de Diffraction de Revolution, Gauthier-Villars, Paris, France, 1964.

Z. S. Hegedus, “Annular pupil arrays—application to confocal scanning,” Optica Acta, vol. 32, pp. 815–826, 1985.

C. J. R. Sheppard and A. Choudhury, “Annular pupils, radial polarization, and superresolution,” Applied Optics, vol. 43, no. 22, pp. 4322–4327, 2004.

T. R. M. Sales and G. M. Morris, “Diffractive superresolution elements,” Journal of the Optical Society of America A, vol. 14, no. 7, pp. 1637–1646, 1997.

H. Luo and C. Zhou, “Comparison of superresolution effects with annular phase and amplitude filters,” Applied Optics, vol. 43, no. 34, pp. 6242–6247, 2004.

M. Yun, L. Liu, J. Sun, and D. Liu, “Three-dimensional superresolution by three-zone complex pupil filters,” Journal of the Optical Society of America A, vol. 22, no. 2, pp. 272–277, 2005.

M. Yun, M. Wang, and L. Liu, “Super-resolution with annular binary phase filter in the 4Pi-confocal system,” Journal of Optics A, vol. 7, no. 11, pp. 640–644, 2005.

J. L. Walsh, “A closed set of normal orthogonal functions,” The American Journal of Mathematics, vol. 55, pp. 5–24, 1923.

K. G. Beauchamp, Walsh Functions and Their Applications, Academic, New York, NY, USA, 1985.

C. J. R. Sheppard, “Binary phase filters with a maximally-flat response,” Optics Letters, vol. 36, no. 8, pp. 1386–1388, 2011.

H. F. Harmuth, Transmission of Information by Orthogonal Functions, Springer, Berlin, Germany, 1972.

H. C. Andrews, Computer Techniques in Image Processing, Academic, New York, NY, USA, 1970.

M. De and L. N. Hazra, “Walsh functions in problems of optical imagery,” Optica Acta, vol. 24, pp. 221–234, 1977.

L. N. Hazra and A. Guha, “Farfield diffraction properties of radial Walsh filters,” Journal of the Optical Society of America A, vol. 3, pp. 843–846, 1986.

L. N. Hazra and A. Banerjee, “Application of Walsh function in generation of optimum apodizers,” Journal of Optics, vol. 5, pp. 19–26, 1976.

L. N. Hazra, “A new class of optimum amplitude filters,” Optics Communications, vol. 21, no. 2, pp. 232–236, 1977.

L. N. Hazra, “Walsh filters for tailoring of resolution in microscopic imaging,” Micron, vol. 38, no. 2, pp. 129–135, 2007.

M. Born and E. Wolf, Principles of Optics, Pergamon, Oxford, UK, 1980.

H. H. Hopkins, “Canonical and real space coordinates used in the theory of image formation,” in Applied Optics and Optical Engineering, R. R. Shannon and J. C. Wyant, Eds., vol. 9, p. 307, Academic, New York, NY, USA, 1983.

E. Fernández, A. Márquez, S. Gallego, R. Fuentes, C. García, and I. Pascua, “Hybrid ternary modulation applied to multiplexing holograms in photopolymers for data page storage,” Journal of Lightwave Technology, vol. 28, pp. 776–783, 2010.

A. Au, C.-S. Wu, S.-T. Wu, and U. Efron, “Ternary phase and amplitude modulations using a twisted nematic liquid crystal spatial light modulator,” Applied Optics, vol. 34, no. 2, pp. 281–284, 1995.

S. Mukhopadhyay, S. Sarkar, K. Bhattacharya, and L. N. Hazra, “Polarisation phase shifting interferometric technique for phase calibration of a reflective phase spatial light modulator,” Optical Engineering, vol. 52, Article ID 035602, 2 pages, 2013.

R. Piestun and J. Shamir, “Synthesis of three-dimensional light fields and applications,” Proceedings of the IEEE, vol. 90, no. 2, pp. 222–244, 2002.