Brevet FR3104736A1 — publié le 20210618

Procédé ultrasonore pour quantifier l’élasticité non linéaire par ondes de cisaillement d’un milieu, et dispositif pour mettre en œuvre ce procédé Procédé ultrasonore pour quantifier l’élasticité non linéaire par ondes de cisaillement d’un milieu, le procédé comprenant les étapes suivantes : A1. - la collecte d’une succession temporelle de données d’élasticité par ondes de cisaillement du milieu, A2. - une application d’une déformation successivement changeante au milieu selon une séquence de déformations prédéterminée pendant la collecte des ondes de cisaillement, A3. - l’observation l’évolution de déformation réelle, et B. - la quantification de l’élasticité non linéaire du milieu en fonction de la succession temporelle de données et l’évolution de déformation. Figures pour l’abrégé : Fig. 2A-2C. Procédé ultrasonore pour quantifier l’élasticité non linéaire par ondes de cisaillement d’un milieu, et dispositif pour mettre en œuvre ce procédé La présente invention 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 quantifier la non-linéarité du milieu, plus notamment pour détecter une région cancéreuse molle dans le milieu (c’est-à-dire non-indurée en comparaison de milieu non-cancéreux). Plus particulièrement, la présente invention 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. Des images ultrasonores peuvent être obtenues par des procédés différents. Par exemple, les images obtenues peuvent être de type élasticité statique (« strain » en langue anglaise), obtenues par le biais de déformation des tissus observés, ou de type élasticité par ondes de cisaillement (« shear wave elastography, SWE » en langue anglaise). Il est connu de détecter des cancers indurés par une imagerie ultrasonore de l’élastographie du milieu par ondes de cisaillement (« shear wave elastography, SWE » en langue anglaise). Cependant, cette technique ne permet pas de détecter des cancers mous, car les cancers mous présentent une élasticité similaire aux milieux non-cancéreux sur une image l’élastographie classique. D’ailleurs, il est connu de déterminer des caractéristiques de non-linéarité d’un milieu, cf. par exemple : H. Latorre-Ossa, JL Gennisson, E. De Brosses, M Fink, Quantitative imaging of nonlinear shear modulus by combining static elastography and shear wave elastography, IEEE Trans Ultrason Ferroelectr Freq Control. 2012 Apr; 59(4):833-9, et M. Bernal, F. Chamming’s, M Couade, J Berfcoff, M. Tanterm JL Gennisson, In Vivo Quantification of the Nonlinear Shear Modulus in Breast Lesions: Feasibility Study. IEEE Trans Ultrason Ferroelectr Freq Control. 2016 Jan; 63(1):101-9. doi: 10.1109/TUFFC.2015.2503601. Epub 2015 Nov 24. Cependant, pour extraire l’information de la non-linéarité, plusieurs problèmes doivent être résolus, en particulier trouver une expression mathématique qui représente la non-linéarité du module de cisaillement comme une fonction du module de cisaillement et de la déformation, de déterminer/mesurer des valeurs qui sont nécessaires pour calculer cette expression mathématique, mettre en place un procédé de mesure et calculs permettant l’obtention de cette information. La présente invention a donc pour but d’améliorer la détection et le diagnostic du cancer non-induré. A cet effet, la présente invention propose un procédé pour quantifier l’élasticité non linéaire par ondes de cisaillement d’un milieu, le procédé comprenant les étapes suivantes : A1. la collecte d’une succession temporelle de données d’élasticité par ondes de cisaillement du milieu, A2. une application d’une déformation successivement changeante au milieu selon une séquence de déformations prédéterminée pendant la collecte des ondes de cisaillement, A3. l’observation de l’évolution de déformation réelle (estimée), et B. la quantification de l’élasticité non linéaire du milieu en fonction de la succession temporelle de données et l’évolution de déformation. Grâce à ces dispositions, la détection ou le diagnostic du cancer mou, par exemple du sein, par image ultrason peut être amélioré. Les cancers mous présentent une élasticité similaire aux tissus sains ou aux lésions bénignes molles sur une image élastographie standard SWE (sans compression). Or il a été découvert qu’en comprimant les tissus (milieux) pendant une élastographie SWE l’élasticité varie différemment sur les tissus sains et sur les lésions malignes. L’imagerie NL-SWE a pour but de qualifier et quantifier ce phénomène. Par exemple, en se basant sur la théorie de l’acousto-élasticité, la non linéarité du module de cisaillement du tissus (« share modulus non linearity » en anglais) ou le coefficient élastique non linéaire de cisaillement ou le durcissement (« stiffening » en anglais) peuvent être déduits des images « strain » et SWE. Pour être précis, dans ce mode, en compressant ou en décompressant le tissu avec la sonde, l'opérateur humain (ou une machine, ou n’importe quoi qui peut tenir ou manipuler une sonde ultrason), peut acquérir une série d’images « SWE » et « strain » utilisées pour déduire une carte qui représente la non linéarité du module de cisaillement ou le taux de durcissement du milieu. Par ailleurs, la présente invention peut résoudre en particulier les problèmes suivants, notamment avec l’objectif d’extraire l’information de la non linéarité : 1. Trouver une expression mathématique qui représente la non linéarité du module de cisaillement comme une fonction du module de cisaillement et de la déformation. 2. Développer une séquence ultrason qui peut acquérir simultanément des images mode B, des images d’élastographie par ondes de cisaillement (SWE) et des images d’élastographie strain (par déformation). 3. Développer une technique pour agréger et recaler les images « strain » et SWE pendant la compression ; ou en fin de compression ou lors d’étapes dédiées prévues durant le processus de compression. 4. Développer une interface utilisateur qui guide l'opérateur pour la compression/décompression. 5. Créer un accessoire pour la sonde qui permet de guider la compression/décompression. Dans divers modes de réalisation du procédé selon la présente invention, on peut éventuellement avoir recours en outre à l'une et/ou à l'autre des dispositions suivantes. Selon un aspect, l’observation l’évolution de déformation réelle est effectuée en temps réel ou en quasi temps réel. Ainsi, grâce à cette observation en temps réel ou en quasi temps réel, la déformation peut être effectuée de maniere continue et progressive (c’est-à-dire sans pause). Selon un aspect, la déformation successivement changeante est appliquée par l’intermédiaire d’une sonde ultrasonore en successivement compressant ou en successivement décompressant le milieu par la sonde, la sonde étant simultanément utilisée pour collecter la succession temporelle de données et/ou pour observer l’évolution de déformation réelle. Cette méthode peut faire appel à un algorithme d’estimation du mouvement, incluant (mais non limité à), une méthode basée sur l’effet Doppler, le tracking du speckle, et/ou les flux optiques. L’observation de l’évolution de déformation réelle peut être représentée ou mesurée par le « strain » du milieu (comme décrit ci-dessous). Selon un aspect, la compression successive et/ou la décompression successive est effectuée de maniere continue et progressive. Par exemple, la compression et/ou la décompression est un mouvement continu sans pause ou arrêt. Selon un aspect, l’étape d’observation de l’évolution de déformation réelle comprend la collecte d’une succession temporelle de données ultrasons du milieu, notamment pour permettre des calculs de déformation (« strain »), par exemple une collecte d’une succession temporelle de données d’élasticité statique (« strain »). Selon un aspect, un niveau de déformation réelle estimée et un niveau de déformation cible sont notifiés sur une interface utilisateur, dans lequel la déformation réelle estimée correspond à la déformation observée, notamment en temps-réel, et la déformation cible change (notamment de manière continue et progressive) selon la séquence de déformation prédéterminée, notamment en temps-réel. Il est noté que le niveau de déformation réelle peut être seulement estimé, c’est-à-dire ce niveau peut avoir une certaine imprécision. Selon un aspect, le niveau de déformation réelle estimé et le niveau de déformation cible sont affichés sur un écran de l’interface utilisateur. Selon un aspect, les données d’élasticité par ondes de cisaillement produites à l’étape (A1) sont générées par les sous-étapes suivantes : A1.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, A1.2. une étape d'observation au cours de laquelle on observe la propagation de l'onde de cisaillement en acquérant une succession temporelle de données ultrasons du milieu, A1.3. une étape de traitement au cours de laquelle on détermine les données d’élasticité à partir desdites données ultrasons du milieu et d’un modèle de propagation d’onde de cisaillement. Ce procédé peut utiliser des ondes planes mais cette étape peut aussi être réalisée par utilisation d’ondes focalisées. Selon un aspect, la collecte d’une succession temporelle de données et l’observation des changements successifs de la déformation sont mises en place par une séquence ultrasons en intercalant des tirs l'élasticité par ondes de cisaillement et des tirs d’élasticité statique, et/ou lesdites données d’élasticité statiques et lesdites de données d’élasticité par ondes de cisaillement sont temporellement entrelacées. Selon un aspect, l’élasticité statique peut être utilisée uniquement pour déduire les mesures de déformation. Selon un aspect, le procédé peut comprendre en outre l’étape C. de détermination d’une image et/ou un score et/ou un symbole montrant le niveau d’élasticité non linéaire du milieu. Le score et/ou symbole peut être un icone, une valeur versus un seuil, des phrases, un tableau de valeurs, etc. Selon un aspect, le procédé peut comprendre en outre les étapes suivantes : A4. la collecte d’une succession temporelle de données mode B du milieu pour calculer une compensation de mouvement du milieu relative à la sonde, A5. l’application de la compensation aux données d’élasticité. Selon un aspect, lesdites données mode B, lesdites données d’élasticité statique et lesdites données d’élasticité par ondes de cisaillement sont temporellement entrelacées. Selon un aspect, les étapes (A1) à (A3) ou (A1) à (A5) sont effectuées simultanément, et/ou les étapes (A1) à (A3) ou (A1) à (A5) sont effectués en temps-réel ou en quasi temps-réel. Selon un aspect, pour chaque instant de la succession temporelle de données ou au moins pour quelques-uns des instants, la valeur d’un paramètre de non linéarité des ondes de cisaillement est déterminée en effectuant un ajustement de régression linéaire selon l’équation : dans laquelle : E(t) est l’élasticité à chaque instant de déformation lors de l’évolution de déformation, E 0 est l’élasticité à un premier instant (t0), ε est la déformation, et A est le paramètre de non linéarité de cisaillement. Selon un aspect, le procédé comprend en outre une étape (0) préalable à l’étape (A1), dans laquelle au moins une image ultrasons ou IRM (imagerie par résonance magnétique) ou mammographie ou image d’élasticité par ondes de cisaillement du milieu est déterminée, sans application d’une déformation selon une séquence de déformation prédéterminée, et cette étape permet de déterminer une zone d’intérêt (ROI region of interest) pour laquelle on quantifie l’élasticité non linéaire selon les étapes (A1) à (B). Selon un aspect, une carte de confiance comportant des niveaux de fiabilités est affichée en parallèle de l’affichage de l’image d’élasticité non linéaire, afin de permettre en particulier au praticien de connaitre le niveau évalué de confiance de chacun des points de l’image « résultat » présentée. Cette carte de confiance peut être une carte (ou une image) qui indique le niveau de fiabilité de chaque partie de l’image montrant le niveau d’élasticité non linéaire du milieu. Le niveau de fiabilité peut être déduit par exemple de statistiques issues du résultat de la régression linéaire et/ou peut tenir compte des caractéristiques du nouvement réel appliqué sur les tissus (vitesse uniforme, mouvement uniaxial, ….). Cette carte peut être implémentée au moyen de matrices, images, codes couleurs, et/ou tableau. Selon un aspect, pour chaque donnée ou pour chaque pixel de l’image, une valeur d’élasticité non linéaire en regard d’une valeur de référence prédéfinie est affichée. Cette valeur de référence prédéfinie peut être une valeur prédéfinie dans la littérature scientifique et/ou peut être calculée pour ce même patient lors d'un précédent examen, ou le résultat de calculs statistiques et optionnellement affichée graphiquement. Par exemple, la valeur peut être affichée à la place ou en plus de la valeur d’élasticité non linéaire (par exemple pour chaque pixel ou groupe de pixels). L’affichage de la valeur prédéfinie peut avantageusement permettre une analyse et/ou un diagnostic fiable, par exemple d’un milieu cancéreux mou et/ou indiquer le type de cancer observé. La présente invention a également pour objet un procédé ultrasonore pour quantifier l’élasticité non linéaire par ondes de cisaillement d’un milieu en utilisant un dispositif d’imagerie (1) comprenant une sonde ultrasonore (6), un micro-ordinateur (4), et optionnellement une interface utilisateur, le procédé comprenant les étapes suivantes : positionner la sonde à la surface du milieu, activer un premier mode d’imagerie de l’élasticité non linéaire par ondes de cisaillement (NL-SWE), pendant lequel, la sonde appliquant une déformation selon une séquence de déformation prédéterminée, en progressivement compressant ou en progressivement décompressant le milieu : le dispositif collecte une succession temporelle de données ultrasons imageant les ondes de cisaillement générées dans le milieu à l’aide de la sonde permettant de déduire l’élasticité du milieu, et le dispositif collecte une succession temporelle de données ultrasons permettant de déduire une évolution de déformation du milieu appliqué par la sonde, quantifier l’élasticité non linéaire du milieu, qui est calculée en fonction de la succession temporelle de données et l’évolution de déformation observée, optionnellement afficher sur l’interface utilisateur une image de l’élasticité non linéaire du milieu, optionnellement enregistrer la succession temporelle d’images, comprenant des calculs et mesurages prédéfinies des images, optionnellement afficher sur l’interface utilisateur la succession temporelle d’images, optionnellement afficher un résultat d’interprétation de l’image non linéaire. Chacune de ces images peut être de type bidimensionnel ou 3D. Avant de lancer ce mode NL-SWE, une etape de réglage des paramètres de la séquence de déformation prédeterminée, par exemple en fonction du cas clinique (taille du patient, taille de l‘organe considéré, en particulier lorsqu’il s’agit de tissus mammaires, profondeur de la lésion observée, densité des tissus, etc.) peut être mise en œuvre. En outre, après l'affichage de l'image de non linéarité, une étape de réglage des paramètres de l'affichage, par exemple en fonction des préférences de l'utilisateur (superposition avec d'autres images B mode, etc. ), changement des unités , changement du code couleurs, etc. Comme la compression (contrainte) ou décompression et un mouvement progressif, un utilisateur ou un bras articulé ou tout autre operateur, peut l’effectuer de manière continue, sans pause. La sonde peut collecter, lors de chacun des étapes de collecte du procédé, des données ultrasons sur lesquelles peut être appliquée une méthode d’estimation du mouvement (speckle tracking sur B Mode ou effet Doppler) afin d'extraire des informations relatives à la déformation et/ou afin de recaler les données ultrasons collectées. Selon un aspect, la séquence de déformation prédéterminée comporte une sous étape d’arrêt optionnellement automatique avant l’étape de quantification. Cela procure en particulier l’avantage pour l’opérateur d’être notifié par exemple que la session d’acquisition est terminée car nombre de données collectée est suffisante. Selon un aspect, le procédé peut comprendre en outre, notamment avant l’activation du premier mode : positionner la sonde à la surface du milieu, activer un deuxième mode d’imagerie d’élasticité par ondes de cisaillement (SWE), déterminer au moins une image d’élasticité par ondes de cisaillement du milieu en appliquant des ondes de cisaillement dans le milieu, sans application d’une déformation selon une séquence de déformation prédéterminée, enregistrer l’image, optionnellement comprenant des calculs et mesurages prédéfinies de l’image, déterminer une zone d’intérêt dans l’image pour imager l’élasticité non linéaire de la zone. En outre, la présente invention propose un procédé pour déterminer un milieu cancéreux non-induré en comparaison du milieu non-cancéreux (une lésion maligne non-induré en comparaison d’une lésion bénigne) en utilisant le procédé tel que cité ci-dessus, dans lequel le milieu cancéreux non-induré (lésion maligne non-induré ) est détecté en fonction de son niveau d’élasticité non-linéaire. En particulier, la présente invention propose un procédé pour déterminer un milieu cancéreux non-induré en comparaison du milieu non-cancéreux en utilisant le procédé tel qu’exposé ci-dessus, dans lequel le milieu étudié est le tissu biologique mammaire et/ou axillaire. La présente divulgation se rapporte également à un dispositif d’imagerie comprenant une sonde ultrasonore et un micro-ordinateur adaptés pour mettre en œuvre le procédé pour quantifier l’élasticité non linéaire d’un milieu tel que cité ci-dessus. La figure 1 est une vue schématique d’un dispositif d’imagerie par ultrasons selon une forme de réalisation de l’invention. Les figures 2A, 2B et 2C montrent de manière schématique un exemple du procédé selon la présente invention, notamment d’une application au milieu d’une déformation schématique successivement changeante. Les figures 3A, 3B et 3C illustrent de manière schématique l’écran d’une interface utilisateur du dispositif lors de l’application d’un procédé selon la présente divulgation. Sur les différentes figures, fournies à titre d’illustration, les mêmes références numériques désignent des éléments identiques ou similaires. 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. 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"). 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. La sonde ultrasonore 6 est constituée par un réseau de n transducteurs ultrasonores T1, T2, ..., Ti, ..., Tn, n étant un nombre entier supérieur ou 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 X perpendiculaire à l'axe Z. La sonde en question peut être également un réseau bidimensionnel (plan ou non) de transducteurs. Les transducteurs T1, 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 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, 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 T1-Tn, soit une onde ultrasonore de compression focalisée en un ou plusieurs points du milieu 2. 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. Le dispositif d’imagerie 1 et le procédé selon la présente invention effectuent une quantification l’élasticité non linéaire par ondes de cisaillement d’un milieu. Pour commencer ce procédé, on peut activer un premier mode d’imagerie de l’élasticité non linéaire par ondes de cisaillement (NL-SWE, « Non-Linear Shear Wave Elasticity » en langue anglaise). Le procédé du mode NL-SWE utilise un procédé ultrasonore pour déterminer l’élasticité du milieu par ondes de cisaillement (SWE). Avantageusement, ce procédé (c’est à dire la collecte des données respective ou le « scan ») et combiné simultanément et/ou intercalé avec un procédé ultrasonore pour observer en temps réel ou en quasi temps réel l’évolution de la déformation. Des exemples détaillés de ces trois procédés ultrasonores sont expliqués en détails ci-dessous. Les procédés SWE et strain peuvent être basés sur des ondes ultrasons planes, tandis que le procédé B-mode peut être basé sur des ondes ultrasons non-plans. Alternativement au procédé ultrasonore pour déterminer l’élasticité statique (« strain »), ou en addition, la sonde peut aussi être équipée avec un capteur de pression qui permet d’observer (ou suivre) l’évolution de déformation du milieu. Les figures 2A à 2C illustrent un exemple du procédé selon la présente invention, notamment d’une application d’une déformation schématique successivement changeante au milieu. Dans cet exemple, le procédé commence dans l’état illustré dans la figure 2A et se termine dans l’état illustré dans la figure 2C. Notamment, les figures 2A à 2C montrent l’évolution de déformation du milieu, causée par une application d’une déformation successivement changeante au milieu. Dans le même temps, le dispositif d’imagerie 1 effectue une quantification de l’élasticité non linéaire par ondes de cisaillement du milieu selon la présente invention. Le procédé peut durer quelques secondes, par exemple entre 5-10s. Ce temps est utile pour collecter assez de données lors de l’évolution de déformation du milieu (c’est à dire pour suffisamment quantifier l’élasticité non linéaire par ondes de cisaillement d’un milieu lors de cette évolution). Un compte à rebours peut être notifié à l’utilisateur, afin qu’il soit prêt à exercer la déformation des tissus par de biais de la sonde ou tout autre moyen. La déformation illustrée dans les figures 2A à 2C est seulement schématique. Le niveau de déformation qui est appliqué selon la présente divulgation est prédéterminé par une séquence de déformation prédéfinie, comme expliqué dans le contexte des figures 3A à 3C. La figure 2A représente une utilisation du dispositif d’imagerie 1 à un instant t0 du procédé, 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 I 0 comprend par exemple une inclusion 2i à une profondeur Z1 par rapport à la surface externe 3. En même temps, le procédé ultrasonore pour quantifier l’élasticité non linéaire par ondes de cisaillement du milieu (comme expliqué ci-dessus) démarre. La figure 2B représente une utilisation du dispositif d’imagerie 1 à un instant t1 du procédé, et avec une pression (ou contrainte) externe P’ exercée, supérieure à P, qui déforme la surface externe 3 dans la direction Z vers l’intérieur du milieu 2. Donc, la déformation du milieu a été successivement changée. En figure 2C est représentée une utilisation du dispositif d’imagerie 1 à un instant t2 du procédé, avec une pression externe P’’ encore plus élevée que P’ exercée à t1. La compression successive, comme illustrée dans les figures 2A à 2C, peut être effectuée de maniere continue et progressive. Cela est possible grâce à l’observation d’évolution de déformation réelle en temps réel ou en quasi temps réel. Alternativement une décompression successive peut être appliquée aux tissus composant le milieu. Comme déjà indiqué, l’évolution de la déformation montrée dans les figures 2A à 2C n’est pas aléatoire, mais est prédéterminée par une séquence de déformation prédéfinie. Par exemple, cette séquence peut être illustrée par exemple graphiquement à l’utilisateur sur l’écran 4a d’une interface utilisateur du dispositif 1, permettant ainsi à l’utilisateur d’appliquer une pression correcte à chaque moment du procédé. Par exemple, un affichage d’un curseur indiquant la pression à exercer ainsi que la pression appliquée mesurée sur un même graphe ou des vibrations dans la sonde pour indiquer à l’utilisateur s’il suit bine la séquence prédéfinie .Il est également possible que l’application de la déformation soit automatisée, par exemple en utilisant un bras robotique qui bouge la sonde d’une façon automatisée selon la séquence prédéterminée. Les figures 3A à 3C illustrent un exemple d’écran 4a d’une interface utilisateur du dispositif lors d’un procédé selon la présente divulgation. Par exemple, les figures 3A à 3C peuvent correspondre au instants temporels t0-t2 des figures 2A à 2C. L’écran 4a peut montrer un champ 41 pour illustrer au moins une image ultrasonore. Par exemple cette image peut être une image d’une zone d’intérêt du milieu (par exemple une image d’élasticité par ondes de cisaillement (SWE)) qui a été prise avant de débuter le procédé, c’est-à-dire avant t0. Après l’instant t2, le champ 41 peut montrer une image montrant le niveau d’élasticité non linéaire du milieu selon la présente invention. L’écran 4a montre en outre un champ 42. Ce champ 42 affiche pendant le procédé le niveau actuel de déformation réel estimé (par exemple avec une croix ou point O) et le niveau de déformation cible (par exemple avec un ballon ou cercle C). La déformation réelle estimée O correspond à la déformation actuellement observée, notamment en temps-réel, et la déformation cible C change successivement selon la séquence de déformation prédéterminée, notamment en temps-réel. La cible C peut évoluer successivement dans le champ 42 pendant le procédé, afin de montrer à chaque instant l’évolution de la séquence de déformation prédéterminée et donc, guide l’utilisateur. La figure 3A peut correspondre à l’état de la figure 2A et donc à l’instant temporel t0. Par conséquent, la sonde 6 exerçant une pression externe P faible ou nulle sur la surface externe 3. Donc, le niveau actuel de déformation réel estimé O et le niveau actuel de déformation cible C sont tous les deux à zéro. La figure 3B peut correspondre à l’état de la figure 2B et donc à l’instant temporel t1. Par conséquent la sonde 6 exerçant une pression externe P augmenté sur la surface externe 3. Dans ce scénario décrit à titre d’exemple, il y a une différence de niveau entre la déformation réelle estimée O et le niveau actuel de déformation cible C. Dans cet exemple, le niveau actuel de déformation réelle estimée O est moins élevé que le niveau actuel de déformation cible C. C’est-à-dire l’utilisateur est guidé par la cible C afin d’augmenter plus la pression. La figure 3C peut correspondre à l’état de la figure 2C et donc à l’instant temporel t2. Par conséquent, la sonde 6 exerçant une pression externe P’’ encore plus élevée que P’ à t1. Comme t2 correspond à la fin de la séquence de déformation prédéterminée, le niveau actuel de déformation réel estimé O et le niveau actuel de déformation cible C sont tous les deux au maximum de la séquence de déformation prédéterminée. Optionnellement, l’imagerie et arrêtée ou « bloquée » (« freeze » en anglais) à l’instant temporel t2. Cela permet avantageusement à l’utilisateur de se concentrer sur la séquence de déformation prédéfinie jusqu’ à t2, sans devoir arrêter l’imagerie manuellement. A l’instant t2 ou après, ou pendant l’acquisition des données, le dispositif 1 quantifie l’élasticité non linéaire du milieu, qui est calculée en fonction de la succession temporelle de données et l’évolution de déformation observée lors la période t0 à t2. L’interface utilisateur peut afficher sur le champ 41 une image de l’élasticité non linéaire du milieu, qui est déterminée en fonction de l’élasticité non linéaire quantifiée. Optionnellement le dispositif 1 peut enregistrer la succession temporelle de données (ou images) collectées lors la période t0 à t2. Cette étape peut comprendre des calculs et mesures prédéfinies des données. On outre, l’écran 4a peut afficher sur l’interface utilisateur la succession temporelle d’images. Après l'affichage de l'image, on peut avoir une étape de réglage des paramètres de l'affichage, en fonction des préférences de l'utilisateur (superposition avec d'autres d'images B mode, etc.), changement des unités, changement du code couleurs. Par exemple on peut changer entre un affichage sans unités, qui est plus facile à interpréter et un affichage avec des unités (par exemple kPa), qui peut être avantageux pour être utilisé avec de la littérature scientifique pour interpréter les valeurs. Ces analyses peuvent être réalisées par le biais d’algorithmes statistiques par tout moyen, par exemple intelligence artificielle. Par ailleurs, le dispositif peut impliquer un deuxième mode d’imagerie ultrason, et/ou d’imagerie d’élasticité par ondes de cisaillement (SWE), et/ou de mammographie. Ce mode peut ainsi être activé avant d’utiliser le mode NL-SWE ou lors d’une session indépendante préalable. Dans ce mode, le dispositif peut déterminer au moins une image d’élasticité par ondes de cisaillement du milieu en appliquant des ondes de cisaillement dans le milieu, sans application d’une déformation selon une séquence de déformation prédéterminée. L’image peut être enregistrée, optionnellement comprenant des calculs et mesures prédéfinies de l’image. Ensuite, l’utilisateur peut déterminer une zone d’intérêt (ROI region of interest) dans l’image pour imager l’élasticité non linéaire de la zone selon le premier mode NL-SWE. Détermination d’images du milieu Ik en « mode B » Le procédé décrit ci-dessous est connu comme l’imagerie en « mode B ». La succession temporelle de données (ou image ou données d’image) mode B peut être utilisée pour calculer une compensation de mouvement du milieu relative à la sonde lors de l’application d’une déformation successivement changeant (comme décrit ci-dessus). Pour générer une image du milieu (I) le dispositif d’imagerie effectue par exemple les étapes suivantes : (i1) 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, (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 (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 (I). 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 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 T1, ... , 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 transducteur, respectivement E1, E2, ... En. 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. 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. 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). 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. Par exemple, on peut déterminer un signal S(t) par la formule suivante : dans laquelle : si est le signal brut perçu par le transducteur n° i après le tir d'onde ultrasonore de compression, 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, 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, V est la vitesse moyenne de propagation des ondes acoustiques ultrasonores de compression dans le milieu viscoélastique observé, et αi(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 αi(x,z) = 1). 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). 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). D’autres techniques de génération d’image du milieu (I) existent, comme les techniques d’imagerie synthétique. Toutes technique d’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. Détermination de l’élasticité statique (« strain ») Le procédé décrit ci-dessous est connu comme l’imagerie d’élasticité statique (« strain »). La succession temporelle de données (ou image ou données d’image) d’élasticité statique (« strain ») peut être utilisée pour observer l’évolution de déformation réelle estimée (comme décrit ci-dessus). Les images du milieu (I 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 (I k ) et l’image du milieu du milieu d’indice k-1 (I k-1 ). 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. Au cours de ce processus d'intercorrélation, on maximise une fonction d'intercorrélation k-1 (x,z) ,S k (x,z) > afin de déterminer le déplacement subi 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 dans les documents suivants : "Internal 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 ) "Elastography: a quantitative method for imaging the elasticity of biological tissues", Ophir et al., Ultrasonic imaging., vol. 13, p.111-134, 1991), et "Pyramidal Implementation of the Lucas Kanade Feature Tracker. Description of the algorithm.", J-Y Bouguet, Intel Corp. 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 u k (x,z) en remplaçant la variable de l’instant temporel t par l’indice d’image k. Ces vecteurs 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 u k à l’indice d’image k est : Ce champ de déplacement (ensemble de vecteurs déplacements) est stocké dans la mémoire M ou dans le micro-ordinateur 4. On peut construire une image du champ de déplacement Iu k , et l’ensemble des images du champ de déplacement Iu k 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 gris ou par un niveau chromatique. Détermination d’images de déformation partielle I Δε k On peut alors calculer les déformations partielles Δε 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 Eventuellement, selon un exemple de mise en œuvre, seule la déformation/déformation partielle uni-axiale 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 uni-axiale dans la direction Z est éventuellement calculée. On simplifie ainsi les calculs. Comme pour les vecteurs déplacements, on peut construire une image de déformation partielle (IΔε k ) à partir de l’ensemble des déformations partielles calculées en tout point de coordonnées (x, z) du plan X-Z (plan d’image). Ces vecteurs de déformations partielles ou images de déformation partielles (IΔε k ) peuvent être visualisés, notamment au moyen de l'écran 4a du micro-ordinateur, sous la forme d'un film ralenti. Détermination du mouvement du milieu par la pression externe 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éplacement des particules 5 et des éléments inclus dans le milieu 2. 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 (I k ) et une image du milieu précédente d’indice k-1 (I k-1 ). Le mouvement peut alors être déterminé par le champ de déplacement (u k ) entre une image du milieu d’indice k et une image du milieu précédente d’indice k-1. Notamment, ce mouvement peut être modélisé par une transformation géométrique T k 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, u k (x, z), c’est-à-dire avec un nombre de paramètres très réduit, et par exemple inférieur à dix. Selon une première variante, la transformation géométrique T k comprend une translation, comme cela est représenté entre la figure 2A et la figure 2B. Selon une seconde variante, la transformation géométrique T k comprend une translation et une homothétie. Selon une seconde variante, la transformation géométrique T k comprend une translation, une homothétie et une rotation. Ainsi, dans le plan X-Z, plan d’image, la transformation géométrique peut être mise sous une forme matricielle du type : avec les paramètres de la transformation géométrique (Tk) suivants : Tx, Ty coefficients de translation dans le plan d’image, Hx, Hz coefficient d’homothétie le plan d’image, et θ angle de rotation d’axe perpendiculaire au plan d’image. Ainsi, si l’on prend trois points P1, P2, et P3 de l’image de cordonnées respectives (x1,z1), (x2,z2), (x3,z3), on a les relations suivantes : avec : dx1 = u x1 dz1 = u z1 dx2 = u x2 dz2 = u z2 dx3 = u x3 dz3 = u z3 pour reprendre les notation des composantes X et Z des vecteurs déplacements des points P1, P2 et P3. En outre, on a une relation permettant de calculer la transformation géométrique T k à 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 de l’inverse de la matrice A : L’application de cette relation permet de calculer la transformation géométrique T k à partir de trois points P1, P2, P3. Inversement, une fois la transformation géométrique T k établie, cette relation permet de connaitre les coordonnées (x, z) de tout point P par la relation inverse : 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 i , l’indice i variant de 1 à k. 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. 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. Alors, les paramètres de la transformation géométrique T k sont obtenus par des valeurs médianes des paramètres calculés à partir de la population des groupes de points de l’image. Notamment, cette technique peut être appliquée à une image du champ de déplacement u k 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. Grâce à l’utilisation d’une population de groupes de points, on peut déterminer une transformation géométrique qui représente de 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. 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 de manière plus globale et plus fiable le mouvement entre images. Détermination d’une image d’élasticité par ondes de cisaillement du milieu IE k Le procédé décrit ci-dessous est connu comme l’imagerie d’élasticité par ondes de cisaillement du milieu. La succession temporelle de données (ou image ou données d’image) d’élasticité par ondes de cisaillement peut être utilisé pour quantifier l’élasticité non linéaire du milieu (comme décrit ci-dessus). 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. On procède en plusieurs étapes : (e1.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, (e1.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 II j,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. 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 I k 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 (II j,k ). 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 II j,k générées. On procède alors à : (e1.3) une étape de traitement des images intermédiaires du milieu pour déterminer une image d’élasticité du milieu IE k . Les images intermédiaires du milieu II j,k sont traitées pour calculer un champ de déplacement u j,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 (II j,k ) et l’image intermédiaire du milieu du milieu d’indice j-1 (I j-1,k ). 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 Δε en chaque point du milieu 2 à partir du champ de déplacement. 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. On peut alors construire une image d’élasticité IE k , correspondant à l’ensemble des paramètres des de propagation de l’onde de cisaillement aux divers points du champ d’observation. Cette image d’élasticité IE k 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. 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, ou le module de cisaillement µ, ou le module d'Young E=3µ, ou l'atténuation α des ondes de cisaillement, ou l'élasticité du cisaillement µ1, ou la viscosité de cisaillement µ2, ou le temps de relaxation mécanique τs des tissus du milieu. 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éristiques de la viscosité locale du milieu. 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 : où ρ est la densité des tissus τS 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 : La résolution de cette équation de propagation avec l’ensemble des déplacements u permet d’obtenir les paramètres de propagation (cS, τS) cités ci-dessus. 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 moyennant par exemple les valeurs sur une bande de fréquences. On peut aussi utiliser les déformations partielles à la place des déplacements. On peut également établir des cartographies des paramètres de propagation, i.e. 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. 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 ce processus durant lequel on prend une pluralité d’images du milieu (images intermédiaires) destinées à déterminer une image d’élasticité IE (paramètre de propagation de l’onde de cisaillement). Ainsi, ces images du milieu ne sont pas corrigées par les mouvements du milieu. Par contre, l’image d’élasticité du milieu IE k 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 I 0 , ou pour pouvoir être comparée à l’image de déformation corrigée Iε k *. Détermination d’une image d’un paramètre de non linéarité INL k Dans un domaine élastique linéaire, le milieu 2 se déforme proportionnellement à la contrainte σ et suit la loi de Hooke : où E 0 est le module d’Young du domaine linéaire, et ε est la déformation 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. 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 : où E est le module d’Young, Δσ est la variation de contrainte locale, c’est-à-dire la contrainte partielle, Δε est la variation de déformation locale, c’est-à-dire la déformation partielle. 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 : « 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 « Quantitative Imaging of Nonlinear Shear Modulus by Combining Static Elastography and Shear Wave Elastography », H. Latorre-Ossa et Al., IEEE Transactions on Ultrasonics, Ferroelectrics, and Frequency Control, Vol. 51, no. 4, p. 833-839, 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ù E est le module d’Young ou module d’élasticité, E 0 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 σ est la contrainte et par exemple la contrainte σ z uni-axiale dans la direction Z, verticale, correspondant sensiblement à la compression induite par la pression exercée externe par l’utilisateur de la sonde 6, et A est le paramètre de non linéarité recherché. Par différentiation de la relation R1 précédente, on obtient : En divisant par la définition du module d’Young de l’équation (eq. 1) et par intégration on obtient une relation R2 entre l’élasticité (module d’Young) E et la déformation ε, soit : où E est le module d’Young E 0 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 ln() est la fonction de logarithme népérien ε est la déformation, et A est le paramètre de non linéarité à déterminer. Procédé ultrasonore pour quantifier l’élasticité non linéaire par ondes de cisaillement d’un milieu, le procédé comprenant les étapes suivantes : A1. la collecte d’une succession temporelle de données d’élasticité par ondes de cisaillement du milieu, A2. une application d’une déformation successivement changeante au milieu selon une séquence de déformations prédéterminée pendant la collecte des ondes de cisaillement, A3. l’observation de l’évolution de déformation réelle, et B. la quantification de l’élasticité non linéaire du milieu en fonction de la succession temporelle de données et l’évolution de déformation. Procédé selon la revendication 1, dans lequel l’observation de l’évolution de déformation réelle est effectuée en temps réel ou en quasi temps réel. Procédé selon l’une des revendications 1 et 2, dans lequel la déformation successivement changeante est appliquée par l’intermédiaire d’une sonde ultrasonore en successivement compressant ou en successivement décompressant le milieu par la sonde, la sonde étant simultanément utilisée pour collecter la succession temporelle de données et/ou pour observer l’évolution de déformation réelle. Procédé selon l’une quelconque des revendications 1 à 3, dans lequel la compression successive et/ou la décompression successive est effectuée de manière continue et progressive. Procédé selon l’une quelconque des revendications 1 à 4, dans lequel l’étape d’observation de l’évolution de déformation réelle comprend la collecte d’une succession temporelle de données ultrasons du milieu. Procédé selon l’une quelconque des revendications 1 à 5, dans lequel un niveau de déformation réelle estimée et un niveau de déformation cible sont notifiés sur une interface utilisateur, et dans lequel la déformation réelle estimée correspond à la déformation observée, notamment en temps-réel, et la déformation cible change selon la séquence de déformation prédéterminée, notamment en temps-réel. Procédé selon la revendication 6, dans lequel le niveau de déformation réelle estimé et le niveau de déformation cible sont affichés sur un écran de l’interface utilisateur. Procédé selon l’une quelconque des revendications 1 à 7, dans lequel les données d’élasticité par ondes de cisaillement produite à l’étape (A1) sont générées par les sous-étapes suivantes : A1.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, A1.2. une étape d'observation au cours de laquelle on observe la propagation de l'onde de cisaillement en acquérant une succession temporelle de données ultrasons du milieu, A1.3. une étape de traitement au cours de laquelle on détermine les données d’élasticité à partir desdites données ultrasons du milieu et d’un modèle de propagation d’onde de cisaillement. Procédé selon l’une quelconque des revendications 1 à 8, dans lequel la collecte d’une succession temporelle de données et l’observation du changement successive de la déformation sont mises en place par une séquence ultrasons en intercalant des tirs l'élasticité par ondes de cisaillement et des tirs d’élasticité statique, et/ou lesdites données d’élasticité statiques et lesdites de données d’élasticité par ondes de cisaillement sont temporellement entrelacées. Procédé selon l’une quelconque des revendications 1 à 9, comprenant en outre l’étape : C. de détermination d’une image et/ou un score et/ou un symbole montrant le niveau d’élasticité non linéaire du milieu. Procédé selon l’une quelconque des revendications 1 à 10, comprenant en outre les étapes : A4. la collecte d’une succession temporelle de données mode B du milieu pour calculer une compensation de mouvement du milieu relative à la sonde, A5. l’application de la compensation aux données d’élasticité. Procédé selon la revendication 11, dans lequel lesdites données mode B, lesdites données d’élasticité statique et lesdites données d’élasticité par ondes de cisaillement sont temporellement entrelacées. Procédé selon l’une quelconque des revendications 1 à 12, dans lequel les étapes (A1) à (A3) ou (A1) à (A5) sont effectuées simultanément, et/ou les étapes (A1) à (A3) ou (A1) à (A5) sont effectués en temps-réel ou en quasi temps-réel. Procédé selon l’une quelconque des revendications 1 à 13, dans lequel, pour chaque instant de la succession temporelle de données ou au moins pour quelques-uns des instants, la valeur d’un paramètre de non linéarité des ondes de cisaillement est déterminée en effectuant un ajustement de régression linéaire selon l’équation : dans laquelle : E(t) est l’élasticité à chaque instant de déformation lors de l’évolution de déformation, E 0 est l’élasticité à un premier instant (t0), ε est la déformation, et A est le paramètre de non linéarité de cisaillement. Procédé selon l’une quelconque des revendications 1 à 14, comprenant en outre : une étape (0) préalable à l’étape (A1), dans laquelle au moins une image ultrasons ou IRM (imagerie par résonance magnétique) ou mammographie ou image d’élasticité par ondes de cisaillement du milieu est déterminée, sans application d’une déformation selon une séquence de déformation prédéterminée, et déterminer une zone d’intérêt (ROI region of interest) pour laquelle on quantifie l’élasticité non linéaire selon les étapes (A1) à (B). Procédé selon l’une quelconque des revendications 1 à 15, dans lequel une carte de confiance est affichée en parallèle de l’affichage de l’image d’élasticité non linéaire. Procédé selon l’une quelconque des revendications 1 à 16, dans lequel pour chaque donnée ou pour chaque pixel de l’image, une valeur d’élasticité non linéaire en regard d’une valeur de référence prédéfinie est affichée. Procédé ultrasonore pour quantifier l’élasticité non linéaire par ondes de cisaillement d’un milieu en utilisant un dispositif d’imagerie (1) comprenant une sonde ultrasonore (6), un micro-ordinateur (4), et optionnellement une interface utilisateur, le procédé comprenant les étapes suivantes : positionner la sonde à la surface du milieu, activer un premier mode d’imagerie de l’élasticité non linéaire par ondes de cisaillement (NL-SWE), pendant lequel, la sonde appliquant une déformation selon une séquence de déformation prédéterminée, en progressivement compressant ou en progressivement décompressant le milieu : le dispositif collecte une succession temporelle de données ultrasons imageant les ondes de cisaillement générées dans le milieu à l’aide de la sonde permettant de déduire l’élasticité du milieu, et le dispositif collecte une succession temporelle de données ultrasons permettant de déduire une évolution de déformation du milieu appliqué par la sonde, quantifier l’élasticité non linéaire du milieu, qui est calculée en fonction de la succession temporelle de données et l’évolution de déformation observée, optionnellement afficher sur l’interface utilisateur une image de l’élasticité non linéaire du milieu, optionnellement enregistrer la succession temporelle d’images, comprenant des calculs et mesurages prédéfinies des images, optionnellement afficher sur l’interface utilisateur la succession temporelle d’images, et optionnellement afficher un résultat d’interprétation de l’image non linéaire. Procédé selon la revendication 18, dans lequel la séquence de déformation prédéterminée comporte une sous étape d’arrêt optionnellement automatique avant l’étape de quantification. Procédé selon l’une des revendications 18 et 19, comprenant en outre, notamment avant l’activation du premier mode : positionner la sonde à la surface du milieu, activer un deuxième mode d’imagerie d’élasticité par ondes de cisaillement (SWE), déterminer au moins une image d’élasticité par ondes de cisaillement du milieu en appliquant des ondes de cisaillement dans le milieu, sans application d’une déformation selon une séquence de déformation prédéterminée, enregistrer l’image, optionnellement comprenant des calculs et mesurages prédéfinies de l’image, déterminer une zone d’intérêt dans l’image pour imager l’élasticité non linéaire de la zone. Dispositif d’imagerie (1) comprenant une sonde ultrasonore (6) et un micro-ordinateur (4) adaptés pour mettre en œuvre le procédé pour quantifier l’élasticité non linéaire d’un milieu selon l’une quelconque des revendications précédentes.