Methods and systems for suppressing clutter in radar systems

09971027 ยท 2018-05-15

Assignee

Inventors

Cpc classification

International classification

Abstract

Methods and systems for suppressing clutter, for example, ground clutter, in radar systems are provided. The methods and systems can be employed in radar systems having an antenna system and at least two receive beams, for example, a main beam and an auxiliary beam. The methods include receiving data streams from each of the at least two receive beams, where each data stream is associated with range bins and include data representing clutter, and, before or after Doppler filtering, generating an adaptive weight from summations of the data streams for each of the range bins, and applying the generated weight to at least one of the data streams to provide Doppler filtered and spatially nulled data streams that can be used to more accurately identify targets, such as, aircraft.

Claims

1. A method for suppressing clutter when detecting targets of interest with a radar system comprising an antenna system and at least two receive beams, the method comprising: receiving I/Q data streams from each of the at least two receive beams, each of the I/Q data streams associated with range bins and including data representing clutter; Doppler filtering the I/Q data streams at each of the range bins and thereby generating for each of the I/Q data streams multiple Doppler-filtered I/Q data streams for each of the range bins; generating a weight from summations of the I/Q data streams for each of the range bins; applying the generated weight to at least one of the multiple Doppler-filtered I/Q data streams to provide at least one weighted I/Q data stream; combining the at least one weighted I/Q data stream with one of the other multiple Doppler-filtered I/Q data streams to thereby generate at least one weighted Doppler-filtered I/Q data stream representing spatial nulling of clutter for each of the range bins; and using the at least one weighted Doppler-filtered I/Q data stream representing spatial nulling of clutter for each of the range bins to detect targets.

2. The method as recited in claim 1, wherein generating the weight comprises computing the weight based on a summation of the I/Q data streams across pulses for each of the I/Q data streams.

3. The method as recited in claim 1, wherein the at least two receive beams comprise at least a main beam and an auxiliary beam.

4. The method as recited in claim 1, wherein the method further comprises comparing at least one of the summations of the I/Q data streams to a predetermined threshold.

5. The method as recited in claim 4, wherein, when the at least one of the summations exceeds the predetermined threshold, equating the weight to zero.

6. The method as recited in claim 1, wherein the method further comprises comparing an absolute value of the weight to a predetermined threshold.

7. The method as recited in claim 6, wherein generating the weight from the summations of the I/Q data streams comprises generating the weight from a quotient including the summations of the I/Q data streams, and wherein, when the absolute value of the quotient exceeds the predetermined threshold weight, decreasing the quotient.

8. The method as recited in claim 6, wherein generating the weight from the summations of the I/Q data streams comprises equating the weight to a quotient of one of the summations of the I/Q data streams divided by a sum of another of the summations of the I/Q data streams and a predetermined constant, and wherein comparing the absolute value of the weight to the predetermined threshold comprises comparing the absolute value of the quotient to the predetermined threshold.

9. The method as recited in claim 8, wherein, when the absolute value of the quotient is greater than the predetermined threshold, decreasing the quotient.

10. The method as recited in claim 1, wherein generating the weight from the summations of the I/Q data streams comprises generating the weight, W, from the quotient including the summations of the I/Q data streams, .sub.M and .sub.A, by the equation:
W=.sub.M/(.sub.A+); wherein comprises a predetermined constant.

11. A method for suppressing clutter when detecting targets of interest in a radar system comprising an antenna system and at least two receive beams, the method comprising: receiving I/Q data streams from each of the at least two receive beams, each of the I/Q data streams associated with range bins and including data representing clutter; generating a weight from summations of the I/Q data streams for each of the range bins; applying the generated weight to at least one of the I/Q data streams to provide at least one weighted I/Q data stream; combining the at least one weighted I/Q data stream with one of the other I/Q data streams to thereby generate at least one weighted I/Q data stream representing spatial nulling of clutter for each of the range bins; Doppler filtering the at least one weighted I/Q data stream for each of the range bins to generate at least one Doppler-filtered, spatially-nulled data stream having suppressed clutter for the range bins; and using the at least one Doppler-filtered, spatially-nulled data stream having suppressed clutter for the range bins to detect targets.

12. The method as recited in claim 11, wherein generating the weight comprises computing the weight based on a summation of the I/Q data streams across pulses for each of the I/Q data streams.

13. The method as recited in claim 11, wherein the at least two receive beams comprise at least a main beam and an auxiliary beam.

14. The method as recited in claim 11, wherein the method further comprises comparing at least one of the summations of the I/Q data streams to a predetermined threshold.

15. The method as recited in claim 14, wherein, when the at least one of the summations exceeds the predetermined threshold, equating the weight to zero.

16. The method as recited in claim 11, wherein the method further comprises comparing an absolute value of the weight to a predetermined threshold.

17. The method as recited in claim 16, wherein generating the weight from the summations of the I/Q data streams comprises generating the weight from a quotient including the summations of the I/Q data streams, and wherein, when the absolute value of the quotient exceeds the predetermined threshold weight, decreasing the quotient.

18. The method as recited in claim 16, wherein generating the weight from the summations of the I/Q data streams comprises equating the weight to a quotient of one of the summations of the I/Q data streams divided by a sum of another of the summations of the I/Q data streams and a predetermined constant, and wherein comparing the absolute value of the weight to the predetermined threshold comprises comparing the absolute value of the quotient to the predetermined threshold.

19. The method as recited in claim 18, wherein, when the absolute value of the quotient is greater than the predetermined threshold, decreasing the quotient.

20. The method as recited in claim 11, wherein generating the weight from the summations of the I/Q data streams comprises generating the weight, W, from the quotient including the summations of the I/Q data streams, .sub.M and .sub.A, by the equation:
W=.sub.M/(.sub.A+); wherein comprises a predetermined constant.

Description

BRIEF DESCRIPTION OF THE DRAWINGS

(1) The subject matter, which is regarded as the invention, is particularly pointed out and distinctly claimed in the claims at the conclusion of this specification. The foregoing and other aspects, features, and advantages of the invention will be readily understood from the following detailed description of aspects of the invention taken in conjunction with the accompanying drawings in which:

(2) FIG. 1 is a schematic diagram of a radar antenna system that may be adapted to implement aspects of the present invention.

(3) FIGS. 2, 3, and 4 are schematic diagrams illustrating the manipulation of the signals received by an antenna system according to aspects of the invention.

(4) FIG. 5 is a schematic block diagram of a signal processing antenna system that may be used for the antenna system shown in FIG. 1 according to aspects of the present invention.

