WO2018210887A1 - Procede d'evaluation de caracteristiques d'une cible par filtrage d'echos radar - Google Patents

Procede d'evaluation de caracteristiques d'une cible par filtrage d'echos radar Download PDF

Info

Publication number
WO2018210887A1
WO2018210887A1 PCT/EP2018/062628 EP2018062628W WO2018210887A1 WO 2018210887 A1 WO2018210887 A1 WO 2018210887A1 EP 2018062628 W EP2018062628 W EP 2018062628W WO 2018210887 A1 WO2018210887 A1 WO 2018210887A1
Authority
WO
WIPO (PCT)
Prior art keywords
time
target
signal
coefficients
radar
Prior art date
Application number
PCT/EP2018/062628
Other languages
English (en)
Inventor
Jacques PARENT DU CHÂTELET
Mohammed TAHANOUT
Original Assignee
Meteo-France
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Meteo-France filed Critical Meteo-France
Publication of WO2018210887A1 publication Critical patent/WO2018210887A1/fr

Links

Classifications

    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
    • G01S13/00—Systems using the reflection or reradiation of radio waves, e.g. radar systems; Analogous systems using reflection or reradiation of waves whose nature or wavelength is irrelevant or unspecified
    • G01S13/88—Radar or analogous systems specially adapted for specific applications
    • G01S13/95—Radar or analogous systems specially adapted for specific applications for meteorological use
    • G01S13/951—Radar or analogous systems specially adapted for specific applications for meteorological use ground based
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
    • G01S13/00—Systems using the reflection or reradiation of radio waves, e.g. radar systems; Analogous systems using reflection or reradiation of waves whose nature or wavelength is irrelevant or unspecified
    • G01S13/02—Systems using reflection of radio waves, e.g. primary radar systems; Analogous systems
    • G01S13/06—Systems determining position data of a target
    • G01S13/08—Systems for measuring distance only
    • G01S13/10—Systems for measuring distance only using transmission of interrupted, pulse modulated waves
    • G01S13/22—Systems for measuring distance only using transmission of interrupted, pulse modulated waves using irregular pulse repetition frequency
    • G01S13/225—Systems for measuring distance only using transmission of interrupted, pulse modulated waves using irregular pulse repetition frequency with cyclic repetition of a non-uniform pulse sequence, e.g. staggered PRF
    • Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
    • Y02A90/00—Technologies having an indirect contribution to adaptation to climate change
    • Y02A90/10—Information and communication technologies [ICT] supporting adaptation to climate change, e.g. for weather forecasting or climate simulation

