DE102012108353A1 - Klassifikationsvorverarbeitung in der medizinischen ultraschallscherwellenbildgebung - Google Patents

Klassifikationsvorverarbeitung in der medizinischen ultraschallscherwellenbildgebung Download PDF

Info

Publication number
DE102012108353A1
DE102012108353A1 DE102012108353A DE102012108353A DE102012108353A1 DE 102012108353 A1 DE102012108353 A1 DE 102012108353A1 DE 102012108353 A DE102012108353 A DE 102012108353A DE 102012108353 A DE102012108353 A DE 102012108353A DE 102012108353 A1 DE102012108353 A1 DE 102012108353A1
Authority
DE
Germany
Prior art keywords
displacement
tissue
signal
over time
profile
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
DE102012108353A
Other languages
English (en)
Other versions
DE102012108353B4 (de
Inventor
Fan Liexiang
Paul Freiburger
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 DE102012108353A1 publication Critical patent/DE102012108353A1/de
Application granted granted Critical
Publication of DE102012108353B4 publication Critical patent/DE102012108353B4/de
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • 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
    • 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
    • A61B8/0858Detecting organic movements or changes, e.g. tumours, cysts, swellings involving measuring tissue layers, e.g. skin, interfaces
    • 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/52Devices using data or image processing specially adapted for diagnosis using ultrasonic, sonic or infrasonic waves
    • A61B8/5207Devices using data or image processing specially adapted for diagnosis using ultrasonic, sonic or infrasonic waves involving processing of raw data to produce diagnostic data, e.g. for generating an image
    • 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/5223Devices using data or image processing specially adapted for diagnosis using ultrasonic, sonic or infrasonic waves involving processing of medical diagnostic data for extracting a diagnostic or physiological parameter from medical diagnostic data
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00Image analysis
    • GPHYSICS
    • G16INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
    • G16HHEALTHCARE INFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR THE HANDLING OR PROCESSING OF MEDICAL OR HEALTHCARE DATA
    • G16H50/00ICT specially adapted for medical diagnosis, medical simulation or medical data mining; ICT specially adapted for detecting, monitoring or modelling epidemics or pandemics
    • G16H50/30ICT specially adapted for medical diagnosis, medical simulation or medical data mining; ICT specially adapted for detecting, monitoring or modelling epidemics or pandemics for calculating health indices; for individual health risk assessment
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B8/00Diagnosis using ultrasonic, sonic or infrasonic waves
    • A61B8/13Tomography
    • A61B8/14Echo-tomography
    • 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/52019Details of transmitters
    • G01S7/5202Details of transmitters for pulse systems
    • G01S7/52022Details of transmitters for pulse systems using a sequence of pulses, at least one pulse manipulating the transmissivity or reflexivity of the medium
    • GPHYSICS
    • G03PHOTOGRAPHY; CINEMATOGRAPHY; ANALOGOUS TECHNIQUES USING WAVES OTHER THAN OPTICAL WAVES; ELECTROGRAPHY; HOLOGRAPHY
    • G03BAPPARATUS OR ARRANGEMENTS FOR TAKING PHOTOGRAPHS OR FOR PROJECTING OR VIEWING THEM; APPARATUS OR ARRANGEMENTS EMPLOYING ANALOGOUS TECHNIQUES USING WAVES OTHER THAN OPTICAL WAVES; ACCESSORIES THEREFOR
    • G03B42/00Obtaining records using waves other than optical waves; Visualisation of such records by using optical means
    • G03B42/06Obtaining records using waves other than optical waves; Visualisation of such records by using optical means using ultrasonic, sonic or infrasonic waves

Landscapes

  • Health & Medical Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Public Health (AREA)
  • Medical Informatics (AREA)
  • Biomedical Technology (AREA)
  • General Health & Medical Sciences (AREA)
  • Pathology (AREA)
  • Veterinary Medicine (AREA)
  • Animal Behavior & Ethology (AREA)
  • Surgery (AREA)
  • Molecular Biology (AREA)
  • Biophysics (AREA)
  • Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
  • Radiology & Medical Imaging (AREA)
  • Heart & Thoracic Surgery (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • General Physics & Mathematics (AREA)
  • Physiology (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • Radar, Positioning & Navigation (AREA)
  • Remote Sensing (AREA)
  • Data Mining & Analysis (AREA)
  • Primary Health Care (AREA)
  • Epidemiology (AREA)
  • Databases & Information Systems (AREA)
  • Theoretical Computer Science (AREA)
  • Ultra Sonic Daignosis Equipment (AREA)

Abstract

Es wird eine Klassifikationsvorverarbeitung (30–40) für die medizinische Ultraschallscherwellenbildgebung (46) bereitgestellt. Als Antwort auf Belastung wird die Verlagerung an einem oder mehr Orten in einem Patienten gemessen (32). Die Verlagerung im Zeitverlauf ist eine Kurve, die eine Ortsverschiebung darstellt. Ein oder mehr Merkmale der Kurve, wie etwa Signal-Rausch-Verhältnis und maximale Verlagerung, werden zum Klassifizieren (42) des Orts verwendet. Der Ort wird als Flüssigkeit oder flüssiges Gewebe, festes Gewebe oder nicht ermittelbar klassifiziert. Eine nachfolgende Scherungsbildgebung (46) kann Scherungsinformationen für Orte mit festem Gewebe und nicht an anderen Orten bereitstellen.

Description

  • HINTERGRUND
  • Die vorliegenden Ausführungsformen beziehen sich auf die Ultraschallbildgebung. Insbesondere kann die Ultraschallscherwellenbildgebung verbessert werden.
  • Informationen zur Scherwellengeschwindigkeit können nützlich für die Diagnose sein. Scherwelleninformationen können zusätzlich zu akustischer Impedanz-Bildgebung (z. B. B-Mode) und Doppler-Bildgebung (z. B. Strömungsmodus) Gewebemerkmale anzeigen. Allerdings kann die Komplexität der Scherwellenfortpflanzung im Gewebe zu wesentlichen Fehlern führen. Zum Beispiel wird die Scherwellengeschwindigkeit unter Umständen für Flüssigkeit oder flüssiges Gewebe weniger genau ermittelt, was zu scheinbar willkürlichen Werten führt.
  • Um die Flüssigkeitseffekte aus Scherwellenbildern zu entfernen, können flüssige Regionen von festen Geweberegionen getrennt werden, indem die Strömung nach dem Anregen gemessen wird. Zum Anregen wird Schallenergie verwendet, um eine Strömung in der Flüssigkeit zu erzeugen. Ein Korrelationskoeffizient zwischen aufeinanderfolgenden oder sequenziellen Echosignalen wird berechnet, um festes Gewebe von der angeregten Flüssigkeit zu unterscheiden. Bei der Live-Abtastung kann allerdings die Bewegung des festen Gewebes die gleiche Größenordnung an Dekorrelation auslösen wie Flüssigkeitsbewegung, die durch das Schallenergieanregen hervorgerufen werden. Die Unterscheidungsfähigkeit ist unter Umständen beschränkt.
  • KURZE ZUSAMMENFASSUNG
  • Als Einleitung enthalten die im Folgenden beschriebenen bevorzugten Ausführungsformen Verfahren, Anweisungen und Systeme zur Klassifikationsvorverarbeitung für die medizinische Ultraschallscherwellenbildgebung. Als Antwort auf Belastung wird die Verlagerung an einem oder mehreren Orten in einem Patienten gemessen. Die Verlagerung im Zeitverlauf ist eine Kurve, die die Verschiebung der Position des Orts darstellt. Ein oder mehr Merkmale der Kurve, wie etwa das Signal-Rausch-Verhältnis und die maximale Verlagerung, werden zum Klassifizieren des Orts verwendet. Der Ort wird als Flüssigkeit oder flüssiges Gewebe, festes Gewebe, nicht ermittelbar oder andere Klasse klassifiziert. Eine nachfolgende Scherungsbildgebung kann Scherungsinformationen für Orte mit festem Gewebe und nicht an anderen Orten bereitstellen.
  • In einem ersten Aspekt wird ein Verfahren für die Klassifikationsvorverarbeitung zur Verwendung in der medizinischen Ultraschallscherwellenbildgebung bereitgestellt. Die Verlagerung im Zeitverlauf an einem ersten Ort im Innern eines Patienten als Antwort auf eine Impulsanregung wird mittels Ultraschall. gemessen. Ein Maximalwert der Verlagerung im Zeitverlauf wird ermittelt. Ein Signal-Rausch-Verhältnis der Verlagerung im Zeitverlauf wird ermittelt. Ein Prozessor klassifiziert den ersten Ort als einen ersten einer Vielzahl von Gewebetypen. Die Klassifikation ist eine Funktion des Maximalwerts der Verlagerung im Zeitverlauf und des Signal-Rausch-Verhältnisses der Verlagerung im Zeitverlauf.
  • In einem zweiten Aspekt sind auf einem nichtflüchtigen, computerlesbaren Speichermedium Daten gespeichert, die durch einen programmierten Prozessor ausführbare Anweisungen für die Klassifikationsvorverarbeitung in der medizinischen Ultraschallscherwellenbildgebung darstellen. Das Speichermedium umfasst Anweisungen zum Ermitteln eines Verlagerungsprofils der Antwort in einem Patienten, zum Berechnen eines ersten Merkmals des Verlagerungsprofils, zum Unterscheiden von Flüssigkeit und flüssigem Gewebe von festem Gewebe, was zumindest teilweise auf dem ersten Merkmal basiert, zum Ausmaskieren von Scherwelleninformationen aus einem Bild für Flüssigkeit oder flüssiges Gewebe.
  • In einem dritten Aspekt wird ein System für die Klassifikationsvorverarbeitung in der medizinischen Ultraschallscherwellenbildgebung bereitgestellt. Ein Schallkopf ist zum Senden einer akustischen Impulsanregung in einen Patienten konfiguriert und zum Abtasten einer Region des Patienten mit Ultraschall konfiguriert. Ein Empfangsstrahlformer ist dazu konfiguriert, Daten zu erzeugen, welche die Region zu verschiedenen Zeitpunkten nach der akustischen Impulsanregung darstellen. Die Daten werden aus der Abtastung mittels Ultraschall erzeugt. Ein Prozessor ist dazu konfiguriert, die durch die akustische Impulsanregung ausgelöste Gewebeverlagerung abzuschätzen, mindestens einen Parameter abzuleiten, der ein Merkmal eines Profils der Gewebeverlagerung beschreibt, und das Gewebe in der Region auf der Grundlage des mindestens einen Parameters zu klassifizieren. Eine Anzeige ist zum Anzeigen eines Bilds konfiguriert, das die Scherung für erste Orte darstellt. Die ersten Orte werden als eine Funktion der Klassifikation des Gewebes ermittelt.
  • 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 der Schwerpunkt 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 für die Klassifikationsvorverarbeitung zur Verwendung in der medizinischen Ultraschallscherwellenbildgebung;
  • 2 ist ein beispielhaftes medizinisches Ultraschallbild;
  • 3 ist ein beispielhaftes Schaubild der Verlagerung im Zeitverlauf von zwei Orten, die in 2 dargestellt sind;
  • 4 zeigt beispielhafte geglättete und Rauschverlagerungskurven für die zwei in 2 dargestellten Orte;
  • 5 ist eine graphische Darstellung der Klassifikation auf der Grundlage der maximalen Verlagerung und des Signal-Rausch-Verhältnisses, gemäß einer Ausführungsform;
  • 6 zeigt zwei beispielhafte medizinische Bilder mit einem B-Mode-Bild auf der linken Seite und einem aus dem Signal-Rausch-Verhältnis erzeugten Bild auf der rechten Seite; und
  • 7 ist ein Blockschaltplan einer Ausführungsform eines Systems für die Klassifikationsvorverarbeitung in der medizinischen Ultraschallscherwellenbildgebung.
  • DETAILLIERTE BESCHREIBUNG DER ZEICHNUNGEN UND GEGENWÄRTIG BEVORZUGTE AUSFÜHRUNGSFORMEN
  • Orte in einem Untersuchungsbereich werden als festes Gewebe oder andere Materialtypen klassifiziert. Feste Geweberegionen erlauben unter Umständen eine genauere Schermessung. Andere Typen, wie etwa Flüssigkeit oder flüssiges Gewebe, können unter Umständen weniger genau auf Schergeschwindigkeit gemessen werden. Die mit festem Gewebe in Zusammenhang stehenden Orte werden in die Scherungsbildgebung eingeschlossen, und andere Orte können vorab von der Bildgebung des Scherwellenfortpflanzungsfelds ausgeschlossen werden. Die Klassifikation kann alternativ oder zusätzlich mit Abschattung verbundene Regionen und/oder echoarme Regionen feststellen, um diese aus der Scherwellenbildgebung auszuschließen.
  • Ein Untersuchungsbereich kann Flüssigkeit (z. B. Blut oder Urin), Knochen und Gewebe umfassen. Gewebe kann flüssigkeitsartiges Gewebe (d. h. flüssiges Gewebe) oder festes Gewebe einschließen. Zu festen Geweben zählen Muskeln, Fett, Organe oder andere Gebilde mit einer verhältnismäßig elastischeren Beschaffenheit als Flüssigkeit und Knochen. Zysten und andere Gewebsstrukturen können einen hohen Flüssigkeitsgehalt aufweisen. Um Ungenauigkeiten in der Scherwellenbildgebung zu reduzieren, wird das flüssige Gewebe oder die Flüssigkeit von festem Gewebe segmentiert.
  • Wenn festes Gewebe und flüssiges Gewebe durch eine Impulskraft angeregt werden, kommt es bei beiden Gewebe zur Verlagerung. Das Merkmal des Verlagerungsprofils im Zeitverlauf kann je nach Gewebetyp unterschiedlich sein. Die Verlagerung im festen Gewebe wird von der Scherwellengleichung bestimmt, und die Verlagerung in flüssigem Gewebe wird von der Navier-Stokes-Gleichung bestimmt. Das sich ergebende Verlagerungsprofil in flüssigem Gewebe zeigt ein im Vergleich mit dem Verlagerungsprofil in festem Gewebe signifikantes Rauschen. Das Signal-Rausch-Verhältnis (SRV) des Echosignals trägt ebenfalls zu dem SRV des Verlagerungsprofils bei.
  • Andere Parameter des Verlagerungsprofils können zusätzlich zu dem SRV oder als Alternative zum SRV für die Klassifikation eingesetzt werden. Zum Beispiel kann die maximale Verlagerung im Verlagerungsprofil Flüssigkeit oder flüssiges Gewebe von einer echoarmen Region festen Gewebes unterscheiden. Werden das SRV und die maximale Verlagerung zusammen verwendet, können sie Flüssigkeit und flüssiges Gewebe von weichen und harten festen Geweben abgrenzen.
  • Gewebeverlagerungen, die durch Impulsanregung ausgelöst werden, werden geschätzt. Mindestens ein Parameter, der die Merkmale des Verlagerungsprofils beschreibt, wird geschätzt. Gewebe werden basierend auf dem abgeleiteten Parameter in mindestens eine von einer Vielzahl von Kategorien klassifiziert. In einem Beispiel werden das SRV und der Maximalwert des Verlagerungsprofils zum Klassifizieren von Gewebe in flüssige oder feste Gewebetypen verwendet. In einem anderen Beispiel wird ein anderes Merkmal wie etwa eine Kenngröße (z. B. Mittelwert oder eine höherwertige Kenngröße) eines Echosignals zusammen mit dem SRV und der maximalen Verlagerung zur Klassifikation verwendet.
  • 1 zeigt ein Verfahren für die Klassifikationsvorverarbeitung zur Verwendung in der medizinischen Ultraschallscherwellenbildgebung. Das Verfahren wird durch das in 7 gezeigte System oder ein anderes System implementiert. Es können zusätzliche, verschiedene oder weniger Vorgänge bereitgestellt werden. Zum Beispiel wird Vorgang 30 nicht durchgeführt, und die Belastungsquelle wird vom Körper, manuell, unter Verwendung eines Stoßspannungsgenerators oder durch einen anderen Mechanismus bereitgestellt. Als weiteres Beispiel werden verschiedene Parameter der Verlagerung statt des SRV in Vorgang 36 und/oder statt der maximalen Verlagerung in Vorgang 38 berechnet. Vorgang 40 ist optional. In noch einem weiteren Beispiel werden Vorgänge 44 und/oder 46 nicht bereitgestellt, und die Klassifikation wird für andere Zwecke verwendet. Die Vorgänge werden in der beschriebenen oder gezeigten Reihenfolge durchgeführt, können aber in anderen Reihenfolgen durchgeführt werden.
  • Die Beschreibung des Verfahrens von 1 wird im Kontext eines in 2 dargestellten Beispiels bereitgestellt. 2 zeigt ein B-Mode- oder Echobild einer zweidimensionalen Brustgeweberegion. Das Verfahren kann für andere Organe verwendet werden. Die Region enthält eine Zyste (dunklere Region um Ort B). Es können sich mehr als eine oder keine Zysten in der Region befinden. Die Region kann andere Materialtypen einschließen, wie etwa Knochen oder Flüssigkeit. Es werden zwei repräsentative Orte A und B gezeigt. Ort A liegt innerhalb eines festen Gewebebereichs, und Ort B liegt innerhalb von flüssigem Gewebe. Das Verfahren wird für einen oder mehr Orte durchgeführt. In einer Ausführungsform wird das Verfahren für jeden der B-Mode- oder Abtastungsorte in einem ganzen Sichtfeld oder einem Untersuchungsbereich durchgeführt. Es kann eine weniger dichte oder eine dichtere Abtastung verwendet werden.
  • In Vorgang 30 von 1 wird eine akustische Anregung in einen Patienten gesendet. Die akustische Anregung wirkt als eine Impulsanregung. Zum Beispiel wird eine 400-Zyklensendewellenform mit Leistungs- oder Spitzenamplitudenpegeln, die ähnlich wie bei B-Mode-Sendungen oder höher sind, für die Bildgebung von Gewebe gesendet. In einer Ausführungsform ist die Sendung eine Strahlungskraftsequenz, die auf das Sichtfeld angewendet wird. Es kann jede ARFI-Sequenz (Schallstrahlungskraftbildgebung, Acoustic Radiation Force Imaging) verwendet werden.
  • Die Sendung wird durch Kraft, Amplitude, Zeitgebung oder ein anderes Merkmal konfiguriert, um eine Belastung auf das Gewebe zu verursachen, die zum Verlagern des Gewebes an einem oder mehr Orten ausreicht. Zum Beispiel wird ein Sendebrennpunkt in der Nähe eines unteren Mittelpunkts des Sichtfelds positioniert, um eine Verlagerung durch das gesamte Sichtfeld zu verursachen. Die Sendung kann für unterschiedliche Teilregionen wiederholt werden.
  • Die Anregung wird von einem Ultraschall-Schallkopf gesendet. Die Anregung ist Schallenergie. Die Schallenergie wird gebündelt, wodurch sich ein dreidimensionales Strahlenprofil ergibt. Die Anregung wird unter Verwendung einer phasengesteuerten Anordnung und/oder eines mechanischen Brennpunkts gebündelt. Die Anregung kann in einer Dimension ungebündelt sein, wie etwa der Höhendimension. Die Anregung wird in das Gewebe eines Patienten gesendet.
  • In Vorgang 32 wird ein Verlagerungsprofil der Antwort in dem Patienten ermittelt. Zum Beispiel sind die Verlagerungsprofile für die Orte A und B in 3 gezeigt. Die Anregung verursacht die Verlagerung des Gewebes. Eine Scherwelle wird erzeugt und pflanzt sich von der Brennregion aus fort. Während sich die Scherwelle durch das Gewebe bewegt, wird das Gewebe verlagert. Es können Längswellen oder andere Verlagerungsursachen verwendet werden. Das Gewebe wird gezwungen, sich in dem Patienten zu bewegen.
  • Die durch die Kraft oder Belastung verursachte Verlagerung wird gemessen. Die Verlagerung wird im Zeitverlauf an einem oder mehr Orten gemessen. Die Messung der Verlagerung kann beginnen, bevor die Belastung oder der Impuls endet, wie etwa unter Verwendung einer verschiedenen Frequenz oder Kodierung. Eine andere Möglichkeit besteht darin, dass die Messung der Verlagerung beginnt, nachdem der Impuls endet. Da die Scher-, Längs- oder andere Welle, welche die Verlagerung im Gewebe auslöst, ausgehend von dem Punkt oder der Region der Belastung Zeit braucht, um sich weiter zu bewegen, kann die Verlagerung von einem entspannten oder teilbelasteten Zustand zu einer maximalen Verlagerung und dann zu einem entspannten Zustand gemessen werden, wie dies in 3 dargestellt ist. Eine andere Möglichkeit besteht darin, dass die Verlagerung nur dann gemessen wird, während sich das Gewebe entspannt, um die maximale Ausdehnung zu bilden.
  • Die Messung betrifft den Grad oder die Größenordnung der Verlagerung. Das Gewebe wird in eine beliebige Richtung bewegt. Die Messung kann entlang der Richtung der größten Bewegung erfolgen. Die Größenordnung des Bewegungsvektors wird ermittelt. Eine andere Möglichkeit besteht darin, dass die Messung entlang einer gegebenen Richtung erfolgt, wie etwa senkrecht zu der Abtastlinie, unabhängig davon, ob das Gewebe mehr oder weniger in andere Richtungen verlagert wird.
  • Die Verlagerung wird mittels Ultraschallabtastung erkannt. Eine Region, wie etwa ein Untersuchungsbereich, das gesamte Sichtfeld oder ein Teiluntersuchungsbereich wird mittels Ultraschall abgetastet. Über einen gegebenen Zeitraum wird der Ultraschall zu dem Gewebe oder Untersuchungsbereich gesendet. Es kann jede derzeit bekannte oder später entwickelte Verlagerungsbildgebung verwendet werden. Zum Beispiel werden Impulse mit einer Dauer von 1–5 Zyklen mit einer Intensität von weniger als 720 mW/cm2 verwendet. Es können auch Impulse mit anderen Intensitäten verwendet werden.
  • Echos oder Reflexionen der Sendung werden empfangen. Die Echos sind strahlengeformt, und die strahlengeformten Daten repräsentieren einen oder mehr Orte. Um die Verlagerung zu erkennen, wird Ultraschallenergie zu dem Gewebe gesendet, das die Verlagerung durchläuft, und Reflexionen der Energie werden empfangen. Jede Sende- und Empfangssequenz kann verwendet werden.
  • Durch mehrmaliges Durchführen des Sendens und Empfangens werden Daten empfangen, die eine ein-, zwei- oder dreidimensionale Region zu verschiedenen Zeitpunkten darstellen. Das Senden und Empfangen wird mehrmals durchgeführt, um die durch die Verlagerung bedingte Veränderung zu erkennen. Durch wiederholtes Abtasten mittels Ultraschall wird die Position des Gewebes zu verschiedenen Zeitpunkten ermittelt.
  • Die Echos werden unter Verwendung von B-Mode- oder Doppler-Erkennung erkannt. Die Verlagerung wird aus den Differenzen für jeden räumlichen Ort erkannt. Zum Beispiel wird die Geschwindigkeit, Varianz, Veränderung des Intensitätsmusters (z. B. Speckle-Tracking) oder eine andere Information aus den empfangenen Daten als die Verlagerung erkannt.
  • In einer Ausführungsform unter Verwendung von B-Mode-Daten werden die Daten aus verschiedenen Abtastungen korreliert. Zum Beispiel wird ein aktueller Satz von Daten mehrmals mit einem Bezugsdatensatz korreliert. Es werden unterschiedliche relative Übersetzungen und/oder Rotationen zwischen den zwei Datensätzen durchgeführt. Der Ort eines Teildatensatzes, der an einem gegebenen Ort im Bezugssatz zentriert ist, wird in dem aktuellen Satz festgestellt.
  • Der Bezugspunkt ist ein erster Datensatz oder Daten von einer anderen Abtastung. Der gleiche Bezugspunkt wird für die gesamte Verlagerungserkennung verwendet, oder die Bezugsdaten verändern sich in einem laufenden oder sich bewegenden Fenster.
  • Die Korrelation ist ein-, zwei- oder dreidimensional. Zum Beispiel wird Korrelation entlang einer Abtastlinie vom Schallkopf weg und zum Schallkopf hin verwendet. Für eine zweidimensionale Abtastung erfolgt die Übersetzung an zwei Achsen entlang mit oder ohne Rotation. Für eine dreidimensionale Abtastung erfolgt die Übersetzung an drei Achsen entlang mit oder ohne Rotation um drei oder weniger Achsen. Der Grad der Ähnlichkeit oder Korrelation der Daten an jeder der verschiedenen Versetzungspositionen wird berechnet. Die Übersetzung und/oder Rotation mit einer größten Korrelation stellt den Bewegungsvektor oder die Versetzung für den Zeitpunkt dar, der mit den aktuellen Daten verbunden ist, die mit dem Bezugspunkt verglichen werden.
  • Es kann jede derzeit bekannte oder später entwickelte Korrelation verwendet werden, wie etwa Kreuzkorrelation, Mustervergleich oder Mindestsumme absoluter Differenzen. Gewebestruktur und/oder Speckle werden korreliert. Unter Verwendung der Doppler-Erkennung übergibt ein Störzeichenfilter Informationen, die in Zusammenhang mit sich bewegendem Gewebe stehen. Die Geschwindigkeit des Gewebe wird aus mehreren Echos abgeleitet. Die Geschwindigkeit wird zum Erkennen einer Verlagerung zum Schallkopf hin oder vom Schallkopf weg verwendet. Alternativ können die Relativgeschwindigkeiten oder die Differenz zwischen Geschwindigkeiten an verschiedenen Orten auf eine Belastung oder Verlagerung hinweisen.
  • 3 zeigt zwei beispielhafte Verlagerungsprofile. Es wird die Größenordnung im Abstand des Bewegungsvektors im Zeitverlauf von den Bezugsdaten gezeigt. Der Analysezeitraum erstreckt sich über etwa 10 Millisekunden, kann aber länger oder kürzer sein. Das Verlagerungsprofil für den in Zusammenhang mit festem Gewebe stehenden Ort A hat ein glatteres Erscheinungsbild als das Verlagerungsprofil für den Ort B, der in Zusammenhang mit dem flüssigen Gewebe steht. Andere Verlagerungsprofile sind möglich.
  • In Vorgang 34 wird ein Merkmal des Verlagerungsprofils berechnet. Es kann ein beliebiges Merkmal verwendet werden. Es kann mehr als ein Merkmal berechnet werden. In dem in 1 gezeigten Beispiel werden sowohl das SRV des Verlagerungsprofils (Vorgang 36) als auch die maximale Verlagerung (Vorgang 38) berechnet. Es können unterschiedliche, zusätzliche oder weniger Merkmale berechnet werden.
  • In Vorgang 36 wird das SRV des Verlagerungsprofils ermittelt. Das SRV der Verlagerung im Zeitverlauf oder als Funktion der Zeit wird ermittelt, indem das Rauschen aus dem Profil festgestellt wird. Das Rauschen kann auf beliebige Weise festgestellt werden, wie etwa durch Auswahl hochfrequenter Komponenten des Profils. Zum Ermitteln des hochfrequenten Rauschens kann eine Fourier-Transformation verwendet werden.
  • In einer Ausführungsform wird das Verlagerungsprofil zur Feststellung des Rauschens gefiltert. Ein Tiefpassfilter, wie etwa ein Butterworth-Filter, wird auf das Profil angewendet (d. h. zeitliche Filterung). Das Filter ist ein Filter mit unendlicher Impulsantwort (IIR) oder ein Filter mit endlicher Impulsantwort (FIR). Die gefilterte Verlagerung im Zeitverlauf wird von der Verlagerung im Zeitverlauf vor der Filterung subtrahiert. Die Differenz stellt das Rauschen dar.
  • 4 zeigt die Verlagerungsprofile für Orte A und B vor und nach der Filterung (linke Seite). Die ungefilterten Verlagerungskurven sind ebenfalls in 4 dargestellt, nämlich als die gestrichelten Linien auf der linken Seite. In dem gefilterten Verlagerungsprofil werden die Hochfrequenzinformationen entfernt oder reduziert. Außerdem ist das Rauschen gezeigt, das sich ergibt, wenn man die gefilterten Verlagerungskurven von den ungefilterten Verlagerungskurven subtrahiert (rechte Seite). Das Rauschen ist für Ort B, der in Zusammenhang mit flüssigem Gewebe steht, variabler und hat eine größere Amplitude als für Ort A, der in Zusammenhang mit festem Gewebe steht.
  • Zum Berechnen des SRV werden die Rauschinformationen quantifiziert. Das quadratische Mittel des Rauschsignals wird berechnet, um den Rauschpegel darzustellen. Es können auch andere Berechnungen verwendet werden, wie etwa ein Durchschnitt der Absolutwerte der Spitzen.
  • Zum Ermitteln das SRV wird der Signalpegel berechnet. In einer Ausführungsform wird der Signalpegel durch die Fläche unter dem gefilterten Verlagerungsprofil dargestellt. Das Integral der gefilterten Verlagerung wird berechnet. Es können auch andere Signalpegelmessungen verwendet werden.
  • Das SRV wird bereitgestellt, indem das Signal (z. B. Integral der gefilterten Verlagerung) durch das Rauschen (z. B. quadratisches Mittel des Rauschens) dividiert wird. Es können auch andere Funktionen verwendet werden, u. a. auch andere Variablen. In dem in 4 gezeigten Beispiel wird das SRV für zwei Orte A, B berechnet. Das SRV kann auch für andere Orte berechnet werden.
  • In Vorgang 38 wird die maximale Verlagerung berechnet. Die maximale Verlagerung wird aus dem Verlagerungsprofil berechnet. Der Spitzenwert oder höchste Grad der Bewegung oder Größenordnung der Verlagerung durch das Gewebe entlang einer Linie, innerhalb einer Ebene oder innerhalb eines Volumens wird berechnet. Die geglättete oder gefilterte Verlagerungskurve wird zum Berechnen des Höchstwerts verwendet. In anderen Ausführungsformen kann die rohe oder ungefilterte Verlagerungskurve verwendet werden. Der Maximalwert über das gesamte Profil oder einen Teil des Profils wird festgestellt oder ermittelt. In dem in 4 gezeigten Beispiel tritt die maximale Verlagerung von 0,9 Mikrometer bei etwa 1,7 Millisekunden für Ort A auf, und für Ort B tritt die maximale Verlagerung von 1,3 Mikrometer bei etwa 0,8 Millisekunden auf.
  • In dem optionalen Vorgang 40 wird ein anderer Parameter, also nicht aus dem Verlagerungsprofil, berechnet. In dem in 1 gezeigten Beispiel können zu den anderen Parametern u. a. eine Kenngröße der Echoinformationen, das Signal-Rausch-Verhältnis des Echosignals und/oder der Dekorrelationsmesswert des Echosignals im Zeitverlauf zählen. Zum Beispiel wird der Mittelwert der B-Mode-Informationen für einen Ort über den gleichen Zeitraum, der für das Verlagerungsprofil verwendet wurde, berechnet. Höherwertige Kenngrößen können stattdessen oder zusätzlich berechnet werden. Es kann jede Kenngröße oder jede andere Parameter für die Klassifikation verwendet werden. Die Kenngröße kann räumlich und/oder zeitlich sein. Das SRV des Echosignals kann als Verhältnis zwischen der Empfangssignalintensität mit oder ohne Senden berechnet werden. Der Dekorrelationskoeffizient kann als der 1 minus den normalisierten Korrelationskoeffizienten der Echosignale berechnet werden, die vor und nach dem ARFI- oder einem anderen Anregungsstrahl erkannt werden.
  • In Vorgang 42 werden ein oder mehr Orte klassifiziert. Das Gewebe oder sonstige Material an jedem Ort wird unabhängig oder separat klassifiziert. Die Klassifikation für jeden Ort hängt von den Daten für diesen Ort ab und nicht von den Daten für andere Orte. In anderen Ausführungsformen können räumliche Filterung oder Informationen aus benachbarten Orten verwendet werden, um einen gegebenen Ort zu klassifizieren.
  • Die Klassifikation stellt das Material fest. Es kann jeder beliebige Materialtyp festgestellt werden. In einer Ausführungsform stellt die Klassifikation den Ort als Flüssigkeit (eine Kategorie oder ein Typ, der flüssiges Gewebe einschließt), festes Gewebe oder sonstiges (z. B. nicht ermittelbar) fest. Der Ort wird als einer der Gruppe klassifiziert. Die Gruppe kann mehr oder weniger Optionen oder Klassen umfassen, wie etwa Auswahl zwischen (a) Flüssigkeit oder flüssigem Gewebe und (b) festem Gewebe.
  • Das Klassifizieren unterscheidet festes Gewebe von anderen Materialien, wie etwa von flüssigem Gewebe und einer anderen Klasse. Ein Prozessor führt das Klassifizieren ohne Benutzerauswahl oder -eingabe durch. Der Prozessor verwendet das Merkmal des Verlagerungsprofils mit oder ohne andere Informationen oder Parameter. Zum Beispiel klassifiziert der Prozessor als eine Funktion des SRV, der maximalen Verlagerung oder des SRV und der maximalen Verlagerung. Das Klassifizieren mit dem SRV und der maximalen Verlagerung kann zwischen flüssigem Gewebe, festem Gewebe und sonstigem oder nicht ermitteltem unterscheiden.
  • 5 zeigt eine Ausführungsform der Klassifikation unter Verwendung von zwei Parametern, nämlich SRV und maximale Verlagerung, des Verlagerungsprofils. Der Flackerrauschpegel wird zur Einstufung einer Kategorie ,nicht ermittelt' verwendet. Ist die maximale Verlagerung gering, liegen vielleicht unzureichende Informationen oder Nicht-Gewebe- oder flüssiges Material vor. Bei geringer maximaler Verlagerung erfolgt keine Klassifikation, oder die Klasse ist die Kategorie ,nicht ermittelt'. Bei geringerem SRV und höherer Verlagerung wird das Material an dem Ort als Flüssigkeit oder flüssiges Gewebe klassifiziert. Bei anderen Kombinationen der zwei Parameter wird das Material an dem Ort als festes Gewebe klassifiziert. Es können auch andere als die in 5 gezeigten Funktionen verwendet werden.
  • Die Klassifikationsfunktion basiert auf Experimenten, einer mathematischen Funktion, statistischen Kenngrößen oder anderen Informationen. Zum Beispiel wird die Fuzzylogik zum Klassifizieren verwendet. In diesem Verfahren wird für jeden Parameter (z. B. SRV, maximale Verlagerung) eine Mitgliedschaftsfunktion jedes Gewebetyps empirisch definiert, und die Ausgabewerte jeder Mitgliedschaftsfunktion je nach den eingegebenen Parametern an einem Ort werden summiert, um die Wahrscheinlichkeit eines Gewebetyps zu erzeugen. Der Maximalwert der Summe entspricht dem zuzuordnenden Gewebetyp. Als anderes Beispiel wird maschinelles Lernen aus einer Sammlung von Schulungsmustern oder Daten mit bekannter Grundwahrheit hinsichtlich der Klasse verwendet, um Kenngrößen oder eine Matrixfunktion für das Klassifizieren zu ermitteln. Eine Wahrscheinlichkeitsfunktion zeigt an, wie wahrscheinlich jede Klasse für einen gegebenen Ort ist, und die Klasse mit der höchsten Wahrscheinlichkeit wird gewählt. Die Klassifikation wird unter Verwendung einer Nachschlagetabelle, Fuzzylogikfunktion, programmierten Funktion oder Matrixfunktion durchgeführt.
  • Wenn man nun nochmals 1 heranzieht, stellt die Rückmeldung von Vorgang 42 zu Vorgang 32 eine Wiederholung der Klassifizierung für eine Vielzahl von Orten dar. Zum Beispiel werden das Verlagerungsprofil und die Merkmale für das Verlagerungsprofil für jeden Ort in einem Sichtfeld ermittelt. Die Wiederholung verwendet die gleiche oder eine andere Sendung in Vorgang 30. Wenn der Untersuchungsbereich klein genug ist, wird ein Impuls verwendet. Die Verlagerung an verschiedenen Orten wird mit einem auf jeden der Orte zentrierten Fenster ermittelt. Für jeden Ort wird das Fenster oder der Kern über dem Ort zentriert. Die Daten, die räumliche Orte innerhalb des Fensters darstellen, werden zur Korrelation verwendet. Die' Verlagerung wird für jeden Ort separat ermittelt. In anderen Ausführungsformen wird der Sendevorgang 30 wiederholt. Verlagerungen für einen, einige oder einen Teilsatz eines Untersuchungsbereichs werden als Antwort auf jede Sendung von Vorgang 30 ermittelt.
  • Die Verlagerungsprofile werden über eine Region von beliebiger Größe erkannt. In einer Ausführungsform wird die Verlagerung in einem Untersuchungsbereich erkannt, der das zu diagnostizierende Gewebe wahrscheinlich einschließt, wie etwa ungefähr ein Drittel bis die Hälfte der vollständigen Abtastungsregion für B-Mode-Bildgebung. Es kann ein größerer, kleinerer oder kein Untersuchungsbereich verwendet werden, wie etwa beim Erkennen der Verlagerung über die gesamte Bildgebungsregion. Engere Untersuchungsbereiche ermöglichen es unter Umständen, die Verlagerung mit weniger Wiederholungen der Sendung der Anregungswellenform von Vorgang 30 zu erkennen. Je nach der Zahl der Empfangsstrahlen, die geformt werden können, und der Probendichte können entweder keine, eine oder mehrere Wiederholungen verwendet werden.
  • Eine vollständige Abtastung, wie etwa Abtastung der Verlagerung an jedem B-Mode-Abtastungsort, kann verwendet werden. Es kann auch eine größere oder kleinere (z. B. weniger dichte) Abtastung der Verlagerung in Bezug auf das B-Mode- Abtastraster verwendet werden. 6 zeigt ein B-Mode-Bild und ein SRV-Bild für dieselbe Abtastdichte. Die SRV-Informationen liefern einen besseren Hinweis auf die Zyste. Das SRV-Bild kann dem Benutzer angezeigt werden oder nicht.
  • Die Klassifikation wird zum räumlichen Ausmaskieren von Scherwelleninformationen in Vorgang 44 verwendet. Scherwelleninformationen sollen für Orte mit festem Gewebe angezeigt werden. Für andere Orte, wie etwa für flüssiges Gewebe oder nicht ermittelte Orte, werden die Scherwelleninformationen nicht angezeigt. Die Orte werden maskiert. Auf der Grundlage der Maskierung werden ermittelte Scherwelleninformationen entweder nicht angezeigt, oder die Scherwelleninformationen werden nicht einmal berechnet. Die unzuverlässigen Scherungsinformationen aus Flüssigkeit oder flüssigem Gewebe werden dem Benutzer nicht in dem Bild dargestellt, wodurch Verwirrung vermieden wird.
  • In Vorgang 46 wird die Scherwellenbildgebung durchgeführt. Schergeschwindigkeit, Schermodul oder andere aus der Reaktion des Gewebes auf eine Scherwelle ermittelte Informationen werden erkannt. Es kann eine beliebige Scherungsbildgebung verwendet werden. Aufgrund der Maskierung von Vorgang 44 stellt das angezeigte Bild die Scherwelleninformationen für feste Geweberegionen und nicht für Flüssigkeit oder flüssige Geweberegionen dar.
  • Die Scherungsinformationen werden für die gleichen Orte erkannt wie die Verlagerung. Es werden andere Sendungen und Abtastungen für die Scherungsbildgebung verwendet als für die Berechnung der Verlagerung. Eine andere Möglichkeit besteht darin, die gleichen Sendungen und sogar Verlagerungsprofile für die Scherungsbildgebung zu verwenden wie für die Klassifizierung. Die Zeit, die die Scherwelle braucht, um sich vom Ursprung (z. B. Sendebrennregion) zu dem Ort zu bewegen, wird ermittelt. Die maximale Verlagerung oder ein anderer Teil des Verlagerungsprofils gibt die Ankunftszeit der Scherwelle an. Die Geschwindigkeit der Scherwelle wird aus den Zeitgebungsinformationen berechnet.
  • Für die Scherwellenbildgebung erzeugt eine Impulsanregung eine Scherwelle an einem räumlichen Ort. Wo die Anregung stark genug ist, wird eine Scherwelle erzeugt. Die Scherwelle pflanzt sich durch Gewebe langsamer fort als die Längswelle entlang der Richtung der akustischen Wellenaussendung. Die Scherwelle pflanzt sich in verschiedenen Richtungen fort, u. a. in einer Richtung, die senkrecht zu der Richtung der angewendeten Belastung steht. Die Verlagerung der Scherwellen ist an den Orten größer, die näher an dem Ort sind, an dem die Scherwelle erzeugt wird.
  • Ultraschalldaten werden gewonnen. Zumindest einige der Ultraschalldaten reagieren auf eine Scherwelle. Ein Untersuchungsbereich wird überwacht, um die Scherwelle zu erkennen. Der Untersuchungsbereich hat eine beliebige Größe, wie etwa 6 mm in lateraler und 10 mm in axialer Richtung. Diese Erkennungsregion wird mittels Ultraschall überwacht. Zum Beispiel werden B-Mode-Abtastungen durchgeführt, um die durch die Scherwelle verursachte Gewebeverlagerung zu erkennen. Zur Scherwellenü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 Scherwelle werden wiederholte B-Mode-Sendungen an einer einzelnen Abtastlinie entlang und Empfangsvorgänge an vier benachbarten Abtastlinien entlang durchgeführt. In anderen Ausführungsformen wird nur ein einzelner Empfangsstrahl geformt, oder es wird eine andere Zahl von Empfangsstrahlen als Antwort auf jede Sendung geformt. Es kann eine beliebige Zahl von Wiederholungen verwendet werden, wie etwa ungefähr 120 Mal. Manche der Ultraschalldaten, wie etwa am Beginn oder am Ende der Wiederholungen, reagieren unter Umständen nicht auf die Scherwelle.
  • Während sich die Scherwelle durch die Abtastlinien fortpflanzt, kann die B-Mode-Intensität aufgrund der Verlagerung des Gewebes variieren. Für die überwachten Abtastlinien wird eine Datensequenz bereitgestellt, die ein Zeitprofil der Gewebebewegung darstellt, die sich aus der Scherwelle ergibt. Zum Beispiel werden Daten von einer Vielzahl räumlicher Orte (z. B. an den Abtastlinien entlang) als Funktion der Zeit korreliert. Es kann eine beliebige Elastizitätserkennung verwendet werden. Für jede Tiefe oder jeden räumlichen Ort wird eine Korrelation über eine Vielzahl von Tiefen oder räumlichen Orten durchgeführt (z. B. Kern von 64 Tiefen, wobei die mittlere Tiefe der Punkt ist, für die das Profil berechnet wird). Es kann eine zwei- oder dreidimensionale Verlagerung im Raum verwendet werden. Es kann eine eindimensionale Verlagerung entlang einer anderen Richtung als die der Abtastlinien oder Strahlen verwendet werden.
  • Die räumliche Versetzung mit der höchsten oder ausreichenden Korrelation weist auf Grad und Richtung der Verlagerung hin. Das Ermitteln der Verlagerungen erfolgt für einen bestimmten Ort zu verschiedenen Zeitpunkten. Das zeitliche Profil für einen bestimmten Ort weist auf die Erkennung der Scherwelle hin. Das Profil wird auf Nicht-Rauschen oder einen einzelnen Fall von Variation untersucht. Eine Spitze im Profil, mit oder ohne zeitliches Herausfiltern der niedrigen Frequenzanteile, weist auf das Passieren der Scherwellenfront hin. Die größte Verlagerung wird gewählt, aber es kann auch der Durchschnitt oder eine andere Kenngröße der Verlagerung verwendet werden. Die maximale Scherung an einem bestimmten Ort wird erkannt. Alternativ kann auch eine durchschnittliche oder andere Scherung erkannt werden.
  • Zur Überwachung einer größeren Region werden zusätzliche Empfangsstrahlen als Antwort auf den Überwachungssendestrahl geformt. Eine andere Möglichkeit besteht darin, dass eine andere Scherwelle erzeugt wird und die Sendestrahlen und Empfangsstrahlen in einem anderen Abstand von dem Punkt der Scherwellenerzeugung bereitgestellt werden. In dem 6 mm × 10 mm großen Beispiel oben können 36 Empfangsabtastlinien bereitgestellt werden. Bei vier Empfangsstrahlen pro Sendestrahl wird der Prozess für unterschiedliche laterale Abstände neunmal wiederholt. Für jeden Empfangsstrahlort wird ein Zeitprofil der Bewegungsinformationen, dargestellt durch die Ultraschalldaten, bereitgestellt. Sendungen entlang verschiedener Abtastlinien zur Überwachung derselben Scherwelle werden während der Bildung des zeitlichen Profils vermieden, um eine höhere zeitliche Auflösung bereitzustellen, aber es können verschachtelte oder wechselnde Abtastpositionen bereitgestellt werden.
  • Die obige Beschreibung betrifft eine Tiefe. Die Abtastung kann so angeordnet werden, dass ein ,Gate' bereitgestellt wird, welches das gesamte axiale Ausmaß des Untersuchungsbereichs abdeckt. In einer anderen Ausführungsform werden Abtastpunkte für jeden Empfangsstrahl bei verschiedenen Tiefen erhalten. Ein separates Zeitprofil wird für jede axiale Tiefe sowie für jeden lateralen Ort bereitgestellt. Es kann eine beliebige Zahl von Tiefen verwendet werden, wie etwa ungefähr 200 für 5 mm oder 400 für 10 mm.
  • Ultraschalldaten, die verschiedene Orte in dem Untersuchungsbereich darstellen, werden gewonnen. Die Ultraschalldaten werden in Echtzeit mit der Abtastung gewonnen oder aus einem Speicher gewonnen. Die Bewegungsinformationen stellen für jeden Ort die Antwort zu verschiedenen Zeitpunkten dar. Es können auch andere Abtastungen, Überwachungen oder Techniken verwendet werden, um Ultraschalldaten zum Schätzen der Schergrößenordnung zu gewinnen.
  • Die Schergeschwindigkeit wird für die verschiedenen räumlichen Orte des Gewebes erkannt. Für jeden Ort wird die Verlagerung als Funktion der Zeit ermittelt. Die Schergeschwindigkeit wird gewonnen, indem eine Zeit ab der Erzeugung der Scherwelle bis zur Erkennung der Scherwelle an einem anderen Ort ermittelt wird. Die Zeit und der Abstand zu dem Ort bestimmen die Geschwindigkeit. Der Abstand ist bekannt aus dem Abtastlinienabstand (d. h. der Position des Sendestrahls zum Erzeugen der Scherwelle und der Position des Empfangsstrahls zum Erkennen der Scherwelle). Die Zeit ist bekannt aus der relativen Zeit zwischen Erzeugung und Erkennung der Scherwelle.
  • Es können auch andere Techniken verwendet werden, um die Spitze in dem Profil zu erkennen. Zum Beispiel wird eine Regression angewendet. Da die Scherwellengeschwindigkeit linear ist, kann eine robuste lineare Regression mit automatischer Ausreißererkennung die Scherwellengeschwindigkeit angeben. Die Ultraschalldaten für alle Abtastpunkte im Untersuchungsbereich werden für Abstand als Funktion der Zeit oder nach Zeit und Abstand in einem Diagramm dargestellt. Die lineare Regression wird auf das Diagramm oder die Daten angewendet, wodurch eine an die Daten angepasste Linie bereitgestellt wird. Die Steigung der Linie gibt die Scherwellengeschwindigkeit an.
  • Die Scherwelleninformationen werden für eine Farbüberlagerung oder andere Modulation von Anzeigewerten verwendet. Zum Beispiel werden Scherwelleninformationen über oder mit B-Mode-Informationen angezeigt. Die Verlagerungsdaten sind in einem Anzeigeformat oder können eine in ein Anzeigeformat umgewandelte Abtastung sein. Die Verlagerungsdaten sind Farb- oder Grauskaladaten, können aber auch Daten vor der Abbildung mit Grauskala- oder Farbskala sein. Die Informationen können linear oder nichtlinear zu den Anzeigewerten abgebildet sein.
  • Das Bild stellt die Verlagerungsinformationen, wie etwa Scherung oder Module (z. B. Schermodule), für die verschiedenen Orte dar. Wenn die Werte für alle Rasterpunkte in einem Untersuchungsbereich oder Sichtfeld ermittelt werden, stellen die Pixel (Bildpunkte) der Anzeige die Schergeschwindigkeiten für diese Region dar. Das Anzeigeraster kann sich von dem Abtastraster und/oder dem Raster, für das die Verlagerungen berechnet werden, unterscheiden. Farbe, Helligkeit, Luminanz, Farbton oder ein anderes Merkmal wird als Funktion der Verlagerung moduliert.
  • Das Bild kann andere Daten einschließen. Zum Beispiel werden B-Mode-Daten oder andere Daten, die Gewebe, Flüssigkeit oder Kontrastmittel in der gleichen Region darstellen, eingeschlossen. Die Verlagerungsdaten werden zur Überlagerung von oder zur Kombination mit anderen Daten verwendet. Die anderen Daten unterstützen den Benutzer dabei, den Ort des Strahls im Verhältnis zu dem zu behandelnden Gewebe zu ermitteln.
  • 7 zeigt eine Ausführungsform eines Systems 10 für die Klassifikationsvorverarbeitung in der medizinischen Ultraschallscherwellenbildgebung. 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 Benutzereingabe zur Interaktion des Benutzers mit dem System bereitgestellt.
  • Das System 10 ist ein medizinisches diagnostisches Ultraschallbildgebungssystem. In alternativen Ausführungsformen ist das System 10 ein PC (Personal Computer), ein Computerarbeitsplatz, ein PACS-Arbeitsplatz (Bildarchivierungs- oder Kommunikationssystem) oder eine andere an einem gleichen Ort befindliche oder über ein Netzwerk verteilte Anordnung zur Echtzeit- oder Post-Erfassungs-Bildgebung.
  • 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. Eine Sequenz von Sendestrahlen wird erzeugt, um eine zwei- oder dreidimensionale Region abzutasten. Es können Sektor-, Vector®-, lineare oder andere Abtastformate verwendet werden. Die gleiche Region wird mehrmals abgetastet. Für die Strömungs- oder Doppler-Bildgebung und für die Scherungsbildgebung wird eine Sequenz von Abtastungen verwendet. Bei der Doppler-Bildgebung kann die Sequenz mehrere Strahlen entlang derselben Abtastlinie umfassen, bevor eine benachbarte Abtastlinie abgetastet wird. Bei der Scherungsbildgebung kann eine Abtast- oder Frame-Verschachtelung verwendet werden (d. h. Abtastung der gesamten Region, bevor diese erneut abgetastet wird). In alternativen Ausführungsformen erzeugt der Sendestrahlformer 12 eine ebene Welle oder eine divergierende Welle für eine schnellere Abtastung.
  • Der gleiche Sendestrahlformer 12 erzeugt Impulsanregungen oder elektrische Wellenformen zur Erzeugung von Schallenergie, um die Verlagerung auszulösen. In alternativen Ausführungsformen wird ein unterschiedlicher Sendestrahlformer zum Erzeugen der Impulsanregung bereitgestellt. Der Sendestrahlformer 12 veranlasst den Schallkopf 14 zur Erzeugung von hochintensiven gebündelten Ultraschallwellenformen.
  • Der Schallkopf 14 ist eine Anordnung zum Erzeugen von Schallenergie aus elektrischen Wellenformen. Bei einer Anordnung bündeln die relativen Verzögerungen die Schallenergie. Ein gegebenes Sendeereignis entspricht der Sendung von Schallenergie durch verschiedene Elemente zu einem im Wesentlichen gleichen Zeitpunkt in Anbetracht der Verzögerungen. Das Sendeereignis stellt einen Ultraschallenergieimpuls zum Verlagern des Gewebes bereit. Der Impuls ist eine Impulsanregung. Die Impulsanregung schließt Wellenformen mit vielen Zyklen (z. B. 500 Zyklen) ein, die aber in einem relativ kurzen Zeitraum auftreten, um die Gewebeverlagerung über einen längeren Zeitraum auszulösen.
  • Der Schallkopf 14 ist eine 1-, 1,25-, 1,5-, 1,75- oder 2-dimensionale Anordnung von piezoelektrischen oder kapazitiven Membranelementen. Der Schallkopf 14 umfasst eine Vielzahl von Elementen für das Umformen zwischen akustischen und elektrischen Energien. Die Empfangssignale werden als Antwort auf die Ultraschallenergie (Echos) erzeugt, die auf die Elemente des Schallkopfs 14 auftrifft. Die Elemente stehen in Verbindung mit Kanälen des Sende- und Empfangsstrahlformers 12, 16. Alternativ wird ein einzelnes Element mit einem mechanischen Brennpunkt verwendet.
  • 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 ist durch Hardware oder Software dazu konfiguriert, relative Verzögerungen, Phasen und/oder Apodisierung anzuwenden, um einen oder mehr Empfangsstrahlen als Antwort auf jede Bildgebungssendung zu formen. Empfangsbetrieb tritt unter Umständen nicht für Echos aus der Impulsanregung ein, die zum Verlagern von Gewebe verwendet wird. Der Empfangsstrahlformer 16 gibt Daten aus, die unter Verwendung der Empfangsssignale 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 werden kombiniert, um Informationen auf einem gewünschten Frequenzband zu isolieren, wie etwa einem zweiten harmonischen, kubischen Grundfrequenz- oder anderem Band.
  • In Koordination mit dem Sendestrahlformer 12 erzeugt der Empfangsstrahlformer 16 Daten, welche die Region zu verschiedenen Zeitpunkten darstellen. Nach der akustischen Impulsanregung erzeugt der Empfangsstrahlformer 16 Strahlen, die verschiedene Linien oder Orte zu verschiedenen Zeitpunkten darstellen. Durch das Abtasten des Untersuchungsbereichs mittels Ultraschall werden Daten (z. B. strahlengeformte Abtastpunkte) erzeugt.
  • 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 Verlagerung. Eine andere Möglichkeit besteht darin, dass die B-Mode-Daten ebenfalls zum Ermitteln der Verlagerung verwendet werden. Als weiteres Beispiel werden Daten für auf Verlagerung basierte Klassifikation und Scherungsbildgebung mit einer Reihe von gemeinsam verwendeten Abtastungen durchgeführt, und B-Mode- oder Doppler-Abtastung wird separat oder unter Verwendung einiger der gleichen Daten durchgeführt.
  • 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 zur Anzeige 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 Verlagerung und Berechnen von Verlagerungsprofileigenschaften. Zum Beispiel ist der separate Prozessor durch Software und/oder Hardware zum Durchführen einer beliebigen Kombination einer oder mehrerer der in 1 gezeigten Vorgänge konfiguriert.
  • Der Prozessor 18 ist zum Schätzen der Gewebeverlagerung konfiguriert, die durch die akustische Impulsanregung induziert wird. Mittels Korrelation, Nachverfolgung, Bewegungsdetektion oder anderen Messungen der Verlagerung wird der Grad der Positionsverschiebung des Gewebes geschätzt. Die Schätzung wird im Verlauf eines Zeitraums mehrmals durchgeführt, wie etwa von vor der durch den Impuls ausgelösten Bewegung des Gewebes bis nach dem Zeitpunkt, an dem das Gewebe zum größten Teil oder vollständig wieder zu einem entspannten Zustand zurückgekehrt ist (z. B. bis das Gewebe sich von der Belastung erholt hat, die durch die Impulsanregung verursacht wurde).
  • Der Prozessor 18 ist dazu konfiguriert, mindestens einen Parameter abzuleiten, der ein Merkmal eines Profils der Gewebeverlagerung beschreibt. Zum Beispiel wird das Signal-Rausch-Verhältnis des Verlagerungsprofils abgeleitet. Als ein weiteres Beispiel wird die maximale Verlagerung des Verlagerungsprofils abgeleitet. Der Prozessor 18 kann andere Parameter berechnen, wie etwa eine Kenngröße der Daten in Zeit, Raum oder Zeit und Raum. Zum Beispiel wird ein durchschnittlicher B-Mode- oder akustischer Impedanzwert im Zeitverlauf und/oder im Raum für jeden Ort berechnet.
  • Der Prozessor 18 ist zum Klassifizieren von Gewebe in der Region konfiguriert. Gewebe kann flüssiges Gewebe oder festes Gewebe sein. Die Klassifikation kann zwischen Typen von Gewebe, zwischen Gewebe und Flüssigkeit oder zwischen Typen von Gewebe und einer oder mehr anderen Klassen sein. In einer Ausführungsform klassifiziert der Prozessor 18 das Material an einem Ort als Flüssigkeit/flüssiges Gewebe, festes Gewebe und sonstiges/nicht ermittelbar.
  • Die Klassifikation basiert auf mindestens einem Parameter. Zum Beispiel basiert die Klassifikation auf Werten für ein oder mehr Merkmale des Verlagerungsprofils. Das Signal-Rausch-Verhältnis und die maximale Verlagerung sind zwei solche Merkmale. Es können auch andere Informationen, wie etwa eine Kenngröße anderer Daten, in der Klassifikation verwendet werden.
  • Der Prozessor 18 implementiert Fuzzylogik, eine Wahrscheinlichkeitsfunktion, Nachschlagetabelle oder einen anderen Prozess. Die eingegebenen Merkmale (z. B. Merkmale des Verlagerungsprofils) werden auf den Prozess angewendet, um eine Kategorie zu ermitteln, mit der man den Ort bezeichnen sollte.
  • Der Prozessor 18 wird entsprechend den Anweisungen betrieben, die im Speicher 22 oder einem anderen Speicher für die Klassifikationsvorverarbeitung in der medizinischen Ultraschallscherwellenbildgebung gespeichert sind. Der Prozessor 18 ist dazu programmiert, Rückmeldungen für hochintensiven gebündelten Ultraschall bereitzustellen.
  • 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 ist vom Prozessor 18 oder einer anderen Vorrichtung durch Eingabe der Signale, die als ein Bild angezeigt werden sollen, konfiguriert. Die Anzeige 20 zeigt ein Bild an, das die Scherung für verschiedene Orte in einem Untersuchungsbereich oder einem gesamten Bild darstellt. Die Klassifikation wird zum Segmentieren der Orte verwendet. Orte, die in Zusammenhang mit festem Gewebe stehen, werden moduliert oder schließen Scherungsinformationen ein. Andere Orte werden maskiert, so dass sie keine Scherungsinformationen einschließen. Die anderen Orte können stattdessen einer größeren Filterung der Scherungsinformationen unterworfen sein.
  • 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 (21)

  1. Verfahren für die Klassifikationsvorverarbeitung zur Verwendung in der medizinischen Ultraschallscherwellenbildgebung, wobei das Verfahren umfasst: Messen (32), mittels Ultraschall, der Verlagerung im Zeitverlauf an einem ersten Ort in einem Patienten als Antwort auf eine Impulsanregung; Ermitteln (38) eines Maximalwerts der Verlagerung im Zeitverlauf; Ermitteln (36) eines Signal-Rausch-Verhältnisses der Verlagerung im Zeitverlauf; und Klassifizieren (42), mit einem Prozessor (18), des ersten Orts als ein erster einer Vielzahl von Gewebetypen, wobei das Klassifizieren (42) eine Funktion des Maximalwerts der Verlagerung im Zeitverlauf und des Signal-Rausch-Verhältnisses der Verlagerung im Zeitverlauf ist.
  2. Verfahren nach Anspruch 1, das weiterhin umfasst: Senden (30) einer akustischen Anregung in einen Patienten, wobei die Impulsanregung die akustische Anregung umfasst; wobei das Messen (32) der Verlagerung das wiederholte Abtasten des ersten Orts mittels Ultraschall umfasst.
  3. Verfahren nach Anspruch 1, worin das Messen (32) der Verlagerung das Senden (30) von Ultraschall zu dem Gewebe und Empfangen von Reflexionen aus dem Senden (30), wobei das Senden (30) des Ultraschalls und das Empfangen mehrere Male durchgeführt werden, und das Erkennen der Verlagerung aus den Reflexionen von dem mehrmaligen Empfangen umfasst.
  4. Verfahren nach Anspruch 1, das weiterhin umfasst: Wiederholen des Klassifizierens (42) für eine Vielzahl von Orten; Durchführen der Scherwellenbildgebung (46) des Patienten einschließlich des ersten Orts und der Vielzahl von Orten; und räumliches Maskieren (44) der Scherwellenbildgebung als eine Funktion des Klassifizierens (42) für den ersten Ort und die Vielzahl von Orten.
  5. Verfahren nach Anspruch 1, wobei das Messen (32) der Verlagerung im Zeitverlauf das Messen (32) der Verlagerung, während sich das Gewebe nach der Impulsanregung wieder entspannt, umfasst.
  6. Verfahren nach Anspruch 1, wobei das Ermitteln (36) des Signal-Rausch-Verhältnisses umfasst: Filtern der Verlagerung im Zeitverlauf; Subtrahieren der gefilterten Verlagerung im Zeitverlauf von der Verlagerung im Zeitverlauf; Berechnen eines quadratischen Mittels eines Resultats des Subtrahierens; und Dividieren eines Integrals der gefilterten Verlagerung durch das quadratische Mittel, wobei ein Resultat des Dividierens das Signal-Rausch-Verhältnis umfasst.
  7. Verfahren nach Anspruch 1, wobei das Klassifizieren (42) das Klassifizieren (42) des ersten als einer von einer Gruppe umfasst, die aus Flüssigkeit, festem Gewebe, und nicht ermittelbar besteht.
  8. Verfahren nach Anspruch 1, wobei das Klassifizieren (42) das Klassifizieren (42) des ersten als flüssiges Gewebe oder festes Gewebe umfasst.
  9. Verfahren nach Anspruch 1, wobei das Klassifizieren (42) das Klassifizieren (42) mittels Fuzzylogik oder einer Wahrscheinlichkeitsfunktion umfasst.
  10. Verfahren nach Anspruch 1, das weiterhin umfasst Ermitteln eines Signal-Rausch-Verhältnisses eines Ultraschallsignals; und Ermitteln eines Dekorrelationskoeffizienten eines Ultraschallsignals im Zeitverlauf; wobei das Klassifizieren (42) das Klassifizieren (42) als eine Funktion des Maximalwerts der Verlagerung im Zeitverlauf, des Signal-Rausch-Verhältnisses der Verlagerung im Zeitverlauf, des Signal-Rausch-Verhältnisses des Ultraschallsignals und des Dekorrelationskoeffizienten des Ultraschallsignals im Zeitverlauf umfasst.
  11. In einem nichtflüchtigen, computerlesbaren Speichermedium (22) mit darauf gespeicherten Daten, die durch einen programmierten Prozessor (18) ausführbare Anweisungen für die Klassifikationsvorverarbeitung in der medizinischen Ultraschallscherwellenbildgebung darstellen, wobei das Speichermedium (22) Anweisungen umfasst zum: Ermitteln (32) eines Verlagerungsprofils der Antwort in einem Patienten; Berechnen (34) eines ersten Merkmals des Verlagerungsprofils; Unterscheiden (42) von Flüssigkeit und flüssigem Gewebe von festem Gewebe, welches zumindest teilweise auf dem ersten Merkmal basiert; und Ausmaskieren (44) von Scherwelleninformationen aus einem Bild für Flüssigkeit oder flüssiges Gewebe.
  12. Nichtflüchtiges, computerlesbares Speichermedium (22) nach Anspruch 11, wobei das Ermitteln (32) des Verlagerungsprofils das Empfangen von Daten, die eine Region des Patienten zu verschiedenen Zeitpunkten darstellen, und das Korrelieren der Daten umfasst.
  13. Nichtflüchtiges, computerlesbares Speichermedium (22) nach Anspruch 11, wobei das Ermitteln (32) des Verlagerungsprofils das Ermitteln einer Größenordnung der Verlagerung als Funktion der Zeit nach einem Impuls, der die Verlagerung von Gewebe in dem Patienten verursacht, umfasst.
  14. Nichtflüchtiges, computerlesbares Speichermedium (22) nach Anspruch 11, wobei das Berechnen (34) das Berechnen (36) eines Signal-Rausch-Verhältnisses des Verlagerungsprofils umfasst.
  15. Nichtflüchtiges, computerlesbares Speichermedium (22) nach Anspruch 11, wobei das Berechnen (34) das Berechnen (38) einer maximalen Verlagerung aus dem Verlagerungsprofil umfasst.
  16. Nichtflüchtiges, computerlesbares Speichermedium (22) nach Anspruch 11, wobei das Unterscheiden (42) das Klassifizieren (42) als eine Funktion des ersten Merkmals umfasst.
  17. Nichtflüchtiges, computerlesbares Speichermedium (22) nach Anspruch 11, wobei das Unterscheiden (42) das Klassifizieren (42) zwischen (a) der Flüssigkeit und dem flüssigen Gewebe, (b) festem Gewebe und (c) einer anderen Klasse umfasst.
  18. Nichtflüchtiges, computerlesbares Speichermedium (22) nach Anspruch 11, wobei das Maskieren (44) das Anzeigen (46) des Bilds umfasst, das Scherwelleninformationen für feste Geweberegionen und nicht für Flüssigkeit oder flüssige Geweberegionen darstellt.
  19. System für die Klassifikationsvorverarbeitung in der medizinischen Ultraschallscherwellenbildgebung, wobei das System umfasst: Schallkopf (14), der zum Senden einer akustischen Impulsanregung in einen Patienten konfiguriert ist und der zum Abtasten einer Region des Patienten mittels Ultraschall konfiguriert ist; Empfangsstrahlformer (16), der zum Erzeugen von Daten konfiguriert ist, die die Region zu verschiedenen Zeitpunkten nach der akustischen Impulsanregung darstellen, wobei die Daten aus der Abtastung mittels Ultraschall erzeugt werden; Prozessor (18), der dazu konfiguriert ist, Gewebeverlagerung, die durch die akustische Impulsanregung verursacht ist, zu schätzen, zumindest einen Parameter, der ein Merkmal eines Profils der Gewebeverlagerung beschreibt, abzuleiten und Gewebe in der Region auf der Grundlage von zumindest dem einen Parameter zu klassifizieren; und Anzeige (20), die zum Anzeigen eines Bilds konfiguriert ist, das die Scherung für erste Orte darstellt, wobei die ersten Orte als eine Funktion der Klassifikation des Gewebes ermittelt werden.
  20. System nach Anspruch 19, wobei der Prozessor (18) zum Ableiten des einen Parameters als ein Signal-Rausch-Verhältnis konfiguriert ist, zum Ableiten eines anderen Parameters als eine maximale Verlagerung des Profils konfiguriert ist und zum Klassifizieren auf Grundlage des Signal-Rausch-Verhältnisses und der maximalen Verlagerung konfiguriert ist.
  21. System nach Anspruch 19, wobei der Prozessor (18) zum Ableiten des einen Parameters als ein Signal-Rausch-Verhältnis konfiguriert ist, zum Ableiten eines anderen Parameters als eine maximale Verlagerung des Profils konfiguriert ist, zum Berechnen einer Kenngröße der Daten in Zeit, Raum oder in Zeit und Raum konfiguriert ist und zum Klassifizieren auf der Grundlage der Kenngröße, des Signal-Rausch-Verhältnisses und der maximalen Verlagerung konfiguriert ist.
DE102012108353.3A 2011-09-09 2012-09-07 Klassifikationsvorverarbeitung in der medizinischen ultraschallscherwellenbildgebung Active DE102012108353B4 (de)

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
US13/229,408 US10338203B2 (en) 2011-09-09 2011-09-09 Classification preprocessing in medical ultrasound shear wave imaging
US13/229,408 2011-09-09

Publications (2)

Publication Number Publication Date
DE102012108353A1 true DE102012108353A1 (de) 2013-03-14
DE102012108353B4 DE102012108353B4 (de) 2023-01-19

Family

ID=47740324

Family Applications (1)

Application Number Title Priority Date Filing Date
DE102012108353.3A Active DE102012108353B4 (de) 2011-09-09 2012-09-07 Klassifikationsvorverarbeitung in der medizinischen ultraschallscherwellenbildgebung

Country Status (5)

Country Link
US (1) US10338203B2 (de)
JP (1) JP6140953B2 (de)
KR (1) KR102051293B1 (de)
CN (1) CN102988082B (de)
DE (1) DE102012108353B4 (de)

Families Citing this family (18)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20120259224A1 (en) * 2011-04-08 2012-10-11 Mon-Ju Wu Ultrasound Machine for Improved Longitudinal Tissue Analysis
US9468421B2 (en) * 2012-02-16 2016-10-18 Siemens Medical Solutions Usa, Inc. Visualization of associated information in ultrasound shear wave imaging
US9332962B2 (en) * 2013-03-13 2016-05-10 Siemens Medical Solutions Usa, Inc. Ultrasound ARFI displacement imaging using an adaptive time instance
CN104640506B (zh) * 2013-09-18 2017-06-30 东芝医疗系统株式会社 超声波诊断装置、医用图像处理装置和医用图像处理方法
US10143442B2 (en) * 2013-10-24 2018-12-04 Ge Medical Systems Global Technology, Llc Ultrasonic diagnosis apparatus
JP5730979B2 (ja) * 2013-11-08 2015-06-10 日立アロカメディカル株式会社 超音波診断装置、及び弾性評価方法
US10451587B2 (en) * 2014-07-16 2019-10-22 Duke University Methods, systems and computer program products for estimating shear wave speed using statistical inference
CN106716172B (zh) * 2014-08-14 2020-06-05 皇家飞利浦有限公司 用于流体池检测和识别的声流
ES2939246T3 (es) * 2014-10-28 2023-04-20 Supersonic Imagine Métodos y aparatos de formación de imágenes para realizar la formación de imágenes por elastografía por ondas de corte
CN104605890B (zh) * 2014-12-18 2017-04-05 深圳开立生物医疗科技股份有限公司 一种剪切波峰值波形校正的方法、装置、系统及其应用
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
US10537305B2 (en) 2015-09-30 2020-01-21 General Electric Company Detecting amniotic fluid position based on shear wave propagation
US11006928B2 (en) * 2016-02-10 2021-05-18 Siemens Medical Solutions Usa, Inc. Sound speed imaging using shear waves
US11138723B2 (en) * 2016-07-22 2021-10-05 Canon Medical Systems Corporation Analyzing apparatus and analyzing method
CN106725610B (zh) * 2016-11-29 2019-08-06 深圳大学 基于移动声束相干激励剪切波的弹性测量方法及系统
US11443422B2 (en) * 2017-10-13 2022-09-13 The Cleveland Clinic Foundation Advanced ultrasonic detection of different tissue types
US11006926B2 (en) 2018-02-27 2021-05-18 Siemens Medical Solutions Usa, Inc. Region of interest placement for quantitative ultrasound imaging
CN109907731B (zh) * 2019-01-31 2021-06-01 浙江大学 基于特征空间的光学相干层析的三维血流造影方法

Family Cites Families (25)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5487387A (en) * 1994-06-03 1996-01-30 Duke University Method and apparatus for distinguishing between solid masses and fluid-filled cysts
US5840038A (en) * 1997-05-29 1998-11-24 Marquette Medical Systems, Inc. Method and apparatus for signal averaging and analyzing high resolution P wave signals from an electrocardiogram
FR2778838A1 (fr) * 1998-05-19 1999-11-26 Koninkl Philips Electronics Nv Procede de detection de variations d'elasticite et appareil echographique pour mettre en oeuvre ce procede
US6594595B2 (en) * 2001-04-03 2003-07-15 Advantest Corporation Apparatus for and method of measuring cross-correlation coefficient between signals
FR2844058B1 (fr) * 2002-09-02 2004-11-12 Centre Nat Rech Scient Procede et dispositif d'imagerie utilisant des ondes de cisaillement
JPWO2006038596A1 (ja) 2004-10-05 2008-05-15 アステラス製薬株式会社 涙液及び唾液乾燥症治療用医薬組成物
US8858441B2 (en) * 2005-05-12 2014-10-14 The Trustees Of Columbia University In The City Of New York System and method for electromechanical wave imaging of body structures
FR2899336B1 (fr) * 2006-03-29 2008-07-04 Super Sonic Imagine Procede et dispositif pour l'imagerie d'un milieu viscoelastique
US7894874B2 (en) * 2006-05-08 2011-02-22 Luna Innovations Incorporated Method and apparatus for enhancing the detecting and tracking of moving objects using ultrasound
CA2652742C (en) * 2006-05-26 2016-09-06 Queen's University At Kingston Method for improved ultrasonic detection
US8118744B2 (en) * 2007-02-09 2012-02-21 Duke University Methods, systems and computer program products for ultrasound shear wave velocity estimation and shear modulus reconstruction
US20080242997A1 (en) * 2007-03-26 2008-10-02 Luna Innovations Incorporated Method and apparatus for classifying gaseous and non-gaseous objects
US9839395B2 (en) * 2007-12-17 2017-12-12 Dexcom, Inc. Systems and methods for processing sensor data
US8583204B2 (en) * 2008-03-28 2013-11-12 Dexcom, Inc. Polymer membranes for continuous analyte sensors
US8187187B2 (en) 2008-07-16 2012-05-29 Siemens Medical Solutions Usa, Inc. Shear wave imaging
US9554770B2 (en) 2008-09-29 2017-01-31 Siemens Medical Solutions Usa, Inc. High pulse repetition frequency for detection of tissue mechanical property with ultrasound
JP5679988B2 (ja) * 2008-12-22 2015-03-04 コーニンクレッカ フィリップス エヌ ヴェ 印加された力に対する組織変位のリアルタイムの監視のためのアブレーション制御装置
US20100241001A1 (en) 2009-03-20 2010-09-23 Palmeri Mark L Ultrasound Methods, Systems and Computer Program Products for Imaging Fluids
US20100286520A1 (en) * 2009-05-11 2010-11-11 General Electric Company Ultrasound system and method to determine mechanical properties of a target region
AU2010295484B2 (en) * 2009-09-17 2014-03-20 University Of Virginia Patent Foundation Ultrasound-based method and related system to evaluate hemostatic function of whole blood
JP6148010B2 (ja) 2009-11-25 2017-06-14 コーニンクレッカ フィリップス エヌ ヴェKoninklijke Philips N.V. 焦束されたスキャンラインビーム形成での超音波剪断波撮像
US8715185B2 (en) * 2010-04-05 2014-05-06 Hitachi Aloka Medical, Ltd. Methods and apparatus for ultrasound imaging
US10004474B2 (en) * 2010-10-27 2018-06-26 Siemens Medical Solutions Usa, Inc. Tissue density quantification using shear wave information in medical ultrasound scanning
US8494791B2 (en) * 2010-12-16 2013-07-23 General Electric Company Methods and systems for improved correlation of shear displacement waveforms
US20130102932A1 (en) * 2011-10-10 2013-04-25 Charles A. Cain Imaging Feedback of Histotripsy Treatments with Ultrasound Transient Elastography

Also Published As

Publication number Publication date
KR20130028867A (ko) 2013-03-20
US20130066204A1 (en) 2013-03-14
JP2013059622A (ja) 2013-04-04
DE102012108353B4 (de) 2023-01-19
CN102988082A (zh) 2013-03-27
CN102988082B (zh) 2015-11-18
KR102051293B1 (ko) 2019-12-03
US10338203B2 (en) 2019-07-02
JP6140953B2 (ja) 2017-06-07

Similar Documents

Publication Publication Date Title
DE102012108353B4 (de) Klassifikationsvorverarbeitung in der medizinischen ultraschallscherwellenbildgebung
DE102013010060B4 (de) Messen der akustischen Absorption oder Dämpfung von Ultraschall
DE102013001230B4 (de) Achsbezogene Charakterisierung von Scherwellen mit Ultraschall
DE102011114943B4 (de) Bestimmen von Scherwelleninformationen in der medizinischen Ultraschallbildgebung
DE102012101312B4 (de) Viskoelastizitätsmessung unter Verwendung einer amplitudenphasenmodulierten Ultraschallwelle
DE102016100367B4 (de) Spärliche Verfolgung in Schallstrahlintensitätsimpuls-Bildgebung
DE102014002747A1 (de) Arfi-ultraschall-verlagerungsbildgebung mittels eines adaptiven zeitpunkts
DE102012101395B4 (de) Mehrstrahliger Spektral-Doppler in der medizinischen diagnostischen Ultraschallbildgebung
DE102009033286B4 (de) Scherwellenbildgebung
DE102017211895A1 (de) Gewebecharakterisierung im medizinischen diagnostischen Ultraschall
DE102013011510A1 (de) Automatische Doppler-Gate-Positionierung in der Spektraldopplersonografie
DE3686401T2 (de) Darstellung von stromlinien in inhomogenen medien.
DE102012022010A1 (de) Adaptive Bildoptimierung in der welleninduzierten Ultraschallbildgebung
DE102019202545A1 (de) Platzierung der Region von Interesse für die quantitative Ultraschallbildgebung
DE102017202141B4 (de) Schallgeschwindigkeitsbildgebung unter Verwendung von Scherwellen
DE202004021722U1 (de) Vorrichtung zum Gewinnen eines Volumenscans eines sich periodisch bewegenden Objekts
DE102016106998A1 (de) Quantitative viskoelastische Ultraschallbildgebung
DE102014003105A1 (de) Fettanteilschatzung mittels ultraschall mit scherwellenausbreitung
DE102018218319A1 (de) Viskoelastische Gewebeschätzung von Schergeschwindigkeit in medizinischer Ultraschallbildgebung
DE102016114783A1 (de) Adaptive Bewegungsschätzung bei Schallstrahlungskraftbildgebung
DE102009044028A1 (de) Systeme und Verfahren zum Erfassen von Gebieten veränderter Steifigkeit
DE102013002065B4 (de) Visualisierung von zugehörigen lnformationen bei der Ultraschall-Scherwellenbildgebung
DE112019001425T5 (de) Ultraschallsystem zur dreidimensionalen darstellung von scherwellen
DE10308320A1 (de) Verfahren und Vorrichtung zur Unterdrückungsfilterung hoher Verformungsraten
DE112020005100T5 (de) Zusätzliche Diagnosedaten in der parametrischen Ultraschallbildgebung

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