(5) FIG. 6 is a detailed flow diagram with annotation for the spatial nulling practiced in the spatial nulling module shown in FIG. 5 according to an aspect of the present invention.

(6) FIG. 7 is a duplication of FIG. 4, and shows a representative plot of the signals received corresponding to range bins, N.sub.R, for pulses, N.sub.P, shown in FIG. 3.

(7) FIG. 8 is a representative plot of the sum of the signals received in each range bin for each pulse corresponding to range bins for pulse signals shown in FIG. 7.

(8) FIG. 9 is a schematic block diagram of another signal processing antenna system that may be used for the antenna system shown in FIG. 1 according to another aspect of the present invention.

(9) FIG. 10 is a detailed flow diagram with annotation for the spatial nulling practiced in spatial nulling module shown in FIG. 9 according to an aspect of the present invention.

(10) FIG. 11 is a representative plot of a typical transmitted radar beam and two receive beams according to one aspect of the invention.

(11) FIG. 12 is a representative plot that illustrates the spatial nulling of a received beam shown in FIG. 11 before and after spatial nulling according to aspects of the invention.

DETAILED DESCRIPTION OF THE INVENTION

(12) FIG. 1 is a schematic diagram of a radar antenna system 10, for example, a pulse-Doppler radar system, that may be adapted to implement aspects of the present invention. The structures shown in FIG. 1, and in the other figures herein, are not intended to be illustrated to scale, but are dimensioned as shown to facilitate illustration, description, and operation of aspects of the invention. It is to be understood that like reference numbers in all figures provided herein refer to similar, if not identical, structures, for example, to identical structures to those structures shown in FIG. 1.

(13) As shown in FIG. 1, in one aspect, system 10 includes one or more transmitting and receiving antenna systems 12 adapted to transmit one or more radar beams (not shown) to detect one or more targets 18, for example, an airplane. Radar system 10 may comprise, for example, a commercial air traffic control system or a military target detection system, among other radar applications. As known in the art, the transmitted radar beam reflects from targets 18 and other structures to yield two or more reflected beams 14 and 16, which are received by antenna system 12. The at least two receive beams may typically be simultaneous receive beams stacked in the vertical or elevation dimension, as known in the art. In one aspect, when the cluster of simultaneous beams is electronically steered, the received beams may be stacked vertically in sine-space. As known in the art, sine-space and elevation are similar when the cluster of beams is steered near broadside, but may differ significantly when the beam cluster of simultaneous beams is electronically steered away from broadside. The received beams 14 and 16 may typically be oriented at different directions or elevations, that is, received beams 14 and 16 may be off-set or squinted as known in the art, that is, oriented at different elevation angles with respect to antenna system 12.

(14) FIG. 1 also illustrates representative sources of clutter or ground clutter 20, that is, buildings, structures, terrain, trees, mountains, and/or bodies of water, and the like, that may interfere with the detection of targets 18. According to aspects of the invention, the influences of the presence of such sources of clutter on the desired detection of targets 18 can be suppressed, minimized, and/or eliminated, and thus enhance the accuracy and/or efficiency of the detection of targets 18.

(15) As shown in FIG. 1, off-set received beams 14 and 16 received by antenna system 12 of system 10 may be directed at different elevations: beam 14 may be oriented at a higher elevation and be less exposed to the influence of clutter 20, and beam 16 may be oriented at a lower elevation and be more exposed to the influence of clutter 20. In one aspect, the higher elevation received beam 14 may be referred to as the main beam and the lower elevation received beam 16 may be referred to as the auxiliary beam or the aux beam.

(16) As known in the art, although beams 14 and 16 are represented in FIG. 1 as single lobes to facilitate illustration, beams 14 and 16 typically comprise a series of pulses of electromagnetic (EM) radiation or energy, for example, radio waves at a given frequency. For example, beams 14 and 16 may comprise a plurality of pulses of EM radiation defined by a pulse repetition frequency (PRF), that is, a rate at which pulses are transmitted from antenna system 12. The plurality of pulses emitted by antenna system 12 may be referred to as a burst of pulses, where a burst may comprise a plurality of EM pulses. Correspondingly, the EM of reflected received beams 14 and 16 from target 18 and clutter 20 received by antenna system 12 are also received in the form of pulses and the received pulses associated with bursts. As is known in the art, the power of EM beam received by antenna system 12 is expressed in the form decibels [dB]. According to aspects of the present invention, the pulses of beams 14 and 16 are received by antenna system 12 and are processed to minimize the effect of clutter 20 on the detection of targets 18.

(17) As known in the art, the beams 14 and 16 from one or more targets 18 and clutter 20 are received by antenna system 12, and the received beams 14 and 16 are processed, for example, modulated, amplified, down-converted, and/or passed through detector circuitry to estimate the location and/or speed of targets 18.

(18) FIGS. 2, 3, and 4 are schematic diagrams illustrating some of the processing and/or manipulation of the signals received by antenna system 12 according to aspects of the invention. Antenna system 12 and target 18 shown in FIG. 1 also appear in FIG. 2.

(19) As known in the art, antenna system 12 identifies the location of the coordinates of target 18 by range (that is, the distance from antenna 12) as indicated by arrow 22 in FIG. 2; azimuth (that is, left to right location relative to a datum, for example, the orientation of antennas 12) as indicated by arrow 24 in FIG. 2; and elevation (that is, height above a datum, for example, the elevation above antenna system 12) as indicated by arrow 26 in FIG. 1.

(20) FIG. 2 is a plan view of the relative location of antenna system 12 and target 18 for an azimuth 24 that characterizes antenna system 12, for example, 1 to 2 degrees, and the range (distance) associated with antenna system 12. As shown in FIG. 2, in this aspect, range 22 is subdivided into range intervals or spaces 28.

(21) FIG. 3 is a representative plot of signal power or signal strength 30, for example, in dB, of the EM signal received by antenna system 12 in FIG. 2 over range 22. For the sake of ease of illustration, FIG. 3 illustrates only one signal or trace 32 for a single pulse for what typically comprises a plurality of traces of reflected pulses 34 received by antenna system 12 for the plurality of pulses, or burst, transmitted by antenna system 12. As shown in FIG. 3, the power 30 of trace 32 increases as indicated at for example, at a range corresponding to the target 18.

