FR2971695A1 - Mesure de viscoelasticite en utilisant une onde a ultrason modulee par amplitude de phase - Google Patents

Abstract

Une propriété de viscoélasticité du tissu est mesurée (32) in vivo. Pour récupérer plus d'information et/ou estimer la viscosité, le module de cisaillement, et/ou d'autres caractéristiques de cisaillement, une forme d'onde modulée en amplitude et en phase est transmise (30) au tissu. Le déplacement créé par la forme d'onde au cours du temps comporte des déplacements associés à la réponse à différentes fréquences. En examinant le déplacement dans le domaine de fréquence, une ou plusieurs propriétés de viscoélasticité peut être calculée (40) pour différentes fréquences. La réponse en fréquence peut indiquer le caractère fin du tissu.

Description

MESURE DE VISCOELASTICITE EN UTILISANT UNE ONDE A ULTRASON MODULÉE PAR AMPLITUDE DE PHASE Les modes de réalisation de la présente invention se rapportent à des mesures de viscoélasticité. En particulier, une propriété de viscoélasticité de tissus est mesurée in vivo. Des caractéristiques de module de cisaillement d'un échantillon de tissu peuvent être mesurées à l'aide d'un rhéomètre. Un rhéomètre est un dispositif mécanique capable de soumettre un échantillon à une déformation qui est soit dynamique (par exemple, sinusoïdale) soit statique (par exemple, linéaire). Le rhéomètre peut mesurer le module de cisaillement de stockage et le module de cisaillement de perte pour un domaine de fréquences. Le domaine est limité sur la base de la rigidité de matériau, tel que de 1 à 10 Hz pour des tissus mous. Pour un usage médical, le tissu est extrait d'un patient pour le placer dans le rhéomètre. Le rhéomètre n'est pas utilisé pour des mesures in vivo. Les mesures dépendent de la dimension et de la forme de l'échantillon de tissu extrait, ainsi que des conditions aux limites dues à l'extraction. Des caractéristiques de cisaillement peuvent être mesurées in vivo par ultrason. Par exemple, une détection de vitesse de cisaillement est utilisée dans diverses applications de diagnostic, comme pour évaluer une maladie du foie. Pour une détection d'onde de cisaillement, une impulsion d'ultrason de poussée (par exemple, une impulsion à un cycle) est déclenchée le long d'une ligne de balayage. L'impulsion de poussée produit une onde de cisaillement qui entraîne le déplacement du tissu. Le déplacement est mesuré. Pour détecter la vitesse d'onde de cisaillement, plusieurs impulsions de poussée le long d'une même ligne de balayage et des balayages de détection de déplacement correspondants sont utilisés. Ces mesures de cisaillement peuvent avoir des informations limitées ou être fonction du niveau de compression. A titre d'introduction, les modes de réalisation préférés décrits ci-dessous incluent des procédés, des instructions, et des systèmes pour mesurer une propriété de viscoélasticité de tissus mous in vivo. Pour récupérer plus d'informations et/ou estimer des caractéristiques de viscosité, de module de cisaillement, et/ou d'autres caractéristiques de cisaillement, une forme d'onde modulée par déplacement d'amplitude et de phase est émise vers le tissu. Le déplacement créé par la forme d'onde au cours du temps inclut des déplacements associés à une réponse à différentes fréquences. En examinant le déplacement dans le domaine de fréquences, une ou plusieurs propriétés de viscoélasticité peuvent être calculées pour différentes fréquences. La réponse en fréquence peut indiquer la santé du tissu. Suivant un premier aspect, un support de stockage pouvant être lu par ordinateur de manière non transitoire comporte mémorisées en son sein des données représentant des instructions pouvant être exécutées par un processeur programmé pour mesurer une propriété de viscoélasticité de tissu in vivo. Le support de stockage comporte des instructions pour émettre une forme d'onde modulée en phase et modulée en amplitude comportant des cycles à différentes fréquences à différents instants, les différentes fréquences se trouvant dans un domaine de fréquences, la forme d'onde modulée en phase et modulée en amplitude étant émise vers le tissu dans un patient, pour calculer le déplacement en fonction du temps du tissu en réponse à la forme d'onde modulée en phase et modulée en amplitude, le déplacement calculé à partir de balayages du tissu, pour appliquer une transformée de Fourier au déplacement au cours du temps ; et pour déterminer la propriété de viscoélasticité à partir d'une transformée de Fourier d'une équation d'onde de cisaillement et de la transformée de Fourier du déplacement au cours du temps. Suivant un deuxième aspect, il est visé un procédé pour la mesure de viscoélasticité en utilisant l'ultrason. Une forme d'onde à ultrason modulée en phase et en amplitude est émise dans un patient in vivo. Une propriété de viscoélasticité est mesurée à partir du cisaillement créé par la forme d'onde à ultrason modulée en amplitude et en phase. La propriété de viscoélasticité est mesurée sur un domaine de fréquences associé à la forme d'onde à ultrason modulée en amplitude et en phase et indépendamment d'une quantité de compression appliquée de manière externe au patient pendant l'émission. Suivant un troisième aspect, il est visé un système pour calculer une onde de cisaillement en utilisant l'ultrason. Un formateur de faisceau d'émission peut fonctionner pour produire une forme d'onde modulée en amplitude et en phase. Un transducteur à ultrason est relié au formateur de faisceau d'émission de sorte que le transducteur à ultrason émette de l'énergie acoustique vers du tissu dans un patient en réponse à la forme d'onde modulée en amplitude et en phase. Un formateur de faisceau de réception peut fonctionner pour émettre en sortie des données représentant des emplacements dans l'espace en fonction de signaux acoustiques reçus. Un processeur est configuré pour estimer le déplacement du tissu au cours du temps en fonction des données de sortie et pour calculer l'information de cisaillement en fonction du déplacement du tissu au cours du temps. Un affichage peut fonctionner pour afficher une image qui est fonction de l'information de cisaillement. La présente invention est définie par les revendications qui suivent et rien dans cette section ne doit être considéré comme une limitation de ces revendications. D'autres aspects et avantages de l'invention sont décrits ci-dessous en liaison avec les modes de réalisation préférés et pourront être plus tard revendiqués indépendamment ou en combinaison.

Les composants et les figures ne sont pas nécessairement à l'échelle, l'important étant d'illustrer le principe de l'invention. En outre, aux figures, des références numériques identiques désignent des parties correspondantes dans toutes les différentes vues.

