Two-dimensional window design

Windowing is a process where an index-limited sequence has its maximum energy concentrated in a finite frequency interval. This can be extended to an "N"-dimension where the "N"-D window has the limited support and maximum concentration of energy in a separable or non-separable "N"-D passband. The design of an "N"-dimensional window particularly a 2-D window finds applications in various fields such as spectral estimation of multidimensional signals, design of circularly symmetric and quadrantally symmetric non-recursive 2D filters, design of optimal convolution functions, image enhancement so as to reduce the effects of data-dependent processing artifacts, optical apodization and antenna array design.
Two-dimensional window.
Due to the various applications of multi-dimensional signal processing, the various design methodologies of 2-D windows is of critical importance in order to facilitate these applications mentioned above, respectively.
Consider a two-dimensional window function (or window array) formula_1 with its Fourier transform denoted by formula_2. Let formula_3 and formula_4 denote the impulse and frequency response of an ideal filter and formula_5 and formula_6 denote the impulse and frequency response of a filter approximating the ideal filter, then we can approximate formula_4 by formula_5. Since formula_3 has an infinite extent it can be approximated as a finite impulse response by multiplying with a window function as shown below
and in the Fourier domain
The problem is to choose a window function with an appropriate shape such that formula_6 is close to formula_4 and in any region surrounding a discontinuity of formula_4, formula_6 shouldn't contain excessive ripples due to the windowing.
2-D window function from 1-D function.
There are four approaches for generating 2-D windows using a one-dimensional window as a prototype.
Approach I
One of the methods of deriving the 2-D window is from the outer product of two 1-D windows, i.e., formula_16 The property of separability is exploited in this approach. The window formed has a square region of support and is separable in the two variables. In order to understand this approach, consider 1-D Kaiser window whose window function is given by
then the corresponding 2-D function is given by
where:
The Fourier transform of formula_1 is the outer product of the Fourier transforms of formula_21. Hence formula_22.
Approach II
Another method of extending the 1-D window design to a 2-D design is by sampling a circularly rotated 1-D continuous window function. A function is said to possess circular symmetry if it can be written as a function of its radius, independent of formula_23 i.e. formula_24
If "w"("n") denotes a good 1-D even symmetric window then the corresponding 2-D window function is
(where formula_26 is a constant) and
The transformation of the Fourier transform of the window function in rectangular co-ordinates to polar co-ordinates results in a Fourier–Bessel transform expression which is called as Hankel transform. Hence the Hankel transform is used to compute the Fourier transform of the 2-D window functions.
If this approach is used to find the 2-D window from the 1-D window function then their Fourier transforms have the relation
where:
and
This is the most widely used approach to design the 2-D windows.
2-D filter design by windowing using window formulations obtained from the above two approaches will result in the same filter order. This results in an advantage for the second approach since its circular region of support has fewer non-zero samples than the square region of support obtained from the first approach which in turn results in computational savings due to reduced number of coefficients of the 2-D filter. But the disadvantage of this approach is that the frequency characteristics of the 1-D window are not well preserved in 2-D cases by this rotation method. It was also found that the mainlobe width and sidelobe level of the 2-D windows are not as well behaved and predictable as their 1-D prototypes. While designing a 2-D window there are two features that have to be considered for the rotation. Firstly, the 1-D window is only defined for integer values of formula_32 but formula_33 value isn't an integer in general. To overcome this, the method of interpolation can be used to define values for formula_1 for any arbitrary formula_35 Secondly, the 2-D FFT must be applicable to the 2-D window.
Approach III
Another approach is to obtain 2-D windows by rotating the frequency response of a 1-D window in Fourier space followed by the inverse Fourier transform. In approach II, the spatial-domain signal is rotated whereas in this approach the 1-D window is rotated in a different domain (e.g., frequency-signal).
Thus the Fourier transform of the 2-D window function is given by
The 2-D window function formula_37 can be obtained by computing the inverse inverse Fourier transform of formula_38.
Another way to show the type-preserving rotation is when the relation formula_39 is satisfied. This implies that a slice of the frequency response of 2-D window is equal to that of the 1-D window where the orientation of formula_40 is arbitrary. In spatial domain, this relation is given by formula_41. This implies that a slice of the frequency response formula_38 is the same as the Fourier transform of the one-directional integration of the 2-D window formula_37.
The advantage of this approach is that the individual features of 1-D window response formula_44 are well preserved in the obtained 2-D window response formula_38. Also, the circular symmetry is improved considerably in a discrete system. The drawback is that it's computationally inefficient due to the requirement of 2-D inverse Fourier transform and hence less useful in practice.
Approach IV
A new method was proposed to design a 2-D window by applying the McClellan transformation to a 1-D window. Each coefficient of the resulting 2-D window is the linear combination of coefficients of the corresponding 1-D window with integer or power of 2 weighting.
Consider a case of even length, then the frequency response of the 1-D window of length N can be written as
Consider the McClellan transformation:
which is equivalent to
Substituting the above, we get the frequency response of the corresponding 2-D window
From the above equation, the coefficients of the 2-D window can be obtained.
To illustrate this approach, consider the Tseng window. The 1-D Tseng window of formula_50 weights can be written as
By implementing this approach, the frequency response of the 2-D McClellan-transformed Tseng window is given by
where formula_1 are the 2-D Tseng window coefficients.
This window finds applications in antenna array design for the detection of AM signals.
The advantages include simple and efficient design, nearly circularly symmetric frequency response of the 2-D window, preserving of the 1-D window prototype features. However, when this approach is used for FIR filter design it was observed that the 2-D filters designed were not as good as those originally proposed by McClellan.
2-D window functions.
Using the above approaches, the 2-D window functions for few of the 1-D windows are as shown below. When Hankel transform is used to find the frequency response of the window function, it is difficult to represent it in a closed form. Except for rectangular window and Bartlett window, the other window functions are represented in their original integral form. The two-dimensional window function is represented as formula_54 with a region of support given by formula_55 where the window is set to unity at origin and formula_56 for formula_57 Using the Hankel transform, the frequency response of the window function is given by
where formula_59 is Bessel function identity.
Rectangular window.
The two-dimensional version of a circularly symmetric rectangular window is as given below
The window is cylindrical with the height equal to one and the base equal to 2a. The vertical cross-section of this window is a 1-D rectangular window.The frequency response of the window after substituting the window function as defined above, using the Hankel transform, is as shown below
Bartlett window.
The two-dimensional mathematical representation of a Bartlett window is as shown below
The window is cone-shaped with its height equal to 1 and the base is a circle with its radius 2a. The vertical cross-section of this window is a 1-D triangle window.The Fourier transform of the window using the Hankel transform is as shown below
Kaiser window.
The 2-D Kaiser window is represented by
The cross-section of the 2-D window gives the response of a 1-D Kaiser Window function. The Fourier transform of the window using the Hankel transform is as shown below