(22) FIG. 4 is representative plot of the signals corresponding to range bins 40, for example, in dB, of the EM signal trace 32 shown in FIG. 3 over range 22. As known in the art, the use of range bins 40 is a procedure by which radar signals are electronically isolated into range intervals 28 (FIG. 1) and into which the respective received signal 32 can be summed over the range of pulses 34 received (as discussed below with respect to FIG. 8). The sums accumulated in the range bins 40 may be analyzed, for example, compared to a given value (such as, the average noise value) and the presence of a target 18 determined. For example, if the power of the received signal increases over successive pulses 34, the use of range bins can assist in differentiating targets from noise. In addition, when antenna system 12 identifies received signals in adjacent range bins 40, it is more likely that a target 18 is present at the range of the range bins. In aspects of the present invention, the summation of data in the range bins 40 shown in FIG. 4 are used to calculate weighting factors used in minimizing clutter.

(23) FIG. 5 is a schematic block diagram of a signal processing antenna system 40 that may be used for antenna system 12 shown in FIGS. 1 and 2 according to aspects of the present invention. As shown in FIG. 5, antenna system 40 may include an antenna/aperture 42 adapted to transmit one or more signals 44 and a receiver 41 adapted to receive reflected signals 45 (for example, a main or upper beam) and 46 (for example, an aux or lower beam) from one or more targets and/or clutter, for example, one or more targets 18 shown in FIG. 2, as described herein. Antenna/aperture 42 typically also includes a beamformer, as known in the art.

(24) As known in the art, receiver 41 of antenna/aperture 42 may be adapted to receive reflected beams 45 and 46 and through, for example, comparison of the two received beams 45 and 46 with the transmitted beam 44 generate output data streams 47 and 49 corresponding to the received beams 45 and 46, respectively. As is typical in the art, receiver 41 may typically include at least one analog-to-digital (A/D) converter for converting received analog return signals 45 and 46 to digital form. Accordingly, output data streams 47 and 49 may typically comprise digital signals containing an amplitude and a phase associated with received beams 45 and 46, as expressed by in-phase (I) and quadrature-phase (Q) signals, commonly referred to as I/Q data streams. Since the I/Q data streams 47 and 49 may be expressed as complex numbers, data streams 47 and 49 may typically be referred to as complex data streams 47 and 49.

(25) As is typical in the art, antenna system 40 may include two or more pulse compressor modules 48 and 50 (for example, a main beam pulse compressor module and an aux beam pulse compressor module) associated with the two or more received reflected signals 45 and 46; one or more spatial nulling modules 52 according to aspects of the present invention; one or more Doppler nulling modules (or Doppler Filter Banks (DFB)) 54; and one or more detection processing modules 56. According to this aspect of the invention, the nulling of spatial nulling module 52 is practiced before the nulling of Doppler nulling module 54, therefore, in one aspect, the antenna system 40 shown in FIG. 5 comprises a Pre-Doppler spatial nulling system.

(26) Pulse compression modules 48 and 50 receive data streams 47 and 49, from antenna/aperture 42 and output modulated upper and lower data streams 51 and 53, respectively. Again, according to aspects of the invention, data streams 51 and 53 may be complex, I/Q data streams. In addition, data streams 51 and 53 may be correlated with respect to range bin n, for example, where the data associated with a range bin n for data stream 51 (for example, the main data stream) is related to the data associated with the same range bin n for data stream 53 (for example, the aux data stream). Data streams 51 and 53 are then processed by spatial nulling module 52, according to aspects of the inventionthe details of the spatial nulling module 52 are disclosed and illustrated with respect to FIG. 6 below.

(27) After spatial nulling 52 the spatially-nulled data stream 55 may then be forwarded to Doppler filter module 54 where spatially-nulled data stream 55 may be Doppler-nulled to produce spatially-nulled, Doppler-nulled data stream 57. As shown in FIG. 5, spatially-nulled, Doppler-nulled data stream 57 may then be forwarded to detection processing 56 to, for example, be used to more accurately locate targets detected by antenna system 40.

(28) FIG. 6 is a detailed flow diagram with annotation for the spatial nulling practiced in spatial nulling module 52 shown in FIG. 5 according to an aspect of the present invention. As shown in FIG. 6, the spatial nulling module 52 (represented by the components outlined by a dash line in FIG. 6) receives data stream 51 from compressor module 48 and data stream 53 from compressor module 50. As noted above, data streams 51 and 53 may comprise complex I/Q data streams.

(29) According to aspects of the invention, in addition to other processes, pulse compressor modules 48 and 50 may also sum or add the received pulse samples by range bin, for example, in each range bin for each received beam. In one aspect, these summations may comprise tapered summations, where a taper is an amplitude profile multiplicatively imposed onto the data. In one aspect, the summation of pulses samples per range bin may typically yield two or more coherent, complex (that is, I/Q) data streams. In one aspect, the complex data streams may comprise complex high clutter-to-noise-ratio (CNR) data streams that are representative of the net clutter signals in the received signals (or signal channels).

(30) The sum of the main data stream data may be represented by .sub.M; the sum of the aux data stream data may be represented by .sub.A. This summing is represented schematically in FIGS. 7 and 8. In one aspect of the invention, the sum of the main data stream data represented by, .sub.M, contains the same data or information as the zero frequency Doppler bin (also know as, the zero Doppler bin).

(31) FIG. 7 is a duplication of FIG. 4, and shows a representative plot of the signals (in dB) received corresponding to range bins, N.sub.R, 40 for pulses, N.sub.P, 34 in a burst shown in FIG. 3. FIG. 8 shows a representative plot of the sum 66 (.sub.M or .sub.A) (in dB) of the signals received in each range bin 40, for each pulse, N.sub.P, in the burst corresponding to range bins, N.sub.R, 40 for pulses, N.sub.P, 34 shown in FIG. 7. As shown in FIG. 8, in this example plot, which is presented for illustration only and does not limit the scope of the present invention, the sum (.sub.M or .sub.A) of the data streams 51 or 53 per pulse result in a peak 68 corresponding to the location of target 18 shown in FIG. 2. As also illustrated in FIG. 8, the sums of the data outside peak 62 may be substantially low, as illustrated by peak 70 in FIG. 8, or non-existent, according to aspects of the invention, and thus may not appear in FIG. 8. In other aspects of the invention, other peaks (such as, peaks similar to peak 66 or 70) may also appear in a plot similar to the plot shown in FIG. 8, which may correspond to targets, clutter, or noise in data streams 51 and 53.

(32) In one aspect, each summation (.sub.M or .sub.A) may represent the output of a Doppler filter centered on a specific frequency, for example, a presumed center of the clutter's power spectral density spectrum, such as, at 0 Hz. In one aspect, the taper used in the summation, for example, Blackman weights, may be selected to produce a filter response that has low side lobes at non-desired Doppler frequencies. Consequently, in one aspect, the summations, for example, as shown in FIG. 8, may have negligible contamination from undesired sources.