La Figure 1 est un organigramme d'un mode de réalisation d'un procédé pour mesurer une propriété de viscoélasticité de tissu in vivo. la Figure 2 est une illustration sous forme de graphique d'une transmission à ultrason induisant du 20 cisaillement ; la Figure 3 illustre à titre d'exemple une forme d'onde modulée en amplitude et en phase ; la Figure 4 représente un mode de réalisation de balayages pour la détection de déplacement relative à la 25 forme d'onde modulée en amplitude et en phase de la Figure 3 ; la Figure 5 illustre des emplacements relatifs d'un transducteur pour créer une propagation d'onde de cisaillement à l'intérieur du tissu ciblé et des bobines 30 de radiofréquence pour mesurer le déplacement dans un mode de réalisation à résonance magnétique donné à titre d'exemple ; et la Figure 6 est un mode de réalisation d'un système pour un calcul d'onde de cisaillement en utilisant l'ultrason. Des propriétés mécaniques de tissu sont mesurées in vivo. Les propriétés de n'importe quel structure ou matériau chez un patient, tel qu'un tissu mou, peuvent être mesurées. Une forme d'onde à ultrason modulée en phase (PM) et modulée en amplitude (AM) entraîne un déplacement de tissu. La forme d'onde AM-PM produit une vibration à des fréquences successives à l'intérieur du patient au niveau d'un tissu souhaité. Le module de cisaillement de tissu et la viscosité du tissu sont mesurés à diverses fréquences avec cette excitation, permettant des mesures en moins de deux secondes. La modulation de phase à des fréquences inférieures (par exemple, de 1 à no Hz) peut fournir un rapport signal bruit plus élevé. La forme d'onde AM-PM permet une récupération d'information pendant que l'on induit le cisaillement ou après que l'on ait induit le cisaillement. Un procédé rapide et fiable peut évaluer la rigidité de matériau pour diverses fréquences en une seule mesure. En récupérant des mesures pour différentes fréquences, la connaissance qui en résulte des propriétés mécaniques du matériau (par exemple, le tissu) peut venir en complément des procédures de diagnostic. Cette connaissance peut comporter des paramètres quantitatifs, tels que la contrainte, la viscosité et/ou le module de cisaillement. Les paramètres quantitatifs sont indépendants de la grandeur de compression appliquée au tissu puisque les paramètres sont dérivés du déplacement. La réponse en fréquence peut être indépendante de la compression. La connaissance peut comporter des paramètres qualitatifs, tels que le déplacement. Le déplacement peut être fonction d'une grandeur de pression appliquée au tissu. La connaissance peut comporter des informations dérivées.

La pente de la courbe donnant le module de cisaillement en fonction de la fréquence et/ou la pente de la courbe donnant la viscosité en fonction de la fréquence peut fournir des informations se rapportant à des propriétés mécaniques du tissu liées à une pathologie du tissu. Les pentes sont indépendantes des niveaux de précompression. Les mesures peuvent être répétées pour augmenter la précision. Comme les mesures ne reposent pas sur un contact, les mesures sont indépendantes de condition aux limites, de la dimension des matériaux et de leur forme. La Figure 1 représente un procédé pour une mesure de viscoélasticité en utilisant l'ultrason. Une transmission à ultrason est utilisée pour produire un déplacement de tissu en réponse à différentes fréquences. De cette manière, l'ultrason est utilisé pour calculer une propriété de viscoélasticité. Le procédé est mis en oeuvre par le système de la Figure 6 ou par un système différent. Des actes ou étapes supplémentaires ou différents, ou moins d'actes ou d'étapes peuvent être prévus. Par exemple, le procédé est effectué sans les actes 38, 42, et/ou 44. A titre d'un autre exemple, l'acte 32 est effectué sans l'un ou plusieurs des actes 34 à 40. Les actes sont effectués dans l'ordre décrit ou représenté, mais peuvent être effectués dans d'autres ordres A l'acte 30, une forme d'onde modulée en phase et modulée en amplitude est émise. Un transducteur à ultrason émet une forme d'onde acoustique convertie à partir d'une forme d'onde électrique. L'énergie acoustique avec la forme d'onde modulée en phase et modulée en amplitude est émise vers le tissu chez un patient. L'émission a lieu in vivo. La forme d'onde acoustique est émise pour produire une onde de cisaillement. L'excitation est une impulsion à ultrason de poussée. L'énergie acoustique est focalisée, résultant en un ou plusieurs faisceaux pour chaque émission. L'excitation est focalisée en utilisant un foyer mécanique et/ou un réseau phasé. L'excitation est focalisée à un emplacement pour permettre la détection de l'onde de cisaillement résultante, comme celle focalisée en un emplacement de tissu entourant une tumeur possible et/ou incluant une tumeur possible. Comme représenté à la Figure 2, l'onde 52 de cisaillement est produite à la région 50 de foyer par le. transducteur 54 et se propage latéralement à partir de cette région 50 de foyer. Les flèches sont représentées dans une direction (par exemple, horizontalement), mais l'onde de cisaillement se propage dans des directions multiples. L'onde de cisaillement diminue d'amplitude au. fur et à mesure que l'onde se propage à travers le tissu. Pour produire l'onde de cisaillement, des excitations de puissance ou d'amplitude élevée sont souhaitées. Par exemple, l'excitation a un indice mécanique qui est proche de 1,9 sans cependant le dépasser. Pour être conservateur et prendre en compte une variation de la sonde, un indice mécanique de 1,7 ou un autre niveau peut être utilisé en tant que la limite supérieure. Des puissances plus grandes (par exemple, un MI (index mécanique) dépassant 1,9) ou moindres peuvent être utilisées. L'émission le long de lignes de balayage identiques ou voisines peut faire en sorte que la température du tissu augmente au cours du temps. Des effets biologiques peuvent inclure l'hyperthermie à des températures de tissu d'environ 41-45°C, une dénaturation de protéine à des températures supérieures à 43-45°C, et une nécrose de tissu à des températures supérieures à 50°C. La rigidité du tissu peut être affectée même à des températures inférieures à 43-45°C. A des températures supérieures à 43-45°C, des augmentations de la viscosité et/ou de 1a rigidité peuvent avoir lieu. A des températures supérieures à 50°C, le tissu peut avoir une rigidité élevée et/ou une atténuation élevée. Des effets biologiques sont limités en empêchant une augmentation de température de plus de 2 degrés Celsius. En variante, les émissions peuvent entraîner des effets biologiques. Les formes d'ondes électriques et acoustiques correspondantes comportent des cycles à différentes fréquences à différents instants. Les différentes fréquences se trouvent dans un domaine de fréquences. Le domaine de fréquences peut être rendu optimum sur la base du type de tissu qui est examiné. Par exemple, des tissus du sein peuvent avoir un déplacement maximum en réponse à 50 Hz, de sorte que la gamme comporte le 50 Hz. A titre d'un autre exemple, des tumeurs peuvent être plus dures que du tissu mou, de sorte qu'ils ont une fréquence plus grande à laquelle la réponse en déplacement maximum a lieu. Pour la mesure de tumeur, le domaine de fréquence peut être plus grand. Suivant un mode de réalisation, le domaine de fréquence se trouve entre 15 et 75 Hz pour la modulation de phase. N'importe quelle fonction de variation de fréquence peut être utilisée pour la modulation de phase. Par exemple, un balayage en fréquence chirp ou à impulsion dans le domaine est utilisé. La Figure 3 représente une forme d'onde avec une enveloppe modulée en amplitude et en phase. L'onde modulée en phase commence à la fréquence la plus basse et balaye ou se modifie graduellement jusqu'à la fréquence la plus élevée dans le domaine. D'autres variations de fréquence peuvent être utilisées, telles qu'une variation non-linéaire, aléatoire, ou avec d'autres étapes entres les différentes fréquences dans n'importe quel ordre.