Definitions

  • the invention relates to a method for evaluating the characteristics of a target, in particular for evaluating the characteristics of an atmospheric target, in particular the quantity of rain and the speed of the raindrops, by means of a precipitation radar. .
  • Precipitation radars are measuring instruments deployed in operation for more than twenty years. They can detect precipitation, as well as the horizontal speed of the wind in precipitation, within a radius of a few hundred kilometers around the radar.
  • the data collected by the radars are transmitted in real time to concentration centers to produce precipitation and wind maps at the scale of a country, as with the French ARAMIS network, which includes 34 radars, or at the same time.
  • a continent such as with the US NEXRAD network, or the European OPERA network that gathers in real time the data collected by the radars of the contributing countries.
  • continent-wide precipitation maps are available every 15 minutes or even every 5 minutes for various applications: meteorology, hydrology, tourism, agriculture, air transport, risk prevention, etc.
  • Rainfall radars operate in panoramic mode with low elevation angles between 0 ° and 5 °. They are thus faced with the difficulty of obtaining representative signals, at the lowest altitude, which are not disturbed by the ground returns of the lower part of the radar beam. These ground returns, also called “ground echoes,” can be very powerful, leading to a false assessment of the characteristics of the target, such as the amount of precipitation or the speed of displacement of the precipitation.
  • the common solutions to this problem are to filter the received signal by eliminating the low frequency components, exploiting the fact that the fixed contribution of the ground is zero Doppler, while that of the precipitation presents a Doppler shift due to the speed. the wind.
  • the ground returns are narrow band, while the precipitation signals are wider band due to the natural variability of the target.
  • the measurement of v which is based on the argument of a complex number, becomes ambiguous.
  • An object of the invention is to respond, at least partially, to this need.
  • the invention relates to a method for evaluating a characteristic of a target, preferably chosen from the speed, the density and the homogeneity of the target, the method comprising the following steps: a) transmission towards said target, by radar, during a residence time, an emitted signal consisting of a series of transmitted pulses, any two consecutive pulses being separated by a time base interval selected from a group consisting of at least one constant base interval; b) receiving a received signal returned by said target in response to the transmitted signal, the received signal consisting of a succession of pulses received, at times ti, in response to the transmitted pulses; c) determining a sampled time series S r (ti) of the received signal, S r (ti) being a complex value of the received signal, at time t; ; d) determining, from the time series S r (ti), several autocorrelation coefficients R (Tj), Tj being such that there are at least two pulses of the received signal separated by a duration equal to T
  • R (Tj) ⁇ S r (ti). S ti - Tj)> i, where * denotes the conjugate complex, and ⁇ > i the arithmetic mean over all the products (S r (ti) .s r * (ti - Tj)) computable from said series temporal; e) determining a lower cut-off time T c positive, greater than an estimate of the coherence time of the target T C ibi e ; f) determining an upper breaking time T c 'greater than T c and less than or equal to the residence time; g) analysis of the autocorrelation coefficients R (Tj) for absolute values of Tj greater than the lower cut-off time T c and less than the upper cut-off time T c ', referred to as "upper coefficients", so as to determine a value
  • a method according to the invention is thus particularly well suited when the radar transmits towards a target close to the ground.
  • a method according to the invention also has one or more of the following optional features: the lower cut-off time T c is chosen between 50 ms and 500 ms, preferably depending on the nature of the target; the lower cut-off time T c is measured from a signal neighboring the received signal obtained by orienting the radar towards a part of the chosen target so that the neighboring signal does not contain a ground contribution; in step f), the upper cut-off time T c 'is chosen equal to the residence time of the radar on the target; in step g), R CO h is a weighted average of said upper coefficients; the weighting coefficients decrease with the approach of the cut-off times T c and T c 'and / or increase as the average of these two cut-off times approaches; in step h), subtract R CO h with said lower coefficients; in step a), the time interval between any two consecutive pulses is always equal to the same basic interval; alternatively, in step a), the time interval between any two consecutive pulses is selected from
  • the invention also relates to: - a computer program comprising program code instructions for the implementation of one or more, preferably all the steps c) to j) of a method according to the invention, when said program is executed by a computer;
  • an electronic circuit integrating functions performing one or more, preferably all the steps c) to j), for example an integration in a circuit of type, in English FPGA (Field Programmable Gate Array) or DPS (Digital Signal Processing);
  • a computer medium on which is recorded such a program for example a memory or a CD-ROM.
  • the subject of the invention is also a center comprising a computer containing a computer program according to the invention, and a precipitation radar, preferably several precipitation radars, connected to the computer, so as to transmit information to it. relating to signals transmitted and / or received by the said radar (s).
  • the center may be stationary or mobile, for example a land vehicle, an aircraft or a spacecraft.
  • dwell time refers to the time of illumination of the target by the radar, for example, when the antenna of the radar is rotating, the time during which the Radar beam passes on the target.
  • the pulses of the transmitted radar signal follow one another in a regular rhythm when the radar signal is constituted by the continuous repetition of a pattern consisting of several pulses.
  • the target target coherence time T is conventionally the time interval after which the autocorrelation coefficient of a signal received from the target is divided by two, relative to its value at zero time.
  • FIG. 1 schematically illustrates the operation of a radar according to the invention
  • FIG. 2 illustrates examples of pulses of the signal emitted by the radar
  • Figure 3 illustrates the power spectrum of a received signal
  • FIG. 4 represents the autocorrelation function of a received signal
  • FIG. 5 is a block diagram summarizing the various steps of the method according to the invention.
  • FIG. 1 schematically illustrates step a) of the method according to the invention where a precipitation radar 10 emits towards a target 12, in this case an atmospheric precipitation, for example a cloud of rain, snow or hail N, a signal S e consisting of a series of transmitted pulses 14e. Any two consecutive sent pulses are separated by a time base interval selected from a group consisting of at least one constant base interval, the pulses succeeding each other preferably in a regular rhythm.
  • a precipitation radar 10 emits towards a target 12, in this case an atmospheric precipitation, for example a cloud of rain, snow or hail N, a signal S e consisting of a series of transmitted pulses 14e.
  • a time base interval selected from a group consisting of at least one constant base interval, the pulses succeeding each other preferably in a regular rhythm.
  • Figure 2a shows the simple case PRT ("Pulse Repetition Time” in English), the time interval between any two consecutive sent pulses is always equal to the same base interval AU.
  • the time interval between any two consecutive pulses is chosen from a group consisting of two basic intervals AU and ⁇ 2 .
  • the time interval between any two consecutive pulses is chosen from a group consisting of three basic intervals AU,
  • FIGS. 2b and 2c illustrate pulses succeeding each other in a regular rhythm, the series of base intervals (AU, ⁇ 2 ) and (AU, ⁇ 2 , ⁇ 3 ) respectively repeating continuously.
  • each base interval is preferably between 100 and 5 ms.
  • the emitted signal S e preferably has a constant wavelength, preferably less than
  • the emitted signal S e is preferably sent in a transmission direction ⁇ forming an angle ⁇ , less than 15 °, at 10 °, preferably less than 5 ° with a horizontal plane 16 on which the radar 10 rests.
  • the radar 10 preferably comprises a parabolic antenna 18 and a transmitter 20 disposed at the focus of the antenna.
  • the radar 10 may comprise a flat antenna, of centimeter dimensions.
  • the antenna 18 focuses the electromagnetic energy in a beam 22 whose angular width a is inversely proportional to the diameter of the antenna.
  • the angular width a of the beam 22 is preferably greater than 0.5 ° and less than 5 °, preferably less than 3 °, at 2 °, more preferably less than 1 °.
  • the horizontal and / or vertical dimension of the field of view is preferably between 800m to 50km, better still between 4000m and 250km.
  • the antenna 18 rotates at an azimuthal rotation speed of preferably between 0.2 rpm and 5 rpm.
  • Step b) of the method according to the invention is also illustrated by FIGS. 1 and 3.
  • the radar receives, at a moment to, a received signal S r returned by the drops of water of the cloud N in response to the transmitted signal S e .
  • FIG. 3 represents the power amplitude A of the signal received in logarithmic scale as a function of the Doppler velocity v of this signal, normalized with respect to the maximum observable speed v a .
  • the narrow line S at zero speed corresponds substantially to the contribution of ground echo to the received signal.
  • the wide band P centered around the speed of 0.3 represents substantially the contribution of the echo of the target, in this case the echo of the cloud N.
  • the contribution of the additive noise B is in the form of a continuum.
  • the received signal S r is the sum of several signals whose coherence times are significantly different. More precisely, one classically distinguishes a coherent component, or "coherent part”, substantially constituted by the echo of the immobile obstacles, and in particular the ground echo, and an incoherent component or “incoherent part” substantially constituted of the echo from the atmospheric target.
  • a computer 24 in communication with the radar 10.
  • This computer 24 receives the received signal S r and executes the instructions of a program consisting in implementing the processing steps c) to j) of the process according to the invention.
  • step c) a sampled time series S r (ti) of the received signal is determined, S r (ti) being the complex value of the received signal S r at times ti, which are the moments of reception of the received pulses 14r in response to the pulses emitted 14th.
  • step d) from the time series S r (ti), several autocorrelation coefficients R (Tj) are determined, Tj being a duration such that there are at least two received pulses 14r separated by a duration equal to Tj:
  • R (Tj) ⁇ S r (t i ). S (t i - Tj)> i , where * denotes the conjugate complex, and ⁇ > i the arithmetic mean over all the products (S r (ti) S r * (ti - Tj)) computable from said time series S r (ti).
  • FIG. 4 represents the autocorrelation function R of the received signal S r as a function of time T or "lag". It makes it possible to read the values of R (Tj), Tj being a particular value of T.
  • the first two coefficients R (Ati) and R (2 ⁇ ) have been represented.
  • the autocorrelation function is the inverse Fourier transform of the power spectrum shown in FIG.
  • the contribution of the ground echo results in a continuum over the full extent of the autocorrelation function.
  • the narrow line PB at zero time also contains the contribution of noise, whose coherence time is practically nil.
  • the coherence time of the target Tdbie is estimated.
  • the estimate can be made by a usual value for the target in question.
  • it can be estimated that the coherence time of an atmospheric target is typically between 50ms and 500ms.
  • Another method consists of evaluating Tdbie from a neighboring signal S r 'coming from a region close to the target, but for which it is certain that it contains no coherent component.
  • the neighboring signal is a signal returned by the target in response to an emitted signal identical to S e , but in a transmission direction ⁇ 'forming with the horizontal plane an angle ⁇ ' greater than ⁇ , that is, ie sent to a part of the target at a higher altitude.
  • the contribution of the ground, and therefore the coherent component of the neighboring signal is negligible in such a neighboring signal.
  • the analysis of the neighboring signal S r allows, in a conventional manner to evaluate a coherence time Tdbie for the target.
  • a "lower" cutoff time T c greater than Tdbie is determined, preferably more than twice as much as Tdbie.
  • the autocorrelation coefficient R (Tj) is substantially identical to the autocorrelation coefficient that would have been obtained by analyzing only the coherent component of the received signal S r .
  • the contribution of the incoherent component of the received signal S r can be considered as substantially zero.
  • the cut-off time is typically between 50ms and 500ms. In some applications, however, it may be different.
  • step f an "upper" cut-off time T c ', greater than T c and less than the residence time of the radar on the target is determined.
  • T c ' can be any, provided that it remains lower than said residence time. More T c is large, the better the accuracy of the estimate higher autocorrelation coefficients.
  • step g) the analysis of autocorrelation coefficients R (Tj) for absolute values of Tj greater than the lower cut-off time T c , and lower than the upper cut-off time T c ', called "higher coefficients" , allows to determine a signal that substantially corresponds to the coherent component of the received signal S r .
  • the autocorrelation coefficients R (Tj) for absolute values of Tj greater than the lower cut-off time T c , and less than the upper cut-off time T c ', referred to as "upper coefficients", are therefore autocorrelation coefficients that would have been obtained if the received signal S r had been perfectly coherent, for example if the target had been fixed. From these coefficients, it is possible to determine a representative value Rcoh of the coherent component of the received signal.
  • the weighting coefficients k i decrease as the time of closure T c and T c 'and / or increase near the mean of these two cut-off times, so as to follow an evolution of the type " bell function "according to Tj.
  • step h the autocorrelation coefficients R (Tj) are processed for absolute values of Tj less than the coherence time of the target Tdbie estimated in step e), referred to as "lower coefficients", in R CO h function, so as to obtain corrected autocorrelation coefficients R (Tj) corr.
  • the characteristics of the target are evaluated from the corrected autocorrelation coefficients R (Tj) C orr, preferably according to the well-known "peer-tap" method.
  • the radial velocity v r of the target is determined by the "peer-drawn" method from the phase of the first coefficient of autocorrelation RT i) correlation of the received signal:
  • v r - ⁇ / (4 ⁇ Ti) arg (R (Ti) corr), where ⁇ is the wavelength of the transmitted radar signal.
  • the density and the homogeneity can be evaluated by the classical methods by exploiting the autocorrelation coefficients R Tj) CO rr.
  • a known method for calculating one or more characteristics of the target is applied.
  • the invention provides a simple and accurate method for evaluating one or more characteristics of a target. Simulation checks have been performed which show that soil echo rejection capabilities by such a method are greater than 40 dB.
  • the invention is not limited to the examples that have just been given, and can be used with advantage in the exploitation of radars of other types: radar wind profilers, ocean observation radars , satellite radar observation, etc.
  • a method according to the invention is useful for filtering the effect of ships when the target consists of waves.

