IL300197B1 - Novel method and system for target height estimatiom using multipath corrupted radar measurements - Google Patents

Description

26/11/2023 Novel Method for Target Height Estimation Using Multipath Corrupted Radar Measurements FIELD OF THE INVENTION The invention is directed to a method for mitigating the effects of multipath on radar target tracking, and specifically, to the use of a neural network to improve target height estimation.

BACKGROUND In radar tracking of low-flying targets, the altitude of a target as measured by a monopulse radar is influenced by a phenomenon that is commonly known as "multipath". Multipath is especially severe when the main beam of the radar antenna illuminates both the target and the Earth's surface. When the range difference between the backscattered (direct ray) and one or more specularly reflected rays propagating between the target, the surface, and the radar is smaller than the radar range resolution, phase interference occurs which give rise to the large fluctuations in the apparent target elevation, known as the multi-path effect. These fluctuations depend primarily upon the instantaneous values of the radar wavelength ( ), the radar height above the surface, the true target height, and the slant range from radar to target. These fluctuations behave as a slowly changing bias which can be easily mistaken for target maneuvers.

In the prior art, multipath is generally treated as a source of angle measurement error which is to be reduced or eliminated. One approach to mitigate multipath effects is to use radar frequency diversity to average over, or "smooth", phase fluctuations at the radar. Other approaches include antenna beam null steering and the use of linear or nonlinear Kalman filters. All of these approaches have difficulty in distinguishing between slow variations in the measured target elevation which are caused by specular multipath reflections and those which are caused by actual target maneuvers in elevation.

