Brevet FR3160773A1 — publié le 20251003

L’invention concerne un procédé d’analyse d’un milieu incluant l’acquisition (101) de signaux échos le calcul (102) d’une pluralité de signaux directionnels à partir des signaux échos l’estimation (103) d’un ensemble de signaux inter-directionnels pour différents temps d’intérêt (T), à partir de la pluralité de signaux directionnels, et la détermination (104) d’au moins une caractéristique du milieu en fonction d’une évolution temporelle des signaux inter-directionnels estimés. Figure à publier avec l’abrégé : Fig. 5 PROCEDE D’ANALYSE D’UN MILIEU DOMAINE DE L'INVENTION La présente invention concerne le domaine technique général de l’analyse d’un milieu par propagation d’ondes, et notamment d’ondes sonores ou ultrasonores, ou électromagnétiques. Plus précisément, la présente invention concerne un procédé et un dispositif pour l’analyse d’un objet cible, ou d’un milieu diffus tel qu’un tissu biologique, humain ou animal. Dans la suite, la présente invention sera décrite en référence à l’imagerie médicale par ultrasons, étant entendu que les enseignements décrits ici peuvent être utilisés dans d’autres types d’applications (ultrasons non médicaux, SONAR, RADAR, etc.) utilisant des ondes dont l’amplitude, la fréquence et la phase peuvent être contrôlées ou estimées. ARRIERE PLAN DE L'INVENTION 1. Principe général de la diffusion On connaît différentes solutions d’analyse d’un milieu utilisant le principe de diffusion d’une onde. Afin de simplifier le propos, une description discrète est adoptée ici bien que la réalité de la diffusion soit continue. En référence à la FIG. 1 , la diffusion est un phénomène par lequel une onde U E est déviée et redistribuée dans diverses directions U 1 , U 2 , P, Ech 1 par une interaction avec un « objet diffuseur » Ob 1 contenu dans un milieu. Les solutions d’analyse consistent généralement à : émettre, grâce à un émetteur E, une onde d’excitation UE dans le milieu à analyser, recevoir, grâce à un récepteur R, des échos Ech 1 , Ech 2 suite à la diffusion de l’onde d’excitation incidente UE par un (ou plusieurs) objet(s) diffuseur(s) Ob 1 , Ob 2 contenu(s) dans le milieu et les convertir en signaux de réception, traiter les signaux de réception pour extraire une information sur le milieu. Les signaux de réception sont constitués de deux composantes : une composante dite de « diffusion simple » , et une composante dite de « diffusion multiple » . Plus précisément et comme rappelé précédemment, la diffusion de l’onde induit sa déviation et sa redistribution dans diverses directions, de sorte que certaines parties P de l’onde déviée peuvent interagir avec d’autres objets diffuseurs Ob 2 contenus dans le milieu avant d’être reçues par le récepteur R. Les parties Ech 1 de l’onde déviée qui sont captées par le récepteur R avant d’interagir avec d’autres objets diffuseurs Ob 2 constituent la composante de diffusion simple. Les parties Ech 2 de l’onde déviées qui sont captées par le récepteur R après avoir interagi avec plusieurs objets récepteurs Ob 1 , Ob 2 constituent la composante de diffusion multiple. Ici la description s’est limitée à deux objets diffuseurs mais est bien entendue généralisable à un plus grand nombre de diffuseurs. Les solutions d’analyse traditionnelles d’un milieu utilisent généralement la composante de diffusion simple des signaux de réception : à chaque instant donné de réception d’un écho, la mesure du temps écoulé entre l’émission de l’onde d’excitation et la réception de la composante de diffusion simple associée est représentative de la position de l’objet diffuseur. Toutefois, la composante de diffusion simple étant combinée à la composante de diffusion multiple, ces solutions d’analyse traditionnelles ne sont pas applicables à des milieux incluant un grand nombre d’objets diffuseurs Ob 1 , Ob 2 , ou si les diffuseurs présents dans le milieu ont un pouvoir diffusant élevé. 2. Principe de la rétrodiffusion cohérente dans le domaine de l’optique Ainsi, dans les milieux fortement hétérogènes, comme par exemple les milieux poreux (os, poumon, …), les solutions d’analyse traditionnelles utilisées pour cartographier une région d’intérêt ne sont plus valables. En effet, la relation existant entre : le temps « t » s’écoulant entre l’émission (par un émetteur/récepteur) d’une onde et la réception (par l’émetteur/récepteur) d’une onde réfléchie, et la distance « d » entre l'émetteur/récepteur et le diffuseur n’est plus valable du fait de l’apparition du phénomène de diffusion multiple. Cette diffusion multiple (principal régime de diffusion dans les milieux poreux tels que le poumon) contient néanmoins des informations caractérisant le milieu. En particulier, si la diffusion multiple est suffisamment importante (diffusion forte, nombre élevé de diffuseurs), il est possible de modéliser le parcours d’une onde dans le milieu comme une marche aléatoire obéissant elle-même aux propriétés de la diffusion, avec une constante de diffusion D, elle-même liée au libre parcours moyen de diffusion L. Ce comportement, identifié en optique, donne lieu à un écho contenant une composante de diffusion incohérente (fond chaotique de diffusion multiple) et une composante de diffusion cohérente, concentrée proche du lieu d’émission, qui donne lieu à un phénomène appelé traditionnellement “cône de rétrodiffusion cohérente” . En particulier, il a été démontré en optique que la cohérence des ondes n'était pas totalement perdue grâce à l'effet de "rétrodiffusion cohérente" , même dans un milieu fortement multi-diffusant. En effet, les interférences constructives entre chemins réciproques font, que lorsqu’une moyenne sur plusieurs illuminations du milieu est réalisée, l'intensité rétrodiffusée est deux fois plus importante dans la direction initiale d'illumination du milieu, que dans les autres directions ( FIG. 2 ), ce qui définit un pic (appelé "cône de rétrodiffusion cohérente" , Crc) dont la largeur et l’évolution au cours du temps sont reliées aux propriétés de transport du milieu. 3. Adaptation du principe de rétrodiffusion cohérente dans le domaine des ultrasons Ces principes d’abord observés en optique ont ensuite été utilisés dans le domaine de l’imagerie ultrasonore, notamment dans les documents WO 2020/117789 et WO 2019/104340. Le document WO 2020/117789 décrit un système de détection de lésions et/ou d'anomalies dans un milieu hétérogène en utilisant des ondes ultrasonores. Ce système comprend : un réseau d'éléments transducteurs ultrasonores conçu pour émettre des signaux ultrasonores et pour recevoir des signaux ultrasonores rétrodiffusés, et un processeur apte à exécuter un algorithme de détection de lésion qui traite les signaux ultrasonores rétrodiffusés pour : obtenir une matrice de réponse inter-éléments (IRM), diviser l'IRM en une pluralité de sous-IRM correspondant à des sous-ensembles des éléments transducteurs ultrasonores, et effectuer ensuite un traitement supplémentaire de détection de dépression afin d’identifier une ou plusieurs lésions et/ou anomalies. Le document WO 2019/104340 décrit quant à lui un procédé de mesure de propriétés micro-architecturales d’un réseau vasculaire. Le procédé comprend les étapes consistant à : injecter un agent de contraste ultrasonore dans un milieu, transmettre une pluralité d'impulsions acoustiques dans le milieu à l'aide d'un réseau de transducteurs ultrasonores, chaque impulsion acoustique respective étant transmise à partir d'un ou de plusieurs éléments du réseau de transducteurs ultrasonores, recevoir une pluralité de signaux rétrodiffusés avec le réseau de transducteurs ultrasonores en réponse à chaque impulsion acoustique respective, obtenir une matrice de réponse comprenant les signaux rétrodiffusés extraire une contribution cohérente ou incohérente aux signaux rétrodiffusés à partir de la matrice de réponse, et quantifier une propriété du réseau vasculaire sur la base de la contribution cohérente ou incohérente aux signaux rétrodiffusés. Dans les documents WO 2020/117789 et WO 2019/104340, le réseau de transducteurs ultrasonores est utilisé d'une manière non conventionnelle. Au lieu d'effectuer une formation de faisceau et une imagerie conventionnelle, les impulsions ultrasonores sont transmises en déclenchant des éléments individuels du réseau de transducteurs ultrasonores un par un et en recevant les signaux rétrodiffusés sur chaque élément séparément. Ce faisant, une matrice de réponse inter-éléments est acquise, qui peut être utilisée pour extraire la contribution cohérente ou incohérente aux signaux rétrodiffusés. Par exemple, dans certaines implémentations, la contribution incohérente est extraite. Ainsi, les documents WO 2020/117789 et WO 2019/104340 proposent d’adapter la théorie optique de "rétrodiffusion cohérente" à une approche par point dans le domaine des ultrasons, c’est à dire en émettant une onde ultrasonore à une position donnée (avec un transducteur ultrasonore unique), et en recevant une onde réfléchie à une autre position (avec un autre transducteur ultrasonore unique). Le cône de rétrodiffusion cohérente est alors observé lorsque la position du transducteur ayant émis l’onde ultrasonore (point d’émission) est proche de la position du transducteur ayant reçu l’onde ultrasonore réfléchie par le milieu (point de réception). Toutefois, ce cône de rétrodiffusion est si fin qu’il est la plupart du temps considéré comme n’existant que lorsque les points d’émission et de réception sont confondus. L’étude porte alors sur le grossissement d’un halo de diffusion incohérente, dont l’étalement spatial en fonction du temps est directement lié à la constante de diffusion D. 4. Problématique associée aux procédés existants dans le domaine des ultrasons La technique décrite dans les documents WO 2020/117789 et WO 2019/104340 consistant à : construire un matrice inter-élément à partir d’émission et de réception ultrasonore dans le milieu, et à utiliser cette représentation pour extraire un signal comprenant un cône de rétrodiffusion cohérente et un halo de diffusion incohérente, possède plusieurs limitations qui dégradent fortement les performances de celle-ci en conditions réelles. En particulier : les points d’émission et de réception doivent être en contact direct avec le milieu étudié, et le caractère incohérent des données étudiées implique une forte sensibilité au bruit. 4.1. Nécessité de contact direct avec le milieu étudié La première limitation (points d’émission et de réception en contact avec le milieu) pose problème lorsque le milieu n’est pas directement accessible. Par exemple, lorsque le milieu étudié est un poumon humain, celui-ci n’est pas directement accessible puisqu’il est protégé par la plèvre et par un couche de tissus sous-cutanés. Il n’est donc pas possible de positionner les transducteurs ultrasonores émettant l’onde ultrasonore (point d’émission) et recevant l’onde ultrasonore réfléchie (point de réception) directement en contact avec le poumon à étudier. Pour s’affranchir de cette limitation, une solution peut consister à focaliser l’onde émise par le transducteur directement sur la surface du poumon afin de créer un point d’émission virtuel à la surface du poumon à étudier. Toutefois, cette solution nécessite la connaissance précise de la position du poumon en temps réel, et la capacité d’obtenir une focalisation suffisamment précise malgré les aberrations introduites par la traversée des couches de tissus sous-cutanés. 4.2. Caractère incohérent des données étudiées La deuxième limitation peut également être mitigée en utilisant une focalisation qui concentre l’énergie de plusieurs transducteurs ultrasonores (du réseau de transducteurs) afin de créer un point d’émission virtuel émettant une onde ultrasonore de plus forte intensité. Cependant, cette focalisation n’apporte qu’un gain d’énergie limité, en particulier en champ proche, et doit être surveillée pour éviter d’endommager la surface du milieu qui peut être particulièrement fragile (comme c’est le cas pour la surface du poumon, qui absorbe particulièrement l’énergie ultrasonore et peut être inflammé en cas d’insonification trop puissante). Un but de la présente invention est de proposer un procédé et un dispositif d’analyse d’un milieu permettant de remédier à au moins l’un des inconvénients précités. BREVE DESCRIPTION DE L'INVENTION A cet effet, l’invention propose un procédé d’analyse d’un milieu diffusant, remarquable en ce qu’il comprend les étapes suivantes : l’acquisition, par un système d’acquisition incluant un émetteur et un récepteur, de signaux, en réponse à l’émission d’une onde dans le milieu se propageant selon une direction d’émission souhaitée, le calcul, à partir des signaux acquis, d’une pluralité de signaux directionnels, chaque signal directionnel étant représentatif de l’évolution au cours du temps d’une onde réfléchie par le milieu et se propageant selon une direction choisie de réception, chaque signal directionnel étant associé à un couple respectif de directions d’émission et de réception, l’estimation d’un ensemble de signaux inter-directionnels pour différents temps d’intérêt, à partir de la pluralité de signaux directionnels, chaque signal inter-directionnel étant représentatif de la rétrodiffusion du milieu à un temps d’intérêt respectif, en fonction d’une ouverture angulaire correspondant à un angle entre la direction d’émission et la direction de réception, la détermination d’au moins une caractéristique du milieu en fonction d’une évolution temporelle des signaux inter-directionnels estimés. On entend, dans le cadre de la présente invention, par “ouverture angulaire” , un angle formé entre une direction d’émission d’une onde, et une direction de réception d’une onde, notamment d’une direction de déplacement d’un front d’une onde acoustique. On entend, dans le cadre de la présente invention, par “direction moyenne” , une direction correspondant à une bissectrice entre des directions d’émission et de réceptions d’ondes. Des aspects préférés mais non limitatifs de l’invention sont les suivants : l’étape d’estimation peut comprendre, pour chaque temps d’intérêt considéré, les étapes suivantes : sélectionner un premier sous-groupe de signaux directionnels parmi la pluralité de signaux directionnels, les signaux directionnels du premier sous-groupe étant associés à des couples de directions d’émission et de réception ayant une première direction moyenne identique, la première direction moyenne correspondant à une bissectrice entre les directions d’émission et de réception de chaque couple, sélectionner au moins un deuxième sous-groupe de signaux directionnels parmi la pluralité de signaux directionnels, les signaux directionnels du deuxième sous-groupe étant associés à des couples de directions d’émission et de réception ayant une deuxième direction moyenne identique, la deuxième direction moyenne étant différente de la première direction moyenne, calculer, à partir des signaux directionnels du premier sous-groupe, un premier ensemble de signaux inter-directionnels élémentaires pour différents instants autour du temps d’intérêt, chaque signal élémentaire du premier ensemble étant représentatif de la rétrodiffusion du milieu à la première direction moyenne et à un instant respectif, en fonction de l’ouverture angulaire, calculer, à partir des signaux directionnels du deuxième sous-groupe, un deuxième ensemble de signaux inter-directionnels élémentaires pour différents instants autour du temps d’intérêt, chaque signal élémentaire du deuxième ensemble étant représentatif de la rétrodiffusion du milieu à la deuxième direction moyenne et à un instant respectif, en fonction de l’ouverture angulaire, combiner les signaux inter-directionnels élémentaires des premier et deuxième ensembles afin d’obtenir le signal inter-directionnel pour le temps d’intérêt considéré ; chaque signal inter-directionnel élémentaire peut être un signal complexe obtenu par transformée de Hilbert ou démodulation, la sous-étape de combinaison consistant à : extraire les modules des signaux inter-directionnels élémentaires des premier et deuxième ensembles, et calculer la moyenne desdits modules afin d’obtenir le signal inter-directionnel pour le temps d’intérêt considéré ; chaque signal inter-directionnel élémentaire peut être un signal complexe obtenu par transformée de Hilbert ou démodulation, la sous-étape de combinaison consistant à : aligner les signaux inter-directionnels élémentaires des premier et deuxième ensembles entre eux en multipliant chaque signal inter-directionnel élémentaire par un nombre complexe ayant un argument d’alignement pour obtenir des signaux inter-directionnels élémentaires alignés, calculer la moyenne des signaux inter-directionnels élémentaires alignés afin d’obtenir le signal inter-directionnel pour le temps d’intérêt considéré ; l’argument d’alignement peut être choisi afin de maximiser l’énergie du signal inter-directionnel final pour le temps d’intérêt considéré autour de l’ouverture angulaire nulle ; l’étape d’estimation peut comprendre, pour chaque temps d’intérêt considéré, les étapes suivantes : sélectionner un premier sous-groupe de signaux directionnels parmi la pluralité de signaux directionnels, les signaux directionnels du premier sous-groupe étant associées à des couples de directions d’émission et de réception ayant une première direction d’émission identique, sélectionner au moins un deuxième sous-groupe de signaux directionnels parmi la pluralité de signaux directionnels, les signaux directionnels du deuxième sous-groupe étant associés à des couples de directions d’émission et de réception ayant une deuxième direction d’émission identique, la deuxième direction d’émission étant différente de la première direction d’émission, calculer, à partir des signaux directionnels du premier sous-groupe, un premier ensemble de signaux inter-directionnels élémentaires pour différents instants autour du temps d’intérêt, chaque signal élémentaire étant représentatif de la rétrodiffusion du milieu à la première direction d’émission et à un instant respectif, en fonction de l’ouverture angulaire, calculer, à partir des signaux directionnels du deuxième sous-groupe, un deuxième ensemble de signaux inter-directionnels élémentaires pour différents instants autour du temps d’intérêt, chaque signal élémentaire étant représentatif de la rétrodiffusion du milieu à la deuxième direction d’émission et à un instant respectif, en fonction de l’ouverture angulaire, combiner les signaux inter-directionnels élémentaires des premier et deuxième ensembles afin d’obtenir le signal inter-directionnel pour le temps d’intérêt considéré ; chaque signal inter-directionnel élémentaire est un signal complexe obtenu par transformée de Hilbert ou démodulation, la sous-étape de combinaison consistant à : extraire les modules des signaux inter-directionnels élémentaires des premier et deuxième ensembles, et calculer la moyenne desdits modules afin d’obtenir le signal inter-directionnel pour le temps d’intérêt considéré ; l’étape de détermination d’au moins une caractéristique du milieu comprend les sous-étapes suivantes : pour chaque signal inter-directionnel estimé : extraire une composante gaussienne représentative d’un cône de rétrodiffusion cohérente, calculer une largeur de la composante gaussienne, calculer un facteur de diminution à partir des largeurs calculées pour les différents temps d’intérêt, ledit facteur étant représentatif de la diminution des largeurs calculées pour les différents temps d’intérêt, définir une caractéristique du milieu à partir du facteur de diminution ; l’émetteur et le récepteur peuvent consister en un réseau de transducteur, chaque signal directionnel étant représentatif de l’évolution au cours du temps : d’une onde acoustique réfléchie par le milieu et se propageant selon une direction de réception souhaitée, en réponse à l’émission d’une onde acoustique dans le milieu se propageant selon une direction d’émission souhaitée. D'autres avantages et caractéristiques du procédé selon l’invention ressortiront mieux de la description qui va suivre de plusieurs variantes d’exécution, données à titre d'exemples non limitatifs, à partir des dessins annexés sur lesquels : FIG. 1 est une représentation schématique illustrant le principe de diffusion d’une onde, FIG. 2 est une représentation schématique d’un cône de rétrodiffusion cohérente, FIG. 3 est une représentation schématique d’un dispositif d’imagerie par ultrasons incluant un réseau de transducteurs et une (ou plusieurs) unité(s) de calcul, FIG. 4 est une représentation schématique d’ondes plane et spirale émises par des réseaux de transducteurs plan et courbe, FIG. 5 est une représentation schématique d’un procédé d’analyse selon l’invention, FIG. 6 est une représentation schématique illustrant le principe d’émission d’une onde plane à partir d’un réseau de transducteurs, FIG. 7 , FIG. 7 , FIG. 7 , FIG. 7 sont des représentations schématiques illustrant le principe de détermination d’une matrice inter-directionnelle, FIG. 8 est une représentation schématique de plusieurs matrices inter-directionnelles K déterminées pour différents instants t, FIG. 9 est une représentation schématique illustrant un premier exemple de construction d’un signal inter-directionnel à partir d’une matrice inter-directionnelle, FIG. 10 est une représentation schématique illustrant un deuxième exemple de formation d’un signal inter-directionnel à partir d’une matrice inter-directionnelle, FIG. 11 est une représentation schématique illustrant différents sous-groupes de signaux directionnels sélectionnés dans une matrice inter-directionnelle, FIG. 12 est une représentation schématique de sous-étapes du procédé d’analyse, FIG. 13 , FIG. 14 et FIG. 15 sont des illustrations schématiques de signaux inter-directionnels, FIG. 16 est un diagramme illustrant l’évolution temporelle d’un cône de rétrodiffusion cohérente. DESCRIPTION DETAILLEE DE L'INVENTION On va maintenant décrire différents modes de réalisation du procédé et du dispositif d’analyse d’un milieu selon l’invention en référence aux figures. Dans ces différentes figures, les éléments équivalents sont désignés par la même référence numérique. L’invention sera décrite en référence au domaine de l’imagerie du corps humain par échographie. Il est bien évident pour l’homme du métier que le procédé et le dispositif d’analyse d’un milieu selon l’invention peuvent être utilisés pour d’autres applications, telle que des applications SONAR, RADAR, ou d’autres applications non médicales (sismographie, étude de matériaux tels que des bétons ou des matériaux polycristallin, etc.). 1. Présentation Le procédé selon l’invention permet l’analyse d’un milieu hétérogène à partir d’ondes émises (dans le milieu) et reçues (du milieu) selon des directions d’émission et de réception souhaitées. Dans le domaine de l’imagerie ultrasonore, émettre et recevoir selon des directions peut être réalisé avec le dispositif d’analyse décrit ci-dessous en référence à la FIG. 3 , par exemple en émettant et recevant des ondes planes dans le milieu grâce au dispositif d’analyse décrit ci-dessous. Ce dispositif d’analyse est configuré pour : l’acquisition de données qui sont rassemblées dans un signal inter-directionnel contenant un cône de rétrodiffusion cohérente, et potentiellement un fond incohérent, et le traitement de ce signal inter-directionnel afin de déterminer une (ou plusieurs) caractéristique(s) du milieu. C’est en particulier l’étude du cône de rétrodiffusion cohérente contenu dans ce signal inter-directionnel qui permet de déterminer la (ou les) caractéristique(s) du milieu. En effet, l’évolution temporelle : de l’intensité du cône de rétrodiffusion cohérente, et de l’ouverture angulaire du cône de rétrodiffusion cohérente dépendent de caractéristiques de transport du milieu (comme la constante de diffusion D ou le libre parcours moyen de diffusion). L’utilisation d’ondes plane en émission et en réception répond à la première limitation de la technique décrite dans les documents WO 2020/117789 et WO 2019/104340 puisqu’une onde plane reste plane lorsqu’elle se propage tant qu’elle ne change pas significativement de milieu, et que sa génération peut donc être faite indépendamment de la position de la surface du milieu à étudier. L’utilisation d’onde plane en émission et en réception répond également en partie à la deuxième limitation en réduisant les risques de dégradation de la surface du milieu à analyser. L’étude du cône cohérent (plutôt que l’étude du halo incohérent proposée dans WO 2020/117789 et WO 2019/104340) permet une sensibilité au bruit réduite, et ouvre la possibilité d’un traitement complètement en phase (sans passer à l’intensité) qui augmente significativement le rapport du signal au bruit. 2. Dispositif d’analyse d’un milieu En référence à la FIG. 3 , on a illustré un exemple de dispositif dans lequel le procédé d’analyse d’un milieu décrit dans la suite peut être mis en œuvre. Ce dispositif comprend : un réseau de transducteurs T1-Tn pour l’acquisition de signaux, et une unité de commande et de traitement Uc pour : le pilotage du réseau de transducteurs T1-Tn et le traitement des signaux acquis par le réseau de transducteurs T1-Tn. 2.1. Réseau de transducteurs Le réseau de transducteurs T 1 -T n comprend un jeu de « n » transducteurs ultrasonores ( « n » étant un entier supérieur ou égal à un) disposés linéairement. En variante, les transducteurs T 1 -T n du réseau peuvent être disposés en courbe, ou en cercles concentriques, ou en matrice. Le réseau de transducteurs T 1 -T n permet d’émettre des ondes ultrasonores d’excitation vers un milieu à analyser (organe, tissu biologique, etc.), et de recevoir des échos acoustiques (i.e. ondes ultrasonores réfléchies par le milieu à analyser). Chaque transducteur T 1 -T n consiste par exemple en une plaquette de matériau piézoélectrique de forme rectangulaire revêtue sur ses faces avant et arrière d'électrodes et recouverte sur la face avant de lentilles et de couches d’adaptation d’impédances acoustiques. De tels transducteurs sont connus de l’homme du métier et ne seront pas décrits plus en détails dans la suite. Dans la variante de réalisation illustrée à la FIG. 3 , tous les transducteurs T 1 -T n du réseau sont utilisés à la fois en émission et en réception. Dans d’autres modes de réalisation, des transducteurs distincts peuvent être utilisés pour l'émission et la réception. 2.2. Unité de commande et de traitement L’unité de commande et de traitement U c est reliée au réseau de transducteurs T 1 -T n . Elle permet de piloter les transducteurs T 1 -T n du réseau, et de traiter les données acquises par les transducteurs T 1 -T n du réseau. Plus précisément, l’unité de commande et de traitement U c permet : de commander aux transducteurs T 1 -T n l'émission d'ondes ultrasonores vers le milieu à analyser, de commander aux transducteurs T 1 -T n la réception des échos réfléchis par le milieu à analyser et leur conversion en signaux de réception, de traiter les signaux de réception. L’unité de commande et de traitement U c peut être composée d’une ou de plusieurs entités physiques distinctes, éventuellement distantes du réseau de transducteurs T 1 -T n . L’unité de commande et de traitement U c comprend par exemple : un (ou plusieurs) contrôleur(s) 11, tel qu’un Smartphone, un assistant personnel (ou « PDA » , sigle de l’expression anglo-saxonne « Personal Digital Assistant » ), ou tout type de terminal mobile connu de l’homme du métier ; et un (ou plusieurs) calculateur(s) 12, tel qu’un/des ordinateur(s), un/des micro-ordinateur(s), une/des stations de travail, et/ou d'autres dispositifs connus de l’homme du métier incluant un/des processeur(s), un/des microcontrôleur(s), un/des automate(s) programmable(s), un/des circuit(s) intégré(s) spécifique(s) d'application, et/ou d'autres circuits programmables, une (ou plusieurs) unité(s) de stockage 13 comportant une (ou plusieurs) mémoire(s) qui peu(ven)t être une mémoire ROM/RAM, une clé USB, une mémoire d’un serveur central. Chaque élément de l’unité de commande et de traitement peut être intégré partiellement ou totalement dans la sonde et/ou dans une entité physique distincte de celle-ci (par exemple contrôleur 11 dans la sonde, calculateur 12 et unité de stockage 13 dans une entité distincte). Outre la conservation de données associées à l’analyse d’un milieu, l’unité de stockage 13 permet également de stocker des instructions de code de programmation destinées à exécuter les étapes du procédé d’analyse décrites dans la suite. 3. Procédé d’analyse 3.1. Généralités L’un des aspects avantageux du procédé d’analyse concerne l’utilisation d’ondes directionnelles en émission et en réception, lesdites ondes étant émises et reçues selon des directions d’intérêt souhaitées. Ces ondes directionnelles peuvent être des ondes planes O P telles qu’illustrées à la FIG. 4 . En variante, ces ondes directionnelles peuvent être des ondes spirales O S . On appelle onde plane O P une onde dont le temps d’arrivée du front d’onde à une profondeur z dépend de la position latérale x selon une loi affine t(x, z) = a(z) + bx où a ne dépend que de z et b une constante qui ne dépend ni de z, ni de x. On appelle onde spirale O S l’équivalent d’une onde plane O P en référentiel polaire. Un référentiel polaire est choisi de manière arbitraire, et l’onde spirale O S est définie comme l’onde dont le temps d’arrivée du front d’onde à un rayon r dépend de l’angle selon une loi affine où a ne dépend que de r et b est une constante qui ne dépend ni de r ni de . L’homme du métier remarquera que l’onde plane O P correspond à un cas limite de l’onde spirale où l’origine du repère s’éloigne à l’infini de la sonde. Les ondes directionnelles sont générées en émission ou en réception par le réseau de transducteurs réel T 1 -T n en appliquant à chaque transducteur le retard correspondant au temps d’arrivée du front d’onde à la position dudit transducteur. Il est bien évident pour l’homme du métier que le type d’onde est indépendant de la forme du réseau de transducteurs. Notamment, un réseau de transducteurs plan Tr 1 ou un réseau de transducteurs courbe Tr 2 peut être configuré pour émettre une onde plane (en utilisant une loi de retard adaptée). Il en va de même dans les cas d’ondes directionnelles spirales ou divergentes. Quel que soit le type d’ondes directionnelles (i.e. planes, spirales, divergentes) utilisé en émission et en réception, les directions d’émission et de réception peuvent être caractérisées par des angles d’émission et de réception, ces angles étant définis comme l’angle formé par la direction d’intérêt relativement à la normale au réseau de transducteurs. Plus précisément, on entend, dans le cadre de la présente invention par « angle d’émission » (respectivement « angle de réception » ) α l’angle entre : une direction orthogonale (ci-après dénommée « normale ») à une tangente au réseau de transducteurs T 1 -T n , et une direction normale à la tangente au front d’onde de l’onde émise (respectivement reçue) grâce au réseau de transducteurs T 1 -T n , et ce, au moment d’émission (respectivement de réception) de l’onde directionnelle. 3.2. Etapes du procédé En référence à la FIG. 5 , le procédé d’analyse comprend les étapes suivantes : acquisition 101 d’une pluralité de signaux échos en réponse à l’émission d’une onde dans le milieu se propageant, calcul 102 d’une pluralité de signaux directionnels à partir des signaux échos acquis, estimation 103 d’un ensemble de signaux inter-directionnels pour différents temps d’intérêt « T » , à partir de la pluralité de signaux directionnels, détermination 104 d’au moins une caractéristique du milieu en fonction d’une évolution temporelle des signaux inter-directionnels estimés. Chaque signal inter-directionnel est représentatif de la rétrodiffusion du milieu en fonction d’une ouverture angulaire correspondant à un angle entre une direction d’émission et une direction de réception. 3.3. Etape d’acquisition de la pluralité de signaux directionnels Chaque signal directionnel est représentatif de l’évolution au cours du temps d’une onde réfléchie se propageant selon une direction de réception souhaitée en réponse à l’émission d’une onde incidente se propageant selon une direction d’émission souhaitée. Chaque signal directionnel est donc associé à un couple de directions d’émission et de réception respectif différent des couples de directions des autres signaux directionnels. Ainsi, pour acquérir la pluralité de signaux directionnels : des ondes directionnelles d’émission – se propageant selon différentes directions d’émission – sont émises dans le milieu à étudier grâce au réseau de transducteurs T 1 -T n , et des ondes directionnelles de réception réfléchies par le milieu à étudier sont synthétisées à partir des signaux reçus par le réseau de transducteurs T 1 -T n . Dans la suite, on décrira l’invention en référence à l’utilisation d’ondes directionnelles planes, étant entendu que le procédé peut utiliser d’autres types d’ondes directionnelles telles que des ondes spirales ou divergentes. 3.3.1. Génération d’ondes planes d’émission Pour générer une onde plane d’émission, les transducteurs T 1 -T n du réseau sont activés (en émission) de sorte que les ondes élémentaires El 1 -El n générées par chacun des transducteurs T 1 -T n se combinent pour former une onde plane d’émission ayant une direction d’émission souhaitée. Plus précisément, en fonction des instants d’activation des transducteurs et de l’amplitude des tensions d’excitation appliquées aux transducteurs T 1 -T n par l’unité de commande et de traitement U c , il est possible de contrôler les transducteurs T 1 -T n pour qu’ils produisent des ondes ultrasonores élémentaires El 1 -El n se combinant pour former une onde ultrasonore plane d’émission 14 qui se propage à travers le milieu à analyser selon une direction 15 souhaitée (cf. FIG. 6 ). Cette onde ultrasonore résultante 14 plane peut être émise selon différentes directions d’émission (i.e. différentes directions) en faisant varier les instants d’activation (t, t+Δt, t+2Δt, … t+nΔt) de chaque transducteur T 1 -T n du réseau. Par exemple, pour la génération d’une onde plane d’émission, tous les transducteurs T 1 -T n peuvent être activés : simultanément pour obtenir une onde plane d’émission se propageant selon une direction d’émission orthogonale au réseau de transducteurs T 1 -T n , ou successivement (en fonction d’une loi de retard d’activation) pour obtenir une onde plane d’émission se propageant selon une direction d’émission comprise entre -90° et 90° relativement au réseau de transducteurs T 1 -T n . Dans tous les cas, les transducteurs T 1 -T n sont tous activés pour générer l’onde plane, c’est-à-dire que les transducteurs T 1 -T n sont tous activés en émission pour chaque réception. Le lecteur appréciera que les ondes planes sont générées par une émission d’une durée finie. Cette émission, qui peut être représentée comme une fonction du temps ou de la fréquence, est caractérisée par une fréquence centrale ainsi qu’une largeur de bande. 3.3.2. Synthèse d’ondes planes de réception Dans la suite, on appelle “ondes planes de réceptions“ les signaux synthétisés à partir des signaux reçus par le réseau de transducteurs T 1 -T n (retard et sommation) représentatifs des ondes réfléchies par le milieu et se propageant selon une direction choisie. Pour synthétiser des ondes planes de réception, une solution peut consister à acquérir simultanément des signaux avec les transducteurs T 1 -T n . Ainsi, après l’émission d’une onde plane d’émission présentant une direction d’émission donnée, chaque transducteur est activé en réception pour enregistrer (acquérir et mémoriser) des signaux captés représentatifs de la réverbération par le milieu de l’onde plane d’émission. Pour chaque transducteur T 1 -T n , un signal capté fonction du temps t, {si(t)}0≤i≤n-1 est enregistré. Les signaux captés par les transducteurs T 1 -T n sont ensuite sommés selon une loi de retard temporelle dépendant de l’angle de réception désiré pour l’onde plane de réception. Par exemple pour recevoir une onde plane de réception ayant un angle de réception α à partir des signaux captés {si(t)} 0≤i≤n-1 mesurés par les transducteurs T 1 -T n , l’opération de sommation suivante est effectuée : Où : « c » représente la célérité de l’onde dans le milieu, « x i » représente la position latérale du transducteur i. Le traitement du bloc de signaux captés permet de « réorienter » les réponses enregistrées par les différents transducteurs T 1 -T n pour obtenir les ondes de réception aux différents angles de réception souhaités. On génère ainsi N ondes planes en réception à partir d’une seule émission. 3.3.3. Obtention des signaux directionnels L’étape d’acquisition permet l’obtention d’une pluralité de signaux directionnels, chaque signal directionnel étant associé à un couple respectif de directions d’émission et de réception (différent des couples des autres signaux directionnels de la pluralité de signaux directionnels). Le lecteur appréciera que chaque signal directionnel est un signal réel dont l’amplitude varie en fonction du temps. Ce signal est appelé signal « RF » (pour « Radio-Frequency » ). Afin d’étudier ce signal avec plus de facilité, sa représentation complexe peut être utilisée. Cette représentation, largement connue dans les domaines de la physique des ondes et du traitement du signal, permet la simplification des calculs et de la résolution de l’équation d’onde. Chaque signal directionnel peut être retrouvé en prenant la partie réelle de sa représentation complexe. De plus, le signal complexe peut être caractérisé par son module, représentatif de l’amplitude ou l’enveloppe de l’onde de réception, et par son argument, représentatif de la phase de l’onde de réception. Pour obtenir ce signal complexe à partir du signal RF, plusieurs solutions sont possibles, comme par exemple l’extraction du signal analytique par transformée de Hilbert, ou la démodulation autour de la fréquence centrale suivie par un filtrage passe-bas. On obtient alors un signal IQ (in-phase quadrature), qui peut être modulé ou démodulé. Cette étape peut être réalisée de manière analogique ou digitale. 3.4. Etape d’estimation de l’ensemble de signaux inter-directionnels Une fois les signaux directionnels acquis, ceux-ci sont utilisés pour estimer l’ensemble de signaux inter-directionnels pour différents temps d’intérêt « T » . Pour illustrer la mise en œuvre de cette étape d’estimation, des matrices inter-directionnelles seront utilisées dans les dessins, étant entendu pour l’homme du métier que la construction de telles matrices inter-directionnelles n’est pas nécessaire à l’estimation de l’ensemble de signaux inter-directionnels. 3.4.1. Détermination d’une matrice inter-directionnelle fenêtrée En référence à la FIG. 7 , le principe de construction d’un coefficient K i,j d’une ligne « i » d’une matrice inter-directionnelle K(t) fenêtrée temporellement comprend les étapes suivantes : pour les transducteurs T 1 -T n du réseau, émettre une onde plane d’émission O inc,i ayant une direction d’émission souhaitée (chaque transducteur étant activé selon une loi de retard d’activation prédéterminée pour générer l’onde plane O inc,i ayant la direction d’émission souhaitée), pour les transducteurs T 1 -T n du réseau, recevoir une onde plane de réception O rec,j ayant une direction de réception souhaitée (obtenue en mettant en œuvre soit la première solution, soit la deuxième solution). Le signal de réception obtenu est variable au cours du temps. Ce signal est tronqué (fenêtré) en fenêtres de temps successives (potentiellement partiellement superposées) – de durée Δt – associées chacune à une matrice inter-directionnelle fenêtrée respective. Pour chaque fenêtre de temps, le signal de réception tronqué correspondant est multiplié par une fonction de fenêtrage et enregistré dans le coefficient K i,j (t) de la matrice inter-directionnelle fenêtrée associée K(t). Comme illustré sur les FIG. 7 et 7c, la ligne « i » de chaque matrice inter- directionnelle fenêtrée K(t) est complétée en recevant des ondes planes de réception O rec2 , O rec3 ayant chacune une direction de réception distincte des autres ondes planes de réception et en tronquant les signaux de réception représentatifs des ondes planes de réception O rec2 , O rec3 en fenêtres de temps successives. Ces étapes sont répétées pour chaque ligne de la matrice K(t). Plus précisément, une fois que toutes les ondes planes de réception ont été synthétisées et tronquées pour l’ensemble des directions de réception souhaitées, la ligne « i » de chaque matrice inter- directionnelle fenêtrée K(t) est déterminée. En référence à la FIG. 7 , la construction de la ligne « i+1 » de chaque matrice inter-angles fenêtrée K(t) peut être initiée : en générant une nouvelle onde plane d’émission O inc,i+1 ayant une direction d’émission distincte de la direction d’émission associée à l’onde plane d’émission O inc,i utilisée pour la ligne « i » précédente, en synthétisant des ondes planes de réception O rec1 , etc. ayant les directions de réception souhaitées distinctes, et en tronquant les signaux de réception représentatifs des ondes planes de réception O rec1 , etc. en fenêtres de temps successives. Ainsi dans chaque matrice inter- directionnelle fenêtrée K(t) : chaque ligne « i » est représentative de la direction d’émission d’une onde plane émise, et chaque colonne « j » est représentative d’une direction de réception d’une onde plane reçue. Dans le cadre de la présente invention, les matrices inter-directionnelles fenêtrées K(t) considérées peuvent être acquises avec des directions d’émission et de réception espacées d’un pas constant. L’homme du métier appréciera que la condition de pas constant n’est pas essentielle dans le cadre du procédé d’analyse selon l’invention. Comme illustré à la FIG. 8 , chaque signal directionnel S d (t) correspond alors à un même coefficient K i,j (t) des matrices inter-directionnelles fenêtrées K(t). Ce signal directionnel S d (t) varie en fonction du temps. 3.4.2. Extraction de chaque signal inter-directionnel Pour chaque temps d’intérêt « T » considéré, on construit ensuite un signal inter-directionnel Sid à partir des valeurs au temps « T » des signaux directionnels S d . Différentes méthodes peuvent être mises en œuvre pour l’estimation d’un signal inter-directionnel pour un temps d’intérêt « T » considéré. Notamment, l’estimation de chaque signal inter-directionnel S id peut être réalisée en sélectionnant des signaux directionnels ayant une même direction moyenne et des ouvertures angulaires différentes. Par exemple, la FIG. 9 illustre la formation d’un signal inter-directionnel S id à partir des valeurs des signaux directionnels sélectionnés (signaux directionnels ayant une direction moyenne constante et des ouvertures angulaires différentes) dans une matrice inter-directionnelle K(T) obtenue pour le temps d’intérêt « T » considéré. En variante, l’estimation de chaque signal inter-directionnel S id peut être réalisée en sélectionnant des signaux directionnels ayant une même direction d’émission et des ouvertures angulaires différentes. Par exemple, la FIG. 10 illustre la formation d’un signal inter-directionnel S id à partir des valeurs des signaux directionnels sélectionnés (signaux directionnels ayant une direction d’émission identique et des ouvertures angulaires différentes) dans une matrice inter-directionnelle K(T) obtenue pour le temps d’intérêt « T » considéré. Avantageusement, chaque signal inter-directionnel (associé à un temps d’intérêt « T » considéré) peut être rendu plus robuste en le combinant à d’autres signaux inter-directionnels obtenus : à différents instants « t » (i.e. différentes profondeurs proches de celles d’intérêt) et à différentes directions moyennes (cf. FIG. 11 ) ou à différentes directions d’émission. Par exemple et en référence à la FIG. 12 , lorsque chaque signal inter-directionnel est rendu plus robuste en travaillant sur les directions moyennes, l’étape d’estimation peut comprendre les sous-étapes suivantes : sélectionner 1011 un premier sous-groupe de signaux directionnels parmi la pluralité de signaux directionnels, les signaux directionnels du premier sous-groupe étant associés à des couples de directions d’émission et de réception ayant une première direction moyenne identique (ou respectivement ayant une première direction d’émission identique), calculer 1013, à partir des signaux directionnels du premier sous-groupe, un premier ensemble de signaux inter-directionnels élémentaires (S E1 … S En ) pour différents instants « t » autour du temps d’intérêt « T » , chaque signal inter-directionnel élémentaire du premier ensemble étant représentatif de la rétrodiffusion du milieu à la première direction moyenne (ou respectivement à la première direction d’émission) et à un instant « t » respectif, en fonction de l’ouverture angulaire, sélectionner 1012 un deuxième sous-groupe de signaux directionnels parmi la pluralité de signaux directionnels, les signaux directionnels du deuxième sous-groupe étant associés à des couples de directions d’émission et de réception ayant une deuxième direction moyenne identique mais différente de la première direction moyenne (ou respectivement ayant une deuxième direction d’émission identique, la deuxième direction d’émission étant différente de la première direction d’émission), calculer 1014, à partir des signaux directionnels du deuxième sous-groupe, un deuxième ensemble de signaux inter-directionnels élémentaires (S E1’ … S En ’ ) pour différents instants « t » autour du temps d’intérêt « T » , chaque signal inter-directionnel élémentaire du deuxième (respectivement troisième, quatrième, cinquième …etc.) ensemble étant représentatif de la rétrodiffusion du milieu à la deuxième direction moyenne (ou respectivement à la deuxième direction d’émission) et à un instant « t » respectif, en fonction de l’ouverture angulaire, combiner 1015 les différents signaux élémentaires sélectionnés des premier et deuxième ensembles afin d’obtenir le signal inter-directionnel pour le temps d’intérêt « T » considéré. Comme il sera décrit plus en détail dans la suite, l’étape de combinaison peut inclure une combinaison dans le temps des signaux élémentaires du premier sous-groupe d’une part, et du deuxième sous-groupe d’autre part, ainsi qu’une combinaison des signaux élémentaires des premier et deuxième sous-groupes. Bien entendu, plus de deux sous-groupes de signaux directionnels peuvent être utilisés pour l’obtention du signal inter-directionnel pour le temps d’intérêt « T » considéré. Par exemple, et comme illustré à la FIG. 11 , des troisième P 3 , quatrième P 4 et cinquième P 5 (etc.) sous-groupes peuvent être sélectionnés. Ainsi, la détermination du signal inter-directionnel (associé à un temps d’intérêt « T » considéré) par combinaison des signaux inter-directionnels élémentaires calculés pour les différents sous-groupes de signaux directionnels permet d’améliorer la qualité de celui-ci. Cette combinaison peut être effectuée de deux façons différentes : moyennage en enveloppe (intensité) ; ou moyennage en phase. 3.4.2.1. Moyennage en enveloppe Dans le cas d’un moyennage en enveloppe, les enveloppes des signaux inter-directionnels sont moyennées à travers le temps et les directions moyennes (ou les directions d’émission). Le signal inter-directionnel obtenu (pour le temps d’intérêt « T » considéré) présente la forme illustrée à la FIG. 13 . Il comprend une composante gaussienne représentative du cône de rétrodiffusion cohérente, comme décrit plus en détails au point 3.5 ci-dessous). 3.4.2.2. Moyennage en phase Dans le cas d’un moyennage en phase, les signaux inter-directionnels élémentaires sont combinés dans leur forme complexe. Chaque signal élémentaire complexe peut-être obtenu par transformée de Hilbert ou par démodulation. Chaque signal élémentaire complexe comprend une « partie réelle » et une « partie imaginaire » : l’argument de chaque signal élémentaire complexe correspondant à la phase du signal élémentaire et le module de chaque signal élémentaire complexe correspondant à l’amplitude de chaque signal élémentaire (cette composante est utilisée pour le moyennage en enveloppe décrit en Section 3.4.2.1). En particulier, les signaux complexes modulés ou démodulés (IQ) sont utilisés. Ainsi, on exploite la cohérence d’une partie des signaux de diffusion multiples, qui sont en phase pour un instant donné et une direction moyenne donnée. Plus précisément, chaque signal inter-directionnel élémentaire complexe, lorsqu’il est extrait à direction moyenne constante, présente un cône de rétrodiffusion cohérente dont la phase dépend du temps « t » et de la direction moyenne mais pas de l’ouverture angulaire. Le reste des composantes de ce signal élémentaire (partie incohérente et bruit) a une phase qui peut être considérée comme aléatoire en fonction du temps, de la direction moyenne et de l’ouverture angulaire. Pour faire ressortir le cône de rétrodiffusion cohérente, il est donc possible de moyenner les signaux inter-directionnels élémentaires avec leurs phases. Néanmoins, la phase du cône de rétrodiffusion cohérente dépendant de l’angle moyen et du temps, il convient d’aligner les phases des cônes de rétrodiffusion cohérente des différents signaux élémentaires avant la mise en œuvre de ce moyennage. Ce moyennage a pour but de maximiser l’énergie du cône cohérent obtenu (i.e. l’énergie aux ouvertures angulaires faibles), et est réalisé en multipliant chaque signal élémentaire par une constante complexe ayant une phase dite d'alignement. Cette phase d’alignement est calculée par exemple : en résolvant le problème de maximisation de l’énergie du cône cohérent après moyennage, ou bien en mesurant (ou estimant) la phase de chaque signal élémentaire autour de l’ouverture angulaire nulle et en utilisant l’opposé de cette phase comme phase d’alignement. Le signal inter-directionnel obtenu présente la forme illustrée à la FIG. 14 . 3.4.2.3. En pratique Pour chaque temps d’intérêt « T » considéré, le signal inter-directionnel final est en pratique immédiatement déduit des signaux directionnels, sans passer par la construction de matrices inter-directionnelles ou de groupes de signaux élémentaires. On note ces signaux directionnels , avec : α l’angle de la direction d'émission (angle entre la normale au réseau de transducteur et la normale au front d’onde de l’onde émise), et β l’angle de la direction de réception (angle entre la normale au réseau de transducteur et la normale au front d’onde de l’onde reçue). Dans ce cas, l’angle de la direction moyenne est défini comme étant égale à , et la demi ouverture angulaire est définie comme étant égale à . A chaque paire de directions d'émission et de réception correspond de manière unique une paire de direction moyenne et d’ouverture angulaire. Les signaux directionnels peuvent être exprimés comme suit : Le signal inter-directionnel est obtenu en moyennant les signaux directionnels en fonction de l’angle moyen et du temps t (autour du temps d'intérêt T) : où : M est un opérateur de moyennage, est un ensemble de directions moyennes à disposition, et est une durée sur laquelle le signal inter-directionnel est moyenné. L’opérateur de moyennage M peut être défini de deux façons différentes. Notamment, l’opérateur de moyennage M peut être défini comme un moyennage en enveloppe ; dans ce cas, le signal inter-directionnel peut être écrit comme suit : où est l’opérateur de moyennage arithmétique sur la variable x. Avec un tel opérateur de moyennage M, un signal inter-directionnel de forme gaussienne tel qu’illustré à la FIG. 13 est obtenu. En variante, l’opérateur de moyennage M peut être défini comme un moyennage en phase avec réalignement. Dans ce cas, on commence par forcer la phase de tous les cônes de rétrodiffusion cohérente à 0 en multipliant par un signal d’alignement, typiquement le conjugué du signal correspondant au même instant, au même angle moyen et 0 une ouverture angulaire nulle. On effectue ensuite un moyennage avant de passer à l’intensité. Dans ce cas, le signal inter-directionnel vaut : où est le conjugué. Alternativement, on peut écrire : où est l’opérateur de phase. Les deux expressions ci-dessus permettent chacune d’obtenir un signal inter-directionnel respectif pour un temps d’intérêt « T » considéré : la première expression (expression 1) permet d’obtenir un premier signal inter-directionnel, la deuxième expression (expression 2) permet d’obtenir un deuxième signal inter-directionnel. Ces premier et deuxième signaux inter-directionnels sont différents, puisqu’ils présentent des expressions théoriques légèrement différentes (le premier signal inter-directionnel est égal au deuxième signal inter-directionnel multiplié par une constante), mais dans les deux cas, les phases sont réalignées avant le moyennage. Quelle que soit l’expression utilisée (expression 1 ou expression 2), on obtient un signal inter-directionnel de la forme illustrée à la FIG. 14 . La différence observée entre les deux types de moyennages (moyennage en enveloppe tel qu’illustré à la FIG. 13 / moyennage en phase tel qu’illustré à la FIG. 14 ) vient de l'interférence destructive du fond incohérent (constant en fonction de l’ouverture angulaire) lors du moyennage en phase, ce qui ne laisse que le cône de rétrodiffusion cohérente. 3.5. Détermination d’une (ou plusieurs) caractéristiques à partir des signaux inter-directionnels Pour chaque temps d’intérêt « T » considéré, le signal inter-directionnel obtenu peut être décomposé en deux parties : une partie incohérente, qui correspond à un fond chaotique de diffusion multiple ; cette partie incohérente est présente dès le début sur toutes les ouvertures angulaires de manière homogène, comme illustré à la FIG. 15 , une partie cohérente, appelée “cône de rétrodiffusion cohérente” , qui correspond à des signaux encore en phase qui interfère de manière constructive entre eux ; cette partie cohérente est concentrée autour des faibles ouvertures angulaires. Pour le signal inter-directionnel obtenu après moyennage en enveloppe, les deux parties qui le composent peuvent être identifiées facilement, par exemple en utilisant la différence entre l’ouverture angulaire nulle et les grandes ouvertures angulaires ou en calculant la fonction de la forme f(x) = 1 + exp(-ax²) la plus proche des données. Pour le signal inter-directionnel obtenu après moyennage en phase, le moyennage a pour effet de faire interférer constructivement les cônes de rétrodiffusion cohérente et destructivement la partie incohérente du signal. Il ne reste donc plus que la partie cohérente. Dans tous les cas, le signal inter-directionnel obtenu dépend fortement du temps. D’abord, son intensité diminue, théoriquement en t (-3/2) (si l’on se rend indépendant des effets de diffusion, de directivité et de masquage spatial). Cette tendance théorique dépend de l’hypothèse d’une diffusion forte dans le milieu, avec un nombre de diffuseurs suffisamment élevé pour engendrer une marche gaussienne de l’onde acoustique, et ce de manière homogène. L’étude de l’intensité du signal inter-directionnel (partie cohérente ou incohérente, les deux étant égales), permet donc de valider les hypothèses faites sur le milieu, ainsi que de localiser les zones à forte diffusion. De plus, l’extension angulaire de la partie cohérente du signal inter-directionnel diminue dans le temps. La théorie donne que la partie cohérente du signal inter-directionnel se comporte comme une gaussienne en fonction de l’ouverture angulaire, gaussienne dont la variance évolue avec le temps (1/t²). L’étude de l’évolution temporelle de l’extension angulaire du cône de rétrodiffusion cohérente permet de remonter à des paramètres physiques du milieu, comme la constante de diffusion D ou le libre parcours de diffusion moyen L (proportionnel à la constante de diffusion). Ces derniers sont liés à la pente de l’extension angulaire en fonction du temps de la composante gaussienne du signal inter-directionnel (partie cohérente). En effet, la théorie nous dit : pour le moyennage en enveloppe (intensité), et pour le moyennage en phase. Si on choisit la gaussienne qui correspond le mieux au cône cohérent et qu’on reporte l’inverse de sa variance angulaire au cours du temps, on obtient en théorie pour le moyennage en phase : . En pratique, on vérifie grâce à des simulations qu’on obtient effectivement l’inverse de la déviation standard linéaire en fonction du temps, comme illustré à la FIG. 16 . 4. Conclusions Le procédé décrit ci-dessus permet une amélioration du rapport signal sur bruit des signaux inter-directionnels calculés pour différents temps d’intérêt « T » . L’utilisation d’ondes directionnelles permet de s’affranchir de la nécessité d’émettre et de recevoir des ondes à la surface du milieu à étudier. Par ailleurs, l’étude de l’évolution temporelle du cône de rétrodiffusion cohérente (qui s’amincit avec la profondeur) plutôt que du halo incohérent (tel que proposé dans les documents WO 2020/117789 et WO 2019/104340) permet de caractériser le milieu de manière plus robuste. Dans le cas d’un moyennage en phase des signaux inter-directionnel déterminés pour différents temps d’intérêt « T » , l’interférence constructive des cônes de diffusion cohérente permet un gain de SNR supplémentaire ainsi qu’une facilité à estimer son extension angulaire. Dans la description qui précède, l’invention était décrite dans le cas de l’utilisation d’onde planes. Le lecteur appréciera que l’invention peut être mise en œuvre en utilisant d’autres types d’ondes fortement directives, tel que des ondes spirales. De même dans la description qui précède, l’invention était décrite en référence à un réseau de transducteurs T 1 -T n ayant une géométrie linéaire. Il est bien évident pour l’homme du métier que le réseau de transducteurs T 1 -T n peut présenter d’autres formes telles qu’une forme courbe, ou matricielle. Dans le cas d’une sonde matricielle, donc bidimensionnelle, la méthode précédente se généralise en définissant des ondes planes ou spirales bidimensionnelles. Ces ondes planes ou spirales ne sont rien d’autres que la combinaison d’une onde plane ou spirale selon un axe avec une onde plane ou spirale selon l’autre axe, donnant ainsi des lois de retard définies dans des repères cartésiens, cylindriques ou polaires.   Procédé d’analyse d’un milieu diffusant, caractérisé en ce que le procédé comprend les étapes suivantes : l’acquisition (101), par un système d’acquisition incluant un émetteur et un récepteur, de signaux, en réponse à l’émission d’une onde dans le milieu se propageant selon une direction d’émission souhaitée, le calcul (102), à partir des signaux acquis, d’une pluralité de signaux directionnels, chaque signal directionnel étant représentatif de l’évolution au cours du temps d’une onde réfléchie par le milieu et se propageant selon une direction choisie de réception, chaque signal directionnel étant associé à un couple respectif de directions d’émission et de réception, l’estimation (103) d’un ensemble de signaux inter-directionnels pour différents temps d’intérêt (T), à partir de la pluralité de signaux directionnels, chaque signal inter-directionnel étant représentatif de la rétrodiffusion du milieu à un temps d’intérêt (T) respectif, en fonction d’une ouverture angulaire correspondant à un angle entre la direction d’émission et la direction de réception, la détermination (104) d’au moins une caractéristique du milieu en fonction d’une évolution temporelle des signaux inter-directionnels estimés. Procédé selon la revendication 1, dans lequel l’étape d’estimation comprend, pour chaque temps d’intérêt (T) considéré, les étapes suivantes : sélectionner (1011) un premier sous-groupe (P1) de signaux directionnels parmi la pluralité de signaux directionnels, les signaux directionnels du premier sous-groupe (P1) étant associés à des couples de directions d’émission et de réception ayant une première direction moyenne identique, la première direction moyenne correspondant à une bissectrice entre les directions d’émission et de réception de chaque couple, sélectionner (1012) au moins un deuxième sous-groupe (P2) de signaux directionnels parmi la pluralité de signaux directionnels, les signaux directionnels du deuxième sous-groupe étant associés à des couples de directions d’émission et de réception ayant une deuxième direction moyenne identique, la deuxième direction moyenne étant différente de la première direction moyenne, calculer (1013), à partir des signaux directionnels du premier sous-groupe (P1), un premier ensemble de signaux inter-directionnels élémentaires (SE1, SEn) pour différents instants (t) autour du temps d’intérêt (T), chaque signal élémentaire du premier ensemble étant représentatif de la rétrodiffusion du milieu à la première direction moyenne et à un instant (t) respectif, en fonction de l’ouverture angulaire, calculer (1014), à partir des signaux directionnels du deuxième sous-groupe (P2), un deuxième ensemble de signaux inter-directionnels élémentaires (SE1’, SEn’) pour différents instants (t) autour du temps d’intérêt (T), chaque signal élémentaire du deuxième ensemble étant représentatif de la rétrodiffusion du milieu à la deuxième direction moyenne et à un instant (t) respectif, en fonction de l’ouverture angulaire, combiner (1015) les signaux inter-directionnels élémentaires des premier et deuxième ensembles afin d’obtenir le signal inter-directionnel pour le temps d’intérêt (T) considéré. Procédé selon la revendication 2, dans lequel chaque signal inter-directionnel élémentaire est un signal complexe obtenu par transformée de Hilbert ou démodulation, la sous-étape de combinaison consistant à : extraire les modules des signaux inter-directionnels élémentaires (SE1-SEn, (SE1’-SEn’) des premier et deuxième ensembles, et calculer la moyenne desdits modules afin d’obtenir le signal inter-directionnel pour le temps d’intérêt (T) considéré. Procédé selon la revendication 2, dans lequel chaque signal inter-directionnel élémentaire est un signal complexe obtenu par transformée de Hilbert ou démodulation, la sous-étape de combinaison consistant à : aligner les signaux inter-directionnels élémentaires des premier et deuxième ensembles entre eux en multipliant chaque signal inter-directionnel élémentaire par un nombre complexe ayant un argument d’alignement pour obtenir des signaux inter-directionnels élémentaires alignés, calculer la moyenne des signaux inter-directionnels élémentaires alignés afin d’obtenir le signal inter-directionnel pour le temps d’intérêt (T) considéré. Procédé selon la revendication 4, dans lequel l’argument d’alignement est choisi afin de maximiser l’énergie du signal inter-directionnel final pour le temps d’intérêt (T) considéré autour de l’ouverture angulaire nulle. Procédé selon la revendication 1, dans l equel l’étape d’estimation comprend, pour chaque temps d’intérêt (T) considéré, les étapes suivantes : sélectionner un premier sous-groupe de signaux directionnels parmi la pluralité de signaux directionnels, les signaux directionnels du premier sous-groupe étant associées à des couples de directions d’émission et de réception ayant une première direction d’émission identique, sélectionner au moins un deuxième sous-groupe de signaux directionnels parmi la pluralité de signaux directionnels, les signaux directionnels du deuxième sous-groupe étant associés à des couples de directions d’émission et de réception ayant une deuxième direction d’émission identique, la deuxième direction d’émission étant différente de la première direction d’émission, calculer, à partir des signaux directionnels du premier sous-groupe, un premier ensemble de signaux inter-directionnels élémentaires pour différents instants (t) autour du temps d’intérêt (T), chaque signal élémentaire étant représentatif de la rétrodiffusion du milieu à la première direction d’émission et à un instant (t) respectif, en fonction de l’ouverture angulaire, calculer, à partir des signaux directionnels du deuxième sous-groupe, un deuxième ensemble de signaux inter-directionnels élémentaires pour différents instants (t) autour du temps d’intérêt (T), chaque signal élémentaire étant représentatif de la rétrodiffusion du milieu à la deuxième direction d’émission et à un instant (t) respectif, en fonction de l’ouverture angulaire, combiner les signaux inter-directionnels élémentaires des premier et deuxième ensembles afin d’obtenir le signal inter-directionnel pour le temps d’intérêt (T) considéré. Procédé selon la revendication 6, dans lequel chaque signal inter-directionnel élémentaire est un signal complexe obtenu par transformée de Hilbert ou démodulation, la sous-étape de combinaison consistant à : extraire les modules des signaux inter-directionnels élémentaires des premier et deuxième ensembles, et calculer la moyenne desdits modules afin d’obtenir le signal inter-directionnel pour le temps d’intérêt (T) considéré. Procédé selon l’une quelconque des revendications 1 à 7, dans lequel l’étape de détermination d’au moins une caractéristique du milieu comprend les sous-étapes suivantes : pour chaque signal inter-directionnel estimé : extraire une composante gaussienne représentative d’un cône de rétrodiffusion cohérente, calculer une largeur de la composante gaussienne, calculer un facteur de diminution à partir des largeurs calculées pour les différents temps d’intérêt (T), ledit facteur étant représentatif de la diminution des largeurs calculées pour les différents temps d’intérêt (T), définir une caractéristique du milieu à partir du facteur de diminutione. Procédé selon l’une quelconque des revendications 1 à 8, dans lequel l’émetteur et le récepteur consistent en un réseau de transducteur, chaque signal directionnel étant représentatif de l’évolution au cours du temps : d’une onde acoustique réfléchie par le milieu et se propageant selon une direction de réception souhaitée, en réponse à l’émission d’une onde acoustique dans le milieu se propageant selon une direction d’émission souhaitée.