Abstract

La présente invention concerne un procédé d'évaluation d'une caractéristique d'une cible comprenant : a) la détermination d'une série temporelle échantillonnée Sr(ti) d'un signal reçu; b) la détermination de coefficients d'autocorrélation R(Tj) à partir de Sr(ti), Tj étant la durée séparant au moins deux impulsions du signal reçu; c) la détermination d'un temps de coupure Tc, supérieur à une estimation du temps de cohérence de la cible; d) la détermination d'un temps de coupure Tc' supérieur à Tc; e) l'analyse des coefficients d'autocorrélation pour des valeurs absolues de Tj supérieures à Tc et inférieures à Tc' de manière à déterminer une valeur Rcoh; f) le traitement des coefficients d'autocorrélation pour des valeurs absolues de Tj inférieures à Tc en fonction de Rcoh, de manière à obtenir des coefficients d'autocorrélation corrigés R(Tj)corr; g) l'évaluation de ladite caractéristique en fonction des coefficients d'autocorrélation corrigés.

Description

PROCEDE D'EVALUATION DE CARACTERISTIQUES D'UNE CIBLE PAR
FILTRAGE D'ECHOS RADAR
Domaine technique
L'invention se rapporte à un procédé d'évaluation de caractéristiques d'une cible, en particulier pour évaluer les caractéristiques d'une cible atmosphérique, notamment la quantité de pluie et la vitesse des gouttes de pluie, au moyen d'un radar précipitation.
Etat de la technique
Les radars précipitations sont des instruments de mesure déployés en opérationnel depuis plus d'une vingtaine d'années. Ils permettent de détecter les précipitations, ainsi que la vitesse horizontale du vent dans les précipitations, dans un rayon de quelques centaines de kilomètres autour du radar. Les données collectées par les radars sont transmises en temps réel à des centres de concentration pour produire des cartes des précipitations et du vent à l'échelle d'un pays comme avec le réseau français ARAMIS qui comporte 34 radars, ou à l'échelle d'un continent, comme avec le réseau US NEXRAD, ou le réseau européen OPERA qui rassemble en temps réel les données collectées par les radars des pays contributeurs. C'est ainsi que des cartes de précipitations à l'échelle d'un continent sont disponibles, toutes les 15mn ou même toutes les 5mn, pour diverses applications : météorologie, hydrologie, tourisme, agriculture, transport aérien, prévention des risques, etc.
Les radars précipitations fonctionnent en mode panoramique avec des angles d'élévation faibles, entre 0° et 5°. Ils sont ainsi confrontés à la difficulté d'obtenir des signaux représentatifs, à l'altitude la plus basse, qui ne soient pas perturbés par les retours de sol de la partie inférieure du faisceau radar. Ces retours de sol, appelés aussi « échos de sol », peuvent être très puissants, amenant à une fausse évaluation des caractéristiques de la cible, comme la quantité de précipitation ou la vitesse de déplacement de la précipitation. Les solutions courantes à ce problème sont de filtrer le signal reçu en éliminant les composantes de basse fréquence, en exploitant le fait que la contribution du sol, fixe, est de Doppler nul, alors que celle des précipitations présente un décalage Doppler dû à la vitesse du vent. De plus, les retours de sol sont de bande étroite, alors que les signaux de précipitations sont de plus large bande en raison de la variabilité naturelle de la cible.
Cependant, lorsque la vitesse radiale de la précipitation est faible, ces techniques de filtrage détruisent une partie du signal utile. Des méthodes d'interpolation dans le domaine fréquentiel sont parfois mises en œuvre pour corriger cette difficulté. Ces méthodes sont toutefois complexes et peu précises.
Par ailleurs, pour une longueur d'onde radar λ, la vitesse maximale observable va est liée à la distance maximale observable ra (portée en distance du radar) par la relation va ra = cA/4, où c est la vitesse de la lumière. Lorsque la vitesse v est supérieure à cette valeur limite va, la mesure de v, qui est basée sur l'argument d'un nombre complexe, devient ambiguë. Compte tenu des longueurs d'onde radar utilisées, des vitesses rencontrées dans l'atmosphère et des portées requises, ces problèmes d'ambiguïté sont considérables et des techniques, dites de « PRT multiples - PRT pour Puise Répétition Time », ont été imaginées et mises en place aux USA (double PRT) et dans les pays européens (triple PRT en France par exemple). Mais les méthodes de filtrage classiques, évoquées plus haut pour éliminer les retours de sol, ne sont pas compatibles avec ces techniques.
Il existe donc un besoin pour un procédé d'évaluation de caractéristiques d'une cible, en particulier d'une précipitation atmosphérique, ne présentant pas ces inconvénients.
Un but de l'invention est de répondre, au moins partiellement, à ce besoin. Résumé de l'invention
L'invention concerne un procédé d'évaluation d'une caractéristique d'une cible, choisie de préférence parmi la vitesse, la densité et l'homogénéité de la cible, le procédé comprenant les étapes suivantes : a) émission, vers ladite cible, au moyen d'un radar, pendant un temps de résidence, d'un signal émis constitué d'une série d'impulsions émises, deux impulsions consécutives quelconques étant séparées par un intervalle de base temporel choisi dans un groupe constitué d'au moins un intervalle de base constant ; b) réception d'un signal reçu renvoyé par ladite cible en réponse au signal émis, le signal reçu étant constitué d'une succession d'impulsions reçues, à des instants ti, en réponse aux impulsions émises; c) détermination d'une série temporelle échantillonnée Sr(ti) du signal reçu, Sr(ti) étant une valeur complexe du signal reçu, à l'instant t; ; d) détermination, à partir de la série temporelle Sr(ti), de plusieurs coefficients d'autocorrélation R(Tj), Tj étant tel qu'il existe au moins deux impulsions du signal reçu séparées par une durée égale à Tj :
R(Tj) = < Sr(ti ) . S ti - Tj) >i, où * désigne le complexe conjugué, et < >i la moyenne arithmétique sur l'ensemble des produits (Sr(ti ) . Sr *(ti - Tj)) calculables à partir de ladite série temporelle ; e) détermination d'un temps de coupure inférieur Tc positif, supérieur à une estimation du temps de cohérence de la cible TCibie; f) détermination d'un temps de coupure supérieur Tc' supérieur à Tc et inférieur ou égal au temps de résidence ; g) analyse des coefficients d'autocorrélation R(Tj) pour des valeurs absolues de Tj supérieures au temps de coupure inférieur Tc et inférieures au temps de coupure supérieur Tc', dits « coefficients supérieurs », de manière à déterminer une valeur
Figure imgf000005_0001
h) traitement des coefficients d'autocorrélation R(Tj) pour des valeurs absolues de Tj inférieures au temps de coupure inférieur Tc, dits « coefficients inférieurs », en fonction de RCOh , de manière à obtenir des coefficients d'autocorrélation corrigés
Figure imgf000005_0002
i) évaluation de ladite caractéristique de ladite cible en fonction des coefficients d'autocorrélation corrigés, de préférence suivant la méthode « puise pair ». Les inventeurs ont constaté qu'un tel procédé permet d'évaluer une caractéristique d'une cible d'une façon plus simple et plus précise que les méthodes habituelles de l'état de l'art.
En particulier, sans être liés par cette théorie, les inventeurs expliquent ce résultat du fait que lorsque le signal reçu Sr comporte une composante cohérente, résultant en particulier de la réflexion du signal émis sur des obstacles fixes, et notamment sur le sol, RCOh est représentative de cette composante cohérente. Le procédé permet de filtrer l'effet de cette composante cohérente.
Un procédé selon l'invention est ainsi particulièrement bien adapté lorsque le radar émet vers une cible proche du sol.
De préférence, un procédé selon l'invention présente encore une ou plusieurs des caractéristiques optionnelles suivantes : le temps de coupure inférieur Tc est choisi entre 50ms et 500ms, de préférence en fonction de la nature de la cible ; le temps de coupure inférieur Tc est mesuré à partir d'un signal voisin du signal reçu obtenu en orientant le radar vers une partie de la cible choisie de façon que le signal voisin ne contienne pas de contribution de sol ; à l'étape f), le temps de coupure supérieur Tc' est choisi égal au temps de résidence du radar sur la cible ; à l'étape g), RCOh est une moyenne pondérée desdits coefficients supérieurs ; les coefficients de pondération diminuent à l'approche des temps de coupure Tc et Tc' et/ou augmentent à l'approche de la moyenne de ces deux temps de coupure ; à l'étape h), on soustrait RCOh auxdits coefficients inférieurs ; à l'étape a), l'intervalle de temps entre deux impulsions consécutives quelconques est toujours égal au même intervalle de base ; alternativement, à l'étape a), l'intervalle de temps entre deux impulsions consécutives quelconques est choisi dans un groupe constitué par plusieurs, en particulier par deux ou trois intervalles de base ; chaque intervalle de base est compris entre 100 μ& et 5 ms ; - la cible est constituée d'un ensemble de diffuseurs distribués aléatoirement dans l'espace et en mouvement à des vitesses variables selon les diffuseurs considérés ; la cible est une couche de l'atmosphère ; à l'étape a), les impulsions sont envoyées suivant une direction formant un angle inférieur à 5° avec un plan horizontal sur lequel repose ledit radar ; - le signal émis a une longueur d'onde constante inférieure à 10 cm, notamment inférieure à 5 cm ; le radar comporte une antenne parabolique de diamètre compris entre 2 m et 8 m ; l'antenne du radar est fixe ou tournante ; l'antenne tourne à une vitesse de rotation azimutale comprise entre 0,2 tours par minute et 5 tours par minute ; la vitesse de rotation de l'antenne du radar est modifiée en fonction du temps de coupure inférieur Tc, de façon à ajuster le temps de résidence du radar sur la cible pour augmenter la valeur de Tc'.
L'invention concerne également : - un programme d'ordinateur comprenant des instructions de code de programme pour la mise en œuvre d'une ou plusieurs, de préférence toutes les étapes c) à j) d'un procédé selon l'invention, lorsque ledit programme est exécuté par un ordinateur ;
- un circuit électronique intégrant des fonctions réalisant une ou plusieurs, de préférence toutes les étapes c) à j), par exemple une intégration dans un circuit de type, en anglais FPGA (Field Programmable Gâte Array) ou DPS (Digital Signal Processing) ;
- un support informatique sur lequel est enregistré un tel programme, par exemple une mémoire ou un CD-ROM.
L'invention a également pour objet un centre comportant un ordinateur contenant un programme d'ordinateur selon l'invention, et un radar précipitations, de préférence plusieurs radars précipitations, connecté(s) à l'ordinateur, de manière à lui transmettre des informations relatives à des signaux émis et/ou reçus par le(s) dit(s) radar(s).
Le centre peut être fixe ou mobile, par exemple être un véhicule terrestre, un aéronef ou un engin spatial.
Définitions
- Le temps de résidence d'un radar sur une cible, en anglais « dwell time », désigne le temps d'éclairement de la cible par le radar, par exemple, lorsque l'antenne du radar est tournante, le temps pendant lequel le faisceau radar passe sur la cible.
- On dit que les impulsions du signal radar émis se succèdent suivant un rythme régulier quand le signal radar est constitué par la répétition continue d'un motif constitué de plusieurs impulsions.
- Le temps de cohérence T cible de la cible est classiquement l'intervalle de temps au bout duquel le coefficient d'autocorrélation d'un signal reçu de la cible est divisé par deux, par rapport à sa valeur au temps nul.
- On appelle classiquement « valeur complexe d'un signal » un signal ayant la forme A.exp(icp) avec A l'amplitude du signal, φ sa phase, exp() la fonction exponentielle et i le nombre imaginaire tel que i2 = -1.
Brève description des figures
L'invention pourra être mieux comprise à la lecture de la description qui va suivre, d'exemples non limitatifs de mise en œuvre de celle-ci, et à l'examen du dessin annexé, sur lequel : la figure 1 illustre schématiquement le fonctionnement d'un radar selon l'invention ; la figure 2 illustre des exemples d'impulsions du signal émis par le radar ; la figure 3 illustre le spectre de puissance d'un signal reçu ; - la figure 4 représente la fonction d'autocorrélation d'un signal reçu ; et la figure 5 est un schéma bloc résumant les différentes étapes du procédé selon l'invention.
Description détaillée
La figure 1 illustre schématiquement l'étape a) du procédé selon l'invention où un radar précipitation 10 émet vers une cible 12, en l'occurrence une précipitation atmosphérique, par exemple un nuage de pluie, de neige ou de grêle N, un signal Se constitué d'une série d'impulsions émises 14e. Deux impulsions émises consécutives quelconques sont séparées par un intervalle de base temporel choisi dans un groupe constitué d'au moins un intervalle de base constant, les impulsions se succédant de préférence suivant un rythme régulier.
On a représenté à la figure 2a le cas simple PRT (« Puise Répétition Time » en anglais), l'intervalle de temps entre deux impulsions émises consécutives quelconques est toujours égal au même intervalle de base AU.
Pour le cas double PRT, représenté à la figure 2b, l'intervalle de temps entre deux impulsions consécutives quelconques est choisi dans un groupe constitué de deux intervalles de base AU et Δί2.
Pour le cas triple PRT, représenté à la figure 2c, l'intervalle de temps entre deux impulsions consécutives quelconques est choisi dans un groupe constitué de trois intervalles de base AU,
Les exemples représentés aux figures 2b et 2c illustrent des impulsions se succédant suivant un rythme régulier, la série d'intervalles de base (AU, Δί2) et (AU, Δί2, Δί3) respectivement se répétant en permanence.
Bien entendu, pour le cas simple PRT, les impulsions émises se succèdent également suivant un rythme régulier. Chaque intervalle de base est de préférence compris entre 100 et 5 ms.
Le signal émis Se a de préférence une longueur d'onde constante, de préférence inférieure à
10 cm, notamment inférieure à 5 cm.
Le signal émis Se est de préférence envoyé suivant une direction d'émission Δ formant un angle β, inférieur à 15°, à 10°, de préférence inférieur à 5° avec un plan horizontal 16 sur lequel repose le radar 10.
Le radar 10 comporte de préférence une antenne parabolique 18 et un émetteur 20 disposé au foyer de l'antenne. Alternativement, le radar 10 peut comporter une antenne plate, de dimensions centimétriques.
L'antenne 18 focalise l'énergie électromagnétique dans un faisceau 22 dont la largeur angulaire a est inversement proportionnelle au diamètre de l'antenne. Le diamètre de l'antenne est de préférence compris entre 2m pour la bande X (9 MHz, λ = 3,3cm) et 8m pour la bande S (2,5 GHz, longueur d'onde λ = 12cm), un diamètre de 4,5m étant bien adapté pour la bande C (5 GHz, λ = 6cm). Ces fréquences sont imposées par la taille des objets renvoyant le signal, par exemple des gouttes d'eau constituant la précipitation, et par la réglementation.
La largeur angulaire a du faisceau 22 est de préférence supérieure à 0,5° et inférieure à 5°, de préférence inférieure à 3°, à 2°, mieux inférieure à 1°. La dimension horizontale et/ou verticale du champ d'observation est de préférence comprise entre 800m à 50km, mieux entre 4000m et 250km.
L'antenne 18 tourne à une vitesse de rotation azimutale de préférence comprise entre 0,2 tours par minute et 5 tours par minute.
L'étape b) du procédé selon l'invention est aussi illustrée par les figures 1 et 3. Le radar reçoit, à un instant to, un signal reçu Sr renvoyé par les gouttes d'eau du nuage N en réponse au signal émis Se.
La figure 3 représente l'amplitude A de puissance du signal reçu en échelle logarithmique en fonction de la vitesse Doppler v de ce signal, normalisée par rapport à la vitesse maximale observable va.
La raie étroite S à la vitesse nulle correspond sensiblement à la contribution de l'écho du sol au signal reçu. La bande large P centrée autour de la vitesse de 0,3 représente sensiblement la contribution de l'écho de la cible, en l'occurrence de l'écho du nuage N. La contribution du bruit additif B est sous la forme d'un continuum.
De manière connue, le signal reçu Sr est la somme de plusieurs signaux dont les temps de cohérence sont signifïcativement différents. Plus précisément, on distingue classiquement une composante cohérente, ou « partie cohérente », substantiellement constituée de l'écho des obstacles immobiles, et en particulier de l'écho de sol, et une composante incohérente ou « partie incohérente » substantiellement constituée de l'écho provenant de la cible atmosphérique.
Le reste des étapes du procédé est réalisé par un ordinateur 24 en communication avec le radar 10. Cet ordinateur 24 reçoit le signal reçu Sr et exécute les instructions d'un programme consistant à mettre en œuvre les étapes de traitement c) à j) du procédé selon l'invention.
A l'étape c), une série temporelle échantillonnée Sr(ti) du signal reçu est déterminée, Sr(ti) étant la valeur complexe du signal reçu Sr aux instants ti, qui sont les instants de réception des impulsions reçues 14r en réponse aux impulsions émises 14e.
A l'étape d), on détermine, à partir de la série temporelle Sr(ti), plusieurs coefficients d'autocorrélation R(Tj), Tj étant une durée telle qu'il existe au moins deux impulsions reçues 14r séparées par une durée égale à Tj :
R(Tj) = < Sr(ti ) . S (ti - Tj) >i, où * désigne le complexe conjugué, et < >i la moyenne arithmétique sur l'ensemble des produits (Sr(ti ) . Sr *(ti - Tj)) calculables à partir de ladite série temporelle Sr(ti).
Par « calculable à partir de la série temporelle », on entend qu'on ne considère que les instants ti pour lesquels une valeur de Sr(ti) et une valeur de Sr(ti - Tj) ont été déterminées à l'étape c).
Tj est égal à une somme de multiples des intervalles de base séparant les impulsions du signal émis: Tj = m.Ati + η2.Δΐ2 + η3.Δΐ3 + ni, n2 et n3 étant des nombres entiers ou nuls, au moins un des coefficients ni, n2 et n3 étant strictement positif.
Par exemple, pour le cas PRT simple, Tj e { Δίι, 2 AU, 3 AU, 4 AU,■■■ } .
Pour le cas double PRT, Tj e { Ati, At2, AU+ At2, 2AU+ At2, AU+ 2At2, 2AU+ 2At2, 3AU+ 2At2, ... } . l<j<J , La valeur maximale Tj est limitée par la durée de la série temporelle disponible. La figure 4 représente la fonction d'autocorrélation R du signal reçu Sr en fonction du temps T ou « lag ». Elle permet de lire les valeurs de R(Tj), Tj étant une valeur particulière de T . A titre d'exemple, les deux premiers coefficients R(Ati) et R(2 Δίι) ont été représentés.
La fonction d'autocorrélation est la transformée de Fourier inverse du spectre de puissance illustré à la figure 3.
La contribution de l'écho de sol se traduit par un continuum sur toute l'étendue de la fonction d'autocorrélation. La contribution de l'écho de la cible, par exemple de l'écho de pluie, au signal reçu Sr, ou « signal de cible », se traduit par une bande large autour de zéro. La raie étroite PB au temps zéro contient également la contribution du bruit, dont le temps de cohérence est pratiquement nul.
A l'étape e), on estime le temps de cohérence de la cible Tdbie. L'estimation peut être faite par une valeur habituelle pour la cible considérée. Par exemple, on peut estimer que le temps de cohérence d'une cible atmosphérique est typiquement compris entre 50ms et 500ms. Une autre méthode consiste à évaluer Tdbie à partir d'un signal Sr' voisin provenant d'une région voisine de la cible, mais pour lequel on est certain qu'il ne contient pas de composante cohérente. De préférence, le signal voisin est un signal renvoyé par la cible en réponse à un signal émis identique à Se, mais suivant une direction d'émission Δ' formant avec le plan horizontal un angle β' supérieur à β, c'est-à-dire émis vers une partie de la cible à une plus grande altitude. Avantageusement, la contribution du sol, et donc la composante cohérente du signal voisin, est négligeable dans un tel signal voisin. L'analyse du signal voisin Sr' permet, de manière classique d'évaluer un temps de cohérence Tdbie pour la cible.
On détermine ensuite un temps de coupure « inférieur » Tc supérieur à Tdbie, de préférence plus de deux fois supérieur à Tdbie. Pour les durées Tj supérieures au temps de coupure inférieur Tc, le coefficient d'autocorrélation R(Tj) est sensiblement identique au coefficient d'autocorrélation qui aurait été obtenu en n'analysant que la composante cohérente du signal reçu Sr. Autrement dit, pour ces coefficients d'autocorrélation R(Tj), la contribution de la composante incohérente du signal reçu Sr peut être considérée comme sensiblement nulle. Le temps de coupure est typiquement compris entre 50ms et 500ms. Dans certaines applications, il peut être cependant différent. Par exemple, il peut être dix fois inférieur si le radar émet autour de 95 GHz, ou bien dix fois supérieur si le radar émet autour de 50 MHz. A l'étape f), on détermine un temps de coupure « supérieur » Tc', plus grand que Tc et inférieur au temps de résidence du radar sur la cible.
La valeur de Tc' peut être quelconque, pourvu qu'elle reste inférieure audit temps de résidence. Plus Tc' est grand, meilleure sera la précision de l'estimation des coefficients d'autocorrélation supérieurs.
A l'étape g), l'analyse des coefficients d'autocorrélation R(Tj) pour des valeurs absolues de Tj supérieures au temps de coupure inférieur Tc, et inférieures au temps de coupure supérieur Tc', dits « coefficients supérieurs », permet de déterminer un signal qui correspond sensiblement à la composante cohérente du signal reçu Sr.
Les coefficients d'autocorrélation R(Tj) pour des valeurs absolues de Tj supérieures au temps de coupure inférieur Tc, et inférieures au temps de coupure supérieur Tc', dits « coefficients supérieurs » sont donc des coefficients d'autocorrélation qui auraient été obtenus si le signal reçu Sr avait été parfaitement cohérent, par exemple si la cible avait été fixe. A partir de ces coefficients, on peut déterminer une valeur représentative Rcoh de la composante cohérente du signal reçu.
Rcoh est de préférence une moyenne pondérée des coefficients supérieurs : RCOh =∑ kj . R(Tj), pour les Tj tels que Tc<Tj<Tc', avec∑ kj = 1.
De préférence, les coefficients de pondération kj diminuent à l'approche des temps de coupure Tc et Tc' et/ou augmentent à l'approche de la moyenne de ces deux temps de coupure, de sorte à suivre une évolution du type « fonction cloche » en fonction de Tj.
A l'étape h), on procède au traitement des coefficients d'autocorrélation R(Tj) pour des valeurs absolues de Tj inférieures au temps de cohérence de la cible Tdbie estimé à l'étape e), dits « coefficients inférieurs », en fonction de RCOh , de manière à obtenir des coefficients d'autocorrélation corrigés R(Tj) corr.
De préférence, pour obtenir les coefficients d'autocorrélation corrigés R(Tj)COrr, on soustrait Rcoh auxdits coefficients inférieurs : R(Tj)COrr = R(Tj) - RCOh, au moins pour des valeurs de Tj telles que |Tj| < Tdbie. A l'étape i), on procède à l'évaluation des caractéristiques de la cible à partir des coefficients d'autocorrélation corrigés R(Tj)Corr, de préférence suivant la méthode « puise pair », bien connue.
Par exemple, en simple PRT, la vitesse radiale vr de la cible, c'est-à-dire sa vitesse en direction du radar, est déterminée par la méthode « puise pair » à partir de la phase du premier coefficient d'autocorrélation R T i)corr du signal reçu:
vr = - λ/(4π Ti) arg(R(Ti)corr ), où λ est la longueur d'onde du signal radar émis.
Quant à la densité et l'homogénéité, elles peuvent être évaluées par les méthodes classiques en exploitant les coefficients d'autocorrélation R Tj)COrr.
De façon plus générale, une fois que le signal reçu est filtré par le procédé selon l'invention, on applique, à la dernière étape, une méthode connue pour calculer une ou plusieurs caractéristiques de la cible.
Comme cela apparaît clairement à présent, l'invention apporte une méthode simple et précise pour l'évaluation d'une ou plusieurs caractéristiques d'une cible. Des vérifications par simulations ont été effectuées qui montrent que les capacités de rejet des échos de sol par une telle méthode sont supérieures à 40 dB.
Bien entendu, l'invention n'est pas limitée aux exemples qui viennent d'être donnés, et peut être utilisée avec profit dans l'exploitation de radars d'autres natures : radar profileurs de vent, radars d'observation de l'océan, observation radar par satellite, etc.
En particulier, un procédé selon l'invention est utile pour filtrer l'effet de navires lorsque la cible est constituée de vagues.

