DE102012101312A1 - Viskoelastizitätsmessung unter Verwendung einer amplitudenphasenmodulierten Ultraschallwelle - Google Patents

Viskoelastizitätsmessung unter Verwendung einer amplitudenphasenmodulierten Ultraschallwelle Download PDF

Info

Publication number
DE102012101312A1
DE102012101312A1 DE102012101312A DE102012101312A DE102012101312A1 DE 102012101312 A1 DE102012101312 A1 DE 102012101312A1 DE 102012101312 A DE102012101312 A DE 102012101312A DE 102012101312 A DE102012101312 A DE 102012101312A DE 102012101312 A1 DE102012101312 A1 DE 102012101312A1
Authority
DE
Germany
Prior art keywords
tissue
displacement
amplitude
waveform
function
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
DE102012101312A
Other languages
English (en)
Other versions
DE102012101312B4 (de
Inventor
Caroline Maleke
Liexiang Fan
Kevin Michael Sekins
Roee Lazebnik
John Benson
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 DE102012101312A1 publication Critical patent/DE102012101312A1/de
Application granted granted Critical
Publication of DE102012101312B4 publication Critical patent/DE102012101312B4/de
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N33/00Investigating or analysing materials by specific methods not covered by groups G01N1/00 - G01N31/00
    • G01N33/48Biological material, e.g. blood, urine; Haemocytometers
    • G01N33/483Physical analysis of biological material
    • 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
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/05Detecting, measuring or recording for diagnosis by means of electric currents or magnetic fields; Measuring using microwaves or radio waves 
    • A61B5/055Detecting, measuring or recording for diagnosis by means of electric currents or magnetic fields; Measuring using microwaves or radio waves  involving electronic [EMR] or nuclear [NMR] magnetic resonance, e.g. magnetic resonance imaging
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B8/00Diagnosis using ultrasonic, sonic or infrasonic waves
    • A61B8/08Detecting organic movements or changes, e.g. tumours, cysts, swellings
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B8/00Diagnosis using ultrasonic, sonic or infrasonic waves
    • A61B8/52Devices using data or image processing specially adapted for diagnosis using ultrasonic, sonic or infrasonic waves
    • A61B8/5215Devices using data or image processing specially adapted for diagnosis using ultrasonic, sonic or infrasonic waves involving processing of medical diagnostic data
    • A61B8/5238Devices using data or image processing specially adapted for diagnosis using ultrasonic, sonic or infrasonic waves involving processing of medical diagnostic data for combining image data of patient, e.g. merging several images from different acquisition modes into one image
    • A61B8/5261Devices using data or image processing specially adapted for diagnosis using ultrasonic, sonic or infrasonic waves involving processing of medical diagnostic data for combining image data of patient, e.g. merging several images from different acquisition modes into one image combining images from different diagnostic modalities, e.g. ultrasound and X-ray
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N29/00Investigating or analysing materials by the use of ultrasonic, sonic or infrasonic waves; Visualisation of the interior of objects by transmitting ultrasonic or sonic waves through the object

Landscapes

  • Health & Medical Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Physics & Mathematics (AREA)
  • Engineering & Computer Science (AREA)
  • General Health & Medical Sciences (AREA)
  • Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
  • Pathology (AREA)
  • Biomedical Technology (AREA)
  • Molecular Biology (AREA)
  • Biophysics (AREA)
  • Animal Behavior & Ethology (AREA)
  • Heart & Thoracic Surgery (AREA)
  • Surgery (AREA)
  • Radiology & Medical Imaging (AREA)
  • Medical Informatics (AREA)
  • Public Health (AREA)
  • Veterinary Medicine (AREA)
  • Chemical & Material Sciences (AREA)
  • General Physics & Mathematics (AREA)
  • Analytical Chemistry (AREA)
  • Biochemistry (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Immunology (AREA)
  • High Energy & Nuclear Physics (AREA)
  • Hematology (AREA)
  • Urology & Nephrology (AREA)
  • Food Science & Technology (AREA)
  • Medicinal Chemistry (AREA)
  • Ultra Sonic Daignosis Equipment (AREA)
  • Magnetic Resonance Imaging Apparatus (AREA)
  • Investigating Or Analyzing Materials By The Use Of Ultrasonic Waves (AREA)

Abstract

Eine viskoelastische Gewebeeigenschaft wird in vivo gemessen (32). Zur Gewinnung von weiteren Informationen und/oder zum Schätzen der Viskosität, des Schubmoduls, und/oder anderer Schubmerkmale wird eine amplituden- und phasenmodulierte Wellenform zu dem Gewebe gesendet (30). Die durch die Wellenform im Zeitverlauf verursachte Verschiebung umfasst Verschiebungen in Zusammenhang mit der Antwort auf verschiedene Frequenzen. Durch Untersuchung der Verschiebung in dem Frequenzbereich können eine oder mehr viskoelastische Eigenschaften für verschiedene Frequenzen berechnet werden (40). Der Frequenzgang kann die Gesundheit des Gewebes anzeigen.

Description

  • HINTERGRUND
  • Die vorliegenden Ausführungsformen beziehen sich auf Viskoelastizitätsmessungen. Insbesondere wird eine viskoelastische Gewebeeigenschaft in vivo gemessen.
  • Die Schubmodulmerkmale einer Gewebeprobe können mit einem Rheometer gemessen werden. Ein Rheometer ist eine mechanische Vorrichtung, die in der Lage ist, eine Probe entweder einer dynamischen (z. B. sinusförmigen) oder einer statischen (z. B. linearen) Verformung auszusetzen. Das Rheometer kann den Speicherschubmodul und den Verlustschubmodul für einen Frequenzbereich messen. Der Bereich ist bedingt durch die Materialsteifigkeit begrenzt, wie etwa 1 bis 10 Hz für weiches Gewebe. Für die medizinische Anwendung wird das Gewebe einem Patienten entnommen und in das Rheometer eingelegt. Das Rheometer wird nicht für In-vivo-Messungen verwendet. Die Messungen richten sich nach der Größe und Gestalt der entnommenen Gewebeprobe sowie nach den Randbedingungen, die durch die Entnahme bedingt sind.
  • Schubmerkmale können in vivo unter Verwendung von Ultraschall gemessen werden. Zum Beispiel wird die Erkennung der Schubgeschwindigkeit in verschiedenen diagnostischen Anwendungen eingesetzt, wie etwa bei der Untersuchung von Lebererkrankungen. Zur Erkennung der Schubwelle wird ein Ultraschallschiebeimpuls (z. B. ein Einzyklus-Impuls) entlang einer Abtastlinie gesendet. Der Schiebeimpuls erzeugt eine Schubwelle, was zur Verschiebung des Gewebes führt. Die Verschiebung wird gemessen. Zur Erkennung der Schubwellengeschwindigkeit werden mehrere Schiebeimpulse entlang der gleichen Abtastlinie und entsprechende Abtastungen zur Verschiebungserkennung verwendet. Diese Schubmessungen können begrenzte Informationen bieten oder vom Kompressionsgrad abhängen.
  • KURZE ZUSAMMENFASSUNG
  • Als eine Einleitung umfassen die unten beschriebenen bevorzugten Ausführungsformen Verfahren, Anweisungen und Systeme zur In-vivo-Messung einer viskoelastischen Eigenschaft von Weichteilen. Zur Erhebung von weiteren Informationen und/oder zur Abschätzung der Viskosität, des Schubmoduls und/oder anderer Schubmerkmale wird eine amplituden- und phasenmodulierte Wellenform zu dem Gewebe gesendet. Die durch die Wellenform im Zeitverlauf verursachte Verschiebung umfasst Verschiebungen, die in Zusammenhang mit der Antwort auf verschiedenen Frequenzen stehen. Durch die Untersuchung der Verschiebung in dem Frequenzbereich können eine oder mehrere viskoelastische Eigenschaften für verschiedene Frequenzen berechnet werden. Der Frequenzgang kann die Gesundheit des Gewebes anzeigen.
  • In einem ersten Aspekt sind in einem nichtflüchtigen, computerlesbaren Speichermedium Daten gespeichert, die Anweisungen darstellen, die von einem programmierten Prozessor zur In-vivo-Messung einer viskoelastischen Gewebeeigenschaft ausführbar sind. Das Speichermedium umfasst Anweisungen zum Senden einer amplitudenmodulierten und phasenmodulierten Wellenform, einschließlich Zyklen in verschiedenen Frequenzen zu verschiedenen Zeitpunkten, wobei die verschiedenen Frequenzen in einem Bereich von Frequenzen liegen und die amplitudenmodulierte und phasenmodulierte Wellenform zu dem Gewebe in einem Patienten gesendet wird, zum Berechnen der Gewebeverschiebung als Funktion der Zeit als Antwort auf die amplitudenmodulierte und phasenmodulierte Wellenform, wobei die Verschiebung aus Abtastungen des Gewebes berechnet wird, zum Anwenden einer Fourier-Transformation auf die Verschiebung im Zeitverlauf und zum Ermitteln der viskoelastischen Eigenschaft aus einer Fourier-Transformation einer Schubwellengleichung und der Fourier-Transformation der Verschiebung im Zeitverlauf.
  • In einem zweiten Aspekt wird ein Verfahren zur viskoelastischen Messung unter Verwendung von Ultraschall bereitgestellt. Eine amplituden- und phasenmodulierte Ultraschallwellenform wird in vivo in einen Patienten gesendet. Aus dem Schub, der durch die amplituden- und phasenmodulierte Ultraschallwellenform verursacht wird, wird eine viskoelastische Eigenschaft gemessen. Die viskoelastische Eigenschaft wird über einen Bereich von Frequenzen gemessen, der in Zusammenhang mit der amplituden- und phasenmodulierten Ultraschallwellenform steht und unabhängig von einem Grad an Kompression ist, die während des Sendens von außen auf den Patienten angewendet wird.
  • In einem dritten Aspekt wird ein System zur Berechnung einer Schubwelle unter Verwendung von Ultraschall bereitgestellt. Ein Sendestrahlformer ist betrieblich einsetzbar, um eine amplituden- und phasenmodulierte Wellenform zu erzeugen. Ein Ultraschall-Schallkopf ist so mit dem Sendestrahlformer verbunden, dass der Ultraschall-Schallkopf als Antwort auf die amplituden- und phasenmodulierte Wellenform akustische Energie zu dem Gewebe in einem Patienten aussendet. Ein Empfangsstrahlformer ist betrieblich einsetzbar, um Daten auszugeben, die räumliche Orte als Funktion der empfangenen akustischen Signale darstellen. Ein Prozessor ist dazu konfiguriert, die Verschiebung des Gewebes im Zeitverlauf als Funktion der ausgegebenen Daten abzuschätzen und Schubinformationen als Funktion der Verschiebung des Gewebes im Zeitverlauf zu berechnen.
  • Eine Anzeige ist betrieblich einsetzbar, um ein Bild anzuzeigen, das eine Funktion der Schubinformationen ist.
  • Die vorliegende Erfindung wird durch die folgenden Ansprüche definiert, und nichts in diesem Abschnitt darf als Einschränkung dieser Ansprüche aufgefasst werden. Weitere Aspekte und Vorteile der Erfindung werden nachstehend in Verbindung mit den bevorzugten Ausführungsformen beschrieben und können später unabhängig voneinander oder in Kombination miteinander beansprucht werden.
  • KURZE BESCHREIBUNG DER ZEICHNUNGEN
  • Die Komponenten und Figuren sind nicht notwendigerweise maßstabsgerecht, stattdessen wurde das Gewicht auf die Erklärung der Grundgedanken der Erfindung gelegt. Außerdem bezeichnen in den Figuren gleiche Referenznummern übereinstimmende Teile in den verschiedenen Ansichten.
  • 1 ist ein Flussdiagramm einer Ausführungsform eines Verfahrens zur In-vivo-Messung einer viskoelastischen Gewebeeigenschaft;
  • 2 ist eine graphische Darstellung einer schubauslösenden Ultraschallsendung;
  • 3 veranschaulicht eine beispielhafte amplituden- und phasenmodulierte Wellenform;
  • 4 zeigt eine Ausführungsform von Abtastungen für die Erkennung der Verschiebung bezüglich der amplituden- und phasenmodulierten Wellenform von 3;
  • 5 veranschaulicht die relativen Orte eines Schallkopfs zum Auslösen einer Schubwellenfortpflanzung innerhalb des Zielgewebes und Hochfrequenzspulen zur Messung der Verschiebung in einer beispielhaften Magnetresonanz-Ausführungsform; und
  • 6 ist eine Ausführungsform eines Systems zur Schubwellenberechnung unter Verwendung von Ultraschall.
  • DETAILLIERTE BESCHREIBUNG DER ZEICHNUNGEN UND GEGENWÄRTIG BEVORZUGTE AUSFÜHRUNGSFORMEN
  • Mechanische Eigenschaften von Gewebe werden in vivo gemessen. Die Eigenschaften eines Gebildes oder Materials in einem Patienten, wie etwa Weichteile, können gemessen werden. Eine amplitudenmodulierte (AM) und phasenmodulierte (PM) Ultraschallwellenform verursacht eine Verschiebung des Gewebes. Die eine AM-PM-Wellenform produziert eine Schwingung in aufeinanderfolgenden Frequenzen im Innern des Patienten an einem gewünschten Gewebe. Der Schubmodul und die Viskosität des Gewebes werden mit dieser einen Anregung bei verschiedenen Frequenzen gemessen, was eine Messung in weniger als zwei Sekunden erlaubt. Phasenmodulation bei niedrigeren Frequenzen (z. B. 1 bis 100 Hz) kann ein höheres Signal-Rausch-Verhältnis bieten. Die AM-PM-Wellenform erlaubt das Erfassen von Informationen während oder nach dem Induzieren des Schubs.
  • Ein rasches und zuverlässiges Verfahren kann die Materialsteifigkeit für verschiedene Frequenzen in einer einzigen Messung beurteilen. Durch die Erhebung von Messwerten für verschiedenen Frequenzen können die daraus hervorgehenden Erkenntnisse über die mechanischen Eigenschaften eines Materials (z. B. Gewebe) diagnostische Verfahren ergänzen. Die Erkenntnisse können quantitative Parameter umfassen, wie etwa mechanische Beanspruchung, Viskosität, und/oder Schubmodul. Die quantitativen Parameter sind unabhängig von dem Grad der auf das Gewebe angewandten Kompression, da die Parameter von der Verschiebung abgeleitet sind. Der Frequenzgang kann unabhängig von der Kompression sein. Die Erkenntnisse können qualitative Parameter umfassen, wie etwa Verschiebung. Die Verschiebung kann von einem Grad des auf das Gewebe angewendeten Drucks abhängig sein. Die Erkenntnisse können abgeleitete Informationen umfassen. Die Steigung von Schubmodul gegenüber Frequenz und/oder die Steigung von Viskosität gegenüber Frequenz kann Informationen über die mechanische Gewebeeigenschaft liefern, die in Zusammenhang mit der Gewebepathologie stehen. Die Steigungen sind unabhängig von den Vorkompressionsgraden. Die Messungen können wiederholt werden, um die Genauigkeit zu steigern. Da die Messungen nicht auf eine Kontaktherstellung angewiesen sind, sind die Messungen unabhängig von Randbedingung, Materialgröße und Gestalt.
  • 1 zeigt ein Verfahren für eine viskoelastische Messung unter Verwendung von Ultraschall. Eine Ultraschallstrahlung wird verwendet, um eine Gewebeverschiebung als Antwort auf verschiedene Frequenzen zu erzeugen. Auf diese Weise wird Ultraschall verwendet, um eine viskoelastische Eigenschaft zu berechnen. Das Verfahren wird durch das in 6 gezeigte System oder ein anderes System implementiert. Es können zusätzliche, verschiedene oder weniger Vorgänge bereitgestellt werden. Zum Beispiel wird das Verfahren ohne Vorgänge 38, 42 und/oder 44 durchgeführt. Als weiteres Beispiel wird Vorgang 32 ohne einen oder mehr der Vorgänge 34 bis 40 durchgeführt. Die Vorgänge werden in der beschriebenen oder gezeigten Reihenfolge durchgeführt, können aber in anderen Reihenfolgen durchgeführt werden.
  • In Vorgang 30 wird eine amplitudenmodulierte und phasenmodulierte Wellenform gesendet. Ein Ultraschall-Schallkopf sendet eine akustische Wellenform aus, die von einer elektrischen Wellenform umgewandelt wurde. Die akustische Energie wird mit der amplitudenmodulierten und phasenmodulierten Wellenform zu dem Gewebe in einem Patienten gesendet. Die Sendung erfolgt in vivo.
  • Die akustische Wellenform wird gesendet, um eine Schubwelle zu erzeugen. Die Anregung ist ein Ultraschall-Schiebeimpuls. Die akustische Energie wird gebündelt, wodurch sich ein oder mehr Strahlen für jede Sendung ergeben. Die Anregung wird unter Verwendung einer phasengesteuerten Anordnung und/oder eines mechanischen Brennpunkts gebündelt. Die Anregung wird an einem Ort gebündelt, um die Erkennung der sich ergebenden Schubwelle zu ermöglichen, wie zum Beispiel an einem Gewebeort gebündelt, die einen möglichen Tumor umgibt und/oder einen möglichen Tumor umfasst.
  • Wie in 2 dargestellt, wird die Schubwelle 52 in der Brennregion 50 von dem Schallkopf 54 erzeugt und pflanzt sich von der Brennregion 50 ausgehend in lateraler Richtung weiter fort. Die Pfeile sind in einer Richtung gezeigt (z. B. horizontal), aber die Schubwelle bewegt sich in mehrere verschiedene Richtungen. Die Amplitude der Schubwelle vermindert sich, während die Welle sich durch das Gewebe bewegt.
  • Zum Erzeugen der Schubwelle sind Anregungen mit hoher Amplitude oder von hoher Kraft erwünscht. Zum Beispiel hat die Anregung einen mechanischen Index von annähernd 1,9, der aber 1,9 nicht überschreitet. Als zurückhaltender Ansatz und um der Sondenvariation Rechnung zu tragen, kann ein mechanischer Index von 1,7 oder ein anderer Wert als Obergrenze verwendet werden. Es kann mehr Kraft (z. B. ein MI von über 1,9) oder weniger Kraft verwendet werden. Das Senden entlang der gleichen oder benachbarten Abtastlinien kann dazu führen, dass die Temperatur des Gewebes im Zeitverlauf ansteigt. Biologische Wirkungen können u. a. sein: Hyperthermie bei einer Gewebetemperatur von etwa 41 bis 45°C, Proteindenaturierung bei Temperaturen über 43 bis 45°C und Gewebenekrose bei Temperaturen über 50°C. Die Gewebesteifigkeit kann sogar bei Temperaturen unter 43 bis 45°C beeinflusst werden. Bei Temperaturen über 43 bis 45°C kann die Viskosität und/oder die Steifigkeit erhöht sein. Bei Temperaturen über 50°C kann das Gewebe eine höhere Steifigkeit und/oder eine hohe Dämpfung aufweisen. Biologische Wirkungen werden begrenzt, indem ein Temperaturanstieg von mehr als 2 Grad Celsius verhindert wird. Ansonsten können die Sendungen biologische Wirkungen auslösen.
  • Die elektrischen und entsprechenden akustischen Wellenformen umfassen Zyklen in verschiedenen Frequenzen zu verschiedenen Zeitpunkten. Die verschiedenen Frequenzen liegen in einem Bereich von Frequenzen. Der Frequenzbereich kann auf der Grundlage des untersuchten Gewebetyps optimiert werden. So weist zum Beispiel Brustgewebe unter Umständen eine maximale Verschiebungsantwort von 50 Hz auf, und daher umfasst der Bereich 50 Hz. Als weiteres Beispiel sind Tumoren unter Umständen härter als Weichteile, und somit haben sie eine größere Frequenz, bei der die maximale Verschiebungsantwort auftritt. Bei der Tumormessung kann der Frequenzbereich größer sein. In einer Ausführungsform liegt der Frequenzbereich für die Phasenmodulation zwischen 15 und 75 Hz.
  • Für die Phasenmodulation kann eine beliebige Frequenzvariationsfunktion verwendet werden. Zum Beispiel wird eine Chirp-Frequenzmodulation innerhalb des Bereichs verwendet. 3 zeigt eine Wellenform mit einer amplituden- und phasenmodulierten Hüllkurve. Die Phasenmodulation beginnt bei der niedrigsten Frequenz und wird moduliert oder allmählich verändert bis hin zur höchsten Frequenz in dem Bereich. Es kann auch eine andere Frequenzvariation verwendet werden, wie etwa nicht-linear, zufällig, oder andere Schritte zwischen den verschiedenen Frequenzen in beliebiger Reihenfolge.
  • Die Amplitudenmodulation ist eine beliebige Funktion. Die Amplitudenmodulation trennt einzelne Frequenzen oder Gruppen von Frequenzen voneinander.
  • In einer Ausführungsform wird die amplituden- und phasenmodulierte Wellenform erzeugt, indem ein sinusförmiger Träger (im MHz-Bereich) mit einer sinusförmigen Amplitudenmodulation mit einem phasenvariierenden Term (im Hz-Bereich) multipliziert wird. Die Amplitudenmodulation mit dem phasenvariierenden Term definiert eine Hüllkurve mit Amplituden- und Phasenmodulation. Die Wellenform kann dargestellt werden als: x(t) = Asin(ωct) × sin((ωm + Δωmt)t) wobei x(t) die Wellenform als Funktion der Zeit ist, A ein Amplitudengewicht ist, ωc eine Mittenfrequenz oder Trägerfrequenz ist und Δωm eine Modulationsfrequenz ist, die sich im Zeitverlauf ändert. In einer Ausführungsform variiert ωm + Δωm zwischen 15 und 75 Hz, aber es können auch andere Bereiche verwendet werden. 3 zeigt eine Wellenform mit einem Bereich zwischen 15 und 75 Hz.
  • Die Wellenform weist eine beliebige Zahl von Hüllkurven- und Trägerzyklen auf. In einer Ausführungsform weist die Wellenform genügend Zyklen auf, um 1333 Millisekunden lang zu sein, aber es kann auch eine längere oder kürzere Zeit verwendet werden. Die Zahl der Trägerzyklen beträgt Zehner, Hunderte, Tausende oder mehr. Die Zahl der Hüllkurvenzyklen oder Modulationszyklen beträgt zwei oder mehr.
  • Die Gewebeantwort ist von der Wellenform, x(t), und von den Gewebemerkmalen abhängig. Die Verschiebung, y(t), des Gewebes im Zeitverlauf kann als eine Faltung der Wellenform, x(t), und der Gewebemerkmale oder Gewebeantwort, h(t), ausgedrückt werden: y(t) = x(t)·h(t). Die Gewebeantwort spiegelt die viskoelastischen Eigenschaften des Gewebes wider.
  • In Vorgang 32 werden eine oder mehrere viskoelastische Eigenschaften aus dem durch die amplituden- und phasenmodulierte Ultraschallwellenform verursachten Schub gemessen. Um die viskoelastischen Eigenschaften in Vorgang 32 zu messen, werden die Verschiebung, y(t), des Gewebes im Zeitverlauf als Antwort auf die amplituden- und phasenmodulierte Wellenform in Vorgang 34 gemessen. Vorgang 32 ist so gezeigt, dass er Vorgänge 34, 36, 38 und 40 einschließt. Es können verschiedene, zusätzliche oder weniger Vorgänge bereitgestellt werden. In anderen Ausführungsformen werden andere Prozesse verwendet, um die viskoelastische Eigenschaft zu messen.
  • In Vorgang 34 wird die Verschiebung als Funktion der Zeit berechnet. Das Gewebe wird mehrere Male abgetastet, um die Verschiebung zu ermitteln, wie etwa ein mindestens dreimaliges Abtasten einer Region zu zwei verschiedenen Zeitpunkten, um die Verschiebung zu ermitteln.
  • Die Verschiebung des Gewebes erfolgt als Antwort auf die amplitudenmodulierte und phasenmodulierte Wellenform. Die durch die Schubwelle verursachte Verschiebung des Gewebes wird im Zeitverlauf ermittelt. Wenn die Schubwelle eine bestimmte Stelle durchläuft, wird das Gewebe um einen Grad oder eine Entfernung verschoben, die sich zu einem Höchstgrad steigert und dann wieder abnimmt, wenn das Gewebe wieder in die Ruheposition zurückkehrt. Da die AM-PM-Wellenform viele Zyklen umfasst, kann das Gewebe fortlaufend verschoben werden. Die Amplitudenmodulation verursacht eine Varianz in der im Zeitverlauf eintretenden Verschiebung, während das Gewebe beginnt, wieder zu einem Normalzustand zurückzukehren, wenn sich die Amplitude verringert. Die Phasenmodulation führt zur Variation des Grads der im Zeitverlauf eintretenden Verschiebung.
  • Das Gewebe wird unter Verwendung einer beliebigen Bildgebungsmodalität abgetastet, die in der Lage ist, während der Reaktion des Gewebes auf die schiebende Wellenform eine Abtastung auf Verschiebung durchzuführen, wie etwa während oder nach der Anwendung der AM-PM-Wellenform. Die Abtastung findet statt, bevor das Gewebe wieder zu einem entspannten oder normalen Zustand oder zu einer entspannten oder normalen Position zurückkehrt, kann aber Abtastungen des Gewebes zu Zeitpunkten umfassen, zu denen das Gewebe wieder in der Ruheposition ist. Ultraschall und Magnetresonanztomographie sind zwei mögliche Modalitäten zum Berechnen der Verschiebung aus Abtastungen des Gewebes.
  • Für die Ultraschallabtastung wird die Schubwelle 52 an Orten erkannt, die an die Brennregion für die AM-PM-Wellenform angrenzen und/oder in einem Abstand von dieser angeordnet sind. Zum Erkennen der Verschiebung wird Ultraschallenergie zu dem Gewebe gesendet, in dem eine Verschiebung stattfindet, und dann werden die Echos der Energie empfangen. Zum Erkennen der Gewebeantwort auf Schubwellen in einem Untersuchungsbereich erfolgen Sendungen an andere Brennregionen, und die Erkennung wird um die anderen Brennregionen herum durchgeführt. Diese anderen Sendungen dienen zum Erkennen der Schubwellen, nicht zum Auslösen der Schubwelle. Die Sendungen zur Erkennung können weniger Kraft und/oder kürzere Impulse aufweisen (z. B. 1 bis 5 Trägerzyklen) und die gleiche oder andere Abtastlinie verwenden wie die AM-PM-Wellenform. Die Sendungen zum Erkennen können ein breiteres Strahlenprofil entlang mindestens einer Dimension, wie etwa lateral, aufweisen, um gleichzeitig Empfangsabtastpunkte entlang einer Vielzahl von Abtastlinien bilden zu können. Die Schubwelle kann in einer, zwei oder mehr Richtungen überwacht werden.
  • Ein Untersuchungsbereich wird überwacht, um die Schubwelle zu erkennen. Der Untersuchungsbereich hat eine beliebige Größe, wie etwa 6 mm in lateraler und 10 mm in axialer Richtung um den Brennpunkt der AM-PM-Wellenform. Diese Erkennungsregion wird mittels Ultraschall überwacht. Zum Beispiel werden B-Mode-Abtastungen durchgeführt, um die durch die Schubwelle verursachte Gewebeverschiebung zu erkennen. Zur Schubwellenüberwachung kann Doppler-, farbkodierter Strömungs- oder ein anderer Ultraschallmodus verwendet werden.
  • Die Überwachung erfolgt für eine beliebige Zahl von Abtastlinien. Zum Beispiel werden vier Empfangsstrahlen als Antwort auf jede Sendung geformt. Nach dem Senden der Anregung zum Erzeugen der Schubwelle werden wiederholte B-Mode-Sendungen entlang einer oder mehrerer der Sende-Abtastlinien sowie Empfangsvorgänge entlang der entsprechenden Empfangs-Abtastlinien durchgeführt. In anderen Ausführungsformen wird nur ein einzelner Empfangsstrahl geformt oder eine andere Anzahl von Empfangsstrahlen werden als Antwort auf jede Sendung geformt. Manche der Ultraschalldaten, wie etwa am Beginn oder am Ende der Wiederholungen, reagieren unter Umständen nicht auf die Schubwelle.
  • Das Senden und Empfangen zur Erkennung erfolgt mehrere Male, um eine Veränderung zu ermitteln, die sich aufgrund der Verschiebung im Zeitverlauf ergibt. Es kann eine beliebige Sende- und Empfangssequenz verwendet werden. Die Erkennung der Verschiebung kann mit anderen Abtastungen verschachtelt sein, wie etwa separaten Abtastungen unterschiedlicher Regionen auf Verschiebung.
  • Das Senden und Empfangen zur Erkennung sind mit der amplituden- und phasenmodulierten Wellenform verschachtelt. Zum Beispiel werden die Sendung(en) und der Empfang (die Empfangsvorgänge) zum einmaligen Abtasten einer Region jedes Mal durchgeführt, wenn die Amplitudenmodulation Null oder nahezu Null erreicht. Als ein anderes Beispiel ist das Abtasten des Gewebes mit dem regelmäßigen Senden von Vorgang 30 verschachtelt, wie etwa alle 1,67 Millisekunden. Die Verschachtelung vermeidet das Auftreten von Interferenzen, wie etwa ein Anhalten der Amplituden- und Phasenmodulations-Wellenform während der Abtastung, sodass Echos von der Wellenform auf ein Mindestmaß verringert werden. 4 zeigt ein Beispiel für die Verschachtelung von Abtastungen der Ebenen 56 mit der amplituden- und phasenmodulierten Wellenform. Der Verschachtelungszeitraum kann mehr oder weniger häufig sein. Wenn die Amplituden- und Phasenmodulation nach dem Anhalten für die Verschiebungsabtastung wieder fortgesetzt wird, beginnt die Wellenform von einem Punkt in der Wellenform, an dem das Senden beendet wurde. In alternativen Ausführungsformen erfolgt die Abtastung auf einer anderen Frequenz oder mit einer anderen Kodierung als die amplituden- und phasenmodulierte Wellenform. Sowohl das Schieben der Wellenform als auch die Abtastung können gleichzeitig erfolgen, und die Frequenz oder die Kodierung wird verwendet, um die Echos voneinander zu unterscheiden.
  • Die obige Beschreibung betrifft die Tiefe oder den Ort. Die viskoelastische Eigenschaft wird für einen Ort gemessen. Um eine größere Region zu messen, werden die Vorgänge 30 bis 40 für andere Orte wiederholt. Für jeden Empfangsstrahlort wird ein Zeitprofil der Bewegungsinformationen bereitgestellt, das durch die Ultraschalldaten dargestellt wird. Ein separates Zeitprofil wird für jede axiale Tiefe sowie für jeden lateralen Ort bereitgestellt.
  • Die Verschiebung wird aus den Ultraschall-Abtastdaten berechnet. Das Gewebe bewegt sich zwischen zwei Abtastungen. Die Daten einer Abtastung werden bezüglich der Daten in der anderen Abtastung in eine, zwei oder drei Dimensionen übersetzt. Für jede mögliche relative Position wird ein Grad der Ähnlichkeit berechnet. Der Ähnlichkeitsgrad wird mit Korrelation ermittelt, wie etwa einer Kreuzkorrelation. Eine Mindestsumme absoluter Differenzen oder eine andere Funktion kann verwendet werden. Die räumliche Versetzung mit der höchsten oder ausreichenden Korrelation weist auf Grad und Richtung der Verschiebung hin.
  • Das Ermitteln der Verschiebungen erfolgt für einen bestimmten Ort zu verschiedenen Zeitpunkten und somit in Verbindung mit sequenziellen Abtastungen. Die Verschiebung wird in Bezug auf einen anfänglichen oder Bezugs-Frame von Abtastdaten ermittelt (d. h. kumulative Verschiebung). Eine andere Möglichkeit besteht darin, dass die Verschiebung aus dem unmittelbar vorhergegangenen Frame von Abtastdaten ermittelt wird, sodass der vorhergegangene Frame jeweils auf laufender Grundlage zum Bezugspunkt wird (d. h. inkrementelle Verschiebung). Das Zeitprofil für einen bestimmten Ort zeigt die durch die Schubwelle im Zeitverlauf und als Antwort auf verschiedene Anteile der amplituden- und phasenmodulierten Wellenform verursachten Verschiebung an.
  • Für die Abtastung mittels Magnetresonanz wird ein ähnlicher Prozess verwendet. 5 zeigt den Schallkopf 54 beim Senden an eine Brennregion, die zumindest teilweise von einem Magnetresonanztunnel umgeben ist. Der Tunnel umfasst Hochfrequenzspulen. Es wird ein gleichförmiges Magnetfeld erzeugt. Unter Verwendung eines Puls-Hochfrequenz-Gradientenverstärkers werden die Drehungen der Moleküle in den Patienten geändert. Die Spulen dienen zur Erkennung der Änderung, wobei k-Raum-Daten erfasst werden.
  • Statt der Korrelation zwischen Frames in einer Sequenz wird ein Satz von Verschiebungsdaten im Zeitverlauf erfasst. K-Raum-Magnetresonanzdaten, die das Gewebe des Patienten darstellen, werden zu verschiedenen Zeitpunkten oder in einer Sequenz erfasst. Ein Satz (Bezugssatz) wird erfasst, ohne dass das Gewebe der amplituden- und phasenmodulierten Wellenform ausgesetzt ist. Ein weiterer Satz wird erfasst, während die amplituden- und phasenmodulierte Wellenform angewendet wird. Da es keine Interferenz zwischen dem Ultraschall der Wellenform und den k-Raum-Daten gibt, können die Abtastdaten ohne Verschachtelung erfasst werden. Eine andere Möglichkeit besteht darin, dass eine Verschachtelung verwendet wird. Die Abtastsequenz ist periodisch (z. B. alle 1,67 ms), je nach den Niederamplitudenanteilen der amplitudenmodulierten Wellenform, oder einer anderen Funktion.
  • Die k-Raum-Daten werden synchronisiert mit einem oder mehreren physiologischen Zyklen erfasst, wie etwa Herz- oder Atmungszyklus. Die Abtastungen der Bezugssequenz werden zu den selben Zeitpunkten im Zyklus erfasst wie der Verschiebungssatz. Während die k-Raum-Daten verwendet werden können, ohne zu einem Bild verarbeitet zu werden, werden in alternativen Ausführungsformen Bilddaten verwendet.
  • Die Magnetresonanz-Bezugsdaten, die das Gewebe darstellen, während es frei von Antworten auf das Senden ist, wird von den Magnetresonanzdaten subtrahiert, die die Gewebeantwort auf das Senden darstellen. Eine andere Möglichkeit besteht darin, dass die Daten, die die Gewebeantwort darstellen, von den Bezugsdaten subtrahiert werden. Die Abtastungen von bezüglich des physiologischen Zyklus ähnlichen Zeitpunkten werden subtrahiert. Die Subtraktion isoliert Differenzen. Da die Schubwelle eine Differenz verursacht, bleibt das versetzte oder verschobene Gewebe übrig, und die anderen Signale heben sich in der Subtraktion auf.
  • Nach einer Tiefpassfilterung oder einer anderen Verarbeitung werden der Grad, die Richtung oder der Grad und die Richtung der Verschiebung zwischen den sequenziellen Frames ermittelt. Ein Ort, an dem ein Punkt oder eine Region mit maximaler Intensität liegt, wird in den Daten für jeden Subtraktionsdaten-Frame festgestellt. Die Differenz von Orten zwischen den sequenziellen Frames oder zwischen einem anderen Frame und einem Bezugs-Frame wird als die Verschiebung berechnet. Die Verschiebung variiert im Zeitverlauf aufgrund der Differenzen in der amplituden- und phasenmodulierte Wellenform, die zu den verschiedenen Zeitpunkten angewendet wird.
  • In Vorgang 36 wird die Verschiebung im Zeitverlauf in den Frequenzbereich umgewandelt. Eine Fourier-Transformation wird auf die Verschiebung im Zeitverlauf angewendet. Es kann eine beliebige Transformation verwendet werden. Da die Verschiebung im Zeitverlauf auf die verschiedenen Frequenzen der amplituden- und phasenmodulierten Wellenform reagiert, stellt die Transformation in den Frequenzbereich verschiedene Antwortstufen als Funktion der Frequenz bereit. Die Fourier-Transformation der Verschiebung, y(t), kann als Y(Ω) dargestellt werden, wobei Ω die Frequenz ist, d. h. die Darstellung der Verschiebung im Frequenzbereich.
  • Die Transformation der Verschiebung oder Gewebeantwort wird dargestellt als:
    Figure 00160001
  • In Vorgang 38 wird die Verschiebung im Frequenzbereich gefiltert. Es kann eine beliebige Filterung verwendet werden, wie etwa Bandpassfilterung. Der Frequenzbereich des Bandpassfilters wird je nach dem Frequenzbereich der gesendeten amplituden- und phasenmodulierten Wellenform festgelegt. Die akustische Stärke ist abhängig von der Quadratzahl des Drucks. Infolgedessen beträgt die Verschiebung oder Vibration, die durch die Wellenform induziert wird, das Zweifache der Frequenz der Wellenform. Zum Beispiel umfasst die Wellenform eine Frequenzvariation der Hüllkurve zwischen 15 und 75 Hz, sodass die reagierenden Verschiebungsinformationen zwischen 30 und 150 Hz eintreten. Die Breite des Frequenzbereichs ist das Zweifache des Bereichs der Wellenform, und die niedrigen und hohen Frequenzen im Bereich betragen das Zweifache der niedrigen und hohen Frequenzen des Bereichs der Wellenform. Die Filterung isoliert die Informationen in diesem Bereich oder einem Teilbereich. Frequenzdaten außerhalb des Bereichs werden nicht verwendet, und Daten innerhalb des Bereichs werden verwendet, was zu einer Bandpassfilterung führt. Es können auch eine andere Bandpassfilterung, andere Filterung, andere Verarbeitung oder keine Datenmodifikationsprozesse verwendet werden.
  • In Vorgang 40 wird die viskoelastische Eigenschaft aus der Fourier-Transformation der Verschiebung im Zeitverlauf ermittelt. Die Schubwelle kann dargestellt werden als:
    Figure 00170001
    Es können auch anderen Darstellungen der Schubwelle verwendet werden.
  • In dem Frequenzbereich wird die Schubwellengleichung dargestellt als:
    Figure 00170002
  • Diese Gleichung stellt die Fourier-Transformation der Schubwellengleichung dar. Je nach der verwendeten Schubwellendarstellung und Fourier-Transformation können auch andere Darstellungen verwendet werden. Die Terme auf der linken Seite der Gleichung sind bekannt oder gemessen, wie etwa die sich ergebende Bewegung (Verschiebung) im Frequenzbereich, Y(Ω), die in einer Frequenz (ω) gleich dem Zweifachen der AM-PM-Frequenz, z. B. 2 × (15 bis 75 Hz) = 30 Hz bis 150 Hz, schwingt. Die rechte Seite stellt die imaginären und reellen Teile der transformierten Schubwellengleichung dar.
  • Die viskoelastische Eigenschaft wird aus der Schubwellengleichung in dem Frequenzbereich ermittelt. Ein beliebiger Modul-, Viskositäts- oder Schubwert kann geschätzt werden. Gewebemodulwerte stellen die Härte oder Steifigkeit des Gewebes dar. Zum Beispiel wird der Schubmodul des Gewebes geschätzt. In alternativen Ausführungsformen wird der Elastizitätsmodul geschätzt. In anderen Ausführungsformen werden andere Schubwerte geschätzt, und zwar entweder quantitative oder qualitative Werte.
  • In einer Ausführungsform wird die Viskosität ermittelt. Die Viskosität wird als Funktion der Ergebnisse der Anwendung der Fourier-Transformation berechnet. Der imaginäre Teil der Fourier-Transformation der Schubwellengleichung kann verwendet werden, um die Viskosität zu ermitteln. Der imaginäre Teil wird durch die Frequenz, ω, dividiert, um die Viskosität zu erhalten. Dies wird dargestellt als:
    Figure 00180001
  • In einer alternativen oder zusätzlichen Ausführungsform wird der Schubmodul ermittelt. Der Schubmodul wird als Funktion der Ergebnisse der Anwendung der Fourier-Transformation berechnet. Der reelle Teil der Fourier-Transformation der Schubwellengleichung kann verwendet werden, um den Schubmodul zu ermitteln. Der reelle Teil selbst stellt den Schubmodul dar. Dies wird dargestellt als:
    Figure 00180002
  • Die viskoelastische Eigenschaft, wie etwa die Viskosität oder der Schubmodul, kann über einen Bereich von Frequenzen (w) ermittelt werden. Zum Beispiel liegen die Verschiebungsdaten im Bereich zwischen 30 und 150 Hz. Die Viskosität wird über den gleichen Bereich ermittelt. Die Viskosität für Gruppen von Frequenzen oder für Teilbänder kann gemittelt werden, wie etwa durch Bereitstellen von Werten für jede 5 Hz. Die Eigenschaft wird nur als Antwort auf die amplitudenmodulierte und phasenmodulierte Wellenform und Abtastungen des Gewebes während der Anwendung der Wellenform ermittelt. Die sequenzielle Durchführung der Vorgänge 30 bis 40 wird nicht bereitgestellt, kann aber bereitgestellt werden. Durch die Verwendung einer Wellenform mit einem Bereich von Frequenzen kann die frequenzbezogene Antwort des Gewebes während eines kurzen Zeitraums gemessen werden, wie etwa ausreichend zum Verschachteln der Abtastungen und der Sendung der Wellenform.
  • Die Viskosität und der Schubmodul sind quantitative Eigenschaften. Die Werte sind von einem Grad der angewendeten Kompression unabhängig. Der während des Sendens auf den Patienten ausgeübte Druck von der akustischen Wellenform oder von einer externen Quelle wird nicht benötigt, um quantitative Eigenschaften zu ermitteln. In alternativen Ausführungsformen wird eine qualitative Eigenschaft verwendet, wie etwa die Verschiebung.
  • In Vorgang 42 wird die Steigung der viskoelastischen Eigenschaft als Funktion der Frequenz festgestellt. Die Eigenschaft als Funktion der Frequenz kann in einem Diagramm dargestellt werden, mindestens über einen Teil oder den gesamten Bereich der Frequenzen (z. B. 30 bis 150 Hz). Eine Linie wird an das Diagramm oder die Daten angepasst, wodurch sich eine Steigung ergibt. Zum Beispiel werden die Steigungen des Schubmoduls und der Viskosität als Funktion der Frequenz ermittelt.
  • Die Steigung kann die Gesundheit des Gewebes anzeigen oder die Diagnose unterstützen. Zum Beispiel kann gesundes Gewebe eine waagrechtere Viskositätssteigung haben als ein Tumor. Der Frequenzgang des Gewebes kann anzeigen, ob das Gewebe gesund ist oder nicht. Es können auch andere Funktionen der Eigenschaft als Funktion der Frequenz berechnet werden, wie etwa die Varianz, der Grad der Änderung oder die Krümmung.
  • In Vorgang 44 wird ein Bild erzeugt. Das Bild stellt die Eigenschaft dar. Es kann eine Zahlen- oder Textangabe der Eigenschaft angezeigt werden, entweder bei einer bestimmten Frequenz oder berechnet aus Daten über verschiedene Frequenzen. In anderen Ausführungsformen werden das Diagramm und/oder die Linienanpassung und der Steigungswert ausgegeben. Die viskoelastische Eigenschaft wird dem Anwender in dem Bild mitgeteilt.
  • Das Bild kann ein Schaubild sein, wie etwa ein Diagramm der Werte als Funktion der Frequenz.
  • Das Bild kann zusätzlich eine ein-, zwei- oder dreidimensionale Darstellung der Eigenschaft oder andere Schubinformationen als Funktion von Raum oder Ort umfassen. Zum Beispiel wird die Schubgeschwindigkeit innerhalb einer Region angezeigt. Die Schubgeschwindigkeit moduliert die Farbe für die Pixel (Bildpunkte) in einer Region in einem grauskala-modulierten B-Mode-Bild. Das Bild kann Verschiebungsinformationen, wie etwa Schub, oder Modulen (z. B. die Schubmodulen) für verschiedene Orte darstellen. Das Anzeigeraster kann sich von dem Abtastraster und/oder dem Raster, für das die Verschiebungen berechnet werden, unterscheiden. Farbe, Helligkeit, Luminanz, Farbton oder ein anderes Merkmal wird als Funktion der Schubinformationen moduliert.
  • Die Vorgänge werden für andere Abtastlinien und/oder andere Tiefen wiederholt. Zum Beispiel werden die Vorgänge 30 bis 40 erneut für jeden Ort in einer ein-, zwei- oder dreidimensionalen Region durchgeführt.
  • 6 zeigt eine Ausführungsform eines Systems 10 zur Schubwellenberechnung unter Verwendung von Ultraschall. Ultraschall erzeugt eine Schubwelle, und die Abtastdaten, die auf das Gewebe reagieren, das auf den Ultraschall reagiert, werden verwendet, um eine Eigenschaft zu ermitteln. Das System 10 implementiert das Verfahren von 1 oder andere Verfahren. Das System 10 umfasst einen Sendestrahlformer 12, einen Schallkopf 14, einen Empfangsstrahlformer 16, einen Bildprozessor 18, eine Anzeige 20 und einen Speicher 22. Es können zusätzliche, verschiedene oder weniger Komponenten bereitgestellt werden. Zum Beispiel wird eine Anwendereingabe zum manuellen oder unterstützten Kennzeichnen eines Untersuchungsbereichs bereitgestellt, für welche Schubinformationen gewonnen werden sollen. Als weiteres Beispiel wird ein zusätzlicher HIFU-Schallkopf zur Behandlung des Gewebes bereitgestellt. Das System 10 ist ein medizinisches diagnostisches Ultraschall-Bildgebungssystem.
  • In alternativen Ausführungsformen ist das System 10 ein Computerarbeitsplatz (PC), eine sogenannte Workstation, eine PACS-Einrichtung (PACS bedeutet Bildarchivierungs- oder Kommunikationssystem) oder eine andere an einem gleichen Ort befindliche oder über ein Netzwerk verteilte Anordnung zur Echtzeit- oder verzögerten Erfassungs-Bildgebung.
  • In weiteren alternativen Ausführungsformen ist das System 10 Teil eines Magnetresonanzsystems. Zum Beispiel werden der Sendestrahlformer 12 und der Schallkopf 14 zum Senden der Wellenform bereitgestellt, um die Schubwelle zu erzeugen, aber der Empfangsstrahlformer wird nicht bereitgestellt. Stattdessen werden Magnetresonanzspulen und -magnete, wie etwa in 5 gezeigt, mit dem Prozessor 18, dem Speicher 22 und der Anzeige 20 für Abtastungen bereitgestellt.
  • Der Sendestrahlformer 12 ist ein Ultraschallsender, ein Speicher, ein Impulsgeber, eine analoge Schaltung, eine digitale Schaltung oder eine Kombination derselben. Der Sendestrahlformer 12 ist betrieblich einsetzbar, um Wellenformen für eine Vielzahl von Kanälen mit verschiedenen oder relativen Amplituden, Verzögerungen und/oder Phasenverschiebung zu erzeugen. Nach dem Senden von akustischen Wellen von dem Schallkopf 14 als Antwort auf die erzeugten Wellen werden ein oder mehr Strahlen geformt. Die Wellenformen sind jeweils eine amplituden- und phasenmodulierte Wellenform, aber mit relativen Verzögerungen und Apodisierung, um die gewünschte akustische Wellenform (z. B. siehe 3) in der Brennregion bereitzustellen.
  • Zum Abtasten der Gewebeverschiebung wird eine Sequenz von Sendestrahlen erzeugt, um eine ein-, zwei- oder dreidimensionale Region abzutasten. Es können Sektor-, Vector®-, lineare oder andere Abtastformate verwendet werden. Die gleiche Region wird mehrmals abgetastet, wie etwa in 4 dargestellt. Die Abtastung durch den Sendestrahlformer 12 ist verschachtelt oder synchronisiert mit dem Senden der amplituden- und phasenmodulierten Wellenform durch den Sendestrahlformer 12. Die gleichen Elemente des Schallkopfs 14 werden sowohl für die Abtastung als auch für die Erzeugung von Schubwellen verwendet, aber es können verschiedene Elemente, Schallkopf und/oder Strahlformer verwendet werden.
  • Der Schallkopf 14 ist eine 1-, 1,25-, 1,5-, 1,75- oder 2-dimensionale Anordnung von piezoelektrischen oder kapazitiven Membranelementen. Alternativ wird ein einzelnes Element mit einem mechanischen Brennpunkt bereitgestellt. Der Schallkopf 14 umfasst eine Vielzahl von Elementen für das Umformen zwischen akustischen und elektrischen Energien. Zum Beispiel ist der Schallkopf 14 eine eindimensionale PZT-Anordnung mit etwa 64 bis 256 Elementen.
  • Der Schallkopf 14 steht in Verbindung mit dem Sendestrahlformer 12, um elektrische Wellenformen in akustische Wellenformen umzuwandeln, und steht in Verbindung mit dem Empfangsstrahlformer 16, um akustische Echos in elektrische Signale umzuwandeln. Der Schallkopf 14 sendet akustische Energie mit einer amplituden- und phasenmodulierten Wellenform aus. Die Wellenform hat ihren Brennpunkt in einem zu untersuchenden Gewebebereich oder Gewebeort im Patienten. Die akustische Wellenform wird als Antwort auf die Anwendung der elektrischen Wellenform auf die Elemente des Schallkopfs erzeugt.
  • Für die Abtastung mittels Ultraschall zur Erkennung einer Verschiebung sendet der Schallkopf akustische Energie und empfängt Echos. Die Empfangssignale werden als Antwort auf die Ultraschallenergie (Echos) erzeugt, die auf die Elemente des Schallkopfs 14 auftrifft.
  • Der Empfangsstrahlformer 16 umfasst eine Vielzahl von Kanälen mit Verstärkern, Verzögerungen und/oder Phasendrehern und einen oder mehr Summierer. Jeder Kanal steht in Verbindung mit einem oder mehr Schallkopf-Elementen. Der Empfangsstrahlformer 16 wendet relative Verzögerungen, Phasen und/oder Apodisierung an, um einen oder mehr Empfangsstrahlen als Antwort auf jede Sendung zur Erkennung zu formen. Der Empfangsstrahlformer 16 gibt Daten aus, die unter Verwendung der empfangenen akustischen Signale räumliche Orte darstellen. Die Strahlformung wird durch relative Verzögerungen und/oder Phasenverschiebung und Summierung von Signalen aus verschiedenen Elementen bereitgestellt. In alternativen Ausführungsformen ist der Empfangsstrahlformer 16 ein Prozessor zum Erzeugen von Abtastpunkten unter Verwendung der Fourier- oder anderer Transformationen.
  • Der Empfangsstrahlformer 16 kann einen Filter umfassen, wie etwa einen Filter zum Isolieren von Informationen auf einem zweiten harmonischen oder anderen im Verhältnis zu dem Sende-Frequenzband stehenden Frequenzband. Derartige Informationen umfassen eher Informationen über das gewünschte Gewebe, Kontrastmittel und/oder die Strömung. In einer anderen Ausführungsform umfasst der Empfangsstrahlformer 16 einen Speicher oder Puffer und einen Filter oder Addierer. Zwei oder mehr Empfangsstrahlen sind kombiniert, um Informationen auf einem gewünschten Frequenzband zu isolieren, wie etwa einem zweiten harmonischen, kubischen Grundfrequenz- oder anderem Band.
  • Der Empfangsstrahlformer 16 gibt summierte Strahlendaten aus, die räumliche Orte darstellen. Es werden Daten für einen einzelnen Ort, Orte entlang einer Linie, Orte für ein Gebiet oder Orte für ein Volumen ausgegeben. Es kann eine dynamische Fokussierung bereitgestellt werden. Die Daten können unterschiedlichen Zwecken dienen. Zum Beispiel können für B-Mode- oder Gewebedaten andere Abtastungen durchgeführt werden als für die Schubwellenerkennung. Eine andere Möglichkeit besteht darin, dass die B-Mode-Daten ebenfalls zum Ermitteln der durch eine Schubwelle verursachten Verschiebung verwendet werden.
  • Der Prozessor 18 ist ein B-Mode-Detektor, ein Doppler-Detektor, ein gepulster Doppler-Detektor, ein Korrelationsprozessor, ein Fourier-Transformationsprozessor, eine anwendungsspezifische integrierte Schaltung, ein allgemeiner Prozessor, ein Steuerprozessor, ein Bildprozessor, ein Field Programmable Gate Array, ein digitaler Signalprozessor, eine analoge Schaltung, eine digitale Schaltung, Kombinationen derselben oder eine andere jetzt bekannte oder später entwickelte Vorrichtung zur Erkennung und Verarbeitung von Informationen aus strahlengeformten Ultraschallabtastpunkten.
  • In einer Ausführungsform umfasst der Prozessor 18 einen oder mehr Detektoren und einen separaten Prozessor. Der separate Prozessor ist ein Steuerprozessor, ein allgemeiner Prozessor, ein digitaler Signalprozessor, eine anwendungsspezifische integrierte Schaltung, ein Field Programmable Gate Array, ein Netzwerk, ein Server, eine Gruppe von Prozessoren, ein Datenpfad, Kombinationen derselben oder eine andere jetzt bekannte oder später entwickelte Vorrichtung zum Ermitteln von Verschiebung und Berechnen von Gewebeeigenschaften. Zum Beispiel führt der separate Prozessor eine beliebige Kombination einer oder mehrerer der in 1 gezeigten Vorgänge durch. Der Prozessor 18 ist durch Software und/oder Hardware zum Durchführen der Vorgänge konfiguriert.
  • In einer Ausführungsform schätzt der Prozessor 18 die Gewebeverschiebung im Zeitverlauf als Funktion der vom Empfangsstrahlformer 16 ausgegebenen Daten. Die Verschiebungen werden als ein Profil oder als Daten geschätzt, die eine Kurve der Größenordnung der Verschiebung als Funktion der Zeit darstellen. Das Verschiebungsprofil kann gewonnen werden, indem der Grad der Ähnlichkeit korreliert oder auf andere Weise ermittelt wird, der zwischen den Bezugsdaten und den während des Sendens der amplituden- und phasenmodulierten Schiebe-Wellenform gewonnenen oder damit verschachtelten Daten besteht.
  • Der Prozessor 18 ist zum Umwandeln der Verschiebung des Gewebes im Zeitverlauf zum Frequenzbereich konfiguriert. Eine Fourier-Transformation, wie etwa die schnelle Fourier-Transformation, wird auf die Verschiebungsdaten angewendet.
  • Der Prozessor 18 ist zum Berechnen der Schubinformationen als Funktion der Verschiebung des Gewebes im Zeitverlauf konfiguriert. Zum Beispiel wird eine Schubgeschwindigkeit aus der Verschiebung im Zeitverlauf berechnet. Die Geschwindigkeit wird durch den Grad der Verschiebung, dividiert durch die Zeit, bereitgestellt. In einer Ausführungsform berechnet der Prozessor 18 die Viskosität als Funktion der Frequenz oder bei einer Frequenz. Die Viskosität wird aus der Verschiebung in dem Frequenzbereich berechnet. Die transformierte Verschiebung wird verwendet, um die Viskosität in Anhängigkeit von der Frequenz zu ermitteln. Die Frequenzen liegen in einem Bereich, der in Zusammenhang mit der amplituden- und phasenmodulierten Wellenform steht. Die Viskosität auf einer Frequenz, eine durchschnittliche Viskosität oder Viskositäten bei verschiedenen Frequenzen werden in dem Frequenzbereich unter Verwendung einer Schubdarstellung berechnet, können aber auch alternativ im Zeitbereich berechnet werden. Der Prozessor 18 kann andere Eigenschaften berechnen, wie etwa den Schubmodul.
  • Der Prozessor 18 erzeugt Bild- oder Anzeigewerte, die von der viskoelastischen Eigenschaft auf der Anzeige 20 abgebildet werden. Zum Beispiel wird die Schubviskosität, der Schubmodul oder ein anderer Wert ermittelt. Eine Text- oder Zahlenangabe der Eigenschaft wird dem Anwender angezeigt. Ein Schaubild der Viskosität oder einer anderen Eigenschaft im Zeitverlauf oder Frequenzverlauf kann angezeigt werden. In einer Ausführungsform passt der Prozessor 18 der Eigenschaft eine Gerade als Funktion der Frequenz an, und die Steigung wird als Wert oder graphisch anzeigt.
  • In zusätzlichen oder alternativen Ausführungsformen werden die Schubinformationen in Abhängigkeit vom Ort angezeigt. Die Größenordnung der Schubwerte moduliert die Farbe, den Farbton, die Helligkeit und/oder andere Anzeigemerkmale für die verschiedenen Pixel, die eine Geweberegion darstellen. Der Prozessor 18 ermittelt einen Pixelwert (z. B. RGB) oder einen Skalarwert, der in einen Pixelwert umgewandelt wird. Das Bild wird als Skalar- oder Pixelwerte erzeugt. Das Bild kann an einen Videoprozessor, eine Nachschlagetabelle, eine Farbkarte oder direkt an die Anzeige 20 ausgegeben werden.
  • Der Prozessor 18 wird entsprechend den Anweisungen betrieben, die im Speicher 22 oder einem anderen Speicher gespeichert sind. Die Anweisungen konfigurieren den Prozessor 18 zum Betrieb, indem diese geladen und/oder ausgeführt werden. Der Prozessor 18 ist zur In-vivo-Messung einer viskoelastischen Gewebeeigenschaft programmiert. Der Speicher 22 ist ein nichtflüchtiges, computerlesbares Speichermedium. Die Anweisungen zum Implementieren der im vorliegenden Dokument besprochenen Prozesse, Verfahren und/oder Techniken werden auf dem computerlesbaren Speichermedium oder den Speichern, wie etwa einem Cache, Puffer, RAM, Wechselmedium, einer Festplatte oder anderen computerlesbaren Speichermedien, bereitgestellt. Die computerlesbaren Speichermedien umfassen verschiedene Arten von flüchtigen und nichtflüchtigen Speichermedien. Die Funktionen, Vorgänge oder Aufgaben, die in den Figuren veranschaulicht oder im vorliegenden Dokument beschrieben sind, werden. als Reaktion auf einen oder mehrere Sätze von Anweisungen ausgeführt, die in oder auf computerlesbaren Speichermedien gespeichert sind. Die Funktionen, Vorgänge oder Aufgaben sind unabhängig von dem speziellen Typ des Anweisungssatzes, der Speichermedien, des Prozessors oder der Verarbeitungsstrategie und können durch Software, Hardware, integrierte Schaltungen, Firmware, Mikrocode und dergleichen ausgeführt werden, die allein oder in Kombination arbeiten. Ebenso können Verarbeitungsstrategien Multiprocessing, Multitasking, Parallelverarbeitung und dergleichen umfassen. In einer Ausführungsform sind die Anweisungen auf einer Wechsel-Speichermedieneinrichtung gespeichert, um von lokalen oder fernen Systemen gelesen werden zu können. In anderen Ausführungsformen sind die Anweisungen an einem fernen Ort gespeichert, um sie über ein Computernetzwerk oder über Telefonleitungen zu übertragen. In noch anderen Ausführungsformen sind die Anweisungen in einem bestimmten Computer, einer CPU, GPU oder einem System gespeichert.
  • Die Anzeige 20 ist eine Kathodenstrahlröhre, ein LCD-Bildschirm, ein Projektor, ein Plasmabildschirm oder eine andere Anzeige zum Anzeigen von zweidimensionalen Bildern oder dreidimensionalen Darstellungen. Die Anzeige 20 zeigt ein oder mehrere Bilder an, die Schubinformationen darstellen. Das Bild ist ein Schaubild, eine numerische, Text- und/oder zweidimensionale Darstellung einer Region. Zum Beispiel wird ein Viskositätswert oder ein Schaubild der Viskosität als Funktion der Frequenz als das Bild angezeigt.
  • Während die Erfindung oben mit Bezug auf verschiedene Ausführungsformen beschrieben wurde, ist es selbstverständlich, dass viele Änderungen und Abwandlungen vorgenommen werden können, ohne vom Umfang der Erfindung abzuweichen. Es ist daher beabsichtigt, dass die obige detaillierte Beschreibung als erläuternd und nicht als Einschränkung betrachtet wird, und dass verstanden wird, dass es die folgenden Ansprüche einschließlich aller Äquivalente sind, die bezwecken, den Erfindungsgedanken und den Umfang dieser Erfindung zu definieren.

Claims (20)

  1. Nichtflüchtiges, computerlesbares Speichermedium mit darauf gespeicherten Daten, die durch einen programmierten Prozessor (18) ausführbare Anweisungen für die In-vivo-Messung (32) einer viskoelastischen Gewebeeigenschaft darstellen, wobei das Speichermedium Anweisungen umfasst zum: Senden (30) einer amplitudenmodulierten und phasenmodulierten Wellenform, die Zyklen in verschiedenen Frequenzen zu verschiedenen Zeitpunkten einschließt, wobei die verschiedenen Frequenzen in einem Bereich von Frequenzen liegen, wobei die amplitudenmodulierte und phasenmodulierte Wellenform zu dem Gewebe in einem Patienten gesendet wird; Berechnen (34) der Gewebeverschiebung als Funktion der Zeit als Antwort auf die amplitudenmodulierte und phasenmodulierte Wellenform, wobei die Verschiebung aus Abtastungen des Gewebes berechnet wird; Anwenden (36) einer Fourier-Transformation auf die Verschiebung im Zeitverlauf; und Ermitteln (40) der viskoelastischen Eigenschaft aus einer Fourier-Transformation einer Schubwellengleichung und der Fourier-Transformation der Verschiebung im Zeitverlauf.
  2. Computerlesbares Speichermedium nach Anspruch 1, wobei die amplitudenmodulierte und phasenmodulierte Wellenform eine Chirp-Frequenzmodulation innerhalb des Bereichs umfasst.
  3. Computerlesbares Speichermedium nach Anspruch 1, wobei die amplitudenmodulierte und phasenmodulierte Wellenform einen sinusförmigen Träger umfasst, der mit einer sinusförmigen Amplitudenmodulation mit einem phasenvariierenden Term multipliziert wird.
  4. Computerlesbares Speichermedium nach Anspruch 1, wobei das Berechnen (34) der Verschiebung das Abtasten des Gewebes mit Ultraschall, verschachtelt mit dem Senden (30), und das Berechnen (34) der Verschiebung des Gewebes zwischen den Abtastungen als Funktion der Korrelation zwischen Daten aus den Abtastungen umfasst.
  5. Computerlesbares Speichermedium nach Anspruch 1, wobei das Berechnen (34) der Verschiebung umfasst: das Erfassen von k-Raum-Magnetresonanzdaten, die das Gewebe zu verschiedenen Zeitpunkten darstellen, während das Gewebe auf das Senden (30) reagiert, das Subtrahieren der k-Raum-Magnetresonanz-Bezugsdaten, die das Gewebe darstellen, während es frei von Antworten auf das Senden (30) ist, und das Berechnen (34) einer Verschiebung aus den Ergebnissen der Subtraktion.
  6. Computerlesbares Speichermedium nach Anspruch 1, wobei das Anwenden (36) der Fourier-Transformation das Transformieren der Verschiebung im Zeitverlauf in einen Frequenzbereich umfasst, wobei die Verschiebung zu verschiedenen Zeitpunkten auf die verschiedenen Frequenzen reagiert.
  7. Computerlesbares Speichermedium nach Anspruch 1, wobei das Ermitteln (40) der viskoelastischen Eigenschaft das Ermitteln der Viskosität umfasst.
  8. Computerlesbares Speichermedium nach Anspruch 1, wobei das Senden (30) und Ermitteln (40) das Ermitteln (40) der viskoelastischen Eigenschaft über den Bereich der Frequenzen lediglich als Antwort auf die amplitudenmodulierte und phasenmodulierte Wellenform und Abtastungen des Gewebes während der Anwendung der Wellenform umfasst.
  9. Computerlesbares Speichermedium nach Anspruch 1, wobei das Ermitteln (40) das Ermitteln (40) sowohl der Viskosität als Funktion eines imaginären Teils und des Schubmoduls als Funktion eines reellen Teils der Fourier-Transformation der Schubwellengleichung umfasst.
  10. Computerlesbares Speichermedium nach Anspruch 1, welches weiterhin das Auffinden (42) einer Steigung der viskoelastischen Eigenschaft als Funktion der Frequenz in dem Bereich umfasst.
  11. Computerlesbares Speichermedium nach Anspruch 1, welches weiterhin die Bandpassfilterung (38) der Ausgabe der Anwendung der Fourier-Transformation der Verschiebung im Zeitverlauf umfasst, wobei der Bandpass etwa den zweifachen Bereich an Breite und hohen und niedrigen Frequenzen aufweist.
  12. Verfahren zur viskoelastischen Messung unter Verwendung von Ultraschall, wobei das Verfahren umfasst: Senden (30) einer amplituden- und phasenmodulierten Ultraschallwellenform zu einem Patienten in vivo; und M essen (32) einer viskoelastischen Eigenschaft aus dem durch die amplituden- und phasenmodulierte Ultraschallwellenform verursachten Schub, wobei die viskoelastische Eigenschaft über einen Bereich von Frequenzen gemessen wird, der in Zusammenhang mit der amplituden- und phasenmodulierten Ultraschallwellenform steht und unabhängig von einem Grad an Kompression ist, die während des Sendens von außen auf den Patienten angewendet wird (30).
  13. Verfahren nach Anspruch 12, wobei das Senden (30) das Multiplizieren eines sinusförmigen Trägers bei einer Ultraschallfrequenz mit einer Hüllkurve umfasst, die die Amplituden- und Phasenmodulation aufweist, sodass der Bereich der Frequenzen zwischen 15 Hz und 75 Hz liegt.
  14. Verfahren nach Anspruch 12, wobei das Messen (32) umfasst: Abtasten des Patienten mit Ultraschall, verschachtelt mit dem Senden (30); und Berechnen (34) der durch das Senden (30) verursachten Verschiebung von Gewebe des Patienten zwischen den Abtastungen des Abtastens, wobei das Berechnen (34) als Funktion der Korrelation zwischen Daten aus den Abtastungen erfolgt.
  15. Verfahren nach Anspruch 12, wobei das Messen (32) umfasst: Erfassen von k-Raum-Magnetresonanzdaten, die das Gewebe des Patienten zu verschiedenen Zeitpunkten darstellen, während das Gewebe auf das Senden (30) reagiert, das Subtrahieren der k-Raum-Magnetresonanz-Bezugsdaten, die das Gewebe darstellen, während es frei von Antworten auf das Senden (30) ist, und das Berechnen (34) einer Verschiebung aus den Ergebnissen der Subtraktion.
  16. Verfahren nach Anspruch 12, wobei das Messen (32) das Anwenden (36) einer Fourier-Transformation auf Daten, die die Verschiebung als Funktion der Zeit darstellen, und das Berechnen (34) der Viskosität als Funktion von Ergebnissen der Anwendung (36) der Fourier-Transformation umfasst.
  17. Verfahren nach Anspruch 16, wobei das Berechnen (34) das Berechnen (34) als Funktion eines imaginären Teils einer Fourier-Transformation der Schubwellengleichung umfasst.
  18. Verfahren nach Anspruch 12, die weiterhin das Auffinden (42) einer Steigung der viskoelastischen Eigenschaft als Funktion der Frequenz in dem Bereich umfasst.
  19. System zur Schubwellenberechnung unter Verwendung von Ultraschall, wobei das System umfasst: einen Sendestrahlformer (12), der betrieblich einsetzbar ist, um eine amplituden- und phasenmodulierte Wellenform zu erzeugen; einen Ultraschall-Schallkopf (14), der so mit dem Sendestrahlformer (12) in Verbindung steht, dass der Ultraschall-Schallkopf akustische Energie als Antwort auf die amplituden- und phasenmodulierte Wellenform in einem Patienten sendet; einen Empfangsstrahlformer (16), der betrieblich einsetzbar ist, um Daten auszugeben, die räumliche Orte als Funktion der empfangenen akustischen Signale darstellen; einen Prozessor (18), der dazu konfiguriert ist, die Verschiebung in einem Patienten im Zeitverlauf als Funktion der ausgegebenen Daten abzuschätzen und Schubinformationen als Funktion der Verschiebung in dem Patienten im Zeitverlauf zu berechnen. eine Anzeige (20), die betrieblich einsetzbar ist, um ein Bild anzuzeigen, das eine Funktion der Schubinformationen ist.
  20. System nach Anspruch 19, wobei der Prozessor (18) konfiguriert ist, um die Verschiebung in dem Patienten im Zeitverlauf in einen Frequenzbereich zu transformieren und die Viskosität und den Schubmodul als Funktion der Frequenz aus der Verschiebung in dem Frequenzbereich zu berechnen, wobei die Frequenzen in einem Bereich liegen, der in Zusammenhang mit der amplituden- und phasenmodulierten Wellenform steht.
DE102012101312.8A 2011-02-17 2012-02-17 Viskoelastizitätsmessung unter Verwendung einer amplitudenphasenmodulierten Ultraschallwelle Active DE102012101312B4 (de)

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
US13/029,369 US8469891B2 (en) 2011-02-17 2011-02-17 Viscoelasticity measurement using amplitude-phase modulated ultrasound wave
US13/029,369 2011-02-17

Publications (2)

Publication Number Publication Date
DE102012101312A1 true DE102012101312A1 (de) 2012-08-23
DE102012101312B4 DE102012101312B4 (de) 2022-08-25

Family

ID=46605148

Family Applications (1)

Application Number Title Priority Date Filing Date
DE102012101312.8A Active DE102012101312B4 (de) 2011-02-17 2012-02-17 Viskoelastizitätsmessung unter Verwendung einer amplitudenphasenmodulierten Ultraschallwelle

Country Status (6)

Country Link
US (1) US8469891B2 (de)
JP (1) JP5984417B2 (de)
KR (1) KR101398948B1 (de)
CN (1) CN102641137B (de)
DE (1) DE102012101312B4 (de)
FR (1) FR2971695B1 (de)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114224382A (zh) * 2021-12-17 2022-03-25 重庆医科大学 粘弹性测量方法及其系统

Families Citing this family (53)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102123666B (zh) * 2008-08-15 2013-09-25 皇家飞利浦电子股份有限公司 用于采集材料的声弹性成像数据和超声数据的换能器装置和方法
US8715185B2 (en) * 2010-04-05 2014-05-06 Hitachi Aloka Medical, Ltd. Methods and apparatus for ultrasound imaging
EP2676143B1 (de) 2011-02-15 2023-11-01 Hemosonics, Llc Charakterisierung von hämostase- und sauerstofftransportparametern
US20120253194A1 (en) * 2011-03-30 2012-10-04 Tadashi Tamura Methods and apparatus for ultrasound imaging
US10667791B2 (en) * 2011-08-19 2020-06-02 The University Of British Columbia Elastography using ultrasound imaging of a thin volume
US9239373B2 (en) * 2011-11-16 2016-01-19 Siemens Medical Solutions Usa, Inc. Adaptive image optimization in induced wave ultrasound imaging
US9125615B2 (en) * 2012-09-14 2015-09-08 Mayo Foundation For Medical Education And Research System and method for non-invasive measurement of carpal tunnel pressure
US9642600B2 (en) * 2012-10-01 2017-05-09 Mayo Foundation For Medical Education And Research Shear wave attenuation from k-space analysis system
CN103800038B (zh) * 2012-11-12 2016-09-21 通用电气公司 改善的系统和装置以用于确定目标组织的机械特性
EP2965694A4 (de) * 2013-03-05 2016-11-16 Hitachi Ltd Ultraschalldiagnosevorrichtung und sende-/empfangsverfahren
US9332962B2 (en) * 2013-03-13 2016-05-10 Siemens Medical Solutions Usa, Inc. Ultrasound ARFI displacement imaging using an adaptive time instance
EP2979645A4 (de) 2013-03-28 2016-12-28 Canon Kk Ultraschalldiagnosevorrichtung und ultraschalldiagnoseverfahren
US9965852B2 (en) * 2013-04-18 2018-05-08 The United States Of America, As Represented By The Secretary Of The Navy System and method for evaluating anisotropic viscoelastic properties of fibrous structures
JP6447879B2 (ja) * 2013-04-22 2019-01-09 ソニー株式会社 超音波処理装置および方法、並びにプログラム
WO2014201020A1 (en) * 2013-06-10 2014-12-18 Mayo Foundation For Medical Education And Research System and method for acoustic radiation force creep-recovery and shear wave propagation for elasticity imaging
US10512450B2 (en) * 2013-09-24 2019-12-24 Siemens Medical Solutions Usa, Inc. Shear wave estimation from analytic data
US10143442B2 (en) 2013-10-24 2018-12-04 Ge Medical Systems Global Technology, Llc Ultrasonic diagnosis apparatus
JP6263372B2 (ja) * 2013-11-28 2018-01-17 ジーイー・メディカル・システムズ・グローバル・テクノロジー・カンパニー・エルエルシー 超音波診断装置
US9332963B2 (en) 2014-01-21 2016-05-10 Siemens Medical Solutions Usa, Inc. Swept focus for acoustic radiation force impulse
JP6293578B2 (ja) * 2014-05-28 2018-03-14 ジーイー・メディカル・システムズ・グローバル・テクノロジー・カンパニー・エルエルシー 超音波診断装置及びプログラム
EP3215018B1 (de) * 2014-10-29 2023-12-06 Mayo Foundation for Medical Education and Research Verfahren zur ultraschallelastografie durch kontinuierliche vibration eines ultraschallwandlers
DE102015101891A1 (de) * 2015-02-10 2016-08-11 Endress + Hauser Gmbh + Co. Kg Vorrichtung zur Bestimmung und/oder Überwachung zumindest einer Prozessgröße eines Mediums
US9726647B2 (en) 2015-03-17 2017-08-08 Hemosonics, Llc Determining mechanical properties via ultrasound-induced resonance
US10376242B2 (en) * 2015-04-16 2019-08-13 Siemens Medical Solutions Usa, Inc. Quantitative viscoelastic ultrasound imaging
US9814446B2 (en) * 2015-04-22 2017-11-14 Siemens Medical Solutions Usa, Inc. Method and system for automatic estimation of shear modulus and viscosity from shear wave imaging
JP6608232B2 (ja) * 2015-09-30 2019-11-20 キヤノンメディカルシステムズ株式会社 医用画像診断装置、医用画像処理装置および医用情報の表示制御方法
WO2017062553A1 (en) * 2015-10-08 2017-04-13 Mayo Foundation For Medical Education And Research Systems and methods for ultrasound elastography with continuous transducer vibration
CN106802969B (zh) * 2015-11-26 2020-08-07 英业达科技有限公司 阻尼材料动态特性的验证系统及其验证方法
CN105455851B (zh) 2015-12-24 2018-03-13 无锡海斯凯尔医学技术有限公司 粘弹性介质的粘弹性参数检测方法和设备
US11717256B2 (en) * 2016-01-08 2023-08-08 Siemens Medical Solutions Usa, Inc. Motion independence in acoustic radiation force impulse imaging
US10376233B2 (en) * 2016-04-08 2019-08-13 Siemens Medical Solutions Usa, Inc. Diffraction source compensation in medical diagnostic ultrasound viscoelastic imaging
US10675007B2 (en) 2016-04-19 2020-06-09 Siemens Medical Solutions Usa, Inc. Frequency compounding in elasticity imaging
CN106175831B (zh) * 2016-09-28 2018-10-26 重庆理工大学 一种基于剪切波幅度和相位检测的组织粘弹性测量方法
CN108872903B (zh) * 2017-05-10 2024-06-14 冯原 头部磁共振弹性成像检测方法及成像驱动装置
US20180310918A1 (en) * 2017-04-27 2018-11-01 Siemens Medical Solutions Usa, Inc. Variable focus for shear wave imaging
WO2018224602A1 (en) * 2017-06-09 2018-12-13 Technische Universiteit Eindhoven Shear wave viscoelasticity imaging using local system identification
CN107505232B (zh) 2017-07-21 2019-09-03 无锡海斯凯尔医学技术有限公司 运动信息获取方法及装置
CN107478723A (zh) * 2017-07-21 2017-12-15 无锡海斯凯尔医学技术有限公司 介质粘弹性的测量方法和装置
US11644440B2 (en) 2017-08-10 2023-05-09 Mayo Foundation For Medical Education And Research Shear wave elastography with ultrasound probe oscillation
JP7253540B2 (ja) * 2017-10-12 2023-04-06 コーニンクレッカ フィリップス エヌ ヴェ 患者適応型のせん断波生成による超音波せん断波イメージング
US11154277B2 (en) * 2017-10-31 2021-10-26 Siemens Medical Solutions Usa, Inc. Tissue viscoelastic estimation from shear velocity in ultrasound medical imaging
CN107753058B (zh) * 2017-11-22 2021-03-02 深圳中科乐普医疗技术有限公司 一种剪切波动态滤波方法
WO2019131811A1 (ja) * 2017-12-27 2019-07-04 キヤノンメディカルシステムズ株式会社 解析装置
CN116869567A (zh) * 2018-05-15 2023-10-13 深圳迈瑞生物医疗电子股份有限公司 一种剪切波弹性测量方法及剪切波弹性成像系统
US11467075B1 (en) * 2018-09-21 2022-10-11 The United States Of America As Represented By The Secretary Of The Navy Method of estimating material properties of an elastic plate
US11980502B2 (en) 2018-11-29 2024-05-14 The Penn State Research Foundation Harmonic shear wave imaging
CN113229846B (zh) * 2018-12-24 2023-10-13 深圳迈瑞生物医疗电子股份有限公司 超声成像方法与设备
CN110547825B (zh) * 2019-09-03 2022-04-26 上海交通大学 一种基于多频剪切波幅值分析的超声弹性成像技术
CN112438751A (zh) * 2019-09-04 2021-03-05 通用电气精准医疗有限责任公司 用于剪切波弹性成像的方法和系统以及存储有对应程序的介质
CN110742644B (zh) * 2019-09-29 2022-07-08 深圳大学 弹性成像系统、弹性成像方法及存储介质
CN112507922B (zh) * 2020-12-16 2023-11-07 平安银行股份有限公司 人脸活体检测方法、装置、电子设备及存储介质
CN116528769A (zh) * 2021-01-04 2023-08-01 深圳迈瑞生物医疗电子股份有限公司 粘弹性测量方法和超声成像系统
CN113281368B (zh) * 2021-05-19 2024-04-02 成都鸣石峻致科技有限公司 一种磁共振弹性测量方法、装置、计算机设备、系统及存储介质

Family Cites Families (28)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2845583B2 (ja) * 1990-07-10 1999-01-13 株式会社資生堂 生体表面部の力学特性の測定方法と装置並びにそのためのセンサ内蔵形加振器の構造
US5810731A (en) 1995-11-13 1998-09-22 Artann Laboratories Method and apparatus for elasticity imaging using remotely induced shear wave
US5606971A (en) 1995-11-13 1997-03-04 Artann Corporation, A Nj Corp. Method and device for shear wave elasticity imaging
JPH1176176A (ja) * 1997-09-11 1999-03-23 Nec Corp 生体組織粘弾性特性計測器とその方法および生体組織粘弾性特性算出プログラムを記録した媒体
US6246895B1 (en) 1998-12-18 2001-06-12 Sunnybrook Health Science Centre Imaging of ultrasonic fields with MRI
WO2000070362A1 (en) 1999-05-14 2000-11-23 Koninklijke Philips Electronics N.V. Mr elastography method
CA2428872C (en) * 2000-11-28 2013-01-08 Allez Physionix Limited Systems and methods for making non-invasive physiological assessments
US6735461B2 (en) 2001-06-19 2004-05-11 Insightec-Txsonics Ltd Focused ultrasound system with MRI synchronization
US6879155B2 (en) 2002-04-10 2005-04-12 Mayo Foundation For Medical Education And Research Magnetic resonance acoustography
US6984209B2 (en) 2003-07-02 2006-01-10 The Brigham And Women's Hospital, Inc. Harmonic motion imaging
US7753847B2 (en) 2003-10-03 2010-07-13 Mayo Foundation For Medical Education And Research Ultrasound vibrometry
JP4189840B2 (ja) * 2003-10-20 2008-12-03 独立行政法人産業技術総合研究所 超音波を利用した軟組織の粘弾性推定装置およびプログラム
WO2006044997A2 (en) * 2004-10-15 2006-04-27 The Trustees Of Columbia University In The City Of New York System and method for localized measurement and imaging of viscosity of tissues
GB0427686D0 (en) 2004-12-17 2005-02-02 Univ Basel Method for detection and imaging of synchronous spin and charged particle motion
US7307423B2 (en) 2005-05-05 2007-12-11 Wisconsin A.Umni Research Foundation Magnetic resonance elastography using multiple drivers
US7414705B2 (en) * 2005-11-29 2008-08-19 Navisense Method and system for range measurement
WO2008139245A1 (en) * 2007-05-16 2008-11-20 Super Sonic Imagine Method and device for measuring a mean value of visco-elasticity of a region of interest
US8137275B2 (en) * 2007-06-28 2012-03-20 Siemens Medical Solutions Usa, Inc. Tissue complex modulus and/or viscosity ultrasound imaging
JP5349173B2 (ja) * 2008-07-14 2013-11-20 株式会社東芝 超音波診断装置及び超音波診断装置制御プログラム
EP2310829B1 (de) * 2008-07-30 2016-04-06 Centre Hospitalier De L'Universite de Montreal System und verfahren zur erkennung, charakterisierung und abbildung von heterogenitäten mithilfe einer durch scherwellen hervorgerufenen resonanz
US9364194B2 (en) * 2008-09-18 2016-06-14 General Electric Company Systems and methods for detecting regions of altered stiffness
US8394026B2 (en) * 2008-11-03 2013-03-12 University Of British Columbia Method and apparatus for determining viscoelastic parameters in tissue
US8398550B2 (en) 2008-12-01 2013-03-19 The Board Of Trustees Of The University Of Illinois Techniques to evaluate mechanical properties of a biologic material
WO2010104863A1 (en) 2009-03-09 2010-09-16 Mayo Foundation For Medical Education And Research Method for ultrasound vibrometry using orthogonal basis functions
US8343050B2 (en) * 2009-05-04 2013-01-01 Siemens Medical Solutions Usa, Inc. Feedback in medical ultrasound imaging for high intensity focused ultrasound
US10172527B2 (en) * 2009-07-31 2019-01-08 Supersonic Imagine Method and apparatus for measuring a physical parameter in mammal soft tissues by propagating shear waves
CN101699280B (zh) * 2009-10-15 2011-08-17 北京索瑞特医学技术有限公司 超声无损检测粘弹性介质弹性的方法及其装置
US9986973B2 (en) * 2010-04-23 2018-06-05 Mayo Foundation For Medical Education And Research Method for shear wave ultrasound vibrometry with interleaved push and detection pulses

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114224382A (zh) * 2021-12-17 2022-03-25 重庆医科大学 粘弹性测量方法及其系统
CN114224382B (zh) * 2021-12-17 2023-09-15 重庆医科大学 粘弹性测量方法及其系统

Also Published As

Publication number Publication date
CN102641137B (zh) 2014-11-12
FR2971695B1 (fr) 2022-08-05
KR101398948B1 (ko) 2014-05-27
KR20120102510A (ko) 2012-09-18
US8469891B2 (en) 2013-06-25
JP2012170823A (ja) 2012-09-10
DE102012101312B4 (de) 2022-08-25
CN102641137A (zh) 2012-08-22
US20120215101A1 (en) 2012-08-23
FR2971695A1 (fr) 2012-08-24
JP5984417B2 (ja) 2016-09-06

Similar Documents

Publication Publication Date Title
DE102012101312B4 (de) Viskoelastizitätsmessung unter Verwendung einer amplitudenphasenmodulierten Ultraschallwelle
DE102013010060B4 (de) Messen der akustischen Absorption oder Dämpfung von Ultraschall
DE102016100367B4 (de) Spärliche Verfolgung in Schallstrahlintensitätsimpuls-Bildgebung
DE102012108353B4 (de) Klassifikationsvorverarbeitung in der medizinischen ultraschallscherwellenbildgebung
DE102017211895A1 (de) Gewebecharakterisierung im medizinischen diagnostischen Ultraschall
EP3215018B1 (de) Verfahren zur ultraschallelastografie durch kontinuierliche vibration eines ultraschallwandlers
DE102013001230B4 (de) Achsbezogene Charakterisierung von Scherwellen mit Ultraschall
DE102011114943B4 (de) Bestimmen von Scherwelleninformationen in der medizinischen Ultraschallbildgebung
DE69011686T2 (de) Verfahren und vorrichtung zur messung und darstellung der gewebekompressibilität oder nachgiebigkeit.
DE102017202141B4 (de) Schallgeschwindigkeitsbildgebung unter Verwendung von Scherwellen
DE102010018857A1 (de) Feedback in der medizinischen Ultraschallbildgebung für hochintensiven fokussierten Ultraschall
US10743814B2 (en) Fat fraction estimation using ultrasound with shear wave propagation
DE102012022010A1 (de) Adaptive Bildoptimierung in der welleninduzierten Ultraschallbildgebung
DE102016106998A1 (de) Quantitative viskoelastische Ultraschallbildgebung
DE102014002747A1 (de) Arfi-ultraschall-verlagerungsbildgebung mittels eines adaptiven zeitpunkts
DE102013005038A1 (de) Fusion von parametrischen Magnetresonanz- und Ultraschallbildern
DE102018218319A1 (de) Viskoelastische Gewebeschätzung von Schergeschwindigkeit in medizinischer Ultraschallbildgebung
DE102014003105A1 (de) Fettanteilschatzung mittels ultraschall mit scherwellenausbreitung
DE102019202545A1 (de) Platzierung der Region von Interesse für die quantitative Ultraschallbildgebung
US9173629B2 (en) Ultrasonic diagnostic apparatus and ultrasonic image processing apparatus
DE102009033286A1 (de) Scherwellenbildgebung
DE102016114783A1 (de) Adaptive Bewegungsschätzung bei Schallstrahlungskraftbildgebung
DE102011014301A1 (de) Motion synchronized destruction for three-dimensional reperfusion mapping in medical diagnostic ultrasound imaging
DE102019107227A1 (de) Frequenzdurchlauf für akustische Strahlungskraftimpulse
DE102017206463A1 (de) Frequenz-compounding bei der elastizitätsbildgebung

Legal Events

Date Code Title Description
R012 Request for examination validly filed
R016 Response to examination communication
R016 Response to examination communication
R018 Grant decision by examination section/examining division
R020 Patent grant now final