FR2965934A1 - Support de memorisation, procede et systeme de resolution d'onde transversale en imagerie medicale par ultrason - Google Patents

Support de memorisation, procede et systeme de resolution d'onde transversale en imagerie medicale par ultrason Download PDF

Info

Publication number
FR2965934A1
FR2965934A1 FR1103043A FR1103043A FR2965934A1 FR 2965934 A1 FR2965934 A1 FR 2965934A1 FR 1103043 A FR1103043 A FR 1103043A FR 1103043 A FR1103043 A FR 1103043A FR 2965934 A1 FR2965934 A1 FR 2965934A1
Authority
FR
France
Prior art keywords
transverse
function
tissue
displacements
locations
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
FR1103043A
Other languages
English (en)
Other versions
FR2965934B1 (fr
Inventor
Liexiang Fan
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Siemens Medical Solutions USA Inc
Original Assignee
Siemens Medical Solutions USA Inc
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Siemens Medical Solutions USA Inc filed Critical Siemens Medical Solutions USA Inc
Publication of FR2965934A1 publication Critical patent/FR2965934A1/fr
Application granted granted Critical
Publication of FR2965934B1 publication Critical patent/FR2965934B1/fr
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B8/00Diagnosis using ultrasonic, sonic or infrasonic waves
    • A61B8/48Diagnostic techniques
    • A61B8/485Diagnostic techniques involving measuring strain or elastic properties
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B8/00Diagnosis using ultrasonic, sonic or infrasonic waves
    • A61B8/13Tomography
    • A61B8/14Echo-tomography
    • A61B8/145Echo-tomography characterised by scanning multiple planes
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B8/00Diagnosis using ultrasonic, sonic or infrasonic waves
    • A61B8/54Control of the diagnostic device
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
    • G01S15/00Systems using the reflection or reradiation of acoustic waves, e.g. sonar systems
    • G01S15/88Sonar systems specially adapted for specific applications
    • G01S15/89Sonar systems specially adapted for specific applications for mapping or imaging
    • G01S15/8906Short-range imaging systems; Acoustic microscope systems using pulse-echo techniques
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
    • G01S7/00Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
    • G01S7/52Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S15/00
    • G01S7/52017Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S15/00 particularly adapted to short-range imaging
    • G01S7/52023Details of receivers
    • G01S7/52036Details of receivers using analysis of echo signal for target characterisation
    • G01S7/52042Details of receivers using analysis of echo signal for target characterisation determining elastic properties of the propagation medium or of the reflective target
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B8/00Diagnosis using ultrasonic, sonic or infrasonic waves
    • A61B8/13Tomography
    • A61B8/14Echo-tomography
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B8/00Diagnosis using ultrasonic, sonic or infrasonic waves
    • A61B8/48Diagnostic techniques
    • A61B8/488Diagnostic techniques involving Doppler signals

Landscapes

  • Engineering & Computer Science (AREA)
  • Health & Medical Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Physics & Mathematics (AREA)
  • Remote Sensing (AREA)
  • Radar, Positioning & Navigation (AREA)
  • Medical Informatics (AREA)
  • Public Health (AREA)
  • Biomedical Technology (AREA)
  • Heart & Thoracic Surgery (AREA)
  • Pathology (AREA)
  • Molecular Biology (AREA)
  • Surgery (AREA)
  • Animal Behavior & Ethology (AREA)
  • General Health & Medical Sciences (AREA)
  • Radiology & Medical Imaging (AREA)
  • Veterinary Medicine (AREA)
  • Biophysics (AREA)
  • Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
  • General Physics & Mathematics (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • Acoustics & Sound (AREA)
  • Ultra Sonic Daignosis Equipment (AREA)

Abstract

Système pour résoudre de l'information d'onde transversale dans de l'imagerie médicale par ultrason, le système étant caractérisé en ce qu'il comprend : un formeur (16) de faisceau de réception pouvant fonctionner pour sortir des données représentant des emplacements dans l'espace en fonction des signaux acoustiques reçus ; un processeur (18) configuré pour estimer des déplacements de tissu en fonction des données de sortie, calculer les coefficients de corrélation pour au moins une paire des déplacements de tissu à partir d'emplacements différents dans l'espace, résoudre un temps de propagation sur la base des coefficients de corrélation et produire une image en fonction du temps de propagation, et un affichage (20) pouvant fonctionner pour afficher l'image.

Description

Résolution d'onde transversale dans de l'imagerie médicale par ultrason Arrière plan technologique La présente invention se rapporte à l'imagerie par ultrason. En particulier, on détermine de l'information d'onde transversale pour de l'imagerie par ultrason. On utilise la détection de la vitesse transversale dans 10 diverses applications de diagnostic, telles que l'évaluation d'une maladie du foie. La vitesse transversale peut indiquer les propriétés du tissu en permettant une détection de tumeurs ou d'autres régions. Mais la précision de détection peut être insuffisante 15 pour être à l'aise dans d'autres applications. Pour effectuer une détection d'onde transversale on amorce une impulsion de poussée d'ultrason le long d'une ligne de balayage. L'impulsion de poussée produit une onde transversale provoquant un déplacement du tissu. On 20 détecte le déplacement. Pour détecter la vitesse de l'onde transversale, on peut utiliser de multiples impulsions de poussée le long d'une même ligne de balayage. Même où de multiples impulsions de poussée ne sont pas amorcées le long d'une même ligne de balayage, 25 les impulsions de poussée sont amorcées séquentiellement le long de lignes de balayage voisines. Des additions d'ondes transversales et de résonances de tissu provoquées par des poussées multiples au même emplacement ou en des emplacements voisins peuvent se traduire par 30 des artéfacts ou par un manque de précision suffisante pour certaines applications.5 Résumé succinct A titre d'introduction, le mode de réalisation préféré décrit ci-dessous englobe des procédés, des instructions et des systèmes de résolution de l'information d'onde transversale dans l'imagerie médicale par ultrason. Plutôt que de déterminer une information d'onde transversale pour chaque emplacement sur la base d'une distance à l'origine de l'onde transversale et du temps jusqu'au déplacement de pic pour l'emplacement, on met des profils de déplacement provenant de diverses combinaisons d'origine et d'emplacement de détection en corrélation. On détecte une information transversale en utilisant des déplacements à partir d'une combinaison spatiale d'emplacements de transmission d'emplacements de détection. On utilise la combinaison de corrélation dans une même fonction pour résoudre simultanément les retards d'emplacements latéraux multiples. La diversité spatiale et la corrélation en couche peuvent fournir une estimation d'onde transversale plus précise. Suivant une première facette de l'invention, un support de mémoire non transitoire, pouvant être lu par un ordinateur et ayant mémorisé en son sein des instructions représentant des données et pouvant être exécutées par un processeur programmé pour résoudre de l'information d'onde transversale en imagerie médicale par ultrason, est caractérisé en ce que le support de mémoire comprend des instructions pour déterminer un premier déplacement de tissu en fonction du temps provoqué par une première onde transversale en réaction à une première transmission acoustique le long d'une première ligne de balayage, le premier déplacement étant sur une quatrième ligne de 2 balayage espacé de la première ligne de balayage par au moins des deuxième et troisième lignes de balayage, déterminer un deuxième déplacement de tissu en fonction du temps provoqué par une deuxième onde transversale en réaction à une deuxième transmission acoustique le long de la troisième ligne de balayage, le deuxième déplacement étant sur la quatrième ligne de balayage, déterminer un troisième déplacement de tissu en fonction du temps provoqué par une troisième onde transversale en réaction à une troisième transmission acoustique le long de la première ligne de balayage, le troisième déplacement étant sur la deuxième ligne de balayage, calculer pour chacun d'une pluralité d'emplacements latéraux un retard correspondant à un maximum de corrélation des premier et deuxième déplacements, du premier et du troisième déplacements et du deuxième et du troisième déplacements, calculer une vitesse transversale dans du tissu pour chacun de la pluralité des emplacements latéraux en fonction du retard respectif.
De préférence :
- déterminer le premier, le deuxième et troisième déplacements en fonction du temps comprend suivre une 25 amplitude du mouvement du tissu en fonction du temps.
- les première, deuxième et troisième et quatrième lignes de balayage sont des lignes de balayage voisines.
30 - calculer le retard correspondant à un maximum de corrélations comprend calculer en fonction des emplacements latéraux et d'une pluralité de retards possibles. calculer le retard correspondant au maximum des corrélations comprend calculer une somme des corrélations des premier et deuxième déplacements, des premier et troisième déplacements et des deuxième et troisième déplacements.
- calculer une vitesse transversale comprend calculer la vitesse transversale en fonction du retard et de l'emplacement latéral. - calculer le retard comprend résoudre simultanément pour une pluralité de lignes de balayage en même temps.
- le support comprend en outre des instructions pour 15 répéter les actes de détermination et de calcul pour d'autres lignes de balayage et d'autres profondeurs.
- le support comprend en outre des instructions pour transmettre les premières, deuxième et troisième 20 transmissions acoustiques dans un ordre évitant une transmission séquentielle immédiatement le long de ligne de balayage voisine.
Suivant une deuxième facette, il est prévu un procédé de 25 résolution de l'information d'onde transversale dans l'imagerie médicale par ultrason, le procédé étant caractérisé en ce qu'on détecte des premières ondes transversales dans du tissu en un même emplacement provoquées par des transmissions en des emplacements 30 différents. On détecte des deuxièmes ondes transversales dans du tissu en des emplacements différents provoquées par des transmissions en un même emplacement. On calcule une vitesse transversale en fonction à la fois des première et deuxième ondes transversales détectées. On10 produit une image représentant le tissu en fonction de la vitesse transversale.
De préférence : - détecter les première et deuxième ondes transversales comprend déterminer des déplacements provoqués par les transmissions. - calculer la vitesse transversale comprend calculer la vitesse transversale pour chacune d'une pluralité de positions latérales par corrélation des déplacements à partir des emplacements différents. - il comprend en outre entrelacer spatialement les transmissions à travers les lignes de balayage et dans lequel calculer comprend répartir spatialement les ondes transversales et la durée dans le temps des ondes transversales pour chacun d'une pluralité d'emplacements latéraux. - calculer la vitesse transversale comprend calculer des coefficients de corrélation pour au moins une paire de profils temporels de déplacement d'onde transversale, la au moins une paire de profils temporels de déplacement d'onde transversale étant l'une des premières ondes transversales avec l'une des deuxièmes ondes transversales, deux des premières ondes transversales, ou deux des deuxièmes ondes transversales, le calcul étant une fonction d'une distance entre des lignes de balayage. - il comprend en outre détecter des troisièmes ondes transversales dans du tissu en des emplacements différents provoqués par des transmissions à d'autres emplacements différents, dans lequel calculer est une fonction des premières, deuxièmes et troisièmes ondes transversales détectées.
Suivant une troisième facette, il est prévu un système pour résoudre de l'information d'onde transversale dans de l'imagerie médicale par ultrason, le système étant caractérisé en ce qu'il comprend un formeur de faisceau de réception pouvant fonctionner pour sortir des données représentant des emplacements dans l'espace en fonction de signaux acoustiques reçus. Un processeur configuré pour estimer des déplacements de tissu en fonction des données de sortie, calculer les coefficients de corrélation pour au moins une paire des déplacements de tissu à partir d'emplacements différents dans l'espace, résoudre pour un temps de propagation sur la base des coefficients de corrélation et produire une image en fonction du temps de propagation, et un affichage pouvant fonctionner pour afficher l'image. De préférence : - le processeur est configuré pour résoudre pour les différents emplacements dans l'espace en même temps. le processeur est configuré pour résoudre pour les divers emplacements dans l'espace avec un jeu de sorties déterminé à partir d'une fonction utilisant des données des divers emplacements dans l'espace, chaque valeur du jeu de sortie étant une fonction des déplacements du tissu pour les divers emplacements dans l'espace. - le processeur est configuré pour calculer les coefficients de corrélation pour au moins trois paires de déplacements de tissu, une première paire étant des déplacements dans du tissu en un même emplacement provoqués par des transmissions en des emplacements différents, une deuxième paire étant des déplacements dans du tissu en des emplacements différents provoqués par des transmissions en un même emplacement et une troisième paire étant des déplacements dans du tissu à des emplacements différents provoqués par des transmissions à d'autres emplacements différents. - le processeur est configuré pour calculer une vitesse transversale à partir du temps de propagation, l'image étant produite pour représenter la vitesse transversale. D'autres caractéristiques et avantages de l'invention sont passés en revue ci-dessous en liaison avec les modes de réalisation préférés. Description succincte des dessins Les éléments et les figures ne sont pas nécessairement à l'échelle, l'accent étant mis au contraire sur l'illustration des principes de l'invention. En outre, dans les figures de mêmes repères désignent des parties correspondantes dans les diverses vues.
La figure 1 est un schéma synoptique d'un mode de réalisation d'un procédé de résolution de l'information d'onde transversale en imagerie médicale par ultrason ; La figure 2 est une illustration graphique d'une séquence donnée à titre d'exemple d'une séquence de transmission 20 et de détection de déplacement ; Les figures 3A à 3C illustrent des séquences de transmission différentes et la figure 3D est un graphe donné à titre d'exemple montrant les effets de températures des différentes séquences de transmission 25 des figures 3A à 3C ; La figure 4 est un graphique montrant un exemple d'échantillonnage de déplacement axial, latéral et temporel ; 2965934 s La figure 5 est un diagramme donné à titre d'exemple ou représentant des relations de plages et de retards ; et La figure 6 illustre un mode de réalisation d'un système de résolution d'une information d'onde transversale en 5 imagerie médicale par ultrason. Description détaillée des dessins et des modes de réalisation préférés actuellement On détecte de l'information transversale en utilisant des déplacements à partir d'une combinaison spatiale 10 d'emplacements de transmission et d'emplacements de détection. Pour éviter une ambiguïté sur l'onde transversale introduite par une excitation répétée (impulsion de poussée) sur le même emplacement spatial, on utilise une séquence d'entrelacement. En raison de 15 l'entrelacement, la séquence peut aussi réduire le chevauchement de faisceau en réduisant l'élévation de température. La séquence peut avoir une intensité de signal d'onde transversale remarquable par rapport à des transmissions multiples le long de la même ligne de balayage. L'intensité peut permettre une détection qui est plus fiable et plus robuste. Au lieu de détecter l'instant de pic de l'onde transversale, on détermine le temps de propagation en utilisant une corrélation de la réaction de l'onde transversale à partir des emplacements qui sont associés à la stratégie de séquençage. On traite le problème d'estimation du temps de propagation comme un problème d'optimisation dans l'espace de coefficient de corrélation, ce qui se traduit par une solution optimum par rapport à une recherche du temps jusqu'au pic ou par rapport à la résolution de l'équation d'Helmholtz.
Dans un mode de réalisation, on détermine au moins une valeur de distance pour un type donné de tissu et un dispositif donné de système d'imagerie. On utilise un entrelacement spatial des impulsions d'excitation pour minimiser l'effet thermique du tissu. On diversifie les emplacements de répartition spatiale d'une impulsion de détection d'onde transversale. On acquière des données en utilisant les divers emplacements avec des transmissions entrelacées, ce qui permet une détermination par corrélation du retard ou du temps de propagation de l'onde transversale. Dans un autre mode de réalisation, on calcule les coefficients de corrélation pour au moins une paire de profils temporels de déplacement d'onde transversale. On estime le temps de propagation de l'onde transversale entre chaque distance donnée à partir des coefficients de corrélation. On inverse le temps de propagation et on multiplie le résultat par la distance pour obtenir la vitesse transversale. D'autres modes de réalisation sont possibles. La figure 1 représente un procédé de résolution d'une onde transversale en imagerie médicale par ultrason. Le procédé est mis en oeuvre par le système de la figure 6 ou par un système différent. Des actes supplémentaires différents ou moins nombreux peuvent être prévus. C'est ainsi par exemple que l'acte 35 n'est pas effectué. Au lieu de cela on détermine une autre caractéristique d'onde transversale. Comme autre exemple, l'acte 36 n'est pas effectué. L'itération ou la répétition représentée par les flèches de l'acte 35 à l'acte 30 peut ne pas être prévue dans d'autres modes de réalisation. Les actes sont effectués dans l'ordre décrit ou représenté mais peuvent être effectués dans des ordres différents.
Dans l'acte 30, des impulsions acoustiques sont transmises pour produire une onde transversale. Les excitations sont transmises à partir d'un transducteur d'ultrason. Des excitations sont de l'énergie acoustique.
L'énergie acoustique est focalisée, ce qui se traduit en un faisceau ou en plusieurs faisceaux pour chaque transmission. Les excitations sont focalisées en utilisant un réseau à commande de phase et/ou un foyer mécanique. On peut défocaliser les excitations dans une seule dimension telle que la dimension d'élévation. On émet les excitations dans du tissu d'un patient. On focalise les excitations en un emplacement pour permettre de détecter l'onde transversale qui s'ensuit, tel qu'un emplacement de tissu entourant et incluant une tumeur éventuelle. On focalise les transmissions en des emplacements différents. Pour balayer une région d'un patient pour produire des ondes transversales, on focalise une transmission ou plusieurs transmissions en un même emplacement. On focalise d'autres transmissions à d'autres emplacements. On produit l'onde transversale dans la région focale et on la propage latéralement à partir de la région focale. On détecte l'onde transversale en des emplacements voisins et/ou à distance de la région focale pour toute transmission donnée. L'onde transversale se réduit en amplitude au fur et à mesure que l'onde se propage dans le tissu. Pour détecter une réaction du tissu à des ondes transversales dans une région à laquelle on s'intéresse, on effectue des transmissions vers d'autres régions focales et on effectue une détection autour des autres régions focales. Les régions de détection de région focale différentes de transmission se chevauchent. C'est ainsi par exemple que la figure 2 montre dix lignes de balayage allant de n à n+10. On peut prévoir plus de lignes de balayage ou moins de lignes de balayage. Deux ou plusieurs transmissions pour produire des ondes transversales sont effectuées le long de chaque ligne de balayage comme représenté par les flèches pour les impulsions de poussée en 40 et 42. Les lignes de balayage pour des impulsions de poussée sont les mêmes que les lignes de balayage pour la détection. Dans des variantes de réalisation, une ou plusieurs des lignes de balayage d'impulsion de poussée sont différentes des lignes de balayage de réception ou de détection. Pour produire une onde transversale on souhaite avoir des excitations de grandes amplitudes ou puissances. C'est ainsi par exemple que chaque excitation a un indice mécanique proche de 1,9 mais sans le dépasser. Pour être prudent et prendre en considération des variations d'échantillons, on peut utiliser un indice mécanique de 1,7 comme limite supérieure. On peut utiliser des puissances plus grandes (par exemple indice mécanique dépassant 1,9) ou plus petites. En transmettant séquentiellement le long des mêmes lignes de balayages ou de lignes de balayage voisines, il peut se faire que le tissu voit sa température s'élever avec le temps. Des effets biologique peuvent comprendre une hyperthermie à une température du tissu d'environ 41 à 45°C, une dénaturation de protéine à des températures supérieures à 45°C, et une nécrose du tissu à des températures supérieures à 50°C. La rigidité du tissu peut être affectée même à des températures inférieures à 45°C. A des températures supérieures à 45°C une augmentation de viscosité et/ou de rigidité peut se produire. A des températures supérieures à 50°C le tissu peut avoir une grande rigidité et/ou une grande atténuation. Les transmissions peuvent provoquer des effets biologiques.
En variante, des effets biologiques sont limités en empêchant une élévation de température de plus de 2 degrés Celsius. Les figures 3A à 3C représentent trois séquences différentes de transmission pour produire des ondes transversales en des emplacements différents. La figure 3A représente une transmission séquentielle le long de lignes de balayage voisines. Les lignes de balayage sont représentées sous la forme de flèches pointant vers le bas. L'ordre de transmission est numéroté de 1 à 5. On peut utiliser moins de lignes de balayage ou davantage de lignes de balayage. Les figures 3B et 3C présentent des séquences de transmission donnant moins d'élévation de température avec le temps. On effectue des transmissions séquentielles dans un ordre évitant des transmissions séquentielles de manière immédiate le long de lignes de balayage voisines en entrelaçant spatialement les transmissions à travers des lignes de balayage. La figure 3B illustre des excitations d'amorçage le long de cinq lignes de balayage mais en le faisant dans une séquence gauche/droite. Le champ de vue est subdivisé en moitiés, les moitiés droite et gauche. La première transmission s'effectue le long d'une ligne de balayage du côté gauche de la moitié gauche. La deuxième transmission s'effectue le long d'une ligne de balayage du côté gauche de la moitié droite. Cette configuration se continue en se terminant par une avant dernière transmission du côté droit de la moitié gauche et une dernière transmission du côté droit de la moitié droite. Dans toute la séquence, les transmissions voisines temporellement dans la séquence sont espacées d'une moitié du champ de vue pour une détection transversale (par exemple région à laquelle on s'intéresse). La figure 3C représente une autre séquence possible. Dans cet exemple le champ de vue est subdivisé en quartiers.
Chaque quartier balayé séquentiellement suivant un cycle entrelacé. On peut utiliser d'autres séquences de balayage avec d'autres divisions, répartitions au hasard ou configurations temporelles/spatiales. La figure 3D illustre une élévation de température associée aux trois séquences différentes des figures 3A à 3C. L'élévation de la température du tissu est représentée dans le champ proche où les champs acoustiques de faisceaux se chevauchent beaucoup. La séquence 1 (ligne en tiret) de la figure 3A a une élévation de température plus grande en fonction du temps. Les séquences 2 et 3 (lignes en traits mixtes et en traits pleins respectivement) des figures 3B et 3C ont une élévation de température réduite. On peut régler, en variante ou supplémentairement, la température dans l'estimation de l'onde transversale en résolvant la vitesse de l'onde transversale en utilisant une corrélation supplémentaire. Pour procurer cette corrélation supplémentaire permettant une résolution simultanée du retard en des emplacements latéraux multiples, on utilise une séquence d'acquisition en phase avec un nombre limité, tel que 2, 3 ou 4, d'impulsions de poussée au même emplacement. En 40 et 42, la figure 2 représente deux phases d'impulsions de poussée ou des transmissions pour produire des ondes transversales. Au lieu de pousser plusieurs fois (par exemple cinq, dix ou plus) au même emplacement avec un intervalle de répétition des impulsions court (PRI), on utilise deux ou plusieurs phases d'acquisition. Chaque phase 40, 42 balaie les lignes de balayage suivant une configuration quelconque telle que représentée aux figures 3A à 3C. En utilisant seulement deux phases, l'excitation PRI pour le même emplacement spatial peut être très petite (de l'ordre de 100ms). D'autres transmissions se produisent. Ces autres transmissions sont destinées à détecter les ondes transversales plutôt qu'à provoquer l'onde transversale.
Les transmissions de détection peuvent avoir une puissance plus petite et/ou des impulsions courtes et utiliser les mêmes lignes de balayage que les impulsions de poussée ou des lignes de balayage différentes pour produire les ondes transversales. Les transmissions pour la détection peuvent avoir un profil de faisceau plus large le long d'au moins une dimension, tel que latéralement, pour former des échantillons de réception le long d'une pluralité de lignes de balayage. Dans l'acte 32, on détermine un déplacement du tissu.
L'onde transversale fait que le tissu se déplace. La figure 2 représente le déplacement pour contrôler l'onde transversale sous la forme d'une flèche en 44. Les flèches sont représentées dans une seule direction, mais l'onde transversale se propage dans des directions multiples. L'onde transversale peut être contrôlée dans une direction, deux directions, ou plusieurs directions. Le déplacement du tissu provoqué par l'onde transversale est déterminé en fonction du temps. Lorsque l'onde transversale est produite, du tissu voisin n'est pas déplacé jusqu'à ce que la contrainte transversale se propage. Comme l'onde transversale passe en un emplacement donné, le tissu se déplace d'une quantité ou d'une distance qui augmente vers une quantité de pic, puis diminue alors que le tissu revient au repos. La quantité de pic est associée à un retard ou à un temps de propagation à partir de la production de l'onde transversale au point focal où le pic de l'onde transversale passe à un autre emplacement. Ce profil de déplacement en fonction du temps est déterminé en suivant l'amplitude du mouvement du tissu en fonction du temps. On détecte le déplacement par un balayage par ultrason.
Pour détecter le déplacement, on transmet de l'énergie ultrasonore au tissu subissant un déplacement et on reçoit des réflexions de l'énergie. La figure 2 représente cette détection en 46. On effectue la transmission et la réception pour la détection plusieurs fois pour déterminer un changement dû au déplacement. On peut utiliser toute séquence de transmission et de réception. On peut entrelacer la détection du déplacement avec d'autres balayages, tel qu'un balayage de région différente pour un déplacement séparément.
On contrôle une région à laquelle on s'intéresse pour détecter l'onde transversale. La région à laquelle on s'intéresse peut avoir n'importe quelle dimension, tel que 6 mm latéralement et 10 mm axialement. On contrôle cette région de détection par ultrason. On effectue par exemple, des balayages en mode B pour détecter un déplacement du tissu provoqué par l'onde transversale. On peut utiliser un mode Doppler, de flux de couleur, ou d'autres modes ultrasonores pour contrôler l'onde transversale.
On effectue le contrôle pour n'importe quel nombre de lignes de balayage. On forme par exemple, quatre faisceaux de réception en réaction à chaque transmission.
Après avoir transmis l'excitation pour produire l'onde transversale, on effectue des transmissions en mode B d'une manière répétitive le long d'une ligne de balayage unique et des réceptions le long de quatre lignes de balayage voisines. Dans d'autres modes de réalisation, on forme seulement un seul faisceau de réception ou d'autres nombres de faisceaux de réception en réaction à chaque transmission. On peut utiliser n'importe quel nombre de répétitions, tel que d'environ 10 à 120 fois. Certaines des données ultrasonores, tel que le début ou la fin des répétitions, peuvent ne pas être sensibles à l'onde transversale. Au fur et à mesure que l'onde transversale se propage dans les lignes de balayage, l'intensité en mode B peut varier. Pour les lignes de balayage contrôlées, il est prévu une séquence de données qui représente un profil dans le temps du mouvement du tissu provenant de l'onde transversale. Par exemple, des données d'une pluralité d'emplacements dans l'espace (par exemple le long des lignes de balayage) sont mises en corrélation avec une trame de référence de donnée ou d'image en fonction du temps. On peut utiliser toute détection d'élasticité. Pour chaque profondeur ou emplacement dans l'espace, on effectue une corrélation sur une pluralité de profondeurs ou d'emplacements (par exemple noyau de 64 profondeurs, la profondeur centrale étant le point où on calcule le profil). On peut utiliser des déplacements dans l'espace en deux ou en trois dimensions. On peut utiliser un déplacement en une seule dimension le long d'une direction différente des lignes de balayage. Le décalage spatial ayant la corrélation la plus grande ou une corrélation suffisante à un instant donné indique la quantité de déplacement. On détermine des déplacements pour un emplacement donné à des instants différents. Le profil temporel pour un emplacement donné indique une détection de l'onde transversale. On examine le profil en ce qui concerne l'absence de bruit ou une instance unique de variation. Un pic dans le profil avec ou sans filtrage passe-bas temporel indique le passage du front de l'onde transversale. On choisit le déplacement le plus grand, mais on peut utiliser statistiquement le déplacement moyen ou un autre déplacement. La figure 2 représente en 48 la réaction temporelle de déplacement détectée. Pour contrôler une région plus grande, on forme des faisceaux de réception supplémentaires en réaction à la même impulsion de poussée. En variante, on produit une autre onde transversale et on fournit les faisceaux de transmission et faisceaux de réception à une distance différente du point de production de l'onde transversale. En utilisant la solution décrite ici, on obtient les déplacements associés à des emplacements latéraux multiples en utilisant seulement deux impulsions de poussée ou plusieurs impulsions de poussée par ligne de balayage. Pour chaque emplacement de faisceau de réception, un profil en fonction du temps d'une information de mouvement est fourni, en étant représenté par les données d'ultrason. On évite des transmissions le long de lignes de balayage différentes pour contrôler une même onde transversale pendant la formation du profil temporel pour obtenir une résolution temporelle plus grande, mais on peut prévoir des positions de balayage entrelacées ou décalées.
Ce que l'on indique ci-dessus s'applique à une seule profondeur. L'échantillonnage peut être conçu pour fournir une grille recouvrant toute l'étendue axiale de la région à laquelle on s'intéresse. Dans un autre mode de réalisation, on obtient des échantillons à des profondeurs multiples pour chaque faisceau reçu. Un profil de temps distinct est prévu pour chaque profondeur axiale ainsi que chaque emplacement latéral. On peut utiliser n'importe quel nombre de profondeurs, tel que environ 200 pour 5 mm ou 400 pour 10 mm. On obtient des données d'ultrason représentant des emplacements différents dans la région à laquelle on s'intéresse. On obtient les données d'ultrason en temps réel avec le balayage ou on les obtient à partir d'une mémoire. Pour chaque emplacement, l'information de déplacement représente la réaction à des instants différents. D'autres balayages, contrôles ou techniques, peuvent être utilisés pour obtenir les données d'ultrason afin d'estimer une amplitude transversale. Utiliser des déplacements dans le temps et un emplacement latéral permet d'estimer une vitesse transversale avec aussi peu que deux impulsions de poussée amorcées par ligne de balayage. On détecte par exemple des ondes transversales dans du tissu en un même emplacement provoquées par des transmissions provenant d'emplacements différents ; on détecte des ondes transversales en des emplacements différents provoquées par des transmissions en un même emplacement ; et on détecte des ondes transversales dans du tissu à des emplacements différents provoquées par des transmissions à d'autres emplacements différents. Cette diversité spatiale et temporelle permet une résolution pour l'onde transversale avec moins d'ambiguïté.
La figure 2 représente l'utilisation de cette diversité d'emplacement. Dans cette stratégie d'impulsion, une force d'excitation est appliquée deux fois à chaque emplacement d'indice (par exemple ligne de balayage) pour produire des ondes transversales, et la séquence d'impulsions de détection est appliquée aussi deux fois pour détecter l'onde transversale à d'autres emplacements. Cependant, chaque séquence d'impulsions de détection correspond à des emplacements différents, ce qui se traduit par un système ayant trois jeux possibles de solutions de détection d'onde transversale. Un jeu correspond au séquençage détecté au même emplacement, mais excité en des emplacements différents, le deuxième jeu correspond au séquençage détecté en des emplacements différents, mais excité à partir du même emplacement et le troisième jeu correspond à la détection et à l'excitation, tous deux à des emplacements différents.
L'idée est illustrée à la figure 2 sous la forme des paramètres p et q de décalage d'indice et de la liaison entre deux formes de tube en 46. On peut utiliser une diversité temporelle et spatiale plus grande ou plus petite, tels qu'en prévoyant des emplacements et/ou des temps et des combinaisons associés supplémentaires. Une zone hétérogène est illustrée à la figure 2 sous la forme de l'ovale des réactions de déplacement détectées en 48. Les réactions temporelles de déplacement sont différentes en des emplacements à l'intérieur et à l'extérieur de la zone hétérogène. L'hétérogénéité du milieu fait qu'aucun de chacun des jeux uniques de la solution ci-dessus n'est parfait. La figure 4 illustre un exemple d'emplacements spatiaux et temporels de réactions temporelles de déplacement et leurs excitations correspondantes. A, B, et C sont des zones transversales produites à des emplacements d'excitation correspondants (lignes en tiret). A', B', et C' sont des zones transversales telles que détectées en des emplacements de détection. Les emplacements de détection sont sur des lignes de balayage équidistantes, mais peuvent avoir d'autres répartitions. Cet exemple est celui de trois impulsions de poussée pour des ondes A, B, et C transversales et des déplacements A', B', et C' détectés associés de tissu en fonction du temps. La figure 4 représente le processus à une épaisseur. L'axe horizontal représente une distance latérale lorsque des impulsions de poussée pour des ondes A et C transversales sont transmises le long d'une première ligne de balayage et impulsion de poussée pour l'onde B transversale est transmise le long d'une troisième ligne de balayage. Un autre axe horizontal représente le temps en montrant des impulsions de poussée pour des ondes A et C transversales transmises le long de la première ligne de balayage se produisant à des instants différents. C'est ainsi par exemple que des impulsions de poussée pour des ondes A et C transversales correspondent à des impulsions 40, 42 de poussée le long d'une même ligne n+3 de balayage de la figure 2. Une impulsion de poussée pour l'onde B transversale est représentée comme étant transmise au même moment, mais cela est relatif. La transmission réelle se produit en séquence le long de la troisième ligne de balayage (par exemple ligne n+6 de balayage à la figure 2). La synchronisation de la détection du déplacement en réaction à des impulsions A et B de poussée est représentée comme se produisant en même temps pour mieux indiquer la relation des déplacements A' et B' détectés sur la même quatrième ligne de balayage.
Dans l'exemple des figures 2 et 4, le déplacement A' provoqué par l'onde A transversale en réaction à la transmission acoustique le long de la ligne n+3 de balayage est détecté. Ce déplacement A' est sur la ligne n+6 de balayage représenté par la flèche à partir de la deuxième impulsion de poussée à la figure 2. On effectue la détection en faisant intervenir deux lignes de balayage. Il n'est pas prévu d'autres lignes de balayage. Dans d'autres modes de réalisation, on utilise seulement une ligne de balayage ou plus de deux lignes de balayage qui interviennent. On détecte le déplacement B' provoqué par l'onde B transversale en réaction à la transmission acoustique le long de la ligne n+5 de balayage. Ce déplacement B' sur la ligne n+6 de balayage, immédiatement voisine de la ligne n+5 de balayage, est utilisé pour l'impulsion de poussée pour produire l'onde B transversale. On peut prévoir des lignes de balayage supplémentaires qui interviennent. Le déplacement C' provoqué par l'onde C transversale en réaction à la transmission acoustique le long de la même ligne n+3 de balayage en tant qu'impulsion de poussée pour l'onde A transversale est détecté. On détecte ce déplacement C' le long de la ligne n+4 de balayage immédiatement voisine, mais on peut prévoir une ou plusieurs lignes de balayage qui interviennent. La figure 2 représente une flèche allant de l'impulsion 40 de poussée sur la ligne n+3 de balayage à l'impulsion 46 de détection sur la ligne n+4 de balayage pour représenter C et C'. On peut prévoir d'autres relations spatiales et temporelles. Les figures 2 et 4 en sont des exemples. On peut utiliser une autre impulsion de poussée le long de la ligne n+3 et/ou n+5 de balayage. Une détection de déplacement supplémentaire le long de la même ligne de balayage ou de la ligne de balayage différente, avec ou sans transmission d'impulsion de poussée supplémentaire, peut être utilisée (par exemple détection C" sur n+5 en plus de C' sur n+4 en réaction à la même impulsion C de poussée). Les figures 2 à 4 représentent une détection se produisant sur des lignes de balayage à distance de la droite des transmissions, mais tout écartement (par exemple vers la gauche en deux dimensions ou avant/arrière en trois dimensions) peut être utilisé. A la figure 4, AA' et BB' forment une paire pour détecter le paramètre de temps de propagation dans la plage _ 1. AA' et CC' forment une autre paire pour détecter le paramètre de temps de propagation dans la plage _ 2, respectivement. BB' et CC' forment une autre paire pour limiter la solution. Les paramètres de plage représentent une distance connue entre des lignes de balayage. Le temps de propagation sur la plage _ 1 et la plage 2 peut être déterminé comme les retards 1 et 2 à partir des paires spatiales/temporelles. On peut utiliser l'information de retard pour déterminer la vitesse ou le temps de propagation de l'onde transversale. Dans l'acte 34, on utilise la diversité spatiale pour calculer les retards. On met les déplacements en corrélation pour calculer le retard. En considérant une plage de retard possible différent, on utilise des profils de déplacement pour la résolution des retards associés à des emplacements différents dans un même calcul. Pour un emplacement donné, on calcule le retard correspondant à une corrélation maximum de divers déplacements dans l'espace. Plutôt que d'effectuer une seule corrélation, le maximum est une fonction de corrélations multiples, telles qu'associées à des paires différentes de déplacements. La vitesse transversale ou une autre information transversale est calculée en fonction des ondes transversales détectées en utilisant une corrélation. Lorsque les déplacements sont déterminés en utilisant une corrélation, tel que pour un suivi de tâches, on utilise une corrélation supplémentaire ou en couche pour déterminer une information transversale à partir des déplacements. On détermine un temps de propagation d'onde classique par la fonction suivante : z- _ max cet-) , d (m, n, t)d (u, v, t + 2) (1) art{T} dans laquelle r est le retard, c est une corrélation, t est un temps, d est le déplacement, et m, n, u, et v sont des coordonnées spatiales en deux dimensions. C'est ainsi par exemple que les deux fonctions de déplacement sont les courbes ou profils A' et B' de déplacement. On peut effectuer cette corrélation de maximum d'une manière adéquate en un rapport signal à bruit relativement grand et pour des formes d'onde de signal moins déformées, mais il peut ne pas être adéquat dans d'autres situations.
Pour compenser ce qui se passe dans la zone transitoire, telle que représentée à la figure 2 comme étant la région de tissu hétérogène, on utilise une diversité spatiale. On utilise des paires multiples de corrélations pour déterminer le retard. On peut utiliser la diversité spatiale et temporelle pour déterminer une répartition dans l'espace des ondes transversales et la durée dans le temps des ondes transversales pour chacune d'une pluralité d'emplacements latéraux. On utilise un comportement de groupe de corrélation croisée pour chercher la corrélation maximum. Pour chaque emplacement latéral, y, une plage de retard possible est utilisée pour déterminer le retard associé à la corrélation maximum. On peut obtenir la fonction de corrélation de toute l'étendue latérale à la même profondeur, sous la forme d'une matrice en deux dimensions c (y, r). Le retard pour chacun d'une pluralité d'emplacements latéraux est résolu simultanément. On peut résoudre le retard pour la vitesse transversale pour une pluralité de lignes de balayage en même temps comme partie d'une même fonction. Une corrélation des déplacements à partir des emplacements différents est calculée pour au moins une paire de profils temporels de déplacement d'ondes transversales. On utilise une fonction pour sortir une matrice de retard pour des emplacements correspondants.
On dénote la réaction temporelle de déplacement par d(i,j,t), les indices i et j représentant l'indice d'emplacement d'excitation (par exemple la ligne de balayage) et l'indice d'emplacement de suivi (par exemple ligne de balayage), respectivement. On utilise une corrélation pour détecter le retard d'impulsion entre deux paires de courbes d(i ;j ;t) de déplacement. C'est-à-dire {d(i,i+p,t), d(i,i+q),t)}, {d(i,i+p,t), d(i+pq,i+p,t)}. Une autre paire {d(i,i+q),t), d(i+p-q,i+p,t)} contribue à la détermination de la zone hétérogène. Dans l'exemple de la figure 4, une paire est A'B', une autre paire est A'C', et la paire contribuant à la détermination de la zone hétérogène est B'C'.
On peut utiliser d'autres paires, telles que là où l'on effectue des transmissions et/ou des détections d'impulsion de poussée supplémentaires. On peut utiliser toutes paires diverses temporellement, spatialement ou spatialement ou temporellement de déplacement. Le calcul est une fonction d'une distance entre des lignes de balayage, tel que fourni en détectant un déplacement d'impulsion de poussée le long d'une même ligne de balayage à des lignes de balayage différent. On utilise la plage entre des lignes de balayage pour déterminer la plage. On utilise la plage avec le retard pour déterminer la vitesse de l'onde transversale.
Une fonction à titre d'exemple calcule une somme des corrélations de trois paires différentes de déplacement. On met des ondes transversales détectées provenant de transmissions différentes le long d'une même ligne de balayage et/ou de lignes de balayage de détection différente en corrélation les unes avec les autres. En faisant l'hypothèse que T. conserve une variation continue dans la même profondeur, la valeur de r pour la même profondeur et le long de toute l'étendue latérale dénotée comme ligne frontière ou courbe zT est obtenue en résolvant la fonction de coût défini par : 20 avec: r-~ r (Y, 'r') d.(i,i+p,t)cad(i+p ,i+ ,~n La courbe est représentée sous la forme d'une matrice de retards r et d'emplacement y latéral, dans laquelle T.« y) 25 est une première, deuxième, troisième ou autre dérivé de T. en fonction de y, a,p,y étant des paramètres de pondération. Ces paramètres de pondération peuvent être fixés à un nombre réel 1.0, fixés à d'autres nombres, ou peuvent être une fonction du T. On peut utiliser toute 7 L~(%,a+p,t)d(i+ p.._ f>i+ 4i -tP,1)dr(i., o -r) (3) plage et nombre de valeurs pour le retard, tels que 32 valeurs possibles de r sur une plage d'unité considérée comme allant de 0 à 100. On résout la courbe de coefficient de corrélation maximum pour chaque profondeur en combinant toute l'information provenant de la propagation de l'onde transversale disponible. Dans des variantes, on peut utiliser des termes ou des corrélations supplémentaires de paire. On peut utiliser différentes combinaisons de diversité.
En utilisant l'équation (2) ou une autre fonction de corrélation ayant des termes multiples de paires diverses, on peut réduire la variance due à l'excitation à deux emplacements différents ayant un module différent ; on peut traiter d'une manière adaptative une hétérogénéité du milieu en introduisant le terme c3 de corrélation (par exemple corrélation de B'C') ; et on peut réduire l'erreur en appliquant les contraintes de voisinage.
Comme représentée à la figure 4, l'estimation du retard r est différente pour la paire de donnée utilisée dans le calcul de la fonction de corrélation. La donnée de paire {d(i,i+p,t), d(i,i+q),t)} (corrélation A'B') estime le retard moyen de temps dans la plage _ 1, et la paire {d(i,i+p,t), d(i+p-q,i+p,t)} (corrélation A'C') estime le retard moyen de temps dans la plage _ 2. Il y a une plage partagée dans ces deux estimations. L'estimation (r) fine entre chaque groupe de lignes voisines est obtenue en développant zi dans la pleine étendue latérale, les valeurs connues aux deux extrémités étant données. La pleine étendue latérale peut s'étendre sur toute la région à laquelle on s'intéresse (par exemple dix lignes de balayage), mais peut s'étendre dans un sous-jeu de la région. On obtient les valeurs connues par le séquençage légèrement différent des impulsions d'excitation et de détection qui est spécifique aux extrémités. Aux deux extrémités du champ de vue pour l'imagerie, la détection C' est au même emplacement que C, et l'excitation B est au même emplacement que l'impulsion B' de suivi, de sorte que le temps estimé de propagation est la valeur moyenne allant de A à A'. On utilise les équations suivantes pour résoudre la valeur de r de l'intervalle entre chaque groupe de faisceaux : ém. té
l:a extrémité Le séquençage aux deux extrémités donne l'estimation de 15 ro et de n_1.
La figure 5 représente le milieu de l'équation 4. Le procédé pour développer les valeurs de retard r étant donné les valeurs connues au point ro et TN-1 d'extrémité 20 peut être décrit comme des équations de solution récursives. ro étant donné dans l'ordre de bas en haut, les équations suivantes donnent : a =b(ri) birr -1) 1) 0 _1), On peut calculer tous les a(n) et b(n), de sorte que l'équation (4) se réduit à deux variables inconnues. Ce (5) 25 a(n) système d'équations étant donné et le fait que 2.N_1 est connu, on peut résoudre l'équation (5) à partir de l'ordre de l'indice de haut en bas, ce qui donne T. pour chaque groupe de faisceaux voisins. On utilise le même procédé pour résoudre le T individuel pour une acquisition de faisceaux parallèles. Par exemple, lorsque quatre faisceaux parallèles indexés de 0 à 4 sont utilisés dans des impulsions de suivi, chacun des faisceaux indexés peut être résolu par les équations (4) à (6). Le coefficient pondéré dans l'équation (4) est modifié en prenant en compte la distance différente entre chaque faisceau indexé.
Dans l'acte 34, on calcule la vitesse transversale en fonction du retard et de l'emplacement latéral. On détecte la vitesse transversale pour les divers emplacements dans l'espace du tissu. On calcule la vitesse transversale dans le tissu pour chacun de la pluralité des emplacements latéraux en fonction du retard respectif.
On obtient la vitesse transversale à partir du temps de propagation ou du retard (T) et de la distance connue (par exemple, plage) entre des lignes de balayage. Le temps et la distance à l'emplacement déterminent la vitesse. On connaît la distance par l'intervalle entre les lignes de balayage (c'est-à-dire, la position du faisceau de transmission pour produire l'onde transversale et la position du faisceau de réception pour détecter l'onde transversale). On peut estimer tout module ou valeur transversale en variante ou en plus de la vitesse transversale. Des valeurs de module de tissu représentent la dureté ou la rigidité aux emplacements. C'est ainsi par exemple que l'on estime le module transversal du tissu. Dans des variantes, on estime un module d'Young. Dans d'autres modes de réalisation, on estime d'autres valeurs transversales, qu'elles soient quantitatives ou qualitatives.
Le module transversal est donné par g = pv2s , dans laquelle p est la masse volumique et Vs est une vitesse transversale estimée. Dans un mode de réalisation, les modules du tissu ou l'information transversale, tel que le module transversal, sont déterminés en fonction de la contrainte ou du déplacement ou de l'information transversale. C'est ainsi par exemple que l'on détermine le module transversal pour chaque emplacement d'échantillon en résolvant par itération une équation de diffusion. En faisant l'hypothèse d'un rapport de Poisson de 0.5 ou en utilisant un rapport de Poisson connu, on calcule par itération le module transversal en des emplacements différents sous la forme d'une fonction du champ de contrainte à des instants différents ou sous une contrainte différente pour les emplacements différents et le module transversal différent.
On répète les actes pour d'autres lignes de balayage et/ou d'autres profondeurs. C'est ainsi par exemple qu'on effectue les actes 30, 32, 34 et 34 à nouveau pour chacune d'une de plusieurs profondeurs. On peut effectuer à nouveau les actes pour d'autres emplacements à distance latéralement. La figure 4 représente par exemple une partie de la séquence globale. On fournit les diverses corrélations pour divers emplacements latéraux en résolvant l'équation (3). En variante, on traite d'autres parties latérales séparément, tels que dans des groupes uniques de quatre lignes de balayage. En variante, on peut utiliser l'impulsion de poussée et les paires de déplacement détectés ou d'information pour d'autres solutions ou parties de la région à laquelle on s'intéresse. On peut par exemple utiliser B et B' comme paire C et C' pour une solution.
Dans l'acte 36, on produit une image. L'image représente le tissu en fonction de la vitesse transversale ou d'autres informations transversales. La région pour l'information transversale peut être un sous-jeu ou une région à laquelle on s'intéresse dans une image représentant une zone plus grande ou un volume plus grand du patient. La vitesse transversale module par exemple une couleur pour des pixels dans une région d'une image en mode B modulée en échelle de gris. L'image peut représenter une information de déplacement, tels. que des modules (par exemple les modules transversaux) pour des emplacements différents. 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. On module la couleur, la brillance, la luminance, la teinte ou d'autres caractéristiques en fonction de l'information transversale.
La figure 6 représente un mode de réalisation d'un système 10 pour résoudre l'information d'onde transversale dans une imagerie médicale par ultrason. Le système 10 met en oeuvre le procédé de la figure 1, de la figure 2, de la figure 4 ou d'autres procédés. Le système 10 comprend un formeur 12 de réseau de transmission, un transducteur 14, un formeur 16 de faisceau de réception, un processeur 18 d'image, un affichage 20, et une mémoire 22. On peut prévoir des éléments supplémentaires différents ou moins d'éléments. Une entrée utilisateur est prévue par exemple pour désigner manuellement ou d'une manière assistée une région à laquelle on s'intéresse et pour laquelle il faut obtenir une information transversale. Comme autre exemple, il est prévu un transducteur 24 HIFU supplémentaire pour traiter le tissu. Le système 10 est un système d'imagerie de diagnostic médical par ultrason. Dans des variantes, le système 10 est un ordinateur personnel, un poste de travail, un poste PACS ou d'autres agencements en un même endroit ou répartis sur un réseau d'imagerie en temps réel ou en post-acquisition.
Le formeur 12 de faisceau de transmission est un transmetteur à ultrason, une mémoire, un pulseur, un circuit analogique, un circuit numérique, ou leurs combinaisons.. Le formeur 12 de faisceau de transmission peut fonctionner pour produire des formes d'onde pour une pluralité de canaux ayant des amplitudes différentes ou relatives, des retards et/ou un déphasage. Par transmission d'ondes acoustiques à partir du transducteur 14 en réaction aux ondes produites, il se forme un ou plusieurs faisceaux. Une séquence de faisceaux de transmission est produite pour balayer une région en deux dimensions ou en trois dimensions. On peut utiliser des formats Sector, Vector (marque de fabrique), linéaire ou d'autres formats de balayage. On peut balayer la même région plusieurs fois. Pour une imagerie de flux ou Doppler et pour une imagerie de contrainte, on utilise une séquence de balayage. Dans l'imagerie Doppler, la séquence peut comprendre des faisceaux multiples le long d'une même ligne de balayage avant de balayer une ligne de balayage voisine. Pour une imagerie transversale, on peut utiliser n'importe quelle séquence, telle que représentée aux figures 3A à 3C. On peut former un, deux, ou plusieurs faisceaux de transmission pour produire des ondes transversales le long de chaque ligne de balayage avant de balayer une autre ligne de balayage ou on peut les entrelacer. La séquence d'acquisition est par exemple séparée du calcul de corrélation. Une mémoire mémorise des données pour fournir l'information de déplacement pour la corrélation. Dans des variantes, le formeur 12 de faisceau de transmission produit une onde plane ou une onde divergente pour un balayage plus rapide.
Le transducteur 14 est un réseau à une dimension 1,25, GYM dimension 1,5, une dimension 1,75 ou une dimension 2 d'éléments piézoélectriques ou à membrane capacitive. Le transducteur 14 comprend une pluralité d'éléments pour une transduction entre des énergies acoustique et électrique. Le transducteur 14 est par exemple un réseau PZT en deux dimensions (par exemple ayant environ 3.000 éléments). Les signaux de réception sont produits en réaction à de l'énergie ultrasonore (échos) arrivant sur les éléments du transducteur 14. Les éléments relient des canaux des formeurs 12, 16 de faisceau de transmission et de réception.
Le formeur 16 de faisceau de réception comprend une pluralité de canaux ayant des amplificateurs, des retards et/ou des rotateurs de phase et un ou plusieurs additionneurs. Chaque canal se relie à l'un ou à plusieurs éléments de transducteur. Le formeur 16 de faisceau de réception applique des retards relatifs, des phases et/ou une apodisation pour former un ou plusieurs faisceaux de réception en réaction à chaque transmission pour une détection. Le former 16 de faisceau de réception sort des données représentant des emplacements dans l'espace en utilisant les signaux acoustiques reçus. Des retards relatifs et/ou des déphasages et une sommation de signaux provenant d'éléments différents procurent une formation de faisceau. Dans des variantes le formeur 16 de faisceau de réception est un processeur pour produire des échantillons utilisant une transformée de Fourier ou d'autres transformées.
Le formeur 16 de faisceau de réception peut comprendre un filtre, tel qu'un filtre pour isoler de l'information à un deuxième harmonique ou à une autre bande de fréquence par rapport à la bande de fréquence de transmission. Une information de ce genre peut comprendre plus probablement du tissu souhaité, de l'agent de contraste et/ou une information de flux. Dans une autre variante, le formeur 16 de faisceau de réception comprend une mémoire ou un tampon et un. filtre ou un additionneur. On combine deux ou plusieurs faisceaux reçus pour isoler de l'information à une bande de fréquence souhaitée, telle qu'un deuxième harmonique, un fondamental cubique ou une autre bande.
Le formeur 16 de faisceau de réception sort une donnée sommée de faisceau représentant des emplacements dans l'espace. Une donnée pour un emplacement unique, des emplacements le long d'une ligne, des emplacements pour une zone ou des emplacements pour un volume sont sortis.
Une focalisation dynamique peut être prévue. La donnée peut être destinée à des buts différents. On effectue par exemple des balayages différents pour une donnée en mode B ou de tissu autre que pour une détection d'angle transversal. En variante, la donnée en mode B est utilisée aussi pour déterminer un déplacement provoqué par une onde transversale.
En variante, le mode B est utilisé aussi pour déterminer un déplacement provoqué par une onde transversale.
Le processeur 18 est un détecteur en mode B, un détecteur Doppler, un détecteur Doppler à onde pulsée, un processeur de corrélation, un processeur à transformée le Fourier, un circuit intégré spécifique à une application, un processeur général, un processeur de commande, un processeur d'images, une matrice de grille programmable par champ, un processeur à signal numérique, un circuit analogique, un circuit numérique, leurs combinaisons ou d'autres dispositifs connus maintenant ou à développer ultérieurement pour détecter et traiter de l'information pour l'affichage à partir d'échantillons à ultrason formé par faisceau.
Dans un mode de réalisation, le processeur 18 comprend un ou plusieurs détecteurs et un processeur distinct. Le processeur distinct est un processeur de commande, un processeur général, un processeur du signal numérique, un circuit intégré spécifique à une application, une matrice de portes programmable par champ, un réseau, un serveur, un groupe de processeurs, un trajet de données ou leurs combinaisons ou d'autres dispositifs connus maintenant ou qui seront développés ultérieurement pour déterminer un déplacement et calculer les propriétés du tissu. Le processeur distinct effectue par exemple toute combinaison d'un ou plusieurs des actes représentés à la figure 1. Le processeur 18 est configuré par logiciel et/ou par matériel pour effectuer les actes.
Dans un mode de réalisation, le processeur 18 estime des déplacements du tissu en fonction de la donnée de sortie du formeur 16 de faisceau de réception. On estime les déplacements sous la forme d'un profil ou de données représentant une courbe d'amplitude du déplacement en fonction du temps. On peut obtenir un profil de déplacement en mettant en corrélation ou en déterminant d'une autre façon un niveau de similarité entre une donnée de référence et une donnée obtenue pendant le passage de l'onde transversale. Le processeur 18 utilise l'information de déplacement provenant d'emplacements multiples. On détermine des coefficients de corrélation pour au moins une paire de déplacements du tissu à partir de différents emplacements dans l'espace. On utilise les coefficients de corrélation provenant de paires différentes pour résoudre un temps de propagation. On résout le temps de propagation pour des emplacements différents dans l'espace en même temps, par exemple en les déduisant d'une équation ayant les mêmes entrées. Le temps de propagation pour des emplacements différents dans l'espace est envoyé dans un jeu de sortie déterminé à partir d'une fonction utilisant une donnée à partir des différents emplacements dans l'espace. Chaque valeur du jeu de sortie est une fonction des déplacements du tissu pour une pluralité d'emplacements dans l'espace y compris des emplacements de l'espace ne définissant pas la distance pour laquelle le temps de propagation est calculé. Les déplacements pour des emplacements différents dans l'espace contribuent à un temps de propagation donné. On calcule par exemple les coefficients de corrélation pour au moins trois paires de ces déplacements du tissu : une première paire étant des déplacements du tissu au même emplacement provoqué par des transmissions en des emplacements différents, une deuxième paire étant des déplacements du tissu en des emplacements différents provoqués par des transmissions en un même emplacement et une troisième paire étant des déplacements du tissu en des emplacements différents provoqués par des transmissions à d'autres emplacements différents.
Le processeur 18 calcule une vitesse transversale à partir du temps de propagation. On peut déterminer une autre information transversale. On détermine les retards associés à des lignes de balayage particulières. On calcule la vitesse transversale à partir des retards et des distances ou plages correspondantes.
Le processeur 18 produit et sort une image ou des valeurs d'affichage cartographiées à partir des propriétés du tissu pour les envoyer à l'affichage 20. On détermine par exemple la vitesse transversale, le module transversal ou d'autres valeurs pour chaque emplacement. L'amplitude des valeurs module la couleur, la teinte, la brillance et/ou d'autres caractéristiques d'affichage. Le processeur 18 détermine une valeur de pixel (par exemple rouge, vert, bleu) ou une valeur scalaire transformée en une valeur de pixel. On produit l'image sous la forme des valeurs scalaires ou de pixel. On peut sortir l'image sur un processeur vidéo, une table, une carte en couleur ou directement sur l'affichage 20.
Le processeur 18 opère suivant des instructions mémorisées dans la mémoire 22 ou dans une autre mémoire. Le processeur 18 est programmé pour résoudre l'information d'onde transversale en imagerie médicale par ultrason. La mémoire 22 est un support de mémoire non transitoire et pouvant être lu par un ordinateur. Les instructions pour la mise en oeuvre des processus, des procédés et/ou techniques mentionnés dans le présent mémoire sont prévus sur les supports de mémoire pouvant être lus par ordinateur ou sur des mémoires tel qu'un cache, un tampon, une RAM, des supports amovibles ou d'autres supports de mémoire à disque dur ou autres pouvant être lus par ordinateur. Des supports de mémoire pouvant être lus par d'autres ordinateurs englobent divers types de support de mémoire rémanente et non rémanente. Les fonctions, actes ou les tâches illustrés dans les figures ou décrits dans le présent mémoire sont exécutés en réaction à un ou à plusieurs jeux d'instructions mémorisés dans un ou sur un support de mémoire pouvant être lu par un ordinateur. Les fonctions, actes ou tâches sont indépendantes du type particulier de jeu d'instruction, du support de mémoire, du processeur de stratégie de traitement et peuvent être effectués par logiciel, par matériel, par circuit intégré, par produit de firme, par.microcode et analogue, opérant seul ou en combinaison. De même, des stratégies de traitement peuvent comprendre un multitraitement, une multitâche, un traitement parallèle et analogue. Dans un mode de réalisation, les instructions sont mémorisées sur un dispositif de support amovible pour lecture par des systèmes locaux ou éloignés. Dans d'autres modes de réalisation, les instructions sont mémorisées en un emplacement éloigné pour transfert par l'intermédiaire d'un réseau informatique ou par des lignes de téléphone.
Dans d'autres modes de réalisation encore les instructions sont mémorisées dans un ordinateur donné, un CPU, un GPU ou un système.
L'affichage 20 est un CRT, un LCD, un projecteur, un plasma, ou un autre affichage pour afficher des images en deux dimensions ou des représentations en trois dimensions. L'affichage 20 affiche une ou plusieurs images représentant une information transversale, telle que la vitesse transversale.

Claims (20)

  1. REVENDICATIONS1. Support (22) de mémoire non transitoire, pouvant être lu par un ordinateur et ayant mémorisé en son sein des instructions représentant des données et pouvant être exécutées par un processeur (18) programmé pour résoudre de l'information d'onde transversale en imagerie médicale par ultrason, caractérisé en ce que le support (22) de mémoire comprend des instructions pour : déterminer (32) un premier déplacement de tissu en fonction du temps provoqué par une première onde transversale en réaction à une première transmission acoustique le long d'une première ligne de balayage, le premier déplacement étant sur une quatrième ligne de balayage espacée de la première ligne de balayage par au moins des deuxième et troisième lignes de balayage ; déterminer (32) un deuxième déplacement de tissu en fonction du temps provoqué par une deuxième onde transversale en réaction à une deuxième transmission acoustique le long de la troisième ligne de balayage, le deuxième déplacement étant sur la quatrième ligne de balayage ; déterminer (32) un troisième déplacement de tissu en fonction du temps provoqué par une troisième onde transversale en réaction à une troisième transmission acoustique le long de la première ligne de balayage, le troisième déplacement étant sur la deuxième ligne de balayage ; calculer (34) pour chacun d'une pluralité 30 d'emplacements latéraux un retard correspondant à un maximum de corrélation des premier et deuxièmedéplacements, du premier et du troisième déplacements et du deuxième et du troisième déplacements ; calculer (35) une vitesse transversale dans du tissu pour chacun de la pluralité des emplacements latéraux en 5 fonction du retard respectif.
  2. 2. Support (22) de mémorisation non transitoire et pouvant être lu par un ordinateur suivant la revendication 1, caractérisé en ce que déterminer (32) le 10 premier, le deuxième et troisième déplacements en fonction du temps comprend suivre une amplitude du mouvement du tissu en fonction du temps.
  3. 3. Support (22) de mémorisation non transitoire et 15 pouvant être lu par un ordinateur suivant la revendication 1 ou 2, caractérisé en ce que les première, deuxième et troisième et quatrième lignes de balayage sont des lignes de balayage voisines. 20
  4. 4. Support (22) de mémorisation non transitoire et pouvant être lu par un ordinateur suivant l'une des revendications précédentes, caractérisé en ce que calculer (34) le retard correspondant à un maximum de corrélations comprend calculer (34) en fonction des 25 emplacements latéraux et d'une pluralité de retards possibles.
  5. 5. Support (22) de mémorisation non transitoire et pouvant être lu par un ordinateur suivant la 30 revendication 1 à 3, caractérisé en ce que calculer (34) le retard correspondant au maximum des corrélations comprend calculer (34) une somme des corrélations des premier et deuxième déplacements, des premier ettroisième déplacements et des deuxième et troisième déplacements.
  6. 6. Support (22) de mémorisation non transitoire et pouvant être lu par un ordinateur suivant l'une des revendications précédentes, caractérisé en ce que calculer (35) une vitesse transversale comprend calculer (35) la vitesse transversale en fonction du retard et de l'emplacement latéral.
  7. 7. Support suivant l'une quelconque des revendications précédentes, caractérisé en ce que calculer (34) le retard comprend résoudre simultanément pour une pluralité de lignes de balayage en même temps.
  8. 8. Support (22) de mémorisation non transitoire et pouvant être lu par un ordinateur suivant l'une des revendications précédentes, caractérisé en ce qu'il comprend en outre répéter les actes de détermination (32) et de calcul (34) pour d'autres lignes de balayage et d'autres profondeurs.
  9. 9. Support (22) de mémorisation non transitoire et pouvant être lu par un ordinateur suivant l'une des revendications précédentes, caractérisé en ce qu'il comprend en outre transmettre (30) les première, deuxième et troisième transmissions acoustiques dans un ordre évitant une transmission séquentielle immédiatement le long de ligne de balayage voisine.
  10. 10. Procédé de résolution de l'information d'onde transversale dans l'imagerie médicale par ultrason, le procédé étant caractérisé en ce que :on détecte (32) des premières ondes transversales dans du tissu en un même emplacement provoquées par des transmissions en des emplacements différents ; on détecte (32) des deuxièmes ondes transversales 5 dans du tissu en des emplacements différents provoquées par des transmissions en un même emplacement ; on calcule (35) une vitesse transversale en fonction à la fois des premières et deuxièmes ondes transversales détectées ; 10 on produit (36) une image représentant le tissu en fonction de la vitesse transversale.
  11. 11. Procédé suivant la revendication 10, caractérisé en ce que détecter (32) les première et deuxième ondes 15 transversales comprend déterminer des déplacements provoqués par les transmissions.
  12. 12. Procédé suivant la revendication 11, caractérisé en ce que calculer (35) la vitesse transversale comprend 20 calculer (35) la vitesse transversale pour chacune d'une pluralité de positions latérales par corrélation des déplacements à partir des emplacements différents.
  13. 13. Procédé suivant l'une quelconque des revendications 25 10 à 12, caractérisé en ce qu'il comprend en outre entrelacer spatialement les transmissions à travers les lignes de balayage et dans lequel calculer (34) comprend répartir spatialement les ondes transversales et la durée dans le temps des ondes transversales pour chacun d'une 30 pluralité d'emplacement latéraux.
  14. 14. Procédé suivant l'une quelconque des revendications 10 à 13, caractérisé en ce que calculer (35) la vitesse transversale comprend calculer des coefficients decorrélation pour au moins une paire de profils temporels de déplacement d'onde transversale, la au moins une paire de profils temporels de déplacement d'onde transversale étant l'une des premières ondes transversales avec l'une des deuxièmes ondes transversales, deux des premières ondes transversales, ou deux des deuxièmes ondes transversales, le calcul étant une fonction d'une distance entre des lignes de balayage.
  15. 15. Procédé suivant l'une quelconque des revendications 10 à 14, caractérisé en ce qu'il comprend en outre détecter (32) des troisièmes ondes transversales dans du tissu en des emplacements différents provoqués par des transmissions à d'autres emplacements différents, dans lequel calculer (35) est une fonction des premières, deuxièmes et troisièmes ondes transversales détectées.
  16. 16. Système pour résoudre de l'information d'onde transversale dans de l'imagerie médicale par ultrason, le système étant caractérisé en ce qu'il comprend : un formeur (16) de faisceau de réception pouvant fonctionner pour sortir des données représentant des emplacements dans l'espace en fonction des signaux acoustiques reçus ; un processeur (18) configuré pour estimer des déplacements de tissu en fonction des données de sortie, calculer les coefficients de corrélation pour au moins une paire des déplacements de tissu à partir d'emplacements différents dans l'espace, résoudre un temps de propagation sur la base des coefficients de corrélation et produire une image en fonction du temps de propagation, et un affichage (20) pouvant fonctionner pour afficher l'image.
  17. 17. Système suivant la revendication 16, caractérisé en ce que le processeur (18) est configuré pour résoudre pour les différents emplacements dans l'espace en même temps.
  18. 18. Système suivant la revendication 17, caractérisé en ce que le processeur (18) est configuré pour résoudre pour les divers emplacements dans l'espace avec un jeu de sortie déterminé à partir d'une fonction utilisant des données des divers emplacements dans l'espace, chaque valeur du jeu de sortie étant une fonction des déplacements du tissu pour les divers emplacements dans l'espace.
  19. 19. Système suivant l'une des revendications 16 à 18, caractérisé en ce que le processeur (18) est configuré pour calculer les coefficients de corrélation pour au moins trois paires déplacement de tissu, une première paire étant des déplacements dans du tissu en un même emplacement provoqué par des transmissions des emplacements différents, une deuxième paire étant des déplacements dans du tissu en des emplacements différents provoqués par des transmissions en un même emplacement et une troisième paire étant des déplacements dans du tissu à des emplacements différents provoqués par des transmissions à d'autres emplacements différents.
  20. 20. Système suivant l'une des revendications 16 à 19, caractérisé en ce que le processeur (18) est configuré pour calculer une vitesse transversale à partir du temps de propagation, l'image étant produite pour représenter la vitesse transversale.
FR1103043A 2010-10-06 2011-10-06 Support de memorisation, procede et systeme de resolution d'onde transversale en imagerie medicale par ultrason Active FR2965934B1 (fr)

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
US12/898,957 US8961418B2 (en) 2010-10-06 2010-10-06 Solving for shear wave information in medical ultrasound imaging

Publications (2)

Publication Number Publication Date
FR2965934A1 true FR2965934A1 (fr) 2012-04-13
FR2965934B1 FR2965934B1 (fr) 2021-10-29

Family

ID=45872583

Family Applications (1)

Application Number Title Priority Date Filing Date
FR1103043A Active FR2965934B1 (fr) 2010-10-06 2011-10-06 Support de memorisation, procede et systeme de resolution d'onde transversale en imagerie medicale par ultrason

Country Status (5)

Country Link
US (1) US8961418B2 (fr)
JP (2) JP5773833B2 (fr)
KR (2) KR101586998B1 (fr)
DE (1) DE102011114943B4 (fr)
FR (1) FR2965934B1 (fr)

Families Citing this family (31)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20110245668A1 (en) * 2010-04-05 2011-10-06 Tadashi Tamura Methods and apparatus for ultrasound imaging
JP6067590B2 (ja) * 2011-02-25 2017-01-25 メイヨ フォンデーシヨン フォー メディカル エジュケーション アンド リサーチ 非合焦超音波による超音波振動法
US20120253194A1 (en) * 2011-03-30 2012-10-04 Tadashi Tamura Methods and apparatus for ultrasound imaging
US8532430B2 (en) * 2011-07-28 2013-09-10 General Electric Company Methods for reducing motion artifacts in shear wave images
US9239373B2 (en) * 2011-11-16 2016-01-19 Siemens Medical Solutions Usa, Inc. Adaptive image optimization in induced wave ultrasound imaging
JP6058295B2 (ja) * 2012-06-19 2017-01-11 東芝メディカルシステムズ株式会社 超音波診断装置、医用画像処理装置、医用画像処理方法、および医用画像処理プログラム
JP5771758B2 (ja) * 2012-12-25 2015-09-02 日立アロカメディカル株式会社 超音波診断装置
US9332962B2 (en) 2013-03-13 2016-05-10 Siemens Medical Solutions Usa, Inc. Ultrasound ARFI displacement imaging using an adaptive time instance
US10743814B2 (en) 2013-03-15 2020-08-18 Siemens Medical Solutions Usa, Inc. Fat fraction estimation using ultrasound with shear wave propagation
CN104640506B (zh) * 2013-09-18 2017-06-30 东芝医疗***株式会社 超声波诊断装置、医用图像处理装置和医用图像处理方法
US10512450B2 (en) 2013-09-24 2019-12-24 Siemens Medical Solutions Usa, Inc. Shear wave estimation from analytic data
US10278671B2 (en) 2013-09-30 2019-05-07 Siemens Medical Solutions Usa, Inc. Shear wave detection in medical ultrasound imaging
JP6243249B2 (ja) * 2014-02-19 2017-12-06 ジーイー・メディカル・システムズ・グローバル・テクノロジー・カンパニー・エルエルシー 超音波診断装置及びプログラム
JP6289225B2 (ja) * 2014-04-07 2018-03-07 キヤノンメディカルシステムズ株式会社 超音波診断装置及び制御プログラム
WO2015168400A1 (fr) * 2014-04-30 2015-11-05 The Cleveland Clinic Foundation Mesure des propriétés biomécaniques dans une image d'oct
JP6266133B2 (ja) * 2014-05-16 2018-01-24 コーニンクレッカ フィリップス エヌ ヴェKoninklijke Philips N.V. 超音波剪断波エラストグラフィにおける自己相関ガイド相互相関
WO2015180069A1 (fr) 2014-05-28 2015-12-03 深圳迈瑞生物医疗电子股份有限公司 Système et procédé d'imagerie échographique
EP3160358B1 (fr) * 2014-06-30 2020-08-05 Koninklijke Philips N.V. Élastographie par onde transversale échographique comprenant une surveillance de thérapie
US11058397B2 (en) * 2015-01-02 2021-07-13 Esaote S.P.A. Method for quantifying the elasticity of a material by ultrasounds
US10327738B2 (en) * 2015-03-31 2019-06-25 Samsung Medison Co., Ltd. Ultrasound imaging apparatus and method of processing ultrasound image thereof
US10376242B2 (en) * 2015-04-16 2019-08-13 Siemens Medical Solutions Usa, Inc. Quantitative viscoelastic ultrasound imaging
US10729404B2 (en) * 2015-08-03 2020-08-04 Koninklijke Philips N.V. Ultrasound system and method for measurement using shear wave
US11717256B2 (en) * 2016-01-08 2023-08-08 Siemens Medical Solutions Usa, Inc. Motion independence in acoustic radiation force impulse imaging
US11006928B2 (en) * 2016-02-10 2021-05-18 Siemens Medical Solutions Usa, Inc. Sound speed imaging using shear waves
JP6604272B2 (ja) * 2016-06-06 2019-11-13 コニカミノルタ株式会社 超音波診断装置及び超音波信号処理方法
CN109073612B (zh) 2016-07-29 2022-03-29 霍尼韦尔国际公司 光电离检测器自动化零级校准
US10646202B2 (en) * 2017-01-23 2020-05-12 Siemens Medical Solutions Usa, Inc. Sheer speed imaging using coherence
US11553901B2 (en) 2017-04-06 2023-01-17 Siemens Medical Solutions Usa, Inc. Liver disease activity estimation with ultrasound medical imaging
US11523774B2 (en) 2017-04-06 2022-12-13 Siemens Medical Solutions Usa, Inc. Tissue property estimation with ultrasound medical imaging
US11452503B2 (en) 2018-05-18 2022-09-27 Siemens Medical Solutions Usa, Inc. Shear wave imaging based on ultrasound with increased pulse repetition frequency
CN210376230U (zh) * 2019-02-28 2020-04-21 霍尼韦尔国际公司 用于监测目标气体的气体检测器

Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20100138163A1 (en) * 2008-10-28 2010-06-03 Caterina Gallippi Methods, systems, and computer readable media for monitored application of mechanical force to samples using acoustic energy and mechanical parameter value extraction using mechanical response models

Family Cites Families (13)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH08191834A (ja) * 1995-01-17 1996-07-30 Hitachi Medical Corp 超音波計測装置
US5810731A (en) 1995-11-13 1998-09-22 Artann Laboratories Method and apparatus for elasticity imaging using remotely induced shear wave
FR2791136B1 (fr) 1999-03-15 2001-06-08 Mathias Fink Procede et dispositif d'imagerie utilisant les ondes de cisaillement
JP2003530941A (ja) * 2000-04-26 2003-10-21 コーニンクレッカ フィリップス エレクトロニクス エヌ ヴィ 剪断波パラメター評価のための超音波方法およびシステム
JP2003010183A (ja) * 2001-07-02 2003-01-14 Matsushita Electric Ind Co Ltd 超音波診断装置
JP4260523B2 (ja) * 2002-04-25 2009-04-30 親良 炭 変位計測装置、歪計測装置、弾性率・粘弾性率計測装置、及び、治療装置
FR2844058B1 (fr) 2002-09-02 2004-11-12 Centre Nat Rech Scient Procede et dispositif d'imagerie utilisant des ondes de cisaillement
US7744535B2 (en) * 2004-07-30 2010-06-29 Wisconsin Alumni Research Foundation Method and apparatus for acoustoelastic extraction of strain and material properties
US8016757B2 (en) * 2005-09-30 2011-09-13 University Of Washington Non-invasive temperature estimation technique for HIFU therapy monitoring using backscattered ultrasound
FR2899336B1 (fr) * 2006-03-29 2008-07-04 Super Sonic Imagine Procede et dispositif pour l'imagerie d'un milieu viscoelastique
JP2010526626A (ja) * 2007-05-16 2010-08-05 スーパー ソニック イマジン 関心領域の粘弾性の平均値を測定するための方法および装置
US8187187B2 (en) 2008-07-16 2012-05-29 Siemens Medical Solutions Usa, Inc. Shear wave imaging
US20100191113A1 (en) 2009-01-28 2010-07-29 General Electric Company Systems and methods for ultrasound imaging with reduced thermal dose

Patent Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20100138163A1 (en) * 2008-10-28 2010-06-03 Caterina Gallippi Methods, systems, and computer readable media for monitored application of mechanical force to samples using acoustic energy and mechanical parameter value extraction using mechanical response models

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
BOUCHARD R ET AL: "Image quality, tissue heating, and frame rate trade-offs in acoustic radiation force impulse imaging", IEEE TRANSACTIONS ON ULTRASONICS, FERROELECTRICS AND FREQUENCY CONTROL, IEEE, US, vol. 56, no. 1, 1 January 2009 (2009-01-01), pages 63 - 76, XP011267404, ISSN: 0885-3010, DOI: 10.1109/TUFFC.2009.1006 *
JEREMY J DAHL ET AL: "A parallel tracking method for acoustic radiation force impulse imaging", IEEE TRANSACTIONS ON ULTRASONICS, FERROELECTRICS AND FREQUENCY CONTROL, IEEE, US, vol. 54, no. 2, 1 February 2007 (2007-02-01), pages 301 - 312, XP011168520, ISSN: 0885-3010, DOI: 10.1109/TUFFC.2007.244 *

Also Published As

Publication number Publication date
JP2015154949A (ja) 2015-08-27
DE102011114943B4 (de) 2023-04-20
JP5773833B2 (ja) 2015-09-02
KR101586998B1 (ko) 2016-01-20
JP6053860B2 (ja) 2016-12-27
US20120089019A1 (en) 2012-04-12
CN102551801A (zh) 2012-07-11
JP2012081269A (ja) 2012-04-26
US8961418B2 (en) 2015-02-24
KR20120035901A (ko) 2012-04-16
FR2965934B1 (fr) 2021-10-29
DE102011114943A1 (de) 2012-04-12
KR20150035834A (ko) 2015-04-07
KR101868381B1 (ko) 2018-06-18

Similar Documents

Publication Publication Date Title
FR2965934A1 (fr) Support de memorisation, procede et systeme de resolution d'onde transversale en imagerie medicale par ultrason
US11678861B2 (en) Ultrasound imaging system memory architecture
US20220071601A1 (en) Systems and methods for improving ultrasound image quality by applying weighting factors
US20230277158A1 (en) Ultrasound imaging system memory architecture
FR3005563B1 (fr) Mesures clairsemees de proprietes du tissu dans l'imagerie medicale par ultrason
FR3031448A1 (fr)
FR3039981A1 (fr)
FR3034975A1 (fr)
FR3003153A1 (fr) Imagerie de deplacement par arfi ultrasonore utilisant une instance temporelle adaptative
JP2012081269A5 (fr)
FR3003154A1 (fr) Estimation de la fraction de matieres grasses en utilisant des ultrasons partir d'une propagation d'onde de cisaillement
CN106529561A (zh) 超声彩色流中的闪光伪像检测
FR2986701A1 (fr) Caracterisation d'onde de cisaillement sur l'axe avec un ultrason
FR3047405A1 (fr)
FR3079406A1 (fr) Balayage de fréquence pour une impulsion à force de rayonnement acoustique
FR3049845A1 (fr)
FR3085760A1 (fr) Angles pour une imagerie par onde de cisaillement à base d’ultrasons
FR3079059A1 (fr) Filtrage adaptatif de fouillis dans une imagerie par ultrasons basée sur une force de rayonnement acoustique
FR3046693A1 (fr)
FR2986960A1 (fr) Procede et systeme de visualisation d'information associee dans une imagerie par onde de cisaillement ultrasonore ainsi que support de stockage lisible par ordinateur
FR3062049A1 (fr) Imagerie par vitesse de cisaillement utilisant une coherence
EP4189380B1 (fr) Procédé de détection de discontinuités et système mettant en oeuvre ce procédé
FR3065811A1 (fr) Foyer variable pour imagerie par ondes de cisaillement
FR3081098A1 (fr) Imagerie à onde de cisaillement reposant sur les ultrasons à intervalle accru de répétition des impulsions

Legal Events

Date Code Title Description
PLFP Fee payment

Year of fee payment: 5

PLFP Fee payment

Year of fee payment: 6

PLFP Fee payment

Year of fee payment: 7

PLFP Fee payment

Year of fee payment: 8

PLFP Fee payment

Year of fee payment: 9

PLFP Fee payment

Year of fee payment: 10

PLFP Fee payment

Year of fee payment: 11

PLFP Fee payment

Year of fee payment: 12

PLFP Fee payment

Year of fee payment: 13