FR3066608B1 - 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 PDFInfo
- Publication number
- FR3066608B1 FR3066608B1 FR1754348A FR1754348A FR3066608B1 FR 3066608 B1 FR3066608 B1 FR 3066608B1 FR 1754348 A FR1754348 A FR 1754348A FR 1754348 A FR1754348 A FR 1754348A FR 3066608 B1 FR3066608 B1 FR 3066608B1
- Authority
- FR
- France
- Prior art keywords
- time
- target
- signal
- coefficients
- radar
- Prior art date
- Legal status (The legal status 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 status listed.)
- Active
Links
- 238000000034 method Methods 0.000 title claims abstract description 44
- 238000001914 filtration Methods 0.000 title description 3
- 238000001556 precipitation Methods 0.000 claims description 17
- 230000004044 response Effects 0.000 claims description 8
- 238000004590 computer program Methods 0.000 claims description 4
- GNFTZDOKVXKIBK-UHFFFAOYSA-N 3-(2-methoxyethoxy)benzohydrazide Chemical compound COCCOC1=CC=CC(C(=O)NN)=C1 GNFTZDOKVXKIBK-UHFFFAOYSA-N 0.000 claims description 3
- 230000005540 biological transmission Effects 0.000 claims 1
- 238000012545 processing Methods 0.000 abstract description 4
- 230000001427 coherent effect Effects 0.000 description 11
- 230000006870 function Effects 0.000 description 8
- 238000013459 approach Methods 0.000 description 4
- 238000005311 autocorrelation function Methods 0.000 description 4
- 230000033764 rhythmic process Effects 0.000 description 3
- FGUUSXIOTUKUDN-IBGZPJMESA-N C1(=CC=CC=C1)N1C2=C(NC([C@H](C1)NC=1OC(=NN=1)C1=CC=CC=C1)=O)C=CC=C2 Chemical compound C1(=CC=CC=C1)N1C2=C(NC([C@H](C1)NC=1OC(=NN=1)C1=CC=CC=C1)=O)C=CC=C2 FGUUSXIOTUKUDN-IBGZPJMESA-N 0.000 description 2
- 239000005442 atmospheric precipitation Substances 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 238000001228 spectrum Methods 0.000 description 2
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 2
- 101100495256 Caenorhabditis elegans mat-3 gene Proteins 0.000 description 1
- 239000000654 additive Substances 0.000 description 1
- 230000000996 additive effect Effects 0.000 description 1
- 230000033228 biological regulation Effects 0.000 description 1
- 238000004891 communication Methods 0.000 description 1
- 238000007796 conventional method Methods 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 238000006073 displacement reaction Methods 0.000 description 1
- 238000002592 echocardiography Methods 0.000 description 1
- 238000005562 fading Methods 0.000 description 1
- 238000005286 illumination Methods 0.000 description 1
- 230000010354 integration Effects 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 230000002265 prevention Effects 0.000 description 1
- 238000004088 simulation Methods 0.000 description 1
- 239000002689 soil Substances 0.000 description 1
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
Landscapes
- Engineering & Computer Science (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Physics & Mathematics (AREA)
- Computer Networks & Wireless Communication (AREA)
- General Physics & Mathematics (AREA)
- Electromagnetism (AREA)
- Radar Systems Or Details Thereof (AREA)
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
Procédé d’évaluation de caractéristiques d’une cible par filtrage d’échos 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 fdtrer 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 fdtrage 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 = cÀ/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 ambigüe. 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 fdtrage 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 ti ; 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 :
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 Rcoh ; 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)corr ; 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, Rœh 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 ps 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. F’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).
Fe centre peut être fixe ou mobile, par exemple être un véhicule terrestre, un aéronef ou un engin spatial. Définitions - Fe 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. - Fe temps de cohérence TCibie 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 Α.βχρ(ίφ) 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 F’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 Δΐι.
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 Δΐι 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 Δΐι, Δΐ2 et Δΐ3.
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 (Δΐι, Δΐ2) et (Δΐι, Δΐ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 ps 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 significativement 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) = <SI(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 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 + mAt3 + ..., m, m et m étant des nombres entiers ou nuis, au moins un des coefficients ni, m et m étant strictement positif.
Par exemple, pour le cas PRT simple, Tj e { Ati, 2 Ati, 3 Ati, 4 Ati, ...
Pour le cas double PRT, Tj e { Ati, At2, Ati+ Δΐ2, 2Δΐι+ Δΐ2, Δίι+ 2Δΐ2, 2Δΐι+ 2Δΐ2, 3Δΐι+ 2Δί2, 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 Ati) 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 TCibie. 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 TCibie à 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 TCibie pour la cible.
On détermine ensuite un temps de coupure « inférieur » Tc supérieur à TCibie, de préférence plus de deux fois supérieur à TCibie. 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 TCibie 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| < TCibie. 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(Ti)COrr du signal reçu: vr = - λ/(4π Τι) 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 profiteurs 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 (10)
- REVENDICATIONS1. 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 (Δΐι, Δΐ2, Δΐΐ); 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 ti ; 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 :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 Rcoh) 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)corr) 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 Tune 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 Tune 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 Tune 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 Tune 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 Tune 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 (Δΐι, Δΐ2, Δΐΐ).
- 8. Procédé selon Tune 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 Tune 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).
Priority Applications (2)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| FR1754348A FR3066608B1 (fr) | 2017-05-17 | 2017-05-17 | Procede d'evaluation de caracteristiques d'une cible par filtrage d'echos radar |
| PCT/EP2018/062628 WO2018210887A1 (fr) | 2017-05-17 | 2018-05-15 | Procede d'evaluation de caracteristiques d'une cible par filtrage d'echos radar |
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 (2)
| Publication Number | Publication Date |
|---|---|
| FR3066608A1 FR3066608A1 (fr) | 2018-11-23 |
| FR3066608B1 true FR3066608B1 (fr) | 2019-07-12 |
Family
ID=60080876
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| FR1754348A Active FR3066608B1 (fr) | 2017-05-17 | 2017-05-17 | 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) |
Family Cites Families (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2005265461A (ja) * | 2004-03-16 | 2005-09-29 | Fujitsu Ten Ltd | レーダ装置 |
| FR2945352A1 (fr) * | 2009-05-05 | 2010-11-12 | Meteo France | Procede et module de determination de la vitesse d'au moins une cible par un radar meteorologique |
-
2017
- 2017-05-17 FR FR1754348A patent/FR3066608B1/fr active Active
-
2018
- 2018-05-15 WO PCT/EP2018/062628 patent/WO2018210887A1/fr active Application Filing
Also Published As
| Publication number | Publication date |
|---|---|
| WO2018210887A1 (fr) | 2018-11-22 |
| FR3066608A1 (fr) | 2018-11-23 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CA2100107C (fr) | Procede et dispositif pour l'evaluation des precipitations sur une zone de terrain | |
| EP3504559B1 (fr) | Methode de traitement d'un signal issu d'un lidar coherent et systeme lidar associe | |
| EP2574957B1 (fr) | Procédé d'estimation de la fréquence doppler non ambigue d'une cible mobile, notamment marine, et radar mettant en oeuvre le procédé | |
| CA2968487C (fr) | Procede de detection et de caracterisation par ultrasons de defauts dans un materiau heterogene | |
| FR2833356A1 (fr) | Radar a ondes de surface | |
| EP3811109B1 (fr) | Procédé de mesure de la hauteur de vagues à l'aide d'un radar aéroporté | |
| EP2455778B1 (fr) | Procédé d'estimation de la position angulaire d'une cible par détection radar et radar mettant en oeuvre le procédé | |
| FR2931952A1 (fr) | Procede d'elimination des echos de sol pour un radar meteorologique | |
| EP2453251B1 (fr) | Procédé pour réaliser une analyse haute résolution d'une zone de l'espace au moyen d'une onde pulsée agile en fréquence | |
| FR2971057A1 (fr) | Procede d'affichage de la surface equivalente radar minimale detectable | |
| EP2544020B1 (fr) | Procédé et dispositif de détection d'une cible masquée par des réflecteurs de forte énergie | |
| EP0901640B1 (fr) | Procede de modelisation fine du fouillis de sol recu par un radar | |
| FR3066608B1 (fr) | Procede d'evaluation de caracteristiques d'une cible par filtrage d'echos radar | |
| EP0786096B1 (fr) | Procede pour evaluer des effets non-euclidiens affectant une image obtenue avec un radar spatial, et satellite pour sa mise en uvre | |
| EP3538916B1 (fr) | Procédé de contrôle de la compatibilité électromagnétique d'un détecteur de radars avec au moins un émetteur de bord de signaux impulsionnels | |
| CN114167419A (zh) | 一种结合卫星遥感图像与河流计数据的河宽提取的方法 | |
| EP2796893A1 (fr) | Procédé d'élimination d'échos parasites en imagerie SAR | |
| FR3027698A1 (fr) | Procede de traitement du signal d'un radar en mode de detection terre/mer ; chaine de traitement et produit programme d'ordinateur associes | |
| EP4330722B1 (fr) | Système radar et procédé associé pour optimiser l'élimination de fouillis radar | |
| FR2550347A1 (fr) | Perfectionnements aux radars doppler a impulsions | |
| EP3538919B1 (fr) | Procédé de contrôle de la compatibilité électromagnétique d'un détecteur de radars avec au moins un émetteur de bord de signaux impulsionnels | |
| Mohammadloo et al. | An uncertainty assessment of the effect of using FM pulses on MBES depth measurements | |
| FR2992068A1 (fr) | Procede de filtrage des echos ambigus d images sar | |
| FR2972055A1 (fr) | Procede de determination d'un rang d'ambiguite en distance d'un signal recu par un radar doppler a impulsions | |
| EP2687863B1 (fr) | Procédé de fitrage d'image radar polluée par des échos ambigus en distance |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PLFP | Fee payment |
Year of fee payment: 2 |
|
| EXTE | Extension to a french territory |
Extension state: PF |
|
| PLSC | Publication of the preliminary search report |
Effective date: 20181123 |
|
| PLFP | Fee payment |
Year of fee payment: 4 |
|
| PLFP | Fee payment |
Year of fee payment: 5 |
|
| PLFP | Fee payment |
Year of fee payment: 6 |
|
| PLFP | Fee payment |
Year of fee payment: 7 |
|
| PLFP | Fee payment |
Year of fee payment: 8 |