(33) In the aspect of the invention shown in FIG. 6, data stream 51 and 53 are also represented by representative matrices or grids or data cubes 58 and 60 to facilitate disclosure of aspects of the invention. Matrices 58 and 60 represent the 2-dimensional nature of data streams 51 and 53. For example, matrices 58 and 60 may comprise summations (.sub.M or .sub.A) of the pulse data for the number of pulses emitted and received, N.sub.P, by receiver 41 (see FIG. 5) for each range bin n for the number of range bins, N.sub.B, generated by receiver 41. In other words, according to one aspect of the invention, data streams 51 and 53 in FIG. 6 may comprise two-dimensional matrices representing the summations (.sub.M or .sub.A) of the amplitudes and phase of the received signals 43 and 45 for each pulse received for each range bin. According to aspects of the invention, at least some of the data streams 51 or 53 contain data related to clutter, as disclosed herein.

(34) As shown in FIG. 6, according to aspects of the present invention, spatial nulling module 52 receives each of the data streams 51 and 53 for each of range bins N.sub.B for each of the pulses N.sub.P in modules 62 and 64, respectively. After receipt, modules 62 and 64 may sum the data streams 51 and 53 for each of range bins N.sub.B for each of the pulses N.sub.P to produce the respective summed data streams .sub.M and/or .sub.A. As shown in FIG. 6, in one aspect of the invention, decision or comparison module 72 compares (for example, gates) the sum of the data stream (.sub.M or .sub.A) at each range bin n of one of the data streams, for example, the sum of the data stream for the main beam, .sub.M, to a predetermined threshold value, T.sub.1, for the sum. In one aspect, the predetermined threshold value, T.sub.1, may be derived from the number of pulses N.sub.P and an estimate of the radar system's phase-noise-induced and pulse-to-pulse jitter Doppler-domain clutter cancellation limit.

(35) According to one aspect of the invention, should the sum of the data for the range bin n be less than (<) the threshold value T.sub.1 (that is, not greater than or equal to T.sub.1), thus indicating that no nulling of this range bin may be necessary, a weighting factor or an adaptive weight, W, is set to zero, as indicated at 74 in FIG. 6. As will be disclosed further below, setting W to zero effectively eliminates the nulling of the data in range cell n, and, for example, no further processing of the nulling of the range cell need be practiced. An illustration of a representative threshold value T.sub.1 with respect to the pulse sums .sub.M or .sub.A for each range bin is shown in FIG. 8.

(36) According to another aspect of the invention, should the sum of the data (.sub.M) for the range bin n be greater than or equal to () the threshold value T.sub.1, thus indicating that nulling of this range bin may be desired, for example, to reduce the single pulse clutter to noise-ratio (CNR), according to an aspect of the invention, an adaptive weighting factor, W, is computed. In one aspect, W may be computed as a function of the summations .sub.M and .sub.A. Specifically, in one aspect, W may be computed as a quotient including the summations .sub.M and .sub.A for each range bin n. In one aspect, the weighting factor, W, is defined as the quotient that appears in Equation 1.
W=.sub.M/(.sub.A+)Equation 1.
In Equation 1, W is the adaptive weighting factor; .sub.M is the sum of one stream of data, for example, the main data stream 51, for range bin n from module 62; .sub.A is the sum of another stream of data, for example, the aux data stream 53, for range bin n from module 64; and 8 is a predetermined constant selected based upon the parameters of the system 40 (see FIG. 5). In one aspect, may be referred to as the null depth control constant.

(37) In one aspect of the invention, the weighting factor W may be applied to one of the aux pulse data stream 53 to generate a spatially-nulled data stream 78 for each range bin n. In one aspect, applying weighting factor W to a data stream may comprise multiplying 80 the data stream, for example, data stream 53, by weighting factor W to yield a spatially-nulled data stream 78.

(38) In another aspect, the weighting factor W may be compared to (for example, gated) with a second threshold value T.sub.2, as indicated at comparison module 82. In one aspect, the absolute value of weighting factor or the magnitude of W (|W|) may be compared with threshold value T.sub.2. According to an aspect of the invention, the comparing of the magnitude of W to T.sub.2 may be practiced to determine whether any unacceptable distortion of the 2-way beam has occurred. According to aspects of the invention, the threshold value T.sub.2 may range from about 0.05 to about 3.0, but is typically in the range of about 0.2 to about 0.5 for example, on the order of about 0.3.

(39) According to this aspect of the invention, should the absolute value of weighting factor W for the range bin n be less than or equal to () the second threshold value T.sub.2, thus indicating that weighting factor W is of sufficient magnitude to provide a valid spatial nulling, weighting factor W may be applied to one of data stream 51 or 53, but typically 53, to generate a spatially-nulled data stream 78 for each range bin n.

(40) In another aspect, should the absolute value of weighting factor W be greater than second threshold value T.sub.2, that is, not less than or equal to () second threshold value T.sub.2, an adjustment may be made to weighting factor W.

(41) In one aspect, the weighting factor W may be reduced. In one aspect, the weighting factor may be reduced to limit both the magnitude of the adaptive weight W, and/or to reduce the resulting spatial cancellation ratio (CR) to no more than necessary to supplement Doppler nulling. In addition, in one aspect, limiting the weight magnitude (|W|) and the CR may minimize the distortion of the effective nulled beam pattern and, consequently, may minimize gain loss in potential target directions.

(42) In one aspect, as indicated at module 84, weighting factor W may be reduced by, first, increasing the value of predetermined constant by a predetermined constant , as indicated by Equation 2.
=+Equation 2.
In one aspect, may be referred to as an increment to . In another aspect, the magnitude of W may be limited by the value of . For example, the value of may affect the receive gain advantage of the auxiliary receive beam, for instance, the aux beam, in the direction of potential surface clutter and be more compliant.

(43) As shown in FIG. 6, with an increase in the value of computed by Equation 2, Equation 1 might be re-evaluated, as indicated in module 76, for a new value of weighting factor W. In one aspect, this new weighting factor W may again be compared to second threshold value T.sub.2, as indicated at 82, and, if less than or equal to T.sub.2, applied to one of the data streams. In one aspect, when the new value of W is not less than or equal to T.sub.2, the value of W may be further decreased, for example, by increasing in module 84, and re-evaluating W in module 76 and, again, comparing W in module 82. This process may continue until weighting factor W complies with the gating of module 82.

