# Far-field Diffraction Properties of Annular Walsh Filters

Authors: Pubali Mukherjee; Lakshminarayan Hazra
Journal: Advances in Optical Technologies
 (2013)

---

## 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.


---

## Body

## 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, Section3 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Ψnε(r) of index n≥0 and argument r over an annular region with ε and 1 as inner and outer radii, respectively, it is necessary to express the integer n in the form
(1)n=∑m=0n~-1km2m,
where km are the bits, 0 or 1 of the binary numeral for n, and (2n~) is the power of 2 that just exceeds n. For all r in (ε,1), Ψnε(r) is defined as
(2)Ψnε(r)=∏m=0n~-1sgn{cos[km2mπ(r2-ε2)(1-ε2)]},
where
(3)sgn(x)={+1,x>0,0,x=0,-1,x<0.
The orthogonality condition implies that
(4)∫ε1Ψmε(r)Ψnε(r)rdr=1-ε22δmn,
where δmn is the Kronecker delta defined as
(5)δmn={0,m≠n,1,m=n.
Figure 1 shows the first four annular Walsh functions for central obscuration ε=0.3 in two dimensions. Figure 2 presents values of the functions Ψn0.3(r), n=0,…,3, along the radius in an azimuthal direction. It should be noted that the order of the functions n is equal to the number of zero crossings, or sign changes of the function in the interval (0.3,1), and locations of the points of zero crossings for members of the set of functions Ψnε(r), n=0,…,3 are given by
(6)ri=[(4-i)ε2+i]4,i=1,2,3.
In general, for the first N, where N=2b and b is a positive integer, and annular Walsh functions Ψnε(r), n=0,1,…,(N-1), the zero crossings are located at(7)ri=[(N-i)ε2+i]N,i=1,2,…,(N-1).
The inner and outer radii of the annulus is r0=ε and rN=1. Note that the set of (N-1) zero crossing locations, ri, i=1,2,…,(N-1), 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Ψnε(r) in two dimensions of order n=0,1,2,3 for central obscuration ε=0.3.Figure 2
Annular Walsh functionsΨnε(r) of order n=0,1,2,3 along radius r for central obscuration ε=0.3.For computational purposes it is often convenient to express an annular Walsh functionΨnε(r) as
(8)Ψnε(r)=∑j=1NhnjζjN(r),
where N=2n~. For a particular value of n, the integer n~ is taken such that N just exceeds n. ζjN(r) are zero-one functions, also known as Walsh block functions, defined as
(9)ζjN(r)={1,rj-1≤r≤rj,0,otherwise.For all annular Walsh functions of ordern<N, the same values of rj, j=0,N are to be used. They are unique for the particular value of N, and their values are given in (7). hnj are the elements of a (N×N) Hadamard matrix whose elements are +1 or −1 [13]. Table 1 gives the Hadamard matrix for N=22. The Hadamard matrices for other values of N can be obtained from the defining relation (2) of annular Walsh functions.Table 1
The Hadamard matrix forN=22.
+1
+1
+1
+1
+1
+1
−1
−1
+1
−1
−1
+1
+1
−1
+1
−1


## 3. Far-Field Amplitude Distribution

Figure3 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 F(p) at a point Q′ 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]:
(10)F(p)=∫01f(r)J0(pr)rdr,
where f(r) is the circularly symmetric pupil function. The variable r is the fractional coordinate for a point P′ on the pupil sphere of center O′ and radius O′E′, where E′ is the axial position of the exit pupil. r is equal to the ratio (r′/h′), where r′ is the perpendicular distance of the point P′ from the axis and h′ is the maximum value of r′, so that 0≤r≤1. p is the reduced diffraction variable given by
(11)p=2πλ(n′sinα′)ξ′,
where α′ is the semiangular aperture of the imaging system, n′ is the refractive index of the image space, 2π/λ is the propagation constant, and ξ′(=O′Q′) is the actual geometrical distance of point Q′ on the image plane from the center O′ 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 f(r) is given by
(12)f(r)={0,0≤r<ε,Ψnε(r)ε≤r<1.
It may be noted that f(r) is binary (value either 0 or +1) only in the case of zero order annular Walsh function Ψ0ε(r); for all other orders, f(r) is ternary with value either 0, +1, or −1.The far-field amplitude patternFnε(p) of an annular Walsh filter Ψnε(r) is given by
(13)Fnε(p)=∫ε1Ψnε(r)J0(pr)rdr=∑j=1NhnjΓj(p),
where
(14)Γj(p)=∫ε1ζjN(r)J0(pr)rdr=[rj2J1(prj)prj-rj-12(prj-1)prj-1].


## 4. Numerical Results and Discussion

Figures4–7 give the far-field amplitude distributions for annular Walsh filters Ψnε(r), with central obscuration ε=0.3,0.5,0.7, 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.Far-field amplitude distributions for annular Walsh filtersΨnε(r) on an annular aperture with central obscuration ε=0.3 for orders (a) n=0, (b) n=1, (c) n=2, and (d) n=3. The red graph in (a) represents amplitude distribution for an unobscured uniform aperture.
(a)
(b)
(c)
(d)Far-field amplitude distributions for annular Walsh filtersΨnε(r) on an annular aperture with central obscuration ε=0.5 for orders (a) n=0, (b) n=1, (c) n=2, and (d) n=3. The red graph in (a) represents amplitude distribution for an unobscured uniform aperture.
(a)
(b)
(c)
(d)Far-field amplitude distributions for annular Walsh filtersΨnε(r) on an annular aperture with central obscuration ε=0.7 for orders (a) n=0, (b) n=1, (c) n=2, and (d) n=3. The red graph in (a) represents amplitude distribution for an unobscured uniform aperture.
(a)
(b)
(c)
(d)Far-field amplitude distributions for annular Walsh filtersΨnε(r) on an annular aperture with central obscuration ε=0.9 for orders (a) n=0, (b) n=1, (c) n=2, and (d) n=3. The red graph in (a) represents amplitude distribution for an unobscured uniform aperture.
(a)
(b)
(c)
(d)It is seen in Figures4(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].Figures4–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 Figures6(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].



---

*Source: 360450-2013-07-24.xml*