System and technique for mitigation of clutter in radar
10955524 ยท 2021-03-23
Assignee
Inventors
Cpc classification
G01S7/2923
PHYSICS
G01S13/0209
PHYSICS
G01S13/0218
PHYSICS
G01S13/5248
PHYSICS
G01S7/026
PHYSICS
G01S13/5246
PHYSICS
International classification
G01S7/41
PHYSICS
G01S13/524
PHYSICS
Abstract
Described are a system and technique to mitigate the impacts of clutter in a radar system. The system and technique require only linear co-polarized measurements can be incorporated into the standard radar signal processing chain without slowing down radar performance.
Claims
1. A method for clutter mitigation in a radar system capable of producing detections using RF signals having two orthogonal polarizations, the method comprising: processing sets of radar return pulses in each of a first polarization and a second, associated orthogonal polarization to form processed signals, with each set of radar return pulses including one or more pulses; performing detections on each processed signal to form a detection set, the detections corresponding to potential targets; associating pairs of detections across the first and second polarizations; and monitoring displacement of the associated pairs of detections across the first and second polarizations across a sequence of sets of radar return pulses.
2. The method of claim 1 further comprising identifying associated pairs of detections that move with respect to each other over time.
3. The method of claim 2 further comprising removing from the detection set those associated pairs of detections that move at different rates of speed with respect to each other over time.
4. The method of claim 3 further comprising processing those associated pairs of detections that do not move with respect to each other over time to track an object.
5. The method of claim 3 wherein identifying associated pairs of detections that move with respect to each other over time comprises identifying associated pairs of detections that move apart over time as false alarms due to clutter.
6. A method for clutter mitigation in a radar system capable of producing detections using RF signals having two orthogonal polarizations, the method comprising: processing sets of radar return pulses in each of a first polarization and a second, associated orthogonal polarization to form processed signals, with each set of radar return pulses including one or more pulses; performing detections on each processed signal to form a detection set; associating pairs of detections across the first and second polarizations; and monitoring displacement of the associated pairs of detections across the first and second polarizations across a sequence of sets of radar return pulses, wherein processing sets of radar return pulses in each of a first polarization and a second, associated orthogonal polarization includes averaging sets of radar return pulses in each of the first polarization and the second, associated orthogonal polarization to form averaged signals, and wherein performing detections on each processed signal to form a detection set includes performing detections on each averaged signal to form a detection set.
7. A method for mitigating sea clutter in a marine radar system capable of producing detections using RF signals having two orthogonal polarizations (H-H and (V-V), the method comprising: processing pulses in a CPI in each of HH and VV to form averaged signals; performing CFAR detection on each averaged signal; associating pairs of CFAR detections across H-H and V-V; and monitoring displacement of the associated pairs of CFAR detections across H-H and V-V across a sequence of CPIs.
8. The method of claim 7 further comprising identifying associated pairs of CFAR detections that move apart over time.
9. The method of claim 8 further comprising filtering out of a detection set those associated pairs of CFAR detections that move apart over time.
10. The method of claim 9 identifying those associated pairs of CFAR detections that move apart over time as false alarms due to sea clutter.
11. The method of claim 10 further comprising providing to a radar tracker those associated pairs of CFAR detections that do not move apart over time.
12. The method of claim 6 further comprising identifying associated pairs of detections that move with respect to each other over time.
13. The method of claim 12 further comprising removing from the detection set those associated pairs of detections that move at different rates of speed with respect to each other over time.
14. The method of claim 13 further comprising processing those associated pairs of detections that do not move with respect to each other over time to track an object.
15. The method of claim 13 wherein identifying associated pairs of detections that move with respect to each other over time comprises identifying associated pairs of detections that move apart over time as false alarms due to clutter.
Description
BRIEF DESCRIPTION OF THE DRAWINGS
(1) The foregoing features may be more fully understood from the following description of the drawings in which:
(2)
(3)
(4)
(5)
(6)
(7)
DETAILED DESCRIPTION
(8) Before describing concepts and techniques directed toward systems and techniques for leveraging polarimetric information to mitigate the impacts of clutter (e.g. sea clutter) on radar detections (e.g. CFAR detection) and related devices and systems, some introductory concepts and terminology are explained.
(9) In accordance with the notation used herein, radio frequency (RF) signals which are horizontally polarized on both transmit and receive, are denoted as H-H or HH signals while RF signals which are vertically polarized on transmit and receive are denoted as V-V or VV signals. Thus, so-called H-H radars produce H-H radar returns (or HH radar data or more simply HH returns). Similarly, so-called V-V radars produce V-V radar returns (or VV radar data or more simply VV returns).
(10) Further, it is recognized that the described concepts, systems and techniques find particular use in maritime radar systems operating at a low grazing-angle (LGA). Thus, in one embodiment described herein, polarimetric radar and a characteristic of sea clutter are used to mitigate the impacts of sea clutter in radar systems. One illustrative embodiment of the technique to mitigate the impacts of sea clutter in radar systems is sometimes referred to herein as Polarimetric Co-location Layering (PCL). Thus, reference is sometimes made herein to an illustrative PCL technique and/or to an illustrative PCL system that is specifically directed toward mitigation of sea clutter in a maritime radar.
(11) Accordingly, while some figures herein are used to help illustrate a particular technique that leverages a characteristic of the Doppler spectra of sea clutter and man-made objects to classify detections, the broad concepts, systems and techniques described herein are not so limited.
(12) For example, it should be appreciated that while in the description herein below reference is made to systems operating with linear HH and VV polarized signals, more generally, the concepts, systems and techniques described herein may be used with any orthogonally polarized signals (e.g. in applications not specifically directed toward mitigation of sea clutter). Thus, while the concepts described herein find application with sea clutter, the general concepts, systems and techniques may also find application with mitigation of other types of clutter or with mitigation of undesired signals in general, in which case it may be desirable or even necessary to use orthogonally polarized signals other than HH, VV (e.g. right hand and left hand circularly polarized signals).
(13) Furthermore, it should be appreciated that reference is sometimes made herein to illustrative embodiments in which averaging of pulses takes place. After reading the disclosure provided herein, those of ordinary skill in the art will readily recognize that while some embodiments may utilize coherent averaging of pulses, it is not required. In some embodiments one could simply compute the sum of the magnitudes rather than the magnitude of the sum and still be successful or leverage other processing as described herein.
(14) Also, with respect to the techniques of
(15) Turning now to
(16) Processing then proceeds to processing block 12 in which detections are performed on each processed signal. If the processing of processing block 10 includes averaging, then detections are performed on each averaged signal. Thus, a set of detections (or detection set) is formed in both the first and second orthogonal polarizations. The detections correspond to potential targets or objects of interest.
(17) The detections are then associated across the orthogonal polarizations (e.g. across HH and VV) to form associated detection pairs as shown in processing block 14. This requires a priori knowledge of known differences in responses to orthogonally polarized data from objects giving rise to clutter returns. For example, the knowledge that CFAR detections on sea clutter in HH radar data move at different speeds across coherent processing intervals (CPIs) relative to CFAR detections on sea clutter in VV.
(18) Once the associated detection pairs are formed, then as shown in processing block 16, the displacement of associated detection pairs is monitored across a sequence of CPIs. Associated detection pairs that move relative to each other over time (e.g. move apart over time) are deemed likely to be false alarms due to clutter (e.g. sea clutter) and are filtered out of the detection set that is passed on to a radar tracker, for example.
(19) Thus, as shown in processing block 18, detections moving at different rates of speed within each associated detection pair are eliminated as potential targets. This again requires a priori knowledge of known differences in responses to orthogonally polarized data from objects giving rise to clutter returns. For example, with the knowledge that CFAR detections on sea clutter in HH radar data move at different speeds across coherent processing intervals (CPIs) relative to CFAR detections on sea clutter in VV and that the same is not true for detections on objects of interest, then it is possible to eliminate detections moving at different rates of speed.
(20) Processing then proceeds to processing block 20 where one or more targets is identified.
(21) Referring now to
(22) First consider range-time intensities (RTIs) from one CPI of this data across HH and VV polarizations. The CPI comprises 100 pulses, slow-time decimated by factor m=7, yielding a pulse repetition frequency (PRF) of approximately 893 Hz. The RTIs for HH and VV are illustrated in
(23)
in which:
(24) p corresponds to the RTI pulse index;
(25) n corresponds to the range index; and
(26) subscripts hh and vv indicate the channel whose RTI is being coherently averaged.
(27) It should again, be appreciated that while coherent averaging is being used in this case, it is not required. As noted above, a system could instead compute the sum of the magnitudes rather than the magnitude of the sum and still be successful or leverage other processing.
(28) The signals s.sub.hh and s.sub.vv are shown in plots 26a, 26b respectively and are identified with reference numerals 27a, 27b. For the purposes of this example, P is fixed at P=100 for all CPIs, regardless of waveform. This is an arbitrary selection and the rationale behind this choice is that it is well known that the transient characteristics of sea clutter decorrelate on timescales less than 10 msec at X-band (Chan, 1987; Antipov, 1998). Consequently, choosing P=100 which comprises about 16 milliseconds at the highest PRF of the so-called Four Eyes radar ensures averaging over timescales at least as long as the usual decorrelation time of transient sea clutter characteristics.
(29) The effect of coherently averaging these P pulses is that one has averaged over the Doppler response of the scatterers in each range bin. That is, the average Doppler frequency of each channel in each bin has been captured.
(30) The operation typically performed by a coherent radar's standard radar signal processing chain on a CPI of pulses is Doppler processing. Doppler processing amounts to computing the discrete Fourier transform of samples across slow-time in each range bin to determine the strength of each of the complex exponentials measured across the CPI time at each of the velocities measured by the radar's slow-time sampling frequency.
(31) The location of a scatterer at range R.sub.0 from the radar as a function of slow-time in either channel (i.e. either hh or vv) can be written as:
(32)
in which:
(33) A captures signal amplitude and attenuation;
(34) c=the speed of light;
(35) p=the pulse number in the CPI pulse sequence;
(36) Rp=the scatterer's range to the radar on the pth pulse; and
(37) T [sec] is the pulse length.
(38) It should be noted that in the above, the stop-and-hop approximation, which assumes target motion stops while the radar pulse is in transit, has been used.
(39) The first exponential term captures a phase shift relative to all terms of the slow-time sample sequence. The second exponential term captures the Doppler frequency of the scatterer.
(40) The radar's slow-time sampling frequency corresponds to the PRF comprised by the second exponential term in Equation (2).
(41) By coherently averaging across slow-time instead of Doppler processing, the signal SNR is improved by a factor of P, but rather than doing so by integrating the strength of the individual Doppler velocities of all scatterers in each range bin, the average Doppler velocity in each range bin in each channel is captured.
(42) The coherent averages of signals s.sub.hh and s.sub.vv are thus understood to capture the strength of the average Doppler response in each range bin. Note again, however, that the average Doppler response in each range bin is captured in noncoherent averages by monitoring the displacement of signal returns relative to the radar across time; hence, noncoherent averaging still works as a means of capturing the average Doppler response.
(43) Passing signals s.sub.hh and s.sub.vv 27a, 27b through a 1-D CFAR detection process, yields two sets of detections 32a, 32b and 36a, 36b as illustrated in plot 30. Detections 32a, 32b and 36a, 36b are functions of range-only in each channel and may be produced, for example, by CA-CFAR processing applied to the signals 27a, 27b.
(44) The probability of false alarm for the CA-CFAR results was set to Pfa=10.sup.6. It should be noted that in this illustrative example the 4 MHz waveform is being used as an instructive case. Thus, conclusions should not be drawn regarding the number of detections that 1-D CFAR operations will yield in the cases of finer resolution waveforms.
(45) As illustrated in plot 30, two pairs of detections 32a, 32b, and 36a, 36b are co-located in range across HH and VV polarization curves 27a, 27b. One of these pairs comprises detections on sea clutter while the other comprises detections on the buoy.
(46) The manner in which those detections corresponding to sea clutter and those detections corresponding to the buoy, may be most easily visualized (and will be described below) with reference to
(47) Briefly, however, it is recognized that the locations of detections on the moving sea surface will change from CPI to CPI, or in the noncoherent case, from noncoherent CPI (NCPI) to NCPI. However, because HH and VV capture different scattering mechanisms of the sea surface with different mean Doppler frequencies and the mean Doppler of HH is more significant than the mean Doppler of VV, it is expected that even those sea clutter detections that persist from CPI to CPI will move across range at different rates with respect to one another. In other words, sea clutter detections that are co-located on one CPI will not remain co-located across a series of CPIs.
(48) If the radar is imagined as an observer floating in the sea at the stationary range equal to the waveform's minimum range, then such detections in HH will wash over the observer or move away from the observer, depending on look direction faster than will the corresponding detections in VV. Thus, one can set a minimum number of CPIs across which a pair of detections must retain their co-location in range in order to be declared detections on objects rather than detections on the sea.
(49) Recognizing also that the strength of an extended object's signature will vary somewhat in range as a function of polarization for higher bandwidth measurements, one can opt to relax the requirement that detections must be precisely co-located in range. This results in a definition of polarimetric co-location (PCL): a pair of detections on the coherently averaged signals s.sub.hh and s.sub.vv are considered polarimetrically co-located if they remain approximately co-located in range across a certain number of CPIs. Approximate co-location in range is captured by an input parameter .sub.n, which is the maximum permitted offset in range that will still allow a detection in HH to be considered co-located with an associated detection in VV. The number of CPIs across which a given detection pair must remain approximately co-located in range is captured by another algorithm input parameter, the PCL CPI criterion, denoted n.sub.CPI.
(50) Referring now to
(51) As the sequence shows, for the 4 MHz waveform, the resolution is low enough such that it suffices to allow no range offset between detection pairs. Moreover, so few detections are generated on the 4 MHz waveform data that three CPIs sufficiently allows sea clutter detection pairs to separate over time due to sea clutter's mean Doppler differential.
(52) On this image sequence, VV and HH signals are indicated with reference numerals 40a, 40b, respectively, and VV and HH detections are indicated with reference numerals 42, 44, respectively. Detection pairs that have remained co-located in range for as long as n.sub.CPI CPIs are indicated as PCL detections.
(53) From the sequence, it is seen that on CPI 3 (
(54) It should be note that a co-located detection pair on sea clutter appears in CPI 10 (
(55) Moreover, approximately co-located pairs appear in CPIs 2, 4, 6, 15, and 18 (
(56) For example, and with reference to
(57) PRF=6250 Hz;
(58) BW=40 MHz (R=3.75 m);
(59) .sub.n=2; and
(60) n.sub.CPI=3
(61) detections 52a-52p (sixteen (16) HH detections) and 54a-54j (ten (10) VV detections) occur in
(62) In
(63) PRF=6250 Hz;
(64) BW=4 MHz (R=37.5 m);
(65) n=0; and
(66) n.sub.CPI=3
(67) only detections 56a-56b (two (2) HH detections) and 58a (one (1) VV detection) occur in
(68) Referring now to
(69) Radar 60 includes an illustrative conventional radar signal processing chain comprising means for receiving radar data (e.g. digitized radar data) 62 and a pulse compression circuit 62 which receives and processes the radar data using any technique known to those of ordinary skill in the art. A processor produces range-time intensities (RTIs) 64 in both HH and VV using the pulse compressed data.
(70) A Doppler processor 66 then performs Doppler processing (again, using any technique known to those of ordinary skill in the art) to produce range-Doppler data (RD maps) 68 in both HH and VV polarizations. A detector 70 (e.g. a conventional CFAR detector) can then be used to identify detections 72 in both HH and VV.
(71) PCL circuits and processors can be integrated into the above described radar signal processing chain as a parallel process. As discussed above, only coherent integration on already pulse-compressed data followed by 1-D CFAR detection operations in HH and VV channels are required to execute the PCL clutter mitigation technique.
(72) To integrate PCL system 61 and the associated technique, the RTIs are provided to a coherent processing 74 which produces coherent average values 76 in both HH and VV polarizations. The CP data is then provided to a PCL processor 78 which processes the data provided there to form associated detection pairs as described above.
(73) A PCL detection filter 80 is coupled to receive both CFAR detection data as well as PCL data. PCL detection filter then filters the detections to eliminate detections moving at different rates of speed within each associated detection pair.
(74) It should be recognized that in the case of a noncoherent radar, the Doppler processing stage is simply omitted, and PCL is incorporated as a serial rather than a parallel process in the detection chain.
(75) The so-filtered detection data is then provided to a tracking and parameter estimator 82 of the radar as is generally known and then subsequently to optional discrimination logic circuit 84 and to an output system which may include, for example, a console/display 86 and/or a system for recording data 88.
(76) It should be noted that a clutter mitigation system and technique (e.g. a PCL sea clutter mitigation technique) is implemented as a parallel process that executes while Doppler processing 66 and CFAR detection (e.g. 2-D CFAR detection) 70 are performing computations. In this illustrative embodiment, it is only necessary to introduce a single processing element serially in the existing radar signal processing chain. In particular, in this example, only detection filter 80 (in this case corresponding to a PCL detection filter 80) which performs clutter mitigation filtering of CFAR detections is included serially in the radar signal processing chain.
(77) It should be appreciated that the output of PCL processor 78 in this structure and for each CPI is a set of range bins containing PCL detections. PCL detection filtering simply filters out 2-D CFAR detections (e.g. as produced by a standard radar signal processing chain) whose ranges are not included in PCL range bin set.
(78) In the case of a noncoherent radar, PCL detection filtering filters out the 1-D CFAR detections in the RTIs whose ranges are not included in the PCL range bin set.
(79) This may be understood by considering a 2-D CFAR that produces detections in all positions indicated with a 1 in the notional binary matrix:
(80)
where the row corresponds to the Doppler bin and the column corresponds to the range bin.
(81) A value of one (1) in the matrix is analogous to 2-D CFAR producing a black circle indicating a detection on the CFAR detection RD images 72 in
(82)
(83) It should, of course, be appreciated that in practical embodiments, there are likely to exist hundreds (or even thousands of false alarms, depending upon radar parameters as discussed above in conjunction with
(84) To integrate PCL for mitigation of sea clutter in parallel with conventional radar processing, it has been found that in an embodiment using linear HH and VV signals additional computation time can be avoided if HH and a VV detections are: (1) within .sub.n range bins of each other and (2) over n.sub.CPI contiguous CPIs.
(85) The sea clutter mitigation system and technique described herein may be applied to both stationary and moving targets. With respect to moving targets, allowance must be made for detection displacement across CPIs. A highest velocity of sea clutter is approximately the wind speed.
(86) Wind speed depends on sea state: v.sub.maxw.sub.max (sea state)
Thus, one can map to the number of range bins a PCL detection may move over the course of one (1) CPI:
(87)
(88) It is, of course, necessary to allow for polarimetric signature variation across HH and VV. For example, in practical systems it is expected that HH and VV will each miss at least one or more detections on extended targets in some CPIs. To overcome this, one can allow propagation of PCL detections
(89) .sub.n: maximum pairwise range offset
(90) n.sub.CPI: CPI criterion
(91) .sub.max: detection displacement limit
(92) n.sub.prop: CPI propagation allowance
(93) Having described preferred embodiments, which serve to illustrate various concepts, structures and techniques, which are the subject of this patent, it will now become apparent to those of ordinary skill in the art that other embodiments incorporating these concepts, structures and techniques may be used. Additionally, elements of different embodiments described herein may be combined to form other embodiments not specifically set forth above.
(94) Accordingly, it is submitted that that scope of the patent should not be limited to the described embodiments but rather should be limited only by the spirit and scope of the following claims.