Spatially variant apodization is a digital image processing technique for eliminating sidelobes produced by Fourier transform of finite data sequences without compromising mainlobe width. This process allows each sample or pixel in an image to receive its own frequency domain aperture amplitude weighting function from an infinite number of possible weighting functions. In its simplest form the weight is a function of the negative of the current sample divided by the sum of the neighboring samples, and the function is limited to a predetermined range such as the range between zero and one half.

FIELD OF THE INVENTION 
This present invention relates to decreasing sidelobes when performing 
signal compression using matched filters or transforms and more 
particularly to the reduction of sidelobes resulting from transform of 
finite data sequences. 
BACKGROUND OF THE INVENTION 
Signal compression is a common operation which is performed in many 
systems, including radar. The compression is often performed as a 
transform of domain, such as from the time domain to the frequency domain. 
The accuracy of the compression is limited by the finite amount of signal 
that can be collected. In the case of imaging radars, a signal consists of 
one or more sine waves in time that must be transformed into the spatial 
domain in order to determine their frequency, magnitude, and sometimes 
phase. The most common method for transformation is the Fourier transform. 
The Fourier transform of a limited duration sine wave produces a waveform 
that can be described by a sinc function (FIG. 1). The sinc function has a 
mainlobe which contains the peak and has a width up to the first zero 
crossing, and a set of sidelobes comprising the oscillating remainder on 
both sides of the mainlobe. In radar and some other fields, the composite 
function of the mainlobe and the sidelobes is termed the impulse response 
(IPR) of the system. The location of the center of the sinc function is 
related to the frequency of the sine wave. If there are more than one sine 
wave present in the signal being analyzed, they will appear in the output 
at other locations. The resolution is related to the width of the 
mainlobe. The presence of sidelobes reduces the ability to discriminate 
between sinc functions. 
Traditionally, the sidelobes of the impulse response have been reduced by 
multiplying the signal prior to compression by an amplitude function that 
is a maximum at the center and tending toward zero at the edges, as 
typified by a Hanning weighting function shown in FIG. 2. Sidelobe 
reduction by amplitude multiplication is called "weighting" or, sometimes, 
"apodization". Unfortunately, employing that kind of apodization to reduce 
sidelobes also results in the broadening of the mainlobe which degrades 
the resolution of the system. 
One family of apodization functions is termed "cosine-on-pedestal". Hanning 
(50% cosine and 50% pedestal, as shown in FIG. 2) and Hamming (54% cosine 
and 46% pedestal) are two of the most popular. Hanning weighting reduces 
the peak sidelobe from -13 dB of the mainlobe's peak to -32 dB but it also 
doubles the mainlobe width (FIG. 3). 
The equivalent of apodization can also be performed in the output domain by 
convolution. In the case of digitally sampled data, convolution is 
performed by executing the following operation on each point in the 
sequence: multiply each sample by a real-valued weight. which is dependent 
on the distance from the point being processed. 
Any of the cosine-on-pedestal family of apodizations is especially easy to 
implement by convolution when the transform is of the same length as the 
data set, i.e., the data set is not padded with zeros before 
transformation. In this case, the convolution weights are non-zero only 
for the sample itself and its two adjacent neighbors. The values of the 
weights vary from [0.5, 1.0, 0.5] in the case of Hanning to [0.0, 1.0, 
0.0] in the case of no apodization. Different cosine-on-pedestal 
apodization functions have different zero crossing locations for the 
sidelobes. The Hanning function puts the first zero crossing at the 
location of the second zero crossing of the unweighted impulse response. 
Not shown in FIGS. 1 and 3, the signs of the IPRs are opposite for all 
sidelobes when comparing unapodized and Hanning apodized signals. 
To improve the process, a method called dual-apodization has been 
developed. In this method, the output signal is computed twice, once using 
no apodization and a second time using some other apodization which 
produces low sidelobes. Everywhere in the output, the two values are 
compared. The final output is always the lesser of the two. In this way 
the optimum mainlobe width is maintained while the sidelobes are generally 
lowered. 
An extension to dual apodization is multi-apodization. In this method, a 
number of apodized outposts are prepared using a series of different 
apodizations, each of which have zero-crossings at different locations. 
The final output is the least among the ensemble of output apodized values 
at each output point. In the limit of an infinite number of apodizations, 
all sidelobes will be eliminated while the ideal mainlobe is preserved. 
The final embodiment of this invention occurred when a method was 
discovered that could compute, for each sample in a sidelobe region, which 
of the cosine-on-pedestal functions provided the zero crossing from among 
the potentially infinite number of possible apodizations. This method is 
called spatially variant apodization (SVA). The method computes the 
optimum convolution weight set for each sample using a simple formula 
based on the value of the sample and two of its neighbors. Under noise 
free conditions, well separated compressed signals show only the 
mainlobes, and all sidelobes are removed. Under the usual noisy 
conditions, the output signal to background ratios are improved and the 
sidelobes are greatly reduced. 
SUMMARY OF THE INVENTION 
When a signal of finite duration undergoes a signal compression via a 
transform, such as a Fourier transform, sidelobes develop that obscure 
details in the output data. This invention is a method for attenuating or 
eliminating the sidelobes without compromising the resolution of the 
signal. The first step is a compression of the signal using little or no 
apodization. The second step is to determine the convolution weights for 
each output sample. The center weight is unity. The outer two are the same 
and are computed as follows: 1) the two adjacent samples are summed, 2) 
the sum is divided into the value of the center sample, and 3) the 
resulting value is limited to a specific range depending on the 
application, e.g. 0 to 0.5. In the final step, the sample is convolved 
using the computed weight set. 
This method has variants depending on the type of compression, the type of 
signal, and the application. Fourier transforms are a standard method of 
compressing sine waves but other transforms are also used, including 
cosine, Hartley and Haddamard. Matched filter compression is also used in 
the cases in which the signal is not a sine wave but some other expected 
waveform. In each compression method, one must search for the convolution 
set that implements a set of apodizations which affect the magnitudes and 
signs of the sidelobes. 
The type of signal can be real or complex, one dimensional or 
multidimensional. For real-valued functions, there is only one channel to 
process. Complexed-valued functions have an in-phase (I) channel and a 
quadrature (Q) channel. Spatially variant apodization can be applied to 
the I and Q channels independently or can, with a slight modification in 
the equation, handle the joint I/Q pair. 
When the signals are two (or higher) dimensional, there are again several 
ways to perform the spatially variant apodization. The first is to apodize 
in one dimension at a time in a serial manner. The second is to apodize 
each dimension, starting from the same unapodized process. The results of 
apodizing in the individual dimensions are combined by taking the minimum 
output among the individual apodizations for each output sample. 
In summary, the family of spatially variant apodization methods select a 
different and optimum apodization at each output position in order to 
minimize the sidelobes arising from signal compressions of finite data.