(44) However, according to some aspects of the invention, for example, in some outlier cases, a particular clutter distribution may result in a weight magnitude (|W|) on the order of one or greater, for instance, despite the limiting effects of the value of . For example, it is envisioned that 1) a single dominant clutter source at an elevation angle that produces a larger response in the main beam than in the auxiliary beam may increase the weight magnitude (|W|) desired, or 2) multiple clutter sources that may coherently sum to a larger response in the main beam than in the auxiliary beam may increase the weight magnitude (|W|) desired. In these and other cases, may be increased to bring the weight magnitude down below the threshold T.sub.2 and processing may proceed with some cancellation of clutter sacrificed to control distortion of the nulled beam shape.

(45) As noted above, should the value of .sub.M not be greater than or equal to threshold value T1, the weighting factor for range bin n may be equated to zero (W=0) and applied to one of the data streams 51 or 53, typically, aux data stream 53. As indicated in FIG. 6 applying W=0 to data stream 53 may comprise multiplying data stream 53 by zero to generate nulled data stream 78 with a value of zero, and thus eliminating at least some spatial nulling of the data stream for the specific range bin n under consideration.

(46) According to aspects of the invention, regardless of the weighting factor W, weighting factor W may be applied to, for example, multiplied into, one of the data streams 51 or 53, for example, to aux data stream 53, to provide a weighted data stream 78. As shown in FIG. 6, in one aspect, weighted data stream 78 may then be applied to one of the data streams 51 or 53, for example, main data stream 51, to produce a spatially nulled data stream 86. In one aspect, applied to may comprise adding weighted data stream 78 to one of the data streams 51 or 53, for example, data stream 51, as indicated at addition module 84, to produce a spatially-nulled data stream 86.

(47) In one aspect, spatially-nulled data stream 86 comprises a pulse-by-pulse, spatially-nulled data stream that is spatially nulled for some or a plurality of range bins, for example, each or all range bins N.sub.R.

(48) According to aspects of the invention, for the aspect shown in FIG. 6 or for any aspect disclosed herein, regardless of the value of weighting factor W, a single, unique weighting factor W may be provided per range bin. For example, in one aspect, a single, unique weighting factor W may be provided for each range bin requiring spatial nulling. In one aspect, a single, unique weighting factor W may be provided for all the Doppler bins at each range bin, N.sub.P.

(49) In the aspect of the invention shown in FIG. 6, spatially-nulled data stream 86 may be represented by a representative matrix or grid or data cube 88 to facilitate disclosure of aspects of the invention. Matrix 88 represents the 2-dimensional nature of data stream 86, for example, comprising pulse data for the number of pulses received, N.sub.P, for each range bin n for the number of range bins, N.sub.B, generated by spatial nulling module 52. In one aspect, matrix or data cube 88 represents the spatially nulled range-pulse data cube for each pulse burst having pulses N.sub.P.

(50) As shown in FIGS. 5 and 6, after the spatially-nulled data stream 86 is generated by spatial nulling module 52, the spatially-nulled data stream 86 may be Doppler-nulled in Doppler-nulling module 54 to generate a Doppler-nulled, spatially-nulled data stream 57 that may then be forwarded to detection processing module 56. Doppler nulling in Doppler-nulling module 54 may be practiced by any Doppler nulling procedure or algorithm to suppress the influence of, for example, any moving targets on the data stream and produce Doppler-nulled, spatially nulled data stream 57. For example, in one aspect, Doppler nulling in module 54 may comprise zero-hertz Doppler nulling or its equivalent.

(51) In the aspect of the invention shown in FIG. 6, Doppler-filtered, spatially-nulled data stream 57 may be represented by a representative matrix or grid or data cube 90 to facilitate disclosure of aspects of the invention. Matrix 90 represents the 2-dimensional nature of data stream 57, for example, comprising Doppler-filtered, spatially-nulled pulse data for each range bin n for the number of range bins, N.sub.R, generated by the number of Doppler filters, N.sub.D, of Doppler nulling module 54. As known in the art, the Doppler Filter Bank (DFB) of Doppler filter module 54 may yield a matrix 90 of range bins, N.sub.R, by Doppler filters, N.sub.D.

(52) According to aspects of the invention, detection processing in detection processing module 56 may comprise any detecting process or algorithm. In one aspect, detection processing may be practiced by conventional means, such as by thresh-holding via CFAR (Constant False Alarm Rate) processing, and like methods.

(53) As illustrated in FIGS. 5 and 6, according to aspects of the invention, each range cell, that is, each range-azimuth cell may be independently nulled to supplement Doppler nulling, if needed.

(54) As known in the art, since Doppler nulling and spatial nulling both comprise linear coherent processing, Doppler nulling and spatial nulling are commutative. That is, though the spatial nulling disclosed with respect to FIGS. 5 and 6 above followed the Doppler filtering, spatial nulling disclosed with respect to FIGS. 5 and 6 may commute with Doppler filtering and precede Doppler filtering. In one aspect, though spatial nulling may precede or follow Doppler filtering, pulse compression may be performed first. In some aspects, pulse compression must be performed prior to spatial nulling and/or Doppler filtering.

(55) FIG. 9 is a schematic block diagram of a signal processing antenna system 100 that may be used for antenna system 12 shown in FIGS. 1 and 2 according to another aspect of the present invention. As shown in FIG. 9, in a fashion similar to antenna system 40 shown in FIG. 5, antenna system 100 may include an antenna/aperture 102 adapted to transmit one or more signals 103 and a receiver 101 adapted to receive reflected signals 104 (for example, a main or upper beam) and 106 (for example, an aux or lower beam) from one or more targets and/or clutter, for example, one or more targets 18 shown in FIG. 1, as described herein.

(56) As known in the art, receiver 101 of antenna/aperture 102 may be similar to and have the functionality of receiver 41 of antenna/aperture 42 shown in FIG. 5. Accordingly, receiver 101 may be adapted to receive reflected beams 104 and 106 from transmitted beam 103, respectively, and through comparison of the two received beams 103 104 and 106 with the transmitted beam 103 generate output data streams 107 and 109 corresponding to the received beams 104 and 106, respectively. As is typical in the art, output data streams 107 and 109 may comprise digital I/Q data streams or complex data streams.

(57) As is typical in the art, antenna system 100 may include two or more pulse compressor modules 108 and 110 (for example, a main beam pulse compressor module and an aux beam pulse compressor module) associated with the two or more received reflected signals 104 and 106; one or more Doppler nulling modules (or Doppler Filter Banks (DFB)) 112, 114; one or more spatial nulling modules 116 according to aspects of the present invention; and one or more detection processing modules 118.

