Brevet FR2778838A1 — publié le 19991126

"PROCEDE DE DETECTION DE VARIATIONS D'ELASTICITE ET APPAREIL ECHOGRAPHIQUE POUR ME1TTRE EN CEUVRE CE PROCEDE" Description L'invention concerne un procédé de détection de variations d'élasticité dans un tissu mou soumis à une compression externe de direction axiale déterminée. L'invention concerne également un appareil échographique muni de moyens pour mettre en oeuvre ce procédé. L'invention trouve son application dans l'industrie de l'imagerie médicale Un procédé de reconstruction d'élasticité de tissu utilisant une méthode de perturbation linéaire est déjà connu de la publication intitulée " Tissue Elasticity Reconstruction Using Linear Perturbation Method " par F.Kallel et M. Bertrand, dans IEEE TRANSACTIONS ON MEDICAL IMAGING, VOL. 15, NO. 3, JUNE 1996, pp. 299-313. Cette publication décrit une méthode pour reconstruire le module élastique d'un tissu mou soumis à une compression externe statique à partir des mesures de déplacements. engendrés par cette compression. Ce procédé met en oeuvre un algorithme connu de résolution d'un problème inverse, appelé algorithme de NewtonRaphson, qui utilise une relation directe donnant l'image d'un ensemble de champs de déplacements par un Modèle à Eléments Finis des équations d'élasticité et qui adapte cette relation directe, au sens des moindres carrés, pour fournir la répartition des modules d'élasticité correspondants. L'ensemble de champs de déplacements axiaux de tissus forme les données de départ qui sont préalablement estimées en utilisant une technique de corrélation appliquée à des signaux ultrasonores. Les problèmes relatifs à la matrice qui permet de résoudre le problème inverse selon l'algorithme de Newton-Raphson sont pris en compte en utilisant une technique connue de régularisation de Tikhonov. Cette technique de régularisation fournit un compromis entre la fidélité aux données observées et l'information A PRIORI de la solution. En utilisant un modèle de formation d'image échographique, cette publication enseigne que l'algorithme converge avec seulement environ 10 itérations. Des figures de la publication citée représentent des images de répartition de module élastique obtenues à titre d'exemple avec 15 itérations par reconstruction à partir de données bruitées en utilisant l'algorithme de Newton-Raphson régularisé à chaque itération par le terme de Tikhonov. Dans le domaine médical de diagnostic d'anormalités dans les tissus mous, les tumeurs malignes se distinguent des tissus sains et des tumeurs bénignes par une différence d'élasticité de ces tissus, due à leur différence de structure. Dans le domaine du diagnostic du cancer du sein, il est particulièrement intéressant de disposer d'une méthode de mesure précise et fiable pour détecter les tumeurs malignes, aussi bien les tumeurs situées à fleur de peau que les tumeurs profondes de quelques centimètres qui sont très difficiles détecter. 2 2778838 Un problème est que la distribution du module d'élasticité dans un tissu mou n'est pas directement mesurable. Seul est mesurable le champ des déplacements dus à une compression du tissu mou. Mais ces déplacements sont extrêmement faibles d'o il résulte que leurs mesures ne sont pas directement exploitables par un praticien. Par contre la répartition des modules d'élasticité forme un ensemble de données tout à fait intéressantes du fait que ces données renseignent bien sur la nature des tissus et ont un contraste important facilement exploitable qui est lié aux mesures de déplacements dues à une compression du tissu mou. Un autre problème est que la répartition du module élastique n'est pas liée aux données du champ de déplacement par une relation directe. En revanche les données du champ de déplacement sont liées à la répartition du module élastique par une relation directe d'o il résulte que la répartition du module élastique dans un tissu doit être calculée par une méthode de résolution du problème inverse à partir des données du champ de déplacement. Une méthode de résolution d'un problème inverse est déjà décrite dans la publication citée. Un problème est que cette méthode connue est trop approximative pour s'appliquer à la détection de tumeurs malignes par la détermination de la répartition du module élastique dans un tissu, du fait que les données mesurées du champ de déplacement sont très faibles et que la participation du bruit lors de l'acquisition de ces données est extrêmement forte. C'est pourquoi la présente invention propose un procédé de détection de variation d'élasticité dans un tissu mou soumis à une compression externe de direction axiale déterminée, comprenant des étapes de: estimation d'un champ de déplacements axiaux dans le tissu, détermination d'un estimateur de module élastique incluant une opération de minimisation d'une distance entre une image d'une répartition de modules élastiques élémentaires par un Modèle à Eléments Finis et le champ de déplacements axiaux et régularisation de la solution de l'estimateur par une matrice diagonale (R) dont les coefficients afi sont des fonctions des déplacements axiaux (d) et sont respectivement appliqués aux modules élastiques élémentaires (e) pour rendre leur variance uniforme. Ce procédé fournit des images de reconstruction des variations d'élasticité dans le tissu qui sont beaucoup moins bruitées, beaucoup plus contrastées et qui permettent donc de détecter des très petites variations d'élasticité correspondant à des inhomogénéités du tissu et de très bien localiser ces défauts. Il en résulte une possibilité meilleure d'appliquer ce procédé à la détection des très petites tumeurs. Ce procédé est mis en oeuvre avantageusement par un appareil échographique qui constitue un moyen non invasif d'aide au diagnostic. L'invention est décrite ci-après en détail en référence aux figures schématiques annexées parmi lesquelles: 3 2778838 la FIG.1 représente par un diagramme de blocs les systèmes mettant en oeuvre les étapes du procédé de détection d'inhomogénéités dans un tissu; la FIG.2 représente par un diagramme de blocs les systèmes mettant en oeuvre les étapes fonctionnelles pour reconstruire la répartition du module d'élasticité à partir du champ de déplacement; la FIG.3 représente des lignes isodéplacement d'un champ de déplacements calculés à partir de signaux échographiques obtenus lors de la compression d'un tissu insonifié; la FIG.4 représente une image des variations du module élastique dans un tissu montrant deux inhomogénéités; les FIGs. 5A, 5B représentent une image des variations du module d'élasticité obtenue par la méthode linéaire de Newton-Raphson régularisée selon une méthode de régularisation de l'invention dans le cas de rapports signal/bruit respectivement égaux à et à 20, avec une seule itération et correspondant à la FIG.4; les FIGs.6A, 6B représentent une image des variations du module d'élasticité obtenue par la méthode de Newton-Raphson régularisée selon la méthode de Tikhonov dans le cas de rapports signal/bruit respectivement égaux à 50 et 20 et correspondant à la FIG.4, pour comparer aux FIGs.5A, 5B. L'invention concerne un procédé de traitement d'un signal échographique pour détecter des inhomogénéités dans un milieu insonifié. L'invention est fondée sur le traitement du signal échographique pour déterminer la répartition du module élastique dans le milieu insonifié en utilisant une méthode de résolution inverse régularisée, à partir de données d'un champ de déplacements mesurées dans ce milieu. L'invention trouve une application très intéressante dans la détection des tumeurs malignes par la détermination de la répartition du module élastique dans des tissus, ainsi que dans la détermination de l'étendue de ces tumeurs. A titre d'exemple, en référence à la FIG.1, le milieu est un tissu mammaire référencé texture 11. Une pression P est appliquée sur le tissu mammaire et une détermination d'un champ des déplacements est effectuée en soumettant une région du tissu à des tirs de signaux ultrasonores et en utilisant une méthode de corrélation appliquée aux signaux ultrasonores. D'une manière particulière, la pression P est appliquée par le moyen d'une compresseur 12 actionné par un moteur 13. La zone 11 du tissu mammaire soumise à la compression présente une face fixe F1, sur laquelle est appliquée une sonde 10 d'un échographe 100, et une face parallèle F2 sur laquelle est appliquée perpendiculairement la pression P du compresseur. Dans le dispositif représenté sur la FIG.1, selon l'invention, la paroi mobile du compresseur 12 appliquée sur la face F2 se déplace à une vitesse donnée v en direction de la face Fl. La sonde 10 4 2778838 appliquée sur la face fixe F1 envoie des signaux RF dans la texture Il suivant des tirs parallèles en direction de la face F2. Les lignes de tirs selon l'axe Z sont perpendiculaires aux faces F1 et F2, et parallèles à la direction de pression du compresseur 12. Les tirs sont générés par l'échographe 100 qui comprend un système 110 de focalisation des signaux RF dans le tissu 11 parallèlement à un axe Z, et un balayage par les tirs de signaux RF linéairement parallèlement à un axe X parallèle à la face Fl. Le balayage selon l'axe X est fait selon un nombre donné N de lignes de tirs avec une période de récurrence T. Lorsque les N tirs de 1 à N parallèlement à Z ont été effectués, la sonde recommence N tirs de 1 à N avec une période NT pour une ligne de tir donnée. La sonde reçoit en retour les signaux échographiques provenant de la texture 11 qui sont transmis à un système de détection de variation d'élasticité 200. Ce système 200 comprend un système 210 de corrélation temporelle pour effectuer la corrélation de chacun des signaux transmis par la sonde provenant d'une localisation déterminée dans la texture 11. Le système de corrélation 210 effectue la corrélation des signaux relatifs à deux tirs successifs de la sonde effectués en une localisation donnée et fournit l'amplitude des déplacements de chaque point de la structure 11 insonifié par la sonde sur chaque ligne de tir. Cette méthode fournit, en continu et en temps réel, durant l'action du compresseur 12, une image du champ de déplacement des structures formant la texture 11. A cet effet l'image du champ de déplacements dans la zone de la texture 11 insonifiée est numérisée. On ne stocke pas les signaux des différents balayages par la sonde. Ces signaux sont directement corrélés par le système de corrélation temporelle 210 qui fournit une image des déplacements de chaque point de la structure par exemple sous la forme de lignes isodéplacements, comme illustré par la FIG.3. Il faut bien noter que les déplacements calculés par le système de corrélation 210 résultent des mesures directes des déplacements des structures de la région insonifiées au moyen de l'échographe. Ces déplacements sont radiaux, parallèles à la pression P et obtenus par l'action de cette pression P externe au tissu 11. De manière préférentielle, plusieurs séries de balayages sont effectuées, avec entre chacune d'elles une compression égale ou un déplacement égal à v x NT de la face F2 du compresseur suivies d'un moyennage des données fournies par le système de corrélation 210. Cette méthode permet d'obtenir une grande précision dans cette étape d'acquisition 210 du champ de déplacements. Les données du champ des déplacements sont calculées dans cette étape 210 sous forme de déplacements élémentaires dl,... dû....d, relatifs à chaque pixel de l'image d'ultrasons numérisée obtenue par la méthode de corrélation. Les données de déplacements élémentaires constituent un vecteur d de données de déplacements. Ce vecteur déplacement d est lié à un vecteur appelé vecteur module élastique _e (ou module d'oung) par une relation directe F qui s'écrit selon (1) 2778838 d = F(e) (1) Le vecteur module élastique _ est constitué de modules élastiques élémentaires el,... el,....e, correspondant respectivement aux déplacements élémentaires relatifs à chaque pixel pi,. Pi,....pn de l'image numérisée. Il est connu que F(e) est l'image réalisée au moyen d'un Modèle à Eléments Finis de la répartition des modules élastiques élémentaires c'est-à-dire du vecteur e. Un Modèle à Eléments Finis est une technique de..DTD: calcul applicable aux équations d'élasticité et connu de l'homme du métier. Le problème qui doit être effectivement résolu est la détermination du vecteur e donnant la répartition du module élastique dans la région insonifiée ayant fourni le champ de déplacement. La détermination du vecteur e permet de construire une image de répartition du module élastique formée de modules élémentaires relatifs à chaque pixel de l'image numérisée. La valeur du vecteur _e est obtenu à partir de la relation (1) par la résolution du problème inverse selon la relation (2) ê (e) = arg min [lidF(e)j12] (2a) (2a) o ê (e) est un estimateur de _e. L'estimateur du module élastique e minimise la distance l Id - F(e)l 12 entre les données d acquises par des mesures dans l'étape 210 de la FIG.1 et l'image F(e) réalisée au moyen du Modèle à Eléments Finis. L'estimateur ê (e) permet de déterminer le vecteur e qui, une fois injecté dans le Modèle à Eléments Finis donne une valeur de vecteur d aussi proche que possible du vecteur mesuré d dans l'étape 1. La relation (2a) est connue de la publication citée au titre d'état de la technique et est référencée algorithme de Newton-Raphson. Cet algorithme permet de calculer le minimum de la distance li d - F(e)l 12 de manière itérative, par exemple, selon l'état de la technique avec 10 ou 15 itérations. L'algorithme de Newton-Raphson peut s'écrire sous la forme suivante: ê (e) = arg min [d - F(e)]T [d- F(e)] (2b) Cet algorithme peut être linéarisé de manière connue en faisant une approximation selon laquelle: d = S e + b (3) o S est une fonction de linéarisation sous la forme d'une matrice appelée matrice de sensitivité S dont une méthode de calcul est décrite dans la publication citée au titre d'état de la technique et qui est une fonction de d, et o b est un bruit blanc gaussien. Il en résulte une nouvelle écriture de l'algorithme de Newton-Raphson selon la relation suivante (4): ê = arg min [d - S e]T [d- S e] (4) Cet algorithme connu montre des inconvénients reconnus qui sont principalement des problèmes de convergence de l'algorithme et des problèmes d'instabilité de la solution liés au bruit qui est généralement important dans l'image de 6 2778838 répartition des déplacements obtenue par ultrasons. Un problème inverse dont la solution est instable à cause du bruit s'appelle un problème mal posé. Selon la publication citée au titre de l'état de la technique, ce problème est résolu par la méthode de régularisation de Tikhonov qui consiste à appliquer un terme de régularisation à l'algorithme de Newton- Raphson (4) selon la relation (5): ê(_e) = arg min [(d - S _e)T (d - Se) + A _eT I e_] (Sa) Le terme de régularisation de Tikhonov qui s'écrit À e IT e est une contrainte de appliquée à l'algorithme de Newton-Raphson pour diminuer l'effet du bruit. Ce terme aplanit la solution de l'estimateur en affectant à tous les modules élastiques élémentaires e1,... el,... en des coefficients de régularisation tous égaux. Dit autrement, dans le terme de régularisation, la matrice de régularisation est IT: la matrice (I) identité Transpose (T) affectée de coefficients tous égaux. Cela signifie que, dans l'équation (3) o l'estimateur ê (e) calcule un minimum, à la fois la distance.et le terme de régularisation doivent être aussi faibles que possible. Le terme de régularisation de Tikhonov peut être écrit selon (7a) de manière connue de l'homme du métier: i l A e IT _ = [ el2 + e2 +...ej2 +... + e2]=A ei (7a) o les e, sont les modules élémentaires e1 à en du vecteur e. L'algorithme de Newton- Raphson, auquel à chaque itération on applique la régularisation de Tikhonov sélectionne la solution o tous les el sont poussés avec la même force vers zéro. Cette méthode de régularisation suppose que l'objet dont on veut déterminer la répartition du module élastique est à l'origine homogène, ou encore que l'acquisition des données d a été faite de manière absolument homogène. On ne se trouve malheureusement jamais dans l'une de ces situations idéales. C'est pourquoi selon l'invention, on propose une nouvelle méthode de régularisation pour appliquer à l'algorithme de Newton-Raphson. Parmi les problèmes que la présente invention cherche à résoudre, se trouvent le fait que la mesure du champ des déplacements liés au vecteur module élastique par la relation directe (1) est affectée de bruit; et le fait que cette mesure dépend de la zone de la région insonifiée o ils ont été sont acquis: ainsi on ne peut pas traiter de façon indifférenciée les données provenant de zones insonifiées qui présentent une forte variabilité de données vis-à-vis des données de régions adjacentes, et les données provenant de zones o la variabilité des données est faible. La méthode de régularisation de l'invention repose donc sur le constat qu'il existe dans la région insonifiée des zones qui contiennent des informations à un niveau plus élevé que d'autres zones. En conséquence la méthode de régularisation selon l'invention est basée sur une détermination de coefficients ou poids à affecter à chaque 7 2778838 élément du vecteur module élastique pour tenir compte de la variabilité des mesures de déplacement et de la confiance que l'on accorde à ces informations dû à leur niveau résultant de ces mesures et dû à leur localisation dans la région insonifiée. Dans la formule (Sa) selon l'état de la technique, l'algorithme de newton- Raphson régularisé par le terme de Tikhonov effectue à chaque itération la minimisation d'une fonction de coût C qui s'écrit: C = [(d - S e)T (d Se) + À _T I e] (6) et qui permet d'écrire l'algorithme de Newton- Raphson régularisé par le terme de Tikhonov sous une nouvelle forme: ê (e) = arg min [STS + AI]' STd (5b) Le terme de régularisation de Tiknonov qui était xI (8a) dans la formule êi (e) (5b) de l'algorithme connu de newton-Raphson est transformé en une nouvelle matrice de régularisation notée. R (8b) Selon l'invention, le terme de régularisation est formé en utilisant une matrice R différente de la matrice identité I de Tikhonov multipliée par la constante À et en calculant explicitement les coefficients notés ai de cette nouvelle matrice de régularisation (R), ceci est différent de la régularisation de Tikhonov o les coefficients de la matrice de régularisation I montrés par la relation (6) sont tous égaux. Ainsi, selon l'invention, les coefficients a, de la matrice de régularisation R sont potentiellement tous différents pour être adaptés aux données di du champ de déplacement mesurées en chaque pixel de la région insonifiée, en sorte que chaque point ei de l'image de reconstruction de la répartition du module élastique est formé d'un élément a, ei du vecteur module élastique qui est spécifiquement adapté à une donnée di correspondante du champ de déplacement. Il en résulte que l'algorithme de Newton-Raphson régularisé selon l'invention s'écrit selon la formule (9a) suivante: êR (e) = arg min [(d S x e)T (d - Se) + _T R e] (9a) qui peut s'écrire selon (9b): êR(e) = arg min [STS + R]-' STd (9b) Selon l'invention, le terme de régularisation R peut ainsi être écrit selon la relation (7b) en utilisant les moyens d'écriture de la relation (7a) d'o il résulte que i-, an eT Re = e + ae, +...e...+,= (7b) 8 2778838 Cette relation (7b) exprime que les éléments diagonaux de la matrice R constituent les coefficients a, appliqués respectivement sur chaque éléments el du vecteur module élastique e. En référence à la FIG.2, I'invention propose un système 200 pour mettre en oeuvre la présente invention avec une méthode avantageuse de détermination de ces coefficients cl pour former des coefficients ou poids à appliquer aux éléments e1. Cette méthode est mise en oeuvre par exemple par un microprocesseur 200 qui met déjà en oeuvre l'opération 210 de corrélation temporelle et détermination du vecteur déplacement d. Cette méthode comporte d'abord: une opération 210 d'acquisition du champ de déplacement ou vecteur déplacement d une opération 221 de calcul de la matrice de linéarisation ou sensitivité S comme connu par exemple de la publication citée; une opération de décomposition tronquée en valeurs singulières de la matrice de sensitivité S. On rappelle que l'opération de décomposition d'une matrice en valeurs singulières consiste à écrire la matrice en question dans une base de vecteurs propres comme connu de l'homme du métier. Une telle décomposition définit des nouvelles matrices dont une matrice de changement de repère notée V et une matrice diagonale notée A formée d'éléments notés 8j. Selon l'invention, le niveau de troncature pour opérer la décomposition tronquée a été déterminé empiriquement et trouvé avantageusement de 2 à 10 inclus. La matrice diagonale A est utilisée pour construire une nouvelle matrice diagonale K dont les éléments diagonaux sont calculés par la formule (8) suivante: Ki = (Ai/cte))- AX (lia) une opération 222 de calcul de la matrice de régularisation R selon l'invention qui est une matrice diagonale dont les éléments diagonaux cç,, aussi notés c, sont les éléments diagonaux d'une matrice J non nécessairement diagonale obtenue par la formulation (9) suivante: J = V KVT (llb) Selon l'invention, le terme de régularisation adoucit la répartition des éléments de module élastique au moyen d'un vecteur formé des coefficients diagonaux de la matrice R qui force certains pixels de l'image à avoir des éléments e1 qui tendent plus fortement vers zéro que d'autres. La fonction de régularisation selon l'invention délimite un intervalle uniforme et force chaque valeur e1 à rester dans cet intervalle centré sur une valeur moyenne eMi spécifique de cet el. Ainsi, la définition de la matrice de régularisation R selon l'invention permet d'obtenir que la variance de e_, qui est 9 2778838 constituée d'un vecteur des écarts à la moyenne des e,, reste uniforme. Ce procédé de régularisation permet d'obtenir une image reconstruite de la répartition du module e de l'objet au moyen d'un estimateur de e qui réalise que les écarts à la moyenne de chacune des composantes e, à e, de e soient substantiellement égaux. Selon l'invention le vecteur de variance de e est fixé uniforme et la matrice de régularisation R est dérivée de la matrice de sensitivité S par les calculs matriciels simples décrits plus hauts. Ces calculs simples sont effectués dans un processeur associé à l'échographe comme montré sur la FIG.1. A l'issue de l'opération 222 de calcul de la matrice R, le processeur associé à I'échographe effectue: une opération 223 de calcul de la matrice M de l'estimateur de e à partir de la relation (9b) selon la relation: M = [STS + R]- ST (10) puis une opération 224 de multiplication matricielle de la matrice M par le vecteur d. Selon l'invention une seule itération suffit pour obtenir une image de reconstruction du module élastique e à partir du vecteur déplacement d par le procédé représenté par des blocs fonctionnels sur la FIG.1 et mis en oeuvre par le microprocesseur 200 associés à l'échographe 10, 100 de la FIG.1. Le terme de régularisation est donc appliqué à cette seule itération. D'o la simplicité de la mise en oeuvre de l'invention, comme illustré par la FIG.2 associée à la FIG.1. Les FIGs.4, SA, 5B et 6A, 6B illustrent les résultats obtenus par la méthode de l'invention. Ces FIGs peuvent être affichées au moyen d'un système d'affichage 6 comportant un moniteur, avec éventuellement des moyens d'enregistrement, associés à l'échographe 100. Ce système d'affichage est par exemple le même que celui qui permet l'affichage et l'enregistrement des images échographiques par l'intermédiaire du système de formation d'images 120. La FIG.4 montre, dans un tissu une zone insonifiée numérisée de largeur selon X de 5 cm et d'épaisseur selon Z de 5 cm. Ce tissu comprend deux inclusions ou inhomogénéités qui résultent en des variations du module élastique notées &e en kilopascals et qui apparaissent sou forme de deux reliefs réguliers. Les FIGs.5A et 5B sont des reconstructions avec la matrice de régularisation R de l'invention à partir de données de déplacement déterminées selon la méthode échographique illustrée par les FIGs. 1 et 2 de variations de du module élastique, avec les mêmes unités, dans le cas o en outre les données de déplacement contiennent du bruit. La FIG.5A a été établie pour un rapport signal/bruit égal à 50 et la FIG.5B pour un rapport signal/bruit égal à 20. 2778838 Les FIGs 6A et 6B sont des reconstructions avec le terme de régularisation de Tikhonov, réalisées dans les mêmes autres conditions que les FIGs.5A et 5B. En comparant, dans les FIGs SA et 6A (rapport signal/bruit égal à 50) les zones de tissu homogène entre les inhomogénéités, on peut constater que la méthode selon l'invention permet de complètement débruiter ces zones, alors que ces zones sont très bruitées par la méthode de régularisation de Tikhonov. Le résultat est que les images reconstruites selon l'invention sont mieux contrastées et que les variations d'élasticité dues à des inclusions ou à des inhomogénéités du tissu sont mieux localisées et plus facilement détectées. De même, dans les FIGs 5B et 6B (rapport signal/bruit égal à 20), en comparant les zones de tissus homogène entre les inhomogénéités, on peut constater que la méthode de l'invention offre encore un excellent contraste entre les différentes zones, alors que la méthode de régularisation de Tikhonov offre très peu de contraste, ce qui est un désavantage dans l'application à la détection de tumeurs. Ainsi un appareil échographique muni des moyens de mise en ceuvre du procédé illustré par les FIGs 1 et 2, constitue un excellent moyen non invasif de détection de tumeurs, et particulièrement de tumeurs mammaires. il 2778838 REVENDICATIONS: 1. Procédé de détection de variations d'élasticité dans un tissu mou soumis à une compression externe de direction axiale déterminée, comprenant des étapes de estimation d'un champ de déplacements axiaux dans le tissu, détermination d'un estimateur de module élastique incluant une opération de minimisation d'une distance entre une image d'une répartition de modules élastiques élémentaires par un Modèle à Eléments Finis et le champ de déplacements axiaux et régularisation de la solution de l'estimateur par une matrice diagonale (R) dont les coefficients ai sont des fonctions des déplacements axiaux (d) et sont respectivement appliqués aux modules élastiques élémentaires (e,) pour rendre leur variance uniforme. 2.. Procédé selon la revendication 1, dans lequel l'estimateur de module élastique est calculé en utilisant une matrice (M) formulée par: M = [ST S + R]-' ST o R est la matrice diagonale de régularisation et o S est une matrice de linéarisation du Modèle à Eléments Finis qui est une fonction des déplacements axiaux, et en effectuant la multiplication matricielle de ladite matrice (M) par le vecteur de champ de déplacements (_). 3. Procédé selon la revendication 2, o le tissu mou est couplé à une sonde ultrasonore associée à un échographe, cette sonde émettant et recevant des signaux échographiques parallèlement à la direction axiale de compression du tissu, et o l'estimation du vecteur champ des déplacements axiaux dans le tissu est effectuée en utilisant une technique de corrélation temporelle des signaux échographiques. 4. Appareil échographique pour détecter des variations d'élasticité dans un tissu mou incluant une sonde ultrasonore et des moyens de focalisation et balayage de la sonde qui est couplée à une surface de référence du tissu, et qui émet et reçoit des signaux échographiques parallèlement à une direction axiale du tissu, cet appareil comprenant: des moyens pour appliquer une compression externe au tissu parallèlement à la direction radiale, des moyens pour estimer un champ de déplacements axiaux dans le tissu, un système estimateur de module élastique pour effectuer une opération de minimisation d'une distance entre une image de répartition de modules élastiques élémentaires par un Modèle à Eléments Finis et le champ de déplacement axiaux, et un système de régularisation de la solution de l'estimateur par une matrice diagonale (R) dont les coefficients (a*,) sont des fonctions des déplacements 12 2778838 axiaux (di) et sont respectivement appliquées aux modules élastiques élémentaires (e) pour rendre leur variance uniforme. 5. Appareil selon la revendication 4, comprenant un processeur pour calculer l'estimateur de module élastique en utilisant une matrice (M) formulée par: M = [Ss + R]'1 sT o R est la matrice diagonale de régularisation et o S est une matrice de linéarisation du Modèle à Eléments Finis qui est une fonction des déplacements axiaux, et pour effectuer une multiplication matricielle de ladite matrice (M) par le vecteur de champ de déplacements (d). 6. Appareil selon la revendication 5 comprenant un processeur de traitement des signaux échographiques émis et reçus par la sonde pour effecteur une corrélation temporelle de ces signaux et pour fournir un champ des déplacements axiaux dans le tissu mou soumis à la compression externe axiale. 7. Appareil selon la revendication 6 dans lequel la sonde est une barrette de détecteurs linéaires émettant des tirs récurrents. 8. Appareil selon l'une des revendications 4 à 7 comprenant un système de formation d'images échographiques pour construire des images médicales à partir des signaux échographiques, un système d'affichage d'image pour afficher des images médicales du tissu et pour afficher des images de mesures de variations des déplacements du tissu soumis à la compression et des images de reconstruction de la répartition du module élastique dans le tissu pour visualiser les variations d'élasticité dans ce tissu. 9. Appareil selon la revendication 8, pour l'aide au diagnostic de tumeurs dans les tissus, comprenant des moyens pour localiser les zones des images de reconstruction du module élastique présentant un substantiellement fort contraste des variations d'élasticité, pour détecter des tumeurs potentielles associées à ces zones. 10. Appareil selon la revendication 9, pour l'aide au diagnostic du cancer du sein, o le compresseur est muni de moyens pour compresser le sein en appliquant une pression radiale sur une face du tissu, avec une autre face du tissu parallèle servant de référence aux déplacements radiaux.