DETAILED DESCRIPTION OF THE INVENTION 
Spatially variant apodization (SVA) allows each pixel in an image to 
receive its own frequency domain aperture amplitude weighting function 
from an infinity of possible weighting functions. In the case of synthetic 
aperture radar (SAR), for example, SVA effectively eliminates 
finite-aperture induced sidelobes from uniformly weighted data while 
retaining nearly all of the good mainlobe resolution and clutter texture 
of the unweighted SAR image. 
FIG. 1 depicts the graph of a sinc function waveform. This serves to model 
the impulse response of performing a Fourier transform on a set of 
finite-aperture data. The mainlobe 10 carries the information from the 
original signal. To maintain the resolution of the image, the mainlobe 10 
must not be widened during the apodization of the image. The sidelobes 12 
do not carry any information about the original signal. Instead, they 
serve to obscure the neighboring details which have weaker signal 
strengths than the sidelobes. 
Spatially variant apodization was developed for synthetic aperture radar in 
response to the problems inherent to finite aperture systems as described 
above. However, there are many different embodiments for spatially variant 
apodization in areas of imagery, digital signal processing, and others. 
FIG. 4 is a simplified block diagram of a synthetic aperture radar system 
utilizing spatially variant apodization. The system can be broken into 
five smaller sections: data acquisition 14, data digitizing 16, digital 
image formation processing 18, detection 20 and display 22. 
Data acquisition 14 for synthetic aperture radar comprises a transmitter 26 
to generate a radio frequency signal to be broadcast by an antenna 24. The 
reflected radio signals returning to the antenna 24 are sent to the 
receiver, where a complex pair of signals are formed and sent to an analog 
to digital converter 16. 
The analog to digital converter 16 samples and digitizes each signal and 
passes the data to the digital processor 18. In the digital processor 18, 
the first function performed is that of motion compensation 30. Since this 
type of system is used in moving aircraft to survey surface features, the 
motion of the plane must be taken into consideration so that the image is 
not distorted. After motion compensation 30, the signals are processed by 
polar formatting circuitry or algorithms to format the data in such a 
manner so that a coherent two dimensional image can be formed by a Fourier 
transform. The next step in digital processing is to transform the data 
from the frequency domain to the space domain via a Fast Fourier Transform 
(FFT) 34. It is at this step that sidelobes are produced in the image. The 
final step in the digital processor 18 is to perform spatially variant 
apodization 36 on the complex data sets. 
Following the digital processing 18, detection 20 takes place to form the 
final signal which drives the display 22. Detection 20 comprises 
determining the magnitude of the complex image. From this data a two 
dimensional image can be displayed on a CRT or on film. 
It is well known that cosine-on-pedestal frequency domain weighting 
functions can be implemented using a 3-point convolver on complex, Nyquist 
sampled imagery. The family of cosine-on-pedestal weightings range from 
uniform weighting to Hanning weighting. Hamming weighting is a special 
case of cosine-on-pedestal which nulls the first sidelobe. Similarly, any 
unweighted aperture sinc function sidelobe can be nulled using one of the 
family of cosine-on-pedestal weighting functions. 
The spatially variant apodization algorithm has a number of different 
forms. The derivation for the most basic form is as follows. Let g(i) 
denote the current sample or pixel, g(i-1) denote the previous sample or 
pixel, and g(i+1) denote the following sample in one dimension of either 
the real (I) or imaginary (Q) parts of a uniformly weighted 
Nyquist-sampled image. Using a 3-point convolver to achieve a given 
cosine-on-pedestal aperture weighting, g(i) is replaced by g'(i) as 
follows: 
EQU g'(i)=w(i)g(i-1)+g(i)+w(i)g(i+1). (1) 
This equation is the convolution of g(i) with the sampled impulse response 
due to a raised-cosine aperture-weighting function. As w(i) varies from 0 
to 1/2, the frequency domain amplitude weighting varies from 
cosine-on-zero pedestal (Harming) at w(i)=1/2 to uniform weighting at 
w(i)=0. The center convolver weight is always unity because it is 
desirable to normalize the peak point-target responses for the family of 
cosine-on-pedestal weightings. 
It is desired therefore, to find the w(i) which minimizes 
.vertline.g'(i).vertline. subject to the constraints 
0.ltoreq.w(i).ltoreq.1/2. The unconstrained w(i) that gives the minimum 
is: 
##EQU1## 
If w(i) in Eq.(2) is substituted into Eq.(1) then g'(i)=0 is the 
unconstrained solution. Therefore, applying the constraints, g'(i)=0 is 
the solution wherever 0.ltoreq.w(i).ltoreq.1/2. However, g'(i) can be 
nonzero wherever w(i)&lt;0 or w(i)&gt;1/2, i.e., 
##EQU2## 
The process given by Eqs. (1)-(4) can be applied according to the block 
diagram in FIG. 5. The real horizontal convolver 38 and the imaginary 
horizontal convolver 40 can each be performed sequentially using two 
parallel processors to produce two modified data sets. Next, the real 
vertical convolver 42 and the imaginary vertical convolver 44 can each be 
performed sequentially on the modified data using two parallel processors 
to create further modified data sets. 
Spatially variant apodization can also be formulated by minimizing (I.sup.2 
+Q.sup.2) jointly. The result is one weight applied jointly to both I and 
Q. Equation (2) becomes: 
##EQU3## 
The block diagram for applying this form of SVA is depicted in FIG. 6. 
First a horizontal convolver 46 is applied to both components I, Q of the 
current sample i using the weighting having the form of equation (5) to 
obtain the modified sample value (I'.sub.i, Q'.sub.i). The convolver 
equation is 
EQU (I'.sub.i,Q'.sub.i)=w(i)*(I.sub.i-1,Q.sub.i-1)+(I.sub.i,Q.sub.i)+w(i)* 
(I.sub.i+1,Q.sub.i+1). 
Then a vertical convolver 48 is applied to both components of the resultant 
(I'.sub.i,Q'.sub.i) again using a weighting having the form of equation 
(5) and the same form of convolver equation. Another approach to 
processing the two dimensional data is to operate both the horizontal and 
vertical convolvers on the current sample value to calculate two modified 
values and then select the minimum of the two modified values. 
The upper and lower limits for weights as set forth heretofore are 
theoretical ideals which should deliver the optimal results. However, 
actual systems which may utilize spatially variant apodization are not 
ideal systems, and the designers may prefer some result other than the 
theoretical ideal result. For this reason the invention shall not be 
limited to only those theoretical limits set forth herein. 
The sampling rate for SAR systems is typically not limited to the Nyquist 
rate. Any integer or non-integer sampling rate may be chosen as a matter 
of design choice. If the sampling rate is an N integer multiple of the 
Nyquist rate then the convolver and weight algorithms should utilize the 
current pixel or sample, g(i), and the Nth neighboring pixels or samples 
on each side, g(i-N) and g(i+N). For cases calculating I and Q jointly, 
the weight is determined using the following expression: 
##EQU4## 
denoted w(i) where I.sub.i and Q.sub.i represent the current real and 
imaginary components, respectively, of the current sample, I.sub.i-N and 
Q.sub.i-N represent the real and imaginary components, respectively, of 
the Nth sample preceding the current sample, and I.sub.i+N and Q.sub.i+N 
represent the real and imaginary components, respectively, of the Nth 
sample following the current sample. The usual limits apply to this case 
except that, in order to obtain good results when using sampling rates 
that are non-integer multiples of the Nyquist rate, the upper limit for 
the convolver weight must be modified. Let p be the percent zero-padding 
in the system. The required maximum for w(i) can be calculated in the 
following manner: 
##EQU5## 
Finally, a new complex value, (I.sub.i ',Q.sub.i '), for the current 
sample, (I.sub.i,Q.sub.i), is determined according to the expression: 
EQU (I.sub.i ',Q.sub.i ')=w(i)*(I.sub.i-N,Q.sub.i-N)+(I.sub.i,Q.sub.i)+w(i)* 
(I.sub.i+N,Q.sub.i+N). 
The techniques that have been addressed here have been two dimensional. The 
algorithm is applied to only one dimension at a time. As few as one 
dimension may be processed using spatially variant apodization. The 
techniques may also be expanded to as include as many dimensions as is 
necessary. As stated above, the convolving process for more than one 
dimension can be employed serially with each convolution operating on the 
resultant of the previous convolution, or alternatively each convolution 
can operate on the initial samples and the minimum value selected. 
FIG. 7 illustrates the effect of the SVA algorithm on a data set having two 
peaks. The solid line is the sum of two sincs separated by 3.5 samples. 
The output of the SVA algorithm is shown in the dashed line which reveals 
the two distinct peaks with no sidelobes and no broadening of the 
mainlobes. The same result was reached using either the independent 
treatment of I and Q, or the joint treatment. 
Another embodiment of the method of minimizing I and Q jointly, uses 
complex values of the data in the convolver equation and achieves better 
suppression of the sidelobes. In this method the weight is expressed by 
##EQU6## 
subject to the constraints if w(i)&lt;W.sub.min ; 
g'(i)=[I(i)+jQ(i)][1-w.sub.min /w(i)]; 
if w.sub.min .ltoreq.w(i).ltoreq.w.sub.max ; g'(i)=0; and 
if w(i)&gt;w.sub.max ; g'(i)=[I(i)+jQ(i)][1-w.sub.max /w(i)]. 
Thus the form of the convolver equation to determine the modified point g' 
as well as the value of the weight w(i) depends on the range of the 
weight. In the ideal case the upper weight limit w.sub.max =1/2 and the 
lower weight limit w.sub.min =0 as in the previous embodiments.