FR 2844058 A1 20040305 FR 0210838 A 20020902 Procédé et dispositif d'imagerie utilisant des ondes de cisaillement. La présente invention est relative aux procédés et 5 dispositifs d'imagerie utilisant des ondes de cisaillement. Plus particulièrement, l'invention concerne un procédé d'imagerie utilisant des ondes de cisaillement pour observer un milieu viscoélastique diffusant qui contient des particules réfléchissant les ondes ultrasonores de com10 pression, procédé comprenant: (a) une étape d'excitation au cours de laquelle on génère une onde élastique de cisaillement dans le milieu viscoélastique, (b) une étape d'observation au cours de laquelle on 15 observe la propagation de l'onde de cisaillement simultanément en une multitude de points d'un champ d'observation dans le milieu viscoélastique, cette étape d'observation comprenant les sous-étapes suivantes: (bl) faire émettre dans le milieu viscoélasti20 que, par un réseau de transducteurs commandés indépendamment les uns des autres, une succession de tirs d'ondes ultrasonores de compression non focalisées à une cadence d'au moins 500 tirs par seconde, (b2) faire détecter et enregistrer en temps réel 25 des signaux acoustiques reçus du milieu viscoélastique 2, comprenant les échos générés par les ondes ultrasonores de compression non focalisées en interagissant avec les particules réfléchissantes dudit milieu viscoélastique, (c) et au moins une étape de traitement au cours de 30 laquelle: (cl) on traite les signaux acoustiques successifs reçus du milieu viscoélastique au cours de la sousétape (b2) pour déterminer des images de propagation successives de l'onde de cisaillement, (c2) et pour déterminer au moins un paramètre de mouvement du milieu viscoélastique en différents points du champ d'observation. On obtient ainsi un "film" illustrant clairement la 5 propagation de l'onde de cisaillement dans le milieu viscoélastique, qui peut permettre une analyse qualitative et/ou quantitative pour repérer des zones de dureté différente du reste du milieu viscoélastique ou des zones ayant un temps de relaxation différent du reste du milieu viscoé10 lastique. Le document WO-A-00/55 616 décrit un exemple d'un tel procédé, dans lequel on génère l'onde de cisaillement à la surface du milieu viscoélastique. Ce procédé donne toute satisfaction notamment pour imager des zones situées 15 relativement près de la surface du milieu viscoélastique. Mais ce procédé connu ne permet pas d'observer certaines zones du milieu viscoélastique, notamment: - des zones suffisamment profondes pour ne pas pouvoir être atteintes par des ondes de cisaillement géné20 rées en surface (les ondes de cisaillement s'atténuent rapidement), - et des zones d'ombre masquées par des obstacles (notamment des parties du squelette d'un patient ou des zones liquides telles que des kystes liquides) qui gênent la 25 propagation des ondes de cisaillement. De plus, si le champ d'observation est partiellement dans une zone d'ombre, il peut être nécessaire de déplacer le dispositif générateur des ondes de cisaillement en cours d'observation, ce qui est fastidieux pour l'utili30 sateur. Enfin, le dispositif générateur des ondes de cisaillement est relativement lourd et complique le dispositif. La présente invention a notamment pour but de pal- lier ces inconvénients. A cet effet, selon l'invention, un procédé du genre en question est caractérisé en ce qu'au cours de l'étape d'excitation (a), on fait générer l'onde élastique de ci5 saillement en faisant émettre au moins une onde ultrasonore focalisée dans le milieu viscoélastique par ledit réseau de transducteurs, la focalisation et la chronologie de ladite onde ultrasonore focalisée, ainsi que la chronologie desdites ondes ultrasonores non focalisées, étant adaptées pour 10 qu'au moins certaines desdites ondes ultrasonores non focalisées parviennent dans le champ d'observation lors de la propagation de l'onde de cisaillement dans ce champ d'observation. Ainsi, le même réseau de transducteurs permet à la 15 fois d'engendrer l'onde élastique de cisaillement de façon choisie dans le champ d'observation, et d'observer ensuite cette propagation, grâce au fait que le dispositif d'imagerie est adapté pour générer soit des ondes ultrasonores focalisées permettant de générer l'onde élastique de 20 cisaillement, soit des ondes ultrasonores non focalisées permettant d'observer la propagation de l'onde de cisaillement, et grâce au choix judicieux: - de la chronologie de ces différentes émissions, - et du ou des points de focalisation de l'onde 25 ultrasonore focalisée. Le procédé d'imagerie selon l'invention est donc aisé à mettre en oeuvre pour un utilisateur, avec un dispositif relativement simple et léger. L'invention est donc d'un très faible cot de revient en comparaison de techni30 ques concurrentes telles que 1'IRM, et elle permet le cas échéant de constituer des systèmes d'imagerie ambulatoires utilisables tant en imagerie préopératoire qu'en imagerie postopératoire, voire en imagerie per-opératoire. A titre d'exemple, dans des applications médicales, le procédé selon l'invention peut permettre de repérer efficacement des zones cancéreuses dans les tissus d'un patient. La propagation des ondes de cisaillement s'y déroule en effet très différemment des zones voisines. Ce repérage 5 s'effectue beaucoup plus facilement que par une observation classique par simple échographie ultrasonore, puisque la propagation des ondes de cisaillement est fonction du module de cisaillement du milieu, lui-même très variable entre une zone de tissus sains et une zone de tissus cancé10 reux: le module de cisaillement varie typiquement dans un rapport de 1 à 30 entre une zone saine et une zone cancéreuse, alors que le module de compression, qui régit la propagation des ondes acoustiques de compression utilisées dans l'échographie ultrasonore, varie seulement de l'ordre 15 de 5% entre un tissu sain et un tissu cancéreux. De même, on peut ainsi repérer des zones nécrosées au sein d'un tissu, par exemple des zones tumorales ayant subi un traitement par hyperthermie aux ultrasons, notamment pour évaluer l'efficacité de ce traitement par hyper20 thermie. Une autre application possible de l'invention concerne l'évaluation quantitative du degré de fibrose du foie, qui est un paramètre important dans les pathologies du foie, notamment l'hépatite C. On notera que l'invention permet de générer l'onde de cisaillement, et d'observer sa propagation, y compris à travers une zone liquide ou une barrière osseuse masquant totalement ou partiellement le champ d'observation (crâne, cage thoracique, etc.), puisqu'il est possible de focaliser 30 des ondes ultrasonores au travers de telles barrières (voir notamment le document WO-A-02/32316 ou la demande de brevet français n 02 10682 du 28 aot 2002). Dans des modes de réalisation préférés du procédé selon l'invention, on peut éventuellement avoir recours en outre à l'une et/ou à l'autre des dispositions suivantes: - au cours de la sous-étape (b2), pour déterminer ledit paramètre de mouvement, on compare plusieurs images successives de propagation (par exemple par corrélation, 5 Doppler, etc.) avec une même image de référence du milieu viscoélastique, cette image de référence était déterminée en tirant au moins une onde ultrasonore de compression non focalisée dans ledit milieu viscoélastique puis en détectant et en enregistrant des échos générés par ladite onde 10 ultrasonore de compression non focalisée lorsqu'elle interagit avec les particules réfléchissantes du milieu viscoélastique (on améliore ainsi la précision de la mesure du paramètre de mouvement (par exemple le déplacement) du milieu viscoélastique, notamment pour de faibles amplitudes 15 de mouvement (typiquement moins de 30 Mm avec la technique d'excitation par ondes acoustiques utilisée ici); - l'étape (a) est précédée d'une étape d'observation initiale (aO) au cours de laquelle on tire au moins une onde ultrasonore de compression non focalisée puis on 20 détecte et on enregistre des échos générés par ladite onde ultrasonore de compression non focalisée en interagissant avec les particules réfléchissantes du milieu viscoélastique, ces échos correspondant (directement ou indirectement) à une image initiale du milieu viscoélastique, et au cours 25 de la sous-étape (b2) ladite image initiale constitue ladite image de référence pour traiter au moins certaines des images successives de déplacement; - au cours de l'étape d'observation initiale (aO), on tire successivement plusieurs ondes ultrasonores de com30 pression non focalisées puis on détecte et on enregistre des échos générés par chaque onde ultrasonore de compression non focalisée en interagissant avec les particules réfléchissantes du milieu viscoélastique, ces échos correspondant (directement ou indirectement) à plusieurs images préliminaires successives du milieu viscoélastique, et on détermine ladite image initiale du milieu viscoélastique en combinant lesdites images préliminaires successives; - ledit paramètre de mouvement est un déplacement du milieu viscoélastique; - l'onde ultrasonore focalisée émise au cours de l'étape d'excitation (a) présente une fréquence f comprise entre 0,5 et 15 MHz et est émise pendant une durée k/f en secondes, o k est un entier compris entre 50 et 5000 et f 10 est exprimé en Hz; - l'onde ultrasonore focalisée émise au cours de l'étape d'excitation (a) présente une fréquence comprise entre 0,5 et 15 MHz et est émise pendant une succession de périodes d'émission séparées par des périodes de repos, les 15 périodes d'émission se succédant à une cadence comprise entre 10 et 1000 émissions par seconde; - l'onde ultrasonore focalisée émise au cours de l'étape d'excitation (a) est une combinaison linéaire (notamment une somme) de deux signaux monochromatiques de fré20 quences respectives fi et f2 telles que 20 Hz émettant des ondes ultrasonores différentes au cours des étapes d'excitation (a) successives, puis on combine les différentes cartographies obtenues au cours des étapes de 10 cartographie (d) successives, pour calculer une cartographie synthétique du champ d'observation. Par ailleurs, l'invention a également pour objet un dispositif d'imagerie utilisant des ondes de cisaillement pour observer un milieu viscoélastique diffusant qui 15 contient des particules réfléchissant les ondes ultrasonores de compression, ce dispositif comprenant un réseau de transducteurs commandés indépendamment les uns des autres par au moins une unité centrale électronique adaptée pour: - faire générer au moins une onde élastique de ci20 saillement dans le milieu viscoélastique, observer la propagation de l'onde de cisaillement simultanément en une multitude de points d'un champ d'observation dans le milieu viscoélastique, en faisant émettre dans le milieu viscoélastique, par ledit réseau de 25 transducteurs, une succession de tirs d'ondes ultrasonores de compression non focalisées, à une cadence d'au moins 500 tirs par seconde, puis en faisant détecter en temps réel par ledit réseau de transducteurs, et pour enregistrer en temps réel, des signaux acoustiques reçus du milieu viscoé30 lastique, comprenant les échos générés par les ondes ultrasonores de compression non focalisées en interagissant avec les particules réfléchissantes dudit milieu viscoélastique, - et traiter les signaux acoustiques successifs reçus du milieu viscoélastique pour déterminer des images de propagation successives de l'onde de cisaillement, puis déterminer au moins un paramètre de mouvement du milieu viscoélastique en différents points du champ d'observation, caractérisé en ce que l'unité centrale électronique est 5 adaptée pour faire générer l'onde élastique de cisaillement en faisant émettre au moins une onde ultrasonore focalisée dans le milieu viscoélastique par ledit réseau de transducteurs, la focalisation et la chronologie de ladite onde ultrasonore focalisée, ainsi que la chronologie desdites on10 des ultrasonores non focalisées, étant adaptées pour que lesdites ondes ultrasonores non focalisées parviennent dans le champ d'observation lors de la propagation de l'onde de cisaillement dans ce champ d'observation. D'autres caractéristiques et avantages de l'inven15 tion apparaîtront au cours de la description suivante d'une de ses formes de réalisation, donnée à titre d'exemple non limitatif, en regard du dessin joint. Sur le dessin, la figure 1 est une vue schématique d'un dispositif d'imagerie par ondes de cisaillement selon 20 une forme de réalisation de l'invention. Le dispositif d'imagerie 1 représenté sur la figure 1 est destiné à étudier la propagation des ondes élastiques de cisaillement dans un milieu viscoélastique 2 qui est diffusant vis à vis des ondes ultrasonores de compression, 25 et qui peut être par exemple: - un corps inerte, notamment dans le cas du contrôle de qualité pour des applications industrielles, ou un corps vivant, par exemple une partie du corps d'un patient, dans le cas des applications médicales. 30 Ces mouvements sont suivis par exemple au moyen d'un micro-ordinateur 4 (comprenant au moins une interface d'entrée 4a telle qu'Tun clavier ou autre, et une interface de sortie telle qu'un écran ou autre) ou toute autre unité centrale électronique, qui fait envoyer dans le milieu 2, à partir de sa surface extérieure 3, des ondes ultrasonores de compression qui interagissent avec les particules diffusantes 5 contenues dans le milieu 1, lesquelles particules sont réfléchissantes pour les ondes ultrasonores de com5 pression. Les particules 5 peuvent être constituées par toute hétérogénéité du milieu 1, et notamment, lorsqu'il s'agit d'une application médicale, par des particules de collagène présentes dans les tissus humains (ces particules forment sur les images échographiques des points connus 10 sous le terme "speckle"). Pour observer la propagation de l'onde de cisaillement, on utilise une sonde ultrasonore 6 disposée contre la surface extérieure 3 du milieu observé 1. Cette sonde envoie, selon un axe X, des impulsions d'ondes ultrasonores 15 de compression du type de celles couramment utilisées en échographie, à une fréquence comprise par exemple entre 0,5 et 100 MHz et de préférence entre 0,5 et 15 MHz, par exemple de l'ordre de 4 MHz. La sonde ultrasonore 6 est constituée par un réseau 20 de n transducteurs ultrasonores Tl, T2,..., Ti,..., Tn, n étant un nombre entier n au moins égal à 1. Cette sonde 6 peut se présenter par exemple sous la forme d'une barrette linéaire pouvant comprendre par exemple n = 128 transducteurs alignés selon un axe Y perpendi25 culaire à l'axe X. Mais la sonde en question pourrait être également un réseau bidimensionnel (plan ou non) de transducteurs. Les transducteurs T1, T2,... Tn sont commandés indépendamment les uns des autres par le micro-ordinateur 4, 30 éventuellement par l'intermédiaire d'une unité centrale CPU qui est contenue par exemple dans une baie électronique 7 reliée par un câble souple à la sonde 6. Les transducteurs T1-Tn peuvent ainsi émettre sélectivement: - soit une onde ultrasonore de compression "plane" (c'est à dire en l'occurrence une onde dont le front d'onde est rectiligne dans le plan X, Y) ou tout autre type d'onde non focalisée éclairant l'ensemble du champ d'observation dans le milieu 2, par exemple une onde générée en faisant 5 émettre des signaux acoustiques aléatoires par les différents transducteurs T1-Tn, - soit une onde ultrasonore de compression focalisée en un ou plusieurs points du milieu 2. Pour observer la propagation de l'onde de cisaille10 ment dans le milieu 2, on procède en plusieurs étapes successives: (a) une étape d'excitation au cours de laquelle le micro-ordinateur 4 fait générer une onde élastique de cisaillement dans le milieu viscoélastique 2, en faisant 15 émettre au moins une onde ultrasonore focalisée dans le milieu viscoélastique par la sonde 6, (b) une étape d'observation au cours de laquelle on observe la propagation de l'onde de cisaillement simultanément en une multitude de points du champ d'observation dans 20 le milieu viscoélastique 2, cette étape d'observation comprenant les sous- étapes suivantes: (bl) le micro-ordinateur 4 fait émettre dans le milieu viscoélastique, par la sonde 6 une succession de tirs d'ondes ultrasonores de compression non focalisées à 25 une cadence d'au moins 500 tirs par seconde (la focalisation et la chronologie de l'onde ultrasonore focalisée émise à l'étape (a), ainsi que la chronologie desdites ondes ultrasonores non focalisées, sont adaptées pour qu'au moins certaines desdites ondes ultrasonores non focalisées 30 parviennent dans le champ d'observation lors de la propagation de l'onde de cisaillement dans ce champ d'observation, pour au moins certaines émissions d'onde ultrasonore non focalisée), (b2) le micro-ordinateur 4 fait détecter par la sonde 6 et enregistrer en temps réel des signaux acoustiques reçus du milieu viscoélastique 2, comprenant les échos générés par les ondes ultrasonores de compression non focalisées en interagissant avec les particules réfléchissantes 5 5 du milieu viscoélastique, ces échos correspondant (directement ou indirectement) à des images successives de déplacement du milieu viscoélastique, (c) et au moins une étape de traitement au cours de laquelle: (cl) le micro- ordinateur 4 traite les signaux acoustiques successifs reçus du milieu viscoélastique 2 au cours de la sous-étape (b2) pour déterminer des images de propagation successives, (c2) et le micro-ordinateur 4 détermine au moins 15 un paramètre de mouvement du milieu viscoélastique 2 en différents points du champ d'observation. L'onde ultrasonore focalisée émise au cours de l'étape d'excitation (a) peut être une onde monochromatique de fréquence f comprise entre 0,5 et 15 MHz, par exemple 20 environ 4 MHz, émise pendant une durée k/f en secondes, o k est un entier compris entre 50 et 5000 (par exemple de l'ordre de 500) et f est exprimé en Hz. Eventuellement, une telle onde peut être émise pendant une succession de périodes d'émission séparées par des périodes de repos, les pé25 riodes d'émission se succédant à une cadence comprise entre et 1000 émissions par seconde. En variante, l'onde ultrasonore focalisée émise au cours de l'étape d'excitation (a) est une combinaison linéaire (notamment une somme) de deux signaux monochromati30 ques de fréquences respectives fl et f2 telles que 20 Hz Ifl - f21 Par ailleurs, l'onde ultrasonore focalisée émise au cours de l'étape d'excitation (a) peut éventuellement être focalisée, simultanément ou non, en plusieurs points, de façon que l'onde de cisaillement générée présente la forme souhaitée (par exemple, on peut ainsi générer une onde de cisaillement plane, ou au contraire une onde de cisaille5 ment focalisée) et illumine les zones souhaitées du milieu 2. Au cours de l'étape (bl), qui peut durer par exemple moins d'une seconde, on peut émettre les ondes ultrasonores de compression non focalisées, à une cadence comprise 10 entre 500 et 10 000 tirs par seconde et de préférence comprise entre 1000 et 5 000 tirs par seconde (cette cadence est limitée par le temps d'aller-retour de l'onde de compression dans le milieu 2, donc par l'épaisseur du milieu 2 dans la direction X: il faut en effet que tous les échos 15 générés par l'onde de compression aient été reçus par la sonde 6 avant d'envoyer une nouvelle onde de compression). Chaque onde ultrasonore de compression non focalisée se propage dans le milieu 2 avec une vitesse de propagation beaucoup plus élevée que les ondes de cisaillement 20 (par exemple de l'ordre de 1500 m/s dans le corps humain), et interagit avec les particules réfléchissantes 5, ce qui génère des échos ou autres perturbations analogues du signal, connus en soi sous le nom de "bruits de speckle" dans le domaine de l'échographie. Ces "bruits de speckle" sont captés par les transducteurs Tl,...,Tn au cours de la sous-étape (b2), après chaque tir d'onde ultrasonore de compression non focalisée. Le signal sij(t) ainsi capté par chaque transducteur Ti après le tir n j est tout d'abord échantillonné à haute 30 fréquence (par exemple de 30 à 100 MHz) et numérisé en temps réel (par exemple sur 12 bits) par un échantillonneur appartenant à la baie 7 et relié à ce transducteur, respectivement El, E2,... En. Le signal sij(t) ainsi échantillonné et numérisé est ensuite mémorisé, également en temps réel, dans une mémoire Mi appartenant à la baie 7 et propre au transducteur Ti. Chaque mémoire Mi présente par exemple une capacité 5 de l'ordre de 128 Mo, et contient l'ensemble des signaux sij(t) reçus successivement pour les tirs j = 1 à p. En temps différé, après la mémorisation de tous les signaux sij (t) correspondant à une même propagation d'onde de cisaillement, l'unité centrale CPU fait retraiter ces 10 signaux par un circuit sommateur S appartenant à la baie 7 (ou bien elle effectue elle-même ce traitement, ou encore ledit traitement peut être effectué dans le microordinateur 4), par un processus classique de formation de voies correspondant à la sous-étape (cl). On génère ainsi des signaux Sj (x, y) qui correspondent chacun à l'image du champ d'observation après le tir n j (dans le cas o l'onde ultrasonore non focalisée est une onde plane). Par exemple, on peut déterminer un signal Sj (t) 20 par la formule suivante: n Sj(t) = a, (x, y).s/[t(x, y) + d, (x, y) / V] i=1 o: - sij est le signal brut perçu par le transducteur n i après le tir d'onde ultrasonore de compression n j, - t(x,y) est le temps mis par l'onde ultrasonore de compression pour atteindre le point du champ d'observation de coordonnées (x,y), avec t = 0 au début du tir n j, - di(x,y) est la distance entre le point du champ 30 d'observation de coordonnées (x,y) et le transducteur n i, ou une approximation de cette distance, - V est la vitesse moyenne de propagation des on- des acoustiques ultrasonores de compression dans le milieu viscoélastique observé. - et ai(x,y) est un coefficient de pondération tenant compte de lois d'apodisation (en pratique, on pourra dans de nombreux cas considérer que ai(x,y)=l). La formule ci-dessus s'applique mutatis mutandis lorsque le champ d'observation est à 3 dimensions (réseau bidimensionnel de transducteurs), en remplaçant les coordonnées spatiales (x,y) par (x,y, z). Après l'étape éventuelle de formation de voies, l'unité centrale CPU mémorise dans une mémoire centrale M appartenant à la baie 7 les signaux d'images Sj(x,y) ou Sj(x) ou Sj(x,y,z), qui correspondent chacun au tir n j. Ces signaux peuvent également être mémorisés dans le micro15 ordinateur 4 lorsqu'il effectue lui-même le traitement d'image. Ces images sont ensuite traitées en temps différé à la sous-étape (c2), par corrélation et avantageusement par intercorrélation soit deux à deux, soit de préférence avec 20 une image de référence qui peut être: - soit une image de déplacement précédemment déterminée comme expliqué ci-dessus et utilisée comme image de référence pour les images de déplacement ultérieures (ou pour un nombre limité d'images de déplacement ultérieures, 25 par exemple 30 images de déplacement), - soit déterminée au cours d'une étape préliminaire d'observation initiale (aO), comme les images successives de déplacement susmentionnées, en faisant émettre une ou plusieurs ondes ultrasonores non focalisées par la sonde 30 6 avant l'étape d'excitation (a) qui génère l'onde de cisaillement (lorsque plusieurs ondes ultrasonores de compression non focalisées sont ainsi émises avant la phase d'excitation, on enregistre des échos générés par chaque onde ultrasonore de compression non focalisée en interagis- sant avec les particules réfléchissantes du milieu viscoélastique, ces échos correspondant à plusieurs images préliminaires successives du milieu viscoélastique, et on détermine ladite image initiale du milieu viscoélastique en com5 binant lesdites images préliminaires successives et notamment en moyennant les valeurs des pixels desdites images préliminaires). L'intercorrélation susmentionnée peut être réalisée par exemple dans un circuit électronique spécialisé DSP ap10 partenant à la baie 7, ou être programmée dans l'unité centrale CPU ou dans le microordinateur 4. Au cours de ce processus d'intercorrélation, on maximise une fonction d'intercorrélation afin de déterminer le déplacement subi 15 par chaque particule 5 donnant lieu à un écho ultrasonore. Des exemples de tels calculs d'intercorrélation sont donnés dans l'état de la technique, notamment par O'Donnell et al. ("Internal displacement and strain imaging using speckle tracking", IEEE transactions on ultrasonic, 20 ferroelectrics, and frequency control, vol. 41, n 3, mai 1994, p. 314-325) et par Ophir et al. ("Elastography: a quantitative method for imaging the elasticity of biological tissues", Ultrasonic imag. , vol. 13, p.111-134, 1991). On obtient ainsi un ensemble de vecteurs déplace25 ments u(?,t) engendrés par les ondes de cisaillement en chaque position r du milieu 2 sous l'effet de l'onde de cisaillement (ces vecteurs déplacements peuvent éventuellement être réduits à une seule composante dans l'exemple considéré ici). Cet ensemble de vecteurs déplacements est stocké dans la mémoire M ou dans le micro-ordinateur 4 et peut par exemple être visualisé, notamment au moyen de l'écran 4a du micro-ordinateur, sous la forme d'un film ralenti o la valeur des déplacements est illustrée par un paramètre opti- que tel que par un niveau de gris ou par un niveau chromatique. On visualise ainsi parfaitement les différences de propagation de l'onde de cisaillement entre les zones de 5 caractéristiques différentes du milieu 2, par exemple les tissus sains et les tissus cancéreux dans le cas d'une application médicale. Ce film de propagation de l'onde de cisaillement est en outre superposable avec une image échographique 10 classique, qui peut être générée également par le dispositif 1 décrit ci-dessus. Par ailleurs, il est également possible de calculer non pas les déplacements de chaque point du milieu observé 2, mais les déformations du milieu 2, c'est à dire des vec15 teurs dont les composantes sont les dérivées des composantes des vecteurs déplacements respectivement parrapport aux variables d'espace (coordonnées selon X, Y dans l'exemple considéré). Ces vecteurs de déformation sont utilisables comme les vecteurs déplacements pour visualiser clai20 rement la propagation de l'onde de cisaillement sous la forme d'un film, et présentent en outre l'avantage de s'affranchir des déplacements de la sonde 6 par rapport au milieu observé 2. A partir des champs de déplacements ou de déforma25 tions, le microordinateur 4 peut avantageusement procéder ensuite à une étape de cartographie (d) au cours de laquelle, à partir de l'évolution du paramètre de mouvement (déplacement ou déformation) au cours du temps dans le champ d'observation X, Y (ou X, Y, Z dans le cas d'un ré30 seau bidimensionnel de transducteurs), on calcule au moins un paramètre de propagation de l'onde de cisaillement, soit en certains points du champ d'observation choisis par l'utilisateur à partir du micro-ordinateur 4, soit dans tout le champ d'observation. Le paramètre de propagation de l'onde de cisaillement qui est calculé au cours de l'étape de cartographie est choisi par exemple parmi: la vitesse Cs des ondes de cisaillement, le module de cisaillement pi ou le module 5 d'Young E=3pt, l'atténuation a des ondes de cisaillement, l'élasticité de cisaillement i1, la viscosité de cisaillement t2, ou le temps de relaxation mécanique ts des tissus. Par exemple, on peut calculer en différents points du champ d'observation: - la valeur de la célérité Cs de l'onde de cisaillement, qui donne accès à la dureté des tissus, - la valeur du temps de relaxation mécanique -s des tissus, caractéristique de la viscosité locale du milieu. Pour cela, on utilise l'équation de propagation (1) suivante, à laquelle obéissent les déplacements u(,t) engendrés par les ondes de cisaillement en chaque position r du milieu: a 2(?, t) 2 o ôt 2=> cs2 (1+v-rs).V2u(rt) (1) o p est la densité des tissus, TS est le temps de relaxation mécanique des tissus et cS est la célérité de l'onde de cisaillement, directement reliée au module d'Young E des tissus par la relation: es (2) Dans le domaine de Fourier, l'équation d'onde (1) ci-dessus peut s'écrire: o2pU(F, o) = cs2(1 +jor).AU(,) (3) o U(QV)est la transformé de Fourier du champ de déplacement u(?,t) mesuré en chaque point et AU(?,w)est la transformée de Fourier du laplacien spatial de ce champ u(?,t) Etant donné que wTS " 1, on a donc accès à une expression simplifiée: c 2 U(,2o) S A U(?, co) S AU(r,@(4) s = 1 tan(( U(, Co)) AU(,) (5) o U(x) est la phase de la variable complexe x. Les fonctions U(o.)et AU(7,o) étant connues en chaque point de l'image échographique, il est donc possible de mesurer en tout point de l'espace le module d'Young et le 15 temps de relaxation mécanique des tissus, en établissant ainsi une cartographie de ces deux paramètres. De plus, les équations (4) et (5) étant vraies à chaque fréquence, le calcul de cS et TS peut être avantageusement moyenné sur l'ensemble de la bande des fréquences 20 portées par l'onde de cisaillement, améliorant ainsi fortement la qualité de la cartographie réalisée. A cet effet, on peut utiliser les formules suivantes: 2 1 12 U(F,w) d) Co -0 OO AU(?,W) (6) T o = |-ltan(Â( U(F,_))) dco I-@oSo AU(r,) (7) o 0 et X 1 sont les fréquences minimales et maximales que 30 porte l'onde de cisaillement Le mode de calcul serait le même en utilisant non plus les déplacements, mais les déformations du milieu observé 2. Par ailleurs, on peut avantageusement établir suc5 cessivement plusieurs cartographies des paramètres de propagation souhaités, par exemple cS et TS, en générant successivement des ondes de cisaillement différentes, obtenues par exemple en émettant des ondes ultrasonores de compression focalisées successivement en plusieurs points ou ayant 10 des formes d'ondes différentes. On peut ensuite combiner les différentes cartographies obtenues, par exemple en les moyennant, de façon à obtenir une cartographie synthétique plus riche et plus précise. REVENDICATIONS 1. Procédé d'imagerie utilisant des ondes de cisaillement pour observer un milieu viscoélastique diffusant 5 (2) qui contient des particules (5) réfléchissant les ondes ultrasonores de compression, procédé comprenant: (a) une étape d'excitation au cours de laquelle on génère une onde élastique de cisaillement dans le milieu viscoélastique (2), (b) une étape d'observation au cours de laquelle on observe la propagation de l'onde de cisaillement simultanément en une multitude de points d'un champ d'observation dans le milieu viscoélastique (2), cette étape d'observation comprenant les sous-étapes suivantes: (bl) faire émettre dans le milieu viscoélastique (2), par un réseau de transducteurs (6) commandés indépendamment les uns des autres, une succession de tirs d'ondes ultrasonores de compression non focalisées à une cadence d'au moins 500 tirs par seconde, (b2) faire détecter et enregistrer en temps réel des signaux acoustiques reçus du milieu viscoélastique (2), comprenant les échos générés par les ondes ultrasonores de compression non focalisées en interagissant avec les particules réfléchissantes (5) dudit milieu viscoélastique, (c) et au moins une étape de traitement au cours de laquelle: (cl) on traite les signaux acoustiques successifs reçus du milieu viscoélastique (2) au cours de la sous-étape (b2) pour déterminer des images de propagation 30 successives de l'onde de cisaillement, (c2) et on détermine au moins un paramètre de mouvement du milieu viscoélastique (2) en différents points du champ d'observation, caractérisé en ce qu'au cours de l'étape d'excitation (a), on fait générer l'onde élastique de cisaillement en faisant émettre au moins une onde ultrasonore focalisée dans le milieu viscoélastique (2) par ledit réseau de transducteurs (6), la focalisation et la chronologie de ladite onde ul5 trasonore focalisée, ainsi que la chronologie desdites ondes ultrasonores non focalisées, étant adaptées pour qu'au moins certaines desdites ondes ultrasonores non focalisées parviennent dans le champ d'observation lors de la propagation de l'onde de cisaillement dans ce champ d'observation, 10 pour au moins certaines émissions d'onde ultrasonore non focalisée. 2. Procédé selon la revendication 1, dans lequel, au cours de la sous-étape (b2), pour déterminer ledit paramètre de mouvement, on compare plusieurs images successives 15 de propagation avec une même image de référence du milieu viscoélastique (2), cette image de référence était déterminée en tirant au moins une onde ultrasonore de compression non focalisée dans ledit milieu viscoélastique puis en détectant et en enregistrant des échos générés par ladite 20 onde ultrasonore de compression non focalisée lorsqu'elle interagit avec les particules réfléchissantes (5) du milieu viscoélastique. 3. Procédé selon la revendication 2, dans lequel l'étape (a) est précédée d'une étape d'observation initiale 25 (aO) au cours de laquelle on tire au moins une onde ultrasonore de compression non focalisée puis on détecte et on enregistre des échos générés par ladite onde ultrasonore de compression non focalisée en interagissant avec les particules réfléchissantes (5) du milieu viscoélastique, ces 30 échos correspondant à une image initiale du milieu viscoélastique, et au cours de la sous-étape (b2) ladite image initiale constitue ladite image de référence pour traiter au moins certaines des images successives de déplacement. 4. Procédé selon la revendication 3, dans lequel, au cours de l'étape d'observation initiale (aO), on tire successivement plusieurs ondes ultrasonores de compression non focalisées puis on détecte et on enregistre des échos générés par chaque onde ultrasonore de compression non fo5 calisée en interagissant avec les particules réfléchissantes (5) du milieu viscoélastique, ces échos correspondant à plusieurs images successives du milieu viscoélastique, et on détermine ladite image initiale du milieu viscoélastique en combinant lesdites images successives. 5. Procédé selon l'une quelconque des revendications précédentes, dans lequel ledit paramètre de mouvement est un déplacement du milieu viscoélastique (2). 6. Procédé selon l'une quelconque des revendications précédentes, dans lequel l'onde ultrasonore focalisée 15 émise au cours de l'étape d'excitation (a) présente une fréquence f comprise entre 0,5 et 15 MHz et est émise pendant une durée k/f en secondes, o k est un entier compris entre 50 et 5000 et f est exprimé en Hz. 7. Procédé selon l'une quelconque des revendica20 tions 1 à 5, dans lequel l'onde ultrasonore focalisée émise au cours de l'étape d'excitation (a) présente une fréquence comprise entre 0,5 et 15 MHz et est émise pendant une succession de périodes d'émission séparées par des périodes de repos, les périodes d'émission se succédant à une cadence 25 comprise entre 10 et 1000 émissions par seconde. 8. Procédé selon l'une quelconque des revendications 1 à 5, dans lequel l'onde ultrasonore focalisée émise au cours de l'étape d'excitation (a) est une combinaison linéaire (notamment une somme) de deux signaux monochroma30 tiques de fréquences respectives fl et f2 telles que 20 Hz 9. Procédé selon l'une quelconque des revendications précédentes, dans lequel l'onde ultrasonore focalisée émise au cours de l'étape d'excitation (a) est focalisée simultanément en plusieurs points. 10. Procédé selon l'une quelconque des revendications précédentes, dans lequel l'étape (c) de traitement d'images est suivie par une étape (d) de cartographie au 5 cours de laquelle, à partir d'une évolution du paramètre de mouvement au cours du temps, on calcule au moins un paramètre de propagation de l'onde de cisaillement en au moins certains points du champ d'observation pour déterminer ainsi une cartographie dudit paramètre de propagation dans le 10 champ d'observation. 11. Procédé selon l'une quelconque des revendications précédentes, dans lequel le paramètre de propagation de l'onde de cisaillement qui est calculé au cours de l'étape de cartographie (d) est choisi parmi la vitesse des 15 ondes de cisaillement, le module de cisaillement, le module d'Young, l'atténuation des ondes de cisaillement, l'élasticité de cisaillement, la viscosité de cisaillement et le temps de relaxation mécanique. 12. Procédé selon la revendication 11, dans lequel 20 on répète successivement les étapes (a) à (d) en émettant des ondes ultrasonores focalisées différentes au cours des étapes d'excitation (a) successives, puis on combine les différentes cartographies obtenues au cours des étapes de cartographie (d) successives, pour calculer une cartogra25 phie synthétique du champ d'observation. 13. Dispositif d'imagerie pour la mise en oeuvre d'un procédé selon l'une quelconque des revendications précédentes, utilisant des ondes de cisaillement pour observer un milieu viscoélastique diffusant (2) qui contient des 30 particules (5) réfléchissant les ondes ultrasonores de compression, ce dispositif comprenant un réseau de transducteurs (6) commandés indépendamment les uns des autres par au moins une unité centrale électronique (4, CPU) adaptée pour: - faire générer au moins une onde élastique de cisaillement dans le milieu viscoélastique (2), - observer la propagation de l'onde de cisaillement simultanément en une multitude de points d'un champ 5 d'observation dans le milieu viscoélastique (2), en faisant émettre dans le milieu viscoélastique, par ledit réseau de transducteurs (6), une succession de tirs d'ondes ultrasonores de compression non focalisées, à une cadence d'au moins 500 tirs par seconde, puis en faisant détecter en 10 temps réel par ledit réseau de transducteurs (6), et pour enregistrer en temps réel, des signaux acoustiques reçus du milieu viscoélastique (2), comprenant les échos générés par les ondes ultrasonores de compression non focalisées en interagissant avec les particules réfléchissantes (5) dudit 15 milieu viscoélastique, - et traiter les signaux acoustiques successifs reçus du milieu viscoélastique (2) pour déterminer des images de propagation successives de l'onde de cisaillement, puis déterminer au moins un paramètre de mouvement du mi20 lieu viscoélastique (2) en différents points du champ d'observation, caractérisé en ce que l'unité centrale électronique (4) est adaptée pour faire générer l'onde élastique de cisaillement en faisant émettre au moins une onde ultrasonore focalisée 25 dans le milieu viscoélastique par ledit réseau de transducteurs (6), la focalisation et la chronologie de ladite onde ultrasonore focalisée, ainsi que la chronologie desdites ondes ultrasonores non focalisées, étant adaptées pour que lesdites ondes ultrasonores non focalisées parviennent dans 30 le champ d'observation lors de la propagation de l'onde de cisaillement dans ce champ d'observation.