Brevet FR2958409A1 — publié le 20111007

Procédé et système d'échographie ultrasonore anharmonique La présente invention concerne un procédé et un système d'exploration d'un milieu par échographie ultrasonore comportant l'émission vers ledit milieu d'au moins un signal ultrasonore par au moins un émetteur ultrasonore et la réception par au moins un récepteur ultrasonore d'au moins un signal d'écho ultrasonore anharmonique, généré lors de la propagation dudit signal ultrasonore dans ledit milieu, et dont la forme générale peut être exprimée par x(t) = xo + x, cos4(t)), où I(t) est la phase dudit signal d'écho. Elle s'applique en particulier au domaine de l'échographie médicale et à l'analyse des signaux d'échographie reçus par la sonde échographique après propagation, diffusion et réflexion de ces signaux dans le milieu exploré. L'échographie est basée sur l'utilisation d'ondes ultrasonores, qui sont des ondes de pression se propageant dans un milieu élastique. Lorsqu'une telle onde rencontre une interface entre deux milieux d'impédances acoustiques différentes, une première partie de son énergie est transmise, et une deuxième partie est réfléchie par l'interface, formant un écho de l'onde d'origine. Ainsi, en échographie, les ondes ultrasonores émises par la sonde se propagent dans les tissus, sont réfléchies par de telles interfaces, puis reçues par la sonde. Une image du tissu balayé par la sonde peut alors être reconstituée sur un écran, en affectant à chaque point de l'écran une luminosité ou un niveau de gris dépendant de l'amplitude de l'écho reçu en provenance d'une zone correspondante du tissu. Le problème général rencontré en échographie ultrasonore consiste le plus souvent à obtenir une image aussi exacte que possible du milieu analysé, et d'une résolution suffisante pour visualiser les détails de ce milieu. La résolution des images obtenues par échographie est d'autant plus fine que la fréquence des ondes réfléchie est élevée. Cependant, l'absorption des ultrasons par les tissus augmente également avec la fréquence des ondes, de telle sorte que l'échographie à haute fréquence est généralement limitée à l'exploration d'organes superficiels. Ainsi, la fréquence des ondes ultrasonores émises dans les tissus est généralement inférieure à 5 MHz, la résolution des images étant ainsi au maximum de l'ordre de quelques microns. Afin d'obtenir des images de meilleure résolution tout en limitant l'absorption des ondes par les tissus, il est connu de recourir à la technique d'échographie harmonique, qui exploite les propriétés non-linéaires de la propagation des ondes ultrasonores. En effet, les tissus biologiques sont des matériaux non-linéaires, qui, lorsqu'ils transmettent une onde ultrasonore, déforment cette onde au cours de sa propagation dans les tissus. Ainsi, une onde ultrasonore émise par la sonde à une fréquence f peut être reçue, après propagation, diffusion et réflexion dans les tissus, sous la forme d'un signal anharmonique, c'est-à-dire non-linéaire, qui peut être décomposé en un signal fondamental, de fréquence f, et plusieurs harmoniques, de fréquences multiples de la fréquence f. L'échographie harmonique consiste ainsi à émettre dans les tissus un signal ultrasonore à une fréquence fondamentale f, et à recueillir le signal réfléchi à une fréquence harmonique, généralement l'harmonique double, de fréquence 2f, le signal fondamental étant éliminé soit par un filtrage, soit par une technique d'inversion de phase des impulsions. Cette technique permet d'obtenir des images de meilleure résolution, tout en limitant l'absorption de l'onde par les tissus, puisque les signaux harmoniques apparaissent lors de la propagation de l'onde, et ne traversent pas les tissus deux fois, de la sonde vers l'interface de réflexion puis de cette interface vers la sonde. L'échographie harmonique permet également d'augmenter le contraste entre les organes à visualiser et la graisse. Cet effet non-linéaire est généralement accentué par l'utilisation d'agents de contrastes ultrasonores, qui comportent des microbulles de gaz encapsulées, de diamètre de l'ordre de quelques microns, injectés par voie intraveineuse. Ces agents de contraste se comportent comme des oscillateurs non-linéaires et amplifient ainsi la déformation de l'onde et la génération d'harmoniques. Cependant, l'échographie harmonique n'exploite que partiellement le signal réfléchi par les tissus, puisque seule l'harmonique double est analysée. De plus, cette harmonique double n'est caractérisée que par sa seule amplitude. Ainsi, la non-linéarité, ou anharmonicité, du signal réfléchi par les tissus n'est mesurée qu'au moyen de cette amplitude, qui n'est pas facilement interprétable en terme de non-linéarité. Par ailleurs, il est connu de décomposer un signal périodique de fréquence f, harmonique ou non-harmonique, au moyen des séries de Fourier, introduites en 1822 par Joseph Fourier. Cette technique consiste à décomposer le signal en une somme infinie de fonctions sinusoïdales de fréquences multiples de f. Tout signal périodique x(t) de période T=1/f peut ainsi être exprimé comme une somme de fonctions sinusoïdales, du type : ;2,r- t x(t)= 1cn(x)e n=~ Les coefficients cn, appelés coefficients de Fourier, qui sont définis par la formule : T 2 en (x) = T J x(t)e i2~cT t T 2 (1) (2) constituent un codage du signal x(t). A partir de cette analyse, la donnée des coefficients en permet de synthétiser le signal x(t). Afin d'augmenter la compacité de ce codage, il convient de limiter le nombre de coefficients cn, à priori infini, et de ne conserver que les premiers termes de l'expression (1). Ces termes doivent cependant être en nombre suffisant pour que le signal synthétisé à partir du codage soit aussi proche que possible du signal x(t) d'origine. La décomposition en séries de Fourier est universelle, en ce qu'elle permet de décomposer tout type de signal périodique. Cependant, lorsque le signal analysé n'est pas harmonique ou quasi harmonique, la décomposition de Fourier nécessite de conserver un grand nombre de coefficients, coefficients auxquels il est difficile de donner un sens physique. Le but de l'invention est donc de permettre une mesure plus précise et plus pertinente de l'anharmonicité du signal non-linéaire généré lors de la propagation d'une onde ultrasonore, en permettant l'analyse et la représentation de ce signal au moyen d'un petit nombre de paramètres, par rapport au nombre de paramètres nécessaires au codage par série de Fourier, lesdits paramètres étant porteurs d'un sens physique et constituant une signature simple et explicite de la dynamique et de la structure interne du milieu qui produit ce signal. A cet effet, l'invention a pour objet un procédé d'exploration d'un milieu par échographie ultrasonore du type précité, caractérisé en ce que ledit procédé comporte l'émission vers ledit milieu d'au moins un signal ultrasonore par au moins un émetteur ultrasonore et la réception par au moins un récepteur ultrasonore d'au moins un signal d'écho ultrasonore anharmonique, généré lors de la propagation dudit signal ultrasonore dans ledit milieu, et dont la forme générale peut être exprimée par x(t) = xo + x, cos4(t)), où I(t) est la phase dudit signal d'écho, caractérisé en ce qu'il comporte en outre les étapes suivantes : - analyse dudit signal d'écho, comportant les étapes de détermination d'une expression d'une équation de phase F((13) = , et de détermination d'une expression de la phase I(t) en fonction de paramètres mesurant l'anharmonicité dudit signal d'écho et sa morphologie, à partir des fonctions pcosn et psinn définies par : k pcosä(t,r)= Ecos(kt)kn et psinä(t,r)= Esin(kt) rk , k=1 k=1 - génération d'au moins une image permettant de visualiser la valeur d'au moins un desdits paramètres dudit signal d'écho ou une combinaison de ces paramètres en fonction de sa provenance dans ledit milieu. dt 1ùr2 dans laquelle r, variant dans [0,1], est un paramètre mesurant l'anharmonicité dudit signal d'écho ; - l'onde x(t) est exprimée au moyen de deux paramètres r et c Do, sous la forme : x(t) = xo + a1 h sin(t, r)+ b1 hcos(t, r) où a1 = x1 cos»o) et b1 = ùx1 sin(cDo ), les fonctions hsin et hcos étant définies par : hcos : (t, r) (1 + r2)cos(t)+ 2r et hsin : (t, r) (1û r2)sin(t) 1 + r2 û 2rcos(t) 1 + r2 û 2rcos(t) - l'équation de phase est exprimée sous la forme : F(d)) = P(~) , Q(iv) dans laquelle P(cD) et Q(cl)) sont des polynômes trigonométriques ; - l'expression de la phase t(t) est déterminée sous la forme : n t( )= +Eakpsin»ùpk,rk)ùbkpcos»ùPk,rk) k=1 dans laquelle les fonctions psin1 et pcos1 sont définis par : pcos1(t,r)=Ecos(kt)-k et psin1(t,r)=Esin(kt)-k , k=1 k k=1 k - l'émission vers ledit milieu d'au moins un signal ultrasonore comporte la focalisation dudit signal ultrasonore sur une zone cible dudit milieu. Ainsi réalisé, le procédé selon l'invention permet d'analyser l'ensemble de l'onde non-linéaire reçue par la sonde, et de caractériser cette onde au moyen d'un petit nombre de paramètres, par rapport aux procédés de décomposition selon l'état de la technique. De plus, ces paramètres ont un sens physique, ce qui permet de comparer plusieurs ondes en comparant directement les paramètres issus de leur analyse. Selon un autre aspect, l'invention a également pour objet un système d'exploration d'un milieu par échographie ultrasonore comportant au moins un émetteur ultrasonore 4 Le procédé selon l'invention comporte également les caractéristiques suivantes, prises séparément ou en combinaison : - l'équation de phase est exprimée sous la forme : dcl) 1+r2+2rcos(4) apte à émettre vers ledit milieu au moins un signal ultrasonore, au moins un récepteur ultrasonore, apte à recevoir au moins un signal d'écho ultrasonore anharmonique, généré lors de la propagation dudit signal ultrasonore dans ledit milieu, et dont la forme générale peut être exprimée par x(t) = xo + x, cos(d)(t)), où t(t) est la phase dudit signal d'écho, caractérisé en ce qu'il comporte en outre : - des moyens pour analyser ledit signal d'écho, comportant des moyens pour déterminer une expression d'une équation de phase F(.:13)= d , et des moyens pour déterminer une expression de la phase t(t) en fonction de paramètres mesurant l'anharmonicité dudit signal d'écho et sa morphologie, à partir des fonctions pcosn et psinn définies par : k pcosä(t, r) = Ecos(kt) kn et psinä(t, r) = E sin(kt) rk , k=1 k=1 - des moyens pour générer au moins une image permettant de visualiser la valeur d'au moins un desdits paramètres ou toute combinaison de ces paramètres dudit signal d'écho en fonction de sa provenance dans ledit milieu. - il comporte des moyens pour exprimer l'onde x(t) au moyen de deux paramètres r et cDo, sous la forme : x(t) = xo + a, h sin(t, r)+ b, hcos(t, r) où a, = x, cos4o) et b, = ùx, sin(d3o ), les fonctions hsin et hcos étant définies par : hcos : (t, r) (1 + r2)cos(t)+ 2r et hsin : (t, r) (1ù r2 )sin(t) 1 + r2 ù 2rcos(t) 1 + r2 ù 2rcos(t) - il comporte des moyens pour exprimer l'équation de phase sous la forme : F(d)) = P(~) , Q(iv) Le système selon l'invention comporte également les caractéristiques suivantes, prises séparément ou en combinaison : - il comporte des moyens pour exprimer l'équation de phase sous la forme : cet:. 1+r2 +2rcos((D ) dt 1ùr2 20 dans laquelle r, variant dans [0,1], est un paramètre mesurant l'anharmonicité dudit signal d'écho ; 25 30 6 dans laquelle P(cD) et Q(t) sont des polynômes trigonométriques ; - il comporte des moyens pour exprimer la phase t(t) sous la forme : n t( )= cl) +Eakpsin1(cl) ûpk,rk)ûbkpcoslûPk,rk) k=1 dans laquelle les fonctions psin1 et pcos1 sont définis par : pcos1(t,r)=Ecos(kt)-k et psin1(t,r)=Esin(kt)-k , k=1 k k=1 k - il comporte des moyens pour focaliser ledit signal ultrasonore émis sur une zone cible dudit milieu. L'invention sera mieux comprise au regard d'un exemple de réalisation de l'invention qui va maintenant être décrit en faisant référence aux figures annexées parmi lesquelles : la figure 1 est un schéma représentant un système d'échographie ultrasonore selon un mode de réalisation de l'invention ; et - la figure 2 est un schéma synoptique illustrant le procédé d'échographie ultrasonore selon un mode de réalisation de l'invention. Le système 1 d'échographie ultrasonore représenté sur la figure 1 comprend une sonde échographique 3, reliée à une unité 5 de commande et d'analyse, elle-même reliée à un écran 6. La sonde échographique 3 comporte une pluralité de transducteurs piézoélectriques, et elle est apte à émettre et à recevoir des signaux ultrasonores. L'unité 5 de commande et d'analyse est apte à commander l'émission par la sonde échographique de signaux ultrasonores, à analyser les signaux ultrasonores reçus par la sonde, à générer à partir de ces signaux des images représentatives des tissus biologiques 7 analysés, et à transmettre ces images à I `écran 6. La figure 2 est un schéma synoptique illustrant le procédé d'échographie ultrasonore selon un mode de réalisation de l'invention. Lors d'une étape 10, l'unité 5 de commande et d'analyse commande l'émission par la sonde 3 d'un signal ultrasonore périodique linéaire, à une fréquence f, et pendant une durée de l'ordre 1 à 2 millisecondes, vers les tissus 7 à analyser, éventuellement dans lesquels ont préalablement été injectés des agents de contraste ultrasonore. Cette étape d'émission comprend une étape de focalisation du signal sur une zone de taille restreinte du milieu analysé, située par exemple à une interface entre deux milieux d'impédances acoustiques différentes. Avantageusement, cette focalisation est réalisée au moyen d'une technique de focalisation électronique, ou d'une technique de focalisation par renversement temporel. Lors d'une étape 14, le signal ultrasonore se propage dans les tissus 7, et des échos sont engendrés par réflexion sur la zone cible et par diffusion du signal par les tissus, en direction de la sonde 3. Lors de la propagation du signal, celui-ci subit des déformations, de la part des tissus biologiques et des microbulles de l'agent de contraste, de telle sorte que le signal que la sonde reçoit est un signal anharmonique, c'est-à-dire non-linéaire. La sonde reçoit donc en 16, après un délai td suivant l'émission 10, ce signal anharmonique, et produit un signal électrique anharmonique x(t), envoyé vers l'unité 5 de commande et d'analyse. Puis, lors d'une étape 18, l'unité 5 analyse le signal périodique anharmonique x(t) reçu, de manière à en extraire des paramètres caractéristiques de la zone cible. Tout signal périodique simple, c'est-à-dire possédant un maximum et un minimum par période, peut être décrit sous la forme suivante : x(t) = xo + x, cos(d)(t)) (3) dans laquelle toute la dépendance temporelle est contenue dans la fonction de phase cl). Or, dans un signal périodique anharmonique, la principale contribution à l'anharmonicité provient de la brisure de symétrie de la dynamique de phase. Ainsi, toute l'information dynamique pertinente est exprimée par la dynamique de phase. Lors de l'analyse du signal x(t), il convient donc d'étudier cette phase 4)(t), et en particulier la dynamique de phase exprimée par la fonction F, dérivée de la fonction cl) par rapport au temps t : Ainsi, la morphologie du signal est complètement déterminée par la connaissance 25 de F. L'étape d'analyse 18 du procédé selon l'invention consiste donc à décrire cette fonction F au moyen d'un très petit nombre de paramètres. On entendra par petit nombre de paramètres un nombre de paramètres réduit par rapport au nombre de paramètres nécessaire à la décomposition de la même fonction, au moyen des séries de Fourier, 30 avec un niveau de précision équivalent. Cette étape d'analyse 18 comprend ainsi une première étape consistant à exprimer la phase cl), et en particulier la fonction F, dérivée de cl) par rapport au temps. dt (4) Dans le cas le plus simple, et pour un signal de période 27t, la dynamique de phase peut être écrite sous la forme : F(dcl) - = 1 + r2 + 2rcos(4) dt -r 2 appelée équation de phase. La fonction F présente dans ce cas une symétrie de réflexion par rapport à l'axe 43=0. Cette expression de la dynamique de phase ne contient qu'un seul paramètre, r, qui varie dans l'intervalle [0,1]. La limite r=0 correspond à un signal harmonique, la limite r=1 à un signal infiniment anharmonique. Le signal x(t), qui peut être écrit sous la forme : x(t) = xo + x, cos(d)(t, r)ù ) (8) où cl)o est une origine de phase, est décomposé et réécrit sous une forme faisant intervenir les paramètres r et cl)o : x(t) = xo + alh sin(t, r)+ b,hcos(t, r) (9) avec a, = x, cos4o) et b, = ùx, sin(d3o ), et dans laquelle on a défini les fonctions hcos et hsin suivantes : hcos (t, r) (1 + r2 )cos(t)+ 2r (10) : 1 + r2 ù 2rcos(t) hsin (t, r) (1ù r2)sin(t) (11) : 1+ r2 ù 2rcos(t) Ainsi, la décomposition du signal x(t) fait intervenir seulement deux paramètres, r et 4 . r, appelé paramètre d'anharmonicité, mesure le degré d'anharmonicité du signal, la limite r=0 correspond à un signal harmonique, la limite r=1 à un signal infiniment anharmonique. Par ailleurs, le paramètre cDo, qui définit la composition du signal dans les deux fonctions hcos et hsin, est un paramètre de morphologie, qui correspond à l'angle de symétrie de réflexion de la dynamique de phase. Dans le cas général, c'est-à-dire pour un signal périodique quelconque, l'équation de phase peut s'écrire sous la forme : F(`1,) - Pn (`:) (6) Qm (d') dans laquelle Pn et Qm sont des polynômes trigonométriques de degrés respectifs n et m. La forme générale d'un polynôme trigonométrique de degré n est : (5) n Pn((3)= ao +Enak cos(k(l))+bk sin(k(l)) k=1 L'analyse du signal x(t) consiste alors à déterminer une expression de cl) faisant intervenir un petit nombre de paramètres, ce qui permet de déterminer une expression du signal x(t) en fonction de ces paramètres. Avantageusement, l'équation de phase (6) peut être réécrite sous la forme : 1 dt 0,7,4) (12) F(^:13) el) Pn ((D) La factorisation du polynôme Pn(t) permet de transformer F( iv) en une somme de termes simples, ce qui permet de réécrire l'équation de phase sous la forme : dt ù a + ak cos(cl) ù pk )+ bk sin(cl) + pk) L (13) c o s ( c l ) dcl) dans laquelle les paramètres rk, compris entre 0 et 1, mesurent l'anharmonicité du signal x(t), et les paramètres pkcaractérisent sa morphologie. La période T du signal peut être déterminée en intégrant cette équation par rapport à4,entre 0et2n: =2z de r a T= f =22r ao+L k k F(~) k 1ù rk A partir de ce résultat, et des contraintes selon lesquelles la période est égale à 27t et le signal est harmonique lorsque les coefficients rk sont tous nuls, l'équation de phase peut être exprimée ainsi : (7) (14) dt =1+EDk(cl)ùpk) dl) k=1 n (15) où la fonction Dk est définie par : rk(akcos(c)+bksin(c))ùak) Dk cl) (1 + rk ù 2rk cos()) et vérifie : f Dk(cl-))dcl-) =0 ,13=0 La définition des fonctions des fonctions polycos et polysin, notées pcosn et psinn, qui s'expriment par : k pcosn (t, r) = E cos(kt) - n k=1 k (16) (17) (18) k psinn(t,r)= Esin(kt) kn k=1 et possèdent entre autres les propriétés suivantes : coso (t, r) = r(cos(t) ù r) p 1 + 1-2 ù 2r cos(t) psino (t, r) = r sin(t) 1 + r2 ù 2rcos(t) pcosl(t, r) = ù 2 In(1 + r2 ù 2r cos(t)) psinl(t, r) = tan-1 r sin(t) ~1ù rcos(t), permet de réécrire l'équation de phase sous la forme : dt =1+Eakpcos0(Dù pk,rk)+bkpsin0(Dùpk,rk) (24) k=1 La résolution de cette équation permet d'accéder à une expression analytique de t(1), qui s'exprime par : n t( )= cl) +Enakpsinlùpk,rk)ùbkpcosl(cl) ù pk,rk) (25) k=1 Le temps t est donc exprimé en fonction de la phase cl), et de manière duale la phase cl) est exprimée en fonction du temps t, à l'aide de paramètres indépendants clairement définis, qui mesurent l'anharmonicité (paramètres r ou rk), et la morphologie (paramètres 4 ou Pk). Ainsi, lors de l'étape 18 d'analyse, le processeur code le signal x(t) au moyen d'un petit nombre de paramètres. Selon un mode de réalisation, le signal anharmonique x(t) est décrit de manière quasi-exacte par une période, une amplitude, une harmonicité r et une morphologie c Do. Selon un autre mode de réalisation, le signal anharmonique x(t) est décrit de manière encore plus précise par deux couples de paramètres (ri, p1) et (r2, P2), complétés de leurs poids respectifs. La zone cible, d'où provient le signal réfléchi x(t), est donc caractérisée par un nombre restreint de paramètres, notamment par le délai td, caractéristique de la profondeur de la zone cible dans les tissus, et par les paramètres d'harmonicité et de morphologie, qui permettent de caractériser l'interface de réflexion. L'étape 10 d'émission d'une onde ultrasonore est réitérée, en modifiant la zone cible de convergence de cette onde, avec une fréquence de l'ordre de plusieurs milliers d'émissions par seconde. Ainsi, à l'issue de nouvelles étapes 14, 16 et 18, chacune des (19) (20) (21) (22) (23) zones cibles du tissu analysé est caractérisée par ces paramètres d'harmonicité et de morphologie. Le processeur génère alors dans une étape 20 plusieurs images, représentant chacune une cartographie du tissu analysé pour un paramètre donné, et ces images sont affichées par l'écran 6. Par exemple, selon un mode de réalisation, une première image représente, en niveau de gris, la valeur du paramètre r, comprise entre 0 et 1, en chacune des zones analysées, une telle image permettant de visualiser le comportement non-linéaire des tissus analysés, et une deuxième image représente la valeur du paramètre cDo, représentatif de la morphologie de la réponse de ces tissus. Le procédé selon l'invention permet ainsi d'analyser l'ensemble des signaux d'échographie, et d'extraire de ces signaux un nombre restreint de paramètres, permettant une représentation compacte et pertinente des tissus analysés. Il devra toutefois être compris que l'exemple de réalisation présenté ci-dessus n'est pas limitatif. Notamment, les signaux ultrasonores peuvent être émis par la sonde sans focalisation. Bien entendu, d'autres modes de réalisation encore peuvent être envisagés. REVENDICATIONS1.- Procédé d'exploration d'un milieu (7) par échographie ultrasonore comportant une émission (10) vers ledit milieu d'au moins un signal ultrasonore par au moins un émetteur ultrasonore (3) et la réception (16) par au moins un récepteur ultrasonore (3) d'au moins un signal d'écho ultrasonore anharmonique, généré lors de la propagation (14) dudit signal ultrasonore dans ledit milieu (7), et dont la forme générale peut être exprimée par x(t) = xo + x, cos(cl)(t)), où cl:^(t) est la phase dudit signal d'écho, caractérisé en ce qu'il comporte en outre les étapes suivantes : - analyse (18) dudit signal d'écho, comportant les étapes de détermination d'une expression d'une équation de phase F(4) _ - , et de détermination d'une expression de la phase cl:^(t) en fonction de paramètres (r, rk, 4 , Pk) mesurant l'anharmonicité dudit signal d'écho et sa morphologie, à partir des fonctions pcosn et psinn définies par : k pcosä(t, r) _ Ecos(kt) kn et psinä(t, r) _ E sin(kt) rk , k=1 k=1 - génération (20) d'au moins une image permettant de visualiser la valeur d'au moins un desdits paramètres (r, rk, 4 , Pk) dudit signal d'écho ou une combinaison de ces paramètres en fonction de sa provenance dans ledit milieu. 2.- Procédé d'exploration d'un milieu par échographie ultrasonore selon la revendication 1, caractérisé en ce que l'équation de phase est exprimée sous la forme : dl) 1+r2+2rcos(4) dt 1ùr2 dans laquelle r, variant dans [0,1], est un paramètre mesurant l'anharmonicité dudit signal d'écho. 3.- Procédé d'exploration d'un milieu par échographie ultrasonore selon la revendication 2, caractérisé en ce que l'onde x(t) est exprimé au moyen de deux paramètres r et cDo, sous la forme : x(t) = xo + a, h sin(t, r)+ b, hcos(t, r) où a, = x, cos(cDo) et b, = ûx, sin(cDo ), les fonctions hsin et hcos étant définies par : hcos : (t, r) (1 + r2)cos(t)+ 2r et hsin : (t, r) (1û r2)sin(t) 1 + r2 û 2rcos(t) 1 + r2 û 2rcos(t) 4.- Procédé d'exploration d'un milieu par échographie ultrasonore selon la revendication 1, caractérisé en ce que l'équation de phase est exprimée sous la forme : F(cl)) = P(13) 5.- Procédé d'exploration d'un milieu par échographie ultrasonore selon la revendication 4, caractérisé en ce que l'expression de la phase cl:^(t) est déterminée sous la forme : n t()= .13+ ùpk,rk)ùbkpcosl(cl) ùpk,rk) k=1 dans laquelle les fonctions psin1 et pcos1 sont définis par : 10 pcos1(t,r)=Ecos(kt)rk et psin1(t,r)='sin(kt)rk k=1 k k=1 k 6.- Procédé d'exploration d'un milieu par échographie ultrasonore selon l'une des revendications précédentes, caractérisé en ce que l'émission (10) vers ledit milieu d'au moins un signal ultrasonore comporte la focalisation dudit signal ultrasonore sur une zone 15 cible dudit milieu (7). 7.- Système d'exploration d'un milieu (7) par échographie ultrasonore comportant au moins un émetteur (3) ultrasonore apte à émettre vers ledit milieu (7) au moins un signal ultrasonore, au moins un récepteur (3) ultrasonore, apte à recevoir au moins un signal 20 d'écho ultrasonore anharmonique, généré lors de la propagation dudit signal ultrasonore dans ledit milieu (7), et dont la forme générale peut être exprimée par x(t) = xo + x1 cos(cl)(t)), où cl:^(t) est la phase dudit signal d'écho, caractérisé en ce qu'il comporte en outre : - des moyens (5) pour analyser ledit signal d'écho, comportant des moyens pour 25 déterminer une expression d'une équation de phase F(4) = , et des moyens pour déterminer une expression de la phase cl:^(t) en fonction de paramètres (r, rk, 4 , Pk) mesurant l'anharmonicité dudit signal d'écho et sa morphologie, à partir des fonctions pcosn et psinn définies par : pcosn(t,r)= Ecos(kt) kn et psinn(t,r)= Esin(kt) kn , k=1 k=1 Q(cD) , dans laquelle P(cl:^) et Q(D) sont des polynômes trigonométriques.5 14 - des moyens (5) pour générer au moins une image permettant de visualiser la valeur d'au moins un desdits paramètres (r, rk, 4 , pk) ou une combinaison de ces paramètres dudit signal d'écho en fonction de sa provenance dans ledit milieu. 8.- Système d'exploration d'un milieu (7) par échographie ultrasonore selon la revendication 7, caractérisé en ce qu'il comporte des moyens (5) pour exprimer l'équation de phase sous la forme : 1+r2+2rcos(4) dt 1ùr2 dans laquelle r, variant dans [0,1], est un paramètre mesurant l'anharmonicité dudit signal 10 d'écho. 9.- Système d'exploration d'un milieu (7) par échographie ultrasonore selon la revendication 8, caractérisé en ce qu'il comporte des moyens (5) pour exprimer l'onde x(t) au moyen de deux paramètres r et cDo, sous la forme : x(t) = xo + a1 h sin(t, r)+ b1 hcos(t, r) où a1 = x1 cos(cDo) et b1 = -x1 sin(cDo ), les fonctions hsin et hcos étant définies par : hcos : (t, r) (1 + r2)cos(t)+ 2r et hsin : (t, r) (1- r2)sin(t) 1 + r2 - 2rcos(t) 1 + r2 - 2rcos(t) 10.- Système d'exploration d'un milieu (7) par échographie ultrasonore selon la 20 revendication 7, caractérisé en ce qu'il comporte des moyens (5) pour exprimer l'équation de phase sous la forme : F(D) = «cl)) , Q(cD) dans laquelle P(cl:^) et Q(') sont des polynômes trigonométriques. 25 11.- Système d'exploration d'un milieu (7) par échographie ultrasonore selon la revendication 10, caractérisé en ce qu'il comporte des moyens (5) pour exprimer la phase cl:^(t) sous la forme : n t()= .13+ k=1 dans laquelle les fonctions psin1 et pcos1 sont définis par : 30 pcos1(t, r) = E- cos(kt)rk et psin1(t, r) = E sin(kt)-k k=1 k k=1 k 15 12.- Système d'exploration d'un milieu (7) par échographie ultrasonore selon l'une des revendications 7 à 11, caractérisé en ce qu'il comporte des moyens pour focaliser ledit signal ultrasonore émis sur une zone cible dudit milieu (7).5