(58) According to this aspect, pulse compression modules 108 and 110 may be similar in functionality to pulse compressors 48 and 50 shown in FIG. 5. Accordingly, pulse compression modules 108 and 110 may receive data streams 107 and 109, from antenna/aperture 102 and output modulated upper and lower data streams 111 and 113. According to aspects of the invention data streams 111 and 113 may be I/Q data streams. Data streams 111 and 113 are then processed by Doppler filter modules 112 and 114, respectively, to produce Doppler-nulled data stream 115 and 117.

(59) Doppler nulling in Doppler-nulling modules 112 and 114 may be practiced by any Doppler nulling process or algorithm to suppress the influence of, for example, any moving targets in the data stream, and produce Doppler-nulled data streams 115 and 117, respectfully. For example, in one aspect, Doppler nulling in modules 112 and 114 may comprise zero-hertz Doppler nulling or its equivalent.

(60) As shown in FIG. 9, data streams 111 and 113 and Doppler-nulled data streams 115 and 117 may then be forwarded to spatial nulling module 116, according to aspects of the invention to yield a Doppler-nulled and spatially-nulled data stream 119. The details of the spatial nulling module 116 are illustrated and described with respect to FIG. 10 below. The Doppler-nulled and spatially-nulled data stream 119 from spatial nulling module 116 may then be forwarded to detection processing 118 to, for example, be used to more accurately locate targets detected by antenna system 100.

(61) According to this aspect of the invention, the nulling of spatial nulling module 116 may be practiced after the nulling of Doppler nulling modules 112 and 114. Therefore, in one aspect, the antenna system 100 shown in FIG. 9 comprises a Post-Doppler spatial nulling system.

(62) FIG. 10 is a detailed flow diagram with annotation for the spatial nulling practiced in spatial nulling module 116 shown in FIG. 9 according to an aspect of the present invention. In order to facilitate disclosure of aspects of the invention, compressors 108 and 110, DFBs 112 and 114, and detection processing module 118 shown in FIG. 9 also appear in FIG. 10.

(63) As shown in FIG. 10, the spatial nulling module 116 (represented by the components outlined by a dash line in FIG. 10) receives data stream 111 from compressor module 108, data stream 113 from compressor module 110, Doppler-nulled data stream 115 from Doppler nulling module 112, and Doppler-nulled data stream 117 from Doppler nulling module 114. As noted above, data streams 111, 113, 115, and 117 may comprise complex I/Q data streams.

(64) According to aspects of the invention, in addition to other processes, pulse compressor modules 108 and 110 may also sum or add the received pulse samples by range cell, for example, in each range cell for each received beam. In one aspect, these summations may comprise tapered summations. In one aspect, the summation of pulses samples per range cell may typically yield two or more complex (that is, I/Q) data streams. In one aspect, the complex data streams may comprise complex high clutter-to-noise-ratio (CNR) data streams that are representative of the net clutter signals in the received signals (or signal channels).

(65) In a manner similar to the process disclosed in FIG. 6, the sum of the main data stream data may be represented by .sub.M in FIG. 10, and the sum of the aux data stream data may be represented by .sub.A. This summing is represented schematically in FIGS. 7 and 8, as discussed above.

(66) In the aspect of the invention shown in FIG. 10, as in other aspects disclosed herein, data stream 111, 113, 1115, and 117 are also represented by representative matrices or grids or data cubes 120, 122, 124, and 126, respectively, to facilitate disclosure of aspects of the invention. Matrices 120, 122, 124, and 126, represent the possible 2-dimensional nature of data streams 111, 113, 115, and 117. For example, data streams 111 and 113 may comprise pulse data for the number of pulses emitted and received, N.sub.P, by receiver 101 (see FIG. 9) for each range bin for the number of range bins, N.sub.B, generated by receiver 101. In other words, according to one aspect of the invention, data streams 111 and 113 in FIG. 10 may comprise two-dimensional matrices representing the summations (.sub.M or .sub.A) of the amplitudes or power of the received signals 104 and 106 for each pulse received for each range bin. Matrices 124 and 126 represent the 2-dimensional nature of data stream 115 and 117, for example, comprising Doppler-filtered pulse data for each range bin n for the number of range bins, N.sub.R, generated by the number of Doppler filters, N.sub.D, of Doppler nulling modules 112 and 114. As known in the art, the Doppler Filter Banks (DFB) of Doppler filter modules 112 and 114 may yield matrices 124 and 126 of range bins, N.sub.R, by Doppler filters, N.sub.D. According to aspects of the invention, at least some of the data streams 111, 113, 115, and 117 contain data related to clutter, as disclosed herein.

(67) In one aspect of the invention, the data processing practiced in spatial nulling module 116 may be similar, if not substantially identical, to the data processing described for nulling module 52 shown in and described with respect to FIG. 6. In one aspect, the distinction between the processing shown in and described with respect to FIG. 10 differs from that shown in and described with respect to FIG. 6 in the relative timing of the spatial nulling with respect to the Doppler nulling. For example, in FIG. 6, the spatial nulling is practiced before the Doppler nulling; in FIG. 10, the spatial nulling is practiced after the Doppler nulling.

(68) As shown in FIG. 10, according to aspects of the present invention, spatial nulling module 116 receives each of the data streams 111 and 113 for each of range bins N.sub.B for each of the pulses N.sub.P for each of range bins N.sub.B in modules 132 and 134, respectively. After receipt, modules 132 and 134 may sum the data streams 111 and 113 for each of range bins N.sub.B for each of the pulses N.sub.P to produce the respective summed data streams .sub.M and/or .sub.A. In one aspect of the invention, decision or comparison module 142 compares (for example, gates) the sum of the data stream at each range bin n of one of the data streams, for example, the sum of the data stream for the main beam, .sub.M, to a predetermined threshold value, T.sub.1, for the sum. According to this aspect of the invention, should the sum of the data for the range bin n be less than (<) the threshold value T.sub.1, thus indicating that no nulling of this range bin is necessary, a weighting factor or an adaptive weight, W, is set to zero, as indicated at 144, effectively eliminating the nulling of the data in range cell n under consideration. An illustration of threshold T.sub.1 with respect to the pulse sums .sub.M or .sub.A for each range bin is shown in FIG. 8. The value of threshold T.sub.1 may be substantially the same as the value of threshold value T.sub.1 shown in and described with respect to FIG. 6.

(69) As shown at 142 in FIG. 10, according to another aspect of the invention, should the sum of the data for the range bin n be greater than or equal to () the threshold value T.sub.1, thus indicating that nulling of this range bin is desired, a weighting factor, W, is computed as a function of the summations .sub.M and .sub.A, for example, a function of a quotient having the sums .sub.M and .sub.A for each range bin n. In one aspect, the weighting factor, Equation 1 shown in and described with respect to FIG. 6, and repeated here defines W.
W=.sub.M/(.sub.A+)Equation 1.
The variables in Equation 1 and the range of values may be the same as those described with respect to FIG. 6.