Claims

REVENDICATIONS
1. Procédé d'évaluation d'une caractéristique d'une cible (12), choisie de préférence parmi la vitesse, la densité et l'homogénéité de la cible (12), le procédé comprenant les étapes suivantes : a) émission, vers ladite cible (12), au moyen d'un radar (10), pendant un temps de résidence, d'un signal émis (Se) constitué d'une série d'impulsions émises (14e), deux impulsions consécutives quelconques étant séparées par un intervalle de base temporel choisi dans un groupe constitué d'au moins un intervalle de base constant (Ah, Δί2, Δΐ3); b) réception d'un signal reçu (Sr) renvoyé par ladite cible (12) en réponse au signal émis (Se), le signal reçu (Sr) étant constitué d'une succession d'impulsions reçues (14r), à des instants ti, en réponse aux impulsions émises (14e); c) détermination d'une série temporelle échantillonnée Sr(ti) du signal reçu, Sr(ti) étant une valeur complexe du signal reçu, à l'instant t; ; d) détermination, à partir de la série temporelle Sr(ti), de plusieurs coefficients d'autocorrélation R(Tj), Tj étant tel qu'il existe au moins deux impulsions du signal reçu séparées par une durée égale à Tj :
R(Tj) = < Sr(ti ) . Sr *(ti - Tj) >i, où * désigne le complexe conjugué, et < >i la moyenne arithmétique sur l'ensemble des produits (Sr(ti ) . Sr *(ti - Tj)) calculables à partir de ladite série temporelle ; e) détermination d'un temps de coupure inférieur Tc positif, supérieur à une estimation du temps de cohérence de la cible TCibie ; f) détermination d'un temps de coupure supérieur Tc' supérieur à Tc et inférieur ou égal au temps de résidence ; g) analyse des coefficients d'autocorrélation R(Tj) pour des valeurs absolues de Tj supérieures au temps de coupure inférieur Tc et inférieures au temps de coupure supérieur Tc', dits « coefficients supérieurs », de manière à déterminer une valeur
RcohJ h) traitement des coefficients d'autocorrélation R(Tj) pour des valeurs absolues de Tj inférieures au temps de coupure inférieur Tc, dits « coefficients inférieurs », en fonction de RCOh , de manière à obtenir des coefficients d'autocorrélation corrigés
R(Tj)corrJ i) évaluation de ladite caractéristique de ladite cible en fonction des coefficients d'autocorrélation corrigés, de préférence suivant la méthode « puise pair ».
2. Procédé selon la revendication précédente, dans lequel le temps de coupure inférieur Tc est choisi entre 50ms et 500ms, de préférence en fonction de la nature de la cible (12).
3. Procédé selon l'une quelconque des revendications précédentes, dans lequel le temps de coupure inférieur Tc est estimé à partir d'un signal voisin du signal reçu Sr, un signal voisin étant renvoyé par la cible en réponse à un signal identique au signal émis (Se), mais suivant une direction d'émission (Δ') formant avec un plan horizontal (16) un angle (β') supérieur à l'angle (β ) que forme le signal émis (Se) avec ce plan.
4. Procédé selon l'une quelconque des deux revendications précédentes, dans lequel, à l'étape g), RCOh est une moyenne pondérée desdits coefficients supérieurs.
5. Procédé selon l'une quelconque des revendications précédentes, dans lequel, à l'étape h), on soustrait RCOh auxdits coefficients inférieurs : R(Tj)COrr = R(Tj) - RCOh.
6. Procédé selon l'une quelconque des revendications précédentes, dans lequel, à l'étape a), l'intervalle de temps entre deux impulsions consécutives quelconques est toujours égal au même intervalle de base (Δίι).
7. Procédé selon l'une quelconque des revendications 1 à 4, dans lequel, à l'étape a), l'intervalle de temps entre deux impulsions consécutives quelconques est choisi dans un groupe constitué de deux (Δίι, Δί2) ou trois intervalles de base (Δίι, Δί¾ Δί3).
8. Procédé selon l'une quelconque des revendications précédentes, dans lequel la cible est une couche de l'atmosphère et/ou dans lequel, à l'étape a), les impulsions sont envoyées suivant une direction formant un angle inférieur à 5° avec un plan horizontal sur lequel repose ledit radar (10).
9. Programme d'ordinateur comprenant des instructions de code de programme pour la mise en œuvre d'une ou plusieurs, de préférence toutes les étapes c) à j) d'un procédé selon l'une quelconque des revendications précédentes, lorsque ledit programme est exécuté par un ordinateur.
10. Centre comportant un ordinateur contenant un programme d'ordinateur selon la revendication précédente, et un radar précipitations (10), de préférence plusieurs radars précipitations (10), connecté(s) à l'ordinateur, de manière à lui transmettre des informations relatives à des signaux émis et/ou reçus par le(s) dit(s) radar(s) (10).
PCT/EP2018/062628 2017-05-17 2018-05-15 Procede d'evaluation de caracteristiques d'une cible par filtrage d'echos radar WO2018210887A1 (fr)

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
FR1754348 2017-05-17
FR1754348A FR3066608B1 (fr) 2017-05-17 2017-05-17 Procede d'evaluation de caracteristiques d'une cible par filtrage d'echos radar