La modulation d'amplitude se fait suivant n'importe quelle fonction. La modulation d'amplitude sépare des fréquences individuelles ou des groupes de fréquences individuelles les unes des autres.

Suivant un mode de réalisation, la forme d'onde modulée en phase et en amplitude est produite en multipliant une porteuse sinusoïdale (dans le domaine du MHz) par une modulation d'amplitude sinusoïdale avec un terme variant en phase (dans le domaine des Hz). La modulation d'amplitude avec le terme variant en fonction de la phase définit une enveloppe avec la modulation d'amplitude et de phase. La forme d'onde peut être représentée sous la forme : x(t) =Asin(coc t)X Sin«a),n +acorn t)t} où x(t) est la forme d'onde en fonction du temps, A est une pondération d'amplitude, w, est une fréquence centrale ou une fréquence de porteuse, et Awm est une fréquence de modulation qui varie au cours du temps. Suivant un mode de réalisation wm + Awm varie entre 15 et 75 Hz, mais d'autres domaines peuvent être utilisés. La Figure 3 représente une forme d'onde avec le domaine de 15 à 75 Hz. La forme d'onde comporte n'importe quel nombre de cycles d'enveloppe et de porteuse. Suivant un mode de réalisation, la forme d'onde a des cycles suffisants pour être de 1333 millisecondes, mais des longueurs de temps plus grandes ou moindres peuvent être utilisées. Le nombre de cycles de la porteuse est en dizaine, centain milliers, ou plus. Le nombre de cycles de l'enveloppe ou de la modulation est de deux ou plus. La réponse du tissu est fonction de la forme d'onde, x(t) et des caractéristiques du tissu. Le déplacement, y(t), du tissu au cours du temps peut être exprimé sous la forme d'une convolution de la forme d'onde, x(t) et des caractéristiques ou réponse de tissu, h(t): y(t)=x(t)*h(t). La réponse du tissu reflète les propriétés de viscoélasticité du tissu. Une ou plusieurs propriétés de viscoélasticité sont mesurées à partir du cisaillement créé par la forme d'onde à ultrason modulée en amplitude et en phase à l'acte 32. Afin de mesurer les propriétés de viscoélasticité à l'acte 32, le déplacement, y(t), du tissu au cours du temps en réponse à la forme d'onde modulée en amplitude et en phase est mesurée à l'acte 34. L'acte 32 est représentée comme incluant les actes 34, 36, 38 et 40. Des actes différents, supplémentaires, ou moins d'actes peuvent être prévus. Dans d'autres modes de réalisation, d'autres processus sont utilisés pour mesurer la propriété de viscoélasticité. A l'acte 34, le déplacement est calculé en fonction du temps. Le tissu est balayé plusieurs fois pour déterminer le déplacement, par exemple en balayant une région au moins trois fois pour déterminer le déplacement à deux instants différents. Le déplacement du tissu se fait en réponse à la forme d'onde modulée en phase et modulée en amplitude. Le déplacement du tissé créé par l'onde de cisaillement est déterminé au cours du temps. Au fur et à mesure que l'onde de cisaillement passe devant un emplacement donné, le tissu se déplace d'une quantité ou distance qui augmente jusqu'à une valeur pic et ensuite diminue au fur et à mesure que le tissu revient au repos. Comme la forme d'onde AM-PM comporte de nombreux cycles, le tissu peut être déplacé de manière continue. La modulation en amplitude entraîne une variation du déplacement au cours du temps au fur et à mesure que le tissu commence à revenir à un état normal au fur et à mesure que l'amplitude diminue. La modulation de phase résulte en une variation de la quantité du déplacement au cours du temps. Le tissu est balayé en utilisant n'importe quelle modalité d'imagerie capable de balayer un déplacement pendant la réponde du tissu à la forme d'onde de poussée, comme par exemple pendant ou après l'application de la forme d'onde AM-PM. Le balayage a lieu avant que le tissu ne revienne à un état normal ou relâché, mais peut comporter des balayages du tissu à des instants où le tissu est revenu au repos. Une imagerie par résonance magnétique ou à ultrason sont deux modalités possibles pour calculer le déplacement à partir de balayages du tissu. Pour un balayage à ultrason, l'onde 52 de cisaillement est détectée à des emplacements voisins de la région de foyer et/ou à distance de la région de foyer pour la forme d'onde AM-PM. Pour détecter le déplacement, une énergie à ultrason est émise vers les tissus subissant le déplacement, et des réflexions de l'énergie sont reçues. Pour détecter la réponse du tissu à des ondes de cisaillement dans un région d'intérêt, des émissions sont faites vers d'autres régions de foyer, et la détection est effectuée autour des autres régions de foyer. Ces autres transmissions servent à détecter les ondes de cisaillement plutôt que de créer l'onde de cisaillement. Les émissions pour la détection peuvent avoir une puissance plus faible et/ou des impulsions courtes (par exemple, de 1 à 5 cycles de porteuse) et utiliser la ligne de balayage différente ou la même ligne de balayage que la forme d'onde AM-PM. Les émissions pour la détection peuvent avoir un profil de faisceau plus large le long d'au moins une dimension, par exemple latéralement, pour former en simultané des échantillons de réception le long d'une pluralité de lignes de balayage. L'onde de cisaillement peut être surveillée dans une, deux, ou plus de deux directions. Une région d'intérêt est surveillée pour détecter l'onde de cisaillement. La région d'intérêt est de n'importe quelle dimension, par exemple 6 mm dans la direction latérale et 10 mm en direction axiale entourant l'emplacement de foyer de la forme d'onde AM-PM. Cette région de détection est surveillée par ultrason. Par exemple, des balayages en mode B sont effectués pour détecter un déplacement de tissu créé par l'onde de cisaillement. Des modes Doppler, d'écoulement en couleur, ou d'autres modes à ultrason peuvent être utilisés pour surveiller l'onde de cisaillement. La surveillance est effectuée pour n'importe quel nombre de lignes de balayage. Par exemple, quatre faisceaux de réception sont formés en réponse à chaque transmission. Après avoir émis l'excitation pour produire l'onde de cisaillement, des émissions en mode B sont effectuées de manière répétée le long d'une ou de plusieurs lignes de balayage d'émission et des réceptions le long de lignes de balayage de réception correspondantes. Dans d'autres modes de réalisation, uniquement un seul faisceau de réception ou un autre nombre de faisceaux de réception sont formés en réponse à chaque émission. Certaines des données d'ultrason, telles qu'au début ou à la fin des répétitions, peuvent ne pas être une réponse à l'onde de cisaillement. L'émission et la réception pour la détection sont effectuées une pluralité de fois pour déterminer des variations dues au déplacement au cours du temps. N'importe quelle séquence d'émission et de réception peut être utilisée. La détection de déplacement peut être entrelacée avec un autre balayage, tel que le balayage de différentes régions pour un déplacement de manière séparée. L'émission et la réception pour la détection sont entrelacées avec la forme d'onde modulée en amplitude et en phase. Par exemple, la ou les émissions et la ou les réceptions pour balayer une région une fois sont effectuées chaque fois que la modulation d'amplitude se trouve à un niveau zéro, ou à proximité d'un niveau zéro. A titre d'un autre exemple, le balayage du tissu est entrelacé avec l'émission de l'acte 30 sur une base périodique, telle que chaque 1,67 millisecondes. L'entrelacement évite l'interférence, telle que le cas où la forme d'onde de modulation en phase et en amplitude qui serait arrêtée pendant le balayage, de sorte que des échos provenant de la forme d'onde sont minimisés. La Figure 4 représente un exemple d'entrelacement de balayages des plans 56 avec la forme d'onde modulée en amplitude et en phase. La durée d'entrelacement peut être plus ou moins fréquente. Lorsque la modulation d'amplitude et de phase se poursuit après avoir cesser le balayage de déplacement, la forme d'onde part d'un point dans la forme d'onde où l'émission a cessé. Dans des modes de réalisation en variante, le balayage est effectué à une fréquence différente ou avec un codage différent de la forme d'onde modulée en amplitude et en phase. La forme d'onde de poussée et le balayage peuvent tous les deux avoir lieu simultanément et la fréquence ou le codage est utilisé pour distinguer des échos de chacun. La description ci-dessus se rapporte à une profondeur ou emplacement. La propriété de viscoélasticité est mesurée pour un emplacement. Pour surveiller une région plus grande, des actes 30 à 40 sont répétés pour d'autres emplacements. Pour chaque emplacement de faisceau de réception, un profil temporel d'information de déplacement est fourni, représenté par les données à ultrason. Un profil temporel distinct est fourni pour chaque profondeur axiale ainsi que chaque emplacement latéral.