(70) As shown in FIG. 10, in one aspect of the invention, the weighting factor W may be applied to the aux pulse data stream 117 to generate a spatially-nulled data stream 148 for each range bin n. In one aspect, applying weighting factor W to a data stream may comprise multiplying 150 the data stream, for example, data stream 117 by weighting factor W to yield a spatially-nulled data stream 148.

(71) In another aspect, the weighting factor W in FIG. 10 may be compared to (for example, gated) with a second threshold value T.sub.2, as indicated at comparison module 152. In one aspect, the absolute value weighting factor W (|W|) may be compared with threshold value T.sub.2. According to this aspect of the invention, should the absolute value of weighting factor W for the range bin n be less than or equal to () the second threshold value T.sub.2, thus indicating that weighting factor W is of sufficient magnitude to provide a valid spatial nulling, weighting factor W may be applied to one of the Doppler-nulled data stream 115 or 117 to generate a spatially-nulled data stream 148 for each range bin n. The value of threshold T.sub.2 may be substantially the same as the value of threshold value T.sub.2 shown in and described with respect to FIG. 6.

(72) In another aspect, should the absolute value of weightings factor W be greater than second threshold value T.sub.2, that is, not less than or equal second threshold value T.sub.2, an adjustment may be made to weighting factor W. In one aspect, the weighting factor W is reduced. For example, in one aspect, as indicated at module 154, weighting factor W is reduced by, first, increasing the value of predetermined constant by a predetermined constant , for example, in a fashion substantially identical to that illustrated in FIG. 6, for example, by Equation 2, which is reproduced here.
=+Equation 2.
In one aspect, value of may be substantially identical to the values of disclosed with respect to FIG. 6.

(73) As shown in FIG. 10, with the increase in the value of computed by Equation 2, Equation 1 is re-evaluated, as indicated in module 146, for a new value of weighting factor W. In one aspect, this new weighting factor W may again be compared with second threshold value T.sub.2, and, if less than or equal to T.sub.2, applied to one of the data streams. In one aspect, when the new value of W is not less than or equal to T.sub.2, the value of W may be decreased, for example, by again increasing in module 154, and re-evaluating W in module 146 and comparing W in module 152. This process may continue until weighting factor W complies with the gating of module 152.

(74) As noted above, should the value of .sub.M not be greater than or equal to threshold value T1, the weighting factor for range bin n is equated to zero (W=0) as shown at 144 and applied to one of the Doppler-nulled data streams 115 or 117, typically, aux data stream 117. As indicated in FIG. 10, applying W=0 to data stream 117 may comprise multiplying data stream 117 by zero, thus eliminating any spatial nulling of the data for the specific range bin n under consideration.

(75) According to aspects of the invention, regardless of the weighting factor W, weighting factor W is then applied to, for example, multiplied into as indicated at 150, one of the Doppler-nulled data streams 115 or 117, for example, to aux data stream 117, to provide a weighted data stream 148. As shown in FIG. 10, in one aspect, weighted data stream 148 may then be applied to one of the Doppler-nulled data streams 115 or 117, for example, data stream 115, to produce a Doppler-nulled and spatially-nulled data stream 158. In one aspect, applied to may comprise adding weighted data stream 148 to one of the Doppler-nulled data streams 115 or 117, for example, data stream 115, as indicated at addition module 156, to produce a Doppler-nulled and spatially-nulled data stream 158.

(76) In one aspect, Doppler-nulled and spatially-nulled data stream 158 may comprise a pulse-by-pulse, Doppler-nulled, spatially-nulled data stream that is Doppler-nulled and spatially-nulled for some or a plurality of range bins, for example, for each or all range bins N.sub.R.

(77) In the aspect of the invention shown in FIG. 10, Doppler-nulled and spatially-nulled data stream 158 may be represented by a representative matrix or grid or data cube 160 to facilitate disclosure of aspects of the invention. Matrix 160 represents the possible 2-dimensional nature of Doppler-nulled and spatially-nulled data stream 158, for example, comprising pulse data for the number of pulses received, N.sub.P, for each range bin n for the number of range bins, N.sub.B, generated by spatial nulling module 116. In one aspect, matrix or data cube 160 represents the spatially nulled range-pulse data cube for each pulse burst having pulses N.sub.P.

(78) As shown in FIGS. 9 and 10, after Doppler-nulled and spatially-nulled data stream 158 is generated by spatial nulling module 116, the Doppler-nulled and spatially-nulled data stream 158 may be forwarded to detection processing module 118. According to aspects of the invention, detection processing in detection processing module 118 may comprise any detecting process or algorithm. In one aspect, detection processing may comprise any one of the detection processes disclosed above with respect to detection processing module 56 in FIG. 6.

(79) As illustrated in FIGS. 9 and 10, according to aspects of the invention, each range bin, that is, each range-azimuth bin may be independently nulled to supplement Doppler nulling, if needed.

(80) FIG. 11 is a representative plot 200 of a typical unweighted transmitted radar beam represented by curve 202 and two weighted receive beams 204 and 206, according to one aspect of the invention. The two weighted receive beams 204 and 206 may be a first receive beam represented by curve 204 corresponding to a main beam, and a second received beam represented by curve 206, corresponding to an aux beam, according to aspects of the invention. The second received beam or aux beam 206 is a 30-dB-weighted receive beam suitable for use as the auxiliary beam, according to aspects of the invention.

(81) The abscissa 208 of plot 200 represents the elevation angle of the beams (in degrees) and the ordinate 210 represents the antenna gain (in decibels (dB)). As indicated in FIG. 11, the elevation angles of the transmitted main beam 202 and of first received main beam 204 are in the range of about 2-3 degrees, and the elevation of second received aux beam 206 is about 0 degrees. Accordingly, in this example, the transmitted main beam 202 and the first received main beam 204 at about the same elevation form a two-way pair of transmitted and received beams for target detection. In this example, the transmitting/receiving antenna comprises an array-type antenna having 32 rows of antenna elements and 32 elements per row.

(82) In the representative example shown in FIG. 11, the ground clutter is assumed to be located at about 0 degrees elevation. Since the two-way pair of beams 202 and 204 are elevated from the assumed location of ground clutter, the first receive beam 204 has reduced gain on the ground clutter at 0 degreesthat is, since there was little or no reflectance from the clutter. In contrast, the second receive beam 206 directed at about 0 degrees exhibits a relatively high gain, that is, indicating relatively high reflectance from ground clutter as indicated at 212, for example, a maximum gain (36 dB).

