DE102016116658A1 - Sparkleartefakt-erkennung im ultraschall-farbfluss - Google Patents

Sparkleartefakt-erkennung im ultraschall-farbfluss Download PDF

Info

Publication number
DE102016116658A1
DE102016116658A1 DE102016116658.8A DE102016116658A DE102016116658A1 DE 102016116658 A1 DE102016116658 A1 DE 102016116658A1 DE 102016116658 A DE102016116658 A DE 102016116658A DE 102016116658 A1 DE102016116658 A1 DE 102016116658A1
Authority
DE
Germany
Prior art keywords
locations
sparkle
values
doppler
color flow
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.)
Pending
Application number
DE102016116658.8A
Other languages
English (en)
Inventor
Chi Hyung Seo
King Yuen Wong
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 DE102016116658A1 publication Critical patent/DE102016116658A1/de
Pending legal-status Critical Current

Links

Images

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F18/00Pattern recognition
    • G06F18/20Analysing
    • G06F18/22Matching criteria, e.g. proximity measures
    • 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/5269Devices using data or image processing specially adapted for diagnosis using ultrasonic, sonic or infrasonic waves involving detection or reduction of artifacts
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B8/00Diagnosis using ultrasonic, sonic or infrasonic waves
    • A61B8/06Measuring blood flow
    • 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
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
    • G01S15/00Systems using the reflection or reradiation of acoustic waves, e.g. sonar systems
    • G01S15/88Sonar systems specially adapted for specific applications
    • G01S15/89Sonar systems specially adapted for specific applications for mapping or imaging
    • G01S15/8906Short-range imaging systems; Acoustic microscope systems using pulse-echo techniques
    • G01S15/8979Combined Doppler and pulse-echo imaging systems
    • G01S15/8988Colour Doppler imaging
    • 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
    • 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/52025Details of receivers for pulse systems
    • G01S7/52026Extracting wanted echo signals
    • 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/5205Means for monitoring or calibrating
    • 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/52053Display arrangements
    • G01S7/52057Cathode ray tube displays
    • G01S7/52071Multicolour displays; using colour coding; Optimising colour or information content in displays, e.g. parametric imaging
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T5/00Image enhancement or restoration
    • G06T5/70Denoising; Smoothing
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00Image analysis
    • G06T7/20Analysis of motion
    • G06T7/254Analysis of motion involving subtraction of images
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06VIMAGE OR VIDEO RECOGNITION OR UNDERSTANDING
    • G06V10/00Arrangements for image or video recognition or understanding
    • G06V10/20Image preprocessing
    • G06V10/30Noise filtering
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B8/00Diagnosis using ultrasonic, sonic or infrasonic waves
    • A61B8/46Ultrasonic, sonic or infrasonic diagnostic devices with special arrangements for interfacing with the operator or the patient
    • A61B8/461Displaying means of special interest
    • A61B8/463Displaying means of special interest characterised by displaying multiple images or images and diagnostic data on one display
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B8/00Diagnosis using ultrasonic, sonic or infrasonic waves
    • A61B8/48Diagnostic techniques
    • A61B8/488Diagnostic techniques involving Doppler signals
    • 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
    • G06T2207/00Indexing scheme for image analysis or image enhancement
    • G06T2207/10Image acquisition modality
    • G06T2207/10132Ultrasound image
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06VIMAGE OR VIDEO RECOGNITION OR UNDERSTANDING
    • G06V2201/00Indexing scheme relating to image or video recognition or understanding
    • G06V2201/03Recognition of patterns in medical or anatomical images
    • G06V2201/031Recognition of patterns in medical or anatomical images of internal organs

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Health & Medical Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Remote Sensing (AREA)
  • Radar, Positioning & Navigation (AREA)
  • General Physics & Mathematics (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • General Health & Medical Sciences (AREA)
  • Molecular Biology (AREA)
  • Veterinary Medicine (AREA)
  • Public Health (AREA)
  • Animal Behavior & Ethology (AREA)
  • Surgery (AREA)
  • Biophysics (AREA)
  • Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
  • Pathology (AREA)
  • Radiology & Medical Imaging (AREA)
  • Biomedical Technology (AREA)
  • Heart & Thoracic Surgery (AREA)
  • Medical Informatics (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Theoretical Computer Science (AREA)
  • Acoustics & Sound (AREA)
  • Multimedia (AREA)
  • Hematology (AREA)
  • Data Mining & Analysis (AREA)
  • Bioinformatics & Computational Biology (AREA)
  • Artificial Intelligence (AREA)
  • Bioinformatics & Cheminformatics (AREA)
  • General Engineering & Computer Science (AREA)
  • Evolutionary Biology (AREA)
  • Evolutionary Computation (AREA)
  • Ultra Sonic Daignosis Equipment (AREA)
  • Physiology (AREA)

Abstract

Es wird ein Sparkle bei der Farbfluss-Bildgebung erkannt. Farbflussdaten werden mit unterschiedlicher Impulswiederholungsfrequenz (PRF) geschätzt. Durch Korrelieren der mit unterschiedlichen PRFs geschätzten Farbflussdaten wird ein Sparkle identifiziert. Farbflussbilder können gefiltert werden, um eine Bewegung zu reduzieren, während die Sparkle-Region beibehalten wird (z. B. Nierenstein-Abbildung), oder um die Sparkle-Region zu reduzieren, während die Bewegung beibehalten wird (z. B. Entfernen des Sparkle als Systemrauschen).

Description

  • HINTERGRUND
  • Die vorliegenden Ausführungsformen betreffen die Sparkle-Erkennung bei der Farbfluss-Bildgebung. Die Farbfluss-Bildgebung ist für Artefakte anfällig, wenn ein Korrelationsverlust in den zurückgestreuten Echos auftritt. Ein Korrelationsverlust kann zugrunde liegende Physiologien, wie z. B. Nierensteine, aufgrund einer Oberflächenrauheit anzeigen. Ein Korrelationsverlust kann stattdessen auch an Systemunvollkommenheiten liegen, wie z. B. Phasenrauschen oder Echorauschen. Normalerweise wird die Verstärkung oder die Sendeleistung zurückgefahren, um dieses Sparkle-Artefakt zu reduzieren, allerdings auf Kosten der Empfindlichkeit. Unter Verwendung der räumlichen Varianz des Farbflusses kann das Sparkle identifiziert werden, jedoch kann ein Alias-Fluss oder eine Turbulenz ebenso wie Artefakte verringert werden.
  • KURZDARSTELLUNG
  • Einführend umfassen die bevorzugten Ausführungsformen ein Verfahren, ein System, ein computerlesbares Medium und Befehle zur Sparkle-Erkennung bei der Farbfluss-Bildgebung. Farbflussdaten werden mit unterschiedlicher Impulswiederholungsfrequenz (Pulse Repetition Frequency, PRF) geschätzt. Durch korrelieren der mit unterschiedlichen PRFs geschätzten Farbflussdaten wird ein Sparkle identifiziert. Farbflussbilder können gefiltert werden, um eine Bewegung zu verringern, während die Sparkle-Region beibehalten wird (z. B. Nierenstein-Abbildung), oder um die Sparkle-Region zu verringern, während die Bewegung beibehalten wird (z. B. Entfernen des Sparkle als Systemrauschen).
  • In einer ersten Erscheinungsform wird ein Verfahren zum Erkennen von Sparkle-Artefakten beim Farbfluss bereitgestellt. Es werden erste Farbflussdaten erzeugt, welche Stellen in einem Patienten repräsentieren. Die ersten Farbflussdaten werden mit einer ersten Impulswiederholungsfrequenz erzeugt. Es werden zweite Farbflussdaten erzeugt, welche die Stellen in dem Patienten repräsentieren. Die zweiten Farbflussdaten werden mit einer zweiten Impulswiederholungsfrequenz erzeugt. Für jede der Stellen wird ein Maß an Ähnlichkeit zwischen den ersten und zweiten Farbflussdaten bestimmt. Ein Farbflussbild wird gefiltert, wobei die Filterung auf dem Maß an Ähnlichkeit als Funktion der Stelle basiert. Das gefilterte Farbflussbild wird angezeigt.
  • In einer zweiten Erscheinungsform sind auf einem permanenten computerlesbaren Speichermedium Daten gespeichert, welche Befehle repräsentieren, die durch einen programmierten Prozessor zur Sparkle-Erkennung bei der Doppler-Bildgebung ausführbar sind. Das Speichermedium umfasst Befehle zum: Abtasten mehrerer Stellen mit Ultraschall, wobei das Abtasten jeder der Stellen eine Fluss-Abtastprobenanzahl von zurückkehrenden Abtastproben liefert; Schätzen erster Doppler-Werte für die Stellen unter Verwendung der zurückkehrenden Abtastproben der Fluss-Abtastprobenanzahl; Schätzen zweiter Doppler-Werte für die Stellen unter Verwendung einer Teilabtastung der zurückkehrenden Abtastproben der Fluss-Abtastprobenanzahl; Multiplizieren der ersten Doppler-Werte mit einem Faktor, der eine Funktion der Teilabtastung ist; Korrelieren der zweiten Doppler-Werte mit den multiplizierten ersten Doppler-Werten und Erkennen des Sparkle aus Ergebnissen der Korrelation.
  • In einer dritten Erscheinungsform wird ein System zur Sparkle-basierten Verarbeitung in Fließbildern bereitgestellt. Ein Messaufnehmer und ein Strahlformer werden zum Abtasten einer Abtastregion bereitgestellt. Eine Doppler-Schätzeinheit ist dafür konfiguriert, aus der Abtastung erste Bewegungswerte zu schätzen, welche Stellen der Abtastregion repräsentieren, und zweite Bewegungswerte zu schätzen, welche die Stellen der Abtastregion repräsentieren. Die ersten Bewegungswerte werden mit einer anderen Anzahl von Abtastproben aus der Abtastung geschätzt als die zweiten Bewegungswerte. Ein Prozessor ist dafür konfiguriert, anhand eines Vergleichs der ersten und zweiten Bewegungswerte zu identifizieren, welche der Stellen ein Sparkle aufweisen.
  • Die vorliegende Erfindung wird durch die folgenden Patentansprüche definiert und nichts in diesem Abschnitt sollte als Beschränkung jener Patentansprüche angesehen werden. Weitere Erscheinungsformen und Vorteile der Erfindung werden nachstehend in Verbindung mit den bevorzugten Ausführungsformen erörtert.
  • KURZBESCHREIBUNG DER ZEICHNUNGEN
  • Die Komponenten und die Figuren sind nicht notwendigerweise maßstabsgetreu, stattdessen liegt die Betonung darauf, die Prinzipien der Erfindung zu veranschaulichen. Überdies sind in den Figuren überall in den verschiedenen Ansichten entsprechende Teile durch gleiche Bezugszahlen gekennzeichnet.
  • 1 ist ein Ablaufplan einer Ausführungsform eines Verfahrens zur Sparkle-Erkennung bei der Bewegungs-Bildgebung;
  • 2A bis 2D sind beispielhafte Schaubilder von Gewichtungen für die Filterung als eine Funktion der Ähnlichkeit;
  • 3A und 3B sind beispielhafte Geschwindigkeitsbilder, erzeugt mit unterschiedlicher PRF, 3C ist ein binäres Filter oder eine Maske zum Entfernen von Sparkle und 3D ist das Fließbild der 3A nach dem Anwenden des Filters der 3C;
  • 4A und 4B sind beispielhafte Geschwindigkeitsbilder, erzeugt mit unterschiedlicher PRF, 4C ist ein binäres Filter oder eine Maske und 4D ist das Fließbild der 4A nach dem Filtern zum Verstärken von Sparkle-Regionen und
  • 5 ist ein Blockschaubild einer Ausführungsform eines Systems zur Sparkle-basierten Verarbeitung in Fließbildern.
  • DETAILLIERTE BESCHREIBUNG DER ZEICHNUNGEN UND DERZEIT BEVORZUGTE AUSFÜHRUNGSFORMEN
  • Es werden Sparkle-Artefakte bei der Ultraschall-Bildgebung erkannt. Sparkle, manchmal auch als Funkeln bezeichnet, ist zufälliger als eine Fluid- oder Gewebebewegung. Durch Variieren der Impulswiederholungsfrequenz (PRF) können unterschiedliche Informationen für Sparkle-Artefakte erzeugt werden. Eine aus dem Blutfluss geschätzte Geschwindigkeit ist im Allgemeinen unabhängig von der PRF. Durch Erzeugen von zwei Bildern mit unterschiedlichen PRFs und Vornehmen einer normalisierten Kreuzkorrelation kann eine Gewichtungsmatrix angewendet werden, um Bilder mit reduziertem Sparkle (z. B. nur Fluss oder saubereres Fließbild) oder mit reduziertem Fluss (z. B. nur Artefakt oder verstärkter Stein) zu erzeugen. Bei Verwendung der PRF-Variation können Sparkle-Artefakte erkannt werden, ohne die Empfindlichkeit zu beeinträchtigen.
  • 1 zeigt eine Ausführungsform eines Verfahrens zur Sparkle-Artefakt-Erkennung beim Farbfluss. Ein Farbfluss wird verwendet, um eine Raumbewegungs-Bildgebung anzuzeigen, wie z. B. eine Fluid- oder Gewebebewegung. „Farbe” wird verwendet, um von der Spektral-Doppler-Bildgebung zu unterscheiden, wobei das Leistungsspektrum für ein Bereichstor geschätzt wird. Die Farb-„Fluss”-Daten müssen nicht von einem Fluid sein (z. B. können sie von einer Gewebebewegung sein) und/oder müssen keine Farbe repräsentieren (z. B. können sie ein Skalar sein). Das Sparkle-Artefakt wird erkannt, um ein zurückkehrendes Signal von Steinen zu verstärken oder um ein Systemrauschen für die Farbfluss-Bildgebung zu verringern.
  • Das Verfahren wird durch das Ultraschall-Bildgebungssystem 10 der 5, den Prozessor 24 oder ein anderes System, ein anderes Filter und/oder einen anderen Prozessor durchgeführt. Zum Beispiel führt das Ultraschall-Bildgebungssystem 10 die Schritte aus. Als ein anderes Beispiel steuert der Prozessor 24 im Schritt 30 einen Strahlformer zum Abtasten, bewirkt in den Schritten 32 bis 36 die Erzeugung von Bewegungsdaten durch eine Doppler-Schätzeinheit und bewirkt im Schritt 46 die Erzeugung des Bildes durch einen Abtastungswandler, einen Grafikspeicher und/oder eine Anzeigevorrichtung, bestimmt aber selbst die Ähnlichkeit im Schritt 40, erkennt im Schritt 42 ein Sparkle und filtert im Schritt 44. In noch einem anderen Beispiel führt ein separates Filter den Schritt 44 aus.
  • Die Schritte der 1 werden in der dargestellten Reihenfolge oder einer anderen Reihenfolge ausgeführt. Zum Beispiel werden die Schritte 34 und 36 in der dargestellten Reihenfolge oder einer entgegengesetzten Reihenfolge ausgeführt. Als ein anderes Beispiel werden die Schritte 30 und 32 für eine PRF wiederholt, nachdem die Schritte 30 und 32 für eine andere PRF ausgeführt worden sind.
  • Es können zusätzliche, andere oder weniger Schritte angewendet werden, als in 1 dargestellt. Zum Beispiel wird der Schritt 42 nicht separat ausgeführt, sondern stattdessen wird das Filtern des Schritts 44 durchgeführt, wobei die Ähnlichkeit die Sparkle-Stellen anzeigt. In einem anderen Beispiel wird im Schritt 46 kein Bild erzeugt. Die Schritte 34 und 36 sind eine Ausführungsform zum Ausführen des Schritts 32, es können jedoch auch andere Ausführungsformen bereitgestellt werden.
  • Im Schritt 30 werden verschiedene Stellen innerhalb einer Abtastregion des Patienten mit Ultraschall abgetastet. In einer Ausführungsform wird unter Verwendung eines Ultraschallsystems ein Patient oder eine Region in Echtzeit mit der Bildgebung abgetastet. Bei der abgetasteten Region handelt es sich um ein Inneres eines Objekts, z. B. des Patienten. Die Abtastung ist von einer Volumen-, Ebenen- oder Linienregion. Das Abtasten einer Ebene liefert Daten, die verschiedene Stellen oder Abtastproben der Ebene repräsentieren. Die Daten, welche die Region repräsentieren, werden aus einer räumlichen Abtastung des Objekts gebildet. Die räumlichen Abtastproben sind für Stellen, die in einem akustischen Abtastgitter verteilt sind.
  • Die Region für die Farbflussabtastung ist eine Region von Interesse, die kleiner ist als ein Gesichtsfeld, oder für das gesamte Gesichtsfeld. Das Ultraschallsystem kann das Gesichtsfeld durch B-Modus-Bildgebung abtasten. Die Farbflussregion ist eine Teilmenge dieses Gesichtsfelds. Der Benutzer oder ein Prozessor bestimmt die Region von Interesse, in welcher die Farbflussabtastung erfolgt. Alternativ ist die Farbflussregion das vollständige Gesichtsfeld.
  • Es werden räumliche Abtastproben entlang einer oder mehreren Abtastlinien empfangen. Wo der Sendestrahl nur eine Empfangs-Abtastlinie beschallt, werden zurückkehrende oder Echo-Abtastproben entlang dieser Abtastlinie empfangen. Wo der Sendestrahl mehrere Abtastlinien beschallt, können Abtastproben entlang den mehreren Abtastlinien empfangen werden. Um die Abtastproben für verschiedene Empfangsstrahlen gleichzeitig zu erzeugen, wird eine parallele Empfangsstrahlformung durchgeführt. Zum Beispiel kann ein System in der Lage sein, zwei oder mehr, mehrere zehn oder hunderte von Empfangsstrahlen parallel zu formen. Alternativ werden von den Elementen empfangene Signale gespeichert und nacheinander verarbeitet. Räumliche Abtastproben werden für mehrere Empfangslinien in der Region von Interesse in Reaktion auf einen und/oder in Reaktion auf aufeinander folgende Sendestrahlen erhalten.
  • Das Abtasten kann mehrfach durchgeführt werden, um die Region abzudecken. Die Schritte werden wiederholt, um verschiedene Abschnitte der Region von Interesse abzutasten. Alternativ liefert eine einfache Ausführung die Daten für die gesamte Region von Interesse.
  • Die vollständige Region von Interesse wird mehrere Male nacheinander abgetastet. Ein Abtasten zu verschiedenen Zeiten nacheinander liefert räumliche Abtastproben, die mit einer Bewegung verbunden sind. Es können beliebige derzeit bekannte oder später entwickelte Impulssequenzen verwendet werden. Es wird eine Sequenz von mindestens zwei (Fluss-Abtastprobenanzahl) Aussendungen entlang jeder Abtastlinie bereitgestellt. Zum Beispiel beträgt die Fluss-Abtastprobenanzahl 10 bis 20, was zu 10 bis 20 Abtastproben für jede Stelle führt. Es können beliebige Impulswiederholungsfrequenzen (d. h. Abtastgeschwindigkeiten für eine Stelle), Fluss-Abtastprobenanzahlen (d. h. Anzahlen an Abtastproben für eine Stelle oder Anzahlen an Abtastproben, die zum Schätzen verwendet werden) und Impulswiederholungsintervalle (d. h. Zeiten zwischen den jeweiligen Abtastprobenerwerben für eine Stelle) verwendet werden.
  • Die Echo-Antworten auf die Aussendungen der Sequenz oder die zurückkehrenden Abtastproben werden verwendet, um Geschwindigkeit, Energie (Leistung) und/oder Varianz zu einer gegebenen Zeit zu schätzen. Die Aussendungen entlang einer Linie (Linien) können mit Aussendungen entlang einer anderen Linie (anderen Linien) verschachtelt werden. Mit oder ohne Verschachtelung werden die räumlichen Abtastproben für eine gegebene Zeit unter Verwendung von Aussendungen von verschiedenen Zeiten erhalten. Die Schätzungen von verschiedenen Abtastlinien können nacheinander erhalten werden, aber schnell genug, um aus einer Benutzerperspektive eine gleiche Zeit zu repräsentieren. Es können mehrere Abtastungen durchgeführt werden, um Schätzungen für verschiedene Zeiten zu erhalten.
  • Um Daten in Reaktion auf verschiedene PRFs oder Fluss-Abtastprobenanzahlen zu erzeugen, können die Abtastproben mit einer Abtastung bei einer PRF oder Fluss-Abtastprobenanzahl erhalten werden. Eine Teilgruppe dieser Abtastproben kann verwendet werden, um bei der niedrigeren PRF zu schätzen. Alternativ wird die Abtastung wiederholt. Für eine erste Wiederholung werden die zurückkehrenden Abtastproben bei einer PRF erhalten. Für eine anschließende Wiederholung werden die zurückkehrenden Abtastproben bei der anderen PRF erhalten. Durch ein nacheinander erfolgendes Aussenden bei den verschiedenen PRFs, um zurückkehrende Abtastproben in verschiedenen Fluss-Abtastprobenanzahlen zu erhalten, werden Abtastproben zum Schätzen eines Flusses mit verschiedenen PRFs erhalten.
  • In alternativen Ausführungsformen werden die zurückkehrenden Abtastproben oder Fluss-Schätzungen durch Übertragung über ein Netzwerk und/oder Laden aus einem Speicher erhalten. Es werden Daten erhalten, die zuvor durch Abtasten erhalten wurden.
  • Im Schritt 32 erzeugt eine Schätzeinheit oder ein Detektor Farbflussdaten, welche Stellen in dem Patienten repräsentieren. Farbflussdaten umfassen Schätzungen der Geschwindigkeit, Energie (z. B. Leistung) und/oder Varianz. In einer Ausführungsform wird zumindest die Geschwindigkeit geschätzt. Die Farbflussdaten können für Fluid oder Gewebe sein. Es können Schätzungen der Geschwindigkeit, Energie und/oder Varianz von Gewebebewegung erzeugt werden. Es können beliebige Bewegungsdaten erhalten werden, entweder aus Fluid- oder Gewebebewegung. Farbflussdaten werden in den folgenden Beispielen als Fluidbewegung verwendet, können aber alternativ oder zusätzlich auch Gewebebewegungsdaten sein.
  • Die empfangenen räumlichen Abtastproben können Clutter-gefiltert werden. Das Clutter-Filter lässt Frequenzen durch, die zu einer Fluid- und nicht zu einer Gewebebewegung gehören oder zu einer Gewebebewegung und nicht zu einer Fluidbewegung gehören. Das Clutter-Filtern erfolgt an Signalen in der Impulssequenz zum Schätzen einer Bewegung zu einer gegebenen Zeit (z. B. Abtastproben der Fluss-Abtastprobenanzahl). Ein gegebenes Signal kann für Schätzungen verwendet werden, die verschiedene Zeiten repräsentieren, z. B. verbunden mit einem sich bewegenden Fenster für die Clutter-Filterung und Schätzung. Es werden verschiedene Filterausgaben verwendet, um eine Bewegung für eine Stelle zu verschiedenen Zeiten zu schätzen.
  • Aus den räumlichen Abtastproben werden Farbflussdaten erzeugt. Es kann eine Doppler-Verarbeitung, wie z. B. Autokorrelation, angewendet werden. In anderen Ausführungsformen kann eine zeitliche Korrelation angewendet werden. Es kann ein anderes Verfahren angewendet werden, um die Farbflussdaten zu schätzen. Farbdoppler-Parameterwerte (z. B. Geschwindigkeits-, Energie- oder Varianzwerte) werden aus den räumlichen Abtastproben geschätzt, die zu verschiedenen Zeiten erhalten werden. Die Änderung der Frequenz (z. B. Dopplerverschiebung) zwischen zwei Abtastproben für die gleiche Stelle zu verschiedenen Zeiten zeigt die Geschwindigkeit an. Es kann eine Sequenz (Fluss-Abtastprobenanzahl) von zwei oder mehr Abtastproben verwendet werden, um die Farbdoppler-Parameterwerte zu schätzen. Schätzungen werden für verschiedene Gruppierungen empfangener Signale gebildet, z. B. vollständig separate oder unabhängige Gruppierungen oder überlappende Gruppierungen. Die Schätzungen für jede Gruppierung repräsentieren die räumliche Stelle zu einer gegebenen Zeit.
  • Die Schätzung wird für die verschiedenen abgetasteten räumlichen Stellen durchgeführt. Zum Beispiel werden aus Echos in Reaktion auf die Abtastung Geschwindigkeiten für die verschiedenen Stellen in einer Ebene geschätzt. Es können mehrere Rahmen von Farbflussdaten erhalten werden, um entsprechend die Region von Interesse zu verschiedenen Zeiten zu repräsentieren.
  • Die Schätzungen können Schwellenwerte aufweisen. Schwellenwerte werden für die Geschwindigkeiten und/oder Leistungen angewendet. Zum Beispiel wird ein unterer Geschwindigkeits-Schwellenwert angewendet. Geschwindigkeiten unterhalb des Schwellenwerts werden entfernt oder auf einen anderen Wert gesetzt, z. B. Null. Als ein anderes Beispiel wird, wo die Energie unterhalb eines Schwellenwerts liegt, der Geschwindigkeitswert für dieselbe räumliche Stelle entfernt oder auf einen anderen Wert gesetzt, z. B. Null. Alternativ werden die geschätzten Geschwindigkeiten ohne Schwellenwerte verwendet.
  • Bei den erhaltenen Bewegungs- oder Farbflussdaten handelt es sich um einen Datenrahmen oder ein Bild, welches den Patienten zu einer gegebenen Zeit repräsentiert, obwohl sie aus empfangenen Signalen über die Fluss-Abtastprobenanzahl geschätzt sind. Es können auch andere Daten erzeugt werden, z. B. B-Modus-Daten. Ein B-Modus-Bild kann mit einer beinhalteten Region von Interesse überlagert sein oder eine solche aufweisen, welche die Farbdoppler-Geschwindigkeiten zeigt. Innerhalb der Region von Interesse sind Stellen ohne Fluss als B-Modus-Daten dargestellt.
  • Um ein Sparkle zu erkennen, werden zwei oder mehr Datenrahmen erzeugt, welche dieselben Stellen repräsentieren, aber bei unterschiedlicher PRF. Durch Erzeugen der Rahmen aus denselben Daten sind die Stellen dieselben. Wo die Rahmen mit sequentieller Abtastung erzeugt werden, sind dieselben Stellen aufgrund einer Messaufnehmer- und/oder Gewebebewegung möglicherweise nicht genau dieselben, sind aber ungefähr dieselben.
  • Die Schritte 34 und 36 zeigen eine Ausführungsform, wobei die zwei Rahmen mit unterschiedlicher PRF unter Verwendung eines Faktors 2 erzeugt werden. Im Schritt 34 wird ein Datenrahmen für eine PRF erzeugt. Zum Beispiel werden Doppler-Werte unter Verwendung der zurückkehrenden Abtastproben einer Fluss-Abtastprobenanzahl (z. B. 16) geschätzt. Geschwindigkeiten oder andere Bewegungswerte werden für jede der Stellen unter Verwendung der PRF und einer entsprechenden Fluss-Abtastprobenanzahl geschätzt.
  • Die wahre Geschwindigkeit vt eines Flusspixels kann wie folgt ausgedrückt werden:
    Figure DE102016116658A1_0002
    wobei z die Farb-Ensemble-Abtastproben sind, vd die Dopplergeschwindigkeits-Schätzung ist, i ein Index der Abtastproben ist, n eine ganze Zahl ist und –π ≤ vd < π. Das Geschwindigkeitsbild oder der Rahmen von Geschwindigkeitswerten kann ausgedrückt werden als velimage = vtmod 2π = vd
  • Es können auch andere Ausdrücke des Rahmens von Geschwindigkeitswerten verwendet werden.
  • Wo die Werte aus derselben Fluss-Abtastprobenanzahl zu schätzen sind, ist die PRF für einen der Rahmen für die vollständigen oder alle der Fluss-Abtastproben, was die maximale PRF bei gegebenen erhaltenen Abtastproben liefert. Alternativ wird eine Teilabtastung angewendet (z. B. Erhalten von 20 Abtastproben für die Fluss-Abtastprobenanzahl, aber Verwenden von nur 16).
  • Für den Datenrahmen mit der höheren PRF werden die geschätzten Geschwindigkeiten mit einem Faktor multipliziert. Der Faktor bezieht sich auf das Maß Teilabtastung für die andere PRF oder ist genauso hoch wie dieses oder ist eine Reflexion der Differenz der PRFs. Zum Beispiel werden die resultierenden Geschwindigkeitswerte bei einer Teilabtastung der Fluss-Abtastprobenanzahl mit einem Faktor 2 verdoppelt. Für einen genaueren Vergleich werden die Geschwindigkeitsschätzungen unter Verwendung der vollständigen oder einer größeren Fluss-Abtastprobenanzahl (d. h. doppelte PRF) mit 2 multipliziert. {2vd}mod 2π ist einfach das Originalbild, multipliziert mit 2 und anschließend eingewickelt („wrapped”) mit 2π.
  • Es kann jede Multiplikation angewendet werden, z. B. eine Multiplikation der geschätzten Geschwindigkeiten mit dem Faktor. In einer Ausführungsform werden die Geschwindigkeiten mit 2 multipliziert, wobei eine Arithmetik endlicher Genauigkeit verwendet wird, so dass die Einwicklung während des Overflow inherent ist (d. h., wenn es sich bei den Originalgeschwindigkeiten um 8-Bit-2s-Komplement handelt, wird durch Durchführen der Multiplikation mit 2 in einer 2s-Komplement-Arithmetik automatisch die gewünschte Multiplikation erreicht). Wo der Faktor ein anderer als zwei ist, kann die Multiplikation wie folgt ausgedrückt werden: veltimesδ = {δvd}mod 2π
  • Alternativ werden die geschätzten Geschwindigkeiten nicht mit dem Faktor multipliziert. Die Korrelation zum Erkennen von Sparkle beruht auf Variation nach Stelle, so dass es in einigen Fällen ohne Einwicklung möglicherweise nicht notwendig ist, für jede Stelle ähnliche Stärken zu haben.
  • Zum Vergleich mit den Bewegungsdaten mit einer PRF werden im Schritt 36 Bewegungsdaten mit einer anderen PRF erzeugt. Die Doppler-Schätzeinheit erzeugt einen zweiten Rahmen von Farbflussdaten, welche die Stellen in dem Patienten repräsentieren. Der zweite Farbflussdatenrahmen wird mit einer zweiten Impulswiederholungsfrequenz erzeugt. Es wird eine andere Fluss-Abtastprobenanzahl verwendet. Es wird die gleiche Gruppe zurückkehrender Signale oder eine andere Gruppe zurückkehrender Signale verwendet. Zum Schätzen von Dopplerwerten aus einer anderen Gruppe weist die andere Gruppe eine andere Fluss-Abtastprobenanzahl und/oder PRF auf. Zum Schätzen von Dopplerwerten aus einer gleichen Gruppe wird eine Teilabtastung angewendet. Es werden dieselben zurückkehrenden Abtastproben teilabgetastet, um die PRF zu senken und die Fluss-Abtastprobenanzahl zu verringern.
  • Es kann eine beliebige Verringerung oder Teilabtastung angewendet werden. Zum Beispiel wird die Fluss-Abtastprobenanzahl um die Hälfte oder um einen Faktor 2 verringert. Alles andere (z. B. werden für einen Rahmen die Abtastproben 1 bis 16 verwendet und für den anderen Rahmen werden nur Abtastproben mit gerader oder ungerader Nummer verwendet) oder eine andere Gruppierung führt zu einer Verwendung der Hälfte der zurückkehrenden Abtastproben. Es kann jeder dritte oder eine andere Gruppierung mit Verwendung von weniger als allen, weniger als der Hälfte und oder weniger als den für den anderen Rahmen verwendeten Geschwindigkeitswerten verwendet werden. Es kann ein beliebiger ganzzahliger Teilabtastungsfaktor oder Teilabtastungsfaktor mit reellem Wert verwendet werden.
  • In einer Darstellung wird ein teilabgetastetes Bild velsubsampled erzeugt als:
    Figure DE102016116658A1_0003
  • Für echten Blutfluss ist zu erwarten, dass das teilabgetastete Bild unter Verwendung eines Faktors 2 das Doppelte der wahren Geschwindigkeit beträgt, dargestellt als:
    Figure DE102016116658A1_0004
  • Bei Anwendung dieser Näherung auf das teilabgetastete Bild ergibt sich: velsubsampled ≈ {2vt}mod 2π = {4nπ + 2vd}mod 2π = {2vd}mod 2π, was anzeigt, dass die Geschwindigkeit um den Teilabtastungsfaktor höher ist als die wahre Geschwindigkeit. Die mit der doppelten PRF geschätzten Geschwindigkeiten führen zu Geschwindigkeiten, die doppelt so hoch sind wie Geschwindigkeiten, die mit der PRF geschätzt worden sind.
  • Da eine Teilabtastung mit anderen Faktoren als 2 angewendet werden kann, kann mehr als ein Rahmen der teilabgetasteten Geschwindigkeitsdaten erzeugt werden. Zum Beispiel werden Bewegungsdaten mit Teilabtastungsfaktoren 2 und 3 erzeugt. Die resultierenden Rahmen können gemittelt oder auf andere Weise kombiniert werden. Es werden mehrere Versionen von velsubsampled gemittelt, um die Genauigkeit der wahren Geschwindigkeit zu erhöhen, bevor die Korrelation durchgeführt wird. Zum Beispiel können argΣiz(2i) × z*(2i – 2) und arg Σiz(2i – 1) × z*(2i + 1) gemittelt werden, um velsubsampledAVG zu erzeugen. Der Mittelwert ist allgemeiner auszudrücken als: velsubsampledAVG ≈ AVG({δvd}mod 2π, {δvd'}mod 2π} wobei δ der Teilabtastungsfaktor ist.
  • Im Schritt 40 bestimmt ein Prozessor ein Maß an Ähnlichkeit zwischen den Farbflussdaten, die mit der unterschiedlichen PRF geschätzt werden, und Fluss-Abtastprobenanzahlen. Um die Ähnlichkeit zu bestimmen, werden zwei Datenrahmen verglichen, z. B. Geschwindigkeitsrahmen.
  • Für veltimes2 = {2vd}mod 2π sollten dann für Bereiche, die den Blutfluss repräsentieren, velsubsampled und veltimes2 eine hohe räumliche Korrelation aufweisen. Für Artefakte zufälliger Natur, z. B. Sparkle, ist die räumliche Korrelation niedrig.
  • Die Ähnlichkeit wird als eine Funktion der Stelle bestimmt. Für jede der Stellen wird die Ähnlichkeit zwischen Daten aus verschiedenen Rahmen bestimmt. Die Ähnlichkeit kann eine Differenz sein. Alternativ basiert die Ähnlichkeit auf Kernen beliebiger Größe (z. B. 5 axial und 3 lateral oder die 2D-Kerngröße ist 2A + 1 Abtastproben axial × 2B + 1 Abtastproben lateral), zentriert an der jeweiligen Stelle. Die Ähnlichkeit wird als eine Funktion der Stelle bestimmt.
  • Es kann ein beliebiges Maß an Ähnlichkeit verwendet werden. In einer Ausführungsform wird eine minimale Summe absoluter Differenzen verwendet. In einer anderen Ausführungsform wird eine Korrelation, z. B. eine Kreuzkorrelation, verwendet. Zum Beispiel werden die Dopplerwerte aus der teilabgetasteten Schätzung mit den Dopplerwerten aus der vollständigen Schätzung, multipliziert mit dem Teilabtastungsfaktor, korreliert. Die Kreuzkorrelation oder andere Ähnlichkeit kann normalisiert werden. Das Maß an Ähnlichkeit kann durch normalisierte Kreuzkorrelation (Normalized Cross-Correlation, NCC) zwischen den zwei Signalen velsubsampledAVG und veltimesδ quantifiziert werden. Eine 2D normalisierte Kreuzkorrelation unter Verwendung von Datensegmenten entlang der axialen Richtung und lateralen Richtung bei null Verzögerung wird wie folgt durchgeführt:
    Figure DE102016116658A1_0005
  • Der Koeffizient der normalisierten Kreuzkorrelation ρ(i, j) bei null Verzögerung wird für jede Abtastprobe berechnet. Es können andere Korrelationen oder Ähnlichkeitsmaße zwischen den Datenrahmen, die mit unterschiedlicher PRF geschätzt werden, oder Fluss-Abtastprobenanzahlen verwendet werden.
  • Im Schritt 42 erkennt der Prozessor das Sparkle aus der Ähnlichkeit. Stellen mit einem hohen Korrelationsgrad zwischen Daten mit unterschiedlicher PRF sind kein Sparkle, sondern sind Bewegungen. Bei Stellen mit einem geringen Korrelationsgrad handelt es sich um Sparkle. Die Ergebnisse der Korrelation identifizieren die Stellen mit geringerer Korrelation. Es kann ein beliebiger Schwellenwert verwendet werden, um Sparkle von Bewegung zu unterscheiden, z. B. ein empirisch entwickelter Schwellenwert für eine gegebene Bildgebungsanwendung.
  • In alternativen Ausführungsformen erkennt der Prozessor nicht speziell das Sparkle oder Stellen mit Sparkle. Die Ähnlichkeit als Funktion der Stelle kann die Stellen ohne spezielles Identifizieren einer gegebenen Stelle als Sparkle anzeigen. Die Erkennung erfolgt durch Erzeugung der Ähnlichkeit als eine Funktion der Stelle.
  • Im Schritt 44 filtert der Prozessor oder ein Filter ein Farbflussbild. Bei dem Farbflussbild handelt es sich um Geschwindigkeits- oder andere Bewegungsdaten, die verwendet werden, um ein Bild zu erzeugen (z. B. skalare Schätzung oder Dopplerwerte), oder um RGB- oder andere Daten für die Anzeige. Das Farbflussbild wird aus den Schätzungen erzeugt, die mit der höheren PRF erzeugt werden, z. B. den Geschwindigkeiten, die mit der vollständigen Abtastprobenanzahl erzeugt werden. Alternativ wird das Farbflussbild aus der Teilabtastung erzeugt. In noch einer anderen Alternative wird das Farbflussbild aus Abtastproben von einer anderen Abtastung erzeugt, also nicht aus den Daten erzeugt, die für die Ähnlichkeitsbestimmung verwendet werden.
  • Die Filterung basiert auf dem Ähnlichkeitsgrad als Funktion der Stelle. Um das Sparkle zu unterdrücken, verringert die Filterung stärker Werte an Stellen geringerer Ähnlichkeit. Um das Sparkle zu verstärken, verringert die Filterung stärker Werte an Stellen stärkerer Ähnlichkeit. Durch die Filterung wird vermieden, einige Dopplerwerte zu ändern, es werden aber andere verringert. Alternativ führt die Filterung zu einer Änderung der Höhe aller Werte.
  • Für Geschwindigkeitswerte bedeutet die Verringerung, die Geschwindigkeit näher an Null zu bringen. Die negativen oder positiven Geschwindigkeiten werden verringert, indem negative oder positive Geschwindigkeiten abgesenkt werden, während das Vorzeichen beibehalten wird. Für Energiewerte weisen die Schätzungen kein Vorzeichen auf. Durch die Verringerung werden die Schätzungen näher an Null gebracht.
  • Um zu verringern oder zu filtern, werden die Bewegungsdaten mit Gewichtungen multipliziert. Die Gewichtungen sind Bruchwerte von 0 bis 1. Es kann auch eine andere Gewichtung angewendet werden, z. B. Addieren von, Subtrahieren von oder Dividieren durch Gewichtungen. Das Farbflussbild oder die Gewichtungen werden mit den Gewichtungen als Funktion der Stelle gewichtet.
  • Die Gewichtungen werden aus dem Ähnlichkeitsgrad abgebildet. Es wird eine Pixel-für-Pixel- oder stellenspezifische Gewichtungsmatrix erzeugt. Der Ähnlichkeitsgrad für jede Stelle wird auf eine Gewichtung für diese Stelle abgebildet. Als ein Ergebnis ist die Gewichtung zielabhängig. Die Gewichtungen variieren als Funktion des Ähnlichkeitsniveaus, was für eine Verringerung als Funktion des Sparkle sorgt, z. B. um Bewegungssignale durchzulassen und Sparkle-Artefakte zu entfernen oder zu dämpfen.
  • Es kann ein beliebiges Maß an Verringerung vorgesehen sein. Für eine 2D-Kreuzkorrelation können die Koeffizienten der normalisierten Kreuzkorrelation im Bereich von –1 bis 1 oder von 0 bis 1 liegen. Zwei Signale sind identisch, wenn der Kreuzkorrelationskoeffizient 1 ist, und sie werden als unkorreliert angesehen, wenn der Koeffizient nahe oder unter Null liegt. Es kann ein beliebiger Schwellenwert innerhalb des Bereichs normalisierter Werte verwendet werden. In einer Ausführungsform wird der Schwellenwert für binäre Gewichtungen verwendet. Oberhalb des Schwellenwerts ist die Gewichtung 1 oder 0. Unterhalb des Schwellenwerts ist die Gewichtung das andere aus 0 oder 1. Gleich dem Schwellenwert kann entweder auf 0 oder 1 abgebildet werden. Das Ergebnis ist eine binäre Maske als Filter. Werte, die zu einem Sparkle gehören, sind die einzigen verbleibenden oder die einzigen entfernten nach dem Filtern mit der binären Maske. Die Werte, die nicht zu einem Sparkle gehören, werden entfernt oder bleiben.
  • In einer Ausführungsform wird dann, wenn der Koeffizient der normalisierten Kreuzkorrelation größer oder gleich einem eingestellten Schwellenwert ε > 0 ist, der Abtastprobenwert mit 0 multipliziert. Wenn der Koeffizient kleiner als der Schwellenwert ε ist, wird der Abtastprobenwert mit 1 multipliziert. Dieser Fall kann zum verstärkten Darstellen von Nierensteinen oder Gallensteinen verwendet werden, wobei ein Sparkle gewünschte Informationen anzeigt. Die Stellen, die zu einer größeren Ähnlichkeit gehören (d. h. mit Bewegung), werden entfernt oder verringert.
  • In einer anderen Ausführungsform wird der Abtastprobenwert dann mit 1 multipliziert, wenn der Koeffizient größer oder gleich einem eingestellten Schwellenwert ε > 0 ist. Wenn der Koeffizient kleiner als der Schwellenwert ε ist, wird der Abtastprobenwert mit 0 multipliziert. Die Bewegungsdaten für Stellen, die zu einer geringeren Ähnlichkeit gehören (d. h. mit Sparkle), werden entfernt oder verringert, während Bewegungsdaten für Bewegungsstellen beibehalten werden.
  • Es kann auch eine andere Gewichtungsabbildung als eine binäre Maske verwendet werden. Die Gewichtungsmatrix kann verallgemeinert werden, anstatt eine Maske von 0 und 1 zu verwenden. Einige Beispiele für eine Abbildung der Ähnlichkeit auf eine Gewichtung sind in 2A bis 2D dargestellt. Der Koeffizient der normalisierten Kreuzkorrelation wird auf Gewichtungen abgebildet. Die Abbildungen der 2A bis 2D sind für eine Verringerung der Geschwindigkeits- oder anderen Bewegungsdaten an Sparkle-Stellen. Durch Anwenden der Gewichtungen der 2A bis 2D werden Stellen mit größerer Ähnlichkeit stärker gewichtet als Stellen mit geringerer Ähnlichkeit, wodurch die Dopplerwerte für die Stellen des Sparkle verringert werden. Es können auch umgekehrte Abbildungen verwendet werden, wobei Stellen mit geringerer Ähnlichkeit stärker gewichtet werden als Stellen mit höherer Ähnlichkeit, um die Dopplerwerte für andere Stellen als Sparkle zu verringern. Es kann auch eine andere lineare oder nichtlineare Abbildung der Ähnlichkeit auf Gewichtungen verwendet werden.
  • Im Schritt 46 wird das gefilterte Farbflussbild angezeigt. Das Ultraschallsystem verarbeitet den gefilterten Datenrahmen, um das Bild zu erzeugen. Es wird eine räumliche Filterung, eine zeitliche Filterung, eine Abtastungsumwandlung oder eine andere Bildverarbeitung durchgeführt. Die Skalarwerte werden auf Anzeigewerte abgebildet, z. B. auf Farbwerte unter Verwendung einer Geschwindigkeitsskala. Das resultierende Bild wird zur Anzeige gepuffert. Die Anzeigewerte werden der Anzeigevorrichtung von dem Puffer bereitgestellt.
  • Es wird ein Farbfluss-(z. B. Doppler-Energie- oder Doppler-Geschwindigkeits-), ein Doppler-Gewebebewegungs- oder ein anderes Bewegungsbild erzeugt. Das Bild kann auch andere Informationen umfassen. Zum Beispiel ist das Bild eine Überlagerung von B-Modus-daten durch die Farbflussdaten. Für Nicht-Gewebe-Stellen oder Stellen, die zu einem ausreichenden Fluss gehören, werden die Farbflussdaten (z. B. Geschwindigkeiten) verwendet, um eine Farbe für die Anzeige zu bestimmen. Für Gewebestellen oder Stellen mit wenig/keinem Fluss werden die B-Modus-Daten verwendet.
  • Das Bild umfasst Bewegungswerte (z. B. Geschwindigkeiten oder Energie), welche verändert worden sind, um das Sparkle-Artefakt oder Bewegungsinformationen, die nicht zu einer rauen Fläche gehören, zu entfernen oder zu verringern. Das Farbflussbild, z. B. ein Geschwindigkeitsbild, wird nach dem Filtern auf der Grundlage der Ähnlichkeit aus den Dopplerwerten erzeugt. Als ein Ergebnis sind in dem Farbflussbild Steine oder andere Ziele mit rauer Fläche ohne andere Bewegungsinformationen hervorgehoben oder es wird ein Fluss mit weniger Sparkle-Artefakt hervorgehoben. Zum Beispiel wird ein Nierenstein- oder Gallensteinbild angezeigt. Als ein anderes Beispiel wird ein Fluidgeschwindigkeitsbild angezeigt.
  • 3A bis 3D und 4A bis 4C zeigen zwei Beispiele. 3A bis 3B und 4A bis 4B zeigen Geschwindigkeitsbilder in Kombination mit B-Modus. Die Geschwindigkeitsinformationen werden farblich abgebildet. Die Geschwindigkeitsteile dieser Bilder repräsentieren Rahmen von Bewegungsdaten. Für 3A und 4A werden für die Bilder Schätzungen aus einer vollständigen Fluss-Abtastprobenanzahl oder einer höheren PRF verwendet. Für 3B und 4B werden für die Bilder Schätzungen aus einer Teilabtastung derselben zurückkehrenden Abtastproben mit einem Faktor 2 verwendet (z. B. Schätzung der Geschwindigkeit unter Verwendung jeder anderen Abtastprobe der vollständigen Fluss-Abtastprobenanzahl). 3B und 4B sind aus Bewegungsdaten, die mit einer geringeren PRF geschätzt werden.
  • 3C zeigt zwei weiße Regionen ohne erkanntes Sparkle. Flussregionen mit Korrelation unterhalb eines Schwellenwerts werden mit 0 gewichtet, was zu der Maske zum Durchlassen eines fehlerfreien Flusses der 3C führt. Durch Anwenden der Maske der 3C auf das Bild der 3A entsteht ein Geschwindigkeitsbild der 3D. In dem Bild ist das Sparkle entfernt oder verringert. Für Nicht-Bewegungsstellen können B-Modus-Informationen bereitgestellt werden.
  • 4A zeigt ein Farbflussbild. Eine binäre Maske (z. B. 4C), welche Stellen durchlässt (z. B. Gewichtung 1), die zu einem Sparkle gehören, und Flussinformationen entfernt, die nicht zu einem Sparkle gehören, wird durch Korrelieren der Bewegungsinformationen der 4A und 4B erzeugt. Durch Anwenden der resultierenden binären Maske auf 4A werden die Farbfluss- oder Bewegungsinformationen für Stellen, die zu einem Sparkle gehören, beibehalten, während Bewegungsinformationen für andere Stellen entfernt werden. Die resultierende 4D zeigt Farbinformationen als Sparkle, welche einen Nierenstein ohne Fluss aus Blut- oder Gewebebewegung repräsentieren. Für nicht-Bewegungsstellen können B-Modus-Informationen bereitgestellt werden.
  • 5 zeigt eine Ausführungsform eines Systems 10 für eine Sparkle-basierte Verarbeitung in Flussbildern (d. h. bei einer Farbfluss-, Gewebebewegungs- oder anderen Bewegungsbildgebung). Das System 10 realisiert das Verfahren der 1 oder ein anderes Verfahren. Durch Schätzen von Bewegungsdaten mit unterschiedlicher PRF können die Bewegungsdaten verglichen werden, um ein Sparkle anzuzeigen. Die Ergebnisse des Vergleichs werden für eine Verarbeitung, Filterung oder andere Bildgebung auf der Grundlage erkannter Sparkle-Stellen verwendet.
  • Das System 10 umfasst einen Sendestrahlformer 12, einen Messaufnehmer 14, einen Empfangsstrahlformer 16, einen Speicher 18, ein Filter 20, eine Fluss-Schätzeinheit 22, einen anderen Speicher 28, einen Prozessor 24 und eine Anzeigevorrichtung 27. Außerdem können andere oder weniger Komponenten bereitgestellt werden. Zum Beispiel umfasst das System einen B-Modus-Detektor. Als ein anderes Beispiel werden die Fluss-Schätzeinheit 22 und der Prozessor 24 ohne die Vorfeldkomponenten wie z. B. den Sende- und Empfangsstrahlformer 12, 16 bereitgestellt. In noch einem anderen Beispiel sind die Speicher 18 und 28 eine Komponente.
  • In einer Ausführungsform ist das System 10 ein medizinisches diagnostisches Ultraschallsystem. In einer alternativen Ausführungsform ist das System 10 ein Computer oder eine Arbeitsplatzstation. In noch einer anderen Ausführungsform ist die Fluss-Schätzeinheit 22 Teil eines medizinischen diagnostischen Ultraschallsystems oder eines anderen medizinischen Bildgebungssystems und der Prozessor 24 ist Teil einer separaten Arbeitsplatzstation oder eines entfernten Systems, aus welchem das Ultraschall-Bildgebungssystem aufgebaut ist.
  • Der Messaufnehmer 14 ist eine Matrix mehrerer Elemente. Die Elemente sind piezoelektrische oder kapazitive Membranelemente. Die Matrix ist als eine eindimensionale Matrix, eine zweidimensionale Matrix, eine 1,5D-Matrix, eine 1,25D-Matrix, eine 1,75D-Matrix, eine Ringmatrix, eine mehrdimensionale Matrix, eine Wobbelmatrix, Kombinationen dieser oder eine beliebige andere derzeit bekannte oder später entwickelte Matrix konfiguriert. Die Messaufnehmerelemente wandeln zwischen akustischer und elektrischer Energie. Der Messaufnehmer 14 ist über einen Sende/Empfangs-Schalter mit dem Sendestrahlformer 12 und dem Empfangsstrahlformer 16 verbunden, in anderen Ausführungsformen können aber auch separate Verbindungen verwendet werden.
  • Der Sende- und Empfangsstrahlformer 12, 16 sind Strahlformer zum Abtasten einer Region des Patienten mit dem Messaufnehmer 14. Der Sendestrahlformer 12 sendet unter Benutzung des Messaufnehmers 14 einen oder mehrere Strahlen zum Abtasten einer Region. Es können Vector®, Sektor-lineare oder andere Abtastformate verwendet werden. Die Empfangslinien und/oder Sendestrahlen werden in der Abtastregion verteilt. Der Empfangsstrahlformer 16 tastet in verschiedenen Tiefen mit den Empfangsstrahlen ab. Durch Abtasten derselben Stellen zu verschiedenen Zeiten wird eine Sequenz für die Flussschätzung erhalten.
  • Bei dem Sendestrahlformer 12 handelt es sich um einen Prozessor, eine Verzögerungseinheit, ein Filter, einen Signalformgenerator, einen Speicher, einen Phasendreher, einen Digital-analog-Wandler, einen Verstärker, Kombinationen dieser oder eine beliebige andere derzeit bekannte oder später entwickelte Sendestrahlformer-Komponente. In einer Ausführungsform erzeugt der Sendestrahlformer 12 digital Enveloppe-Abtastproben. Durch Filterung, Verzögerungen, Phasendrehung, Digital-analog-Wandlung und Verstärkung wird die gewünschte Sendesignalform erzeugt. Es können auch andere Signalformgeneratoren verwendet werden, wie z. B. Schaltimpulsgeber oder Signalformspeicher.
  • Der Sendestrahlformer 12 ist zu mehreren Kanälen zum Erzeugen elektrischer Signale einer Sendesignalform für jedes Element einer Sendeapertur auf dem Messaufnehmer 14 konfiguriert. Die Signalformen sind unipolare, bipolare, gestufte, sinusförmige oder andere Signalformen einer gewünschten Mittenfrequenz oder eines Frequenzbands mit einem, mehreren und/oder einer Bruchzahl an Zyklen. Die Signalformen weisen eine relative Verzögerung und/oder Phase und Amplitude zum Fokussieren der akustischen Energie auf. Der Sendestrahlformer 12 umfasst eine Steuerung zum Ändern einer Apertur (z. B. der Anzahl aktiver Elemente), eines Apodisationsprofils (z. B. Typ oder Schwerpunkt) über die mehreren Kanäle, eines Verzögerungsprofils über die mehreren Kanäle, eines Phasenprofils über die mehreren Kanäle, einer Mittenfrequenz, eines Frequenzbands, einer Signalform, einer Anzahl an Zyklen und/oder von Kombinationen dieser. Auf der Grundlage dieser Strahlformungsparameter wird ein Sendestrahlfokus erzeugt.
  • Bei dem Empfangsstrahlformer 16 handelt es sich um einen Vorverstärker, ein Filter, einen Phasendreher, eine Verzögerungseinheit, einen Summierer, ein Basisbandfilter, einen Prozessor, Puffer, einen Speicher, Kombinationen dieser oder andere derzeit bekannte oder später entwickelte Empfangsstrahlformer-Komponenten. Der Empfangsstrahlformer 16 ist zu mehreren Kanälen zum Empfangen elektrischer Signale konfiguriert, welche Echos oder akustische Energie repräsentieren, die auf den Messaufnehmer 14 treffen. Ein Kanal von jedem der Elemente der Empfangsapertur innerhalb des Messaufnehmers 14 verbindet mit einem Verstärker und/oder einer Verzögerungseinheit. Ein Analog-digital-Wandler digitalisiert das verstärkte Echosignal. Die digitalen Funkfrequenz-Empfangsdaten werden auf eine Basisbandfrequenz demoduliert. Anschließend werden von dem Verstärker und/oder der Verzögerungseinheit beliebige Empfangsverzögerungen, z. B. dynamische Empfangsverzögerungen und/oder Phasendrehungen, angewendet. Ein digitaler oder analoger Summierer kombiniert Daten von verschiedenen Kanälen der Empfangsapertur, um zurückkehrende Abtastproben für einen oder mehrere Empfangsstrahlen zu bilden. Bei dem Summierer handelt es sich um einen einzelnen Summierer oder einen Kaskadensummierer. In einer Ausführungsform ist der Strahlformsummierer dafür konfiguriert, gleichphasige und gegenphasige Kanaldaten in komplexer Weise zu summieren, so dass Phaseninformationen für den gebildeten Strahl bewahrt werden. Alternativ summiert der Strahlformsummierer Amplituden- oder Intensitätsdaten, ohne die Phaseninformationen zu bewahren.
  • Der Empfangsstrahlformer 16 ist dafür konfiguriert, Empfangsstrahlen in Reaktion auf die Sendestrahlen zu bilden. Zum Beispiel empfängt der Empfangsstrahlformer 16 einen, zwei oder mehr Empfangsstrahlen in Reaktion auf jeden Sendestrahl. Die Empfangsstrahlen sind kollinear, parallel mit Abstand oder nichtparallel zu den entsprechenden Sendestrahlen. Der Empfangsstrahlformer 16 gibt räumliche Abtastproben aus, welche verschiedene räumliche Positionen einer abgetasteten Region repräsentieren. Sobald die Kanaldaten strahlgeformt oder auf andere Weise kombiniert werden, um räumliche Positionen entlang den Abtastlinien 11 zu repräsentieren, werden die Daten aus der Kanaldomäne in die Bilddatendomäne umgewandelt. Die Phasendreher, Verzögerungseinheiten und/oder Summierer können für eine parallele Empfangsstrahlformung wiederholt werden. Einer oder mehrere der parallelen Empfangsstrahlformer können Teile von Kanälen gemeinsam benutzen, z. B. die Anfangsverstärkung gemeinsam benutzen.
  • Zum Abbilden einer Bewegung, z. B. Gewebebewegung oder Fluidgeschwindigkeit, werden für jede von mehreren im Wesentlichen gleichen räumlichen Positionen mehrere Sendungen und entsprechende Empfänge durchgeführt. Phasenänderungen zwischen den verschiedenen Empfangsereignissen für jede gegebene Position zeigen die Geschwindigkeit des Gewebes oder Fluids an. Eine Geschwindigkeits-Abtastprobengruppe entspricht mehreren Sendungen für jede von mehreren Abtastlinien 11. Die Anzahl der Male, wie oft eine im Wesentlichen gleiche räumliche Position, z. B. eine Abtastlinie 11, innerhalb einer Geschwindigkeits-Abtastprobengruppe abgetastet wird, ist die Geschwindigkeits- oder Fluss-Abtastprobenanzahl. Die Sendungen für verschiedene Abtastlinien 11, verschiedene Geschwindigkeits-Abtastprobengruppierungen oder verschiedene Bildgebungstypen können verschachtelt sein. Die Zeitdauer zwischen Sendungen an eine im Wesentlichen gleiche Abtastlinie 11 innerhalb der Geschwindigkeits-Abtastprobenanzahl ist das Impulswiederholungsintervall. Aus dem Impulswiederholungsintervall ergibt sich die Impulswiederholungsfrequenz oder umgekehrt.
  • Bei dem Speicher 18 handelt es sich um einen Video-Direktzugriffsspeicher, einen Direktzugriffsspeicher, ein entfernbares Medium (z. B. eine Diskette oder eine CD), eine Festplatte, eine Datenbank, einen Corner-Turning-Speicher oder eine andere Speichereinheit zum Speichern von Daten oder Videoinformationen. In einer Ausführungsform ist der Speicher 18 ein Corner-Turning-Speicher einer Bewegungsparameter-Schätzstrecke. Der Speicher 18 ist dafür konfiguriert, Signale in Reaktion auf mehrere Sendungen entlang einer im Wesentlichen gleichen Abtastlinie zu speichern. Der Speicher 22 ist dafür konfiguriert, Ultraschalldaten zu speichern, die in einem akustischen Gitter, in einem kartesischen Gitter, sowohl in einem kartesischen Koordinatengitter als auch in einem akustischen Gitter formatiert sind, oder Ultraschalldaten, die ein Volumen in einem 3D-Gitter repräsentieren. Es werden die zurückkehrenden Abtastproben der Fluss-Abtastprobenanzahl für jede von mehreren Positionen gespeichert.
  • Bei dem Filter 20 handelt es sich um ein Clutter-Filter, ein Filter mit endlicher Impulsantwort, ein Filter mit unendlicher Impulsantwort, ein analoges Filter, ein digitales Filter, Kombinationen davon oder andere derzeit bekannte oder später entwickelte Filter. In einer Ausführungsform umfasst das Filter 20 einen Mixer zum Verschieben von Signalen auf das Basisband und eine programmierbare Tiefpassfilter-Reaktion zum Entfernen oder Minimieren von Informationen bei Frequenzen, die vom Basisband entfernt sind. In anderen Ausführungsformen ist das Filter 20 ein Tiefpass-, Hochpass- oder Bandpassfilter. Das Filter 20 verringert Geschwindigkeiten von Fluiden oder verringert alternativ den Einfluss von Daten von Gewebe, während die Geschwindigkeitsinformationen von Fluiden beibehalten werden. Das Filter 20 weist eine eingestellte Reaktion auf oder kann programmiert werden, z. B. auf eine Änderung der Operation als eine Funktion der Signalrückführung oder auf ein anderes adaptives Verfahren. In noch einer anderen Ausführungsform sind der Speicher 18 und/oder das Filter 20 Teil der Fluss-Schätzeinheit 22.
  • Die Doppler- oder Fluss-Schätzeinheit 22 ist ein Doppler-Prozessor oder Kreuzkorrelationsprozessor zum Schätzen der Farbflussdaten. In alternativen Ausführungsformen kann eine andere derzeit bekannte oder später entwickelte Einheit zum Schätzen der Geschwindigkeit, Leistung (z. B. Energie) und/oder Varianz aus beliebigen oder verschiedenen Eingabedaten bereitgestellt werden. Die Fluss-Schätzeinheit 22 empfängt mehrere Signale, die zu einer im Wesentlichen gleichen Position gehören, zu verschiedenen Zeiten und schätzt eine Dopplerverschiebungsfrequenz auf der Grundlage einer Änderung oder einer mittleren Änderung der Phase zwischen aufeinanderfolgenden Signalen von derselben Stelle. Aus der Dopplerverschiebungsfrequenz wird die Geschwindigkeit berechnet. Alternativ wird die Dopplerverschiebungsfrequenz als eine Geschwindigkeit verwendet. Die Leistung und die Varianz können ebenfalls berechnet werden.
  • Farbflussdaten (z. B. Geschwindigkeit, Leistung und/oder Varianz) werden für räumliche Positionen in der Abtastregion aus den Abtastproben der strahlgeformten Abtastung geschätzt. Zum Beispiel repräsentieren die Flussdaten mehrere verschiedene Positionen in einer Ebene. Die Farbflussdaten sind Bewegungsdaten für Gewebe und/oder Fluid.
  • Die Fluss-Schätzeinheit 22 kann einen oder mehrere Schwellenwerte anwenden, um ausreichende Bewegungsinformationen zu identifizieren. Zum Beispiel werden Geschwindigkeits- und/oder Leistungsschwellenwerte zum Identifizieren von Geschwindigkeiten angewendet. In alternativen Ausführungsformen wendet ein separater Prozessor oder ein separates Filter Schwellenwerte an. In anderen Ausführungsformen werden die Schwellenwerte nach einer Bewegungsunterdrückung angewendet, z. B. durch den Prozessor 24.
  • Die Fluss-Schätzeinheit 22 gibt Rahmen von Daten aus, welche die Abtastregion zu verschiedenen Zeiten repräsentieren. Die Abtastproben der strahlgeformten Abtastung für eine gegebene Fluss-Abtastprobenanzahl werden verwendet, um für ein Mal zu schätzen. Ein sich bewegendes Fenster mit Überlappen der Daten wird verwendet, um für andere Male zu schätzen. Es werden Geschwindigkeiten für jede Position zu verschiedenen Zeiten ausgegeben.
  • Aus der Abtastung für dieselben Positionen werden zwei oder mehr Rahmen geschätzt. Die Rahmen sind für eine gleiche oder eine verschiedene Zeit. Für sequentielle Abtastung repräsentieren die Positionen zu verschiedenen Zeiten. Wenn die Abtastproben für dieselbe Fluss-Abtastprobenanzahl verwendet werden, um die Rahmen zu schätzen, dann repräsentieren die Rahmen eine gleiche Zeit.
  • Die zwei oder mehr Rahmen werden unterschiedlich geschätzt. Die Anzahl der Abtastproben, die bei der Schätzung für jeden Rahmen verwendet werden ist verschieden. Die PRF oder Fluss-Abtastprobenanzahl bei der Schätzung variiert von Rahmen zu Rahmen. Zum Beispiel werden zum Schätzen der Geschwindigkeit in einem Rahmen 16 Abtastproben erhalten und verwendet. Alle anderen Abtastproben (z. B. 8 Abtastproben) werden verwendet, um die Geschwindigkeit in dem anderen Rahmen zu schätzen. Der Teilabtastungsfaktor für den anderen Rahmen ist 2. Um die Geschwindigkeit der zwei Rahmen zu schätzen, wird eine unterschiedliche Anzahl an Abtastproben verwendet.
  • Die Fluss-Schätzeinheit 22 oder der Prozessor 24 ist dafür konfiguriert, die Bewegungsdaten der zwei oder mehr Rahmen zu normalisieren. Aufgrund einer Teilabtastung kann die Geschwindigkeitsschätzung eine andere sein, z. B. führt die Teilabtastung mit dem Faktor 2 zu den Ergebnissen, dass die Autokorrelation im Vergleich zu dem Rahmen, der aus der vollständigen Fluss-Abtastprobenanzahl geschätzt wird, verdoppelt wird. Die Bewegungsdaten aus der vollständigen Fluss-Abtastprobenanzahl werden zum Normalisieren multipliziert. Alternativ wird keine Multiplikation oder Normalisierung durchgeführt und die Rahmen werden verwendet, obwohl sie unterschiedliche Größenbereiche aufweisen.
  • Bei dem Prozessor 24 handelt es sich um einen digitalen Signalprozessor, einen allgemeinen Prozessor, eine anwendungsspezifische integrierte Schaltung, ein feldprogrammierbares Gate-Array, einen Steuerprozessor, ein digitales Schaltungssystem, ein analoges Schaltungssystem, eine Graphikverarbeitungseinheit, ein Filter, Kombinationen davon oder eine andere derzeit bekannte oder später entwickelte Einheit zum Realisieren von Berechnungen, Algorithmen, Programmierungen oder anderen Funktionen. Der Prozessor 24 arbeitet nach Befehlen, die in dem Speicher 28 oder einem anderen Speicher für die Sparkle-basierte Bildverarbeitung bereitgestellt werden. Es können zusätzliche oder mehrere Prozessoren benutzt werden. Der Prozessor 24 wird durch Software, Firmware und/oder Hardware konfiguriert.
  • Der Prozessor 24 empfängt Farbflussdaten aus der Fluss-Schätzeinheit 22, dem Speicher 28 und/oder einer anderen Quelle. Der Prozessor 24 ist dafür konfiguriert, unter Verwendung der empfangenen Bewegungsdaten zu identifizieren, welche der Stellen der Abtastregion einem Sparkle in den Schätzungen entsprechen. Der Prozessor 24 ist dafür konfiguriert, die Bewegungswerte aus den verschiedenen Rahmen zu vergleichen. Zum Beispiel wird eine räumliche Kreuzkorrelation zwischen den Daten der Rahmen durchgeführt. Die Kreuzkorrelation zeigt ein Maß an Ähnlichkeit zwischen den verschiedenen PRF-Rahmen an. Die Ähnlichkeit wird für jede der mehreren Stellen bestimmt, die in den Rahmen repräsentiert wird. Die Stellen mit weniger Korrelation zwischen Bewegungswerten, die aus unterschiedlicher PRF geschätzt werden, sind Sparkle-Stellen. Die Stellen mit mehr Korrelation sind tatsächliche Bewegung statt eines Sparkle.
  • Der Prozessor 24 identifiziert speziell Sparkle-Stellen. Alternativ oder zusätzlich verwendet der Prozessor 24 die Anzeige des Sparkle, die durch das Ähnlichkeitsmaß geliefert wird, zum Filtern oder für eine andere Bildverarbeitung. Das Ähnlichkeitsmaß kann auf Filtergewichtungen abgebildet werden. Die Gewichtungen werden auf einen Rahmen von Bewegungswerten abgebildet, z. B. auf den Rahmen, der unter Verwendung der höheren Anzahl an Abtastproben geschätzt wird. Durch das Abbilden und die entsprechenden Gewichtungen werden Stellen reduziert, die nicht zu einem Sparkle gehören, oder Stellen reduziert, die zu einem Sparkle gehören. Andere Stellen bleiben gleich, werden verstärkt oder nicht so stark reduziert. Die Gewichtungsmatrix und/oder Sparkle-Stellen können gefiltert werden, z. B. mit einem Tiefpassfilter, um Ausreißerstellen zu entfernen oder Ausreißergewichtungen zu verringern.
  • Der Prozessor 24 oder eine andere Komponente verwendet die gefilterten Bewegungswerte, um ein Bild zu erzeugen. Der Rahmen wird gemäß der Abtastung umgewandelt und farblich abgebildet. Die resultierenden Farbwerte werden zu einem B-Modus-Bild hinzugefügt, z. B. diesem überlagert, oder allein verwendet. Die Farbwerte werden in einen Anzeigepuffer gegeben, um ein Bild auf der Anzeigevorrichtung 27 anzuzeigen.
  • Bei der Anzeigevorrichtung 27 handelt es sich um eine CRT-, LCD- oder Plasma-Anzeigevorrichtung, einen Projektor, einen Monitor, einen Drucker, einen berührungsempfindlichen Bildschirm oder eine andere derzeit bekannte oder später entwickelte Anzeigevorrichtung. Die Anzeigevorrichtung 27 empfängt RGB-, andere Farbwerte oder andere Bewegungswerte und gibt ein Bild aus. Das Bild kann ein Schwarzweiß- oder Farbbild sein. Das Bild repräsentiert die Region des Patienten, die von dem Strahlformer und Messaufnehmer 14 abgetastet wird, wobei jedoch Rücksignale von Sparkle-Stellen entfernt, verringert oder beibehalten werden, während Rücksignale für andere Stellen verringert werden. Die Anzeigevorrichtung 27 zeigt ein Doppler- oder anderes Farbbild aus den gefilterten Bewegungswerten an. Das resultierende Bild kann raue Flächen darstellen, z. B. von Steinen, ohne oder mit weniger Informationen aus Fluidfluss oder Gewebebewegung. Alternativ kann das resultierende Bild Fluidfluss oder Gewebebewegung darstellen, ohne oder mit weniger Informationen aus Sparkle-Artefakten.
  • Bei dem Speicher 28 handelt es sich um einen Video-Direktzugriffsspeicher, einen Direktzugriffsspeicher, ein entfernbares Medium (z. B. eine Diskette oder eine CD), eine Festplatte, eine Datenbank oder eine andere Speichereinheit zum Speichern von Farbfluss- oder anderen Bewegungsdaten. Die gespeicherten Daten liegen in einem Format in Polarkoordinaten oder kartesischen Koordinaten vor. Der Speicher 28 wird vom Prozessor 24 für die verschiedenen Filter-, Erkennungs-, Identifizierungs-, Korrelations-, Rechen- oder anderen Schritte verwendet.
  • Die Befehle zum Realisieren der oben erörterten Prozesse, Verfahren und/oder Techniken werden auf computerlesbaren Speichermedien oder Speichern bereitgestellt, wie z. B. einem Cache, einem Puffer, einem RAM, einem entfernbaren Medium, einer Festplatte oder anderen computerlesbaren Speichermedien. Der Speicher 28 oder andere Speicher speichert die Befehle zur Sparkle-Erkennung bei der Doppler-Bildgebung. Permanente computerlesbare Speichermedien umfassen verschiedene Typen flüchtiger und nichtflüchtiger Speichermedien. Die Funktionen, Schritte oder Aufgaben, die in den Figuren veranschaulicht oder hierin beschrieben werden, werden in Reaktion auf einen oder mehrere Sätze von Befehlen ausgeführt, die in oder auf computerlesbaren Speichermedien gespeichert sind. Die Funktionen, Schritte oder Aufgaben sind von dem speziellen Typs des Befehlssatzes, des Speichermediums, des Prozessors oder der Verarbeitungsstrategie unabhängig und können durch Software, Hardware, integrierte Schaltungen, Firmware, Mikrocode und dergleichen ausgeführt werden, die allein oder in Kombination wirken. In ähnlicher Weise können Verarbeitungsstrategien Mehrprozessorbetrieb, Mehrprogrammbetrieb, parallele Verarbeitung und dergleichen umfassen. In einer Ausführungsform werden die Befehle auf einem entfernbaren Medium zum Lesen durch lokale oder entfernte Systeme gespeichert. In anderen Ausführungsformen werden die Befehle an einer entfernten Stelle für eine Übermittlung durch ein Computernetzwerk oder über Telefonleitungen gespeichert. In wiederum anderen Ausführungsformen werden die Befehle innerhalb eines gegebenen Computers, einer gegebenen CPU, einer gegebenen GPU oder eines gegebenen Systems gespeichert.
  • Obwohl die Erfindung vorstehend unter Bezugnahme auf verschiedene Ausführungsformen beschrieben worden ist, versteht es sich, dass viele Änderungen und Modifikationen vorgenommen werden können, ohne vom Umfang der Erfindung abzuweichen. Die vorstehende detaillierte Beschreibung soll daher als veranschaulichend und nicht als beschränkend angesehen werden und es versteht sich, dass die Idee und der Umfang der vorliegenden Erfindung durch die folgenden Patentansprüche einschließlich aller Äquivalente definiert werden.

Claims (20)

  1. Verfahren zur Sparkle-Artefakt-Erkennung im Farbfluss, wobei das Verfahren umfasst: Erzeugen erster Farbflussdaten, welche Stellen in einem Patienten repräsentieren, wobei die ersten Farbflussdaten mit einer ersten Impulswiederholungsfrequenz erzeugt werden; Erzeugen zweiter Farbflussdaten, welche die Stellen in dem Patienten repräsentieren, wobei die zweiten Farbflussdaten mit einer zweiten Impulswiederholungsfrequenz erzeugt werden; Bestimmen eines Maßes an Ähnlichkeit zwischen den ersten und zweiten Farbflussdaten für jede der Stellen; Filtern eines Farbflussbildes, wobei das Filtern auf dem Maß an Ähnlichkeit als Funktion der Stelle basiert; und Anzeigen des gefilterten Farbflussbildes.
  2. Verfahren nach Anspruch 1, wobei das Erzeugen der ersten und zweiten Farbflussdaten aufeinander folgendes Senden mit der ersten bzw. zweiten Impulswiederholungsfrequenz umfasst.
  3. Verfahren nach Anspruch 1, wobei das Erzeugen der ersten Farbflussdaten Senden mit der ersten Impulswiederholungsfrequenz umfasst und wobei das Erzeugen der zweiten Farbflussdaten Teilabtasten von Rücksignalen der Aussendung umfasst.
  4. Verfahren nach Anspruch 3, wobei das Erzeugen des ersten Farbflusses ferner Schätzen einer Geschwindigkeit aus den Rücksignalen und Multiplizieren der geschätzten Geschwindigkeiten mit zwei umfasst und wobei das Teilabtasten Teilabtasten mit einem Faktor zwei umfasst.
  5. Verfahren nach einem der Ansprüche 1 bis 4, wobei das Bestimmen Durchführen einer normalisierten Kreuzkorrelation zwischen den ersten und zweiten Farbflussdaten umfasst.
  6. Verfahren nach einem der Ansprüche 1 bis 5, wobei das Filtern Abbilden von Gewichtungen als Funktion des Maßes an Ähnlichkeit je Stelle und Gewichten des Farbflussbildes als Funktion der Stellen mit den Gewichtungen umfasst.
  7. Verfahren nach einem der Ansprüche 1 bis 6, wobei das Filtern Gewichten von Stellen mit geringerer Ähnlichkeit stärker als Stellen mit größerer Ähnlichkeit umfasst und wobei das Anzeigen Anzeigen eines Nieren- oder Gallensteinbildes umfasst.
  8. Verfahren nach einem der Ansprüche 1 bis 7, wobei das Filtern Gewichten von Stellen mit größerer Ähnlichkeit stärker als Stellen mit geringerer Ähnlichkeit umfasst und wobei das Anzeigen Anzeigen eines Fluidgeschwindigkeitsbildes umfasst.
  9. Verfahren nach einem der Ansprüche 1 bis 8, wobei das Filtern Anwenden von binären Gewichtungen als Funktion des Maßes an Ähnlichkeit umfasst.
  10. Permanentes computerlesbares Speichermedium, auf welchem Daten gespeichert sind, die Befehle repräsentieren, die von einem programmierten Prozessor zur Sparkle-Erkennung bei der Doppler-Bildgebung ausführbar sind, wobei das Speichermedium Befehle umfasst zum: Abtasten mehrerer Stellen mit Ultraschall, wobei das Abtasten jeder der Stellen eine Fluss-Abtastprobenanzahl von zurückkehrenden Abtastproben liefert; Schätzen erster Dopplerwerte für die Stellen unter Verwendung der zurückkehrenden Abtastproben der Fluss-Abtastprobenanzahl; Schätzen zweiter Dopplerwerte für die Stellen unter Verwendung einer Teilabtastung der zurückkehrenden Abtastproben der Fluss-Abtastprobenanzahl; Multiplizieren der ersten Dopplerwerte mit einem Faktor, der eine Funktion der Teilabtastung ist; Korrelieren der zweiten Dopplerwerte mit den multiplizierten ersten Dopplerwerten und Erkennen des Sparkle aus Ergebnissen der Korrelation.
  11. Permanentes computerlesbares Speichermedium nach Anspruch 10, wobei das Schätzen der ersten Dopplerwerte Schätzen erster Geschwindigkeiten mit allen der zurückkehrenden Abtastproben der Fluss-Abtastprobenanzahl umfasst und wobei das Schätzen der zweiten Dopplerwerte Schätzen zweiter Geschwindigkeiten mit der Hälfte oder weniger der zurückkehrenden Abtastproben der Fluss-Abtastprobenanzahl umfasst.
  12. Permanentes computerlesbares Speichermedium nach Anspruch 11, wobei der Faktor zwei ist, wobei das Multiplizieren Multiplizieren der ersten Geschwindigkeiten mit zwei umfasst und wobei das Schätzen der zweiten Geschwindigkeiten Schätzen mit allen anderen der zurückkehrenden Abtastproben umfasst.
  13. Permanentes computerlesbares Speichermedium nach einem der Ansprüche 10 bis 12, wobei das Korrelieren normalisiertes Kreuzkorrelieren umfasst.
  14. Permanentes computerlesbares Speichermedium nach einem der Ansprüche 10 bis 13, wobei das Erkennen Identifizieren der Stellen mit weniger Korrelation als andere Stellen umfasst.
  15. Permanentes computerlesbares Speichermedium nach einem der Ansprüche 10 bis 14, ferner umfassend Verringern der ersten Dopplerwerte für die Sparkle-Stellen und Erzeugen eines Dopplerbildes aus den ersten Dopplerwerten nach dem Verringern.
  16. Permanentes computerlesbares Speichermedium nach einem der Ansprüche 10 bis 15, ferner umfassend Verringern der ersten Dopplerwerte für die anderen Stellen als das Sparkle und Erzeugen eines Dopplerbildes aus den ersten Dopplerwerten nach dem Verringern.
  17. System zur Sparkle-basierten Verarbeitung in Flussbildern, wobei das System umfasst: einen Messaufnehmer und Strahlformer zum Abtasten einer Abtastregion; eine Doppler-Schätzeinheit, welche dafür konfiguriert ist, aus der Abtastung erste Bewegungswerte, welche Stellen der Abtastregion repräsentieren, und zweite Bewegungswerte zu schätzen, welche die Stellen der Abtastregion repräsentieren, wobei die ersten Bewegungswerte mit einer anderen Anzahl an Abtastproben aus der Abtastung geschätzt werden als die zweiten Bewegungswerte; und einen Prozessor, der dafür konfiguriert ist, anhand eines Vergleichs der ersten und zweiten Bewegungswerte zu identifizieren, welche der Stellen Sparkle aufweisen.
  18. System nach Anspruch 17, wobei die Doppler-Schätzeinheit dafür konfiguriert ist, die ersten und zweiten Bewegungswerte als Geschwindigkeiten zu schätzen, wobei die andere Anzahl eine andere Impulswiederholungsfrequenz ist.
  19. System nach Anspruch 17 oder 18, wobei der Prozessor dafür konfiguriert ist, anhand des Vergleichs, der eine Kreuzkorrelation ist, zu identifizieren, dass die Stellen mit dem ersten und zweiten Bewegungswert, die eine geringere Korrelation als andere Stellen aufweisen, die Stellen mit dem Sparkle sind.
  20. System nach einem der Ansprüche 17 bis 19, welches ferner eine Anzeigevorrichtung umfasst, die dafür konfiguriert ist, ein Dopplerbild aus den ersten Bewegungswerten anzuzeigen, die als eine Funktion der Stellen mit Sparkle gefiltert sind.
DE102016116658.8A 2015-09-10 2016-09-06 Sparkleartefakt-erkennung im ultraschall-farbfluss Pending DE102016116658A1 (de)

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
US14/850,707 US11096671B2 (en) 2015-09-10 2015-09-10 Sparkle artifact detection in ultrasound color flow
US14/850,707 2015-09-10

Publications (1)

Publication Number Publication Date
DE102016116658A1 true DE102016116658A1 (de) 2017-03-16

Family

ID=58160924

Family Applications (1)

Application Number Title Priority Date Filing Date
DE102016116658.8A Pending DE102016116658A1 (de) 2015-09-10 2016-09-06 Sparkleartefakt-erkennung im ultraschall-farbfluss

Country Status (5)

Country Link
US (1) US11096671B2 (de)
JP (1) JP6385992B2 (de)
KR (1) KR102041433B1 (de)
CN (1) CN106529561B (de)
DE (1) DE102016116658A1 (de)

Families Citing this family (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP6604272B2 (ja) * 2016-06-06 2019-11-13 コニカミノルタ株式会社 超音波診断装置及び超音波信号処理方法
RU2728836C2 (ru) 2016-10-28 2020-07-31 Ппг Индастриз Огайо, Инк. Покрытия для увеличения расстояния обнаружения до объекта, обнаруживаемого с помощью электромагнитного излучения ближнего инфракрасного диапазона
US11419584B2 (en) 2017-06-28 2022-08-23 Tokitae, LLC Devices, systems, and methods for diagnosis of pulmonary conditions through detection of b-lines in lung sonography
JP7285826B2 (ja) * 2017-08-21 2023-06-02 コーニンクレッカ フィリップス エヌ ヴェ 肺超音波検査におけるbラインの検知、提示及び報告
JP7052313B2 (ja) * 2017-11-20 2022-04-12 コニカミノルタ株式会社 超音波信号処理装置、超音波診断装置、および、超音波信号処理方法
CN111629670B (zh) * 2017-12-20 2024-04-05 韦拉索恩股份有限公司 用于超声***的回波窗口伪影分类和视觉指示符
US11561329B2 (en) 2019-01-07 2023-01-24 Ppg Industries Ohio, Inc. Near infrared control coating, articles formed therefrom, and methods of making the same
EP4031055A1 (de) * 2019-09-20 2022-07-27 Mayo Foundation for Medical Education and Research Nichtmetallische ultraschallnachweisbare markierungen
CN111436971B (zh) * 2020-04-20 2022-06-10 汕头市超声仪器研究所股份有限公司 一种超声混响伪像抑制方法
CN112256474B (zh) * 2020-09-11 2022-12-13 安徽中科新辰技术有限公司 一种针对软件实时同步的双机热备份***

Family Cites Families (16)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5961462A (en) * 1998-05-18 1999-10-05 Atl Ultrasound Ultrasonic doppler imaging at high frame rates of display
US6123672A (en) 1998-12-31 2000-09-26 General Electric Company Color flow imaging for enhancing segmentation and flow dynamics
US7288068B2 (en) 2003-12-15 2007-10-30 Siemens Medical Solutions Usa, Inc. Automatic optimization for ultrasound medical imaging
US7946990B2 (en) * 2005-09-30 2011-05-24 Siemens Medical Solutions Usa, Inc. Ultrasound color flow imaging at high frame rates
JP5443760B2 (ja) * 2005-11-10 2014-03-19 ブラッコ・シュイス・ソシエテ・アノニム 流動力学解析に基づく医療画像アプリケーションにおける固定化した造影剤の検出
JP2008171472A (ja) * 2007-01-09 2008-07-24 Fuji Electric Device Technology Co Ltd 情報記録媒体用ガラス基板の製造方法、磁気記録ディスクの製造方法および磁気記録ディスク
US20100185085A1 (en) * 2009-01-19 2010-07-22 James Hamilton Dynamic ultrasound processing using object motion calculation
US8235905B2 (en) 2009-05-26 2012-08-07 General Electric Company System and method for automatic ultrasound image optimization
CN102113900B (zh) * 2010-01-05 2015-07-15 深圳迈瑞生物医疗电子股份有限公司 彩色血流动态帧相关方法和装置
JP5719017B2 (ja) 2010-04-22 2015-05-13 ザ ユニバーシティ オブ ワシントン スルー イッツ センター フォーコマーシャライゼーションThe University Of Washington Through Its Center For Commercialization 結石を検出し、その除去を促進する超音波ベースの方法及び装置
US8684933B2 (en) * 2010-08-17 2014-04-01 Imsonic Medical, Inc. Handheld ultrasound color flow imaging system with mechanically scanned, mechanically focused multi-element transducers
CN103718365B (zh) 2011-08-09 2016-02-17 丰田自动车株式会社 电池以及电池的制造方法
CN103123721B (zh) 2011-11-17 2016-04-27 重庆海扶医疗科技股份有限公司 一种实时减少图像中伪影的方法以及装置
US9420997B2 (en) * 2012-06-14 2016-08-23 Siemens Medical Solutions Usa, Inc. Motion artifact suppression in ultrasound diagnostic imaging
US9320498B2 (en) 2014-02-17 2016-04-26 Siemens Medical Solutions Usa, Inc. Twinkle artifact suppression in ultrasound color flow
US20150320384A1 (en) * 2014-05-06 2015-11-12 University Of Washington Methods and Systems for Detecting an Object in a Subject with Ultrasound

Also Published As

Publication number Publication date
CN106529561A (zh) 2017-03-22
KR20170035796A (ko) 2017-03-31
JP6385992B2 (ja) 2018-09-05
US20170071577A1 (en) 2017-03-16
CN106529561B (zh) 2020-03-17
JP2017051617A (ja) 2017-03-16
US11096671B2 (en) 2021-08-24
KR102041433B1 (ko) 2019-11-06

Similar Documents

Publication Publication Date Title
DE102016116658A1 (de) Sparkleartefakt-erkennung im ultraschall-farbfluss
DE102016100367B4 (de) Spärliche Verfolgung in Schallstrahlintensitätsimpuls-Bildgebung
DE10058449B4 (de) Verfahren und Vorrichtung zur Visualisierung einer Bewegung bei einer Ultraschallflussabbildung unter Verwendung einer kontinuierlichen Datenerfassung
US11391828B2 (en) Methods and systems for filtering ultrasound image clutter
DE60025438T2 (de) Verfahren und Gerät zum Einstellen eines interessanten Gebiets in einer Abbildung
DE10141516B4 (de) Medizinische Diagnostische Ultraschallabbildungsverfahren für ein erweitertes Betrachtungsfeld
DE102012108353B4 (de) Klassifikationsvorverarbeitung in der medizinischen ultraschallscherwellenbildgebung
DE102013011510A1 (de) Automatische Doppler-Gate-Positionierung in der Spektraldopplersonografie
DE10058452B4 (de) Verfahren und Gerät zur Bewegungsdarstellung bei Ultraschall-Fluss-Bilddarstellung unter Verwendung von Paketdatenerfassung
DE102017211895A1 (de) Gewebecharakterisierung im medizinischen diagnostischen Ultraschall
DE102016108737A1 (de) Wissensbasierte Ultraschallbildverbesserung
DE102005029564A1 (de) Verfahren und System zum Abschätzen einer Zeitverzögerung für die Verwendung in der Ultraschallbildgebung
DE102016104263A1 (de) Kontinuierlich ausgerichtete verbesserte Ultraschallbildgebung eines Teilvolumens
DE102016107222A1 (de) Erfassungssteuerung für Gemischtmodus-Ultraschallbildgebung
DE102012101395A1 (de) Mehrstrahliger Spektral-Doppler in der medizinischen diagnostischen Ultraschallbildgebung
DE10306806A1 (de) Verfahren und Vorrichtung zur automatischen Regelung einer Spektraldopplerabbildung
DE102005034697A1 (de) Kontrastmittelbildgebung mit einer (Kontrast)mittelspezifischen Ultraschalldetektion
DE102011114943A1 (de) Bestimmen von Scherwelleninformationen in der medizinischen Ultraschallbildgebung
DE102016114783A1 (de) Adaptive Bewegungsschätzung bei Schallstrahlungskraftbildgebung
DE102017206463A1 (de) Frequenz-compounding bei der elastizitätsbildgebung
DE112015007245T5 (de) Zusammengesetzte Ultraschall-Strömungsbildgebung mit hoher Auflösung
DE102017205566A1 (de) Beugungsquellen-kompensation bei der medizinisch-diagnostischen viskoelastischen ultraschallbildgebung
US8343054B1 (en) Methods and apparatus for ultrasound imaging
DE102016113783B4 (de) Verbesserte Doppler-Bildgebung
DE102013002065B4 (de) Visualisierung von zugehörigen lnformationen bei der Ultraschall-Scherwellenbildgebung

Legal Events

Date Code Title Description
R012 Request for examination validly filed
R016 Response to examination communication
R082 Change of representative

Representative=s name: BALS & VOGEL PATENTANWAELTE PARTG MBB, DE