Publications (1)

Publication Number Publication Date
WO2018210887A1 true WO2018210887A1 (fr) 2018-11-22

Family

ID=60080876

Family Applications (1)

Application Number Title Priority Date Filing Date
PCT/EP2018/062628 WO2018210887A1 (fr) 2017-05-17 2018-05-15 Procede d'evaluation de caracteristiques d'une cible par filtrage d'echos radar

Country Status (2)

Country Link
FR (1) FR3066608B1 (fr)
WO (1) WO2018210887A1 (fr)

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP1577680A2 (fr) * 2004-03-16 2005-09-21 Fujitsu Ten Limited Dispositif de radar
EP2251711A1 (fr) * 2009-05-05 2010-11-17 Meteo France Procédé et module de détermination de la vitesse d'au moins une cible par un radar météorologique

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP1577680A2 (fr) * 2004-03-16 2005-09-21 Fujitsu Ten Limited Dispositif de radar
EP2251711A1 (fr) * 2009-05-05 2010-11-17 Meteo France Procédé et module de détermination de la vitesse d'au moins une cible par un radar météorologique

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
CHORNOBOY E S ET AL: "Variable-PRI processing for meteorologic Doppler radars", RADAR DEVELOPMENTS IN A CHANGING ECONOMY. RECORD OF THE NATIONAL RADAR CONFERENCE. ATLANTA, MAR. 29 - 31, 19; [PROCEEDINGS OF THE NATIONAL RADAR CONFERENCE], NEW YORK, IEEE, US, 29 March 1994 (1994-03-29), pages 85 - 90, XP010117838, ISBN: 978-0-7803-1438-2, DOI: 10.1109/NRC.1994.328103 *
SACHIDANANDA M ET AL: "CLUTTER FILTERING AND SPECTRAL MOMENT ESTIMATION FOR DOPPLER WEATHER RADARS USING STAGGERED PULSE REPETITION TIME (PRT)", JOURNAL OF ATMOSPHERIC AND OCEANIC TECHNO, AMERICAN METEOROLOGICAL SOCIETY, BOSTON, MA, US, vol. 17, no. 3, 1 March 2000 (2000-03-01), pages 323 - 331, XP009016655, ISSN: 0739-0572, DOI: 10.1175/1520-0426(2000)017<0323:CFASME>2.0.CO;2 *