(83) FIG. 12 is a representative plot 220 that illustrates the spatial nulling of first received signal 204 (shown in FIG. 11), for example, the main channel, before and after spatial nulling according to aspects of the invention. Again, the abscissa 222 of plot 220 represents the elevation angle of the beams (in degrees) and the ordinate 224 represents the antenna gain (in decibels (dB)).

(84) As shown in FIG. 12, the first received signal represented by curve 204 (also shown in FIG. 11) comprises a weighted (for example, 30 dB Taylor weighted) main channel receive beam before spatial nulling. Curve 226 represents the weighted (for example, 30 dB Taylor weighted) receive beam after spatial nulling on clutter (at 0 degrees elevation), according to aspects of the invention.

(85) As shown in FIG. 12, in this example, where the initial quiescent received beam 204 indicates an antenna gain of about 27 dB at 0 degrees elevation (that is, an indication of the clutter located at 0 degrees elevation), the spatially-nulled beam 226 indicates an antenna gain of about 0 dB at 0 degrees elevation. That is, according to aspects of the present invention, the clutter located about 0 degrees elevation (as detected, for example, by the gain associated with first or main receive beam 204 (see FIG. 11) has been suppressed in second received or aux beam 206 (see FIGS. 11 and 12) to yield a clutter-suppressed signal 226.

(86) According to aspects of the invention, the clutter-suppressed signal 226 (along with other similar clutter suppressed signals) can be used to more accurately detect targets 18, for example, with target detection processing module 118 shown in FIG. 10. Moreover, according aspects of the invention, this clutter suppression can be practiced by range bin, for example, at each range bin under consideration by the system or method.

(87) As also indicated in FIG. 12, the perturbation to the pattern is also apparent, where a perturbation is a deviation from its unadapted gain pattern. That is, in one aspect, the amount of the clutter reduction or suppression can be controllable as previously described. Reducing the spatial nulling reduces the pattern perturbation and helps to restore loss of target gain.

(88) According to aspects of the invention disclosed herein, a single, unique weight W may be provided and applied to data streams based upon range bin, that is, based upon individual range bins. It is also envisioned that adaptive weight W may also be computed based upon more than one range bin, for example, two or more range bins. In a further refinement, it envisioned that an adaptive weight W may also be computed from a range of range bins (for example, a contiguous block of range bins) or two or more ranges of range bins (for example, two or more split windows or non-contiguous ranges of range bins). The one or more adaptive weights W computed from two or more range bins may then be applied (for example, multiplied) to one or more range bins, for example, to one or more range bins associated with or not associated with the range bins from which the one or more weights W are computed.

(89) It is believed that providing one or more adaptive weights W based upon one or more range bins, as disclosed herein, may provide advantages over the use of weights based upon single range bins. The potential advantages include: 1) fewer weights W need be computed because the resulting average weights can be applied to multiple range cells; and 2) signals received from targets may inherently be excluded from the weight computations when the weights are based upon one set of range cells and applied to another. However, it is also recognized that there may be disadvantages to providing one or more adaptive weights based upon one or more range bins. These disadvantages may include 1) more computations per weight W; and 2) reduced spatial cancellation performance because the clutter signal can fluctuate by tens of dB from cell to cell while the weight is tuned to an average. However, it is believed that these perceived disadvantages may be overcome by the actual advantages of aspects of the invention.

(90) As described herein, it is envisioned that the spatial nulling according to aspects of the invention, as indicated by FIGS. 11 and 12, may augment or be added to the Doppler-based reduction or suppression of clutter and provide improvements over the existing art. For example, aspects of the invention are believed to provide unique advantages over the prior art. Some of these advantages include: 1) simplicity, because each weight, W, may be simply the function of the ratio of two complex numbers; 2) superior spatial cancellation performance, because each weight may be substantially tuned to the particular clutter scatterer distribution and signal it is intended to suppress; 3) target signals may be inherently excluded from the weight computations due to Doppler offsets combined with heavy Doppler weighting; and 4) any target signals at very low Doppler frequency may generally be much lower than the terrain clutter.

(91) As disclosed herein, aspects of the invention may be uniquely adapted for use in radar systems, for example, in pulse-Doppler radar systems. However, it is envisioned that the systems and methods of the present invention may also be used to minimize or eliminate clutter in any form of remote sensing system, including, but not limited to, sonar systems, navigation systems, and guidance systems.

(92) The embodiments disclosed herein are described by way of example only, and other embodiments for implementing the methods and systems described herein may be contemplated by one of skill in the pertinent art without departing from the intended scope of this disclosure.

(93) In additions, the methods and the systems described herein may include conventional data processing, storage, and user interfaces. For example, the method and systems disclosed herein may include memory devices containing data, which may include instructions, the instructions when executed by a processor or multiple processors, cause the steps of the methods for performing the signal processing and weight generation as set forth herein. Software may be embodied in a non-transitory machine-readable medium upon which software instructions may be stored, the stored instructions when executed by a processor may cause the processor to perform the steps of the methods described herein. Any suitable machine-readable storage medium may be used, including, but not limited to, magnetic or optical disks, for example, CD-ROM, DVD-ROM, floppy disks, and the like. Other storage media also fall within the intended scope of this disclosure, for example, dynamic random access memory (DRAM), random access memory (RAM), read-only memory (ROM), or flash memory may also be used.

(94) While various embodiments have been described above, it should be understood that these embodiments and their many aspects have been presented by way of example, and not limitation. It will be apparent to persons skilled in the relevant art that various changes in form and detail can be made therein without departing from the spirit and scope of the invention. Moreover, it is to be understood that the various embodiments of the invention, although different, are not necessarily mutually exclusive. Furthermore, a particular feature, structure, or characteristic described herein in connection with one embodiment may be implemented within other embodiments without departing from the scope of the invention. In addition, it is to be understood that the location or arrangement of individual elements within each disclosed embodiment may be modified without departing from the scope of the invention. The detailed description presented herein, therefore, is not to be taken in a limiting sense, and the scope of the present invention is defined only by the appended claims, appropriately interpreted, along with the full range of equivalents to which the claims are entitled.

(95) Although the term at least one may often be used in the specification, claims and drawings, the terms a, an, the, said, etc. also signify at least one or the at least one in the specification, claims and drawings.

(96) While several aspects of the present invention have been described and depicted herein, alternative aspects may be affected by those skilled in the art to accomplish the same objectives. Accordingly, it is intended by the appended claims to cover all such alternative aspects as fall within the true spirit and scope of the invention.