US Patent Number 10,209,344 (hereinafter '344) to Dougherty et al., entitled "Methods and Systems for Mitigating Multipath Propagation", and dated February 19, 2019, discloses a system and method for mitigating multipath propagation which may include collecting a plurality of detections of a target, forming a plurality of models each assuming at least one parameter causing multipath propagation, determining which model best fits the detections of the target, using the 26/11/2023 best fit model to approximate the ground conditions, and using the approximated ground conditions to remove the multipath error from the observed signals.

In the field of machine learning, a Sequence-to-Sequence (Seq2Seq) neural network is one which converts sequences from one domain to sequences in another domain. A recurrent neural network (RNN) is a type of Seq2Seq neural network in which there are connections between nodes wherein the output from some nodes affects the subsequent input to the same nodes. This enables an RNN to model the dynamics of time-varying systems and to process input data sequences of variable length.

A paper by T. Feigl et al., entitled "Recurrent Neural Networks on Drifting Time-of-Flight Measurements", published in proceedings of the IEEE 2018 International Conference on Indoor Positioning and Indoor Navigation, teaches a data-driven approach that uses training sequences to derive a near-optimal position estimator. In the presence of multipath propagation, time-of-flight (ToF) measurement errors drift due to small-scale motion. This results in changing phases of the multipath components which cause a drift in the ToF measurements. A Long Short-Term Memory RNN learns to interpret drifting errors in ToF measurements of a tracked dynamic object directly from raw ToF data.

While the above-mentioned paper by Feigl et al. describes a target location system based on RF beacon reception by a set of synchronized antennae in an indoor controlled environment having rigid walls with stationary reflective characteristics, the present invention deals with a completely different system: a moving radar in a non-stationary environment with unknown reflection and target characteristics.

SUMMARY OF THE INVENTION This application presents a paradigm shift in the treatment of multipath for target height estimation. Rather than viewing it as a source of error or noise, the effect of multipath is leveraged based on the understanding that the large deviations around the true line-of-sight caused by multipath are deterministic in nature, with characteristic patterns when plotted vs. radar-to-target range, and contain information related to the true target altitude above the reflective surface. A 26/11/2023 dedicated neural network is constructed, which properly recognizes multipath patterns in a wide variety of radar-target scenarios, thus providing an estimate of the true target height.

For the purpose of clarity, it should be noted that, for a person skilled in the art of artificial intelligence and machine learning, any Sequence-to-Sequence (Seq2Seq) neural network could be used in the implementation.

Furthermore, the term "plot" used hereafter refers to a radar report including the current raw radar measurements of: target slant range, Doppler frequency shift, line-of-sight angles (in azimuth and elevation), and the operating RF frequency (or wavelength).

The method of the invention consists of an off-line training stage followed by a real-time stage. In the training stage, a neural network is trained and constructed to recognize multipath effects as they appear in a wide variety of off-line simulated radar-target scenarios. In the real-time stage, the neural network is used to generate an estimate of the true target height and an estimation confidence factor, based on a sliding window of multipath corrupted (MPC) raw radar data.

In contrast to other classical estimation methods, the present invention does not require any prior assumptions about the target dynamics or other stochastic parameters, such as target RCS and surface reflection coefficients.

Furthermore, in contrast with patent ‘344, in which multipath phenomena are handled as errors and therefore are subtracted from real-time MPC data according to a best-fit scenario [e.g. see column 6 lines 28-33 of patent ‘344 which recites "Once a multipath solution has been identified over a portion of a target's path, the calculated multipath signal … may be output to the signal correction module 18, where it is subtracted from the observed signal"], the present invention leaves the multipath phenomena in place, and leverages the valuable information and characteristics of the multipath signal by using a pre-trained neural network for providing the target true altitude.

According to one aspect of the presently disclosed subject matter, there is provided a method for estimating a true target height from received radar signals scattered from the target and a surface in a multipath scenario. The method includes the steps of: deriving from a sequence of 26/11/2023 radar measurements received over a period of time successive data pairs of a measure of slant range and a multipath corrupted (MPC) target elevation; generating from said data pairs a representation of MPC target elevation with slant range over a range window having a variable length; inputting said representation to a Sequence-to-Sequence (Seq2Seq) neural network, said neural network being trained on a dataset comprising variations of MPC target elevation with slant range for targets at different heights; and generating an estimate of the true target height using said neural network.

According to some aspects, the neural network is a recurrent neural network (RNN).

According to some aspects, the RNN is a gated recurrent unit (GRU) neural network.

According to some aspects, the surface is a sea or land surface.

According to some aspects, the successive data pairs correspond to different radar operating wavelengths.

According to some aspects, the radar signals are those of a radar employing frequency agility techniques.

According to some aspects, the radar signals are those of a radar employing phase or amplitude monopulse techniques.

According to some aspects, the radar signals are those of a radar mounted on a moving platform.

According to some aspects, the radar measurements incorporate Doppler frequency shifts which may be used as a decision aid to differentiate between approaching, receding, or tangential target trajectories.

According to some aspects, the dataset further includes variations of MPC target elevation with slant range for targets executing maneuvers.

According to some aspects, the dataset includes variations of the MPC target elevation with slant range obtained from computer simulations and/or experimental radar recordings. 26/11/2023 According to some aspects, the dataset further includes effects of earth curvature and/or atmospheric refraction.

According to some aspects, the neural network includes interpolation to a predetermined discrete range grid.

According to some aspects, the neural network comprises a sliding window algorithm.

According to some aspects, the sliding window algorithm includes a backwards sliding window with an adaptive history length.

According to some aspects, the method further includes generating a confidence factor corresponding to the estimate of the true height of the target.

According to another aspect of the presently disclosed subject matter, there is provided a computer-readable storage medium including computer-executable instructions that, responsive to execution by a processor, implement the method of claim 1.

According to yet another aspect of the presently disclosed subject matter, there is provided a radar system configured to implement the method of claim 1.

BRIEF DESCRIPTION OF THE DRAWINGS The invention is herein described, by way of example only, with reference to the accompanying drawings. Like reference numerals are used to denote similar or like elements in the drawings.

FIG. 1: A definition sketch of a multipath scattering geometry.

FIG. 2A: An exemplary graph of MPC target elevation measurements vs. slant range for a non-maneuvering target.

FIG. 2B: An exemplary graph of MPC target elevation measurements vs. slant range for a target executing a maneuver. 26/11/2023 FIG. 3: A block diagram of an exemplary algorithm for forming a training database of radar measurements corresponding to a multiplicity of multipath scenarios.

FIG. 4: A block diagram of an exemplary algorithm for training a neural network, in accordance with the invention.

FIG. 5: A block diagram of an exemplary algorithm for real-time estimation of a true target height by a neural network, in accordance with the invention.

DETAILED DESCRIPTION FIG. 1 shows a definition sketch of a multipath scattering geometry. The X-axis is a nominal horizontal axis, representing the reflecting surface of the Earth, which can be land or sea, in a simplified "flat-earth" model that is commonly used by those skilled in the art of radar engineering.

The definitions of the geometric symbols used in FIG. 1 are given in the following table.

Table 1: Symbols used FIG.

R Radar position T Target position R' Radar image position T' Target image position S Specular reflection point hR Radar height hT True target height Rd Slant range to target xT Ground range to target R1, R2 Specularly reflected ray path lengths d Elevation angle of direct ray  Grazing angle of specularly reflected ray 26/11/2023 In FIG. 1, the direct ray is represented by the segment R-T, and specularly reflected rays are denoted by the segments R-S and S-T. To observers positioned at the radar and at the target, the reflected rays appear to emanate from the target image T' and the radar image R', respectively. Note that propagation along paths R-T, R-S and S-T is bi-directional, so that there are four possible ray combinations: direct-direct, reflected-reflected, direct-reflected and reflected-direct. The dashed line 110 corresponds to the boresight of a sum channel of a monopulse radar antenna. For simplicity and clarity of presentation, the dashed line 110 is drawn as being approximately parallel to the reflecting surface represented by the x-axis, and effects due to curvature of the Earth and of electromagnetic rays refracted by the Earth's atmosphere are not shown in the simplified planar geometry of FIG. 1. When the Earth is treated as a curved surface, the term "target height" denotes the distance of the target from the surface along the local vertical direction. Methods to compensate for the effects of Earth curvature are generally well-known to those skilled in the art of radar engineering.

Radar antenna pattern 120 shows the dependence of a monopulse sum channel gain on the elevation angle relative to boresight. In monopulse radars, the radar antenna also includes a difference channel which enables estimation of a target elevation angle by comparing the amplitude and/or phase of signals received in the sum and difference channels, using techniques that are familiar to those skilled in monopulse radar technology. The radar also measures a slant range to the target, Rd, based on a time-of-arrival of the received radar signal, using standard range-gating techniques. When the radar beam illuminates both the target and the reflective surface, the direct and reflected signals combine at the radar. If the range and Doppler frequency shift difference between the direct and reflected signals are simultaneously smaller than the radar resolution cell, interference gives rise to the MPC angular measurements errors.

Target RCS pattern 130 shows the variation of the RCS with target aspect angle. For targets consisting of multiple scattering elements, the RCS pattern 130 may be highly anisotropic. The RCS corresponding to the direct ray in the direction T-to-R is referred to as a "monostatic" RCS; and the RCS corresponding to the reflected ray in the direction T-to-S is referred to as a "bistatic" RCS.

Region 140 represents a Fresnel zone on the reflecting surface which surrounds the specular reflection point S. The complex Fresnel reflection coefficient depends, in amplitude and 26/11/2023 phase, on the frequency-dependent electromagnetic properties of the surface, including, for example, its electrical conductivity and dielectric constant. Surface roughness may give rise to diffuse reflections, which, unlike specular reflections, are amenable to mitigation by conventional averaging techniques.

FIG. 2A shows an exemplary graph of multipath corrupted (MPC) target height measurements vs. target slant range, Rd, for a non-maneuvering target, at fixed values of , hR and hT, and under variable sea state conditions. The dashed line 240 marks the ground truth which, in this example, corresponds to a target flying at a constant altitude of 150 meters above the reflecting surface. The characteristic peaks, 220a and 220b, and valleys, 230a and 230b, are caused by a phase difference, 2| R|/ , where R=2(R1+R2–Rd) is the two-way path length difference, between the direct and reflected rays (for the reflected-reflected ray combination).

The peaks 220a and 220b, and the valleys 230a and 230b, lie on differently colored curves in FIG. 2A. The different colors correspond to different sea states. The peaks and valleys change in amplitude but their positions on the radar-to-target slant range axis are unaltered. A similar deterministic relationship holds true for other stochastic parameters of the scenario, such as the monostatic and bistatic target RCS and the surface reflection coefficients.

FIG. 2A indicates that, for fixed values of , hR, and hT, there exists a deterministic relationship between the peaks and valleys in the MPC target elevation measurement and the slant range Rd. This is a direct consequence of the physical nature of specular reflections.

FIG. 2B shows an exemplary graph of MPC target height measurements vs. target slant range for a target executing a maneuver. The sea state is constant, and the measurements are shown by the curve 250. The true target trajectory consists of level flight at a constant height of 250 meters, shown by dashed line 240a, followed by a maneuver, shown inside the dashed oval 240b, and then followed by level flight at a constant height of 850 meters, shown by line 240c. During the maneuver, the true target height increases linearly from 250 to 850 Inside oval 240b, the MPC measurements undergo rapid oscillations between peaks and valleys. The oscillations are a consequence of a rapid phase wrapping that occurs as the vertical separation increases between the target T and its image T'. 26/11/2023 FIG. 3 shows an exemplary algorithm 300 for forming a training database of radar measurements corresponding to a multiplicity of multipath scenarios. Block 310 incorporates parameters describing a target's trajectory, such as its time dependent position, velocity and acceleration vectors, for a multiplicity of possible target maneuvers. Block 320 incorporates radar parameters, such as the wavelength, polarization, and power of the radar transmitter, the gain patterns of the sum and delta channels of the radar antenna, and the height of the radar platform above the reflecting surface of the Earth. Block 330 incorporates parameters describing the angular dependence of the monostatic and bistatic target RCS. The latter depend also on the wavelength and polarization of the radar's transmission. Block 340 incorporates parameters describing the topography of the reflecting surface, such as sea-state or terrain slope, as well as electromagnetic properties of the reflecting surface, such as Fresnel reflection coefficients which depend upon wavelength, polarization and grazing angle . Blocks 350 and 360 represent software modules which simulate the radar measurements for a given multipath scenario. The output of block 360 is validated by comparison with data collected in radar field tests, which are stored in block 370. Data-flow paths 365a and 365b represent the validation feedback loop connecting blocks 360 and 370. Range discretization block 380 receives as input a time-series of radar measurements from block 360. The radar measurements include the radar-target slant range, Rd, the MPC target elevation and, optionally, the target Doppler frequency shift which is proportional to a time derivative of the slant range. In accordance with the dependence of the multipath bias pattern on range-height geometry, block 380 converts and resamples the time sequence of radar measurements into discrete values corresponding to a predetermined slant range grid. The latter conversion typically requires the use of digital interpolation and re-sampling techniques, which are familiar to those skilled in the art of digital signal processing.

Range discretization block 380 outputs sequences of data triplets, which are stored in a comprehensive training database 390. Database 390 is configured to provide output to the supervised optimization process, where the output 395 includes successive values of the slant range after discretization and resampling, the MPC target height and the associated true target height (hT). 26/11/2023 FIG. 4 shows a block diagram of an exemplary algorithm 400 for training a neural network, in accordance with the invention. RNN 420 receives as input from database 390 a sequence {hT,i , ui , i= 0, 1, … n}, where "i" is an index of the sequence, hT,i is the instantaneous true target height, and ui is a data pair whose two components are a discretized value of slant range and an MPC target height. In RNN 420, the parameters win, wout, and wst denote weighting parameters whose values are to be determined by optimization. Each hidden state hi is calculated as a function of the input data pair ui to generate an estimate of true target height, h_esti. The estimate error, ei = h_esti-hT , is input to optimization block 430. Optimization block 430 calculates a loss function, defined as: n Loss =  ei (equation 1) i= and, using for example a first-order iterative minimization algorithm such as stochastic gradient descent, determines the optimal weights, k (w*in, w*out, w*st) = argmin (1/k)  Lossj (equation 2) (win, wout, wst) j= where j is an index chosen randomly from the dataset, k is the batch size for the stochastic gradient descent process, and "argmin" is a minimization over the weighting parameters win, wout, and wst. RNN 420 is configured to operate on sequences having variable length. An alternative embodiment of the neural network may comprise a sliding window algorithm, and the sliding window algorithm may include a backwards sliding window with an adaptive history. Furthermore, the neural network may be implemented as a gated recurrent unit (GRU) neural network.

FIG. 5 shows a block diagram of an exemplary algorithm 500 for real-time estimation of a target height by a neural network, in accordance with the invention. Radar tracking system 5provides a time sequence of radar measurements 515a and auxiliary parameters 515b. The radar measurements 515a include target slant range and MPC target elevation angle. The auxiliary parameters 515b include instantaneous radar wavelength, , closing velocity (Doppler) and radar height hR. Preprocessing block 520 performs a weighted interpolation of the radar measurements 26/11/2023 to the nearest slant range grid point and outputs a sequence of discretized data pairs to a trained RNN 530. The trained RNN provides a real-time estimate of target height, hT_est, and an associated confidence factor, CF. The computing resources required by the trained RNN block 530 are highly efficient for real-time computation. The values of hT_est and CF are typically fed back to a tracking algorithm which is resident in the radar tracking system 510, as indicated by paths 535a and 535b. It will be appreciated that the above descriptions are intended only to serve as examples, and that many other embodiments are possible within the scope of the present invention as defined in the appended claims.

Claims (18)

26/11/2023 CLAIMS

1. A method for estimating a true height of a target (hT) from received radar signals scattered from the target and a surface in a multipath scenario, the method comprising the steps of: a) deriving from a sequence of radar measurements received over a period of time successive data pairs of a measure of slant range and a multipath corrupted (MPC) target elevation; b) generating from said successive data pairs a representation of MPC target elevation with slant range over a range window having a variable length; c) inputting said representation to a Sequence-to-Sequence (Seq2Seq) neural network, said neural network being trained on a dataset comprising variations of MPC target elevation with slant range for targets at different heights; and d) generating an estimate of the true height of a target using said neural network.