Le déplacement est calculé à partir des données de balayage à ultrason. Le tissu se déplace entre deux balayages. Les données d'un balayage sont traduites en une, deux, ou trois dimensions relatives aux données dans l'autre balayage. Pour chaque position relative possible, une valeur de similarité est calculée. La valeur de similarité est déterminée avec une corrélation, telle qu'une corrélation croisée. Une somme minimum de différences absolues ou une autre fonction peuvent être utilisées. Le décalage dans l'espace avec la corrélation la plus élevée ou la corrélation suffisante indique la quantité et la direction de déplacement. Les déplacements sont déterminés pour un emplacement donné à différents instants, tels qu'associés avec des balayages séquentiels. Le déplacement est déterminé par rapport à une trame initiale ou de référence de données balayée (c'est-à-dire des déplacements cumulatifs). En variante, le déplacement est déterminé à partir de la trame immédiatement précédente des données de balayage, ce qui affecte la trame précédente en tant que référence sur une base qui se poursuit (c'est-à-dire un déplacement par incrément). Le profil temporel pour un emplacement donné indique le déplacement créé par l'onde de cisaillement au cours du temps et en réponse à différentes parties de la forme d'onde modulée en amplitude et en phase. Pour le balayage avec une résonance magnétique, un processus similaire est utilisé. La Figure 5 représente le transducteur 54 qui émet vers une région de foyer entourée au moins en partie par un alésage de résonance magnétique. L'alésage comporte des bobines de fréquence radio. Un champ magnétique uniforme est produit. En utilisant un amplificateur à gradient de fréquence radio d'impulsion, les spins des molécules chez le patient sont modifiés. Les bobines sont utilisées pour détecter la modification, faisant l'acquisition de données d'espace k. Plutôt qu'effectuer une corrélation entre des trames dans une séquence, un ensemble de données de déplacement au cours du temps est acquis. Des données de résonance magnétique d'espace k représentant le tissu du patient sont acquises à différents instants ou en séquence. Un ensemble (ensemble de référence) est acquis sans que le tissu soit soumis à la forme d'onde modulée en amplitude et en phase. Un autre ensemble est acquis pendant l'application de la forme d'onde modulée en amplitude et en phase. Comme l'ultrason de la forme d'onde n'interfère pas avec les données d'espace k, les données de balayage peuvent être acquises sans entrelacement. En variante, un entrelacement est utilisé. La séquence du balayage est périodique (par exemple, chaque 1,67 ms) sur la base des parties d'amplitude basse de la forme d'onde modulée en amplitude, ou une autre fonction. Les données d'espace k sont acquises en synchronisation avec un ou plusieurs cycles physiologiques, tel que le cycle cardiaque ou de respiration. Les balayages de la séquence de référence sont acquis aux mêmes points dans le cycle que l'ensemble de déplacement. Alors que les données d'espace k peuvent être utilisées sans être traitées en une image, des données d'image sont utilisées dans des modes de réalisation en variante. Les données de résonance magnétique de référence représentant le tissu libre de répondre à l'émission sont soustraites des données de résonance magnétique représentant la réponde du tissu à l'émission. En variante, les données représentant la réponse du tissu sont soustraites des données de référence. Les balayages provenant d'instants similaires relatifs au cycle physiologique sont soustraits. La soustraction isole des différences. Comme l'onde de cisaillement crée une différence, le tissu déplacé ou décalé reste et les autres signaux s'annulent à partir de la soustraction.

Après un filtrage passe bas ou un autre traitement, la quantité, la direction, ou la quantité et la direction du déplacement entre des trames séquentielles est déterminée. Un emplacement d'un point d'intensité pic ou une région d'intensité pic est identifié dans les données pour chaque trame de données de soustraction. La différence dans les emplacements entre les trames de séquence ou entre une trame différente et une trame de référence est calculée en tant que le déplacement. Le déplacement varie au cours du temps en raison des différences dans la forme d'onde modulée en amplitude et en phase qui est appliquée aux différents instants. A l'acte 36, le déplacement au cours du temps est transformé en le domaine de fréquence. Une transformée de Fourier est appliquée au déplacement au cours du temps.

N'importe quelle transformée peut être utilisée. Comme le déplacement au cours du temps est en réponse à différentes fréquences de la forme d'onde modulée en amplitude et en phase, la transformée dans le domaine de fréquence fournit différents niveaux de réponse en fonction de la fréquence. La transformée de Fourier du déplacement, y(t), peut être représentée sous la forme Y(Q), où Q est la fréquence, c'est-à-dire la représentation du déplacement dans le domaine de fréquence. La transformée du déplacement ou de la réponse du tissu est représentée sous la forme h(t)< > H(M - X~ YS)

A l'acte 38, le déplacement dans le domaine de fréquence est filtré. Tout filtrage peut être utilisé, tel qu'un filtrage à bande passante. Le domaine de fréquence du filtrage à bande passante est réglé sur la base du domaine de fréquences de la forme d'onde modulée en phase et en amplitude émise. La force acoustique est fonction du carré de la pression. Il en résulte que le déplacement ou la vibration induite par la forme d'onde est deux fois la fréquence de la forme d'onde. Par exemple, la forme d'onde comporte une variation de fréquence de l'enveloppe de 15 à 75 Hz, de sorte que l'information de déplacement en réponse a lieu entre 30 et 150 Hz. La largeur du domaine de fréquence est deux fois la largeur de la forme d'onde, et les fréquences basse et haute dans le domaine sont deux fois les fréquences basse et haute du domaine de la forme d'onde. Le filtrage isole l'information dans ce domaine ou dans un sous-domaine. Des données en fréquence à l'extérieur du domaine ne sont pas utilisées et des données à l'intérieur du domaine sont utilisées, ce qui se traduit par un filtrage à bande passante. D'autres filtrages à bande passante, d'autres filtrages, d'autres traitements ou aucun traitement de modification de données peuvent être utilisés. A l'acte 40, la propriété de viscoélasticité est déterminée à partir de la transformée de Fourier du déplacement au cours du temps. L'onde de cisaillement peut être représentée sous la forme : a2yz _ c2Av = Q at a2yz /' a at2 - ~/1 + 17ât, °yz D'autres représentations d'onde de cisaillement peuvent être utilisées. Dans le domaine de fréquences, l'équation d'onde de 5 cisaillement est représentée sous la forme : -w2Yp) AY(Q) ,~ X1 1 Cette équation représente la transformée de Fourier de l'équation d'onde de cisaillement. D'autres représentations peuvent être utilisées, en fonction de la 10 représentation d'onde de cisaillement et de la transformée de Fourier utilisée. Les termes du côté gauche de l'équation sont connus ou mesurés, tels que le mouvement (déplacement) résultant dans le domaine de fréquences, Y(0), qui oscille à une fréquence (w) égale à 15 deux fois la fréquence AM-PM, par exemple 2x(15 à 75 Hz)=30 Hz à 150 Hz. Le côté droit représente les parties imaginaire et réelle de l'équation d'onde de cisaillement transformée. La propriété de viscoélasticité est déterminée à 20 partir de l'équation d'onde de cisaillement dans le domaine de fréquence. N'importe quel module, viscosité, ou valeur de cisaillement peut être estimé. Des valeurs de module de tissu représentent la dureté ou la rigidité du tissu. Par exemple, le module de cisaillement du tissu 25 est estimé. Dans des modes de réalisation en variante, le module de Young est estimé. Dans d'autres modes de réalisation, d'autres valeurs de cisaillement sont estimées, qu'elles soient quantitatives ou qualitatives. Dans un mode de réalisation, la viscosité est 30 déterminée. La viscosité est calculée en fonction de résultats de l'application de la transformée de Fourier. La partie imaginaire de la transformée de Fourier de l'équation d'onde de cisaillement peut être utilisée pour déterminer la viscosité. La partie imaginaire est divisée par la fréquence, w, pour obtenir la viscosité. Ceci est représenté sous la forme : -w2Y(Q) dY(~2) CO Suivant un mode de réalisation en variante ou supplémentaire, le module de cisaillement est déterminé.

Le module de cisaillement est calculé en fonction des résultats de l'application de la transformée de Fourier. La partie réelle de la transformée de Fourier de l'équation d'onde de cisaillement peut être utilisée pour déterminer le module de cisaillement. La partie réelle soi-même représente le module de cisaillement. Ceci est représenté sous la forme : _02y(Q) AY(M La propriété de viscoélasticité, telle que la viscosité ou le module de cisaillement, peut être déterminée sur un domaine de fréquences (w). Par exemple, les données de déplacement se trouvent à l'intérieur d'un domaine de 30 à 150 Hz. La viscosité est déterminée sur le même domaine. La viscosité pour des groupes de fréquences ou pour des sous-bandes peut être moyennée, par exemple en prévoyant des valeurs tous les 5 Hz. La propriété est déterminée en réponse à uniquement la forme d'onde modulée en phase et modulée en amplitude et des balayages du tissu pendant l'application de la forme d'onde. Une mise en oeuvre séquentielle des actes 30 à 40 n'est pas prévue mais peut l'être. En utilisant une forme d'onde avec un domaine de fréquences, la réponse se rapportant à la fréquence du tissu peut être mesurée pendant un court intervalle de temps, par exemple suffisant pour entrelacer les balayages et la transmission de la forme d'onde. La viscosité et le module d'élasticité sont des propriétés quantitatives. Les valeurs sont indépendantes d'une grandeur de compression appliquée. La pression au patient pendant l'émission de la forme d'onde acoustique ou en provenance d'une. source extérieure n'est pas nécessitée pour déterminer des propriétés quantitatives. Suivant des modes de réalisation en variante, une propriété qualitative est utilisée, telle que le déplacement.

A l'acte 42, la. pente de la propriété de viscoélasticité en fonction de la fréquence est trouvée. La propriété en fonction de la fréquence peut être tracée point par point, au moins sur une partie ou sur tout le domaine de fréquences (par exemple, de 30 à 150 Hz). Une ligne est adaptée aux différents points ou données, pour fournir une pente. Par exemple, les pentes du module de cisaillement et de la viscosité en fonction de la fréquence sont déterminées. La pente peut indiquer la santé du tissu ou aider au diagnostic. Par exemple, des tissus sains peuvent avoir une pente de viscosité plus horizontale qu'une tumeur. La réponse en fréquence du tissu peut indiquer un tissu sain ou non. D'autres fonctions de la propriété en fonction de la fréquence peuvent être calculées, telles que la variance, la quantité de variation, ou la courbure. A l'acte 44, il est produit une image. L'image représente la propriété. Une indication textuelle ou numérique de la propriété à une fréquence donnée ou calculée à partir de données pour différentes fréquences peut être affichée. Dans d'autres modes de réalisation, l'adaptation de points et/ou de ligne et la valeur de pente sont émises en sortie. La propriété de viscoélasticité est communiquée à l'utilisateur dans l'image. L'image peut être un graphique, tel qu'un ensemble de points de valeurs en fonction de la fréquence. L'image peut comporter en outre une représentation à une, deux, ou trois dimensions de la propriété ou d'une autre information de cisaillement en fonction de l'emplacement ou de l'espace. Par exemple, la vitesse de cisaillement dans toute une région est affichée. La vitesse de cisaillement module la couleur pour des pixels dans une région suivant une image en mode B modulée en échelle de gris. L'image peut représenter une information de déplacement, telle qu'un cisaillement ou des modules (par exemple, les modules de cisaillement) pour différents emplacements. La grille d'affichage peut être différente de la grille de balayage et/ou de la grille pour laquelle des déplacements sont calculés. La couleur, la clarté, la luminosité, la teinte, ou une autre caractéristique peut être modulée en fonction de l'information de cisaillement. Les actes sont répétés pour d'autres lignes de balayage et/ou d'autres profondeurs. Par exemple, les actes 30 à 40 sont effectués de nouveau pour chaque emplacement dans une région à une, deux, ou trois dimensions. La Figure 6 représente un mode de réalisation d'un système 10 pour un calcul d'onde de cisaillement utilisant l'ultrason. L'ultrason produit une onde de cisaillement, et des données de balayage en réponse au tissu qui répond à l'ultrason sont utilisées pour déterminer une propriété. Le système 10 met en oeuvre le procédé de 1a Figure 1 ou d'autres procédés. Le système 10 comporte un formateur 12 de faisceau d'émission, un transducteur 14, un formateur 16 de faisceau de réception, un processeur 18 d'image, un dispositif 20 d'affichage, et une mémoire 22. Des composants supplémentaires, différents ou d'autres composants peuvent être prévus. Par exemple, une entrée d'utilisateur est prévue pour une désignation manuelle ou assistée d'une région d'intérêt pour laquelle des informations de cisaillement doivent être obtenues. A titre d'un autre exemple, un transducteur HIFU supplémentaire est prévu pour traiter le tissu. Le système 10 est un système d'imagerie à ultrason pour diagnostic médical.

Suivant des modes de réalisation en variante, le système 10 est un ordinateur personnel, une station de travail, une station PALS, ou d'autres agencements à un même emplacement ou répartis sur un réseau pour une imagerie d'acquisition en temps réel ou post acquisition.

Suivant encore d'autres modes de réalisation en variante, le système 10 fait partie d'un système de résonance magnétique. Par exemple, le formateur 12 de faisceau d'émission et le transducteur 14 sont prévus pour émettre la forme d'onde pour produire l'onde de cisaillement, mais le formateur de faisceau de réception n'est pas prévu. A la place, des bobines et des aimants de résonance magnétique, comme représenté à la Figure 5, sont prévus avec le processeur 18, la mémoire 22 et l'affichage 20 pour le balayage.

Le formateur 12 de faisceau d'émission est un émetteur à ultrason, une mémoire, un dispositif d'impulsion, un circuit analogique, un circuit numérique, ou une combinaison de ceux-ci. Le formateur 12 de faisceau d'émission peut fonctionner pour produire des formes d'onde pour une pluralité de canaux avec des amplitudes différentes ou relatives, des retards, et/ou des déphasages. Lors de la transmission d'ondes acoustiques à partir du transducteur 14 en réponse aux ondes produites, un ou plusieurs faisceaux sont formés. Les formes d'ondes sont chacune une forme d'onde modulée en phase et en amplitude, mais avec des retards relatifs et une apodisation pour fournir la forme d'onde acoustique souhaitée (par exemple, voir la Figure 3) à la région de foyer. Pour un balayage de déplacement de tissu, une séquence de faisceaux d'émission sont produits pour balayer une région à une, deux ou trois dimensions. Des formats de balayage Sector, Vector (marque déposée), linéaire, ou d'autres formats de balayage peuvent être utilisés. La même région est balayée plusieurs fois, comme représenté à la Figure 4. Le balayage par le formateur 12 de faisceau d'émission est entrelacé ou synchronisé avec la transmission de la forme d'onde modulée en amplitude et en phase par le formateur 12 de faisceau d'émission. Les mêmes éléments du transducteur 14 sont utilisés à la fois pour balayer et produire des ondes de cisaillement, mais différents éléments, transducteurs, et/ou formateurs de faisceau peuvent être utilisés. Le transducteur 14 est un réseau à 1, 1,25, 1,5 1,75 ou 2 dimensions d'éléments à membrane capacitifs ou piézoélectriques. En variante, un seul élément ayant un foyer mécanique est prévu. Le transducteur 14 comporte une pluralité d'éléments pour effectuer une transduction entre des énergies acoustique et électrique. Par exemple, le transducteur 14 est un réseau PZT à une dimension comportant environ 64 à 256 éléments. Le transducteur 14 est connecté au formateur 12 de faisceau d'émission pour convertir des formes d'ondes électriques en des formes d'ondes acoustiques, et est relié au formateur 16 de faisceau de réception pour convertir des échos acoustiques en des signaux électriques. Le transducteur 14 émet de l'énergie acoustique avec une forme d'onde modulée en amplitude et en phase. La formé d'onde est focalisée au niveau d'une région de tissu ou à un emplacement d'intérêt chez le patient. La forme d'onde acoustique est produite en réponse à l'application de la forme d'onde électrique aux éléments de transducteur. Pour un balayage avec un ultrason pour détecter un déplacement, le transducteur émet de l'énergie acoustique et reçoit des échos. Les signaux de réception sont produits en réponse à de l'énergie d'ultrason (échos) incidente sur les éléments du transducteur 14. Le formateur 16 de faisceau de réception comporte une pluralité de canaux comportant des amplificateurs, des dispositifs à retard, et/ou des rotateurs de phase, et un ou plusieurs dispositifs de sommation. Chaque canal est relié à un ou plusieurs éléments transducteurs. Le formateur 16 de faisceau de réception applique des retards relatifs, des phases relatives, et/ou une apodisation pour former un ou plusieurs faisceaux de réception en réponse à chaque émission pour la détection. Le formateur 16 de faisceau de réception émet en sortie des données représentant des emplacements dans l'espace en utilisant les signaux acoustiques reçus. Des retards relatifs et/ou des phases relatives et des sommations de signaux provenant de différents éléments fournissent une formation de faisceau. Suivant des modes de réalisation en variante, le formateur 16 de faisceau de réception est un dispositif de traitement pour produire des échantillons en utilisant une transformée de Fourier ou d'autres transformées. Le formateur 16 de faisceau de réception peut comporter un filtre, tel qu'un filtre pour isoler une information à une seconde harmonique ou à une autre bande de fréquence relative à la bande de fréquence d'émission. De telles informations peuvent plus probablement comporter des tissus souhaités, des agents de contraste, et/ou des informations concernant le débit ou l'écoulement. Suivant un autre mode de réalisation, le formateur 16 de faisceau de réception comporte une mémoire ou un tampon et un filtre ou un dispositif d'addition. Deux ou plus de deux faisceaux de réception sont combinés pour isoler de l'information à une bande de fréquence souhaitée, telle qu'une seconde harmonique, à une fréquence fondamentale du troisième ordre, ou une autre bande. Le formateur 16 de faisceau de réception émet en sortie des données sommées de faisceau représentant des emplacements dans l'espace. Des données pour un emplacement unique, pour des emplacements le long d'une ligne, pour des emplacements dans une zone, ou pour des emplacements dans un volume sont émis en sortie. Une focalisation dynamique peut être prévue. Les données peuvent être pour différents buts. Par exemple, différents balayages sont effectués pour des données de tissu ou des données de mode B plutôt que pour une détection d'onde de cisaillement. En variante, les données en mode B sont également utilisées pour déterminer un déplacement créé par une onde de cisaillement. Le processeur 18 est un détecteur de mode B, un détecteur Doppler, un détecteur Doppler à onde pulsée, un processeur de corrélation, un processeur de transformée de Fourier, un circuit intégré spécifique d'application, 25 un processeur général, un processeur de commande, un processeur d'images, un réseau à grille programmable par champ, un processeur de signaux numériques, un circuit analogique, un circuit numérique, des combinaisons de ceux-ci ou n'importe quel autre dispositif connu aujourd'hui ou qui sera développé ultérieurement pour détecter et traiter des informations à partir d'échantillons à ultrason formés sous la forme de faisceau.

Suivant un mode de réalisation, le processeur 18 comporte un ou plusieurs détecteurs et un processeur distinct. Le processeur distinct est un processeur de commande, un processeur général, un processeur de signaux numériques, un circuit intégré spécifique d'application, un réseau de grille programmable sur champ, un réseau, un serveur, un groupe de processeurs, un trajet de données, des combinaisons de ceux-ci ou n'importe quel autre dispositif connu aujourd'hui ou qui sera développé ultérieurement pour déterminer un déplacement et calculer des propriétés de tissu. Par exemple, le processeur distinct effectue n'importe quelle combinaison de l'un ou plusieurs des actes représentés à la Figure 1. Le processeur 18 est configuré par un système logiciel et/ou matériel pour effectuer les actes.

Suivant un mode de réalisation, le processeur 18 estime un déplacement de tissu au cours du temps en fonction des données émises en sortie provenant du formateur 16 de faisceau de réception. Les déplacements sont estimés sous la forme d'un profil ou de données représentant une courbe d'amplitude de déplacement en fonction du temps. Le profil de déplacement peut être obtenu en corrélant ou en déterminant d'une autre manière un niveau de similarité entre des données de référence et des données obtenues pendant la transmission de la forme d'onde de poussée modulée en phase ou en amplitude ou entrelacée avec celle-ci. Le processeur 18 est configuré pour transformer le déplacement du tissu au cours du temps en le domaine de fréquences. Une transformée de Fourier, telle que la transformée de Fourier Rapide, est appliquée aux données de déplacement. Le processeur 18 est configuré pour calculer une information de cisaillement en fonction du déplacement du tissu au cours du temps. Par exemple, une vitesse de cisaillement est calculée à partir du déplacement au cours du temps. La quantité de déplacement divisée par le temps fournit la vitesse. Suivant un mode de réalisation, le processeur 18 calcule la viscosité en fonction de la fréquence ou à une fréquence. La viscosité est calculée à partir du déplacement dans le domaine de fréquences. Le déplacement transformé est utilisé pour déterminer la viscosité en fonction de la fréquence. Les fréquences se trouvent dans un domaine associé à la forme d'onde modulée en amplitude et en phase. La viscosité à une fréquence, une viscosité moyenne, ou des viscosités à différentes fréquences sont calculées dans le domaine de fréquences en utilisant une représentation du cisaillement, mais peuvent en variante être calculées dans le domaine de temps. Le processeur 18 peut calculer d'autres propriétés, telles que le module de cisaillement. Le processeur 18 produit et émet en sortie des images ou affiche des valeurs cartographiées à partir des propriétés de viscoélasticité sur le dispositif 20 d'affichage. Par exemple, la viscosité de cisaillement, le module de cisaillement ou une autre grandeur est déterminé. Un texte ou une indication numérique de la propriété est affiché pour l'utilisateur. Un graphe de la viscosité ou d'une autre propriété au cours du temps ou en fonction de la fréquence peut être affiché. Suivant un mode de réalisation, le processeur 18 adapte une ligne droite à la propriété en fonction de la fréquence et la pente est affichée en tant qu'une valeur ou de manière graphique. Suivant des modes de réalisations en variante ou supplémentaires, une information de cisaillement est affichée en fonction de l'emplacement. L'amplitude des valeurs de cisaillement module la couleur, la teinte, la clarté, et/ou d'autres caractéristiques d'affichage pour différents pixels représentant une région de tissu. Le processeur 18 détermine une valeur de pixel (par exemple, Rouge, Vert, Bleu) ou une valeur scalaire convertie en une valeur de pixel. L'image est produite en tant que la valeur scalaire ou la valeur de pixel. L'image peut être émise en sortie vers un processeur vidéo, une table de consultation, une carte en couleur, ou directement sur le dispositif 20 d'affichage.

Le processeur 18 fonctionne à la suite d'instructions mémorisées dans la mémoire 22 ou une autre mémoire. Les instructions configurent le processeur 18 pour un fonctionnement en étant chargées et/ou exécutées. Le processeur 18 est programmé pour mesurer une propriété de viscoélasticité du tissu in vivo. La mémoire 22 est un support de stockage pouvant être lu par ordinateur de manière non transitoire. Les instructions pour mettre en oeuvre les processus, les procédés et/ou les techniques qui sont ici décrites sont fournis sur le support ou mémoires de stockage pouvant être lu par ordinateur, par exemple un cache, un tampon, une RAM, un support amovible, un disque dur ou n'importe quel autre support de stockage pouvant être lu par ordinateur. Des supports de stockage pouvant être lus par ordinateur comportent divers types de supports de stockage volatiles et non volatiles. Les fonctions, actes, ou tâches illustrées aux figures ou ici décrites sont exécutées en réponse à un ou plusieurs ensembles d'instructions mémorisées dans un ou sur des supports de stockage pouvant être lus par ordinateur. Les fonctions, actes ou tâches sont indépendantes du type particulier de l'ensemble d'instruction, du support de stockage, des stratégies de traitement ou de processeur et peuvent être effectuées par logiciel, matériel, circuits intégrés, firmware, micro code et analogues, fonctionnant seul ou en combinaison. De même, des stratégies de traitement peuvent comporter un multitraitement, des multitâches, un traitement en parallèle et analogue. Suivant un mode de réalisation, les instructions sont mémorisées sur un dispositif de support amovible pour une lecture par des systèmes à distance ou locaux. Suivant d'autres modes de réalisation, les instructions sont mémorisées dans un emplacement à distance pour un transfert par l'intermédiaire d'un réseau d'ordinateurs ou par des lignes téléphoniques. Suivant encore d'autres modes de réalisation, les instructions sont mémorisées à l'intérieur d'un ordinateur donné, d'une CPU, d'une GPU ou d'un système.

Le dispositif 20 d'affichage est un CRT, LCD, projecteur, plasma, ou autre dispositif d'affichage pour afficher des images en deux dimensions ou en trois dimensions. Le dispositif 20 d'affichage affiche une ou plusieurs images représentant des informations de cisaillement. L'image est un graphique, un nombre, un texte, et/ou des représentations en deux dimensions d'une région. Par exemple, une valeur de viscosité ou un graphe de viscosité en fonction de la fréquence est affiché en tant que l'image.

Bien que l'invention ait été décrite ci-dessus en référence à divers modes de réalisation, il convient de comprendre que de nombreuses variations et modifications peuvent être effectuées sans sortir de l'étendue de protection de l'invention. Par conséquent, la description détaillée précédente doit être considérée uniquement à titre d'illustration et non pas de limitation et il va de soi que ce sont les revendications qui suivent, en y incluant touts les modes de réalisation équivalents, qui ont pour but de définir l'esprit et l'étendue de protection de la présente invention.

Claims (20)

1. REVENDICATIONS1. Support de stockage pouvant être lu par ordinateur de manière non-transitoire ayant mémorisées en son sein des données représentant des instructions pouvant être exécutées par un processeur (18) programmé pour mesurer (32) une propriété de viscoélasticité de tissu in vivo, caractérisé en ce que le support de stockage comporte des instructions pour émettre (30) une forme d'onde modulée en phase et modulée en amplitude comportant des cycles à différentes fréquences à différents instants, lesdites différentes fréquences se trouvant dans un domaine de fréquences, la forme d'onde modulée en phase et modulée en amplitude ayant été émise vers le tissu chez un patient ; calculer (34) un déplacement en fonction du temps du tissu en réponse à la forme d'onde modulée en phase et modulée en amplitude, le déplacement étant calculé à partir de balayages du tissu ; appliquer (36) une transformée de Fourier au déplacement au cours du temps ; et déterminer (40) la propriété de viscoélasticité à partir d'une transformée de Fourier d'une équation d'onde de cisaillement et de la transformée de Fourier du déplacement au cours du temps.
2. 2. Support de stockage pouvant être lu par ordinateur suivant la revendication 1, caractérisé en ce que la forme d'onde modulée en phase et modulée en amplitudecomporte un balayage en fréquence à impulsion chirp à l'intérieur du domaine.
3. 3. Support de stockage pouvant être lu par ordinateur suivant la revendication 1, caractérisé en ce que la forme d'onde modulée en phase et modulée en amplitude comporte une porteuse sinusoïdale multipliée par une modulation d'amplitude sinusoïdale comportant un terme variant avec la phase.
4. 4. Support de stockage pouvant être lu par ordinateur suivant la revendication 1, caractérisé en ce que le calcul (34) du déplacement comporte le fait de balayer le tissu par ultrason entrelacé avec l'émission (30) et le calcul (34) du déplacement du tissu entre les balayages en fonction de la corrélation entre des données provenant des balayages.
5. 5. Support de stockage pouvant être lu par ordinateur suivant la revendication 1, caractérisé en ce que l'étape (34) de calcul du déplacement comporte le fait d'acquérir des données de résonance magnétique d'espace k représentant le tissu à différents instants tandis que le tissu répond à l'émission (30), de soustraire des données de résonance magnétique d'espace k de référence représentant le tissu qui ne répond pas à l'émission (30), et de calculer (34) un déplacement à partir des résultats de la soustraction.
6. 6. Support de stockage pouvant être lu par ordinateur suivant la revendication 1, caractérisé par le fait qu'appliquer (36) la transformée de Fourier comporte le fait de transformer le déplacement au cours du temps enun domaine de fréquence où le déplacement aux différents instants est en réponse aux différentes fréquences.
7. 7. Support de stockage pouvant être lu par ordinateur suivant la revendication 1, caractérisé en ce que la détermination (40) de la propriété de viscoélasticité comporte le fait de déterminer la viscosité.
8. 8. Support de stockage pouvant être lu par ordinateur suivant la revendication 1, caractérisé en ce que l'émission (30) et la détermination (40) comportent le fait de déterminer (40) la propriété de viscoélasticité sur le domaine de fréquences en réponse à uniquement la forme d'onde modulée en phase et modulée en amplitude et à des balayages du tissu pendant l'application de la forme d'onde.
9. 9. Support de stockage pouvant être lu par ordinateur suivant la revendication 1, caractérisé en ce que la détermination (40) comporte le fait de déterminer (40) à la fois la viscosité en fonction d'une partie imaginaire et le module de cisaillement en fonction d'une partie réelle de _la transformée de Fourier de l'équation d'onde de cisaillement.
10. 10. Support de stockage pouvant être lu par ordinateur suivant la revendication 1, caractérisé en ce qu'il comporte en outre le fait de trouver (42) une pente de la propriété de viscoélasticité en fonction de la fréquence dans le domaine.
11. 11. Support de stockage pouvant être lu par ordinateur suivant la revendication 1, caractérisé en ce qu'il comporte en outre un filtrage (38) à bande passante de lasortie de l'application de la transformée de Fourier du déplacement au cours du temps, la bande passante étant égale à environ deux fois le domaine en termes de largeur et des fréquences haute et basse.
12. 12. Procédé pour la mesure de viscoélasticité en utilisant l'ultrason, le procédé étant caractérisé en ce qu'il comporte les étapes dans lesquelles : on émet (30) une forme d'onde à ultrason modulée en phase et en amplitude chez un patient in vivo ; et on mesure (32) une propriété de viscoélasticité à partir d'un cisaillement créé par la forme d'onde à ultrason modulée en phase et en amplitude, la propriété de viscoélasticité étant mesurée sur un domaine de fréquences associé à la forme. d'onde à ultrason modulée en amplitude et en phase et étant indépendante d'une grandeur de compression qui est appliquée extérieurement au patient pendant la transmission (30).
13. 13. Procédé suivant la revendication 12, caractérisé en ce que l'émission (30) comporte le fait de multiplier une porteuse sinusoïdale à une fréquence à ultrason avec une enveloppe ayant la modulation en phase et en amplitude de sorte que le domaine de fréquences se trouve entre 15 Hz et 75 Hz.
14. 14. Procédé suivant la revendication 12, caractérisé en ce que l'étape de mesure (32) comporte : le fait de balayer le patient avec un ultrason 30 entrelacé avec l'émission (30) ; et le fait de calculer (34) le déplacement créé par l'émission (30) du tissu du patient entre des balayages, le calcul (34) étant fonction de la corrélation entre des données provenant des balayages.
15. 15. Procédé suivant la revendication 12, caractérisé en ce que la mesure (32) comporte : le fait d'acquérir des données de résonance magnétique d'espace k représentant du tissu du patient à différents instants tandis que le tissu est en train de répondre à l'émission (30), de soustraire des données de résonance magnétique d'espace k de référence représentant le tissu qui ne répond pas à l'émission (30) et de calculer (34) un déplacement à partir des résultats de la soustraction.
16. 16. Procédé suivant la revendication 12, caractérisé en ce que la mesure (32) comporte le fait d'appliquer (36) une transformée de Fourier à des données représentant un déplacement en fonction du temps et de calculer (34) la viscosité en fonction des résultats de l'application (36) de la transformée de Fourier.
17. 17. Procédé suivant la revendication 16, caractérisé par le fait que le calcul (34) comporte le fait de calculer (34) en fonction d'une partie Imaginaire d'une transformée de Fourier de l'équation d'onde de cisaillement.
18. 18. Procédé suivant la revendication 12, caractérisé en ce qu'il comporte en outre le fait de trouver (42) une pente de la propriété de viscoélasticité en fonction de la fréquence dans le domaine.
19. 19. Système pour le calcul d'onde de cisaillement utilisant l'ultrason, le système étant caractérisé en ce qu'il comporte :30un formateur (12) de faisceau d'émission qui peut fonctionner pour produire une forme d'onde modulée en amplitude et en phase ; un transducteur (14) à ultrason qui est relié au formateur (12) de faisceau d'émission de sorte que le transducteur à ultrason émet de l'énergie acoustique chez un patient en réponse à la forme d'onde modulée en amplitude et en phase ; un formateur (16) de faisceau de réception qui peut fonctionner pour émettre en sortie des données représentant des emplacements dans l'espace en fonction de signaux acoustiques reçus ; un processeur (18) qui est configuré pour estimer le déplacement chez le patient au cours du temps en fonction des données de sortie et pour.calculer l'information de cisaillement en fonction du déplacement chez le patient au cours du temps ; et un dispositif (20) d'affichage qui peut fonctionner pour afficher une image, l'image étant fonction de 20 l'information de cisaillement.
20. 20. Système suivant la revendication 19, caractérisé en ce que processeur (18) est configuré pour transformer le déplacement chez le patient au cours du temps en un 25 domaine de fréquence et pour calculer la viscosité et le module de cisaillement en fonction de la fréquence à partir du déplacement dans le domaine de fréquence, les fréquences dans un domaine étant associées à la forme d'onde modulée en phase et en amplitude.