Brevet FR2986620A1 — publié le 20130809

PROCEDE DE LEVEE D'AMBIGUÏTE SUR LES MESURES EN DISTANCE La présente invention concerne le domaine des radars Doppler à 5 impulsions et en particulier le traitement des signaux émis et reçus par de tels radars. La présente invention se rapporte à un procédé de levée des ambiguïtés distance. Un radar Doppler à impulsions permet de caractériser une cible 10 éclairée par le faisceau radar suivant deux paramètres, sa distance D par rapport au radar et sa vitesse de rapprochement Vr. Pour cela, le radar Doppler à impulsions émet une onde radioélectrique composée d'une série d'impulsions périodiques ou train d'impulsions. La distance D se déduit du retard t séparant l'instant d'émission d'une 15 impulsion et l'instant de réception de son écho rétro-diffusé par la cible. Le retard t est donné par la formle : t=2D/c, où c désigne la vitesse de la lumière. La vitesse de rapprochement Vr est estimée en émettant un train de N impulsions périodiques et en effectuant une analyse spectrale (par 20 transformée de Fourier) sur les N signaux reçus et démodulés par la fréquence émise. En effet, après démodulation, le signal constitué par les N échos rétro-diffusés par la cible est un signal périodique dont la fréquence, appelée fréquence Doppler fd, est proportionnelle à la vitesse de rapprochement de la cible. Cette fréquence Doppler est donnée par la 25 formule: fd=2Vr/X, où X désigne la longueur d'onde de radar. Du fait de la discrétisation de ce signal sur n échantillons avec une fréquence d'échantillonnage Fr correspondant à la fréquence de répétition des impulsions émises, seule la valeur « ambiguë » de la fréquence Doppler, 30 c'est à dire la valeur modulo Fr , peut être estimée. Pour travailler avec une valeur d'ambiguïté Doppler et donc d'ambiguïté vitesse la plus haute possible, on privilégie les formes d'onde avec une fréquence de répétition des impulsions, aussi appelée fréquence 35 de récurrence, élevée. À Cependant, cette augmentation de la fréquence de récurrence se fait au détriment de la mesure distance puisqu'elle rend à son tour cette mesure distance ambiguë : il est en effet impossible de faire la distinction entre un écho situé à une distance D et un autre situé à D+Damb, avec Damb=c/(2Fr). Dès lors que les signaux rétrodiffusés peuvent s'étendre sur une plage distance supérieure à la distance ambiguë, on parle de forme d'onde ambiguë en distance, la distance de la cible est seulement connue à un rang d'ambiguïté Distance près. En conséquence, pour un radar Doppler à impulsions, le choix de la fréquence de récurrence résulte d'un compromis. Une fréquence de récurrence relativement haute conduit à une valeur d'ambiguïté Doppler relativement haute, mais à une distance ambiguë relativement faible. A l'inverse, une fréquence de récurrence relativement basse conduit à une distance ambiguë relativement haute, mais à une valeur d'ambiguïté Doppler relativement basse. De façon à lever les ambiguïtés distance, il est connu dans l'art antérieur une méthode consistant à émettre successivement plusieurs trains d'impulsions en faisant varier la fréquence de récurrence (et donc la valeur de la distance ambiguë) d'un train à l'autre. L'ambiguïté est ensuite levée par « dépliement » en dupliquant pour chaque train d'impulsions l'image Distance-Doppler associée, de manière à couvrir entièrement le domaine Distance-Doppler non ambigu. Après superposition de ces K images Distance-Doppler " dépliées " issues des K trains d'impulsions, la position non ambiguë de la cible correspond à la position du pixel où la cible est présente sur les K images dépliées. Cette méthode est notamment décrite dans l'ouvrage "Air and Spaceborne Radar Systems, an introduction" au paragraphe 8.6.6. Un inconvénient de cette méthode est que les différents dépliements dupliquent également tous les échos radars parasites notamment les échos de retour de sol : le fouillis (ou clutter selon le terme anglo-saxon). Afin d'éviter cette duplication des échos parasites en distance, il convient de diminuer la fréquence de répétition des impulsions cependant du fait de l'interdépendance des paramètres radar ceci entraine une durée d'illumination du radar plus grande et donc diminue sa discrétion. Un but de l'invention est notamment de corriger les inconvénients 5 précités en proposant une méthode de levée d'ambigüité distance des cibles limitant les repliements des échos parasites en distance et compatible avec l'utilisation de fréquences de récurrence basses. A cet effet, l'invention a pour objet un procédé de levée d'ambiguïté 10 sur les mesures en distance mis en oeuvre par un système radar utilisant au moins un radar Doppler à impulsions émettant un signal comportant deux trains de N impulsions périodiques, et recevant les signaux rétrodiffusés par une cible mobile, lesdites périodes de récurrence des deux trains d'impulsions émis par le radar étant différentes et le procédé comportant : 15 - Une étape de séparation des deux trains d'impulsions du signal reçu en deux sous-signaux comportant chacun un train de N impulsions, - Une étape de suppression de la contribution des échos parasites du fouillis sur chacun des sous-signaux, 20 - Une étape de test successif des différentes hypothèses de rang d'ambiguïté possible. Selon une particularité de l'invention, la différence ATr entre les périodes de récurrence des deux trains d'impulsions est égale à un nombre entier !ceci strictement supérieur à 1 de case distance et est inférieure à 1/(2N 25 Fdmax) Où Fdmax représente la fréquence Doppler maximale envisagée. Selon une autre particularité l'étape de suppression de la contribution des échos parasites du fouillis comporte - une étape de filtrage simple annulation, - une étape de transposition dans le domaine fréquentiel, 30 - une étape de compensation du filtrage simple annulation, - une étape d'effacement des échos parasites du fouillis, - une étape de transposition dans le domaine temporel. De façon avantageuse, la transposition dans le domaine fréquentiel est réalisée par transformée de Fourier rapide et en ce que la transposition dans le domaine temporel est réalisée par transformée de Fourier rapide inverse. Suivant un mode de mise en oeuvre, le test des différentes hypothèses de rang d'ambiguïté possible comporte pour chaque rang d'ambiguïté testé : - une étape de regroupement des deux sous-signaux, - une étape de pondération du signal regroupé, - une étape de transposition dans le domaine fréquentiel, - une étape de conversion en décibel, - une étape de filtrage des échos issus des rangs d'ambiguïté différents du rang d'ambiguïté sous test. Selon une variante de mise en oeuvre du procédé, l'étape de regroupement pour un rang d'ambiguïté r testé correspond au regroupement 15 des deux sous-signaux après décalage temporel d'un des deux sous-signaux d'une valeur de r.kd cases distances. Selon une autre particularité, l'étape de filtrage s'effectue en calculant la valeur minimale entre un premier signal correspondant à l'interpolation sur 2N points du premier sous-signal, un deuxième signal correspondant à 20 l'interpolation sur 2N points du deuxième sous-signal après décalage suivant l'axe case distance d'une valeur proportionnelle au rang sous test et un troisième signal correspondant au spectre du signal regroupé du rang sous test. Avantageusement, l'interpolation des sous-signaux après suppression 25 des échos parasites du fouillis sur 2N points comprend : une étape (51) de pondération, une étape (52) d'extension par zero-padding (extension par des zéros) une étape (53) de transposition dans le domaine fréquentiel, 30 - une étape (54) de conversion en décibel. Selon une autre particularité, la transposition dans le domaine fréquentiel est réalisée par transformée de Fourier rapide. Un avantage de la présente invention est que le procédé peut être mis en oeuvre sur un système radar classique sans nécessiter d'évolutions importantes du matériel. D'autres particularités et avantages de la présente invention apparaîtront plus clairement à la lecture de la description ci-après, donnée à titre illustratif et non limitatif, et faite en référence aux dessins annexés, dans lesquels : La figure 1 représente des chronogrammes d'exemple d'impulsions radar ; La figure 2 représente un exemple d'organigramme représentant les étapes possibles du procédé selon l'invention ; 15 La figure 3 représente un exemple d'algorithme complet de levée d'ambiguïté distance selon l'invention ; La figure 4 représente un exemple d'algorithme de l'étape de suppression de la contribution des échos parasites du fouillis ; La figure 5 représente un exemple d'organigramme de l'étape 20 d'interpolation des sous-signaux ; La figure 6 représente, par un organigramme, un exemple de calcul de la carte sous l'hypothèse de rang d'ambiguïté 0 ; La figure 6 bis représente un exemple de signal obtenu après regroupement des deux sous-signaux sous l'hypothèse d'un rang 25 d'ambiguïté 0 avant et après passage dans le domaine fréquentiel pour une cible de rang d'ambiguïté 0 ; La figure 7 représente, par un organigramme, un exemple de calcul de la carte sous l'hypothèse de rang d'ambiguïté 1 ; La figure 7 bis représente un exemple de signal obtenu après 30 regroupement des deux sous-signaux sous l'hypothèse d'un rang d'ambiguïté 1 avant et après passage dans le domaine fréquentiel pour une cible de rang d'ambiguïté 0 ; 10 La figure 8 représente l'étape de filtrage du rang d'ambiguïté 0; La figure 8 bis représente les amplitudes des signaux des cartes Tz_b1, Tz_b2, Tr_a1 et de la valeur minimale det_a1 en fonction de la case distance; La figure 9 représente l'étape de filtrage du rang d'ambiguïté 1 ; La figure 9 bis représente les amplitudes des signaux des cartes Tz_b1, Tz_b2, Tr_a2 et de la valeur minimale det_a2 en fonction de la case distance. En référence à la figure 1, afin d'illustrer le principe de l'invention, on considère une forme d'onde cohérente. On rappelle qu'un radar est dit cohérent en phase lorsque ce radar maîtrise la différence de phase entre deux signaux rétro-diffusés par une même cible, c'est à dire qu'il sait en extraire le terme relatif à cette cible. Ce type de radar permet d'exploiter la phase des signaux reçus, en plus de la puissance rétrodiffusée, typiquement pour pouvoir estimer une fréquence Doppler par transformée de Fourier. On parle alors de traitement cohérent et d'intégration cohérente. Le premier chronogramme 11 représente un exemple de signal électromagnétique pouvant être émis par un radar Doppler à impulsions mettant en oeuvre l'invention. Suivant une particularité de l'invention, la forme d'onde utilisée est formée de deux trains de N impulsions identiques périodiques. Les impulsions du premier train sont émises avec une période de répétition ou période de récurrence Tri. Le second train d'impulsions est émis avec une période de récurrence Tr2 différente de Tri telle que Tr2=Tr1-FATr. De façon avantageuse, la différence \Tr entre les deux périodes de récurrence correspond à un nombre entier strictement supérieur à 1 de case distance. Chaque impulsion dans le train est identifiée par un numéro de récurrence n avec n E [0, 2N-1]. On considère le cas de deux cibles situées à des distances Di et D2 30 du radar dans deux rangs d'ambiguïté distance différents. Les chronogrammes 12 et 13 représentent des exemples d'échos reçus par le radar après réflexion de l'onde émise par le radar sur ces cibles. Dans l'exemple représenté par le chronogramme 12, la cible est située à une distance Di du radar et est non ambigüe, c'est-à-dire que le retard induit par cette distance est inférieur à la période de récurrence Tri. Dans ce cas, on parle d'ambiguïté distance de rang zéro. Dans ce cas de figure, chaque impulsion de l'onde émise par le radar donne lieu à un écho qui est reçu avant l'émission de l'impulsion suivante. Le 5 décalage temporel t1 par rapport à la dernière impulsion émise est donc indépendant de la période de récurrence. Le troisième chronogramme 13 représente le cas d'une cible ambigüe en distance. L'écho de l'impulsion rétrodiffusée par la cible est reçu par le radar après l'émission d'une nouvelle impulsion. Dans cet exemple, le 10 premier écho est reçu pendant la période de récurrence correspondant à la deuxième impulsion (n=1 cf Figure 2), son rang d'ambiguïté est égal à un. Contrairement au cas précédent, dans le cas d'une cible ambigüe en distance, le décalage temporel (de valeur t21 ou t22) dépend de la période de récurrence et du rang d'ambiguïté. 15 L'invention utilise de façon avantageuse cette différence de comportement entre une cible non ambiguë et une cible ambiguë pour déterminer le rang d'ambiguïté. Pour la suite, on supposera que le traitement est mono-cible. 20 En référence à la figure 2, les différentes étapes possibles du procédé de levée d'ambiguïté sur les mesures en distance dans le traitement d'un signal d'écho reçu par un radar Doppler à impulsions sont présentées. On suppose que le signal émis par le radar Doppler à impulsions comporte 25 deux trains de N impulsions périodiques chacun. Suivant une particularité de l'invention, les périodes de récurrence Tr puis Tr+ATr de ces deux trains d'impulsions sont différentes. La valeur de ATr est choisie de façon à ce que le nombre de cases distances kid correspondant soit entier et supérieur à 1, c'est à dire 30 3,Tr=kcd.Tech où Tech représente la période d'échantillonnage sur l'axe distance. Pour la suite, de façon à simplifier les explications et la représentation, le nombre de case distance kcd sera pris égal à 2. Cet exemple est nullement limitatif et une extension à un nombre de case distance supérieur est possible avec le même raisonnement La valeur de ATr est également choisie de façon à ce que le déphasage pour la fréquence Doppler maximale fdmax envisagée soit suffisament faible pour permettre une intégration cohérente. Le signal doit rester dans la même case fréquence. Cela peut se traduire par : ATr A la réception, le signal rétrodiffusé par une cible est séparé en deux au cours d'une étape 21 de séparation des deux trains d'impulsions du signal reçu en deux sous-signaux comportant chacun un train de N impulsions. La contribution des échos parasites du fouillis est ensuite supprimée sur chacun des sous-signaux dans une étape 22. Une étape 23 de test successif des différentes hypothèses de rang d'ambiguïté possible est ensuite opérée afin de déterminer le rang d'ambiguïté du signal reçu. La figure 3 présente un exemple d'algorithme complet de levée d'ambiguïté distance selon l'invention. De façon à ne pas surcharger la figure, l'algorithme est présenté à travers un traitement sur deux rangs d'ambiguïtés. Cet exemple n'est nullement limitatif et peut s'étendre aux cas de rangs d'ambiguïtés distance supérieurs. L'invention consiste à balayer un à un l'ensemble des rangs d'ambiguïté 25 distance possible et à filtrer à chaque fois les échos issus des rangs d'ambiguïtés différents du rang d'ambiguïté sous test. En référence à la figure 4, l'étape 22 de suppression de la contribution des échos parasites du fouillis est présentée. Cette étape 22 a pour objectif 30 de préparer chaque sous-signal avant leur regroupement. Elle est réalisée séparément sur les deux sous-signaux issus du signal reçu. Sur chacun de ces sous-signaux, un filtrage simple annulation 221, une transposition dans le domaine fréquentiel 222 et une compensation 223 du filtrage simple annulation sont appliqués. Suivant un mode de mise en 35 oeuvre, la transposition dans le domaine fréquentiel est réalisée par transformée de Fourrier rapide ou FFT pour Fast Fourrier Transform selon la terminologie anglo-saxonne. Le filtre simple annulation (SA) et sa compensation (CSA) ont pour but d'atténuer les échos parasites du fouillis (ou clutter suivant le terme anglo - saxon) et en particulier les lobes dus au passage dans le domaine fréquentiel tout en conservant le signal utile et ce même après un retour dans le domaine temporel. De façon avantageuse, cette opération permet d'obtenir une meilleure détection des cibles. Les échos parasites du fouillis localisés autour du lobe principal sont ensuite supprimés au cours d'une étape 224 d'effacement des échos parasites du fouillis. De façon avantageuse, la suppression de ces échos parasites permet d'éviter la formation de forts lobes secondaires après regroupement des deux sous-signaux. Ces lobes sont générés par le regroupement de retour de sol de deux cases distances différentes. Une transposition dans le domaine temporel 225 du signal obtenu est ensuite opérée. Cette transposition peut être réalisée par transformée de Fourier inverse (ou iFFT pour inverse Fast Fourier Transform). Dans le cas où le procédé selon l'invention est mis en oeuvre par un radar aéroporté, une étape supplémentaire de compensation de phase sera 20 ajoutée avant l'étape de filtrage simple annulation 221. La figure 5 présente un exemple d'organigramme de l'étape d'interpolation des sous-signaux. Cette opération à pour but de préparer les sous-signaux à N impulsions pour 25 pouvoir les comparer au signal à 2N impulsions calculé après le regroupement des deux sous-signaux. Cette interpolation est réalisée en complétant le signal par N zéros supplémentaires (technique de zeropadding selon le terme anglo-saxon). Le traitement comporte une étape 51 de pondération, par exemple par un filtre de Hanning ou de Tchebychev, une 30 étape 52 d'extension du signal obtenu par des zéros (technique de zeropadding), une étape 53 de transposition dans le domaine fréquentiel, par exemple par l'intermédiaire d'une transformée de Fourrier rapide et se termine par une étape 54 de conversion en décibel (PdB). La figure 6 présente un exemple de calcul du test sous l'hypothèse de rang d'ambiguïté 0. La première étape consiste à regrouper 231 les deux sous-signaux après suppression des échos du fouillis. L'hypothèse étant que le rang d'ambiguïté est nul, le deuxième sous-signal n'est pas décalé. Le signal regroupé est ensuite pondéré 232 et subit une transposition dans le domaine fréquentiel 233 par exemple par une transformée de Fourrier rapide. Le signal obtenu est ensuite converti en décibel 234. La figure 7 présente un exemple de calcul du test sous l'hypothèse de rang d'ambiguïté 1. Sous l'hypothèse que la cible possède un rang d'ambiguïté de 1 la case distance de cette cible dans le premier sous-signal est égale à sa case distance dans le second sous-signal plus kid cases. Pour regrouper les contributions de ces deux sous-signaux, il suffit de regrouper le signal de la case distance du premier sous-signal avec le signal de la case distance du second sous-signal décalé de kod cases distances. Le traitement consiste à décaler les cases distances du second sous-signal 230 de kcd cases, à regrouper les récurrences des deux sous-signaux 231, à 20 effectuer une pondération 232, une transposition dans le domaine fréquentiel 233 puis à convertir le signal obtenu en décibel 234. Comme précédemment, si les cibles recherchées ne sont pas ambigües en vitesse, il est possible de compenser le décalage Doppler des impulsions du deuxième sous-signal avant la transposition en fréquence de 25 façon à annuler tout élargissement de lobe. Dans cet exemple nullement limitatif, le nombre d'hypothèses de rang d'ambigüité est pris égal à 2. Dans un cas général où le nombre d'hypothèse serait supérieur, les branches de la Figure 7 sont dupliquées pour chaque rang d'ambiguïté supplémentaire, le décalage des récurrences du deuxième 30 sous-signal, pour un rang r, serait de r.lccd cases distance. En référence aux figures 8 et 9 les étapes de filtrage respectivement du rang d'ambigüité 1 et du rang d'ambigüité 0 vont être présentées. Ce filtrage repose sur le fait que si une cible de rang d'ambiguïté 0 est 35 située dans la case distance d'indice icd dans la carte Tr_a1 de rang d'ambiguïté 0, elle sera située dans la même case dans la carte Tz_b1 issue du premier sous-signal et dans la carte Tz_b2 issue du second sous-signal. Ceci est illustré figure 6 bis où est représenté un exemple de signal obtenu après regroupement des deux sous-signaux et après suppression de l'influence du fouillis sur chaque sous-signal sous l'hypothèse d'un rang d'ambiguïté 0, avant et après passage dans le domaine fréquentiel, pour une cible de rang d'ambiguïté 0. Dans cet exemple, un seul signal est présent au niveau de la case distance 500. Par contre, une cible de rang d'ambiguïté 1 sera partagée sur deux cases distances, ica et icd -kid, dans la carte Tr a1 d'ambiguïté 0. Sur la carte Tz_b1 issue du premier sous-signal elle sera concentrée dans la case icd et sur la carte Tz_b2 issue du second sous-signal elle sera concentrée dans la case icd -kcd. Pour filtrer les cibles d'ambiguïté 1 et conserver les cibles d'ambiguïté 0, le traitement consiste à calculer la valeur minimale det_a1 sur les signaux des trois cartes Tr_a1, Tz_b1 et Tz_b2. La Figure 8 bis présente le résultat (courbe delt_a1) du filtrage du rang d'ambiguïté 0 appliqué à une cible de rang d'ambiguïté 0. Le pic représentant la cible est bien conservé sans perte. Suivant un raisonnement analogue, une cible d'ambiguïté 0 sera partagée sur deux cases distances, icd et icd -kcd, dans la carte Tr_a2 d'ambiguïté 1. Ceci est illustré figure 7 bis où est représenté un exemple de signal obtenu après regroupement des deux sous-signaux, après suppression de l'influence du fouillis, sous l'hypothèse d'un rang d'ambiguïté 1, avant et après passage dans le domaine fréquentiel, pour une cible de rang d'ambiguïté 0. Sur cette représentation, le signal occupe la case distance 500 pour les 32 premiers échantillons et la case distance 502 pour les 32 échantillons suivants. Sur la carte Tz_b1 issue du premier sous-signal, la cible d'ambiguïté 0 sera située dans la case icd seule et sur la carte Tz_bZ issue du second sous-signal, elle sera également située dans la case icd. Si une cible de rang d'ambiguïté 1 est située dans la case distance d'indice icd dans la carte Tr_a2 de rang d'ambiguïté 1, elle sera située dans la même case dans la carte Tz_b1 issue du premier sous-signal, et dans la case distance icd -kid dans la carte Tz_b2 issue du second sous-signal. Pour filtrer les cibles d'ambiguïté 0 et conserver les cibles d'ambiguïté 1, le traitement consiste à calculer la valeur minimale det_a2 sur les 3 cartes Tr_a2, Tz_b1 et la carte Tz_b2 décalée de kcd cases distances. La Figure 9 bis présente le résultat (courbe delt_a2) du filtrage du rang d'ambiguïté 1 appliqué à une cible de rang d'ambiguïté O. Le pic représentant la cible de rang d'ambiguïté 0 est bien supprimé. Comme précédemment, dans un cas général où le nombre 10 d'hypothèses serait supérieur à 2, les branches de la figure 9 sont dupliquées avec un décalage 230 des récurrences du deuxième sous-signal d'une valeur de r.kcd cases distance, pour un rang d'ambiguité r. REVENDICATIONS1. Procédé de levée d'ambiguïté sur les mesures en distance mis en oeuvre par un système radar utilisant au moins un radar Doppler à impulsions émettant un signal comportant deux trains de N impulsions périodiques, et recevant les signaux rétrodiffusés par une cible mobile, caractérisé en ce que les périodes de récurrence des deux trains d'impulsions émis par le radar sont différentes et en ce que le procédé comporte : - Une étape (21) de séparation des deux trains d'impulsions du signal reçu en deux sous-signaux comportant chacun un train de N impulsions, - Une étape (22) de suppression de la contribution des échos parasites du fouillis sur chacun des sous-signaux, - Une étape (23) de test successif des différentes hypothèses de rang d'ambiguïté possible. 2. Procédé selon la revendication précédente caractérisé en ce que la différence ATr entre les périodes de récurrence des deux trains d'impulsions est égale à un nombre entier kcj strictement supérieur à 1 de case distance et est inférieure à 1/(2N Fdmax) où Fdmax représente la fréquence Doppler maximale envisagée. 3. Procédé selon une des revendications précédentes caractérisé en ce que l'étape (22) de suppression de la contribution des échos parasites du fouillis comporte - une étape (221) de filtrage simple annulation, - une étape (222) de transposition dans le domaine fréquentiel, - une étape (223) de compensation du filtrage simple annulation, - une étape (224) d'effacement des échos parasites du fouillis, - une étape (225) de transposition dans le domaine temporel. 4. Procédé selon la revendication précédente caractérisé en ce que la transposition dans le domaine fréquentiel est réalisée par transformée de Fourier rapide et en ce que la transposition dans le domaine temporel est réalisée par transformée de Fourier rapide inverse. 5. Procédé selon une des revendications précédentes caractérisé en ce que le test des différentes hypothèses de rang d'ambiguïté possible comporte pour chaque rang d'ambiguïté testé : - une étape (231) de regroupement des deux sous-signaux, - une étape (232) de pondération du signal regroupé, - une étape (233) de transposition dans le domaine fréquentiel, - une étape (234) de conversion en décibel, - une étape (235) de filtrage des échos issus des rangs d'ambiguïté différents du rang d'ambiguïté sous test. 6. Procédé selon la revendication précédente caractérisé en ce que l'étape de regroupement (231) pour un rang d'ambiguïté r testé correspond au regroupement des deux sous-signaux après décalage temporel (230) d'un des deux sous-signaux d'une valeur de r.kcd cases distances. 7. Procédé selon la revendication 5 ou 6 caractérisé en ce que l'étape de filtrage s'effectue en calculant la valeur minimale entre un premier signal correspondant à l'interpolation sur 2N points du premier sous-signal, un deuxième signal correspondant à l'interpolation sur 2N points du deuxième sous-signal après décalage suivant l'axe case distance d'une valeur proportionnelle au rang sous test et un troisième signal correspondant au spectre du signal regroupé du rang sous test. 8. Procédé selon la revendication précédente caractérisé en ce que l'interpolation des sous-signaux après suppression des échos parasite du fouillis sur 2N points comprend : - une étape (51) de pondération,une étape (52) d'extension par zero-padding (extension par des zéros) une étape (53) de transposition dans le domaine fréquentiel, - une étape (54) de conversion en décibel. 9. Procédé selon une des revendications 5 à 8 caractérisé en ce que la transposition dans le domaine fréquentiel est réalisée par transformée de Fourier rapide.