Brevet FR3095941A1 — publié le 20201120

Description Titre de l'invention : Procédé pour déterminer par ultrasons une image corrigée d'un milieu, et dispositif pour mettre en oeuvre ce procédé Domaine technique 100011 La présente divulgation est relative aux procédés et dispositifs d'imagerie par ondes ultrasonores pour fournir des images d'un milieu viscoélastique. Notamment, le procédé a pour but de corriger les images produites en prenant en compte une variation de pression exercée par une sonde ultrasonore d'imagerie sur une surface externe du milieu. Technique antérieure [0002] Plus particulièrement, la présente divulgation concerne un procédé d'imagerie utilisant des ondes ultrasonores pour observer un milieu viscoélastique diffusant qui contient des particules réfléchissantes desdites ondes ultrasonores de compression. [0003] Les images corrigées par le procédé peuvent être de divers type. Par exemple, les images peuvent être de type déformation (« strain » en langue anglaise) ou de type élasticité (« elasticity » en langue anglaise). [0004] Le document US 2005/252295 propose une technique d'imagerie d'élasticité pour laquelle on génère une onde de cisaillement élastique dans le milieu_ et on observe la propagation de cette onde de cisaillement. élastique en enregistrant une pluralité d'images du milieu On détermine alors des images d'au moins un paramètre de mouvement du milieu (déplacement, déformation partielle), et on en déduit des images d'un paramètre de propagation de l'onde de cisaillement, tel que la vitesse de l'onde de cisaillement, le module de cisaillement., le module d'Young, l'élasticité de cisaillement, etc... [0005] Cette dernière technique est très utilisée. Cependant, il a été constaté que les images du paramètre de mouvement ou les images du paramètre de propagation de l'onde de cisaillement sont sensibles à la pression exercée par la sonde ultrasonore d'imagerie sur la surface du milieu. Or, cette pression de la sonde est exercée par le praticien et peut varier lors d'un examen ou d'un examen à un autre. En effet, le praticien fait varier cette pression externe pour obtenir l'observation désirée sur l'image considérée. [0006] En outre, ces images comprennent des déformations liées aux variations de la pression externe exercée, ce qui génère des imprécisions dans leur représentation, dans leur reproductibilité et dans leur détermination. Exposé de l'invention [0007] La présente invention a pour but de perfectionner les procédés d'imagerie de ce type, 2 notamment pour en améliorer les images. [0008] A cet effet, la présente divulgation propose un procédé pour déterminer par ultrasons une image corrigée d'un milieu caractérisée en ce que : [0009] - la détermination d'une succession temporelle d'images du milieu (Ik), k étant un indice d'image compris entre 0 et N, et la première image de ladite succession étant supposée sans déformation, [0010] - la détermination du mouvement pour chaque image du milieu d'indice k entre 1 et N, le mouvement étant entre une image du milieu d'indice k (Ik) et l'image du milieu précédente d'indice k-1 (Ik_!), [0011] - la détermination d'une succession temporelle d'images de déformation partielle (lAFk) à partir de l'image du milieu d'indice k (Ik) et l'image du milieu précédente d'indice k-1 (Ik_!), [0012] - la détermination d'images de déformation partielle corrigée (IAF.L") par com- pensation desdites images de déformation partielle (IAFk) à partir des mouvements précédents, c'est-à-dire les mouvements associés aux images du milieu d'indices compris entre 1 et k inclus, et [0013] - la détermination d'une image de déformation corrigée (IFL*) par une somme des images de déformation partielle corrigée (MEC) des indices entre 1 et k inclus. [0014] Grâce à ces dispositions, on assure que l'image produite est corrigée. [0015] Dans divers modes de réalisation du procédé selon la présente divulgation, on peut éventuellement avoir recours en outre à l'une et/ou à l'autre des dispositions suivantes. [0016] Selon un aspect, la détermination du mouvement à l'indice k est effectuée par : [0017] - détermination d'un champ de déplacement entre une image du milieu et l'image du milieu précédente, [0018] - détermination d'une transformation géométrique d'image à partir dudit champ de déplacement, ladite transformation géométrique représentant le champ de déplacement avec un nombre de paramètres inférieur à dix. [0019] Selon un aspect, le champ de déplacement est calculé par intercorrélation entre une image du milieu et l'image du milieu précédente, ou par un algorithme de suivi de sous-images entre une image du milieu et l'image du milieu précédente. [0020] Selon un aspect, l'algorithme de suivi est un algorithme de Lucas-Kanade. [0021] Selon un aspect, la transformation géométrique comprend au moins une translation, ou une translation et une homothétie, ou une translation, une homothétie et une rotation. [0022] Selon un aspect, la transformation géométrique Tk comprend une translation, une ho- mothétie et une rotation, et dans lequel la transformation géométrique est mise sous une forme matricielle telle que : [0023] 3 1lx. case 11 -z si ne Tx Hx. sine Hz. case Tz 0 0 1 [0024] avec les paramètres de la transformation géométrique Tk dans le plan d'image suivants : [0025] Tx, Tz coefficients de translation, [0026] Hx, Hz coefficient d'homothétie, et [0027] 0 angle de rotation d'axe perpendiculaire au plan d'image [0028] Selon un aspect, des paramètres de la transformation géométrique sont obtenus par des valeurs médianes d'une population de paramètres calculés à partir de groupes de points du champ de déplacement. [0029] Selon un aspect : [0030] - un groupe comprend entre trois et dix points du champ de déplacement, et [0031] - la population est supérieure à cent groupes. [0032] Selon un aspect, le procédé comprend en outre : [0033] (el) la détermination d'une image (l'élasticité du milieu. et [0034] (e2) la détermination d'une image d'élasticité corrigée par compensation de ladite image d'élasticité à partir des mouvements précédents, c'est-à-dire les mouvements associés aux images du milieu d'indices compris entre 0 et k inclus. [0035] Selon un aspect, l'image d'élasticité produite à l'étape (el) est générée par les sous- étapes suivantes : [0036] (e1.1) une étape d'excitation au cours de laquelle on génère une onde de cisaillement. dans le milieu, en faisant émettre au moins une onde ultrasonore focalisée, [0037] (e1.2) une étape d'observation au cours de laquelle on observe la propagation de l'onde de cisaillement en déterminant une succession temporelle (l'images intermédiaires du milieu, j étant un indice d'image intermédiaire compris entre 0 et M, [0038] (e1.3) une étape de traitement au cours de laquelle on détermine l'image d'élasticité à partir desdites images intermédiaires du milieu et d'un modèle de propagation d'onde de cisaillement. 100391 Selon un aspect, on effectue en outre : 100401 - la détermination d'une pluralité d'images de déformation corrigée, 100411 - la détermination d'une pluralité d'images d'élasticité corrigée, lesdites images de déformation corrigé et lesdites images d'élasticité étant temporellement entrelacées, et 100421 - la détermination d'une image d'un paramètre de non linéarité à partir de la pluralité d'images de déformation corrigée et de la pluralité d'images d'élasticité corrigée. 100431 Selon un aspect, la valeur de chaque pixel de l'image d'un paramètre de non linéarité est déterminée par régression linéaire en déformation de couples de valeurs, la première valeur de chaque couple correspondant à la valeur d'un même pixel d'une image de déformation corrigée de la pluralité et la deuxième valeur du couple correspondant à la valeur d'un même pixel d'une image d'élasticité de la pluralité, ladite image de déformation corrigée de la pluralité et ladite image d'élasticité corrigée de la pluralité étant temporellement successives ou proches, et [0044] dans lequel la régression linéaire en déformation est établie sur la hase de la relation entre la déformation et l'élasticité suivante : [0045] A ln (E) = ln (Bo) 4 . E E° [0046] où [0047] E est le module d'Young du pixel de l'image d'élasticité corrigée considérée, [0048] E1) est le module d'Young du pixel de la première image d'élasticité colligée, [0049] ln() est la fonction de logarithme népérien, [0050] e est la déformation du pixel de l'image de déformation corrigée considérée, et [0051] A est le paramètre de non linéarité déterminé dudit pixel par ladite régression linéaire en déformation. [0052] Selon un aspect, on effectue en outre : [0053] - la détermination d'une succession temporelle d'images de contrainte partielle corrigée par une multiplication de l'image de déformation partielle corrigée d'indice k et de l'image d'élasticité corrigée d'indice k, [0054] - la détermination d'une image de contrainte colligée par une somme des images de déformation partielle corrigée des indices entre 1 et k inclus. [0055] Selon un aspect, on effectue en outre : [0056] - la détermination d'une pluralité d'images de contrainte corrigée, [0057] - la détermination d'une pluralité d'images d'élasticité colligée, lesdites images de déformation corrigé et lesdites images d'élasticité étant temporellement entrelacées, et [0058] - la détermination d'une image d'un paramètre de non linéarité à partir de la pluralité d'images de contrainte corrigée et de la pluralité d'images d'élasticité corrigée. [0059] Selon un aspect, la valeur de chaque pixel de l'image d'un paramètre de non linéarité est déterminée par régression linéaire en contrainte de couples de valeurs, la première valeur de chaque couple correspondant à la valeur d'un même pixel d'une image de contrainte corrigée de la pluralité et la deuxième valeur du couple correspondant à la valeur d'un même pixel d'une image d'élasticité de la pluralité, ladite image de contrainte corrigée de la pluralité et ladite image d'élasticité corrigée de la pluralité étant temporellement successives ou proches, et [0060] dans lequel la régression linéaire en contrainte est établie sur la base de la relation entre la contrainte et l'élasticité suivante : [0061] A In (E) = ln (E0) 4 . £ 3. Et) [0062] où [0063] E est le module d'Young du pixel dc l'image d'élasticité corrigée considérée, [0064] El) est le module d'Young du pixel de la première image d'élasticité corrigée, [0065] ln() est la fonction de logarithme népérien, [0066] s est la déformation du pixel dc l'image dc déformation corrigée considérée, et [0067] A est le paramètre de non linéarité déterminé dudit pixel par ladite régression. [0068] La présente divulgation se rapporte également à un dispositif d'imagerie comprenant une sonde ultrasonore et un micro-ordinateur adaptés pour mettre en oeuvre le procédé pour déterminer par ultrasons une image corrigée d'un milieu tel que cité ci-dessus. Brève description des dessins 100691 D'autres caractéristiques et avantages de la présente divulgation apparaîtront au cours de la description suivante de divers modes de réalisation, donnés à titre d'exemple non limitatif, en regard des dessins joints. 100701 Sur les dessins : 100711 - la figure 1 est une vue schématique d'un dispositif d'imagerie par ultrasons selon une forme de réalisation dc la divulgation ; [0072] - les figures 2A et 2B montrent de manière schématique le mouvement d'une image du milieu entre une pression externe nulle et une pression externe qui déforme la surface externe du milieu ; [0073] - la figure 2C montre de manière schématique un mouvement plus complexe que celui de la figure 2B. ce mouvement incluant une rotation ; [0074] - la figure 3 est un diagramme présentant une forme de réalisation du procédé de calcul d'une image de déformation corrigée ; [0075] - la figure 4 est un diagramme présentant un procédé pour déterminer une image d'élasticité corrigée ; [0076] - la figure 5 est un diagramme présentant une première forme de réalisation d'un procédé pour déterminer une image d'un paramètre de non linéarité corrigé ; [0077] - la figure 6 est un graphique présentant un exemple de mise en oeuvre du procédé de la figure 5 ; [0078] - la figure 7A est un diagramme présentant un procédé pour déterminer une image de contrainte corrigée ; [0079] - la figure 7B est un diagramme présentant une deuxième forme de réalisation d'un procédé pour déterminer une image d'un paramètre de non linéarité corrigé ; et 6 100801 - la figure 8 est un graphique présentant un exemple de mise en oeuvre du procédé de la figure 7B. [0081] Sur les différentes figures, fournies à titre d'illustration, les mêmes références nu- mériques désignent des éléments identiques ou similaires. Description des modes de réalisation Dispositif d'imagerie 1 [0082] Le dispositif d'imagerie 1 représenté sur la figure 1 est destiné à fournir des images d'un milieu viscoélastique 2 qui est diffusant vis à vis des ondes ultrasonores de compression, et qui peut être par exemple un corps vivant, par exemple une partie du corps d'un patient (sein, foie, abdomen, ...), dans le cas des applications médicales. Ce dispositif d'imagerie 1 est également apte à étudier la propagation des ondes élastiques de cisaillement pour fournir des images d'élasticité du milieu 2. [0083] Les images du milieu sont produites par exemple au moyen d'un micro-ordinateur 4 (comprenant au moins une interface d'entrée 4b telle qu'un clavier ou autre, et une interface de sortie 4a 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 2, lesquelles particules sont réfléchissantes pour les ondes ultrasonores de compression. Les particules 5 peuvent être constituées par toute hétérogénéité du milieu 2, 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 sous le terme "speckle"). 100841 Pour observer le milieu 2 et pour générer des images du milieu, on utilise une sonde ultrasonore 6 disposée contre la surface extérieure 3 du milieu observé 2. Cette sonde envoie, selon un axe Z, des impulsions d'ondes ultrasonores 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. 100851 La sonde ultrasonore 6 est constituée par un réseau de n transducteurs ultrasonores Tl, T2, ..., Ti, ..., Tn, n étant un nombre entier supérieur ou au moins égal à 1. 100861 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 X perpendiculaire à l'axe Z. La sonde en question peut être également un réseau bidimensionnel (plan ou non) de transducteurs. 100871 Les transducteurs Tl, T2, Tn peuvent être commandés indépendamment les uns des autres par le micro-ordinateur 4, éventuellement par l'intermédiaire d'une unité centrale CPU qui est contenue par exemple dans une baie électronique 7 reliée par un 7 câble souple à la sonde 6. Les transducteurs Tl-Tn peuvent ainsi émettre sélectivement : [0088] - 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, Z) 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 émettre des signaux acoustiques aléatoires par les différents transducteurs Tl -Tn, [0089] - soit une onde ultrasonore de compression focalisée en un ou plusieurs points du milieu 2. [0090] Le document US 2009/234230 propose une technique d'imagerie synthétique utilisant plusieurs ondes de compression non focalisées, par exemple des ondes de type onde planes d'angles différents, et qui combine les ondes en retour de ces ondes planes pour obtenir de manière très rapide une image du milieu de qualité améliorée. Détermination d'images du milieu Ik [0091] Pour générer une image du milieu (1) le dispositif d'imagerie effectue par exemple les étapes suivantes : [0092] (i 1) le micro-ordinateur 4 fait émettre dans le milieu viscoélastique, par la sonde 6 au moins un tir d'onde ultrasonore de compression non focalisée, [0093] (i2) 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 l'onde ultrasonore de compression non focalisée en interagissant avec les particules réfléchissantes 5 du milieu viscoélastique, et [0094] (i3) une étape de traitement au cours de laquelle le micro-ordinateur 4 traite les signaux acoustiques reçus du milieu viscoélastique 2 au cours de la sous-étape (i2) pour déterminer une ou plusieurs images du milieu (1). [0095] L'onde ultrasonore de compression non focalisée se propage dans le milieu 2 avec une vitesse de propagation très élevée, par exemple de l'ordre de 1500 mis 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 speckde" dans le domaine de l'échographie. [0096] Ces "bruits de speckle" sont captés par les transducteurs Tl, , Tn au cours de la sous-étape (i2), après le tir d'onde ultrasonore de compression non focalisée. Le signal si(t) ainsi capté par chaque transducteur Ti est tout d'abord échantillonné à haute fréquence (par exemple de 30 à 100 MHz) et numérisé en temps réel par un échantillonneur appartenant à la baie 7 et relié à ce tnuisducteur, respectivement El, E2, En. [0097] Le signal si(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. 8 [0098] Chaque mémoire Mi présente par exemple une capacité de l'ordre de 128 Mo, et contient l'ensemble des signaux si(t) reçus. [0099] Après la mémorisation de tous les signaux si(t), l'unité centrale CPU fait retraiter ces 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 micro-ordinateur 4), par un processus classique de formation de voies correspondant à la sous-étape (i3). [0100] Des signaux S(x, z) sont ainsi générés, ces signaux correspondant chacun à l'image du champ d'observation du milieu après le tir de l'onde ultrasonore non focalisée. [0101] Par exemple, on peut déterminer un signal 5(t) par la formule suivante : [0102] S(X, Z, t) z). si[t(x, di(r, z)/V1 É=1 [0103] dans laquelle : [0104] si est le signal brut perçu par le transducteur n° i après le tir d'onde ultrasonore de compression, [0105] t(x,z) est le temps mis par l'onde ultrasonore de compression pour atteindre le point du champ d'observation de coordonnées (x,z), avec t = 0 au début du tir, [0106] di(x,z) est la distance entre le point du champ d'observation de coordonnées (x,z) et le transducteur n° i, ou une approximation de cette distance, [0107] V est la vitesse moyenne de propagation des ondes acoustiques ultrasonores de com- pression dans le milieu viscoélastique observé, et [0108] ai(x,z) est un coefficient de pondération tenant compte de lois d'apodisation (en pratique, on pourra dans de nombreux cas considérer que cci(x,z) = 1). [0109] 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 du plan (x,z) par des coordonnées spatiales (x,y,z). [0110] Après l'étape optionnelle de formation de voies, l'unité centrale CPU mémorise dans une mémoire centrale M appartenant à la baie 7 les signaux d'images S(x,z) correspondent au dernier tir. Ces signaux peuvent également être mémorisés dans le micro-ordinateur 4 afin qu'il effectue lui-même le calcul de l'image du milieu (I). [0111] D'autres techniques de génération d'image du milieu (I) existent, comme les techniques d'imagerie synthétique. Toutes technique (l'imagerie permettant d'obtenir des images du milieu est utilisable. De préférence on utilisera une technique permettant d'obtenir des images à grande cadence. [0112] Le dispositif d'imagerie 1 et le procédé selon la présente divulgation effectuent une succession temporelle (l'images du milieu pour suivre la variation de la pression 9 externe P exercée sur la surface externe 3. Ainsi, on considérera qu'un nombre N+1 d'images du milieu sont prises à des instants temporels successifs, ces instants temporels n'étant pas obligatoirement séparés d'une période constante. Les images du milieu peuvent donc être identifiées par un indice d'image k compris entre 0 et N. [0113] Les images du milieu seront donc identifiées par la notation Ik. [0114] Cet indice k sera utilisé pour toute autre image ou grandeur déterminée à partir de l'image du milieu d'indice k et éventuellement d'images du milieu précédente. Cet indice pourra donc être utilisé pour repérer l'instant temporel associé à cette image du milieu d'indice k. [0115] En outre, la première image de cette succession d'images est notée par convention l'image du milieu d'indice 0. On supposera que pour cette première image, la pression externe P est faible ou nulle, et que la déformation induite de la surface externe 3 est faible ou nulle. Cette première image est par exemple représentée en figure 2A. [0116] Les figures 2A à 2C illustrent le problème résolu par la présente divulgation. [0117] En figure 2A est représentée une utilisation du dispositif d'imagerie 1, la sonde 6 exerçant une pression externe P faible ou nulle sur la surface externe 3. La surface externe 3 demeure sensiblement horizontale (dans la direction X). L'image du milieu Io comprend par exemple une inclusion 2i à une profondeur Z1 par rapport à la surface externe 3. [0118] En figure 2B est représentée une utilisation du dispositif d'imagerie 1 à un instant temporel ultérieur de celui de la figure 2A, et avec une pression externe P' supérieure à P, qui déforme la surface externe 3 dans la direction Z vers l'intérieur du milieu 2. Le dispositif d'imagerie 1 produit une image du milieu II de même taille que l'image Io produite en figure 2A, mais les couches supérieures du milieu au-dessus de l'inclusion 2i sont naturellement plus comprimées que les couches inférieures du milieu 2 en dessous de l'inclusion 2i, de sorte que l'inclusion 2i est dans cette image du milieu II de la figure 2B déplacée vers la sonde 6 en comparaison de la première image L de la figure 2A. Autrement dit, l'inclusion 2i est alors à une profondeur Z2 par rapport à la surface externe 3 inférieure à la profondeur Z1 de la première image Io. [0119] Par conséquent, pendant une variation de la pression externe, les éléments internes du milieu se déplacent dans l'image du milieu à cause de l'élasticité dudit milieu 2. Ainsi, il est nécessaire de connaitre ces mouvements généraux pour effectuer un suivi et une analyse précise de caractéristiques de ces éléments internes. [0120] L'exemple des figures 2A et 2B montre qu'une translation dans la direction Z pourrait donner une correction adéquate. [0121] En figure 2C est représenté une utilisation du dispositif d'imagerie 1 à un autre instant temporel, par exemple encore ultérieur aux précédents instants temporels représentés aux figures 2A et 2B, et avec une pression externe P différente, et surtout un basculement de la sonde 6 par rapport à la direction verticale Z (basculement exagéré sur cette figure 2C). La surface externe 3 est déformée en direction Z et elle est inclinée d'un angle O. Le dispositif d'imagerie 1 produit une image du milieu In de même taille, mais les couches supérieures du milieu sont aussi comprimées et l'inclusion 2i est déplacée latéralement dans l'image, c'est-à-dire déplacée en direction X par l'inclinaison de la sonde 6. [0122] Ainsi, les variations de pression et de position de la sonde 6 sur la surface externe 3 déplacent les éléments internes du milieu dans l'image du milieu et il apparait nécessaire de connaitre ces mouvements généraux (très lents par rapport aux ondes ultrasonores) pour en effectuer un suivi et/ou pour les compenser dans les images du milieu prises afin d'effectuer une analyse précise de caractéristiques des éléments internes. [0123] La sonde 6 peut éventuellement avoir un déplacement plus complexe qu'une simple translation en direction Z ou qu'une combinaison d'une translation et d'une rotation. Le mouvement de correction nécessaire sera alors plus complexe. Les détails de la détermination de ces mouvements seront explicités plus loin dans la présente divulgation. Détermination d'un champ de déplacement uk [0124] Les images du milieu ('k) peuvent être traitées par corrélation et avantageusement par intercorrélation soit deux à deux, c'est-à-dire entre une image du milieu d'indice k (h) et l'image du milieu du milieu d'indice k-1 (I ,). [0125] L'intercorrélation peut être réalisée par exemple dans un circuit électronique spécialisé DSP appartenant à la baie 7, ou être programmée dans l'unité centrale CPU ou dans le micro-ordinateur 4. [0126] Au cours de ce processus d'intercorrélation, on maximise une fonction d'intercorrélation afin de déterminer le déplacement subi par chaque particule 5 donnant lieu à un écho ultrasonore. [0127] Des exemples de tels calculs d'intercorrélation sont donnés dans l'état de la technique, notamment dans les documents suivants : [0128] - "Internai displacement and strain imaging using speckle tracking", O'Donnell et al, IEEE transactions on ultrasonic, ferroelectrics, and frequency control, vol. 41, n° 3, mai 1994, p. 314-325 ), [0129] - "Elastography: a quantitative methodkr imaging die elasticity of biological tissues", Ophir et al., Ultrasonic imaging., vol. 13, p.111-134, 1991), et [0130] - "Pyramidal Implementation of the Lucas Kanade Feature Tracker. Description of the algorithm.", J-Y Bouguet, Intel Corp. [0131] On obtient ainsi un champ de déplacement, c'est à dire un ensemble de vecteurs dé- placements ou u(x,z,t) en chaque position (x,z) du milieu 2, que l'on peut noter uk(x,z) en remplaçant la variable de l'instant temporel t par l'indice d'image k. Ces vecteurs 11 déplacements peuvent éventuellement être réduits à une seule composante ou à deux ou trois composantes. Dans l'exemple considéré, le champ de déplacement uk à l'indice d'image k est : [0132] Uxi u= [uzi [0133] Ce champ de déplacement (ensemble de vecteurs déplacements) est stocké dans la mémoire M ou dans le micro-ordinateur 4. [0134] On peut construire une image du champ de déplacement Iuk, et l'ensemble des images du champ de déplacement Iuk correspondant à la succession temporelle des images du milieu peuvent être visualisés, notamment au moyen de l'écran 4a du micro-ordinateur, par exemple, sous la forme d'un film ralenti où la valeur des déplacements est illustrée par un paramètre optique tel que par un niveau de pis ou par un niveau chromatique. [0135] Détermination d'images de déformation partielle I At k [0136] On peut alors calculer les déformations partielles At en chaque point du milieu 2, c'est à dire des vecteurs dont les composantes sont les dérivées des composantes des vecteurs déplacements respectivement par rapport aux variables d'espace (coordonnées selon X, Z dans l'exemple considéré), c'est-à-dire : [0137] du,- dx du, As = dz [0138] Eventuellement, selon un exemple de mise en oeuvre, seule la déformation/dé- formation partielle uniaxiale en direction Z est calculée. En effet, la compression liée à la pression externe de la sonde est principalement dans cette direction. De la même manière, dans un calcul suivant de contrainte, seul la contrainte/contrainte partielle uniaxiale dans la direction Z est éventuellement calculée. On simplifie ainsi les calculs. [0139] Comme pour les vecteurs déplacements, on peut construire une image de dé- formation partielle (1Ask) à partir de l'ensemble des déformations partielles calcifiées en tout point de coordonnées (x, z) du plan X-Z (plan d'image). [0140] Ces vecteurs de déformations partielles ou images de déformation partielles (lAsk) peuvent être visualisés, notamment au moyen de l'écran 4a du micro-ordinateur, sous la forme d'un film ralenti. [0141] Détermination du mouvement du milieu par la pression externe [0142] La déformation du milieu 2 liée aux variations de la pression externe P exercée sur la surface extérieure 3 est une déformation lente, élastique quasi-uniforme. Cette déformation se propage à l'intérieur du milieu viscoélastique 2, et provoque le dé- 12 placement des particules 5 et des éléments inclus dans le milieu 2. [0143] Ce déplacement est néfaste à l'analyse des caractéristiques des divers points de l'image du milieu, et notamment l'évaluation de la déformation en ces points et/ou de l'élasticité de ces points. [0144] La présente divulgation considère donc de corriger des images pour compenser ces mouvements, ce qui permet de conserver les particules et/ou éléments internes du milieu 2 à des emplacements quasi constants dans l'image, permettant ainsi d'en étudier les caractéristiques plus précisément. [0145] La déformation externe provoquée par la pression externe P peut être considérée comme provoquant un mouvement entre une image du milieu d'indice k (ID et une image du milieu précédente d'indice k-1 (IkA). [0146] Le mouvement peut alors être déterminé par le champ de déplacement (uk) entre une image du milieu d'indice k et une image du milieu précédente d'indice k-1. [0147] Notamment, ce mouvement peut être modélisé par une transformation géométrique T is d'indice k qui représente d'une manière générale et de manière simple l'ensemble des déplacements du champ de déplacement à l'étape d'indice d'image k, uk(x, z), c'est-à-dire avec un nombre de paramètres très réduit, et par exemple inférieur à dix. [0148] Alors, une compensation de tout type d'images à partir des mouvements précédents peut être effectué par application des transformations géométriques précédentes, c'est-à-dire par application des transformations géométriques Ti, l'indice i variant de 1 à k. [0149] Selon une première variante, la transformation géométrique Tk comprend une translation, comme cela est représenté entre la figure 2A et la figure 2B. [0150] Selon une seconde variante, la transformation géométrique T, comprend une translation et une homothétie. [0151] Selon une seconde variante, la transformation géométrique T, comprend une translation, une homothétie et une rotation. [0152] Ainsi, dans le plan X-Z, plan d'image, la transformation géométrique peut être mise sous une forme matricielle du type : [0153] { HX'. cose -Hz.sine Tx 7'x.= Hx. sine Hz.cose Tz 0 0 1 10154] avec les paramètres de la transformation géométrique (Tk) suivants : 10155] Tx, Tz coefficients de translation dans le plan d'image, 10156] Hx, Hz coefficient d'homothétie le plan d'image, et 10157] () angle de rotation d'axe perpendiculaire au plan d'image. 10158] Ainsi, si l'on prend trois points Pl, P2, et P3 de l'image de cordonnées respectives 13 [0159] (xl,z1), (x2,z2), ((3,z3), on a les relations suivantes : [0160] [xl x2 x31 A = 21 z2 z3 1 1 1 xl ± tixl x2 + cix2 x3 + dx3 F I B = zi+ dzl z2 + clz2 z3 + dz3 [0161] avec [0162] chi = u', dz1 = uz1 [0163] dx2 = u'2dz2 = 142 [0164] ch.3 = u'3dz3 = 143 [0165] pour reprendre les notation des composantes X et Z des vecteurs déplacements des points Pl, P2 et P3. [0166] En outre, on a une relation permettant de calculer la transformation géométrique Tk à partir des matrices A et B précédentes, c'est-à-dire une matrice A correspondant aux cordonnées des trois points dans l'image, et une matrice B correspondant aux coordonnées de ces mêmes trois points avec les déplacements de ces points (le mouvement). Ainsi, la matrice de la transformation géométrique est obtenue par le produit matriciel de la matrice B et der inverse de la matrice A : [0167] Tk = B. A-1 [0168] L'application de cotte relation permet de calculer la transformation géométrique Tk à partir de trois point Pl, P2, P3. [0169] Inversement, une fois la transformation géométrique Tk établie, cotte relation permet de connaitrc les coordonnées (x, z) de tout point P par la relation inverse : [0170] B = Tk.A [0171] Alors, une compensation de tout type d'images à partir des mouvements précédents peut être effectué par multiplication matricielle des transformations géométriques précédentes, c'est-à-dire par multiplication matricielle des transformations géométriques T l'indice i variant de 1 à k. 10172] Les relations précédentes établies avec 3 points peuvent être généralisées à un groupe de points de l'image, le groupe de points comprenant 3, 4, 5 ou 6 points de l'image. Le groupe de points comprend entre trois et dix points de l'image. 10173] En outre, dans une variante, on prendra avantageusement une population de trois points (de groupe de points) de l'image, ladite population comprenant un nombre Ng de groupes de points de grande valeur. Par exemple, la taille de cette population, le nombre Ng de groupes de points est supérieur à cent. 14 [0174] Alors, les paramètres de la transformation géométrique Tk sont obtenus par des valeurs médianes des paramètres calculés à partir de la population des groupes de points de l'image. [0175] Notamment, cette technique peut être appliquée à une image du champ de dé- placement uk pour en déduire le mouvement entre une image du milieu d'indice k et une image du milieu précédente d'indice k-1. [0176] Grâce à l'utilisation d'une population de groupes de points, on peut déterminer une transformation géométrique qui représente dc manière plus globale le mouvement entre l'image du milieu d'indice k et une image du milieu précédente d'indice k-1 [0177] Avantageusement, la sélection des points dans chaque groupe se fait de manière aléatoire dans l'image. Grâce à cette disposition, on peut déterminer une transformation géométrique qui représente dc manière plus globale et plus fiable le mouvement entre images. [0178] Détermination d'images de déformation partielle corrigée I Ae k [0179] Les images de déformation partielle IAEk sont alors corrigées pour ramener ces images dans une situation conforme à la première image du milieu L, supprimant ainsi les mouvements provoqués par la déformation de la pression externe P dans ces images de déformation partielles IAFk. [0180] Cette correction est effectuée en utilisant à chaque image de déformation partielle IAFI, d'indice d'image k, l'ensemble des mouvements global déterminés entre chaque image du milieu Plus précisément, il faut en appliquer les mouvements inverses pour revenir dans la situation de la première image L. [0181] Autrement dit, les images de déformation partielle corrigées 1_4E( contiennent des déformations partielles con-igées AEk" et elles sont calculées par compensation desdites images de déformation partielle lAsk à partir des mouvements associés aux images du milieu d'indices compris entre 1 et k inclus. [0182] Détermination d'une ima e de déformation corri ée I F k [0183] La déformation sk(x,z) est obtenue en cumulant les déformations partielles suc- cessives, c'est-à-dire : [0184] ek (X, Z) = Aei (X34 i=1 10185] Dans le cas présent, la déformation sk(x,z) (déformation corrigée) est obtenue en cumulant les déformations partielles successives corrigées AFk. : 10186] Ek (x, z) = z) i=i 101871 On peut construire une image de déformation corrigée (lAskt) à partir de l'ensemble des déformations corrigées sk(x,z) calculées en tout point de coordonnées (x, z) du plan X-Z (plan d'image). [0188] Ces images de déformation corrigées Ir.k* peuvent être visualisés sous la forme d'un film ralenti. Procédé pour déterminer une image corrigée [0189] En résumé, selon une première forme de réalisation du procédé représenté en figure 3 selon la présente divulgation, le procédé mis en oeuvre comprend les étapes suivantes : [0190] (dl) ) la détermination d'une succession temporelle d'images du milieu (10, k étant un indice d'image compris entre 0 et N, et la première image de ladite succession étant supposée sans déformation, [0191] (d2) la détermination du mouvement pour chaque image du milieu d'indice k entre 1 et N, le mouvement étant entre une image du milieu d'indice k (Ik) et l'image du milieu précédente d'indice k-1 (I"), [0192] (d3) la détermination d'une succession temporelle d'images de déformation partielle (lAsk) à partir de l'image du milieu d'indice k (Ik) et l'image du milieu précédente d'indice k-1 (I"), [0193] (d4) la détermination d'images de déformation partielle corrigées (1-Ask') par com- pensation desdites images de déformation partielle (lAsk) à partir des mouvements précédents, c'est-à-dire les mouvements associés aux images du milieu d'indices compris entre 1 et k inclus, [0194] (d5) la détermination d'une image de déformation corrigée (Irr,*) par une somme ou cumul des images de déformation partielle corrigées (IAEr*) des indices compris entre 1 et k inclus. [0195] Par somme ou cumul d'images, on comprend que les valeurs de pixels de même co- ordonnée spatiale desdites images sont sommées (i.e ajoutées les unes aux autres) pour former la valeur du pixel correspondant de l'image résultante. [0196] Grâce à ces dispositions, on obtient une image de déformation corrigée Ie* du milieu 2 qui n'est pas ou peu influencée par les déformation et mouvements induits par la pression externe P, ce qui permet d'obtenir des valeurs de déformation plus exacte et plus précises que dans l'ait antérieur. [0197] Détermination d' une ima e d'élasticité du milieu IE k [0198] Le document US 2005/252295 propose une technique d'imagerie d'élasticité du milieu viscoélastique 2 : Le dispositif d'imagerie 1 étudie la propagation des ondes élastiques de cisaillement dans ce milieu Les mouvements des ondes élastiques de cisaillement sont suivis par les moyens décrits précédemment, et notamment par le micro-ordinateur 4. 16 [0199] On procède en plusieurs étapes : [0200] (c1.1) 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 émettre au moins une onde ultrasonore focalisée dans le milieu viscoélastique par la sonde 6, [0201] (c1.2) 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 le milieu viscoélastique 2, cette étape comprenant des sous-étapes durant lesquelles on génère une pluralité d'images intermédiaires du milieu ni,k temporellement successives, j étant un indice d'image intermédiaire compris entre 0 et M inclus, M+1 étant le nombre d'images intermédiaires produite. [0202] Chaque image intermédiaire du milieu est par exemple générée par le procédé de dé- termination d'une image du milieu précédemment décrit, dans lequel on fait émettre au moins un tir d'onde ultrasonore de compression non focalisée par la sonde 6, et on détecte et enregistre par la sonde 6 les signaux acoustiques reçus, et on traite ces signaux acoustiques pour construire une image intermédiaire du milieu (Hoc). [0203] La focalisation et la chronologie de l'onde ultrasonore focalisée émise à l'étape (e1), ainsi que la chronologie des ondes ultrasonores non focalisées émises à l'étape (e2), sont adaptées pour que des 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, l'onde de cisaillement est visible dans les images intermédiaires du milieu Hi,k générées. [0204] On procède alors à : [0205] (e1.3) une étape de traitement des images intermédiaires du milieu pour déterminer une image d'élasticité du milieu IEL. [0206] Les images intermédiaires du milieu Hoc sont traitées pour calculer un champ de dé- placement ui,k, par exemple par corrélation ou intercorrélation comme précédemment. Par exemple par intercorrélation entre une image intermédiaire du milieu d'indice j (Hi, k) et l'image intermédiaire du milieu du milieu d'indice j-1 u). Les vecteurs de dé- placement déterminés sont utilisables pour visualiser la propagation de l'onde de cisaillement sous la forme d'un film. On calcule éventuellement en outre les déformations partielles As en chaque point du milieu 2 à partir du champ de déplacement. [0207] A partir des champs de déplacements ou de déformations, le micro-ordinateur 4 peut avantageusement procéder ensuite à une étape de cartographie au cours de laquelle, à partir de l'évolution du paramètre de mouvement (déplacement ou déformation partielle) au cours du temps dans le champ d'observation, on calcule au moins un paramètre de propagation de l'onde de cisaillement, soit en certains points du champ d'observation choisis, soit dans tout le champ d'observation. 17 [0208] On peut alors construire une image d'élasticité TEk, correspondant à l'ensemble des paramètres des de propagation de l'onde de cisaillement aux divers points du champ d'observation. [0209] Cette image d'élasticité IEk du milieu peut être visualisée, notamment au moyen de l'écran 4a du micro-ordinateur où la valeur du paramètre de propagation est illustrée par un paramètre optique tel que par un niveau de gris ou par un niveau chromatique. [0210] 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 : [0211] - la vitesse Cs des ondes de cisaillement, ou [0212] - le module de cisaillement Ft, ou [0213] - le module d'Young E=3p, ou [0214] - l'atténuation a des ondes de cisaillement, ou [0215] - l'élasticité du cisaillement pl, ou [0216] - la viscosité de cisaillement Ft2, ou [0217] - le temps de relaxation mécanique TS des tissus du milieu. [0218] Par exemple, on peut calculer en différents points du champ d'observation : [0219] - la valeur de la célérité Cs de l'onde de cisaillement, qui donne accès à la dureté des tissus, [0220] - la valeur du temps de relaxation mécanique -rs des tissus, caractéristiques de la viscosité locale du milieu. [0221] Pour cela, on utilise un modèle de propagation de l'onde de cisaillement, par exemple représenté par l'équation de propagation suivante, à laquelle obéissent les déplacements u engendrés par les ondes de cisaillement en chaque position r du milieu : [0222] 5217(2F,t) ô - p 5t =c3 2 (1 + rs 71).v (F,t) 102231 Où 102241 p est la densité des tissus, 102251 TS est le temps de relaxation mécanique des tissus, et 102261 cS est la célérité de l'onde de cisaillement, directement reliée au module d'Young E des tissus par la relation : 102271 Cs = [0228] La résolution de cette équation de propagation avec l'ensemble des déplacements u permet d'obtenir les paramètres de propagation (cS, tS) cités ci-dessus. [0229] Des variantes de calcul du/des paramètre(s) de propagation sont possibles. Notamment, on peut utiliser l'équation d'onde dans le domaine de Fourier, en 18 moyennant par exemple les valeurs sur une bande de fréquences. On peut aussi utiliser les déformations partielles à la place des déplacements. [0230] On peut également établir des cartographies des paramètres de propagation, i.c. des images d'élasticité, avec des ondes de cisaillement différentes. Il est alors possible de les combiner, par exemple en les moyennant, de façon à obtenir une cartographie plus précise. [0231] La vitesse de propagation de l'onde de cisaillement dans le milieu 2 est suffisamment grande pour considérer qu'il n'y a pas de variation de la pression externe P exercée sur la surface externe 3 pendant cc processus durant lequel on prend une pluralité d'images du milieu (images intermédiaires) destinées à déterminer une image d'élasticité TE (paramètre de propagation de l'onde de cisaillement). Ainsi, ces images du milieu ne sont pas corrigées par les mouvements du milieu. [0232] Par contre, l'image d'élasticité du milieu IL doit être corrigée avec le mouvement à l'instant temporel considéré ou à un instant temporel proche (représenté par l'indice k) pour pouvoir être comparée à l'image du milieu initial Io, ou pour pouvoir être comparée à l'image de déformation corrigée 'pic*. [0233] Détermination d'une image d 'élasticité corrigée I E k * [0234] L'image d'élasticité TEk est alors corrigée dans un procédé représenté en figure 4, dans lequel on cherche à ramener cette image d'élasticité TEk dans une situation conforme à la première image du milieu lo pour supprimer les mouvements provoqués par la déformation de la pression externe P. [0235] Cette correction est effectuée en utilisant à chaque image d'élasticité IEk d'indice d'image k, l'ensemble des mouvements globaux déterminés entre chaque image du milieu. Plus précisément, il faut en appliquer les mouvements inverses pour revenir dans la situation de la première image L. [0236] Autrement dit, le procédé comprend les étapes suivantes : [0237] (el) détermination d'une image d'élasticité (IEk), [0238] (e2) la détermination d'une image d'élasticité corrigée (IEk'') par compensation de l'image d'élasticité (IL) à partir des mouvements précédents associés aux images du milieu d'indices compris entre 1 et k inclus. [0239] Détermination d'un e ima e d 'un aramètre de non linéarité INL k [0240] Dans un domaine élastique linéaire, le milieu 2 se déforme proportionnellement à la contrainte o et suit la loi de Hooke : [0241] = Eo. E [0242] où [0243] ED est le module d'Young du domaine linéaire, et [0244] s est la déformation. 19 [0245] Dans un domaine élastique non linéaire, cette proportionnalité n'est plus valable. Dans la plupart de cas, le module d'Young E du matériau du milieu 2 augmente avec la compression. [0246] On définit alors le module d'élasticité ou module d'Young E comme étant la pente de la courbe contrainte-déformation, c'est-à-dire : [0247] ao- (eq. 1) E = 3k£ [0248] où [0249] E est le module d'Young, [0250] An est la variation de contrainte locale, c'est-à-dire la contrainte partielle, [0251] As est la variation de déformation locale, c'est-à-dire la déformation partielle. [0252] Un paramètre de non linéarité de l'élasticité peut par exemple être le coefficient de Landau appelé le module d'élasticité en cisaillement du troisième ordre A mentionné dans les documents : [0253] - « Acoustoelasticity in soft solids : Assesment of the non-linear shear modulus with the acoustic radiation force », Gennisson et Al., J. Acoust. Soc. Am (122), December 2007, p. 3211-3219, et [0254] - « Quantitative Imaging of Vonlinear Shear Modulus by Combining Static Elas- tography and Shear Wave Elastography », H. Latorre-Ossa et Al., IEEE Transactions on Ultrasonics, Ferroelectrics, and Frequency Control, Vol. 51, no. 4, p. 833-839. [0255] Notamment, l'équation (1) de ce dernier document peut être réécrite sous la forme d'une relation (R1) entre l'élasticité (module d'Young) E et la contrainte o : [0256] A (eq. 2) E = Eo . 3.E0 [0257] où [0258] E est le module d'Young ou module d'élasticité, [0259] E0 est le module d'Young du domaine linéaire, c'est-à-dire le module d'Young du matériau du milieu sans compression, donc la première image d'élasticité corrigée dans la succession temporelle des images, [0260] o est la contrainte et par exemple la contrainte a, uniaxiale dans la direction Z, verticale, correspondant sensiblement à la compression induite par la pression exercée externe par l'utilisateur de la sonde 6, et [0261] A est le paramètre de non linéarité recherché. [0262] Par différentiation de la relation RI précédente, on obtient : [0263] A 21E = - âo- -3.ifo [0264] En divisant par la définition du module d'Young de l'équation (eq. 1) et par in- tégration on obtient une relation R2 entre l'élasticité (module d'Young) E et la déformation e. soit : [0265] A (eq. 3) ln (E) = ln (Ho) 4 E 3.E0 [0266] où [0267] E est le module d'Young, [0268] E, est le module d'Young du domaine linéaire, c'est-à-dire le module d'Young du matériau du milieu sans compression, donc la première image d'élasticité corrigée dans la succession temporelle des images, [0269] ln() est la fonction de logarithme népérien, [0270] e est la déformation, et [0271] A est le paramètre de non linéarité à déterminer. [0272] Selon une première variante du procédé NLa, représenté en figure 5, basé sur des images de déformation et des images d'élasticité, le procédé mis en oeuvre par le micro-ordinateur 4 comprend en outre les étapes suivantes : [0273] (NL1a) la détermination d'une pluralité d'images de déformation corrigée (Iekt), comme présenté précédemment, [0274] (NL2a) la détermination d'une pluralité d'images d'élasticité corrigée (IEkl, comme présenté précédemment, lesdites images de déformation corrigé et lesdites images d'élasticité étant temporellement entrelacées, et [0275] (NL3a) la détermination d'une image d'un paramètre de non linéarité (INLk) à partir de la pluralité d'images de déformation corrigée (Iek*) et de la pluralité d'images d'élasticité corrigée (jEk"). [0276] En fait, un paramètre de non linéarité d'un pixel de l'image est déterminé par l'ensemble (la pluralité) des valeurs de déformation corrigée et d'élasticité corrigée de ce même pixel. [0277] Par exemple, ce paramètre de non linéarité est déterminé par une régression linéaire de cette pluralité de couples de valeurs. Nous appellerons cette régression linéaire, une régression linéaire en déformation. [0278] Chaque couple de valeurs est formé pour un pixel de l'image et comprend : [0279] - une première valeur de ce même pixel d'une image de déformation corrigée, et [0280] - une deuxième valeur de ce même pixel d'une image de d'élasticité corrigée, [0281] les images de déformation corrigée et d'élasticité considérées étant des images tem- porellement suspensives ou temporellement proches. Ces images peuvent présenter des indices d'image identiques ou proche. 21 [0282] Eventuellement, les images de déformation et les images d'élasticité peuvent avoir des indices d'image différents et il est nécessaire de prendre des images d'indices correspondant à des instants temporels successifs ou temporellement proches comme explicités ci-dessus, c'est-à-dire que l'on effectue une correspondance entre le temps et l'indice de chaque image considérée. [0283] La figure 6 est une représentation pour un pixel d'image d'une pluralité de points Pa, chaque point Pa correspondant à un tel couple de valeurs (Ek, Fil Tous ces points Pa sont sensiblement alignés selon une droite Da dans un graphique ayant pour abscisse la valeur de déformation corrigée e (pour ce pixel) et ayant pour ordonnée le logarithme népérien de la valeur du module d'Young E (pour ce même pixel de l'image). [0284] La régression linéaire en déformation utilise alors la relation R2 entre l'élasticité (module d'Young) E et la déformation e (cg. 3) : In (E) = ln (E0) 4 A E .E [0286] où plus précisément dans le cas présent d'application : [0287] E est le module d'Young du pixel de l'image d'élasticité corrigée considérée, [0288] E1) est le module d'Young du pixel de la première image d'élasticité corrigée, [0289] ln() est la fonction de logarithme népérien, [0290] c est la déformation du pixel de l'image de déformation corrigée considérée, et [0291] A est le paramètre de non linéarité déterminé dudit pixel par ladite régression linéaire en déformation. [0292] La résolution de cette régression linéaire en déformation pour un ensemble de couples de valeurs extraits des pluralités d'images de déformation corrigée (Iek*) et d'images d'élasticité corrigée (IEkt) donne les valeurs du paramètre de non linéarité A de chaque pixel de l'image. [0293] Grâce à ces calculs, on obtient une image d'un paramètre de non linéarité INLk du milieu 2 qui utilise des images corrigées. Cette image ne présente pas les artéfacts d'images non corrigées de l'art antérieur. Cette image n'est pas ou peu influencée par les déformations et mouvements induits par la pression externe P. Les valeurs de paramètre de non linéarité A obtenues sont plus exactes et plus précises qu'avec les méthodes de l'art antérieur. [0294] Selon une seconde variante du procédé NLb, représenté en figure 7A et 7B, basé sur des images de contrainte et des images d'élasticité, le procédé mis en oeuvre par le micro-ordinateur 4 comprend : [0295] - une première phase durant laquelle des images de déformation partielles sont. converties en images de contrainte corrigées, et [0285] [0296] - une deuxième phase durant laquelle une image d'un paramètre de non linéarité est déterminée. [0297] Durant la première phase représentée en figure 7A, on effectue : [0298] (cl) la détermination d'une succession temporelle d'images de contrainte partielle corrigée (lAok*) par une multiplication de de déformation partielle corrigée (lApi(*) d'indice k et de l'image d'élasticité corrigée (TEC) d'indice k, et [0299] (c2) la détermination d'une image de contrainte corrigée (lok*) par une somme des images de déformation partielle corrigée (IAnC) des indices entre 1 et k inclus. [0300] Par multiplication d'images, on comprend ici que les valeurs de pixels de même co- ordonnée spatiale desdites images sont multipliées (i.e. multipliée l'une avec l'autre) pour former la valeur du pixel correspondant de l'image résultante. Ici, on multiplie donc une valeur de déformation partielle corrigée Ar', calculée par les étapes du procédé explicité précédemment et illustré en figure 3, par une valeur d'élasticité corrigée F.- (module d'Young corrigé), calculée par les étapes du procédé explicité précédemment et illustré en figure 4, pour obtenir la valeur de contrainte corrigée Ace dudit pixel. [0301] Ainsi, la contrainte partielle est calculée en chaque pixel de l'image par : [0302] M* = . As* 103031 Par somme ou cumul d'images, on comprend que les valeurs de pixels de même co- ordonnée spatiale desdites images sont sommées (i.e ajoutées les unes aux autres) pour former la valeur du pixel correspondant de l'image résultante. [0304] C'est-à-dire que : [0305] arc (x, z) = ​ (x z) i=i [0306] Cette formule est appliquée en tout pixel de l'image pour construire l'image de contrainte corrigée(Icek*). [0307] Durant la deuxième phase représentée en figure 7B, on effectue : [0308] (NLb1) la détermination d'une pluralité d'images de contrainte corrigée (Iok*), [0309] (NLb2) la détermination d'une pluralité d'images d'élasticité corrigée (IEkt), lesdites images de déformation corrigé et lesdites images d'élasticité étant temporellement entrelacées, et [0310] (NLb3) la détermination d'une image d'un paramètre de non linéarité (11.11_,,,) à partir de la pluralité d'images de contrainte corrigée 04'1 et de la pluralité d'images d'élasticité corrigée (IEk"). [0311] En fait, un paramètre de non linéarité d'un pixel de l'image est déterminé par l'ensemble (la pluralité) des valeurs de contrainte corrigée et d'élasticité corrigée de ce 23 même pixel. [0312] Par exemple, ce paramètre de non linéarité est déterminé par une régression linéaire de cette pluralité de couples de valeurs. Nous appellerons cette régression linéaire, une régression linéaire en contrainte. [0313] Chaque couple de valeurs est formé pour un pixel de l'image et comprend : [0314] - une première valeur de ce même pixel d'une image de contrainte corrigée, et [0315] - une deuxième valeur de cc même pixel d'une image de d'élasticité corrigée, [0316] les images de contrainte corrigée et d'élasticité considérées étant des images tempo- rellement suspensives ou temporellement proches. Ces images peuvent présenter des indices d'image identiques ou proche. [0317] Eventuellemcnt, les images de contrainte et les images d'élasticité peuvent avoir des indices d'image différents et il est nécessaire de prendre des images d'indices correspondant à des instants temporels successifs ou temporellement proches comme explicités ci-dessus, c'est-à-dire que l'on effectue une correspondance entre le temps et l'indice de chaque image considérée. [0318] La figure 8 est une représentation pour un pixel d'image d'une pluralité de points Ph, chaque point Ph correspondant à un couple de valeurs (Ek, uk). Tous ces points Pb sont sensiblement alignés selon une droite Db dans un graphique ayant pour abscisse la valeur de contrainte corrigée o (pour ce pixel) et ayant pour ordonnée la valeur du module d'Young E (pour ce même pixel de l'image). [0319] La régression linéaire en contrainte utilise alors la relation RI entre l'élasticité (module d'Young) E et la contrainte ci (eq. 2) : [0320] A E0 u [0321] où plus précisément dans le cas présent d'application : [0322] E est le module d'Young du pixel de l'image d'élasticité con-igée considérée, [0323] E0 est le module d'Young du pixel de la première image d'élasticité colligée, [0324] ln() est la fonction de logarithme népérien, [0325] E est la déformation du pixel de l'image de déformation corrigée considérée, et [0326] A est le paramètre de non linéarité déterminé dudit pixel par ladite régression en contrainte. [0327] La résolution de cette régression linéaire en contrainte pour un ensemble de couples de valeurs extraites des pluralités d'images de contrainte corrigée (Iakt) et d'images d'élasticité corrigée (IEkt) donne les valeurs du paramètre de non linéarité A de chaque pixel de l'image. [0328] Grâce à ces calculs, on obtient une image d'un paramètre de non linéarité INLk du 24 milieu 2 qui utilise des images corrigées. Cette image ne présente pas les artefacts d'images non corrigées de l'art antérieur. Cette image n'est pas ou peu influencée par les déformations et mouvements induits par la pression externe P. On obtient des valeurs de paramètre de non linéarité À plus exacts et plus précis que dans l'art antérieur. [Revendication 1] [Revendication 2] [Revendication 3] [Revendication 4] [Revendication 5] REVENDICATIONSProcédé pour déterminer par ultrasons une image corrigée d'un milieu, le procédé comprenant : - la détermination d'une succession temporelle d'images du milieu (1k), k étant un indice d'image compris entre 0 et N, et la première image de ladite succession étant supposée sans déformation, - la détermination du mouvement pour chaque image du milieu d' indice k entre 1 et N, le mouvement étant entre une image du milieu d'indice k (Ik) et l'image du milieu précédente d' indice k-1 (Ik ,), - la détermination d'une succession temporelle d'images de déformation partielle (iAek) à partir de l'image du milieu d'indice k (Ik) et l'image du milieu précédente d'indice k-1 - la détermination d'images de déformation partielle corrigée (Iàek') par compensation desdites images de déformation partielle (Iàek) à partir des mouvements précédents, c'est-à-dire les mouvements associés aux images du milieu d'indices compris entre 1 et k inclus, et - la détermination d'une image de déformation corrigée (kit) par une somme des images de déformation partielle corrigée (Iàek*) des indices entre 1 et k inclus. Procédé selon la revendication 1 dans lequel la détermination du mouvement à l'indice k est effectuée par : - détermination d'un champ de déplacement (uk) entre une image du milieu (Ik) et l'image du milieu précédente (Ik 1), - détermination d'une transformation géométrique (Tk) d'image à partir dudit champ de déplacement (uk), ladite transformation géométrique représentant le champ de déplacement (uk) avec un nombre de paramètres inférieur à dix. Procédé selon la revendication 2, dans lequel le champ de déplacement est calculé par intercorrélation entre une image du milieu (&) et l'image du milieu précédente (Ik-1), ou par un algorithme de suivi de sous-images entre une image du milieu (Ik) et l'image du milieu précédente (W-1). Procédé selon la revendication 3, dans lequel l'algorithme de suivi est un algorithme de Lucas-Kanade. Procédé selon l'une des revendications 2 à 4, dans lequel la transformation géométrique (Tk) comprend au moins une translation, ou une translation et une homothétie, ou une translation, une homothétie et une 26 [Revendication 6] [Revendication 7] [Revendication 8] [Revendication 9] [Revendication 10] rotation. Procédé selon l'une des revendications 2 à 4, dans lequel la transformation géométrique (Tk) comprend une translation, une homothétie et une rotation, et dans lequel la transformation géométrique est mise sous une forme matricielle telle que : os0 -Hz. sine = sine Hz. cosP T:1 T z 0 0 1 avec les paramètres de la transformation géométrique (Tk) dans le plan d'image suivants : Tx, Tz coefficients de translation, Hx, Hz coefficient d'homothétie, et 0 angle de rotation d'axe perpendiculaire au plan d'image. Procédé selon la revendication 2, clans lequel des paramètres de la transformation géométrique (Tk) sont. obtenus par des valeurs médianes d'une population de paramètres calculés à partir de groupes de points du champ de déplacement (uk). Procédé selon la revendication 7, clans lequel : - un groupe comprend entre trois et dix points du champ de déplacement (uk), et - la population est supérieure à cent groupes. Procédé selon la revendication 1, comprenant en outre : (el) la détermination d'une image d'élasticité (TEk) du milieu. et (e2) la détermination d'une image d'élasticité corrigée (IEk*) par compensation de ladite image d'élasticité (IEk) à partir des mouvements précédents, c'est-à-dire les mouvements associés aux images du milieu d'indices compris entre 0 et k inclus. Procédé selon la revendication 9, dans lequel l'image d'élasticité (1E10 produite à l'étape (el) ) est générée par les sous-étapes suivantes : (e1.1) une étape d'excitation au cours de laquelle on génère une onde de cisaillement dans le milieu, en faisant émettre au moins une onde ultrasonore focalisée, (e1.2) une étape d'observation au cours de laquelle on observe la propagation de l'onde de cisaillement en déterminant une succession temporelle d'images intermédiaires du milieu (11J,k), j étant un indice d'image intermédiaire compris entre 0 et M, (e1.3) une étape de traitement au cours de laquelle on détermine l'image 27 [Revendication 11] [Revendication 12] d'élasticité (In() à partir desdites images intermédiaires du milieu (IT;k) et d'un modèle de propagation d'onde de cisaillement. Procédé selon la revendication 9, dans lequel on effectue en outre : - la détermination d'une pluralité d'images de déformation corrigée (Tr.k *). - la détermination d'une pluralité d'images d'élasticité corrigée (TE:), lesdites images de déformation corrigé et lesdites images d'élasticité étant temporellement entrelacées, et - la détermination d'une image d'un paramètre de non linéarité (TNLk) à partir de la pluralité d'images de déformation corrigée (TEM et de la pluralité d'images d'élasticité corrigée Procédé selon la revendication 11, dans lequel la valeur de chaque pixel de l'image d'un paramètre de non linéarité (TNLk) est déterminée par régression linéaire en déformation de couples de valeurs, la première valeur de chaque couple correspondant à la valeur d'un même pixel d'une image de déformation corrigée de la pluralité et la deuxième valeur du couple correspondant à la valeur d'un même pixel d'une image d'élasticité de la pluralité, ladite image de déformation corrigée de la pluralité et ladite image d'élasticité corrigée de la pluralité étant temporellement successives ou proches, et dans lequel la régression linéaire en déformation est établie sur la base de la relation entre la déformation et l'élasticité suivante : A ln (E) = ln (E0) 4 . E 3 'E° [Revendication 13] où E est le module d'Young du pixel de l'image d'élasticité corrigée considérée, E0 est le module d'Young du pixel de la première image d'élasticité corrigée, ln() est la fonction de logarithme népérien, s est la déformation du pixel de l'image de déformation corrigée considérée, et A est le paramètre de non linéarité déterminé dudit pixel par ladite ré- gression linéaire en déformation. Procédé selon la revendication 9, dans lequel on effectue en outre : - la détermination d'une succession temporelle d'images de contrainte 28 [Revendication 14] [Revendication 15] partielle corrigée (iAc,") par une multiplication de l'image de déformation partielle corrigée (IAFI(*) d'indice k et de l'image d'élasticité corrigée (fEk*) d'indice k, - la détermination d'une image de contrainte corrigée (Iol,*) par une somme des images de déformation partielle corrigée (1-Anij des indices entre 1 et k inclus. Procédé selon la revendication 13, dans lequel on effectue en outre : - la détermination d'une pluralité d'images de contrainte corrigée (ink*), - la détermination d'une pluralité d'images d'élasticité corrigée (IE:), lesdites images de déformation corrigé et lesdites images d'élasticité étant temporellement entrelacées, et - la détermination d'une image d'un paramètre de non linéarité (JNLk) à partir de la pluralité d'images de contrainte corrigée (irk*) et de la pluralité d'images d'élasticité corrigée (IEL'). Procédé selon la revendication 14, dans lequel la valeur de chaque pixel de l'image d'un paramètre de non linéarité (JNLk) est déterminée par régression linéaire en contrainte de couples de valeurs, la première valeur de chaque couple correspondant à la valeur d'un même pixel d'une image de contrainte corrigée de la pluralité et la deuxième valeur du couple correspondant à la valeur d'un même pixel d'une image d'élasticité de la pluralité, ladite image de contrainte corrigée de la pluralité et ladite image d'élasticité colligée de la pluralité étant temporellement successives ou proches, et dans lequel la régression linéaire en contrainte est établie sur la base de la relation entre la contrainte et l'élasticité suivante : A 1 (E) = In (E0) . E Eo où E est le module d'Young du pixel de l'image d'élasticité corrigée considérée, E0 est le module d'Young du pixel de la première image d'élasticité corrigée, ln() est la fonction de logarithme népérien, s est la déformation du pixel de l'image de déformation corrigée considérée, et A est le paramètre de non linéarité déterminé dudit pixel par ladite ré- 29 [Revendication 16] gression. Dispositif d'imagerie (1) comprenant une sonde ultrasonore (6) et un micro-ordinateur (4) adaptés pour mettre en oeuvre le procédé pour déterminer par ultrasons une image corrigée d'un milieu (2) selon l'une des revendications 1 à 15.