FR 3084166 A1 20200124 FR 1856720 A 20180719 PROCÉDÉS ET SYSTÈMES DE CARACTÉRISATION ULTRASONORE NON INVASIVE D'UN MILIEU HÉTÉROGÈNE DOMAINE TECHNIQUE [001] La présente description concerne des procédés et systèmes de caractérisation ultrasonore non invasive d'un milieu hétérogène, et s'applique notamment à l'imagerie médicale ou au contrôle non destructif et plus généralement à tous les domaines dans lesquels l'imagerie ultrasonore peut être utilisée. ETAT DE L'ART [002] Dans le domaine de l'imagerie acoustique, on cherche à caractériser un milieu inconnu en le sondant de manière active avec des ondes ultrasonores. C'est notamment le principe de l'échographe en imagerie médicale. [3] La résolution d'un système d'imagerie acoustique peut être définie comme la capacité à discerner les petits détails d'un objet. En principe, un système d'imagerie acoustique est limité par la diffraction et la résolution théorique est donnée par X. / 2 (où X. est la longueur d'onde du son dans le milieu), ou par l'ouverture angulaire finie du détecteur. En pratique, la résolution est cependant souvent détériorée par les variations de vitesses du son dans le milieu de propagation. [4] En effet, la plupart du temps en imagerie acoustique, le milieu est considéré comme homogène, avec une vitesse du son co constante. Or l'hypothèse d'un milieu homogène n'est pas toujours respectée. Par exemple, dans le cas de l'échographie du foie, la sonde est placée entre les côtes du patient. Les ondes acoustiques traversent une succession de couches de graisse et de muscle avant d'atteindre l'organe visé. Les tissus mous possèdent des propriétés mécaniques différentes. La vitesse du son est donc loin d'être homogène, entre 1450 m/s pour les tissus adipeux et 1600 m/s pour le foie. Les variations de vitesse du son entraînent un déphasage différent des ondes suivant les endroits par lesquels elles se propagent. Il en résulte une aberration du front d'onde acoustique qui conduit à une distorsion de l'image échographique résultante, et donc à une dégradation de sa résolution et de son contraste. Ces aberrations peuvent être telles qu'elles compromettent l'examen médical. [5] Au cours des dernières décennies, différentes méthodes ont été développées pour corriger les aberrations, principalement dans le domaine de l'imagerie médicale, afin d'obtenir 2 une image échographique ayant une résolution optimale lorsque celle-ci est dégradée par les hétérogénéités du milieu. Toutefois, dans le cas d'aberrations complexes, les distorsions du front d'onde varient suivant les zones de l'image. Qualifiées de domaines d'isoplanétisme, ces zones sont souvent impossibles à déterminer a priori et constituent de ce fait une limite 5 fondamentale pour l'application des techniques de l'état de l'art. [6] Comme illustré sur les FIGS lA - 1C, les méthodes d'échographie conventionnelle utilisent un réseau 10 de transducteurs piézo-électriques 11 qui peuvent émettre et recevoir des impulsions ultrasonores de manière indépendante. On repère la position de chacun des transducteurs par le vecteur u. Lorsqu'un tel réseau est placé en vis-à- 10 vis d'un milieu 20 que l'on cherche à étudier, celui-ci peut être insonifié de différentes manières. [7] La première manière pour insonifier le milieu à étudier est d'émettre une impulsion ultrasonore depuis un des transducteurs du réseau dont la position est repérée par uji, (FIG. 1 A, schéma de gauche). Cela donne lieu à une onde incidente cylindrique (ou 15 sphérique) divergente pour un réseau 1D (ou 2D) de transducteurs. Celle-ci est réfléchie par les diffuseurs 21 du milieu 20 et le champ rétrodiffusé est enregistré par chacun des transducteurs 11 en fonction du temps (FIG. 1A, schéma de droite). En répétant cette opération avec chaque transducteur utilisé successivement comme source, on mesure l'ensemble des réponses impulsionnelles Knout, uin, t) entre chaque transducteur, où lama 20 désigne la position du détecteur. Ces réponses forment la matrice de réflexion Ruti(t) exprimée dans la base des transducteurs. L'intérêt d'une telle mesure réside dans le fait que cette matrice contient toute l'information sur le milieu étudié, un ensemble d'opérations matricielles pouvant ensuite lui être appliquées à des fins d'imagerie du milieu, par exemple. En revanche, une telle acquisition suppose que le milieu reste fixe pendant toute la durée des 25 mesures, ce qui peut être très limitant dans le cas d'expérience in-vivo. De plus, l'énergie émise par un seul élément piézo-électrique est faible, ce qui peut induire un mauvais rapport signal à bruit. [8] D'autres manières sont connues pour insonifier le milieu à étudier qu'on connaît sous le nom de « formation de voies » (ou « beamforming » en anglais), ces techniques 30 désignant le contrôle cohérent d'ondes par application de lignes à retard. Parmi ces techniques, on connaît les émissions focalisées, technique largement utilisée par les échographes conventionnels. Comme montré par la FIG. 1B (schéma de gauche), il suffit d'appliquer aux transducteurs 11 un jeu de retards appropriés afin de corriger les temps de parcours pour que toutes les impulsions arrivent ensemble au point focal rh, visé. En raison 3 des limites physiques de la diffraction, on concentre les ultrasons dans une zone délimitée par l'ouverture de la sonde échographique. Notons que dans un échographe, la focalisation se fait également à la réception. L'ensemble des échos captés par les éléments 11 du réseau 10 est alors traité pour simuler l'effet d'une lentille en réception, comme décrit sur la Fig. 1B, 5 schéma de droite. Les signaux reçus par les transducteurs sont remis en phase en les décalant temporellement. Ces délais sont identiques à ceux appliqués à l'émission. Dans la phase d'émission, tous les signaux interfèrent au point nu. En réception, les signaux provenant de ce même point rue = fin interfèrent électroniquement par sommation des signaux au temps balistique t = 211uout - rift Ilic. Cette sommation donne le résultat final de la focalisation en 10 réception. La méthode illustrée sur la FIG.1B, dite méthode confocale à double focalisation à l'émission et à la réception permet d'imager directement la réflectivité du milieu avec une résolution latérale limitée par la diffraction, une excellente résolution axiale seulement limitée par la durée de l'impulsion initiale et un excellent contraste. Toutefois, cette méthode est chronophage car il faut en toute théorie focaliser à l'émission en chacun des points du milieu. 15 [009] Une autre technique de formation de voies, développée plus récemment, consiste à insonifier le milieu avec une série d'ondes planes. La FIG. 1C illustre le principe de cette échographie dite en ondes planes, décrite par exemple dans l'article de G. Montaldo et al. « Coherent plane-wave componnding for vety high.frante rate ultrasonography and transient elaslography" (IEEE Trans. Ultrason., Ferroelect. Freq. Control 56 489-506, 2009). Des 20 retards sont appliqués sur chaque signal à l'émission (FIG. 1C, schéma de gauche) pour la formation d'un front d'onde incliné d'un angle Ob, par rapport au réseau de transducteurs 10. A la réception (FIG. 1 C, schéma de droite), le champ rétrodiffusé par le milieu, R(u our, Gin, t) est mesuré par tous les capteurs uuut pour une série d'ondes planes incidentes dont on fait varier l'angle d'incidence O. L'ensemble de ces réponses forment une matrice de réflexion 25 Rua(t) définie entre la base de Fourier (ou base d'ondes planes) en entrée et la base des transducteurs en sortie. Une fois cette matrice enregistrée, les signaux sont décalés temporellement avant d'être sommés de manière cohérente afin de focaliser numériquement les données à l'émission et à la réception pour chaque point ru, de l'image échographique. Le nombre d'acquisitions nécessaires pour former une image échographique est réduit par 30 rapport à l'échographie standard (émissions focalisées), et ceci pour un même niveau de contraste et de résolution de l'image échographique. [0010] Inspirées des techniques de correction des aberrations optiques utilisées en astronomie dans les années 1950, des techniques de focalisation adaptative ont été introduites pour corriger la distorsion de phase en imagerie ultrasonore (voir par exemple D. J. Philips, S. 4 W. Smith, O. T. Von Ramm et F. L. Thurstone, 'A phase compensation technique for B-Mode echoencephalography' Ultrasound in Medicine, 1, 345-404, 1975). [00111 Comme illustré sur la FIG. 2, cette technique de correction s'applique à toutes les méthodes d'imagerie échographique conventionnelle (FIGS 1A-1C). Cette technique suppose 5 une couche aberratrice 22 fine et qui induit de simples décalages temporels sur les transducteurs 11. Ces derniers peuvent alors être compensés en ajoutant une loi de retard supplémentaire aux délais généralement utilisés en formation de voies. Comme illustré schématiquement sur la FIG. 2, cette technique de correction comprend l'excitation (25) du milieu 20 sans prise en compte des aberrations, l'estimation (26 - 27) de la loi d'aberration et 10 la nouvelle émission/réception (28 - 29) en prenant en compte l'aberrateur. En utilisant une boucle de rétroaction entre l'émetteur et le détecteur, le système est itéré en temps réel afin de trouver les décalages temporels optimaux pour corriger l'aberration. Plusieurs approches ont été proposées pour calculer les retards à appliquer sur chacun des transducteurs du réseau à partir du speckle ultrasonore, c'est-à-dire le signal résultant de la réflexion des ultrasons par 15 une répartition aléatoire de diffuseurs sous-résolus. Une méthode est basée sur des corrélations entre les signaux reçus par des transducteurs voisins (voir par exemple M. O'Donnell et S. W. Flax. « Phase aberration mectsurements in medical ultrasound: liman studios », Ultrasonic Imaging 10, 1-11, 1988). Une autre méthode est basée sur l'optimisation de contraste du speckle ultrasonore (L. Nock, G. E. Trahey et S. W. Smith. « Phase aberration 20 correction in medical ultrasound using speckle brightness as a quality factor », J. Acoust. Soc. Am, Vol. 85, 5, pp 1819-1833, 1989). Cependant, ces techniques de correction reposent sur l'hypothèse d'une mince couche aberratrice 22 qui se situe au plus près du réseau de transducteurs 10. Pour un aberrateur volumique, ces méthodes ne sont donc valables que sur un seul domaine d'isoplanétisme. 25 [0012] La notion de domaines d'isoplanétisme est illustrée sur les FIGS. 3A - 3B. Un aberrateur 22 est présent entre le réseau de transducteurs 10 et le plan focal. Dans l'exemple de la FIG. 3A, des retards sont appliqués aux transducteurs 11 du réseau de transducteurs 10 de telle sorte que le front de l'onde acoustique en sortie de l'aberrateur 22 soit corrigé (onde cylindrique). Il en résulte en un point n0 donné du plan focal une tache focale limitée par 30 diffraction. Si ce même front d'onde corrigé est légèrement angulé (FIG. 3A) (ou translaté) pour permettre une focalisation en un point ew tel que I n. - I L, alors le front d'onde issu de l'aberrateur est distordu et la correction appliquée ne permet pas une focalisation limitée par la diffraction en ein. Les points nu et eh, appartiennent dans ce cas à des domaines 5 d'isoplanétisme différents de la zone insonifiée. [00131 Au-delà d'un champ de correction limité à un seul domaine d'isoplanétisme, les techniques de correction d'aberration mentionnées précédemment ne fonctionnent pas en général pour des signaux ultra large bande (tels que ceux utilisés en imagerie ultrasonore), puisque chaque composante spectrale du champ ondulatoire peut subir une distorsion de 10 phase différente. Pour une correction ultra large bande des aberrations en imagerie ultrasonore, des méthodes plus sophistiquées ont été développées, telles que les méthodes de filtrage inverse ou adapté dérivées du concept de retournement temporel. Elles consistent à convoluer les signaux émis et reçus par des filtres spatio-temporels adaptés permettant d'optimiser la focalisation des ondes dans le milieu et de corriger les aberrations de celui-ci. 15 Toutefois, la connaissance de ce filtre nécessite la présence d'un fort diffuseur (C Prada and M. Finlc, « Eigenmodes of the time reversai operator », Wave Motion, 20, 151-163, 1994) ou d'un réseau de détecteurs au sein du milieu que l'on cherche à imager ce qui s'avère être particulièrement invasif (voir par exemple M. Tanter, J.-L. Thomas and M. Fink « Time reversai and the inverse filter », J. Acoust. Soc. Am. 108, 223-234, 2000). 20 [0014] Une approche plus récente tire bénéfice directement de la matrice de réflexion et sa matrice de corrélation pour corriger les aberrations de manière itérative dans le speckle ultrasonore. Ainsi, dans l'article de J-L. Robert et al. (« Green'sfitnclion estimation in speckle using the decomposition of the lime reversai operator: Application to aberration correction in medical imaging », J. Acoust. Soc. Am., 123, 2 (2008)), il s'agit de mesurer la matrice de 25 réflexion Rur entre une base d'émissions focalisée à l'entrée et la base des transducteurs en sortie. L'ensemble des émissions focalisées en des points fin du milieu peut être vu comme un réseau de transducteurs virtuels placés au sein de ce dernier, leur taille respective étant donnée par chaque tâche focale. Ainsi cette matrice de réflexion contient l'ensemble des réponses impulsionnelles R(u ouurin, t) entre chaque transducteur d'un réseau virtuel (rh,) au sein du 30 milieu étudié et chaque transducteur (uout) du réseau réel. En limitant la taille du réseau virtuel à un domaine d'isoplanétisme, il a été montré que le premier vecteur propre de la matrice de corrélation RurRtur, où Rtur est la matrice transconjuguée (ou matrice transposée de la matrice conjuguée) de R', donne le front d'onde permettant de focaliser sur un diffuseur virtuel (étoile artificielle) dont la taille est donnée par la tâche focale initiale. En 6 angulant ce front d'onde, on peut focaliser en chacun des points du domaine d'isoplanétisme. On peut alors réitérer le processus avec cette nouvelle base d'émissions focalisées, reconstruire une nouvelle matrice de réflexion R, et accéder à une loi de correction des aberrations plus précise etc. L'un des inconvénients de cette approche est la nécessité d'itérer 5 le processus, ce qui réduit singulièrement la cadence d'image. Par ailleurs, comme les autres techniques de l'état de l'art, le champ de correction de l'image est, encore une fois, limité à un seul domaine d'isoplanétisme. [00151 Dans le cadre de la présente description, une méthode originale est proposée qui ne repose ni sur la génération d'une étoile artificielle, ni sur une approche itérative. La 10 méthode proposée permet en outre de corriger simultanément les aberrations dans les différents domaines d'isoplanétisme d'une image de speckle, pouvant contenir également des objets plus déterministes (interfaces, gros diffuseurs, etc.) qui donnent lieu à une réflexion spéculaire de l'onde incidente. RESUME 15 [00161 La présente description a pour objet, selon un premier aspect, un procédé de caractérisation ultrasonore non invasive d'un milieu hétérogène, le procédé comprenant : - une étape de génération, au moyen d'au moins un premier réseau de transducteurs, d'une série d'ondes ultrasonores incidentes dans une zone dudit milieu hétérogène , - une étape d'enregistrement d'une matrice de réflexion expérimentale définie entre une 20 base d'émission en entrée et la base des transducteurs en sortie ; - une étape de détermination, à partir de ladite matrice de réflexion expérimentale, d'au moins une première matrice de réflexion large bande, définie dans une base de focalisation en entrée et en sortie, dans une première bande spectrale ; une étape de détermination d'une première matrice distorsion définie dans ladite 25 première bande spectrale, entre ladite base de focalisation et une base d'observation, ladite première matrice distorsion résultant, directement ou après changement de base, du produit matriciel terme à terme de ladite première matrice de réflexion large bande définie entre ladite base de focalisation et une base de correction des aberrations, avec la matrice conjuguée en phase d'une matrice de réflexion de référence, définie pour un 30 milieu modèle, entre les mêmes bases; une étape de détermination à partir de ladite première matrice distorsion, d'au moins une cartographie d'un paramètre physique de ladite zone insonifiée dudit milieu hétérogène. 7 [0017] Comme cela sera décrit plus en détails dans la suite de la description, un paramètre physique du milieu hétérogène peut comprendre par exemple : une loi d'aberrations associée à au moins un domaine d'isoplanétisme contenu dans la zone insonifiée, un paramètre caractéristique de la réflectivité locale du milieu, un paramètre caractéristique de la 5 vitesse du son dans le milieu, un taux de diffusion multiple. [0018] Un milieu hétérogène au sens de la présente description comprend tout milieu comprenant des diffuseurs présentant une rupture d'impédance acoustique par rapport à celle d'un milieu environnant, le dit milieu environnant pouvant également présenter une vitesse du son spatialement inhomogène (e.g milieu multicouches). Les procédés et systèmes décrits 10 dans la présente description permettent notamment une caractérisation ultrasonore non invasive en profondeur d'un tel milieu hétérogène. [0019] Les aberrations d'un front d'onde ultrasonore se propageant à travers un milieu hétérogène correspondent à l'écart entre le front d'onde issu de ce milieu et celui qui serait attendu dans un cas idéal. Ces aberrations peuvent être liées par exemple aux fluctuations 15 spatiales de la vitesse du son dans le milieu. Par exemple, dans le cas de l'échographie, les ondes acoustiques traversent différents types de tissus avant d'atteindre l'organe visé. Les tissus mous présentant des propriétés mécaniques différentes, la vitesse du son est donc loin d'être homogène. Il en résulte une aberration du front d'onde acoustique qui conduit à une distorsion de l'image échographique résultante, et donc à une dégradation de sa résolution et 20 de son contraste. Ces aberrations peuvent être telles qu'elles compromettent l'examen médical. [0020] Une base de correction des aberrations optimale est une base qui maximise la taille des domaines d'isoplanétisme contenus dans la zone insonifiée. Par exemple, une base de correction des aberrations peut être la base des transducteurs dans le cas d'un aberrateur 25 mince situé à proximité du réseau de transducteurs, ou la base d'ondes planes dans le cas d'un milieu multicouches invariant par translation. Dans le cas général d'un aberrateur volumique spatialement inhomogène, il existe plusieurs domaines d'isoplanétisme quelle que soit la base de correction des aberrations considérée. L'étude de la matrice distorsion pourra permettre alors de déterminer les lois d'aberration associées à chacune de ces zones. 30 [0021] La base de focalisation peut être contenue dans un seul plan focal situé à une profondeur donnée du réseau de transducteurs dans le milieu, ou peut être définie dans un volume en considérant simultanément des points de focalisation situés à différentes profondeurs. 8 [0022] La base d'émission dans laquelle est enregistrée la matrice de réflexion expérimentale peut-être indifféremment la base des transducteurs, la base d'ondes planes (« base de Fourier ») ou une base de focalisation quelconque. [0023] Une matrice de réflexion ou la matrice distorsion peuvent être définies entre 5 différentes bases, le changement de base se faisant soit par l'application de lignes à retard dans le domaine temporel soit par des produits matriciels avec des matrices de propagation dans le domaine fréquentiel. Dans la suite de la description, on appelle matrice de réflexion « focalisée » ou matrice distorsion « focalisée » respectivement une matrice de réflexion ou la matrice distorsion définie dans une base de focalisation en entrée et en sortie. 10 [0024] La matrice de réflexion expérimentale enregistrée est une matrice « réelle », c'est- à-dire composée de coefficients réels dans le domaine temporel, les signaux électriques enregistrés par chacun des transducteurs étant des nombres réels. La matrice de réflexion focalisée large bande est en revanche une matrice « complexe », c'est-à-dire composée de coefficients formés de signaux analytiques complexes. Les signaux analytiques complexes qui 15 forment la matrice de réflexion focalisée large bande peuvent être obtenus en ne conservant que les fréquences positives du spectre des signaux électriques enregistrés ; leur argument coïncide avec la phase des ondes réfléchies par le milieu. [0025] Le milieu modèle ou milieu de référence est par exemple un milieu homogène dans lequel la vitesse du son est constante. La matrice de réflexion de référence (ou plus 20 simplement « matrice de référence ») peut être établie théoriquement pour ce milieu de référence comme une matrice de propagation, c'est-à-dire une matrice de réflexion pour les ondes ultrasonores obtenue pour ce milieu de référence avec, en guise de réflecteur, un miroir plan virtuel dans chaque plan de focalisation. Suivant le degré de connaissance a priori du milieu de propagation, le milieu de référence pourra prendre des formes plus élaborées (e.g 25 milieu multicouches etc.). [0026] Le produit matriciel terme à terme ou « produit d'Hadamard » de la première matrice de réflexion large bande définie entre la base de focalisation et la base de correction des aberrations avec la matrice conjuguée en phase d'une matrice de réflexion de référence, définie pour un milieu modèle, entre les mêmes bases, revient à retrancher à la phase de 30 chaque élément de la première matrice de réflexion large bande, la phase de l'élément correspondant de la matrice de réflexion de référence. En effet, les éléments d'une matrice ou conjuguée en phase (ou plus simplement « conjuguée » dans la présente description) d'une première matrice, sont des nombres complexes présentant le même module que ceux de la première matrice mais de phase opposée. On retranche ainsi à la phase de chaque élément de 9 la première matrice de réflexion large bande une composante balistique attendue (définie par le milieu modèle), ce qui permet d'isoler la composante distordue du champ réfléchi par le milieu pour chaque émission focalisée en entrée. Les déposants ont montré qu'une analyse de la matrice distorsion permet notamment de discriminer les domaines d' isoplanétisme contenus 5 dans la zone insonifiée et de déterminer, dans le plan d'observation, une loi d'aberrations associée à chaque domaine d'isoplanétisme [0027] Selon un ou plusieurs exemples de réalisation, la base d'observation est la base de correction des aberrations. On pourra chercher à déterminer la première matrice distorsion entre la base de focalisation et la base de correction des aberrations, par exemple pour 10 quantifier et/ou corriger les aberrations. [0028] Dans un tel cas, la détermination de la première matrice distorsion peut comprendre : un changement de base en entrée ou en sortie de ladite première matrice de réflexion focalisée large bande pour obtenir ladite première matrice de réflexion large bande 15 déterminée entre la base de focalisation et la base de correction des aberrations ; et le produit matriciel terme à terme de ladite première matrice de réflexion large bande avec la matrice conjuguée en phase de la matrice de réflexion de référence définie entre les mêmes bases. [0029] Selon un ou plusieurs exemples de réalisation, la base d'observation est la base de 20 focalisation et on parle de première matrice distorsion « focalisée ». On pourra chercher à déterminer une matrice distorsion « focalisée », c'est-à-dire définie dans la base de focalisation en entrée et en sortie, par exemple pour déterminer une cartographie de la fonction d'étalement (ou PSF pour « Point Spread Function »), la fonction d'étalement étant définie comme la réponse impulsionnelle spatiale entre les bases de focalisation en entrée et 25 en sortie ; à partir de la PSF, on pourra par exemple quantifier localement la qualité de focalisation (i .e le niveau d'aberrations), la vitesse du son ou le taux de diffusion multiple dans la zone insonifiée du milieu étudié. [0030] La première matrice distorsion focalisée peut être obtenue par changement de base à partir de la première matrice distorsion définie entre la base de focalisation et la base 30 de correction des aberrations. [0031] Alternativement, la première matrice distorsion focalisée peut être construite directement à partir de la première matrice de réflexion focalisée large bande au moyen d'une corrélation spatiale entre chaque ligne et/ou colonne de ladite première matrice de réflexion focalisée large bande et la même ligne et/ou colonne de la matrice de réflexion de référence 10 définie dans la même base. A partir de la première matrice distorsion focalisée, il sera possible par simple changement de base de déterminer la première matrice distorsion définie entre la base de focalisation et la base de correction des aberrations. [0032] Selon un ou plusieurs exemples de réalisation, la détermination de ladite première 5 matrice de réflexion focalisée large bande comprend : un changement de base en entrée et en sortie de ladite matrice de réflexion expérimentale exprimée dans le domaine fréquentiel pour former une matrice de réflexion focalisée; une sommation cohérente sur ladite première bande spectrale de ladite matrice de 10 réflexion focalisée pour former ladite première matrice de réflexion focalisée large bande. [0033] Alternativement, la première matrice de réflexion focalisée large bande peut être obtenue en appliquant des décalages temporels sur la matrice de réflexion expérimentale. Toutefois, une telle approche ne permet pas d'analyser les données expérimentales sur 15 différentes bandes spectrales, ce qui peut être préjudiciable en cas de dispersion fréquentielle des aberrations. [0034] Selon un ou plusieurs exemples de réalisation, la détermination d'au moins une cartographie d'un paramètre physique dudit milieu hétérogène comprend : - la détermination des invariants dans ladite base de focalisation de ladite première 20 matrice distorsion, afin d'identifier au moins un premier domaine d'isoplanétisme dans ladite base de focalisation ; - la détermination pour chaque premier domaine d'isoplanétisme identifié, d'une cartographie d'au moins une première loi d'aberration dans ladite base de correction. [0035] Plusieurs méthodes peuvent être mises en oeuvre pour déterminer les domaines 25 d'isoplanétisme et les lois d'aberration associées à partir de la première matrice distorsion. Selon un ou plusieurs exemples de réalisation, la détermination de ces invariants comprend une décomposition en valeurs singulières de ladite première matrice distorsion, une décomposition en valeurs singulières de ladite première matrice distorsion normalisée, c'est à dire dont le module de chacun des éléments aura été normalisé mais dont la phase aura été 30 préservée, ou une décomposition en valeurs singulières d'une matrice de corrélation normalisée de ladite première matrice distorsion, c'est-à-dire une matrice de corrélation de ladite première matrice distorsion dont le module de chacun des éléments aura été normalisé. [0036] Dans l'hypothèse notamment d'une réflexion spéculaire par le milieu hétérogène, les déposants ont montré qu'une cartographie de la réflectivité locale du milieu hétérogène 11 pour l'ensemble de la zone insonifiée peut être obtenue par une combinaison linéaire des vecteurs singuliers de ladite première matrice distorsion. [0037] Les déposants ont également montré que la décomposition en valeurs singulières de la matrice distorsion permet de filtrer le sous-espace bruit (matrice aléatoire sans 5 corrélation entre ses lignes et colonnes) du sous-espace signal (matrice caractérisée par d'importantes corrélations entre ses lignes et/ou ses colonnes), le sous-espace bruit contenant à la fois le bruit expérimental et la contribution incohérente du champ réfléchi induit par les évènements de diffusion multiple se produisant en amont de la base focalisée. [0038] Selon un ou plusieurs exemples de réalisation, le procédé de caractérisation 10 ultrasonore comprend en outre la détermination dans ladite base d'observation d'une première matrice de réflexion corrigée par la ou lesdites première(s) loi(s) d'aberration. Cela permet, notamment dans l'hypothèse d'une réflexion diffuse par l'échantillon, de déterminer une cartographie de la réflectivité locale du milieu, ou « image », corrigée des aberrations. [0039] Sur la base de ladite première matrice de réflexion corrigée de la zone insonifiée, 15 il est possible selon un ou plusieurs exemples de réalisation, de déterminer une matrice distorsion corrigée des aberrations correspondant, dans ladite base de correction des aberrations, au produit matriciel terme à terme de ladite matrice de réflexion corrigée déterminée dans ladite base de correction des aberrations, avec la matrice conjuguée en phase de ladite matrice de réflexion de référence. La matrice distorsion corrigée des aberrations 20 permet par détermination de ses invariants dans la base de focalisation d'affiner la correction dans le premier domaine d'isoplanétisme identifié. Elle permet également d'identifier au moins un deuxième domaine d'isoplanétisme dans la base de focalisation et de déterminer pour chaque deuxième domaine d'isoplanétisme identifié, une cartographie d'une deuxième loi d'aberration dans la base de correction des aberrations. Le procédé peut être ainsi itéré 25 autant de fois que nécessaire suivant le nombre de domaines d'isoplanétisme contenu dans la zone insonifiée pour obtenir une cartographie de la réflectivité locale du milieu, ou « image », corrigée des aberrations. Ce procédé itératif s'applique plus particulièrement dans l'hypothèse d'une réflexion diffuse ou intermédiaire, c'est-à-dire une réflexion mixte spéculaire et diffuse. [0040] Selon un ou plusieurs exemples de réalisation, la détermination d'au moins une 30 cartographie d'un paramètre physique dudit milieu hétérogène comprend la détermination pour au moins un premier domaine d'isoplanétisme identifié, d'une première fonction d'étalement (PSF) dans ladite base de focalisation. La fonction d'étalement correspond à la transformée de Fourier spatiale de la loi d'aberration mesurée dans le plan de correction des aberrations. Elle est invariante spatialement sur chaque domaine d'isoplanétisme. Elle peut 12 être obtenue selon un exemple à partir de la décomposition en valeurs singulières de la matrice distorsion focalisée, c'est-à-dire définie dans la base de focalisation en entrée et en sortie. [0041] Selon un ou plusieurs exemples de réalisation, la détermination d'au moins une 5 cartographie d'un paramètre physique dudit milieu hétérogène comprend la détermination de la vitesse du son dans la zone insonifiée. Il s'agit pour cela d'étudier par exemple l'évolution de la première valeur singulière de la matrice distorsion en fonction du modèle de vitesse utilisé pour modéliser la matrice de référence. Le modèle de vitesse optimal est celui qui maximise la première valeur singulière de la matrice distorsion. Une alternative est d'étudier 10 la matrice distorsion focalisée. Celle-ci permet en effet d'étudier l'évolution de la fonction d'étalement en fonction du modèle de vitesse utilisé pour le calcul de la matrice de référence. Le modèle de vitesse optimal est celui qui minimise la largeur de la fonction d'étalement normalisée par la longueur d'onde. [0042] Selon un ou plusieurs exemples de réalisation, la détermination d'au moins une 15 cartographie d'un paramètre physique dudit milieu hétérogène comprend la détermination d'un taux de diffusion multiple dans la zone insonifiée. Il s'agit pour cela d'étudier par exemple la matrice distorsion focalisée dont on aura au préalable corrigé les aberrations. Le rapport entre la moyenne de l'intensité incohérente (diffusion multiple) obtenue au-delà de la tâche focale et l'intensité au point de focalisation (ligne centrale de la matrice distorsion) 20 donne accès à un taux local de diffusion multiple. [0043] Selon un ou plusieurs exemples de réalisation, le procédé de caractérisation comprend en outre l'identification et/ou la suppression de la composante spéculaire du champ réfléchi et/ou des réflexions multiples induites entre le réseau de transducteurs et les différentes interfaces ou couches du milieu hétérogène. A cette fin, la matrice distorsion peut 25 être projetée dans la base de Fourier à la fois en entrée et en sortie. Dans cette base, les composantes spéculaires et multiplement réfléchies du champ apparaissent pour des couples d'angles incident et réfléchi précis. Elles peuvent donc être aisément filtrées et seule la composante diffuse (speckle) du champ réfléchi est conservée. Cette discrimination de la composante diffuse permet ensuite d'accéder directement aux lois d'aberration à appliquer en 30 entrée et sortie pour corriger la matrice de réflexion et obtenir une image optimale du milieu, ce qui est n'est pas possible si la composante spéculaire prédomine. En outre, le filtrage de la matrice distorsion dans le plan de Fourier permet d'éliminer les réflexions multiples qui polluent par exemple les images échographiques à faible profondeur (échos multiples entre la sonde et les couches superficielles du corps humain). 13 [0044] Selon un ou plusieurs exemples de réalisation, le procédé de caractérisation comprend en outre : une étape de détermination, à partir de ladite matrice de réflexion expérimentale, d'au moins une deuxième matrice de réflexion focalisée large bande, définie dans une base 5 de focalisation en entrée et en sortie, dans une deuxième bande spectrale différente de ladite première bande spectrale ; une étape de détermination d'au moins une deuxième matrice distorsion définie dans ladite deuxième bande spectrale entre ladite base de focalisation et ladite base d'observation, ladite deuxième matrice distorsion résultant, directement ou après 10 changement de base, du produit matriciel terme à terme de ladite deuxième matrice de réflexion large bande, définie entre ladite base de focalisation et ladite base de correction des aberrations, avec la matrice conjuguée en phase d'une matrice de réflexion de référence, définie pour un milieu modèle, entre les mêmes bases. [0045] Les matrices distorsion définies sur différentes bandes spectrales peuvent être 15 utilisées pour déterminer autant de matrices de réflexion corrigées qui pourront ensuite être sommées pour former une seule matrice de réflexion large bande. Cette analyse spectrale permet de corriger les aberrations chromatiques induites par la dispersion fréquentielle de la vitesse du son dans le milieu. La pluralité des matrices distorsion obtenues sur différentes bandes spectrales pourront également être utilisées pour déterminer une cartographie 20 spectrale de la vitesse du son ou du taux de diffusion multiple. [0046] La présente description concerne, selon un deuxième aspect, un système de caractérisation ultrasonique non invasive d'un milieu hétérogène configuré pour la mise en oeuvre de l'ensemble des exemples de procédés de caractérisation ultrasonique tels que décrits précédemment. Le système de caractérisation ultrasonique selon le deuxième aspect 25 comprend : - Au moins un premier réseau de transducteurs configuré pour générer une série d'ondes ultrasonores incidentes dans une zone dudit milieu hétérogène et enregistrer les ondes ultrasonores rétrodiffusées par ladite zone en fonction du temps; - une unité de calcul, couplée audit au moins un premier réseau de transducteurs, et 30 configurée pour : o enregistrer une matrice de réflexion expérimentale définie entre une base d'émission à l'entrée et la base des transducteurs en sortie ; 14 o déterminer, à partir de ladite matrice de réflexion expérimentale, au moins une première matrice de réflexion large bande, définie dans une base de focalisation en entrée et en sortie, dans une première bande spectrale ; o déterminer une première matrice distorsion définie dans ladite bande spectrale, 5 entre ladite base de focalisation et une base d'observation, ladite première matrice distorsion résultant, directement ou après changement de base, du produit matriciel terme à terme de ladite première matrice de réflexion large bande déterminée entre ladite base de focalisation et ladite base de correction, avec la matrice conjuguée en phase d'une première matrice de réflexion de 10 référence, définie pour un milieu modèle, entre les mêmes bases; o déterminer à partir de ladite première matrice distorsion, au moins une cartographie d'un paramètre physique dudit milieu hétérogène. [0047] Le système de caractérisation selon la présente description peut comprendre au moins un réseau de transducteurs à la fois émetteur et récepteur, ou plusieurs réseaux de 15 transducteurs, certains étant dédiés à l'émission, d'autres à la réception des ondes ultrasonores. BREVE DESCRIPTION DES FIGURES [0048] D'autres avantages et caractéristiques de la technique présentée ci-dessus apparaîtront à la lecture de la description détaillée ci-dessous, faite par référence aux FIGS. 20 dans lesquelles : - Les FIGS. lA - 1C (déjà décrites) illustrent des mécanismes d'émission/réception connus pour l'imagerie et la quantification ultrasonore ; - La FIG. 2 (déjà décrite) illustre une technique de correction des aberrations en imagerie ultrasonore, selon l'art antérieur, basée sur la focalisation adaptative ; 25 les FIGS. 3A - 3B (déjà décrites) illustrent de manière schématique la notion de domaines d'isoplanétisme ; la FIG. 4 illustre un exemple de système de caractérisation ultrasonore pour la mise en oeuvre des procédés de caractérisation ultrasonore selon la présente description ; la FIG. 5 illustre une étape du procédé selon la présente description, visant à 30 déterminer, selon un exemple de réalisation, une première matrice de réflexion focalisée large bande à partir d'une matrice de réflexion expérimentale; la FIG. 6 illustre les opérations matricielles de changement de base pour la projection de matrices de réflexion enregistrée dans la base focalisée ; 15 les FIGS. 7A, 7B illustre une étape du procédé selon la présente description, visant à déterminer, selon un exemple de réalisation, une première matrice distorsion à partir de la matrice de réflexion focalisée large bande ; la FIG. 8 illustre des images échographiques d'un phantom sans aberrations et avec 5 aberrations et les matrices distorsion focalisées correspondantes, déterminées à une profondeur donnée et correspondant à une ligne de l'image ; les FIGS. 9A et 9B illustrent respectivement la distribution des valeurs singulières de la matrice distorsion correspondant au phantom avec aberrations montrée sur la FIG. 8 (FIG. 9A) et la phase du premier vecteur propre associé dans le plan de Fourier (FIG. 10 9B); la FIG. 10 illustre un exemple d'utilisation de la matrice distorsion (cas d'une couche d'eau aberrante sur un phantom donnant lieu à un seul domaine d'isoplanétisme) et illustre plus précisément des images échographiques sans/avec correction, les matrices distorsion focalisées correspondantes déterminée à une profondeur donnée et 15 correspondant à une ligne de l'image et une représentation de la fonction d'étalement avec et sans correction moyennée sur la zone insonifiée et issue de la matrice distorsion focalisée ; La FIG. 11 illustre un exemple d'utilisation de la matrice distorsion (sur un phantom de sein humain générant plusieurs domaines d'isoplanétisme) et illustre plus 20 précisément une image échographique conventionnelle et des images échographiques après correction des aberrations à l'aide du premier vecteur de la matrice distorsion et du deuxième vecteur de la matrice distorsion ; les FIGS. 12A - 12C illustrent un exemple d'utilisation de la matrice distorsion pour le filtrage des réflexions multiples et plus précisément, FIG. 12A, un schéma décrivant 25 la configuration expérimentale (plaque de plexiglas placée entre la sonde échographique et le phantom), FIG. 12B, une image échographique pour un modèle de vitesse homogène (c=1800 m/s), avant et après suppression des composantes spéculaires et de réflexions multiples de l'image, des exemples de matrice distorsion focalisée pour une ligne de l'image échographique, avant et après filtrage, les mêmes 30 matrices distorsion projetées dans le base des ondes planes, et FIG. 12C, un schéma expliquant la signature des composantes spéculaire et de réflexion multiple dans la base des ondes planes ; La FIG. 13 illustre un exemple d'utilisation de la matrice distorsion pour une mesure quantitative de la vitesse du son et illustre plus précisément une image échographique 16 du phantom étudié, un exemple de matrice distorsion focalisée pour une ligne de l'image et un modèle homogène de vitesse (c=1540 m/s), la matrice distorsion correspondante entre le plan focal en entrée et le plan de Fourier, la largeur de la fonction d'étalement normalisée par la longueur d'onde en fonction de la vitesse des 5 ondes e, la première valeur singulière normalisée de la matrice distorsion en fonction de cette même vitesse ; la FIG. 14 illustre un exemple d'utilisation de la matrice distorsion pour établir une cartographie de la vitesse du son dans le milieu et illustre plus précisément une image échographique conventionnelle du phantom à travers la plaque de plexiglas pour un 10 modèle de vitesse homogène (c=1540 m/s), et la même image après filtrage des réflexions multiples et optimisation du modèle de vitesse pour chaque ligne de l'image, le profil de vitesse du son intégré sur la profondeur ; la FIG. 15 illustre un exemple d'utilisation de la matrice distorsion pour la quantification locale de la diffusion multiple et montre plus précisément une image 15 échographique d'un phantom standard, un schéma expliquant l'occurrence de diffusion multiple dans la matrice distorsion, la matrice distorsion focalisée pour une ligne de l'image échographique, un exemple de fonction d'étalement déduite de la précédente matrice mettant en évidence un fond incohérent lié à la composante de diffusion multiple, le taux de diffusion multiple en chaque pixel de l'image 20 échographique ; [0049] Dans les différents modes de réalisation qui vont être décrits par référence aux figures, des éléments semblables ou identiques portent les mêmes références. DESCRIPTION DETAILLEE [0050] Dans la description détaillée qui suit, seuls certains modes de réalisation sont 25 décrits en détail pour assurer la clarté de l'exposé mais ces exemples ne visent pas à limiter la portée générale des principes ressortant de la présente description. [0051] Les différents modes de réalisation et aspects décrits dans la présente description peuvent être combinés ou simplifiés de multiples manières. En particulier, les étapes des différents procédés peuvent être répétées, interverties, exécutées en parallèle, sauf précision 30 contraire. [0052] La FIG. 4 illustre un exemple d'un système 40 de caractérisation ultrasonore pour la mise en oeuvre des procédés de caractérisation ultrasonore d'un milieu hétérogène 20, selon la présente description. Le système 40 comprend au moins un premier réseau 10 de 17 transducteurs 11, par exemple un réseau linéaire ou bidimensionnel ; les transducteurs sont par exemple des transducteurs piézoélectriques ultrasonores pouvant se présenter classiquement sous la forme d'une barrette rigide mise en contact avec le milieu 20. Le réseau de transducteurs fait par exemple partie d'un dispositif de sondage 41; le réseau de 5 transducteurs est connecté à une unité de calcul 42, qui peut elle-même être reliée à un dispositif d'affichage 43; l'unité de calcul émet et enregistre des signaux électriques vers /à partir de chacun des transducteurs 11. Les transducteurs ultrasonores transforment ensuite ces signaux électriques en ondes ultrasonores et inversement. L'unité de calcul 42 est configurée pour la mise en oeuvre d'étapes de calcul ou traitement notamment pour la mise en oeuvre 10 d'étapes de procédés selon la présente description. [00531 Dans la FIG. 4 comme dans la suite de la description, il est fait référence à un réseau de transducteurs pour l'émission et la réception, étant bien entendu que, dans un cas plus général, plusieurs réseaux de transducteurs pourront être utilisés simultanément. Ils pourront être à la fois émetteur et récepteur, ou bien seulement émetteur pour certains et 15 seulement récepteur pour d'autres. [00541 Lorsque dans la présente description, il est fait référence à des étapes de calcul ou traitement pour la mise en oeuvre notamment d'étapes de procédés, il est entendu que chaque étape de calcul ou traitement peut être mis en oeuvre par logiciel, hardware, firmware, microcode ou toute combinaison appropriée de ces technologies. Lorsqu'un logiciel est 20 utilisé, chaque étape de calcul ou traitement peut être mise en oeuvre par des instructions de programme d'ordinateur ou du code logiciel. Ces instructions peuvent être stockées ou transmises vers un support de stockage lisible par un ordinateur (ou unité de calcul) et/ou être exécutées par un ordinateur (ou unité de calcul) afin de mettre en oeuvre ces étapes de calcul ou traitement. 25 [00551 Les différents résultats expérimentaux montrés dans la présente description (FIGS 5 à 15) ont été réalisés sur deux phantoms ultrasonores, i.e des échantillons synthétiques reproduisant les propriétés mécaniques des tissus biologiques (c = 1540 m/s) et générant un speckle ultrasonore caractéristique de ces derniers. Le premier phantom permet de quantifier les performances en termes de résolution et de contraste du système d'imagerie. Il comprend 30 en effet des cibles quasi-ponctuelles situées à différentes profondeurs et positions latérales ainsi qu'une inclusion échogène hétérogène de 5 mm de diamètre (FIGS 5 à 10 et 12 à 15). Le second phantom reproduit la forme et les propriétés des tissus du sein et présente des échos brillants mimant la présence de micro-calcifications en son sein (FIG. 11). 18 [0056] La présente description décrit des procédés et systèmes de caractérisation ultrasonore non invasive d'un échantillon hétérogène. Ces procédés et systèmes sont basés sur la détermination d'au moins une première matrice appelée « matrice distorsion » dans la suite de la description. Les FIGS 5 à 7 illustrent selon des exemples des étapes de procédés pour 5 déterminer la matrice distorsion. [0057] Comme illustré sur la FIG. 5, étape 51, le procédé de caractérisation ultrasonore non invasive comprend une étape de génération, au moyen d'un réseau 10 de transducteurs 11 (FIG. 4), d'une série d'ondes ultrasonores incidentes dans une zone dudit milieu hétérogène et une étape d'enregistrement d'une matrice de réflexion expérimentale Rui(t) définie entre une 10 base d'émission i en entrée et la base des transducteurs u en sortie ; la base d'émission i en entrée est par exemple la base des transducteurs u, la base d'ondes planes 0 ou la base focalisée r, comme décrit au moyen des FIGS 1 A - 1C. Dans l'exemple de la FIG. 5, il s'agit par exemple d'une matrice de réflexion expérimentale 12'0(t) enregistrée entre la base d'ondes planes en entrée et la base des transducteurs en sortie. 15 [0058] Le procédé comprend également une étape de détermination, à partir de ladite matrice de réflexion expérimentale Rui(t), d'au moins une première matrice de réflexion focalisée large bande R'(Ao))), définie dans une base de focalisation r en entrée et en sortie, dans une première bande spectrale Ami. Cette étape préalable à la construction de la « matrice distorsion » est illustrée selon un exemple de réalisation, au moyen des FIGS. 5 et 6. 20 [0059] Dans l'exemple de la FIG. 5, une opération 52 de transformée de Fourier discrète permet d'exprimer la matrice de réflexion expérimentale dans le domaine fréquentiel 14.0(0)). Une étape 53 dite de « focalisation large bande » de la matrice de réflexion expérimentale comprend la projection de la matrice de réflexion 12'o(co) dans la base focale en entrée et en sortie afin d'obtenir une matrice de réflexion focalisée R'((]), puis une sommation cohérente 25 sur la bande spectrale Ami. On obtient alors une matrice de réflexion focalisée large bande dont un exemple est illustré (image 54) à une profondeur z donnée. [0060] La FIG. 6 illustre la projection de la matrice de réflexion Ruo(c)) dans la base focale r en entrée et en sortie. A titre d'exemple, l'image 63 représente la phase de la matrice Ride)) à la fréquence f=3.48 MHz, acquise dans le premier phantom, comprenant des cibles 30 quasi-ponctuelles ainsi qu'une inclusion échogène hétérogène de 5 mm de diamètre et dont une image échographique est montrée sur la FIG. 5 (image 56). Cela revient à effectuer une formation de voies en émission (61) et en réception (62) afin de sommer de façon cohérente les échos venant d'un même diffuseur. Mathématiquement, ces opérations peuvent être décrites à l'aide de matrices de propagation décrivant les trajets aller et retour des ondes dans 19 le milieu. Le passage de la base des ondes planes incidentes à la base de focalisation est ainsi modélisé par la matrice G'0(w) = [G'o (rio, 0,n, w)] , dont les éléments sont donnés par ré) = exp [jk(z,'cos0,' + xt'sin0,')] (1) avec r,' = (x,',z,') le vecteur position des points de focalisation dans le milieu, 0,, l'angle 5 d'incidence de l'onde émise et k = / c le nombre d'onde. Le trajet des échos réfléchis depuis le milieu vers le réseau de transducteurs peut être lui décrit par la matrice de propagation Go (w) = [Go (l'out, nota, (0)] dont les éléments sont donnés par exp Huma-rotaMi zout w)GO(rout, bout, w) (2) 1111out-rout Il Iluoui-rmell avec Font = (xouc, z out) le vecteur position des points de focalisation dans le milieu et uont le 10 vecteur position des transducteurs. La dernière équation est une traduction matricielle de l'intégrale de Rayleigh Sommerfeld. Son premier terme est la fonction de Green 2D de l'équation d'onde, ou plutôt ici sa forme asymptotique. Le terme de droite est le cosinus directeur entre la normale du réseau de transducteur et le vecteur radial (uout - Fout) - 100611 Les opérations de formation de voies à l'émission et la réception peuvent être 15 interprétées comme une rétropropagation des échos vers les diffuseurs qui leur sont associés. Dans le domaine fréquentiel, cela se traduit par une conjugaison en phase des matrices de propagation. La matrice de réflexion focalisée R' peut ainsi être déduite de la matrice de réflexion expérimentale par le produit matriciel suivant : Rrr(w) = GO(oe) x 12,3(w) x Gort(cd) (3) 20 où le symbole x représente le produit matriciel et l'exposant t désigne la matrice transposée conjuguée. Bien entendu, la projection de la matrice de réflexion dans les bases focales en entrée et en sortie pourrait se faire également à partir d'une matrice de réflexion expérimentale exprimée en entrée dans une autre base, par exemple la base des transducteurs. Il suffirait pour cela d'utiliser dans l'équation (3) à la place de la matrice de propagation GO la 25 matrice de propagation Go. [0062] L'étape suivante consiste à effectuer une somme des matrices de réflexion obtenues sur la bande spectrale ​ de telle sorte à obtenir une matrice de réflexion focalisée large bande R'(r) = Eweam, Rrr(w) el" (4) 30 T est le temps d'écho associé au réseau virtuel de transducteurs dans la base de focalisation. Chaque coefficient R(ro,,,t, %or) de la matrice de réflexion focalisée large bande correspond au signal analytique qui serait enregistré par le transducteur virtuel en tout immédiatement après l'émission d'une impulsion depuis un transducteur virtuel en ro, (schéma 55, FIG. 5). 20 La diagonale de la matrice R'(T = 0) au temps balistique r = 0 correspond à au signal confocal avec focalisation à l'émission et la réception sur le même point. Elle donne donc directement une ligne de l'image échographique (image 56) si la base de focalisation se limite à un plan focal situé à une distance z du réseau réel de transducteurs. Si la base de focalisation 5 contient au contraire tous les points de focalisation du volume à imager, la diagonale de R'(T = 0) donne directement l'ensemble de l'image échographique. Dans ce dernier cas, la matrice de réflexion est une matrice par blocs dont seuls les blocs diagonaux sont non nuls. Ces derniers correspondent aux matrices de réflexion focalisée à chaque profondeur z du milieu. Dans la suite de la description, afin d'alléger les notations, la matrice de réflexion 10 large bande sera appelée simplement Rrr, il ne sera plus fait mention du temps r car ce dernier sera toujours fixé au temps balistique r = O. Cependant, pour certaines applications, notamment pour la correction d'aberrations d'ordre élevé, l'étude de la diffusion multiple ou la mesure de la vitesse du son, les matrices réflexion et distorsion pourront être également examinées à des temps T # 0. 15 [00631 Les FIGS. 7A, 7B illustrent une étape du procédé selon la présente description, visant à déterminer, selon un exemple de réalisation, une première matrice distorsion D à partir de la matrice de réflexion focalisée large bande Rrr. [0064] Afin d'estimer la composante aberrée des fronts d'onde, on commence par choisir une base optimale de correction des aberrations. Par exemple, dans le cas d'un aberrateur 20 mince situé à proximité du réseau de transducteurs, la base des transducteurs semble être la plus appropriée. Dans le cas choisi dans cet exemple d'illustration du procédé de caractérisation, le milieu hétérogène est un milieu multicouches invariant par translation et la base d'ondes planes semble plus appropriée. C'est par exemple le cas de l'échographie du foie, où les ondes acoustiques doivent traverser une succession de couches de muscle et de 25 tissus adipeux. Cet exemple est considéré par la suite mais la présente description ne se restreint en aucun cas à une correction des aberrations depuis la base de Fourier. [00651 L'image 71 (FIG. 7A) représente un exemple de matrice de réflexion focalisée large bande R' enregistrée en présence d'aberration générée par l'insertion d'une couche d'eau (c = 1480 m/s) de 10 mm d'épaisseur entre le réseau de transducteurs et le phantom 30 (c = 1540 m/s). Dans un premier temps, un changement de base de la matrice de réflexion Rrr en sortie (ou en entrée) dans la base des ondes planes permet de former la matrice de réflexion large bande Roi., telle que Re, = 'G 1(@) X Rrr (5) 21 où on considère la matrice Go' à la fréquence centrale co, de la bande spectrale Avr. Physiquement, chaque ligne de la matrice Ro, ainsi obtenue correspond au front d'onde rétrodiffusé en champ lointain pour chaque point insonifié n.. Chaque front d'onde possède deux composantes : une composante géométrique, lié à la diffraction de l'onde entre le plan 5 des transducteurs et le plan focal, et une composante distordue, liée aux aberrations subies par le front d'onde à l'aller et au retour. Afin de séparer les deux composantes, la matrice Ro, est comparée à une matrice de référence qu'on obtiendrait dans un milieu modèle, par exemple un milieu homogène sans aberration (possédant uniquement la composante géométrique). Cette matrice de référence dans le cas d'un milieu homogène sans aberrations n'est autre que 10 la matrice de propagation en milieu homogène Go On crée alors une nouvelle matrice (FIG. 7A, étape 72), appelée matrice distorsion dans la présente description, à partir du produit d'Hadamard ou produit matriciel terme à terme entre ladite première matrice de réflexion large bande Ro, déterminée entre la base de focalisation et la base de correction des aberrations, et la matrice conjuguée en phase de la matrice de réflexion de référence, Go'''((oc.), 15 définie entre les mêmes bases à la fréquence centrale w, : Dor = Ro, Gj(o.),) (6) où le symbole o représente le produit matriciel terme à terme (produit matriciel d'Hadamard) et l'exposant t désigne la matrice transposée conjuguée. Pour rappel, les éléments de la matrice Gort conjuguée en phase de la matrice eGio sont des nombres complexes présentant le 20 même module que ceux de `Gij mais de phase opposée. En termes de coefficients matriciels, cette équation s'écrit : D (0out. rtn) = R (006t, rtn) X Go' (ri n, 0 11, - out, (oc) (7) [0066] L'opération effectuée ci-dessus est illustrée par la FIG. 7B. Elle consiste à soustraire au front d'onde mesuré sa composante géométrique déterministe et ainsi isoler sa 25 composante distordue. Chaque colonne de la matrice Dor contient cette composante aberrée du front d'onde vue depuis le champ lointain pour chaque point de focalisation rio Ainsi, en pratique, pour isoler la composante distordue du front d'onde réfléchi, on retranche à la phase du champ réfléchi, la phase du champ qui serait obtenu idéalement en l'absence d'aberrations, pour un milieu modèle. 30 [0067] Dans le cas d'un milieu modèle inhomogène, l'expression de la matrice de réflexion de référence Go' devient plus compliquée que celle donnée en équation (1). Par exemple, dans le cas d'un milieu multicouches, la matrice G;) pourra être calculée analytiquement par un produit matriciel entre les différentes matrices de transmission associées à la propagation de l'onde dans chacune des couches. Pour un milieu modèle plus 22 compliqué, typiquement non invariant par translation (e.g interface courbe), la matrice pourra être déterminée au moyen d'une simulation numérique ou d'un calcul semi-analytique. A noter que si les coefficients de la matrice de réflexion de référence présentent un module inhomogéne, la matrice de référence pourra être normalisée en imposant un module constant 5 pour chacun de ces coefficients mais en préservant leur phase. [0068] L'image 73 sur la FIG. 7A montre la phase de la matrice distorsion Der déduite de la matrice de réflexion dont le module est représenté sur l'image 71. Cette dernière correspond à la ligne de l'image échographique 83 (FIG. 8) réalisée sur le phantom en présence d'une couche d'eau aberrante. Le schéma 74 illustre les fronts d'onde aberrés isolés 10 pour chaque point du plan focal et la manière dont ils sont stockés dans la matrice distorsion. [0069] La matrice distorsion peut également être étudiée dans le plan focal. Une matrice distorsion « plan focal » Drr peut être obtenue à partir de la matrice distorsion exprimée dans le plan d'observation Der par le changement de base suivant : Drr = 6[1(Wc)Der (8) 15 où Cru est une matrice de propagation normalisée dont les éléments sont donnés par ô-W0Fulk)°"t'w) = exp [jk4',sin00'(1 (9) (rnnt, °ont, (-±-)) "d'efre'llt-119netm., L'équation (8) se réécrit donc en termes de coefficients matriciels de la manière suivante : D (l'ore, = Eor,', D (0 out. rin) exp[-jitcfnutsin0 out] = R(rout , ) (10) où k, est le nombre d'onde à la fréquence centrale ai,. Dans le cas présent où la base de 20 correction des aberrations est la base de Fourier, la matrice de distorsion focalisée D' est simplement déduite par une transformée de Fourier spatiale de Der en sortie. Toutefois, il faut garder à l'esprit que la relation entre les deux matrices sera différente [équation (8)] si le plan de correction n'est plus le plan de Fourier. [0070] Les colonnes de Drr dont l'image 76 représente le module correspondent au 25 champ réfléchi dans le plan image recentré sur le point de focalisation ri.. La matrice Drr donne donc l'évolution de la fonction d'étalement en réflexion en chaque point de l'image échographique, la fonction d'étalement étant la réponse impulsionnelle spatiale du système d'imagerie, c'est-à-dire la tache focale recentrée sur le point de focalisation. Comme cela sera décrit par la suite, la fonction d'étalement permet de quantifier et caractériser les aberrations 30 ainsi que la diffusion multiple induites par l'échantillon en amont du plan focal. [0071] A titre d'illustration, la FIG. 8 illustre une image échographique 81 du phantom obtenu sans aberrations, au moyen d'un schéma expérimental rappelé sur le schéma 80. L'image 82 illustre la matrice distorsion focalisée déterminée à une profondeur donnée et 23 correspondant à une ligne de l'image 81; l'image 83 est une image échographique du phantom obtenu en présence d'aberrations (couche d'eau aberrante 22, schéma 85) et l'image 84 illustre la matrice distorsion focalisée déterminée à une profondeur donnée et correspondant à une ligne de l'image 83. Dans l'expérience sans aberrations, toute l'énergie réfléchie par le 5 milieu apparaît sur la ligne centrale de la matrice distorsion (rie= re't) : la largeur 8 de la fonction d'étalement est proche de la limite imposée par la diffraction (8-2,z/D, avec D la taille du réseau de transducteurs) ce qui explique l'excellente qualité en terme de résolution et de contraste de l'image échographique. Au contraire, en présence de la couche d'eau aberrante, l'énergie réfléchie s'étale en dehors de la ligne centrale. La dégradation de la 10 réponse impulsionnelle spatiale implique une image échographique nettement moins bien résolue. Cet exemple montre comment la matrice distorsion « plan focal » permet de quantifier la qualité de la focalisation dans le milieu, indépendamment d'une analyse qualitative de l'image échographique. [0072] Pour aller plus loin dans l'analyse de la matrice distorsion Der, on peut 15 avantageusement considérer d'une part le cas d'une réflexion par l'échantillon principalement spéculaire et d'autre part le cas d'une réflexion principalement diffuse. En pratique, on peut se trouver dans des cas intermédiaires qui combinent réflexion diffuse et réflexion spéculaire. Dans les deux cas, la matrice distorsion Der présente des corrélations entre ses colonnes. Ces corrélations correspondent à la répétition du même motif de distorsion pour les fronts d'onde 20 issus du même domaine d'isoplanétisme. Pour un échantillon induisant une réflexion diffuse, par exemple de type speckle (répartition aléatoire de diffuseurs non résolus), on va pouvoir étudier la matrice de corrélation normalisée de la matrice distorsion Der pour extraire de manière plus efficace ces mêmes corrélations. Dans tous les cas, comme cela sera décrit plus en détails par la suite, la recherche des invariants de la matrice Der, par exemple au moyen 25 d'une décomposition en valeurs singulières (SVD), permet d'extraire la transmittance complexe de l'aberrateur pour chaque point du plan focal et ainsi corriger de manière optimale les aberrations dans chaque domaine d'isoplanétisme contenu dans la zone insonifiée. [0073] Dans l'hypothèse d'un échantillon induisant une réflexion diffuse, on décrit ci- 30 dessous des exemples d'utilisation de la matrice distorsion pour déterminer une cartographie d'une ou plusieurs lois d'aberrations dans la base d'observation, notamment pour établir une cartographie de la réflectivité locale de la zone insonifiée. [00741 Pour cela, on recherche les invariants de la matrice distorsion, autrement dit les lois d'aberrations invariantes spatialement sur les domaines d'isoplanétisme de la zone 24 insonifiée. Différentes méthodes sont connues de l'homme du métier pour rechercher les invariants d'une telle matrice, comme par exemple la décomposition en valeurs singulières (ou « SVD » selon l'abréviation anglo-saxonne « Singular Value Decomposition ») ou l'analyse en composantes principales (« PCA » selon l'abréviation anglo-saxonne « Principal 5 Component Analysis »). [0075] La décomposition en valeurs singulières est un outil puissant pour extraire les corrélations entre les lignes ou les colonnes d'une matrice. Mathématiquement, la SVD de la matrice Dor, de taille N x M, s'écrit de la manière suivante : Dor=UxExVt (11) 10 U et V sont des matrices unitaires de taille N x M dont les colonnes il; et V; correspondent aux vecteurs propres en sortie et en entrée. Chaque vecteur propre en sortie tri est défini dans le plan de Fourier repéré par l'angle 0,,'t. Chaque vecteur propre en entrée Vi est donc défini dans le plan focal repéré par le vecteur no. E est une matrice de taille N x M dont seuls les éléments diagonaux sont non nuls : 15 E = o-2, o-N2) (12) Les éléments diagonaux de la matrice E sont les valeurs singulières u, de la matrice Dor qui sont réelles, positives et classées dans un ordre décroissant : > a2 > > aN 2 (13) Les coefficients D( Bout, ri') de la matrice Dor s'écrivent donc comme la somme suivante : 20 D(Oout, rin) Ei=1 ai Ui (eout)W(rin) (14) avec L = min (M, N). [0076] La SVD décompose principalement une matrice en deux sous-espaces : un sous- espace signal (une matrice caractérisée par d'importantes corrélations entre ses lignes et/ou ses colonnes) et un sous-espace bruit (une matrice aléatoire sans corrélation entre ses lignes et 25 colonnes). Le sous-espace signal est associé aux valeurs singulières les plus grandes alors que le sous-espace bruit est associé aux valeurs singulières les plus faibles. D'une part, la SVD de D va donc permettre de filtrer le sous-espace bruit qui contient à la fois le bruit expérimental et la contribution incohérente du champ réfléchi induit par les évènements de diffusion multiple. D'autre part, chaque état singulier du sous-espace signal va permettre d'extraire 30 suivant le vecteur propre en sortie U; la distorsion subie par l'onde dans la base de correction des aberrations pour chaque domaine d'isoplanétisme de l'image. En particulier, dans le cas d'un seul domaine d'isoplanétisme, on peut montrer que le premier vecteur propre en sortie Ul donne directement la transmittance complexe A de l'aberrateur entre le plan de correction et le plan focal : 25 U1 (001a) = A(Bout) (15) [0077] Les figures suivantes illustrent différentes applications de la matrice distorsion, qu'elle soit définie dans les bases focales ou entre la base focale et la base de correction des aberrations. 5 [0078] Les FIGS. 9A et 9B montrent le résultat de la SVD appliquée à la matrice distorsion Do, représentée sur l'image 84 de la FIG. 8. Le spectre des valeurs singulières o-, de Der (FIG. 9A) montre une domination de la première valeur singulière 0-1. Cela indique que nous sommes en présence d'un seul domaine d'isoplanétisme. Le premier vecteur propre Ui (FIG. 9B) associé montre une phase parabolique, caractéristique d'un défaut de mise au point 10 expliqué par la présence de la couche d'eau dans lequel la vitesse du son (1490 m/s) est sensiblement différente par rapport à celle du phantom (1540 m/s). [0079] Pour corriger la matrice de réflexion large bande R,.,. (71, FIG. 7A), on commence dans cet exemple à projeter la matrice R' dans le plan de correction des aberrations (ici le plan de Fourier) : 15 Roo = %r/coi x R' x Gô (16) La correction des aberrations s'effectue ensuite par application du conjugué en phase de U1 en entrée et sortie de la matrice Roo WelE, = exp(j x arg{U }) o R99 o expo x argfUll) (17) où Woo est la matrice de réflexion corrigée et la fonction argf...1 désigne la phase du vecteur 20 entre accolades. La relation (18) se traduit pour les coefficients matriciels de la façon suivante : fe (Bout, 0 Én) = (Ji (0 out) R (0 out, 0 in)LYin) (18) Physiquement, l'opération de conjugaison de phase consiste à réémettre un front d'onde modulé par une phase opposée à celle de la distorsion mesurée et à appliquer une correction 25 similaire en sortie. Cette opération permet alors de compenser parfaitement les distorsions de phase cumulées par l'onde lors de son trajet à l'aller et au retour. Une matrice corrigée R'' peut alors être déduite de Wolip par changement de base en entrée et sortie depuis le plan pupille vers le plan focal. Une image confocale corrigée /' peut être déduite de la diagonale = ro',) de la matrice R'': 30 "rin) = R(rin,fin) (19) r(rin) est alors un estimateur fidèle de la réflectivité p(rin) de l'échantillon. [0080] La FIG. 10 illustre les bienfaits de cette correction en comparant l'image échographique aberrée 101 et l'image échographique 103 obtenue suivant la correction décrite par l'équation (17). On observe un gain qualitatif entre ces deux images notamment grâce aux 26 cibles échogènes qui apparaissent bien plus contrastées et mieux résolues après correction des aberrations. [0081] En échographie, on n'a toutefois pas la plupart du temps de cibles échogènes pour juger de la résolution de l'image. Pour quantifier cette dernière, on peut avoir recours à la 5 matrice distorsion focalisée. Celle déduite de la matrice de réflexion de départ, D', est présentée sur l'image 102 à une profondeur particulière. Celle déduite de la matrice de réflexion corrigée, D', est présentée sur l'image 104 à la même profondeur. Dans le deuxième cas, l'énergie est bien plus concentrée autour de la ligne centrale de la matrice distorsion focalisée que dans le premier cas. 10 [0082] A partir de la matrice distorsion focalisée, on peut également déterminer la fonction d'étalement. Cette dernière correspond à la transformée de Fourier spatiale de la loi d'aberration mesurée dans le plan de correction des aberrations. Elle est invariante spatialement sur chaque domaine d'isoplanétisme. Elle peut être également obtenue à partir de la décomposition en valeurs singulières de la matrice distorsion focalisée. Si la matrice GO est 15 unitaire (ce qui est le cas si les bases de Fourier et de focalisation sont conjuguées l'une de l'autre), alors la matrice D' présente une décomposition en valeurs singulières similaire à Dor. Les vecteurs propres en sortie, que l'on nommera W1, sont directement reliés aux vecteurs propres U1 par une simple transformée de Fourier : IV; = GO-1(coc) x Ili (20) 20 [0083] Le schéma 105 de la FIG. 10 compare le module des fonctions d'étalement ainsi obtenues avant et après correction. Le gain en termes de résolution est évident et la limite de diffraction est atteinte après correction des aberrations. [0084] Dans l'expérience réalisée, les aberrations mesurées ne dépendent pas de la fréquence car les vitesses du son dans le phantom et dans l'eau varient peu sur la bande de 25 fréquence étudiée. En pratique, ce n'est pas toujours le cas et une analyse fréquentielle de la matrice distorsion peut s'avérer judicieuse pour déterminer le filtre adaptatif spatiofréquentiel (ou de manière équivalente spatio-temporel) corrigeant de manière optimale les aberrations. La matrice distorsion D peut en effet être étudiée sur différentes bandes spectrales aa centrées autour de fréquences centrales distinctes cup. Une loi d'aberration 1.11(cop) 30 pourra être extraite sur chacune de ces bandes spectrales et une matrice de réflexion corrigée R'' (co p) pourra en être déduite [équation (17)]. Une somme cohérente des matrices de réflexion corrigées de manière indépendante sur chacune des bandes spectrales Awi permet d'accéder à la matrice de réflexion corrigée ultra large bande 1:11': Rrr = Loci Rrr (coc,t) (21) 27 Cette correction fréquentielle des aberrations est équivalente en formation de voies à la convolution des signaux focalisés émis et mesurés par le filtre spatio-temporel Ui(r), où Ui Di ((Jijel wc.ir (22) [0085] Jusqu'à maintenant, nous nous sommes placés dans le cas d'une couche aberrante 5 invariante par translation, ne générant ainsi qu'un seul domaine d'isoplanétisme. Nous décrivons maintenant un deuxième exemple d'utilisation de la matrice distorsion dans le cas d'un aberrateur volumique et non invariant par translation, générant ainsi plusieurs domaines d'isoplanétisme. [0086] Dans ce cas, une approche itérative (post-traitement) de la correction des 10 aberrations est privilégiée. L'étape #0 du processus consiste à corriger la matrice de réflexion large bande 120, en sortie par le conjugué de la phase du premier vecteur propre Iii de la matrice distorsion Dur: R81) = exp(j x arg{Ui}) 0129' (23) où Roi(rl. est la matrice corrigée résultante . Cette étape initiale permet d'effectuer une 15 correction globale des aberrations sur l'ensemble de l'image. [0087] L'étape #1 consiste à recalculer une nouvelle matrice distorsion Di.' déduite de Wriu)définie cette fois-ci entre le plan de Fourier en entrée et le plan focal en sortie. La nature aléatoire de l'objet rend plus judicieuse l'étude de la matrice de corrélation de D(rlii) dans le plan de Fourier, à savoir B(1) = tD(R.10)D,(10)*. Une décomposition en valeurs singulières de la 20 phase de cette matrice len, exp[jarg{B(')}], donne un nouvel ensemble de vecteurs propres UP-) associés à chaque domaine d'isoplanétisme de l'image. La matrice de réflexion correspondante peut donc être corrigée cette fois-ci en entrée : leo) = R(r10) o exp (-jargfUill) (24) [0088] Les étapes suivantes consistent à reproduire le même processus en alternant 25 correction des aberrations résiduelles en sortie (itérations paires) et entrée (itérations impaires). En revanche, à chaque étape, la correction s'effectue toujours avec le premier vecteur propre 14n) de la phase de la matrice de corrélation, exp[jarg{B(1011, de la matrice distorsion dans le plan pupille. Le choix du domaine d'isoplanétisme est en effet réalisé à l'étape #L Suivant que la correction soit en sortie ou en entrée, la matrice de corrélation dans 30 le plan pupille BO') est donnée par : B(") -D(On2Det pour n pair (sortie) et BM = tpi(10) Drio)( pour n impair (entrée). A chaque étape du processus, une image du champ de vision peut être déduite de la diagonale de la matrice de réflexion 12('nr) exprimée dans le plan focal. En 28 pratique, quelques itérations suffisent à obtenir une correction optimale pour le domaine d'isoplanétisme sélectionnée. Une image de l'ensemble du champ de vision peut ensuite être obtenue en combinant les corrections déterminées pour chaque domaine d'isoplanétisme. [0089] La FIG.11 illustre les bienfaits de la méthode selon la présente description au 5 moyen d'une expérience sur un phantom de sein générant plusieurs domaines d'isoplanétisme et comportant des microcalcifications entre 20 et 25 mm de profondeur. Sur la FIG. 11 est représentée une image échographique conventionnelle 111, une image échographique 112 après correction des aberrations à l'aide du premier vecteur propre 14(1) de la matrice de corrélation normalisée exp[jarg03(1)11, et un zoom 113 sur une partie de l'image 10 échographique précédente. Est également représentée une image échographique 116 après correction des aberrations à l'aide du second vecteur 141) , et sur l'image 117, un zoom sur une partie de l'image échographique précédente. On observe à partir de ces images que le premier vecteur propre 1_41) est associé au domaine d'isoplanétisme contenant les microcalcifications entre 20 et 25 mm de profondeur. L'image échographique (112) est en 15 effet plus contrastée et mieux résolue à cette profondeur. Au contraire, le second vecteur propre 1.1?) corrige de manière optimale l'image (117) à de plus faibles profondeurs (entre 10 et 15 mm). [0090] La connaissance des lois d'aberrations en chaque point de l'échantillon peut non seulement être mise à profit pour l'imagerie mais également pour focaliser physiquement les 20 ondes ultrasonores en n'importe quel point de l'échantillon. Par exemple, sur le dispositif expérimental de la FIG. 2, il peut s'agir d'appliquer au réseau de transducteurs une convolution du signal focalisé initial par un filtre spatio-temporel déduit des vecteurs propres de la matrice distorsion. En présence d'une seule aire d'isoplanétisme, il s'agira de considérer la transformée de Fourier du premier vecteur propre Ui [équation (22)]. Notons que si la loi 25 d'aberration est stable fréquentiellement, le filtre spatio-temporel consistera simplement en l'application d'une loi de retard At directement proportionnelle à la phase du vecteur propre U1 dans la base des transducteurs at(u) = -arg{t/1(u)}/0.), (25) En présence de plusieurs aires d'isoplanétisme, la détermination des lois d'aberrations 30 s'effectue au moyen d'un processus itératif (cf paragraphe précédent). Il s'agira d'accumuler suivant un vecteur W l'ensemble des corrections d'aberrations obtenues en entrée de l'échantillon : W = exp (-j Earg {UM + nt-1 argflen+111) (26) 29 [0091] En régime de réflexion spéculaire et dans le cas d'une image échographique contenue dans un seul domaine d'isoplanétisme, on peut montrer que le premier vecteur propre en sortie Eh de la matrice distorsion donne la distorsion induite par l'aberrateur cumulée sur l'aller et le retour : 5U1(Bout) = A (Ocnit)A (et») 8(90ut + Ont) (27) avec A(Bout) et A(Oin) les distorsions subies par le front d'onde à l'aller et au retour et projetées dans le plan de Fourier. La distribution de Dirac S dans l'équation précédente traduit le fait, qu'en régime de réflexion spéculaire, une onde plane avec un angle d'incidence Ou, se réfléchit en une onde plane avec un angle Bout = Bin en sortie. Par ailleurs, on peut montrer 10 que le vecteur propre en entrée V1 donne quant à lui directement accès à la réflectivité p de l'objet : Vi (fin) = P(rin) (28) La matrice de réflexion peut être corrigée des aberrations en lui appliquant une seule fois en entrée ou en sortie le conjugué en phase du vecteur propre U1 : 15 = exp(-j x arg{Ui})oltur (29) Dans le cas d'une image échographique contenant plusieurs un nombre P>1 de domaines d'isoplanétisme, la matrice distorsion est de rang P. Les vecteurs propres en entrée V; décomposent l'objet dans le plan focal sur différents domaines d'isoplanétisme et sont associés à des lois d'aberration 1.1; différentes dans le plan de Fourier. La combinaison linéaire 20 des modules des vecteurs propres Vi pondérés par les valeurs propres 6i2 associées donne finalement accès à une image du champ de vision corrigée des aberrations induites en amont IP(rin)12 =Ep 1at2 11A(ri.)12 (30) [0092] Cependant, un réflecteur n'est jamais parfaitement spéculaire et sa composante diffuse ne sera alors pas corrigée de manière optimale par l'équation (29). Par ailleurs, un 25 réflecteur spéculaire peut induire des réflexions multiples entre ses interfaces susceptibles de brouiller l'image échographique. Il est donc nécessaire de pouvoir séparer les composantes spéculaire et diffuse du champ réfléchi par le milieu. Selon un ou plusieurs exemples de réalisation, le procédé de caractérisation comprend donc également l'identification et/ou la suppression de la composante spéculaire du champ réfléchi et/ou des réflexions multiples 30 induites entre le réseau de transducteurs et les différentes interfaces ou couches du milieu hétérogène. A cette fin, la matrice distorsion peut être projetée dans la base de correction des aberrations à la fois en entrée et en sortie. Une matrice distorsion « plan de Fourier » D90 peut être obtenue à partir de la matrice distorsion exprimée dans le plan d'observation De, par le changement de base suivant : 30 Doe = Do, tCo(coc.) (31) [0093] L'équation précédente se réécrit donc en terme de coefficients matriciels de la manière suivante: D (Bout, (9' in) - (Bout, Tin) exp [jk(zincoseh, + xi', sine/in)] = R(0 out. oiu + Bout) (32) 5 [0094] Comme illustré sur la FIG. 13C, les composantes spéculaires et multiplement réfléchies apparaissent pour des couples précis d'angle incident et réfléchi, Oir, et 0out, tel que Oin + °out = a0 (33) où AO correspond à l'angle entre le réflecteur spéculaire et le réseau de transducteurs. AO est 10 par exemple nul si le réflecteur spéculaire est parallèle au réseau de transducteurs. Les composantes spéculaires et multiplement réfléchies peuvent donc être aisément filtrées et seule la composante diffuse (speckle) du champ réfléchi pourra être conservée pour déterminer de manière précise les aberrations induites par la plaque de plexiglas et les corriger. Cette discrimination de la composante diffuse permet ensuite d'accéder directement 15 aux lois d'aberration à appliquer en entrée et sortie pour corriger la matrice de réflexion et obtenir une image optimale du milieu [équation (17)], ce qui est n'est pas possible si la composante spéculaire prédomine [équation (29)]. En outre, le filtrage de la matrice distorsion dans le plan de Fourier permet d'éliminer les réflexions multiples qui polluent par exemple les images échographiques à faible profondeur (échos multiples entre la sonde et les couches 20 superficielles du corps humain). [0095] Les FIGS. 12A - 12C illustrent l'utilisation de la matrice distorsion pour le filtrage des réflexions multiples et son application ensuite à la correction des aberrations et à une cartographie de la vitesse du son. La FIG. 12A représente le schéma de la configuration expérimentale, à savoir une plaque de plexiglas 121 placée entre la sonde échographique 41 et 25 le phantom 20. L'image échographique correspondante 122 est construite à l'aide d'un modèle de vitesse homogène (c=1800 m/s). Cette image est fortement dégradée par les réflexions multiples et les aberrations induite par la plaque de plexiglas. La FIG. 12B représente différentes formes de la matrice distorsion pour une profondeur de l'image échographique (z=25 mm). L'image 124 représente le module de la matrice de distorsion 30 focalisée avec une fonction d'étalement particulièrement dégradée par la présence de la plaque de plexiglas. L'image 125 montre la matrice distorsion dans le plan de Fourier. On voit clairement apparaître une forte composante spéculaire pour Oin -H 0 - out = O. L'image 126 montre la même matrice distorsion après application d'un filtre spéculaire qui consiste à éliminer toutes les composantes du champ pour lesquelles + Bou, I REVENDICATIONS1. Procédé de caractérisation ultrasonore non invasive d'un milieu hétérogène, le procédé comprenant : - une étape de génération, au moyen d'au moins un premier réseau de transducteurs, d'une série d'ondes ultrasonores incidentes dans une zone dudit milieu hétérogène ; - une étape d'enregistrement d'une matrice de réflexion expérimentale (Roi (0) définie entre une base d'émission (i) en entrée et la base des transducteurs (u) en sortie ; une étape de détermination, à partir de ladite matrice de réflexion expérimentale, d'au moins une première matrice de réflexion (R'(Aoi) ) large bande, définie dans une base de focalisation en entrée et en sortie, dans une première bande spectrale (Acod; une étape de détermination d'une première matrice distorsion (Dor, Drr) définie dans ladite première bande spectrale, entre ladite base de focalisation et une base d'observation, ladite première matrice distorsion résultant, directement ou après changement de base, du produit matriciel terme à terme de ladite première matrice de réflexion large bande (Fter(Aoi)), définie entre ladite base de focalisation et une base de correction des aberrations, avec la matrice conjuguée en phase d'une matrice de réflexion de référence, définie pour un milieu modèle, entre les mêmes bases; - une étape de détermination à partir de ladite première matrice distorsion, d'au moins une cartographie d'un paramètre physique dudit milieu hétérogène. 2. Procédé selon la revendication I, dans lequel la détermination de ladite première matrice de réflexion large bande (R,(Acui) ) comprend : - un changement de base en entrée et en sortie, à partir de ladite matrice de réflexion expérimentale exprimée dans le domaine fréquentiel (Ro' (w))pour former une matrice de réflexion (Rrr(co)), définie dans le domaine fréquentiel, dans une base de focalisation en entrée et en sortie; - une sommation cohérente sur ladite première bande spectrale donnée (Aoi) de ladite matrice de réflexion définie dans une base de focalisation en entrée et en sortie pour former ladite première matrice de réflexion large bande (R'(An I)) définie dans une base de focalisation en entrée et en sortie. Procédé selon l'une quelconque des revendications précédentes, dans lequel la détermination de ladite première matrice distorsion comprend : 36 un changement de base, en entrée ou en sortie, de ladite première matrice de réflexion large bande (Rrr) définie dans une base de focalisation en entrée et en sortie, pour obtenir ladite première matrice de réflexion large bande (Roi.) déterminée entre la base de focalisation et la base de correction des aberrations; 5 le produit matriciel terme à terme de ladite première matrice de réflexion large bande (Rer) définie entre ladite base de focalisation et la base de correction des aberrations, avec la matrice conjuguée en phase de la matrice de réflexion de référence définie entre les mêmes bases. 4 Procédé selon l'une quelconque des revendications précédentes, dans lequel : 10 ladite base d'observation est la base de focalisation ; ladite première matrice distorsion est définie dans une base de focalisation en entrée et en sortie, et est obtenue soit directement à partir de ladite première matrice de réflexion large bande (Rrr) définie dans une base de focalisation en entrée et en sortie, soit par changement de base à partir de la première matrice distorsion définie entre la 15 base de focalisation et la base de correction des aberrations. 5. Procédé selon la revendication 4, dans lequel la détermination d'au moins une cartographie d'un paramètre physique dudit milieu hétérogène comprend la détermination, à partir de ladite première matrice distorsion définie dans une base de focalisation en entrée et en sortie, et pour au moins un premier domaine d'isoplanétisme 20 identifié, d'au moins une première fonction d'étalement (PSF). 6. Procédé selon l'une quelconque des revendications précédentes, dans lequel la détermination d'au moins une cartographie d'un paramètre physique dudit milieu hétérogène comprend : la détermination des invariants dans ladite base de focalisation de ladite première 25 matrice distorsion, afin d'identifier au moins un premier domaine d'isoplanétisme dans ladite base de focalisation ; la détermination pour chaque premier domaine d'isoplanétisme identifié, d'une cartographie d'au moins une première loi d'aberration dans ladite base de correction des aberrations. 30 7 Procédé selon la revendication 6, dans lequel la détermination des invariants de ladite première matrice distorsion comprend une décomposition en valeurs singulières d'au moins l'une des matrices du groupe de matrices comprenant : ladite première matrice 37 distorsion, ladite première matrice distorsion normalisée, une matrice de corrélation normalisée de ladite première matrice distorsion. 8. Procédé l'une quelconque des revendications 6 ou 7, dans lequel la détermination d'au moins une cartographie d'un paramètre physique dudit milieu hétérogène comprend en 5 outre la détermination dans ladite base d'observation d'une première matrice de réflexion large bande corrigée par la ou lesdites première(s) loi(s) d'aberration. 9 Procédé selon la revendication 8, comprenant en outre : la détermination d'une matrice distorsion corrigée, définie entre ladite base de focalisation et ladite base de correction des aberrations par le produit matriciel terme à 10 terme de ladite première matrice de réflexion large bande corrigée avec la matrice conjuguée en phase de ladite matrice de réflexion de référence. 10. Procédé selon la revendication 9, comprenant en outre : la détermination des invariants dans ladite base de focalisation de ladite matrice distorsion corrigée, afin d'identifier au moins un deuxième domaine d'isoplanétisme 15 dans ledit plan focal ; la détermination pour chaque deuxième domaine d'isoplanétisme identifié, d'une cartographie d'une deuxième loi d'aberration dans ladite base de correction. 11 Procédé selon l'une quelconque des revendications précédentes, comprenant en outre : une étape de détermination, à partir de ladite matrice de réflexion expérimentale, d'au 20 moins une deuxième matrice de réflexion large bande (12,.,(Ao)2) ), définie dans une base de focalisation en entrée et en sortie, dans une deuxième bande spectrale (4co2) différente de ladite première bande spectrale ; une étape de détermination d'au moins une deuxième matrice distorsion (Der(A02), D'(Acon)) définie dans ladite deuxième bande spectrale entre ladite base de 25 focalisation et ladite base d'observation, ladite deuxième matrice distorsion résultant, directement ou après changement de base, du produit matriciel terme à terme de ladite deuxième matrice de réflexion large bande (Ror(Aoq)), définie entre ladite base de focalisation et ladite base de correction des aberrations, avec la matrice conjuguée en phase d'une matrice de réflexion de référence, définie pour un milieu modèle, entre les 30 mêmes bases. 12. Procédé selon la revendication 11, comprenant en outre : 38 la détermination des invariants dans ladite base de focalisation de ladite au moins une deuxième matrice distorsion (Dor(Aco,), afin d'identifier au moins un premier domaine d'isoplanétisme dans ledit plan focal ; la détermination pour chaque domaine d'isoplanétisme identifié, d'une cartographie 5 d'au moins une première loi d'aberration dans ladite base de correction ; - la détermination dans ladite base d'observation d'au moins une deuxième matrice de réflexion large bande corrigée par la ou lesdites première(s) loi(s) d'aberration ; la sommation de ladite première matrice de réflexion large bande corrigée et de ladite au moins une deuxième matrice de réflexion large bande corrigée pour former une 10 matrice de réflexion corrigée ultra large bande. 13. Procédé selon l'une quelconque des revendications précédentes, dans lequel la détermination d'au moins une cartographie d'un paramètre physique dudit milieu hétérogène comprend la détermination de la vitesse du son dans la zone insonifiée. 14. Procédé selon l'une quelconque des revendications précédentes, dans lequel la 15 détermination d'au moins une cartographie d'un paramètre physique dudit milieu hétérogène comprend la détermination d'un taux de diffusion multiple dans la zone insonifiée. 15. Procédé selon l'une quelconque des revendications précédentes, comprenant en outre l'identification et/ou la suppression des composantes de réflexion spéculaire et/ou de 20 réflexions multiples à partir de ladite première matrice distorsion, exprimée en entrée et en sortie dans la base d'ondes planes. 16. Système (40) de caractérisation ultrasonique non invasive d'un milieu hétérogène (20) comprenant : Au moins un premier réseau (10) de transducteurs configure pour générer une série 25 d'ondes ultrasonores incidentes dans une zone dudit milieu hétérogène et enregistrer les ondes ultrasonores rétrodiffusées par ladite zone en fonction du temps; une unité de calcul (42), couplée audit au moins un premier réseau de transducteurs, et configurée pour : o enregistrer une matrice de réflexion expérimentale (Ru; (t)) définie entre une 30 base d'émission (i) à l'entrée et la base des transducteurs (u) en sortie ; o déterminer, à partir de ladite matrice de réflexion expérimentale, au moins une première matrice de réflexion large bande (R'(Awi) ), définie dans une base de 39 focalisation en entrée et en sortie, dans une première bande spectrale donnée (A01); o déterminer une première matrice distorsion (Dor, D') définie dans ladite première bande spectrale, entre ladite base de focalisation et une base 5 d'observation, ladite première matrice distorsion résultant, directement ou après changement de base, du produit matriciel terme à terme de ladite première matrice de réflexion large bande (Ror(Ao,)), définie entre ladite base de focalisation et une base de correction des aberrations, avec la matrice conjuguée en phase d'une matrice de réflexion de référence, définie pour un 10 milieu modèle, entre les mêmes bases; o déterminer à partir de ladite première matrice distorsion, au moins une cartographie d'un paramètre physique dudit milieu hétérogène.