2. The method of claim 1 wherein said neural network is a recurrent neural network (RNN).

3. The method of claim 2 wherein the RNN is a gated recurrent unit (GRU) neural network.

4. The method of claim 1 wherein the surface is a sea or land surface.

5. The method of claim 1 wherein the successive data pairs correspond to different radar operating wavelengths.

6. According to some aspects, the radar signals are those of a radar employing frequency agility techniques.

7. The method of claim 1 wherein the radar signals are those of a radar employing phase or amplitude monopulse techniques.

8. The method of claim 1 wherein the radar signals are those of a radar mounted on a moving platform.

9. The method of claim 1 wherein the radar measurements incorporate Doppler frequency shifts.

10. The method of claim 1 wherein the dataset further comprises variations of the MPC target elevation with slant range for targets executing maneuvers.

11. The method of claim 1 wherein the dataset comprises variations of the MPC target elevation with slant range obtained from computer simulations and/or experimental radar recordings. 26/11/2023

12. The method of claim 1 wherein the dataset further comprises effects of earth curvature and/or atmospheric refraction.

13. The method of claim 1 wherein the neural network includes interpolation to a predetermined discrete range grid.

14. The method of claim 1 wherein the neural network comprises a sliding window algorithm.

15. The method of claim 14 wherein the sliding window algorithm comprises a backwards sliding window with an adaptive history length.

16. The method of claim 1 further comprising generating a confidence factor associated with the estimate of the true height of the target.

17. A computer-readable storage medium comprising computer-executable instructions that, responsive to execution by a processor, implement the method of claim 1.

18. The method of claim 1 wherein said sequence of radar measurements received over a period of time is provided by a radar system.