Also Published As

Publication number Publication date
FR3066608A1 (fr) 2018-11-23
FR3066608B1 (fr) 2019-07-12

Similar Documents

Publication Publication Date Title
CA2100107C (fr) Procede et dispositif pour l&#39;evaluation des precipitations sur une zone de terrain
EP2100161B1 (fr) Procede de traitement radar passif multivoies d&#39;un signal d&#39;opportunite en fm
EP3504559B1 (fr) Methode de traitement d&#39;un signal issu d&#39;un lidar coherent et systeme lidar associe
FR2833356A1 (fr) Radar a ondes de surface
FR2759462A1 (fr) Systemes radar pour l&#39;observation de phenomenes meteorologiques fournissant des donnees d&#39;observation ameliorees
EP3811109B1 (fr) Procédé de mesure de la hauteur de vagues à l&#39;aide d&#39;un radar aéroporté
EP2574957A1 (fr) Procédé d&#39;estimation de la fréquence doppler non ambigue d&#39;une cible mobile, notamment marine, et radar mettant en oeuvre le procédé
FR2931952A1 (fr) Procede d&#39;elimination des echos de sol pour un radar meteorologique
EP2455778B1 (fr) Procédé d&#39;estimation de la position angulaire d&#39;une cible par détection radar et radar mettant en oeuvre le procédé
EP0884605A1 (fr) Procédé de traitement du signal de récéption d&#39;un radar SAR de type DERAMP
EP2453251B1 (fr) Procédé pour réaliser une analyse haute résolution d&#39;une zone de l&#39;espace au moyen d&#39;une onde pulsée agile en fréquence
FR2971057A1 (fr) Procede d&#39;affichage de la surface equivalente radar minimale detectable
EP2544020B1 (fr) Procédé et dispositif de détection d&#39;une cible masquée par des réflecteurs de forte énergie
FR3066608B1 (fr) Procede d&#39;evaluation de caracteristiques d&#39;une cible par filtrage d&#39;echos radar
EP3538916B1 (fr) Procédé de contrôle de la compatibilité électromagnétique d&#39;un détecteur de radars avec au moins un émetteur de bord de signaux impulsionnels
EP2796893B1 (fr) Procédé d&#39;élimination d&#39;échos parasites en imagerie SAR
CN114167419A (zh) 一种结合卫星遥感图像与河流计数据的河宽提取的方法
EP3605145B1 (fr) Procédé de traitement distance haute résolution
FR2550347A1 (fr) Perfectionnements aux radars doppler a impulsions
EP3538917B1 (fr) Procédé de contrôle de la compatibilité électromagnétique d&#39;un détecteur de radars avec au moins un émetteur de bord de signaux impulsionnels
WO2023174774A1 (fr) Système radar et procédé associé pour optimiser l&#39;élimination de fouillis radar
EP3538919B1 (fr) Procédé de contrôle de la compatibilité électromagnétique d&#39;un détecteur de radars avec au moins un émetteur de bord de signaux impulsionnels
FR3027698A1 (fr) Procede de traitement du signal d&#39;un radar en mode de detection terre/mer ; chaine de traitement et produit programme d&#39;ordinateur associes
FR2992068A1 (fr) Procede de filtrage des echos ambigus d images sar
EP2687863B1 (fr) Procédé de fitrage d&#39;image radar polluée par des échos ambigus en distance

Legal Events

Date Code Title Description
121 Ep: the epo has been informed by wipo that ep was designated in this application

Ref document number: 18722641

Country of ref document: EP

Kind code of ref document: A1

NENP Non-entry into the national phase

Ref country code: DE

122 Ep: pct application non-entry in european phase

Ref document number: 18722641

Country of ref document: EP

Kind code of ref document: A1