Brevet FR2932339A1 — publié le 20091211

Procédé et dispositif de sondage par propagation d'ondes. La présente invention est relative aux procédés et dispositifs de sondage par propagation d'ondes. Plus particulièrement, l'invention concerne un procédé de sondage par propagation d'ondes, comprenant : (a) une étape d'émission au cours de laquelle on fait émettre par un ensemble de transducteurs (c'est-à-dire par tout ou partie des transducteurs de l'ensemble) une onde incidente dans un milieu diffusant pour ladite onde, (b) une étape de mesure au cours de laquelle on fait capter, par ledit ensemble de transducteurs (c'est-à-dire par tout ou partie des transducteurs de l'ensemble), des signaux représentatifs d'une onde réfléchie réverbérée par le milieu à partir de l'onde incidente, lesdits signaux captés comprenant : - une composante de diffusion simple, représentative de trajets d'onde où l'onde réfléchie résulte d'une seule réflexion de l'onde incidente par des diffuseurs appartenant au milieu, et le cas échéant une composante de diffusion multiple, représentative de trajets d'onde où l'onde réfléchie résulte de plusieurs réflexions successives de l'onde incidente sur les diffuseurs du milieu avant d'atteindre l'ensemble de transducteurs, (c) une étape de traitement au cours de laquelle on traite lesdits signaux captés en vue de déterminer des caractéristiques du milieu (les caractéristiques en question peuvent consister en une image du milieu, et/ou une valeur d'un paramètre du milieu, et/ou la présence ou non d'un point singulier tel qu'une hétérogénéité, etc.). L'étape de traitement susmentionnée permet par exemple de mesurer un paramètre caractéristique du milieu, et/ou de détecter un point singulier du milieu, et/ou de réaliser une image du milieu. 2 Les procédés de ce type sont utilisés notamment dans les systèmes de détection et d'imagerie, comme par exemple les sonars, les radars, les échographes etc. Dans les procédés connus de ce type et notamment dans les procédés d'imagerie échographique ou radar, on utilise la composante de diffusion simple des signaux captés : si chaque diffuseur n'interagit qu'une seule fois avec l'onde, il y a en effet une équivalence directe entre le temps d'arrivée de chaque écho et la distance séparant le transducteur et le diffuseur ayant engendré cet écho. La détection d'un écho à un instant donné est le signe de la présence d'un diffuseur à la distance correspondant au temps d'arrivée de l'écho. Une image de la réflectivité du milieu, c'est-à-dire de la position des différents diffuseurs au sein du milieu, peut alors être le cas échéant construite à partir des signaux captés. En revanche, la diffusion multiple n'est pas utilisée dans les méthodes d'imagerie échographiques ou radar. Bien au contraire, ces méthodes d'imagerie reposent sur l'hypothèse que cette diffusion multiple est négligeable. Or, en présence d'une composante importante de diffusion multiple, notamment lorsque les diffuseurs contenus dans le milieu ont un pouvoir diffusant élevé et/ou sont très denses à l'intérieur du milieu, les méthodes d'imagerie classiques sont fortement perturbées et ne sont plus fiables. En effet, dans ce cas, il n'y a plus d'équivalence entre le temps d'arrivée d'un d'écho et distance séparant un transducteur et un diffuseur du milieu, ce qui ne permet pas de bâtir une image du milieu. La présente invention a notamment pour but de perfectionner les procédés de sondage tels que définis ci-dessus, de façon à prendre en compte la composante de diffusion multiple. A cet effet, selon l'invention, au cours de l'étape 35 de traitement, on extrait, par filtrage, au moins une composante choisie parmi la composante de diffusion multiple et la composante de diffusion simple. On peut mettre à profit cette séparation, par exemple : - pour utiliser uniquement la composante de diffusion simple en vue de bâtir une image du milieu, - ou encore pour déterminer un paramètre caractéristique du milieu à partir de la composante de diffusion multiple, etc. On notera que l'invention est applicable lorsque les signaux captés ne comportent pas effectivement de composante de diffusion multiple . l'invention permet en effet alors de déterminer que la diffusion multiple est nulle, ce qui peut fournir une indication utile pour caractériser le milieu sondé. De plus, même en l'absence de contribution de diffusion multiple significative, l'invention permet d'éliminer en bonne partie le bruit et les effets d'aberration dans les signaux captés lorsqu'on extrait la contribution de diffusion simple, ce qui peut s'avérer extrêmement utile. Dans divers modes de réalisation du procédé selon l'invention, on peut éventuellement avoir recours en outre à l'une et/ou à l'autre des dispositions suivantes : - au cours de l'étape de traitement (c), on extrait la composante de diffusion simple et on utilise ladite composante de diffusion simple pour détecter (et généralement localiser) au moins un point singulier dans le milieu ; au cours de l'étape de traitement (c), on 30 réalise une image du milieu à partir de ladite composante de diffusion simple ; - au cours de l'étape de traitement (c), on détermine ladite composante de diffusion simple sous forme d'une matrice temporelle filtrée de réponse inter-éléments 35 HF, représentative des réponses entre transducteurs de 4 l'ensemble de transducteurs, et on réalise l'image du milieu par formation de voie à partir de ladite matrice temporelle filtrée de réponse inter-éléments HF ; - au cours de l'étape de traitement (c), on détermine ladite composante de diffusion simple sous forme d'au moins une matrice de transfert fréquentielle filtrée KF, représentative des réponses entre transducteurs de l'ensemble de transducteurs, et on réalise l'image du milieu par la méthode de décomposition de l'opérateur de retournement temporel KF*.KF, où KF*est la matrice conjuguée complexe de ladite matrice de transfert fréquentielle filtrée KF ; au cours de l'étape de traitement (c), on extrait la composante de diffusion multiple et on calcule à partir de ladite composante de diffusion multiple, un indice de fiabilité représentatif de l'importance de la diffusion simple dans le milieu ; - au cours de l'étape de traitement (c), on réalise une image du milieu et on quantifie la fiabilité de ladite image en fonction dudit indice de fiabilité ; - au cours de l'étape de traitement (c), on calcule en plusieurs points ledit indice de fiabilité et on quantifie la fiabilité de plusieurs parties de l'image correspondant auxdits plusieurs points, en fonction dudit indice de fiabilité ; - au cours de l'étape de traitement (c), on sépare la composante de diffusion multiple de la composante de diffusion simple par filtrage d'au moins une matrice de transfert fréquentielle représentative des réponses entre transducteurs de l'ensemble de transducteurs ; - l'étape de traitement (c) comporte au moins les sous-étapes suivantes : (cl) une sous-étape de détermination de matrice de transfert fenêtrée au cours de laquelle on détermine au 35 moins une matrice de transfert fréquentielle fenêtrée K(T,f)=[ki](T,f)] correspondant à une matrice temporelle fenêtrée de réponse inter-éléments K(T,t)=[ki;(T,t)], ladite matrice temporelle fenêtrée de réponse inter-éléments correspondant, sur une fenêtre temporelle voisine d'un temps T et de durée At, aux réponses temporelles hi;(t) entre transducteurs de l'ensemble de transducteurs, f étant la fréquence, (c3) une sous-étape de filtrage au cours de laquelle on sépare la composante de diffusion multiple de la composante de diffusion simple dans chaque matrice de transfert fréquentielle fenêtrée K(T,f) en fonction de la cohérence des coefficients ki;(T,f) de la matrice de transfert fréquentielle fenêtrée K(T,f) sur chaque antidiagonale de ladite matrice de transfert fréquentielle fenêtrée K(T,f) ; - au cours de la sous-étape de filtrage (c3), on extrait une composante de diffusion simple ksij(T,f) de chaque coefficient ki](T,f) de la matrice K(T,f), telle que les composantes de diffusion simple ksij(T,f)soient cohérentes entre elles sur chaque antidiagonale de ladite matrice K(T,f) ; - au cours de la sous-étape de filtrage (c3), on extrait une composante de diffusion multiple kMii(T,f) de chaque coefficient kii(T,f) de la matrice K(T,f) (cette composante est telle que lesdites composantes de diffusion multiples kMi;(T,f) soient aléatoires) ; - l'ensemble de transducteurs comporte au moins un nombre N de transducteurs et la matrice de transfert fréquentielle fenêtrée K(T,f) est une matrice de dimension N*N ; - au cours de la sous-étape (cl) de détermination de matrice de transfert fenêtrée, on détermine plusieurs matrices de transfert fenêtrées K(T,f) pour plusieurs valeurs de T et on procède à la sous-étape de filtrage (c3) pour chacune des matrices de transfert fenêtrées K(T,f) ; 6 - les fenêtres temporelles des différentes matrices de transfert fenêtrées K(T,f) sont jointives (lesdites fenêtres temporelles peuvent soit être adjacentes, soit se recouvrir mutuellement), - au cours de la sous-étape (cl) de détermination de matrice de transfert fenêtrée, on détermine chaque matrice fréquentielle fenêtrées K(T,f) par transformée de Fourier d'une matrice temporelle fenêtrée K(T,t) correspondant, sur ladite fenêtre temporelle voisine du temps T et de durée Lit, aux réponses temporelles h;_3(t) entre transducteurs de l'ensemble de transducteurs ; - dans le procédé : l'étape de traitement (c) comporte en outre une sous-étape de rotation de données (c2), qui est intermédiaire entre la sous-étape (cl) de détermination de matrice de transfert fenêtrée et la sous-étape (c3) de filtrage, et au cours de laquelle on calcule deux matrices Al(T,f)=[aluv(T,f)] et A2=[a2,v(T,f)] à partir de K(T,f), où : aluv(T, f) =ku+v-1, v-u+2M-1 (T, f) , a2uv (T, f) =ku+v, v-u+2M-1 (T, f) , M=(N+3)/4, N est une plus petite dimension de la matrice de transfert K(T,f) et est choisi pour que M soit entier, au cours de la sous-étape de filtrage (c3), on sépare la composante de diffusion multiple de la composante de diffusion simple dans chacune des matrices Ar, r étant un indice égal à 1 ou 2, et on obtient ainsi au moins deux matrices filtrées ArF correspondant respectivement aux matrices Ar et représentatives chacune, soit de la composante de diffusion simple, soit de la composante de diffusion multiple, l'étape de traitement (c) comporte en outre une sous-étape de rotation de données inverse (c4), qui est postérieure à la sous-étape de filtrage (c3), et au cours de laquelle on calcule une matrice de transfert fenêtrée filtrée KF (T, f) = [kFij(T, f) ] , où : - lorsque i-j est pair : kFi7(T,f)=a1F(i_jti2+M, (1+3)/21 - lorsque i-j est impair : kFi,(T, f) =a1F(i_j_i),2+M, ; 5 - au cours de la sous-étape de filtrage (c3), on calcule les deux matrices ArF par la formule : ArF =S`S`Ar , où : S=[su] est un vecteur colonne, tS* étant la transposée du vecteur conjugué du vecteur S, 10 - les composantes su du vecteur S sont des nombres k 2 complexes égaux à : su = s.exp j yu 2R s est une constante (éventuellement égale à 1), yu =(x;ùxAU , avec u=(i-j)/2+M si r=1 et u=(i-j-1)/2+M si r=2, 15 - xi et x3 sont les abscisses des transducteurs d'indices i et j selon un axe X, l'ensemble de transducteurs s'étendant au moins selon ledit axe X, L=2M-1 pour r=1 et L=2M-2 pour r=2, lesdites matrices ArF étant représentatives de la 20 composante de diffusion simple ; - au cours de la sous-étape de filtrage (c3) : on décompose en valeurs singulières chacune des matrices Ar en déterminant ainsi, pour chacune des matrices Ar, M valeurs singulières i étant un indice compris 25 entre 1 et M, les valeurs singulières étant ordonnées par ordre décroissant pour i croissant, - pour chaque matrice Ar, on attribue à la diffusion simple p valeurs singulières les plus importantes 2'1,r à Àp,r, p étant compris entre 1 et M, et on attribue à la 30 diffusion multiple les éventuelles autres valeurs singulières non nulles ; - la détermination des p valeurs singulières attribuées à la diffusion simple est effectuée de manière itérative, comme suit : (c31) on fixe initialement q = 1 ; (c32) on calcule une valeur singulière normalisée Àq,r à partir de Âq r : 2q,r M 1 V22 M+1ùqm~=4 ~m,r (c33) si2q,rest au moins égal à une valeur seuil prédéterminée 2max, on attribue 2q,r à la diffusion simple et on réitère la sous-étape (c32) en incrémentant q d'une 10 unité, (c34) si Àq,r est inférieur à la valeur seuil on attribue 2q,r et les éventuelles valeurs singulières suivantes à la diffusion multiple ; - lesdites ondes sont choisies parmi les ondes 15 mécaniques et les ondes électromagnétiques. Par ailleurs, l'invention a également pour objet un dispositif pour la mise en oeuvre d'un procédé de sondage tel que défini ci-dessus, comprenant un ensemble de transducteurs adaptés pour émettre une onde incidente dans 20 un milieu diffusant et pour capter des signaux représentatifs d'une onde réfléchie réverbérée par le milieu à partir de l'onde incidente, lesdits signaux captés comprenant : une composante de diffusion simple, 25 représentative de trajets d'onde où l'onde réfléchie résulte d'une seule réflexion de l'onde incidente par chaque diffuseur du milieu, - et le cas échéant une composante de diffusion multiple, représentative de trajets d'onde où l'onde 30 réfléchie résulte de plusieurs réflexions successives de l'onde incidente sur les diffuseurs du milieu avant d'atteindre l'ensemble de transducteurs, ledit dispositif comportant en outre des moyens de traitement adaptés pour traiter lesdits signaux captés en vue de déterminer des caractéristiques du milieu, caractérisé en ce que les moyens de traitement sont adaptés pour extraire, par filtrage, au moins une composante choisie parmi la composante de diffusion multiple et la composante de diffusion simple. D'autres caractéristiques et avantages de l'invention apparaîtront au cours de la description suivante de plusieurs de ses modes de réalisation, donnés à titre d'exemples non limitatifs, en regard des dessins joints. Sur les dessins : - la figure 1 est une vue schématique d'un dispositif pour la mise en oeuvre d'un procédé selon un mode de réalisation de l'invention, - la figure 2 est un schéma représentant deux chemins de diffusion simple contribuant au champ rétrodiffusé par le milieu sondé au temps T=2R/c dans l'exemple de la figure 1, - la figure 3 est un schéma représentant deux chemins de diffusion multiple contribuant au champ rétrodiffusé par le milieu sondé au temps T=2R/c dans l'exemple de la figure 1, - les figures 4 et 5 illustrent une sous-étape de rotation de données appartenant à un mode de réalisation du procédé selon l'invention, au cours de laquelle on extrait deux matrices Al et A2 de chaque matrice fréquentielle de transfert fenêtrée K, - la figure 6 illustre une sous-étape de rotation de données inverse au cours de laquelle on détermine une matrice de transfert filtrée KF à partir des matrices Al et 35 A2 susmentionnées, après filtrage desdites matrices, 10 - la figure 7 illustre un processus itératif utilisable pour le filtrage des matrices Al et A2 par décomposition en valeurs singulières, - la figure 8 illustre un dispositif expérimental 5 utilisé une expérience de détection de cible en milieu diffusant, - les figures 9 et 10 sont des exemples d'image du milieu de la figure 8, obtenues par échographie classique (avec une échelle de gris linéaire pour la figure 9 et une 10 échelle logarithmique pour la figure 10), - les figures 11 et 12 sont des graphes distance-fréquence où l'échelle de gris montre la première valeur singulière normalisée en fonction de la fréquence f et de la distance R par rapport à l'ensemble de transducteurs 15 dans le dispositif de la figure 8, après filtrage par décomposition en valeurs singulières (sans effet de seuil sur la figure 11 et en se limitant aux valeurs singulières supérieures à un seuil prédéterminé sur la figure 12), - la figure 13 est une image du milieu 20 reconstituée par repropagation numérique des vecteurs propres, après décomposition en valeurs singulières de l'opérateur de retournement temporel, - la figure 14 est un graphe donnant l'amplitude normalisée A des pixels de l'image du milieu du dispositif 25 de la figure 8, dans un plan perpendiculaire à l'axe Z et passant par la cible, en échographie classique et après un filtrage par décomposition en valeurs singulières selon un mode de réalisation de l'invention, - la figure 15 est un graphe donnant la courbe 30 mesurée de l'intensité simplement diffusée Is(X=O, T) (représentée par des points en forme de petits cercles) et la courbe théorique de cette intensité pour un libre parcours d'extinction 1ext=50 mm, - la figure 16 est un graphe similaire (en échelle 35 logarithmique) donnant la courbe mesurée de l'intensité 11 multiplement diffusée IM(X=O, T) (représentée par des points en forme de petits cercles) et la courbe théorique de cette intensité pour quatre valeurs de libre parcours moyen élastique 1e, ces courbes étant normalisées par le maximum de IS, - et la figure 17 est une image échographique de tissus humains, couplée à un graphe dont le niveau de gris donne la valeur d'un indice de fiabilité fondé sur l'importance relative de la diffusion multiple. La figure 1 montre un exemple de dispositif adapté pour sonder un milieu 1 par émission et réception d'ondes. L'invention sera décrite plus particulièrement ci-après dans le cas où les ondes en question sont des ondes ultrasonores (par exemple de fréquences comprises entre 2 et 4 MHz), mais l'invention est également applicable à des ondes de toute nature, par exemple des ondes mécaniques autres que des ondes ultrasonores de compression, des ondes électromagnétiques, ou autres. Le milieu 1 est diffusant pour les ondes considérées, c'est-à-dire qu'il est hétérogène et/ou contient des diffuseurs 2 répartis de manière aléatoire et capables de réfléchir les ondes émises dans le milieu 1. Le milieu 1 en question peut être par exemple une partie du corps humain, et les diffuseurs 2 peuvent être notamment des particules de petite taille, non résolues, contenues dans le milieu 1 (en échographie de tels diffuseurs engendrent des images dites de �speckle�). Bien entendu, le milieu 1 à sonder pourrait être autre, par exemple une partie d'un objet industriel dont on souhaite contrôler la structure dans le cadre d'un contrôle non destructif. Le dispositif de sondage représenté sur la figure 1 comporte un ensemble 3 de transducteurs 4, par exemple un réseau linéaire de transducteurs ultrasonores pouvant se présenter classiquement sous la forme d'une barrette rigide 12 mise en contact avec le milieu 1. Cet ensemble 3 comporte un nombre N de transducteurs, N pouvant par exemple être compris entre 100 et 500. On peut par exemple utiliser une barrette d'une centaine de transducteurs 4 ayant chacun une largeur d'environ de largeur 0.39 mm. On notera que les transducteurs 4 sont ici des transducteurs ultrasonores capables de transformer un signal électrique en ondes ultrasonores de compression ou inversement, mais au sens du présent brevet, on appellera transducteur plus généralement tout dispositif capable de transformer un signal électrique en une onde de quelque nature que ce soit (onde mécanique, onde électromagnétique, etc.) et inversement. Chaque transducteur 4 de l'ensemble 3 peut être commandé individuellement par une unité centrale 5 (UC) comprenant par exemple des moyens des moyens numériques de traitement de signaux, cette unité centrale 5 pouvant par exemple être adaptée pour présenter une image du milieu 1 sur un écran 6. Pour sonder le milieu 1, l'unité centrale 5 envoie aux transducteurs 4 des signaux électriques qui sont transformés par lesdits transducteurs en ondes émises dans le milieu 1, en l'occurrence des ondes ultrasonores de compression, et ces ondes sont partiellement réfléchies par les diffuseurs 2 contenus dans le milieu. Une partie des ondes diffusées (ou échos) revient ainsi vers les transducteurs 4 qui les captent et les transforment en signaux électriques de réception traités ensuite par l'unité centrale 5. Ces ondes reviennent vers les transducteurs 4 : - soit après avoir été réfléchies une seule fois par les diffuseurs 2 (par exemple selon le trajet schématisé en traits pointillés sur la figure 1), c'est à dire après avoir subi une diffusion simple, soit après avoir été réfléchies successivement 35 sur plusieurs diffuseurs 2 (par exemple selon le trajet 13 schématisé en trait plein sur la figure 1), c'est-à-dire après avoir subi une diffusion multiple. L'onde totale diffusée par le milieu 1 et renvoyée vers les transducteurs 4, comprend donc deux contributions . - une contribution dite �de diffusion simple�, correspondant aux cas où l'onde incidente ne subit qu'une seule réflexion par chaque diffuseur 2 du milieu 1 avant de revenir vers les transducteurs 4, - une contribution dite �de diffusion multiple�, correspondant aux cas où l'onde incidente subit plusieurs réflexions successives sur les diffuseurs 2 du milieu avant d'atteindre les transducteurs 4. La présente invention permet de séparer ces deux contributions par filtrage, afin d'utiliser seulement l'une d'entre elle ou de les traiter séparément. Par exemple : - la contribution de diffusion simple est utile pour construire une image du milieu, notamment par échographie dans le cas décrit ici (mais il pourrait s'agir aussi d'une image sismographique dans le cas d'ondes sismiques, ou une image radar dans le cas d'ondes électromagnétiques, ou autres) : en effet, puisque chaque diffuseur 2 n'interagit qu'une seule fois avec l'onde incidente, il y a une équivalence directe entre le temps d'arrivée de chaque écho et la distance séparant le capteur et le diffuseur ayant engendré cet écho. La détection d'un signal écho à un instant donné est le signe de la présence d'un diffuseur 2 à la distance correspondant au temps d'arrivée de l'écho. Une image de la réflectivité du milieu, c'est-à-dire de la position des différentes cibles diffusantes au sein du milieu, peut alors être construite à partir des signaux enregistrés ; - la contribution de diffusion multiple peut également être utile à des fins de caractérisation des milieux désordonnés, car elle permet de mesurer des 14 paramètres statistiques comme le libre parcours moyen ou le coefficient de diffusion d'une onde se propageant à travers un milieu aléatoire ; l'importance relative de cette contribution de diffusion multiple permet également de caractériser la fiabilité d'une image du milieu 1, obtenue par exemple par échographie. Afin de séparer les contributions de diffusion simple et de diffusion multiple, on enregistre tout d'abord les réponses inter-éléments de chaque paire de transducteurs 4 de l'ensemble 3. A cet effet, comme représenté sur la figure 2, une onde incidente correspondant à un signal impulsionnel est envoyée depuis chaque transducteur n°i de l'ensemble 3 (d'abscisse xi) au cours d'une étape d'émission (a) et l'onde rétrodiffusée (réverbérée) par le milieu suite à l'émission de cette onde incidente, est enregistrée par chaque transducteur n°j (d'abscisse xj) au cours d'une étape de mesure (b). Au cours de cette opération, chaque transducteur j capte un signal temporel hi3 (t) suite à une émission par le transducteur i. Cette opération est réalisée pour tous les couples émetteurs/récepteurs possibles de l'ensemble 3, c'est-à-dire en pratique pour tous les couples de transducteurs 4. L'ensemble des N2 réponses forme une matrice temporelle de réponse inter-éléments H(t)=[hlj(t)] matrice carrée de dimension N*N, qui est la réponse globale du milieu 1. On notera que la matrice temporelle de réponse inter-éléments H(t) peut éventuellement être acquise plus rapidement, sans faire émettre un signal impulsionnel successivement par chaque transducteur i de l'ensemble 3, en procédant par exemple comme enseigné dans le document WO-A-2004/086557. Les signaux hij(t) sont numérisés (par exemple sur 9 bits, ou autre), échantillonnés (par exemple avec un échantillonnage à 20 MHz pour des ondes ultrasonores) et 15 enregistrés par l'unité centrale 5. L'unité centrale 5 procède ensuite à une étape de traitement (c), comprenant une sous-étape initiale (cl) de détermination de matrice de transfert fenêtrée. Au cours de cette sous-étape (cl), chaque réponse impulsionnelle h,i(t) est d'abord tronquée (fenêtrée) en fenêtres de temps successives, de durée At. On obtient ainsi une série de matrices temporelles fenêtrées K(T,t)=[klj(T,t)] de dimension N*N, où k;i(T,t) 10 est la contribution à hi3(t) correspondant à la fenêtre de temps [T-At/2 ;T+Ot/2], soit : kIJ(T,t)=hi(t)xWR(t-T) (1) avec : WR(t)=1 pour re [T-Ot/2;T+At/2] et WR(t)=0 partout 15 ailleurs, - R=c.T/2 est une distance mesurée depuis l'ensemble 3 de transducteurs (voir la figure 2), selon un axe z perpendiculaire à la direction longitudinale x dudit ensemble de transducteurs (constituant un réseau linéaire 20 dans l'exemple considéré ici), c est la célérité de l'onde dans le milieu 1 (soit dans le cas particulier des ultrasons dans l'eau ou dans des tissus mous du corps humain, c=1.5 mm/us), At est une durée prédéterminée correspondant à 25 une distance AR=c.At parcourue par l'onde, les intervalles de temps susmentionnés [T-At/2 ;T+At/2] pouvant avantageusement être jointifs, c'est-à-dire recouvrir entièrement une période de temps déterminée (autrement dit, le traitement effectué couvre une profondeur de champ, 30 c'est-à-dire une plage de valeurs de R qui est également continue). Au cours de la sous-étape (cl), on procède ensuite à une transformée de Fourier discrète des coefficients de la matrice K(T,t), pour obtenir ainsi pour chaque valeur de 16 T, une matrice de transfert de dimension N*N que l'on appellera matrice de transfert fréquentielle fenêtrée K(T,f)=[ki;(T,f)], où ki3(T,f) est la transformée de Fourier discrète de k(T,t) et f est la fréquence. A partir de ces matrices de transfert fréquentielles fenêtrées K(T,f), on peut séparer les contributions de diffusion simple et multiple, par filtrage, au cours d'une sous-étape ultérieure de filtrage (c3) appartenant à l'étape de traitement (c). En particulier, au cours de cette sous-étape ultérieure de filtrage (c3), on peut séparer la composante de diffusion multiple de la composante de diffusion simple dans chaque matrice de transfert fréquentielle fenêtrée K(T,f) en fonction de la cohérence des coefficients kij(T,f) de la matrice de transfert fréquentielle fenêtrée K(T,f) sur chaque antidiagonale de ladite matrice de transfert fréquentielle fenêtrée K(T,f) (on appelle antidiagonale un alignement de coefficients kij(T,f) de ladite matrice, tel que i+j soit constant). En effet, les ondes simplement diffusées présentent une cohérence particulière selon les antidiagonales de la matrice K(T,f) alors que les ondes multiplement diffusées ont un aspect aléatoire et ne présentent pas de direction de cohérence privilégiée dans ladite matrice K(T,f). En filtrant judicieusement ces antidiagonales, les deux contributions peuvent ainsi être séparées. Cette propriété peut être expliquée comme suit. Chacune des réponses impulsionnelles hi;(t) peut se décomposer sous la forme suivante : hl~ (t) = hJ (t) + hl ' (t) (2) où hjs(t)et hw(t)correspondent respectivement aux signaux issus de la diffusion simple (S) et de la diffusion multiple (M). De la même manière, les coefficients k,10", f) de la 17 matrice de transfert fréquentielle fenêtrée K(T,f) peuvent se décomposer chacun sous la forme ku(T, f)=k; (T, f)+k,'(T, f) , où kW,� est la contribution de diffusion simple et k,.M(T, f) est la contribution de diffusion multiple. Chaque contribution kJJ(T,f) peut être considérée comme la somme d'ondes partielles associées à plusieurs chemins de diffusion simple dont deux exemples (chemins dl et d2) sont représentés sur la figure 2. Les ondes simplement diffusées qui parviennent au réseau de capteurs dans la fenêtre de temps [T-At/2 ; T+At/2] correspondent à la réflexion de l'onde initiale sur des diffuseurs situés dans une zone du milieu d'épaisseur AR=c.At, à la profondeur R=cT/2. On parle de �volume isochrone� pour désigner l'ensemble des points qui, à un instant T donné, contribuent au signal capté sur le réseau. En réalité, le volume isochrone n'est pas exactement une tranche parallèle à la surface de la barrette mais résulte de la superposition d'ellipses ayant pour foyers l'élément émetteur (i) et l'élément récepteur (j). En champ lointain, c'est-à-dire lorsque R est suffisamment grand, on assimile le volume isochrone à une tranche d'épaisseur JR parallèle au réseau et distant de celui-ci de R=cT/2. La réponsek;Jentre les éléments i et j peut se décomposer en une somme d'ondes partielles issues de la réflexion sur les Nd diffuseurs du volume isochrone. A deux dimensions, dans le cadre de l'approximation paraxiale, on peut écrire kus(T,f) sous la forme suivante : k (T,.f)= exp 2RjkR)1Adexp j2 [(xi ùXd)2+(xi ùXd)2 (3) d=l R où l'entier d désigne le dème chemin de diffusion simple contribuant au signal reçu au temps T, Xd est la position transverse du diffuseur d (selon l'axe X), k est le nombre d'onde dans le milieu ambiant et Ad est une amplitude caractérisant la réflectivité du diffuseur d. La contribution de diffusion multiple peut elle aussi se décomposer en ondes partielles correspondant à des chemins de diffusion multiple de longueur comprise dans l'intervalle [R-AR/2 ; R+AR/2], comme représenté sur la figure 3 (pour deux chemins d'indices p et q). On peut écrire k(T,f) sous la forme suivante : k'(T,f~= exp[jk(Ztp)+Zzp>tp exp j k p) (x, -X1(P) )2 exp j kp) (x~ -Xzp)y p=i jZj p)Z(P) 2Zt 2Z2 (4) où l'entier p désigne l'indice du chemin de diffusion multiple considéré. Les couples (Xp),ZI(p�et (Xp),Zzp)) désignent respectivement les coordonnées du premier et du dernier diffuseur du chemin p dans l'exemple représenté sur la figure 3. BP correspond à l'amplitude complexe de l'onde partielle liée au chemin p, depuis le premier événement de diffusion (X, Zip)) jusqu'au dernier (X zp) , z ) . Bien que la répartition des diffuseurs 2 dans le milieu 1 soit totalement aléatoire et sans corrélation entre diffuseurs, le signal associé à la diffusion simple k présente une cohérence particulière, contrairement à la contribution de diffusion multiple. En effet, l'équation (3) peut se réécrire sous la forme : Nd ksi (T, f) = exp(R 'kR)exp j k (x, - xj) x Y Ad exp j k {(xi + xj - 2Xd)211 d=I 4R contribution déterministe contribution aléatoire (5) Le terme apparaissant devant la somme de l'équation (5) est indépendant de la répartition exacte des diffuseurs, il s'agit donc d'une contribution déterministe, 30 qui caractérise la diffusion simple. Le terme de droite est25 19 aléatoire, car il dépend explicitement de la position des diffuseurs. A contrario, le signal lié à la diffusion multiple (Equation 4) ne peut pas s'écrire ainsi. Cette propriété des signaux issus de la diffusion simple se traduit par une cohérence particulière suivant les antidiagonales de chaque matrice K(T,f), comme l'illustre la figure 4. En effet, le long de chaque antidiagonale, c'est-à-dire pour les couples émetteur(i)/récepteur(j) tels que i+j=constante, la contribution aléatoire de ksi est constante. Pour un milieu diffuseur donné, en régime de diffusion simple, il existe donc une relation de phase déterministe entre les coefficients de chaque matrice K(T,f) situés sur la même antidiagonale. En revanche, en régime de diffusion multiple, cette propriété n'est plus vérifiée et la matrice K(T,f) ne présente pas de cohérence particulière : les éléments kmsont indépendants les uns des autres. La présente invention permet de tirer profit de cette propriété pour isoler les contributions de diffusion simple et multiple par filtrage à partir des signaux mesurés expérimentalement, en exploitant la symétrie particulière de la contribution de diffusion simple au sein de chaque matrice K(T,f). On peut ainsi extraire par filtrage : - la composante de diffusion simple ksij(T,f) de chaque coefficient kiôT,f) de la matrice K(T,f), telle que les composantes de diffusion simple ks,j(T,f)soient cohérentes entre elles sur chaque antidiagonale de ladite matrice K(T,f), - et la composante de diffusion multiple kMij(T, f) de chaque coefficient kij(T,f) de la matrice K(T,f), telles que les composantes km,J(T,f) soient aléatoires et non corrélées entre elles. Deux exemples de techniques de filtrage, pouvant être utilisées pour séparer les deux contributions, sont présentés ci-après. Dans ces deux techniques, l'étape de traitement (c) comporte les deux sous-étapes suivantes, à la suite de la sous-étape (cl) de détermination de matrice de transfert fenêtrée : - une sous-étape (c2) de rotation de données, par 10 rotation des données de chaque matrice K(T,f) et création de deux sous-matrices A1(T,f) et A2(T,f), - la sous-étape (c3) susmentionnée de filtrage, au cours de laquelle on filtre les matrices Al et A2, notamment par projection (technique 1) ou par décomposition 15 en valeurs singulières (technique 2), et on obtient alors deux matrices filtrées, notées AÎet Al, - une sous-étape de rotation de données inverse (c4), permettant de reconstituer, à partir de Ai et Al, de matrices de transfert filtrées KF(T,f), contenant soit la 20 contribution de diffusion simple soit la contribution de diffusion multiple, selon le type de filtrage employé. Ces sous-étapes (c2) à (c4) sont détaillées ci-après. 25 Sous-étape (c2) : Rotation des données Au cours de cette sous-étape (c2), l'unité centrale 5 calcule deux matrices A1(T,f)=[aluv(T,f)] et A2=[a2uv(T,f)] à partir de chaque matrice K(T,f), où 30 aluv (T, f) =ku+V 1, v-u+2M-1 (T, f) , a2uv(T, f) =ku+v, v-u+2M-1 (T, f) , M=(N+3)/4. N est choisi de façon que M soit un entier : par exemple, N=125, et M=32. Si le nombre total de 21 transducteurs 3 est tel que M n'est pas un entier, on travaille sur un nombre réduit N de transducteurs, tel que M=(N+3)/4 soit entier (dans l'exemple particulier considéré ici, on peut par exemple utiliser une barrette de 128 transducteurs et travailler sur uniquement N=125 d'entre eux). Les matrices Al et A2 sont des matrices carrées constituées de sous-ensembles de la matrice K(T,f) qui sont pivotés de 45° dans le sens antihoraire. Ces matrices Al et A2 sont schématisées respectivement sur les figures 4 et 5, qui doivent être comprises comme suit : - Al est une matrice carrée de taille (2M-1) x (2M- 1) dont les lignes sont constituées par certaines diagonales de la matrice K(T,f) (une sur deux), comprenant les coefficients de la matrice K(T,f) qui sont entourés sur la figure 4 (les coefficients de K(T,f) sont schématisés par des points sur les figures 4 et 5), - A2 est une matrice carrée de taille (2M-2) x (2M- 2) dont les lignes sont constituées par les autres diagonales de la matrice K(T,f), comprenant les coefficients de la matrice K(T,f) qui sont entourés sur la figure 5. Dans ce qui suit, on appellera Ar=[ar,j] les matrices Al et A2 (r=1 ou 2) et on appellera L la dimension 25 de la matrice Ar (pour la matrice Al, on a donc L =2M -1 et pour la matrice A2, L = 2M - 2) . Du fait de la réciprocité spatiale, la matrice K est symétrique par rapport à sa diagonale principale D (kÿ= k). La matrice Ar présente donc elle aussi une 30 symétrie : chaque ligne sa partie supérieure est identique à une ligne de la partie inférieure, symétrique par rapport à une médiane horizontale correspondant à la diagonale principale D. La partie supérieure de la matrice Ar peut donc être déduite directement de la partie inférieure. Ainsi chaque colonne de la matrice A ne comporte que M coefficients indépendants bien qu'elle soit de dimension L>M. Sous-étape (c3) : Filtrage Au cours de la sous-étape de filtrage (c3), l'unité centrale 5 sépare la composante de diffusion multiple de la composante de diffusion simple dans chacune des matrices Ar, r étant un indice égal à 1 ou 2, et on obtient ainsi au moins deux matrices filtrées ArF correspondant respectivement aux deux matrices Ar et représentatives chacune, soit de la composante de diffusion simple, soit de la composante de diffusion multiple. Ce filtrage peut être réalisé notamment selon la technique 1 susmentionnée ou selon la technique 2 susmentionnée. 1. Technique 1 : Filtrage par projection Dans cette première technique, au cours de la sous-étape de filtrage (c3), l'unité centrale 5 calcule deux matrices filtrées ArF qui sont représentatives de la diffusion simple. Chacune de ces matrices filtrées est calculée par la formule : ArF =S`S*Ar , où : - S=[su] est un vecteur colonne, tS* étant la transposée du vecteur conjugué du vecteur S, - les composantes su du vecteur S sont des nombres k su = s.exp j 2R yû }/' - s est une constante, complexes égaux 23 - yu =(x;-xJY-/î , avec u=(i-j)2+M si r=1 et u=(i-j- 1)/2+M si r=2, - xi et xj sont les abscisses des transducteurs d'indices i et j selon un axe X, l'ensemble de transducteurs s'étendant au moins selon ledit axe X, - L=2M-1 pour r=1 et L=2M-2 pour r=2. Cette formule de calcule se justifie comme suit. Chaque matrice Ar est la somme deux termes, Ars et ArM, désignant respectivement les contributions dues à la diffusion simple et à la diffusion multiple : Ar= Ars+ ArM (6) La rotation des données, i.e. le passage de K(T,f) à Ar, se traduit mathématiquement par le changement de coordonnées (xi, x;) -> (yu, yv) yu = (x, - x; )/N/i, et y, = (x; + xj )/-N/i L'équation (5) se réécrit alors dans cette nouvelle base 2 arä(Z, f) ex = P kR exp R Nd krr où rV=Adexp J2RL d=1k 2 2R x (7) contribution aléatoire ù'XdY]} . contribution déterministe Ainsi pour un milieu 1 donné, chaque colonne de la 20 matrice Ars présente une dépendance suivant l'indice des lignes (u) parfaitement déterminée. En revanche, la contribution de diffusion multiple (Equation 4) ne peut se factoriser aussi simplement. Même après rotation de la matrice, le caractère aléatoire de la 25 position des diffuseurs persiste sur les colonnes comme sur les lignes de la matrice ArM. Le filtrage des signaux simplement diffusés peut donc être réalisé en projetant les colonnes de la matrice totale Ar sur l'espace �caractéristique de la diffusion 30 simple�, engendré par le vecteur S de coordonnées : s ä = s. exp j k Y M (8) , 2R où s est une constante, le cas échéant égale à 1. La présence au dénominateur de assure la normalisation du vecteur S. Le vecteur ligne P, résultat de 5 cette projection s'écrit : P=tS`Ar (9), où tS* est la transposée du vecteur conjugué de S. Les coordonnées du vecteur P sont données par : p s M exp(2 jkR)r . y (10 ) vù ~ suauv ù suauv+~ Su auvù y R y+~suauv u u u u 10 Le terme résiduel lsua m correspond à la projection des u signaux multiplement diffusés sur le vecteur S. Ensuite, la matrice filtrée ArF est obtenue en multipliant le vecteur colonne S par le vecteur ligne P : ArF =SP=S.tS".A (11) 15 Les coordonnées de la matrice ArF s'écrivent alors : ar,, = s,,p, = exp(R R kR)exp{ j y; r + sä s~ ,aM (12) Le premier terme est strictement égal à la composante simplement diffusée (Equation 7). On a donc : aû =a v+sulsu'auv (13) 20 En termes de matrices, l'équation (13) se réécrit de la manière suivante: ArF = Ars + S.t S*.AM diffusion bruit résiduel (14) simple La matrice ArF contient bien la contribution liée à la diffusion simple (As). Mais elle contient également un 25 terme résiduel lié à la présence de diffusion multiple (StS*AM). La persistance de ce terme est due au fait que les signaux de diffusion multiple ne sont pas strictement 25 orthogonaux à l'espace caractéristique de la diffusion simple, engendré par le vecteur S. Le filtrage réalisé n'est donc pas parfait, cependant, on peut évaluer l'importance du bruit résiduel. En effet, comme on l'a vu au paragraphe concernant la rotation des données, chaque colonne de la matrice Ar ne dispose que de M coefficients indépendants ; la contribution de diffusion multiple est donc diminuée d'un facteur suite au filtrage. La contribution de diffusion simple restant inchangée, le gain en rapport �signal-surbruit�, ou plus exactement �diffusion simple-sur-diffusion multiple� est donc de l'ordre de fi-g. La technique de filtrage décrite ci-dessus (technique 1) est à utiliser notamment quand on souhaite extraire une contribution de diffusion simple noyée dans de la diffusion multiple, c'est-à-dire dans le cas de milieux pour lesquels les signaux simplement diffusés sont d'amplitude très faible par rapport aux signaux issus de diffusion multiple. Cela s'applique notamment au cas de la détection de cibles enfouies dans un milieu diffusant. 2. Technique 2 : Filtrage par décomposition en valeurs singulières Cette seconde technique consiste à séparer diffusion simple et multiple en effectuant la décomposition en valeurs singulières ou SVD (�Singular Value Decomposition�) des matrices Al et A2 obtenues après rotation. La SVD a en effet la propriété de décomposer une matrice en deux sous-espaces : un espace dit �signal� (matrice caractérisée par une importante corrélation entre lignes et/ou colonnes de la matrice) et un espace dit �bruit� (matrice d'aspect aléatoire, sans corrélation entre éléments). En appliquant la SVD aux matrices Ar obtenues après rotation, l'espace signal correspond à la matrice Ars (contribution de diffusion simple, caractérisée par une grande corrélation suivant ses colonnes) et l'espace bruit est associé à la matrice ArM (contribution de diffusion multiple), avec Ar= Ars+ ArM (Equation 6 déjà citée au paragraphe concernant la technique 1). La SVD des matrices Ar s'écrit de la manière suivante : L Ar=UA`V*=12,U;(15) U et V sont des matrices carrées unitaires de dimension L. Leurs colonnes respectives Ui et Vi correspondent aux vecteurs propres associées à la valeur singulière A est une matrice carrée diagonale de dimension L dont les éléments diagonaux correspondent aux valeurs singulières a,i,r rangées dans l'ordre décroissant. Dans le paragraphe concernant la rotation de données, une symétrie particulière de la matrice Ar a été mise en évidence : cette matrice ne comprend que M lignes indépendantes, elle est donc de rang M La diffusion simple se caractérisant, après rotation des données, par une grande cohérence le long des colonnes des matrices Ar, la SVD fait ressortir cette contribution dans l'espace signal (la contribution diffusion simple sera associée aux valeurs singulières les plus élevées) alors que la contribution de diffusion multiple sera associée aux valeurs singulières plus faibles. Ici, contrairement à la première technique de filtrage, on ne fait donc pas d'hypothèse a priori sur la forme de la cohérence existant sur les antidiagonales de la matrice K(T,f) dans le cas de la diffusion simple, on suppose simplement que cette cohérence existe. 27 Le problème est de déterminer à quel rang de valeur singulière correspond le seuil séparant l'espace �signal� (associé à la diffusion simple) de l'espace �bruit� (associé à la diffusion multiple). Si l'équation (5) était rigoureusement exacte, seule la première des valeurs singulières correspondrait à l'espace signal. Lorsque les hypothèses aboutissant à l'équation (5) ne sont pas rigoureusement exactes, la contribution de diffusion simple n'est pas de rang 1 et plusieurs valeurs singulières portent la trace de cette contribution. Il faut alors établir un critère de séparation entre les contributions de diffusion simple (espace signal) et de diffusion multiple (espace bruit). Pour ce faire, on utilise les résultats de la théorie des matrices aléatoires. Par convention et souci de simplification, les valeurs singulières Ai,r sont normalisées par leur moyenne quadratique : 2i,r 1 M (17) 'Ym,r M+1ù1 m=i Pour une matrice de grande dimension dont les 20 coefficients sont totalement aléatoires, sans corrélation entre eux, la première valeur singulière 2lne dépasse jamais une valeur 2max (2max =2, dans le cas d'une matrice carrée). Expérimentalement, la contribution de diffusion 25 multiple ne correspond pas tout à fait à une matrice totalement aléatoire, car il existe des corrélations résiduelles entre capteurs (notamment à cause du couplage mécanique ou électrique entre transducteurs voisins de l'ensemble 3), ce qui modifie /lrä On peut établir à ax. 30 partir de [A.M. Sengupta and P.P. Mitra, "Distributions of singular values for some random matrices," Phys. Rev. E, Vol.60(3), pp 3389-3392, 1999.], la nouvelle loi de probabilité pour les valeurs singulières d'une telle matrice et en déduire la valeur de i,,, qui va permettre de définir un critère de séparation objectif entre l'espace 5 signal et l'espace bruit. Après rotation des données expérimentales, la matrice Ar à traiter (cf équation 6) est donc la somme d'une matrice Ars de rang p M Af,M = UitV1 (20) i=p+l La matrice Ars contient alors la contribution de diffusion simple (plus une contribution de diffusion multiple 10 résiduelle) et la matrice ArM est associée à la diffusion multiple. On notera que la technique 2 suppose que la première des valeurs singulières normalisée dépasse le seuil 2max. Elle ne peut pas être utilisée dans les milieux 15 hautement diffusants, c'est-à-dire des milieux pour lesquels la contribution de diffusion multiple est prépondérante par rapport à la diffusion simple. Dans ce cas, on utilisera plutôt la technique de filtrage par projection des antidiagonales sur l'espace de diffusion 20 simple (technique 1) pour extraire la contribution de diffusion simple. Si au contraire la contribution de diffusion simple est prépondérante ou du même ordre de grandeur que la diffusion multiple, on pourra utiliser la technique de filtrage par SVD des matrices A (technique 2) 25 et ainsi extraire la contribution de diffusion multiple. Sous-étape (c3) : rotation de données inverse Au cours de la sous-étape de rotation de données 30 inverse (c4), l'unité centrale 5 effectue une 30 transformation inverse de la transformation décrite à la sous-étape (cl), et calcule ainsi une matrice de transfert fenêtrée filtrée KF(T,f)=[kF (T,f)], où : - lorsque i-j est pair : kF1, (T, f) =alF(i_7)/2+M, (i+;) /2, - lorsque i-j est impair : kF,i (T, f) =a2F(I_3_1)/2+M, (1+3-1)/2. La matrice KF(T,f) est une matrice carrée de taille (2M-1)x (2M-1), contenant des signaux qui sont issus soit de la diffusion simple, soit de la diffusion multiple, selon le type de filtrage appliqué. La procédure de rotation inverse des données est schématisée sur la figure 6, qui montre que la matrice KF(T,f) reprend une partie des coefficients des antidiagonales des matrices Al(T,f) et A2(T,f). Les matrices filtrées KF(T,f) peut ensuite être 15 utilisée de diverses façons : si KF(T,f) correspond à la composante de diffusion simple, elle peut être utilisée notamment pour détecter un point singulier du milieu 1 ou pour bâtir une image du milieu 1. A cet effet, deux méthodes d'imagerie peuvent par 20 exemple être utilisées : - Selon une première méthode d'imagerie, l'unité centrale 5 peut par exemple calculer à partir chaque matrice KF(T,f), par transformée de Fourier inverse, une matrice temporelle fenêtrée filtrée de réponse 25 inter-éléments HF(T,t) puis sommer les différentes matrices temporelles HF(T,t) pour obtenir une matrice temporelle filtrée de réponse inter-éléments HF(t), représentative des réponses entre transducteurs de l'ensemble de transducteurs. L'unité centrale réalise 30 alors l'image du milieu par formation de voies à partir de ladite matrice temporelle filtrée de réponse inter-éléments HF(t). > Selon une deuxième méthode d'imagerie, dite méthode de décomposition de l'opérateur de retournement temporel 35 (�méthode DORT�, définie notamment par Prada et al., 31 �Eigenmodes of the time reversai operator : a solution to selective focusing in multiple-target media�, Wave Motion 20 (1994) 151-163, Elsevier Science B.V.), l'unité centrale 5 peut procéder aux étapes suivantes : (dl) déterminer pour chaque matrice de transfert filtrée KF(T,f) un opérateur de retournement temporel ORT(T,f) = KF*(T, f) KF(T,f), correspondant à la ou chaque matrice de transfert KF(T,f), KF*(T,f) étant la matrice conjuguée complexe de KF(T,f), (d2) déterminer au moins des vecteurs propres et/ou des valeurs propres dudit opérateur de retournement temporel, (d3) détecter au moins une singularité dans le milieu en fonction au moins desdits vecteurs propres et/ou desdites valeurs propres. Au cours de cette étape (d3), on peut notamment déterminer des vecteurs propres Vi(f) de l'opérateur de retournement temporel, puis propager numériquement ces vecteurs propres dans un modèle numérique du milieu 1, pour déterminer une position d'un objet réverbérant dans le milieu 1 correspondant à chaque vecteur propre. - si les matrices KF(T,f) correspondent à la composante de diffusion multiple, elles peuvent être utilisées notamment pour calculer à partir de ladite composante de diffusion multiple, un indice représentatif de l'importance de la diffusion multiple dans le milieu. Dans ce cas, on peut par exemple réaliser une image (échographique ou autre) du milieu par tout moyen connu ou par les méthodes susmentionnées, et on quantifie la fiabilité de ladite image en fonction dudit indice représentatif de l'importance de la diffusion multiple. Avantageusement, on peut calculer cet indice en plusieurs points (notamment à diverses profondeurs R) et on quantifie la fiabilité de plusieurs parties de l'image correspondant 32 auxdits plusieurs points, en fonction dudit indice représentatif de l'importance de la diffusion multiple. Plusieurs exemples d'application possibles du procédé selon l'invention vont maintenant être décrits en détail. Application 1 : détection d'une cible dans un milieu diffusant L'application présentée ici consiste à détecter une cible (un gros réflecteur) enfouie dans un milieu diffusant (milieu très concentré en petits diffuseurs) en isolant la contribution de diffusion simple issue de la cible via la technique de filtrage par projection (technique de filtrage 1 susmentionnée). Le dispositif expérimental est représenté sur la figure 8. L'ensemble 3 de transducteurs tel que décrit ci-dessus en regard de la figure 1, est placé dans un bain d'eau 1 formant le milieu 1. Le milieu 1 contient un grand nombre de diffuseurs 2 disposés par exemple à une distance R d'environ 40 mm de l'ensemble de transducteurs 3 et sur une profondeur de champ d'environ 20 mm. Cette distance R et cette profondeur de champ sont mesurées selon l'axe Z. Les diffuseurs peuvent par exemple être constitués par un ensemble de tiges d'acier de 0.8mm de diamètre réparties aléatoirement avec une concentration d'environ 12 tiges/cm2. Le libre parcours moyen le typique pour ce milieu 1 est de 7 mm, de sorte que ce milieu est fortement diffusant. La cible 7 est disposée par exemple environ 10 mm derrière l'ensemble des diffuseurs 2. Cette cible 7 peut par exemple être une tige d'acier d'un diamètre de 10 mm. La matrice temporelle de réponse inter-éléments H(t) est déterminée comme suit : on fait émettre par les transducteurs un signal balayant linéairement une bande de fréquence de 2 à 4 MHz (�chirp�). Les signaux réverbérés 33 par le milieu 1 sont alors enregistrés par les transducteurs 3 et sont ensuite convolués par le signal d'émission afin d'obtenir la réponse impulsionnelle entre chaque paire de transducteurs. L'ensemble de transducteurs peut se présenter sou la forme d'une barrette échographique de 128 voies avec un pas entre inter-éléments de 0.417mm. La matrice H(t) est donc de taille 128x128. Celle-ci est ensuite découpée temporellement selon l'équation (1) avec des fenêtres temporelles de longueur dt=5.5ps. Chaque fenêtre se recouvre de moitié avec ses voisines. Une série de matrice K(T,t) est obtenue et le passage dans le domaine fréquentiel est ensuite réalisé par une transformée de Fourier discrète, ce qui aboutit à l'obtention d'une série de matrices K(T,f). Dans la suite de l'étude, seule la bande de fréquence [2.5 - 3.5] MHz sera considérée. Les figures 9 et 10 montrent les images du milieu 1 obtenues à partir de la matrice H(t) (avec une échelle linéaire pour l'amplitude du signal sur la figure 9 et avec une échelle logarithmique sur la figure 10) par les méthodes classiques d'échographie, par formation de voies focalisée à l'émission et à la réception. Il est impossible de distinguer détecter la cible, située ici à R=70mm. On détecte clairement les premiers diffuseurs du milieu désordonné (sur environ lcm d'épaisseur, soit une épaisseur comparable au libre parcours moyen) mais ensuite, les signaux multiplement diffusés prennent le dessus et une image de tavelures (speckle) est alors obtenue, sans correspondance physique avec la présence de cibles. Par ailleurs, on a également essayé de traiter les signaux obtenus par la méthode DORT susmentionnée, sans le filtrage selon l'invention. Celle-ci consiste en une décomposition en valeurs singulières de chaque matrice K(T,f) : K=UA`V* (21) , 34 avec les notations déjà définies ci-dessus. Il est connu qu'en l'absence de bruit et pour des cibles ponctuelles en régime de diffusion simple, chaque valeur singulière non nulle est associée à une cible du milieu. A chaque fréquence, la valeur singulière Xi est proportionnelle à la réflectivité de la sème cible et le vecteur propre V; décrit le signal qui permet de refocaliser une onde sur le diffuseur en question. La méthode DORT permet donc de distinguer différentes cibles et de focaliser l'onde sélectivement sur chaque cible. De plus, si une cible est détectable dans une bande de fréquence et pas dans une autre, on peut sélectionner la bande de fréquence dans laquelle on détecte la cible afin de construire son image, alors qu'avec l'échographie traditionnelle, toute la bande de fréquence du signal incident est utilisée pour reconstruire l'image du milieu ce qui dégrade la qualité de l'image finale. La méthode DORT a été testée ici sur la série de matrices K(T,f) : une décomposition en valeurs singulières (SVD) est effectuée à chaque temps d'écho T et à chaque fréquence f. Pour chaque matrice K(T,f), la première valeur singulière est normalisée par la moyenne quadratique des valeurs singulières : /t= (22) -\1KÂ) où représente la lére valeur singulière normalisée. Cette normalisation a pour but de définir un critère de détection de cible : si Âl Xseuir , avec une probabilité de fausse alarme donnée, et dépendant de la valeur ;i,.i : si X, 36 Dans la configuration expérimentale étudiée, nous avons retenu ici la valeur seuil Xseuii=3, qui correspond à une probabilité de fausse alarme inférieure à 0.1%. Sur les figures 11 et 12 est présentée l'évolution de Âlavec la distance R (correspondant au temps d'écho T tel que R=cT/2) et la fréquence f, sans seuillage pour la figure 11 et avec seuillage à la valeur seuil Xseuii=3 pour la figure 12. Contrairement à ce qui a été obtenu dans les mêmes conditions avec la méthode DORT seule (aucune cible n'avait pu être détectée), le filtrage de la diffusion simple permet la détection de la cible puisque la première valeur singulière 'il dépasse la valeur seuil ÂS,f sur une certaine plage de fréquence au temps d'écho de la cible. Une fois le critère de détection appliqué, l'image de la cible peut être obtenue (figure 13) en repropageant numériquement les vecteur propres Vi à des profondeurs R(=cT/2) et des fréquences f pour lesquelles la première valeur singulière satisfait le critère de détection. On obtient ainsi une image de la cible, qui permet de la localiser, ce que ni l'échographie ni la méthode DORT seule n'ont pu faire. Sur la figure 14, on a comparé l'amplitude normalisée du signal de l'image dans le plan de la cible 25 7 : - dans le cas de l'image de la figure 13 (courbe en pointillés), ce qui correspond au filtrage de la contribution de diffusion simple couplé à la méthode DORT), - et dans le cas de l'image de la figure 9 (courbe 30 en trait plein), ce qui correspond à l'échographie classique. On peut constater que le signal obtenu en mettant en oeuvre la présente invention (courbe en pointillés) est nettement meilleur et permet de discriminer aisément la 37 cible 7, contrairement au signal obtenu par échographie classique. On peut montrer en effet que l'amélioration du rapport signal sur bruit sur l'image finale croît en N pour l'échographie et en N3/2 avec la technique de l'invention, N désignant le nombre de transducteurs du réseau 3. L'image échographique (trait plein) présente des lobes secondaires assez importants, dus à l'aberration induite par le milieu diffusant. A partir d'une telle image, on ne peut pas conclure avec certitude que le pic principal soit effectivement lié à la présence d'une cible : statistiquement, il peut s'agir d'une fausse alarme due aux fluctuations inhérentes à l'aspect aléatoire de la diffusion multiple. Pour atteindre une probabilité de fausse alarme inférieure à 0 .1%, on peut montrer qu'il faudrait ici que le pic principal ait une valeur 2,5 fois supérieure à la moyenne de l'amplitude de l'image échographique, ce qui n'est pas le cas. Enfin, le pic principal de l'image échographique ne se retrouve pas exactement à la bonne position ; au contraire, l'image obtenue avec l'invention ne semble pas souffrir d'aberration. En conclusion, la technique développée ici qui allie le filtrage de la diffusion simple avec la méthode DORT améliore les capacités d'un réseau d'émetteurs-25 récepteurs à détecter et imager une cible cachée derrière un écran diffusant multiplement les ondes. D'une part, l'analyse temps fréquence mise en jeu permet de sélectionner les bandes de fréquence favorables à la détection de la cible, sélection impossible avec 30 l'échographie traditionnelle. Une analyse théorique fondée sur la théorie des matrices aléatoires permet de montrer que, dans des conditions identiques, la probabilité de fausses alarmes en milieu diffusant est rendue plus faible par notre technique. D'autre part, cette technique réduit les effets d'aberration (apparition de lobes secondaires, déplacement du lobe principal...) dont souffrent les techniques d'imagerie conventionnelles. Application 2 : étude d'un milieu très faiblement diffusant Contrairement à l'application précédente, nous considérons ici le cas d'un milieu désordonné faiblement 10 diffusant : la contribution de diffusion simple domine alors la diffusion multiple. Dans ce type de milieu, les techniques conventionnelles comme l'échographie fonctionnent bien et permettent de construire des images de réflectivité du milieu car la diffusion multiple étant 15 minoritaire, celle-ci ne perturbe que très peu l'image échographique obtenue. Cependant, on cherche ici à avoir accès à d'autres informations intéressantes dans l'étude du milieu de propagation : notamment des paramètres caractérisant purement la partie multiplement diffusée de 20 l'onde. Ces paramètres sont les différents �libres parcours moyen� au sein du milieu ; une fois la contribution de diffusion multiple isolée grâce à notre technique, on peut mesurer ces paramètres, ce que ne peuvent faire les techniques classiques d'imagerie échographique. 25 Lorsqu'une onde se propage au sein d'un milieu diffuseur aléatoire, elle perd progressivement, de façon exponentielle, sa cohérence . au bout d'une distance L, seule une fraction exp(-L/1eXt) de l'énergie initiale continue de se propager en cohérence avec l'onde initiale. 30 Le paramètre 1ext, libre parcours d'extinction globale, caractérise donc la distance d'extinction de la partie cohérente de l'onde. Cette extinction progressive de l'onde cohérente a deux origines distinctes : la diffusion (à chaque rencontre avec un diffuseur, une partie de la 35 cohérence initiale est perdue), et l'absorption intrinsèque5 du milieu de propagation. A ces deux phénomènes sont rattachées deux autres longueurs caractéristiques : le libre parcours moyen élastique 1e, et le libre parcours moyen d'absorption 1a, de telle sorte que exp - L = exp - - x exp L , et donc 1 = 1 + 1 eext ~ ee) ea) eext ee ea effet des effet des pertes par pertes par diffusion absorption Mesurer la longueur d'extinction globale 1 ext est intéressant, mais ne permet pas de distinguer les pertes par absorption (la) des pertes par diffusion (le) . La technique de séparation diffusion simple/diffusion multiple selon l'invention, en sélectionnant dans la réponse du milieu la contribution de diffusion multiple, permet de mesurer de façon distincte le et 1a et ainsi de caractériser plus complètement le milieu sondé. Pour illustrer la faisabilité et l'intérêt de telles mesures, la technique de séparation diffusion simple/diffusion multiple selon l'invention a été appliquée tout d'abord à un milieu synthétique (gel agar-agar faiblement diffusant). Le montage expérimental est similaire à celui de la figure 1, mais l'échantillon de milieu diffusant est placé dans un bain d'eau à une distance a=50 mm de l'ensemble 3 de transducteurs, cet échantillon ayant une épaisseur L=100 mm selon l'axe Z. On extrait la contribution de diffusion multiple dans une série de matrices KF(T,f) obtenues par la technique 2 décrite ci-dessus. Une fois la contribution de diffusion multiple isolée, on calcule l'intensité multidiffusée moyenne 1'4 en fonction de l'écart émetteur-récepteur X (distance entre le transducteur émetteur i et le transducteur récepteur j) et du temps d'écho T : IMX=xi ùxi, k, (T f 2\ {f,(xäxr )} (24) Rappelons que k;"(T,f)2 désigne les coefficients (ième ligne, jème colonne) de la matrice KM(T, f) . Les symboles représentent une moyenne effectuée d'une part sur toute la bande de fréquence f et d'autre part sur tous les couples 5 source/récepteur (i,j) séparés de X=x3-x2. Par ailleurs, un résultat supplémentaire est fourni par l'examen de l'évolution temporelle des intensités simplement et multiplement diffusées obtenues au point source (X=0) , qu'on notera 18(0,T) et IM(0,T) . En effet, la 10 théorie du transfert radiatif [J. Paasschens, "Solution of the time-dependent Boltzmann equation, "Phys. Rev. E, Vol.56(l), pp 1135-1141, 1997] montre qu'il existe des solutions analytiques prédisant l'évolution temporelle de Is(O,T) et IM(0,T). En particulier, il apparaît que 15 l'évolution temporelle de l'intensité simplement diffusée Is(O,T) ne dépend que du libre parcours d'extinction 1ext- Dans le cas du gel agar-agar étudié ici, un ajustement entre la prévision théorique et le résultat de l'expérience (figure 15) aboutit à la valeur lext=50mm. 20 Pour le vérifier, on a mesuré l'intensité diffusée totale I(0, T) et l'intensité multiplement diffusée IM(0, T) comme expliqué ci-dessus, et on en a déduit l'intensité simplement diffusée Is(O,T) par différence. En variante, on peut mesurer IS(O,T) en appliquant la formule 24 non plus 25 aux coefficients de la matrice KM(T,f) mais aux coefficients de la matrice Ks(T,f) déterminée par exemple selon la méthode 2 décrite ci-dessus. La courbe mesurée de IS(O,T), représentée par des points sur la figure 15, a été normalisée par son maximum (fixé à 1). On a calculé par 30 ailleurs l'évolution temporelle de Is(O,T) prédite par la solution exacte de l'équation de transfert radiatif, pour plusieurs parcours moyens d'extinction. Cette courbe, normalisée également en fixant sa valeur maximale à 1, s'ajuste sur la courbe expérimentale pour un libre parcours d'extinction lext=50mm (courbe en trait plein sur la figure 15), ce qui permet de déterminer que lext=50mm. Si l'on considère à présent l'évolution temporelle de l'intensité multiplement diffusée IM(X=0,T), la théorie montre qu'elle dépend distinctement des libres parcours moyens le et la. Ainsi, le peut être déterminés par ajustement des mesures expérimentales de IM(X=O,T) avec la théorie (Figure 16), ce qui donne le =1000mm. On obtient ce résultat en déterminant quelle est la valeur de le pour laquelle la courbe théorique de IM(X=O,T) s'ajuste sur la courbe expérimentale (courbe EXP matérialisée par des points en forme de petits cercles sur la figure 16). On en déduit la par la formule I = i I ce qui donne la=50 mm. eext ee ea Le fait de pouvoir séparer contribution de diffusion simple et contribution de diffusion multiple permet donc de mesurer de façon distincte les pertes par absorption des pertes par diffusion. Ici, l'échantillon étudié se trouve être beaucoup plus absorbant que diffusant puisqu'on trouve un libre parcours moyen élastique le--1000 mm alors que le libre parcours d'absorption 1a est de l'ordre de 50 mm. Cette expérience montre que l'invention permet de mieux caractériser le milieu diffusant, en mesurant séparément des paramètres diffusifs (le, la) ce que l'échographie traditionnelle ne permet pas de faire. Ici, on a considéré un cas extrême (gel faiblement diffusant) où le rapport IM / IS est particulièrement faible, mais la technique fonctionne également pour des milieux plus diffusants pour lesquels le rapport IM / IS est proche de l'unité. 42 Application 3 : imagerie ultrasonore du corps humain Le montage expérimental de la figure 1 a été utilisé pour imager le sein d'une volontaire, par échographie classique. L'image obtenue est représentée sur la figure 17. En parallèle, on a utilisé la technique de filtrage 2 pour déterminer l'intensité multiplement diffusée IM, comme dans l'application 2 ci-dessus et calculer ensuite un rapport IM/I, où I est l'intensité diffusée totale, en fonction de la distance R mesurée selon l'axe Z (c'est à dire en fonction du temps t d'aller-retour de l'onde simplement diffusée). Les valeurs du rapport IM/I, qui est un indice de fiabilité compris entre 0 et 1 représentatif de la quantité de diffusion multiple -et donc représentatif également de l'importance de la diffusion simple dans le milieu 1- ont été traduites en niveau de gris dans la barrette située en partie supérieure de la figure 17, qui permet ainsi de déterminer le rapport IM/I en fonction de t (Ce rapport pourrait être traduit en niveaux de couleurs ou par une courbe). Sur la figure 17, le trait vertical indique la limite de t (et donc de R) à partir de laquelle le rapport IM/I atteint un seuil critique (ici par exemple 0.5) au-delà duquel on ne peut plus faire confiance à l'image échographique, fondée sur le fait que la diffusion simple est prédominante. On peut donc ainsi fournir une aide aux praticiens en matière d'imagerie médicale ou autre, en leur donnant une information supplémentaire quant à la fiabilité de l'image. En effet, lorsque la diffusion multiple devient prépondérante, l'image échographique n'a plus de sens et la diffusion multiple peut être une source d'erreurs. Autres applications :35 43 Concernant l'extraction de la diffusion simple et sa possible application à la détection d'une cible enfouie dans un milieu diffusant, deux exemples d'applications supplémentaires peuvent être cités : - en contrôle non destructif (CND) : dans des matériaux comme le béton ou encore les aciers à grain, la diffusion multiple est loin d'être négligeable et la détection de fissures ou de défauts peut s'avérer difficile dans de tels milieux. L'extraction de la contribution de diffusion simple selon l'invention peut palier cet inconvénient. - la détection de cibles enfouies dans la terre (par exemple objets solides, cavités, poches de fluides...) : le sol est un milieu diffusant vis à vis des ondes utilisées pour détecter les objets enfouis. Le filtrage de la diffusion simple selon l'invention, couplé à la méthode DORT, donne de bien meilleurs résultats que l'échographie standard dans de tels milieux. - la technique de filtrage de la diffusion simple par projection suivant les antidiagonales peut également s'appliquer à des données tout simplement bruitées. C'est souvent le cas avec les données �réelles� comme en CND ou en acoustique sous marine où les rapports signal/bruit sont très faibles. Dans ce cas, la méthode DORT avec un filtrage préalable des antidiagonales donnera un résultat plus avantageux que celui obtenu en échographie en terme de probabilité de fausses alarmes. Cette méthode réduira également les effets d'aberration auxquels sont confrontées les techniques classiques d'imagerie. Concernant l'extraction de la diffusion multiple et la possibilité d'effectuer une caractérisation du milieu à partir de ses paramètres diffusifs, là encore les possibilités d'application sont variées : - l'imagerie médicale : souvent oubliée ou négligée 35 dans les tissus mous du corps humain, la diffusion multiple 44 des ultrasons existe bien. L'invention permet de l'extraire afin de caractériser les différents tissus sondés dans le corps humain. Cela permet de nouvelles mesures que l'échographie standard ne peut pas donner. Par exemple, en présence d'une tumeur, une vascularisation importante se met en place avec un réseau microscopique de vaisseaux pour alimenter celle-ci. On peut raisonnablement penser qu'une telle vascularisation conduit à une augmentation de l'intensité multiplement diffusée et serait donc détectable avec la technique de l'invention. Cela permettrait alors de quantifier l'angiogénèse des tumeurs et offrirait un meilleur diagnostic que celui qui est obtenu actuellement. Un autre exemple est donné par le suivi de l'ostéoporose. En effet les os, lorsqu'ils deviennent poreux, sont particulièrement diffusants. L'extraction de la diffusion multiple permettrait d'obtenir des mesures de paramètres diffusifs et d'imager la porosité locale de l'os. - le contrôle non destructif : des matériaux comme les aciers à grain, le béton et plus généralement les milieux granulaires, sont le siège de diffusion multiple. L'extraction de la diffusion multiple, même en faible quantité, permettrait alors de caractériser et/ou d'imager les matériaux à partir de leur comportement diffusif et relier cela à la microarchitecture du matériau sondé. - dans le domaine des mesures sismiques et études du sous-sol, la séparation des deux contributions de diffusion simple et multiple est également bénéfique. REVENDICATIONS1. Procédé de sondage par propagation d'ondes utilisant un ensemble (3) de transducteurs, ledit procédé 5 comprenant : (a) une étape d'émission au cours de laquelle on fait émettre par l'ensemble (3) de transducteurs une onde incidente dans un milieu (1) diffusant pour ladite onde, (b) une étape de mesure au cours de laquelle on 10 fait capter, par ledit ensemble (3) de transducteurs, des signaux représentatifs d'une onde réfléchie réverbérée par le milieu (1) à partir de l'onde incidente, lesdits signaux captés comprenant : - une composante de diffusion simple, 15 représentative de trajets d'onde où l'onde réfléchie résulte d'une seule réflexion de l'onde incidente par chaque diffuseur du milieu, - et le cas échéant une composante de diffusion multiple, représentative de trajets d'onde où l'onde 20 réfléchie résulte de plusieurs réflexions successives de l'onde incidente sur les diffuseurs du milieu avant d'atteindre l'ensemble de transducteurs, (c) une étape de traitement au cours de laquelle on traite lesdits signaux captés en vue de déterminer des 25 caractéristiques du milieu, caractérisé en ce qu'au cours de l'étape de traitement (c), on extrait, par filtrage, au moins une composante choisie parmi la composante de diffusion multiple et la composante de diffusion simple. 30 2. Procédé selon la revendication 1, dans lequel au cours de l'étape de traitement (c), on extrait la composante de diffusion simple et on utilise ladite composante de diffusion simple pour détecter au moins un point singulier dans le milieu. 35 3. Procédé selon la revendication 2, dans lequel au 46 cours de l'étape de traitement (c), on réalise une image du milieu à partir de ladite composante de diffusion simple. 4. Procédé selon la revendication 3, dans lequel au cours de l'étape de traitement (c), on détermine ladite composante de diffusion simple sous forme d'une matrice temporelle filtrée de réponse inter-éléments HF, représentative des réponses entre transducteurs de l'ensemble de transducteurs, et on réalise l'image du milieu par formation de voie à partir de ladite matrice temporelle filtrée de réponse inter-éléments HF. 5. Procédé selon la revendication 3, dans lequel au cours de l'étape de traitement (c), on détermine ladite composante de diffusion simple sous forme d'au moins une matrice de transfert fréquentielle filtrée KF, représentative des réponses entre transducteurs de l'ensemble de transducteurs, et on réalise l'image du milieu par la méthode de décomposition de l'opérateur de retournement temporel KF*.KF, où KF* est la matrice conjuguée complexe de ladite matrice de transfert fréquentielle filtrée KF. 6. Procédé selon la revendication 1, dans lequel au cours de l'étape de traitement (c), on extrait la composante de diffusion multiple et on calcule à partir de ladite composante de diffusion multiple, un indice de fiabilité représentatif de l'importance de la diffusion simple dans le milieu. 7. Procédé selon la revendication 6, dans lequel au cours de l'étape de traitement (c), on réalise une image du milieu et on quantifie la fiabilité de ladite image en fonction dudit indice de fiabilité. 8. Procédé selon la revendication 7, dans lequel au cours de l'étape de traitement (c), on calcule en plusieurs points ledit indice de fiabilité et on quantifie la fiabilité de plusieurs parties de l'image correspondant auxdits plusieurs points, en fonction dudit indice de 47 fiabilité. 9. Procédé selon l'une quelconque des revendications précédentes, dans lequel au cours de l'étape de traitement (c), on sépare la composante de diffusion multiple de la composante de diffusion simple par filtrage d'au moins une matrice de transfert fréquentielle représentative des réponses entre transducteurs de l'ensemble de transducteurs. 10. Procédé selon la revendication 9, dans lequel 10 l'étape de traitement (c) comporte au moins les sous-étapes suivantes : (cl) une sous-étape de détermination de matrice de transfert fenêtrée au cours de laquelle on détermine au moins une matrice de transfert fréquentielle fenêtrée 15 K(T,f)=[kii(T,f)] correspondant à une matrice temporelle fenêtrée de réponse inter-éléments K(T,t)=[ki](T,t)], ladite matrice temporelle fenêtrée de réponse inter-éléments correspondant, sur une fenêtre temporelle voisine d'un temps T et de durée At, aux réponses temporelles 20 hij(t) entre transducteurs de l'ensemble de transducteurs, f étant la fréquence, (c3) une sous-étape de filtrage au cours de laquelle on sépare la composante de diffusion multiple de la composante de diffusion simple dans chaque matrice de 25 transfert fréquentielle fenêtrée K(T,f) en fonction de la cohérence des coefficients kii(T,f) de la matrice de transfert fréquentielle fenêtrée K(T,f) sur chaque antidiagonale de ladite matrice de transfert fréquentielle fenêtrée K(T,f). 30 11. Procédé selon la revendication 10, dans lequel, au cours de la sous-étape de filtrage (c3), on extrait une composante de diffusion simple ksij(T,f) de chaque coefficient ki;(T,f) de la matrice K(T,f), telle que les composantes de diffusion simple ksii(T,f)soient cohérentes 35 entre elles sur chaque antidiagonale de ladite matriceK(T,f). 12. Procédé selon la revendication 10, dans lequel, au cours de la sous-étape de filtrage (c3), on extrait une composante de diffusion multiple kMi;(T,f) de chaque coefficient ki;(T,f) de la matrice K(T,f). 13. Procédé selon l'une quelconque des revendications 10 à 12, dans lequel l'ensemble (3) de transducteurs comporte au moins un nombre N de transducteurs et la matrice de transfert fréquentielle fenêtrée K(T,f) est une matrice de dimension N*N. 14. Procédé selon l'une quelconque des revendications 10 à 13, dans lequel, au cours de la sous-étape (cl) de détermination de matrice de transfert fenêtrée, on détermine plusieurs matrices de transfert fenêtrées K(T,f) pour plusieurs valeurs de T et on procède à la sous-étape de filtrage (c3) pour chacune des matrices de transfert fenêtrées K(T,f). 15. Procédé selon la revendication 14, dans lequel les fenêtres temporelles des différentes matrices de 20 transfert fenêtrées K(T,f) sont jointives. 16. Procédé selon l'une quelconque des revendications 10 à 15, dans lequel au cours de la sous-étape (cl) de détermination de matrice de transfert fenêtrée, on détermine chaque matrice fréquentielle 25 fenêtrées K(T,f) par transformée de Fourier d'une matrice temporelle fenêtrée K(T,t) correspondant, sur ladite fenêtre temporelle voisine du temps T et de durée zt, aux réponses temporelles hi;(t) entre transducteurs de l'ensemble de transducteurs. 30 17. Procédé selon l'une quelconque des revendications 10 à 16, dans lequel : - l'étape de traitement (c) comporte en outre une sous-étape de rotation de données (c2) qui est intermédiaire entre la sous-étape (cl) de détermination de 35 matrice de transfert fenêtrée et la sous-étape (c3) de 48filtrage, et au cours de laquelle on calcule deux matrices A1(T,f)=[aluä(T,f)] et A2=[a2,,(T,f)] à partir de K(T,f), où aluv (T, f) v-u+2M-1 (T, f) , a2uv, (T , f) =ku+v, v-u+2M-1 (T , f) , M=(N+3)/4, N est une plus petite dimension de la matrice de transfert K(T,f) et est choisi pour que M soit entier, - au cours de la sous-étape de filtrage (c3), on sépare la composante de diffusion multiple de la composante de diffusion simple dans chacune des matrices Ar, r étant un indice égal à 1 ou 2, et on obtient ainsi au moins deux matrices filtrées ArF correspondant respectivement aux matrices Ar et représentatives chacune, soit de la composante de diffusion simple, soit de la composante de diffusion multiple, - l'étape de traitement (c) comporte en outre une sous-étape de rotation de données inverse (c4), qui est postérieure à la sous-étape de filtrage (c3), et au cours de laquelle on calcule une matrice de transfert fenêtrée filtrée KF(T,f)=[kFii(T,f)], où : - lorsque i-j est pair : kF13 (T, f) =alF(i-3)/2+M, (i+;)/2, - lorsque i-j est impair : kFij (T, f) =a1F(i-D-1)i2+M, (i+j-1) /2 18. Procédé selon la revendication 17, dans lequel au cours de la sous-étape de filtrage (c3), on calcule les deux matrices ArF par la formule : ArF =S`S*Ar , où : - S=[s,] est un vecteur colonne, tS* étant la transposée du vecteur conjugué du vecteur S, - les composantes su du vecteur S sont des nombres complexes égaux à k su = s. exp j Yu }/-^/-É , 2R - s est une constante,- yu=(x;ûx,, avec u=(i-j)/2+M si r=1 et u=(i-j- 1)/2+M si r=2, xi et xj sont les abscisses des transducteurs d'indices i et j selon un axe X, l'ensemble de 5 transducteurs s'étendant au moins selon ledit axe X, L=2M-1 pour r=1 et L=2M-2 pour r=2, lesdites matrices ArF étant représentatives de la composante de diffusion simple. 19. Procédé selon la revendication 17, dans lequel 10 au cours de la sous-étape de filtrage (c3) : - on décompose en valeurs singulières chacune des matrices Ar en déterminant ainsi, pour chacune des matrices Ar, M valeurs singulières Ài,r, i étant un indice compris entre 1 et M, les valeurs singulières étant ordonnées par 15 ordre décroissant pour i croissant, - pour chaque matrice Ar, on attribue à la diffusion simple p valeurs singulières les plus importantes Al,r à Ap,r, p étant compris entre 1 et M, et on attribue à la diffusion multiple les éventuelles autres valeurs 20 singulières non nulles. 20. Procédé selon la revendication 19, dans lequel la détermination des p valeurs singulières attribuées à la diffusion simple est effectuée de manière itérative, comme suit . 25 (c31) on fixe initialement q = 1 ; (c32) on calcule une valeur singulière normalisée à partir de Àq,r : 2q,r rq,r 1 M a M+1ù q q m=q (c33) si Àq,r est au moins égal à une valeur seuil prédéterminée 2m, on attribue 2qr à la diffusion simple 30 et on réitère la sous-étape (c32) en incrémentant q d'une unité, 1 (c34) si /lq r est inférieur à la valeur seuil 2max , on attribue 2q,r et les éventuelles valeurs singulières suivantes à la diffusion multiple. 21. Procédé selon l'une quelconque des 5 revendications précédentes, dans lequel lesdites ondes sont choisies parmi les ondes mécaniques et les ondes électromagnétiques. 22. Dispositif pour la mise en oeuvre d'un procédé de sondage selon l'une quelconque des revendications précédentes, comprenant un ensemble (3) de transducteurs adaptés pour émettre une onde incidente dans un milieu diffusant (1) et pour capter des signaux représentatifs d'une onde réfléchie réverbérée par le milieu à partir de l'onde incidente, lesdits signaux captés comprenant : - une composante de diffusion simple, représentative de trajets d'onde où l'onde réfléchie résulte d'une seule réflexion de l'onde incidente par chaque diffuseur du milieu, - et le cas échéant une composante de diffusion multiple, représentative de trajets d'onde où l'onde réfléchie résulte de plusieurs réflexions successives de l'onde incidente sur les diffuseurs du milieu avant d'atteindre l'ensemble de transducteurs, ledit dispositif comportant en outre des moyens de traitement (5) adaptés pour traiter lesdits signaux captés en vue de déterminer des caractéristiques du milieu (1), caractérisé en ce que les moyens de traitement (5) sont adaptés pour extraire, par filtrage, au moins une composante choisie parmi la composante de diffusion multiple et la composante de diffusion simple.