FR 2718250 A1 19951006 FR 9403808 A 19940331 Procédé de sondage d'un canal La présente invention concerne de manière générale un procédé de sondage d'un canal, tel que canal radioélectrique, ou canal acoustique. Un tel procédé trouve application par exemple dans les techniques de sondage de l'ionosphère. L'ionosphère étant constituée de couches ionisées, elle se comporte comme un milieu dont la constante diélectrique varie en fonction de l'ionisation et du carré de la fréquence d'une onde reçue entraînant ainsi la réflexion de cette onde reçue. L'homme du métier appréciera néanmoins que l'invention est également applicable au sondage de canaux filaires de transmission, à la détection de cibles dans des applications radar ou sonar, ou à la radiographie d'organes dans des applications médicales ou de structures dans des applications industrielles. Toutes les applications précitées sont basées sur l'émission par un émetteur d'un signal qui est réfléchi ou transmis dans le canal, et reçu en un signal de bruit additionné à un signal de sondage caractérisant le milieu dans lequel s'est propagé le signal émis en ce qui concerne par exemple la présence de cibles, de couches ionisées,...etc. Le signal de sondage reçu par réflexion ou transmission possède des caractéristiques correlées aux caractéristiques du signal émis mais avec une amplitude très atténuée et donc trés faible. La très faible amplitude du signal de sondage reçu par réflexion ou transmission n'est pas gênante en soi puisqu'il suffit d'amplifier le signal reçu avant de l'analyser. Toutefois une telle atténuation du signal de sondage reçu par réflexion ou transmission est toujours accompagnée d'un bruit de fond d'origine interne ou externe, tel que bruit atmosphérique. Cela implique que des erreurs peuvent résulter de l'analyse du signal reçu en confondant une pointe de bruit d'amplitude élevée avec le signal de sondage (fausse alarme) et inversement un signal de sondage reçu d'amplitude faible avec du bruit. Toutes les solutions de la technique antérieure remédient partiellement à cet inconvénient en effectuant par exemple des calculs de corrélation entre le signal émis et le signal reçu ou en filtrant le signal reçu. Elles nécessitent toutes des émetteurs de signal de puissance très élevée, de l'ordre de plusieurs centaines de kilowatt, pour un sondage à 1000 km, et sont donc de réalisation très coûteuse. L'invention fournit un procédé remédiant à l'inconvénient précité en ce qui concerne la puissance particulièrement élevée des émetteurs. Pour une portée prédéterminée, un émetteur nécessaire à la mise en oeuvre de l'invention offre une puissance environ de 103 à 106 fois moins élevée que celle des émetteurs de la technique antérieure. A cette fin, un procédé de sondage d'un canal par émission de signaux et traitement de signaux de bruit et sondage reçus du canal en réponse aux signaux émis, comprenant les étapes cycliques suivantes - émission d'une impulsion brève dans le canal séparée d'une impulsion brève précédente par un intervalle de temps prédéterminé - réception et acquisition d'une réponse impulsionnelle durant un intervalle de temps de durée égale à l'intervalle de temps prédéterminé succédant à l'émission d'impulsion brève, et - échantillonnage de la réponse impulsionnelle en une suite de M échantillons initiaux qui sont définis par des rangs temporels et qui sont mémorisés dans une matrice d'échantillon initial, est caractérisé en ce qu'il comprend, en outre, les étapes de: - transformation par changement d'espace de M colonnes d'échantillon constituées chacune de N échantillons initiaux qui ont des rangs temporels respectifs égaux dans N dernières réponses impulsionnelles qui ont été acquises et échantillonnées, en une matrice d'échantillon transformée, - association à la matrice d'échantillon transformée d'une matrice binaire contenant des éléments binaires qui sont respectivement associés aux échantillons transformés dans la matrice d'échantillon transformée, et qui prennent chacun l'un de deux niveaux binaires en dépendance du module dudit échantillon associé transformé par rapport à un seuil, - filtrage spatial de la matrice binaire en une matrice binaire filtrée dans laquelle chaque élément binaire prend un état représentatif de la présence ou de l'absence du signal de sondage dans l'échantillon initial associé, et - calcul d'une réponse impulsionnelle de canal exempte de bruit en fonction seulement des échantillons transformés qui sont respectivement associés à des éléments binaires de la matrice binaire filtrée ayant des niveaux qui sont représentatifs de la présence du signal de sondage. L'étape de filtrage spatial consiste à transformer chacun des éléments binaires dans la matrice binaire en un élément binaire dans la matrice binaire filtrée en fonction d'éléments binaires contigus audit chacun des éléments binaires dans la matrice binaire. De préférence, la transformation par changement d'espace est un filtrage par transformée de type Fourier et la réponse impulsionnelle calculée résulte d'une transformée de type Fourier inverse des échantillons initiaux d'une dernière réponse impulsionnelle acquise qui sont respectivement associés à des éléments binaires dans la matrice binaire filtrée ayant un état représentatif de la présence d'un signal de sondage. Selon l'invention, la transformation dite "de type Fourier signifie toute transformation connue, qu'elle soit du type Fourier, Hadamard, Rademacher, etc.. Selon une variante de l'invention, la réponse impulsionnelle de canal calculée est une réponse impulsionnelle efficace dont le module à un instant donné est égal à la valeur efficace des N échantillons initiaux qui sont définis par des rangs temporels d'échantillonnage égaux dans les N dernières réponses impulsionnelles acquises et échantillonnées, et qui sont associés à des éléments binaires dans la matrice binaire filtrée ayant un état représentatif de la présence d'un signal du sondage. Le procédé selon l'invention est de préférence adaptatif. Pour cela, il est prévu une étape de détermination succédant à l'étape de filtrage spatial, pour déterminer la valeur du seuil pour une étape cyclique suivante d'association, en dépendance d'une probabilité de fausse alarme donnée calculée relativement aux échantillons initiaux qui sont respectivement associés à des éléments binaires dans la matrice binaire filtrée ayant un état représentatif de l'absence du signal de sondage. Selon une caractéristique de l'invention, un procédé est caractérisé en ce que la valeur du seuil est déterminée à partir d'une fonction convexe qui est établie en fonction du seuil et des échantillons de la matrice filtrée associés aux éléments binaires dans la matrice binaire filtrée ayant un état représentatif de l'absence du signal de sondage. D'autres caractéristiques et avantages de la présente invention apparaîtront plus clairement à la lecture de la description suivante, en référence aux dessins annexés correspondants dans lesquels - les figures 1A, 1B, 1C et 1D sont des diagrammes illustrant respectivement des étapes d'émission d'impulsions successives, de réception de réponses impulsionnelles d'un canal sondé aux impulsions, d'échantillonnage de ces réponses impulsionnelles, et de mémorisation d'échantillons résultant de l'échantillonnage, pour la mise en oeuvre du procédé selon l'invention; - la figure 2 est un diagramme montrant la séparation, dans une matrice transformée d'échantillon, entre des échantillons relatifs à du bruit et des échantillons relatifs à du bruit additionné à un signal de sondage; - la figure 3 est un diagramme de moyenne de module en fonction de seuil de module utilisé pour séparer les échantillons relatifs à du bruit et les échantillons relatifs à du bruit additionné à un signal de sondage dans la matrice transformée; - les figures 4A à 4D montrent des exemples de motif d'échantillons contigus pour filtrer spatialement des échantillons séparés; - les figures 5A et 5B montrent respectivement une matrice d'échantillon transformée et une matrice d'échantillon séparée et filtrée; - les figures 6A et 6B sont des diagrammes montrant des réponses impulsionnelles calculées sur la base respectivement d'échantillons initiaux et d'échantillons séparés et filtrés; et - la figure 7 est un algorithme d'étapes du procédé selon l'invention. Dans la description qui suit, le terme "canal" est utilisé pour qualifier le milieu qui est sondé par onde électromagnétique, mécanique ou acoustique. La nature du "canal" n'est pas limitative : le canal peut être un canal filaire ou radioélectrique pour des applications de télécommunications ou de techniques radar ou sonar, ou bien un organe pour des applications médicales ou un produit pour des applications industrielles. Les figures 1A à 1D schématisent les étapes essentielles d'obtention d'une matrice d'échantillon ME montrée à la figure 1D sous la forme d'une mémoire qui mémorise des échantillons élémentaires résultant de l'échantillonnage de réponses impulsionnelles, ou signaux reçus au cours du temps t. La matrice d'échantillon ME est obtenue de la manière suivante. Un émetteur produit des impulsions brèves 10 à 1N-1 successives séparées deux à deux par un intervalle de temps prédéterminé. Ces impulsions sont produites à des instants respectifs lo à tN-1 et sont propagées en des ondes électromagnétiques dans le canal à sonder. En pratique, l'intervalle de temps séparant deux impulsions est directement fonction de la distance maximale de sondage. En réponse à ces impulsions brèves 10 à 1N-1 émises par l'émetteur, sont reçus par réflexion ou transmission dans le canal des signaux respectifs, appelés réponses impulsionnelles de canal RIo à RIZ~1. Chacune de ces réponses impulsionnelles est prélevée entre l'instant d'émission de l'impulsion émise correspondante et l'instant d'émission de l'impulsion émise suivante, et donc durant l'intervalle de temps prédéterminé. Chaque réponse impulsionnelle Rin, l'indice n étant compris entre O et l'entier (N-1), contient un signal de bruit mélangé éventuellement à un signal de sondage associé à l'impulsion émise correspondante Ion . Chacune des N réponses impulsionnelles RIo à RIZ~1 est ensuite échantillonnée à une fréquence prédéterminée comme montré à la figure 1C. L'échantillonnage d'une réponse impulsionnelle Rin fournit M échantillons EnO à En~1. Comme montré à la figure 1D, les M échantillons de l'une Rin des N réponses impulsionnelles RIo à RIZ~1 sont mémorisés respectivement dans M cases-mémoire contiguës de la matrice d'échantillon ME. Les N ensembles d'échantillons respectifs des N dernières réponses impulsionnelles RIo à RIZ~1 sont stockés dans N lignes de cases-mémoire de la matrice d'échantillon ME, et chacune de ces N lignes comprend M casesmémoire. Le contenu de la matrice ME s'interprète de la manière suivante : chaque ligne est associée à une réponse impulsionnelle sous forme échantillonnée, et chaque colonne correspond à une distribution d'amplitude des échantillons ayant un même rang m, avec OSmSM, respectivement dans les N réponses impulsionnelles RIo à RIN-1. Dans la description qui précède de l'obtention de la matrice d'échantillon ME, n'est considéré que le résultat figé de la mémorisation de N réponses impulsionnelles échantillonnées. En pratique, les étapes d'émission, réception, échantillonnage et mémorisation sont continues. Cela signifie qu'une (N + l)ième impulsion émise donne lieu à une (N + l)ième réponse impulsionnelle qui est échantillonnée puis écrite dans la Nième ligne de la matrice d'échantillon ME. Les échantillons mémorisés dans la Nième ligne sont déplacés dans la (N - l)ième ligne de la matrice dont les échantillons qu'elles contient sont eux-mêmes déplacés dans la (N - 2)ième ligne et ainsi de suite par récurrence jusqu'aux échantillons de la 2ième ligne qui sont déplacés dans la lère ligne. Les échantillons de la première ligne sont pour leur part effacés pour chaque mémorisation d'une nouvelle réponse impulsionnelle échantillonnée. Le contenu de la matrice d'échantillonnage est donc modifié en permanence à la fréquence d'acquisition des réponses impulsionnelles. Cela est cohérent avec le but recherché pour un sondage de canal, canal qui ne présente pas des caractéristiques constantes dans le temps. Le terme 1,impulsion" utilisé dans la description de la présente demande n'est pas limité à la notion physique restreinte qui lui est attribué dans le terme "impulsion de Dirac". Une telle impulsion peut consister par exemple en une séquence codée de "sousimpulsions" émise par l'émetteur, et utile aux techniques de compression d'impulsion. De telles techniques de compression d'impulsion visent à réduire la puissance crête de l'émetteur et augmenter le gain en puissance en étalant temporellement la durée de l'impulsion en une séquence codée de sousimpulsions présentant une énergie égale à celle de l'impulsion. Dans le récepteur sont prévus des moyens de corrélation pour corréler le signal reçu avec des répliques de la séquence déphasées temporellement, et ainsi comprimer le signal reçu en un signal proche de la réponse impulsionnelle qui aurait été reçue si une impulsion brève avait été émise. L'effet de la corrélation est de remettre en phase les raies du spectre de la séquence reçue par réflexion ou transmission dans le canal relativement à la phase de la réplique coïncidant avec la séquence reçue. Par définition, chaque réponse impulsionnelle écrite sous forme échantillonnée dans une ligne de la matrice ME ne peut pas être associée à une variable aléatoire puisqu'elle dépend des caractéristiques du canal sondé. Par contre, chaque distribution d'amplitude des composantes réelle et imaginaire des échantillons d'un même rang m mémorisée dans une colonne respective de cette matrice ME est une fonction aléatoire de valeur moyenne et de variance finie définies par les échantillons de ce rang Il est démontré ci-dessous en termes mathématiques qu'une fonction aléatoire quelconque peut être rendue gaussienne. Il est connu, en statistique, qu'une variable aléatoire X(t) qui est égale à la somme de P variables aléatoires indépendantes non gaussiennes xl(t) à xp(t), P étant un entier très grand, telle que: X(t) = x1(t) + x2(t) + + xp(t) + .... + xp(t) est une variable aléatoire qui obéit à la loi normale ou loi de Gauss. Ce résultat est connu par le théorème de la limite centrale. La définition précédente peut être expliquée de la façon suivante un évènement dû au hasard peut être considéré comme le résultat d'un grand nombre P de causes, indépendantes entr'elles, exerçant chacune une petite influence. Chacune de ces causes est purement aléatoire mais l'évènement résultant de ces causes peut être apprécié par une loi de Gauss. Si un changement d'espace par transformation, tel que transformée de Fourier discrète, est opéré sur les échantillons de la matrice ME selon l'opération où E(l,t) = Elt désigne un élément de la matrice ME, E(V ,t) est un élément dans l'espace transformé et F(l,v) est une fonction d'une famille de fonctions orthogonales définie par alors chaque échantillon E(V,t) est une variable aléatoire gaussienne. En revenant à la figure 1D, il a été vu que chaque colonne correspond à une distribution d'amplitude des échantillons de même rang respectivement pour N réponses impulsionnelles successives RIo à RIZ~1. Par nature, une distribution d'amplitude pour une colonne donnée possède un caractère aléatoire. Grâce à une transformation par changement d'espace de (7,t) à (v,t) de chacune de ces distributions d'amplitudes initiales des échantillons de même rang, correspondant à une colonne de la matrice ME (figure 1D), un caractère gaussien sera conféré aux distributions d'amplitudes des signaux issus d'une transformation des distributions d'amplitudes initiales consécutives au changement d'espace. La figure 2 montre une matrice MF de modules des échantillons transformés résultant d'un changement d'espace "en colonne" des échantillons mémorisés dans la matrice ME montrée à la figure 1D. Selon ce qui précède, chaque colonne de cette matrice MF mémorise les modules des échantillons dont les composantes réelle et imaginaire présentent la caractéristique d'une distribution de Gauss. En pratique le changement d'espace peut consister en une transformation de Fourier, une transformation de Hadamard ou une transformation de Rademacher, etc. L'axe des ordonnées, relatif aux cases-mémoire prises suivant les colonnes dans la matrice MF, est l'axe de fréquence qui se substitue à l'axe de temps (T - s) dans la matrice ME de la figure 1D. Il est rappelé que cette opération de changement d'espace a pour objectif de rendre Gaussien chacune des M distributions d'amplitude des composantes réelle et imaginaire des N échantillons respectivement dans les M colonnes. Indépendamment de ce changement d'espace, chaque échantillon d'une réponse impulsionnelle dans la mémoire à la figure 1D est porteur ou non d'un signal écho. Chaque échantillon étant porteur d'un bruit mélangé ou non avec un signal de sondage, il est considéré ci-après des relations mathématiques établissant la présence d'un signal de sondage dans un bruit. Il a été vu précédemment que le changement d'espace avait pour objectif de donner un caractère Gaussien à chacune des M distributions d'amplitude des composantes réelle et imaginaire de N échantillons, respectivement dans les M colonnes. A titre d'exemple, soit z(t) = b(t) + A.cos(a)0t) la somme d'un signal de bruit b(t) et d'un signal de sondage A.cos(Ot), supposé ici être une composante sinusoïdale pure. A des fins de simplification, b(t) peut s'écrire b(t) = x(t)cos(O0t - y(t) sinû)0t, soit z(t) = [x(t) + A] cos#Ot - y(t) sina > 0t En posant x'(t) = x(t) + A, la variable MO(t) associée au module d'un échantillon est définie telle que MO(t)2 = x'(t)2 + y(t)2, soit z(t) = MO(t) cos(#Ot + avec e = Arctg (y(t)/x'(t))=Arctg [y(t)/(x(t) + A)]. x' (t) et y(t) étant des variables indépendantes de type gaussienne, on peut écrire exp-[[(x' - A) + y]] / 2# p(x',y)dx'dy=- Z dx'dy (3), 2## où a2 désigne la variance du signal de bruit b(t), et p(x',y) est la densité de probabilité pour le couple (x',y). En passant aux coordonnées polaires, l'équation (3), avec x' = MO cose' et y = MO sin 0', s'écrit: soit où 10 désigne la fonction de Bessel modifiée de première espèce et d'ordre 0. Les résultats mathématiques indiqués ci-dessus établissent qu'en présence d'un signal de sondage reçu dans un échantillon de la matrice MF, la densité de probabilité du module MO de cet échantillon est donnée par l'équation (4). L'équation (4) signifie en outre qu'en l'absence d'un signal de sondage reçu dans un échantillon de la matrice MF, une densité de probabilité du module de cet échantillon est décrite par la loi de Raleigh donnée par la relation Le principe de détection pour sondage utilisé dans l'invention consiste à effectuer, échantillon par échantillon, dans la matrice transformée MF de la figure 2, une comparaison entre le module de cet échantillon et un seuil S qui discrimine la présence et l'absence d'un signal de sondage dans cet échantillon, présence et absence respectivement définies en termes de densité de probabilité du module de l'échantillon par les équations (4) et (5). Comme montré dans la figure 2, les échantillons dans la matrice transformée MF, qui sont obtenus par une transformation de changement d'espace d'une suite de réponses impulsionnelles RIo à RIN-1 sont séparables en des premiers échantillons contenant un signal de sondage et des seconds échantillons dont le module ne contient aucune composante de signal de sondage. Ces premiers et seconds échantillons sont groupés dans des "zones" SIGNAL+BRUIT et BRUIT dans la figure 2. La zone BRUIT, ne contenant que des échantillons de bruit, est caractérisée par une densité de probabilité définie par l'équation (5). Dans cette zone BRUIT, la probabilité PMOa qu'un échantillon ait un module MO compris entre MOa et (MOa+ AMOa) est égale à La distinction entre échantillons de la zone BRUIT et échantillons de la zone SIGNAL+BRUIT est opérée par comparaison entre le module de chaque échantillon et un seuil de module S. Si le module de l'échantillon est inférieur au seuil S, alors l'échantillon appartient à la zone BRUIT. Par contre, si le module de l'échantillon est supérieur au seuil S, l'échantillon appartient à la zone SIGNAL+BRUIT. Pour un seuil prédéterminé S, la zone BRUIT est définie par une probabilité de fausse alarme Pfa égale à la probabilité pour que le module des échantillons de la zone BRUIT soit compris entre S et +oxo. Il est donc obtenu, en utilisant l'équation (5): = exp(-S2 i 2a2). Dans cette même zone BRUIT, la moyenne m Dans la figure 3 est montré la courbe convexe C de la moyenne m Ki prenant des valeurs discrètes prédéterminées. Chacun des coefficients directeurs Ki des droites du faisceau est représentatif d'une probabilité de fausse alarme Pfa. En effet, dans l'équation de droite (m La figure 3 assiste à une détermination graphique d'un seuil optimal, noté Sop. A partir d'un groupe d'échantillons donné de la matrice transformée MF, la moyenne m KiS), cette valeur de coefficient Ki permet alors de déterminer "graphiquement" la valeur optimale du seuil S, notée Sop. En effet il suffit alors de choisir initialement un seuil minimal très faible S=SO. Pour ce seuil minimal SO est calculée une moyenne correspondante m Cette moyenne conduit à la détermination d'un seuil S=S1 selon l'équation de droite associée à la probabilité de fausse alarme Pfa fixée, laquelle droite associée correspond à Ki= K4 selon l'exemple montré dans la figure 3. Ce seuil S1 détermine alors la valeur correspondante de la moyenne m Cette résolution graphique pour la détermination d'un seuil optimal S=Sop en fonction d'une probabilité de fausse alarme donnée Pfa est optimale eu égard aux moyens de calcul à mettre en oeuvre. Le méthode graphique qui vient d'être décrite et qui fixe le seuil optimal Sop en fonction de la moyenne m Grâce à ce seuil optimal Sop calculé, la zone BRUIT et la zone SIGNAL+BRUIT dans la matrice transformée d'échantillons MF sont séparées par comparaison du module de chaque échantillon de la matrice MF avec le seuil Sop. En pratique, cette séparation de zone est matérialisée par l'association à la matrice d'échantillon transformé MF d'une matrice binaire de même dimension. Chacun des éléments binaires dans la matrice binaire est associé à l'un respectif des échantillons dans la matrice transformée MF, et prend le niveau 1 si l'échantillon associé dans la matrice MF a un module supérieur au seuil Sop, et le niveau Xo" si l'échantillon associé dans la matrice MF a un module inférieur au seuil Sop. La séparation des zones BRUIT et SIGNAL+BRUIT, matérialisée dans la matrice d'élément binaire associée, ne s'effectue qu'à un facteur de probabilité de fausse alarme près. Cela signifie qu'il existe dans la matrice binaire associée des éléments binaires à un niveau "O", qui sont donc associés à des échantillons dans la matrice MF qui sont estimés être des échantillons de bruit, alors qu'ils correspondent à des échantillons de signal de sondage additionné à du bruit, et inversement des éléments binaires à un niveau "1", qui sont associés à des échantillons qui sont estimés être des échantillons de signal, alors qu'il correspondent en fait à des échantillons de bruit. L'invention prévoit, pour remédier à ce problème, un filtrage spatial des éléments binaires de la matrice associée en une matrice binaire associée filtrée. Considérons pour cela un groupe de E éléments binaires contigus de la matrice associée. A titre d'exemple, les figures 4A à 4D montrent respectivement un groupe de E=3 x 3 éléments binaires circonscrits par un carré, un groupe de E = 3 x 4 éléments binaires circonscrits par un rectangle un groupe de E=5 éléments binaires en croix, et un groupe de E=3 éléments binaires pris en ligne ou en colonne dans la matrice binaire associée. Si le groupe de E éléments binaires est associé à E échantillons dans la matrice transformée MF appartenant à la zone BRUIT, la probabilité PB(F) pour qu'un nombre F d'éléments binaires, avec F PB(F) = CE-PfaFo Pfa) (7) F E! où C E = désigne le nombre de combinaisons F!(E - F)! de F éléments binaires parmi les E éléments binaires du groupe. De même, si ce groupe de E éléments binaires est associé à un groupe de E échantillons appartenant à la zone SIGNAL + BRUIT, la probabilité PS(F) pour qu'un nombre F d'échantillons, avec F où Pd dénote une probabilité de détection,Cp désigne le nombre de combinaisons de i échantillons parmi p échantillons, et CEp le nombre de combinaisons de (F-i) échantillons parmi (E-p) échantillons. L'équation (8) résulte du fait que dans la zone SIGNAL+BRUIT, un échantillon donné peut n'être porteur que de bruit ou porteur de signal et de bruit. En désignant par (E-p) le nombre d'échantillons du groupe porteurs uniquement de bruit, la probabilité que F échantillons parmi E aient un module supérieur à S résulte de la probabilité (CpPdi(l - Pd)P i) que i é central est considéré porteur de signal et bruit si et seulement si ledit élément binaire central a un niveau égal à "O" et si, dans le groupe d'éléments binaires, au moins F éléments ont un niveau égal à "1". Cette étape de filtrage spatial est cohérente avec la probabilité très faible prise par PB(F) dans la zone BRUIT; et (b) pour chaque groupe d'éléments binaires à E éléments binaires centré sur l'un central respectif des éléments binaires de la matrice associée, l'échantillon associé audit élément binaire central est considéré porteur de signal et bruit si et seulement si ledit élément binaire central a un niveau égal à "1" et si, dans ledit groupe, au moins (F-l) éléments binaires ont un niveau égal à "1". Cette étape de filtrage spatial est cohérente avec la probabilité relativement élevée prise par PS(F) dans la zone SIGNAL + BRUIT. A l'issue de ces étapes de filtrage (a) et (b) est déduite une matrice binaire associée filtrée spatialement qui résulte de la matrice associée dans laquelle les niveaux des éléments binaires prennent la valeur "1" si les échantillons respectifs associés sont estimés être porteurs de signal selon ces étapes de filtrage spatial, et la valeur "O" si les échantillons respectifs associés sont estimés n'être porteurs que de bruit toujours selon les étapes (a) et (b). Les figures 5A et 5B montrent respectivement une matrice d'échantillon transformée MF résultant d'une transformation de changement d'espace de la matrice ME, et en correspondance avec cette dernière matrice, une matrice d'échantillon traitée MT résultant de la "suppression" de tous les échantillons dans la matrice MF associés respectivement à des éléments binaires de niveau "O" dans la matrice binaire associée filtrée. Dans cette dernière matrice MT, seuls sont donc conservés les échantillons appartenant à la zone SIGNAL+BRUIT, les modules des échantillons de BRUIT ayant été mis à zéro. De cette matrice traitée MT est déduite la réponse impulsionnelle du canal sondé, représentative de cibles, couches ionisées, etc.. en dépendance de la nature de ce canal. Plus précisément sont déduites une réponse impulsionnelle efficace et/ou une réponse impulsionnelle instantanée. En référence à la figure 5B, l'amplitude d'une réponse impulsionnelle efficace à un instant t, associée à un rang m d'échantillonnage, avec m compris entre 0 et (ME1), est égale à la valeur efficace des modules de N échantillons initiaux dans la matrice traitée MT entre des bornes fréquentielles vmin et vmax pour une transformation de type transformée de Fourier, lesdits N échantillons ayant des rangs d'échantillonnage égaux, c'est-à-dire coïncidant avec un décalage temporel égal par rapport aux instants d'émission des impulsions respectives. Si H(t) désigne la réponse impulsionnelle efficace, sa valeur à un instant discret ti s'écrit MO(ti,v) désignant le module d'un échantillon de coordonnées (ti,v) dans la matrice traitée MT. En référence à la figure 1D, la réponse impulsionnelle instantanée est obtenue de la manière suivante. Pour chaque dernière réponse impulsionnelle RIZ~1 acquise, celle-ci est échantillonnée puis mémorisée en les échantillons E Ces échantillons sont traités dans les étapes de transformation par changement d'espace et filtrage spatial selon l'invention par modification de la matrice d'échantillon ME en la matrice d'échantillon traitée MT, et par utilisation de matrices binaires intermédiaires associée et associée filtrée. Selon l'exemple de l'invention, la transformation par changement d'espace consiste en la transformée de Fourier "en colonnes" des échantillons. La réponse impulsionnelle instantanée résulte alors de la transformée de Fourier inverse des échantillons E(N1)0 à E(N-1)M-1 dans la matrice traitée MT après leur traitement par changement d'espace et filtrage spatial. Ces échantillons filtrés sont alors expurgés de la quasi-totalité du bruit qui ne dégrade donc pas la réponse impulsionnelle instantanée. A titre d'exemple, les figures 6A et 6B montrent deux réponses impulsionnelles efficaces obtenues respectivement à partir d'échantillons initiaux dans la matrice d'échantillon MF, et à partir d'échantillons correspondants dans le matrice traitée MT résultant d'une transformation par changement d'espace et d'un filtrage spatial de la matrice ME selon l'invention. On notera que trois cibles sont détectées à des instants to, t1 et t2. Sans les traitements selon l'invention, la réponse impulsionnelle (figure 6A) est mélangée à du bruit qui induit une confusion dans la détection de la cible située en t2. La figure 7 résume les étapes essentielles ET1 à ET6 de la mise en oeuvre du procédé selon l'invention. L'étape ET1 est relative à l'émission des impulsions successives IO, 11, IN-lr la réception et l'acquisition des réponses impulsionnelles respectives RIo, RI1,... RIN-1 renvoyées par le canal sondé en réponse à ces impulsions, à l'échantillonnage et à la mémorisation de ces réponses impulsionnelles comme montrées dans les figures 1A à 1D. A la seconde étape ET2, les réponses impulsionnelles échantillonnées, transcrites dans la matrice d'échantillon ME, subissent une transformation de changement d'espace. Cette transformation est appliquée aux échantillons initiaux en "colonnes" respectifs des réponses impulsionnelles RIo ---RIN-l résultant de l'échantillonnage à des instants présentant un même décalage temporel relativement aux impulsions émises correspondantes, c'est-à-dire aux M colonnes à N échantillons de même rang. Cette transformation, comme on l'a vu, présente la propriété de donner un caractère gaussien à chaque distribution des composantes réelle et imaginaire des échantillons ayant un même rang. A cette matrice transformée MF ainsi obtenue par transformation par changement d'espace est associée une matrice binaire dans laquelle les éléments binaires sont respectivement associés aux échantillons de la matrice MF, et prennent un niveau "1" pour des échantillons respectifs associés de module supérieur à un seuil S, et un niveau "O" pour des échantillons respectifs associés de module inférieur à ce seuil S (étape ET3). Afin de parfaire cette distinction entre échantillons porteurs uniquement de bruit et échantillons porteurs de bruit mélangé au signal de sondage, l'étape suivante consiste en le filtrage spatial des éléments binaires de la matrice associée, précédemment commenté en référence aux figures 4A à 4D. A l'issue de ce filtrage spatial est fournie une matrice binaire associée filtrée qui redéfinit la caractéristique de chaque échantillon dans la matrice d'échantillon ME en ce qu'il est porteur ou pas de signal. Il peut donc être défini un nouveau seuil en fonction de la moyenne m Comme montré par la flèche reliant les deux blocs relatif à l'étape ET5 au bloc relatif à l'étape ET3, le procédé selon l'invention est adaptatif. Ce caractère adaptatif est inhérent au seuil S utilisé pour définir une matrice binaire associée relative à une première séparation des échantillons en ce qu'ils sont porteurs ou pas de bruit. Le seuil S est déduit de la matrice binaire associée filtrée, ou plus précisément de la moyenne m Selon dans le préambule de la description, l'un des objectifs principaux de l'invention consiste en une réduction considérable de la puissance nécessaire à un sondage de canal. L'invention procure alors l'avantage de réaliser plusieurs sondages simultanés à différentes fréquences d'émission, sans recourir à une puissance très élevée de l'émetteur compte tenu de la faible puissance nécessaire à chaque sondage selon l'invention. Cette simultanéité de sondages n'est pas concevable selon la technique antérieure au motif que la puissance nécessaire à chaque sondage est élevée. La mise en oeuvre du procédé dans de tels systèmes de sondage à fréquence variable, dits de sondage "flash", autorise une réduction de la durée de sondage d'un facteur Q égal au nombre Q de sondages simultanés à différentes fréquences respectives, ledit nombre Q étant limité par la puissance disponible de l'émetteur. REVENDICATIONS 1 - Procédé de sondage d'un canal par émission de signaux et traitement de signaux de bruit et sondage reçus dudit canal en réponse auxdits signaux émis, comprenant les étapes cycliques suivantes - émission d'une impulsion brève (I1....., IN-1) dans ledit canal séparée d'une impulsion brève précédente (Io,...IN~2) par un intervalle de temps prédéterminé (#1 - #0,. ., #n-1 - #n-2), - réception et acquisition d'une réponse impulsionnelle (RI1,...RIN-1) durant un intervalle de temps de durée égale audit intervalle de temps prédéterminé succédant à ladite émission d'impulsion brève (Io...IN~1) et - échantillonnage de ladite réponse impulsionnelle (RIo,...RIN~l) en une suite de M échantillons initiaux (EOg à EOM-1,...E(N-1)O à E(N1)M-1) qui sont définis par des rangs temporels (O à (M-i)), et qui sont mémorisés dans une matrice d'échantillon initial (ME), caractérisé en ce qu'il comprend, en outre, les étapes de: - transformation par changement d'espace (ET2) de M colonnes d'échantillon constituées chacune de N échantillons initiaux (EOg à E(N-1)O,...EOM-1 à E(N1)M-1) qui ont des rangs temporels respectifs égaux dans N dernières réponses impulsionnelles qui ont été acquises et échantillonnées, en une matrice d'échantillon transformée (MF), - association (ET3) à ladite matrice d'échantillon transformée (MF) d'une matrice binaire contenant des éléments binaires qui sont respectivement associés aux échantillons transformés dans ladite matrice d'échantillon transformée (MF), et qui prennent chacun l'un de deux niveaux binaires ("1", "0") en dépendance du module dudit échantillon associé transformé par rapport à un seuil (S), - filtrage spatial (ET4) de ladite matrice binaire en une matrice binaire filtrée dans laquelle chaque élément binaire prend un état représentatif de la présence ou de l'absence du signal de sondage dans l'échantillon initial associé, et - calcul (ET6) d'une réponse impulsionnelle de canal sensiblement exempte de bruit en fonction seulement des échantillons transformés qui sont respectivement associés à des éléments binaires de la matrice binaire filtrée ayant des niveaux qui sont représentatifs de la présence du signal de sondage. 2 - Procédé conforme à la revendication 1, caractérisé en ce que ladite étape de filtrage spatial (ET4) consiste à transformer chacun des éléments binaires dans ladite matrice binaire en un élément binaire dans ladite matrice binaire filtrée en fonction d'éléments binaires contigus audit chacun des éléments binaires dans la matrice binaire. 3 - Procédé conforme à la revendication 1 ou 2, caractérisé en ce que ladite transformation par changement d'espace (ET2) est obtenue par transformée de type Fourier et en ce que ladite réponse impulsionnelle calculée résulte d'une transformée de type Fourier inverse des échantillons initiaux d'une dernière réponse impulsionnelle acquise qui sont respectivement associés à des éléments binaires dans la matrice binaire filtrée ayant un état représentatif de la présence du signal de sondage. 4 - Procédé conforme à l'une quelconque des revendications 1 à 3, caractérisé en ce ladite réponse impulsionnelle de canal calculée est une réponse impulsionnelle efficace dont le module à un instant donné (ti) est égal à la valeur efficace des modules (MO) de N échantillons initiaux qui sont définis par des rangs temporels d'échantillonnage (O à M) égaux dans lesdites N dernières réponses impulsionnelles acquises et échantillonnées, et qui sont associés à des éléments binaires dans la matrice binaire filtrée ayant un état représentatif de la présence du signal de sondage. 5 - Procédé conforme à l'une quelconque des revendications 1 à 4, caractérisé par une étape de détermination (ET5) succédant à ladite étape de filtrage spatial (ET4), pour déterminer la valeur dudit seuil (S) pour une étape cyclique suivante d'association, en dépendance d'une probabilité de fausse alarme donnée (Pfa) calculée relativement aux échantillons initiaux qui sont respectivement associés à des éléments binaires dans la matrice binaire filtrée ayant un état représentatif de l'absence du signal de sondage. 6 - Procédé conforme à la revendication 5, caractérisé en ce que la valeur du seuil (S) est déterminée à partir d'une fonction convexe qui est établie en fonction dudit seuil et desdits échantillons de la matrice filtrée (MF) associés auxdits éléments binaires dans la matrice binaire filtrée ayant un état représentatif de l'absence du signal de sondage. 7 - Procédé de sondage d'un canal par émission de signaux et traitement de signaux de bruit et sondage reçus dudit canal en réponse auxdits signaux émis, caractérisé par une pluralité de sondages conformes à l'une quelconque des revendications précédentes, respectivement pour une pluralité de fréquences d'émission d'impulsions brèves.