DE60026239T2 - Verfahren und Anordnung zur Spektral- Doppler- Bilderzeugung mit adaptivem Zeitbereichwandfilter - Google Patents

Verfahren und Anordnung zur Spektral- Doppler- Bilderzeugung mit adaptivem Zeitbereichwandfilter Download PDF

Info

Publication number
DE60026239T2
DE60026239T2 DE60026239T DE60026239T DE60026239T2 DE 60026239 T2 DE60026239 T2 DE 60026239T2 DE 60026239 T DE60026239 T DE 60026239T DE 60026239 T DE60026239 T DE 60026239T DE 60026239 T2 DE60026239 T2 DE 60026239T2
Authority
DE
Germany
Prior art keywords
frequency
clutter
wall
low
spectral
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.)
Expired - Fee Related
Application number
DE60026239T
Other languages
English (en)
Other versions
DE60026239D1 (de
Inventor
Larry Y.L. Waukesha Mo
Jr. Richard Marion Richfield Kulakowski
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.)
General Electric Co
Original Assignee
General Electric Co
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 General Electric Co filed Critical General Electric Co
Application granted granted Critical
Publication of DE60026239D1 publication Critical patent/DE60026239D1/de
Publication of DE60026239T2 publication Critical patent/DE60026239T2/de
Anticipated expiration legal-status Critical
Expired - Fee Related legal-status Critical Current

Links

Classifications

    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01S—RADIO 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/00—Systems using the reflection or reradiation of acoustic waves, e.g. sonar systems
    • G01S15/88—Sonar systems specially adapted for specific applications
    • G01S15/89—Sonar systems specially adapted for specific applications for mapping or imaging
    • G01S15/8906—Short-range imaging systems; Acoustic microscope systems using pulse-echo techniques
    • G01S15/8979—Combined Doppler and pulse-echo imaging systems
    • G01S15/8981—Discriminating between fixed and moving objects or between objects moving at different speeds, e.g. wall clutter filter
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B8/00—Diagnosis using ultrasonic, sonic or infrasonic waves
    • A61B8/06—Measuring blood flow

Landscapes

  • Physics & Mathematics (AREA)
  • Engineering & Computer Science (AREA)
  • Radar, Positioning & Navigation (AREA)
  • Remote Sensing (AREA)
  • Acoustics & Sound (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • General Physics & Mathematics (AREA)
  • Ultra Sonic Daignosis Equipment (AREA)
  • Measurement Of Velocity Or Position Using Acoustic Or Ultrasonic Waves (AREA)

Description

  • Die Erfindung betrifft Ultraschalldiagnosesysteme, die mittels Spektral-Doppler-Techniken die Geschwindigkeit eines Blutstroms messen. Insbesondere betrifft die Erfindung eine kontinuierliche Anzeige derartiger Daten, einschließlich maximaler und mittlerer Blutstromgeschwindigkeiten.
  • Auf dem Dopplereffekt basierende Ultraschallscanner zur Blutstromerfassung sind hinlänglich bekannt. Solche Systeme arbeiten, indem sie ein Ultraschalltransducerarray dazu veranlassen, Ultraschallwellen in das Objekt abzustrahlen und Ultraschallechos zu empfangen, die von dem Objekt rückgestreut werden. Um einen Blutstrom zu messen, werden rückkehrende Ultraschallwellen mit einem Frequenzreferenzwert verglichen, um Frequenzverschiebungen zu ermitteln, die den zurückkehrenden Wellen aufgrund von sich bewegenden Objekten, z.B. durch die Gefäßwände und die innerhalb des Blutgefäßes vorhandenen roten Blutkörperchen, aufgeprägt sind. Diese Frequenzverschiebungen lassen sich in Bewegungsgeschwindigkeiten umsetzen.
  • In Ultraschallscannern nach dem Stand der Technik wird der Dopplerkurvenverlauf der gepulsten oder kontinuierlichen Welle berechnet und als ein Grauskalenspektrogramm der gegen die Zeit aufgetragenen Geschwindigkeit in Echtzeit angezeigt, wobei die Intensität (oder Farbe) der Grauskala durch die spektrale Leistung moduliert ist. Die Daten für jede Spektrallinie basieren auf mehreren Frequenzdatenabschnitten für unterschiedliche Frequenzintervalle, wobei die Daten des Leistungsspektrums in jedem Paket für eine entsprechende Spektrallinie in einem entsprechenden Pixel einer entsprechenden Pixelspalte auf dem Anzeigemonitor an gezeigt werden. Jede Spektrallinie repräsentiert einen momentanen Blutstrommesswert.
  • Im Falle des herkömmlichen Spektral-Dopplermodus wird ein Ultraschalltransducerarray dazu veranlasst, einen Ultraschallsendeimpuls auszusenden, der wiederholt mit einer Impulswiederholfrequenz (PRF = Pulse Repetition Frequency) abgegeben wird. Die PRF liegt gewöhnlich im kHz-Bereich. Die Hochfrequenz-(HF)-Echosignale werden durch die Wandlerelemente erfasst und anschließend durch einen Strahlformer zu einem Empfangsstrahl geformt. Im Falle eines digitalen Systems wird das aufsummierte HF-Signal jeder Impulsabgabe durch einen Demodulator in dessen In-Phase- und Quadratur-(I/Q)-Komponenten demoduliert. Die I/Q-Komponenten werden über ein spezielles Zeitintervall integriert (aufsummiert) und anschließend abgetastet. Das Summierungsintervall und die Sendeimpulslänge definieren gemeinsam die Länge des durch den Benutzer spezifizierten Abtastvolumens. Dieser sogenannte "Summierungs- und Ablage"-Betrieb liefert auf effiziente Weise das von dem Abtastvolumen rückgestreute Dopplersignal. Auf das Dopplersignal wird ein Wandfilter angewandt, das ein Hochpassfilter ist, das in dem Signal jede Clutter-Störung abweist, die einem stationären oder sich sehr langsam bewegenden Gewebe entspricht, z.B. ein Abschnitt der Wand (Wände) eines Gefäßes, der sich möglicherweise innerhalb des Abtastvolumens befindet. Das gefilterte Ausgangssignal wird sodann an eine Spektralanalyseeinrichtung ausgegeben, die gewöhnlich die komplexe schnelle Fourier-Transformation (FFT) an einem sich bewegenden Zeitfenster von 64 bis 256 Abtastwerten ausführt. Die Datenabtastwerte innerhalb eines FFT-Analysezeitfensters werden im Folgenden als ein FFT-Paket bezeichnet. Jedes FFT-Leistungsspektrum wird komprimiert und anschließend mittels einer Grauskalenabbildung auf dem Monitor als eine einzelne Spektrallinie in einem speziellen Zeitpunkt in dem Spektrogramm der Doppler-Geschwindigkeit (Frequenz) gegenüber der Zeit angezeigt.
  • Gewöhnlich werden die I- und Q-Komponenten des Dopplersignals voneinander unabhängig durch identische Wandfilter gefiltert, die entweder als FIR- oder IIR-Filter verwirklicht sein können. Für eine strenge Abweisung niederfrequenter Clutter-Störungen ist ein schmales Übergangsband in der Frequenzantwort des Filters erforderlich. Typischerweise wird die Wandfilter-Eckfrequenz über eine Bedienpultsteuertaste manuell ausgewählt. Für gewöhnlich wird die Wandfilter-Eckfrequenz erhöht, wenn in dem spektralen Bild eine helle, niederfrequente Clutter-Störung beobachtet wird. Bei jeder Änderung der Vorgabewerte der Wandfilter-Grenzfrequenz wird ein entsprechender Satz von Filterkoeffizientenwerten aus einer Referenztabelle (LUT) ausgelesen und in die Wandfilter geladen. Siehe beispielsweise das US-Patent 5 494 037, das ein Ultraschallsystem beschreibt, zu dem Wandfilter gehören, die dazu dienen, Clutter-Störungen in den I- und Q-Komponenten abzuweisen.
  • Der Hauptnachteil des Ansatzes einer manuellen Filterwahl nach dem Stand der Technik basiert darauf, dass das Wandfilter nach einer Voreinstellung der Eckfrequenz unveränderlich bleibt, auch wenn sich die Clutter-Störungsfrequenz und Bandbreite möglicherweise aufgrund einer radialen und/oder seitlichen Bewegung der Gefäßwände über den Herzzyklus hinweg zeitlich ändern. Dementsprechend ist die ausgewählte Filter-Eckfrequenz meist lediglich für einen kleinen Abschnitt des Herzzyklus optimal.
  • Die vorliegende Erfindung betrifft ein Verfahren und eine Vorrichtung zum adaptiven Wand-(Hochpass)-Filtern, die die Nachteile des eine manuelle Filterwahl verwendenden Ansatzes aus dem Stand der Technik beseitigen. Gemäß dem bevorzugten Ausführungsbeispiel der Erfindung wird die Wandfilter-Eckfrequenz automatisch ausgewählt, was den Bedienungskomfort und die Effizienz verbessert. Ein weiterer Vorteil ist, dass es ermöglicht ist, für jedes neue FFT-Paket eine maßgeschneiderte Wandfilter-Eckfrequenz einzusetzen.
  • Gemäß dem bevorzugten Ausführungsbeispiel der Erfindung wird niederfrequente Clutter-Störung in den Doppler-I/Q-Daten vor der FFT-Verarbeitung entfernt. Auf die I/Q-Daten wird ein Tiefpassfilter angewandt, dessen Eckfrequenz bei der für die aktuelle Doppleranwendung antizipierten höchsten Clutter-Störungsfrequenz (beispielsweise 40% der PRF) eingestellt ist. Das Tiefpassfilter weist Fließfrequenzkomponenten ab, die über dem Clutter-Störungsfrequenzbereich liegen. Die Gesamtleistung des Tiefpassfilterausgangssignals, d.h. die Summe von (In 2 + Qn 2) über die Paketgröße M der FFT (oder einer Fraktion der FFT) wird anschließend berechnet.
  • Gemäß den bevorzugten Ausführungsbeispielen wird ein Systemrauschmodell verwendet, um die mittlere Systemrauschleistung in dem Tiefpassfilterausgangssignal vorauszuberechnen. Die durch das Systemrauschmodell vorherberechnete mittlere Systemrauschleistung stellt einen Störschwellenwert bereit, um den Grad der in dem aktuellen FFT-Paket vorhandenen Clutter-Leistung zu messen. Falls keine wesentliche Clutter-Störung vorliegt, wird eine Wandfilteraus wahllogik automatisch aus einer Filterkoeffizienten-LUT die niedrigste Wandfilter-Eckfrequenz auswählen. Falls wesentliche Clutter-Störungsleistung in dem FFT-Paket vorhanden ist, fährt der Algorithmus mit der Berechnung des Mittelwerts und der Varianz der Clutter-Störungsfrequenz über das FFT-Paket hinweg fort. Der geschätzte Mittelwert und die Varianz der Clutter-Störungsfrequenz werden anschließend an die Filterauswahllogik ausgegeben, die die am besten geeignete Filter-Eckfrequenz für das aktuelle Clutter-Signal auswählt. Um eine zu starke Schwankung der Wandfilter-Eckfrequenz von einem FFT-Paket zum nächsten zu vermeiden, kann eine beliebige Persistenzfunktion auf die vorgeschriebene Wandfilter-Eckfrequenz angewandt werden. Nachdem die neue optimale Filter-Eckfrequenz ausgewählt ist, verläuft die übrige Verarbeitung wie im Falle einer herkömmlichen Doppler-Wandfilterung.
  • Dem Fachmann sollte klar sein, dass das Verfahren der Erfindung als Hardware (z.B. in Form eines digitalen Signalverarbeitungschips) und/oder als Software verwirklicht werden kann.
  • Die Erfindung wird nun anhand von Beispielen mit Bezug auf die Zeichnungen eingehender beschrieben:
  • 1 zeigt ein Blockschaltbild der grundlegenden Signalverarbeitungskette in einem herkömmlichen Spektral-Doppler-Bildgebungssystem.
  • 2 zeigt in einem Blockdiagramm ein Wandfilter und zugeordnete Bauelemente eines Spektral-Doppler-Bildgebungssystems der in 1 dargestellten Bauart.
  • 3 zeigt das bevorzugte Ausführungsbeispiel der Erfindung anhand eines Flussdiagramms.
  • 1 zeigt im Wesentlichen ein typisches digitales Echtzeit-Ultraschallbildgebungssystem, das einen Spektral-Doppler-Bildgebungsmodus aufweist. Ein Ultraschalltransducerarray 2 wird durch einen Sender eines Strahlformers 4 veranlasst, Ultraschallstrahlen auszusenden, die auf eine gewünschte Sendefokusposition fokussiert sind. Der Sender stellt einen Ultraschallsendeimpuls bereit, der mit einer Impulswiederholfrequenz (PRF = Pulse Repetition Frequency) wiederholt abgegeben wird. Die PRF liegt gewöhnlich im kH-Bereich. Die HF-Echosignale werden durch die Wandlerelemente erfasst und anschließend durch einen Empfänger eines Strahlformers 4 zu einem Empfangsstrahl geformt. Im Falle eines digitalen Systems wird das aufsummierte (strahlgeformte) HF-Signal jeder Impulsabgabe durch einen Demodulator 6 in dessen In-Phase- und Quadratur-(I/Q)-Komponenten demoduliert. Die I/Q-Komponenten werden über ein spezielles Zeitintervall integriert (aufsummiert) und anschließend durch einen "Summierungs- und Ablage"-Block 8 abgetastet. Das Summierungsintervall und die Sendeimpulslänge definieren gemeinsam die Länge des durch den Benutzer spezifizierten Abtastvolumens. Dieser "Summierungs- und Ablage"-Betrieb liefert auf effiziente Weise das von dem Abtastvolumen rückgestreute Dopplersignal. Das Dopplersignal wird durch ein Wandfilter 10 geleitet, das jede Clutter-Störung in dem Signal abweist, die einem stationären oder sich sehr langsam bewegenden Gewebe entspricht. Das gefilterte Ausgangssignal wird sodann an eine Spektralanalyseeinrichtung 12 ausgegeben, die gewöhnlich die schnelle Fourier-Transformation (FFT) an einem sich bewegenden Zeitfenster von 64 bis 256 Abtastwerten ausführt. Jedes FFT-Leistungsspektrum wird (in Block 14) komprimiert und (in Block 16) auf eine Grauskala abgebildet, um als eine einzelne Spektrallinie zu einem speziellen Zeitpunkt in dem Spektrogramm der Doppler-Geschwindigkeit (Frequenz) gegenüber der Zeit auf einem Monitor 18 angezeigt zu werden.
  • In einem typischen Spektral-Doppler-System werden die I- und Q-Komponenten des Dopplersignals voneinander unabhängig durch identische Wandfilter gefiltert, die entweder als FIR- oder IIR-Filter verwirklicht sein können. Für eine strenge Abweisung niederfrequenter Clutter-Störungen ist ein schmales Übergangsband in der Frequenzantwort des Filters erforderlich. IIR-Filter werden im Allgemeinen als vorteilhafter erachtet, da eine Verwirklichung mittels FIR eine große Filterlänge verlangt. In herkömmlichen Spektral-Doppler-Systemen wird gewöhnlich ein IIR-Hochpassfilter verwendet, und zwar häufig in Form einer Kaskade von drei oder vier Stufen zweiter Ordnung. Das Eingangssignal x(n) und das Ausgangssignal y(n) jeder Stufe zweiter Ordnung genügen der folgenden Differenzengleichung: y(n) = a1y(n – 1) + a2y(n – 2) + b0x(n) + b1x(n – 1) + b2x(n – 2)
  • Für eine vorgegebene Eckfrequenz können die Filterkoeffizienten {a1, a2, b0, b1, b2} für jede Stufe zweiter Ordnung mittels Standard-Filterentwurfsformeln vorausberechnet und in einer LUT 22 gespeichert werden (siehe 2).
  • Typischerweise wird die Wandfilter-Eckfrequenz über eine Bedienpultsteuertaste 20 manuell ausgewählt. Für gewöhnlich wird die Wandfilter-Eckfrequenz erhöht, wenn in dem spektralen Bild eine helle, niederfrequente Clutter- Störung beobachtet wird. Bei jeder Änderung der Vorgabewerte der Wandfilter-Grenzfrequenz wird ein entsprechender Satz von Filterkoeffizientenwerten aus der LUT 22 ausgelesen und in die Wandfilter 10 geladen. Um Einschwingrauschen zu minimieren, können die Werte von y(n – 1) und y(n – 2) für jede Filterstufe unmittelbar nachdem der neue Satz von Filterkoeffizienten geladen ist als Null angenommen werden.
  • Die vorliegende Erfindung stellt eine Verbesserung gegenüber den in 2 gezeigten manuellen Wandfilter-Eckfrequenztechniken dar. Gemäß dem bevorzugten Ausführungsbeispiel der Erfindung wird eine niederfrequente Clutter-Störung in den Doppler I/Q-Daten vor der FFT-Verarbeitung entfernt. Wie in 3 gezeigt, wird auf die I/Q-Daten ein Tiefpassfilter (LPF = Low-Pass Filter) 26 angewandt, dessen Eckfrequenz bei der für die aktuelle Doppleranwendung antizipierten höchsten Clutter-Störungsfrequenz (beispielsweise 40% der PRF) eingestellt ist. Das Tiefpassfilter 26 weist Fließfrequenzkomponenten ab, die über dem Clutter-Störungsfrequenzbereich liegen. Die Gesamtleistung des Tiefpassfilterausgangssignals, d.h. die Summe von (In 2 + Qn 2) über die Paketgröße M der FFT (oder einer Fraktion der FFT) wird anschließend berechnet (Schritt 28).
  • Gemäß den bevorzugten Ausführungsbeispielen wird ein Systemrauschmodell 30 verwendet, um die mittlere Leistung des Systemrauschpegels innerhalb des Durchlassbands des Tiefpassfilters vorauszuberechnen. In dem am meisten bevorzugten Ausführungsbeispiel setzt das Modell einen insgesamt digitalen Scanner voraus, dessen Systemrauschen hauptsächlich auf den in jedem Empfangskanal in dem Strahlformer vorhandenen Vorverstärker zurückzuführen ist. Das Johnsonsche Vorverstärkerrauschen wird häufig als eine Effektiv spannung pro Hz1/2 bei Raumtemperatur spezifiziert (z.B. 10 nV/Hz1/2). Auf diese Weise sollte sich mit Kenntnis der äquivalenten Rauschbandbreiten sämtlicher Filter in dem Doppler-Signalpfad (ausgehend von dem Demodulator bis hin zu dem "Summierungs- und Ablage"-Filter) ein absoluter Effektivspannungsrauschpegel in Abhängigkeit von dem Verstärkungsgrad des Systems berechnen lassen. Jedes auf eine Analog-Digital-Konvertierung in dem Empfänger zurückzuführende Quantisierungsrauschen kann ebenfalls in einer geeigneten Weise hinzugefügt werden. Weiter sollte es mit Kenntnis der Abtastvolumenposition und der Öffnungsstrategie in dem Spektral-Dopplermodus ohne weiteres möglich sein, das Gesamtsystemrauschen durch Summieren über sämtliche aktive Empfangskanäle (einschließlich der Arrayapodisationseffekte) für eine vorgegebene Abtastvolumenposition zu berechnen. Die mittlere Rauschleistung am Ausgang des Tiefpassfilters kann basierend auf der Bandbreite des Tiefpassfilters berechnet werden. Dem Fachmann sollte klar sein, dass ähnliche Rauschmodelle für Scanner entwickelt werden können, deren Doppler-Signalpfade sich von der grundlegenden Konstruktion in 1 unterscheiden. Während von einem Standpunkt der Verwirklichung her ein Systemrauschmodell deutlich das effizienteste ist, kann ferner zur Durchführung derselben Funktion eine LUT mit mehreren Eingängen verwendet werden. Eine derartige LUT lässt sich entweder anhand von Rauschkalibrierungsmesswerten oder durch eine Simulation des Systemrauschmodells aufstellen. In der ersten Alternative wird das System vorkalibriert, indem verschiedene Kombinationen von Verstärkungsgradvorgabewerten probeweise eingesetzt werden, die sich ergebenden Rauschwerte aufgezeichnet und die betreffenden Verstärkungsgradvorgabewerte und entsprechenden Rauschwerte in einer LUT abgespeichert werden. In der zweiten Alternative werden die Rauschmodellwerte vorausberechnet und in einer LUT gespeichert.
  • Die durch das Systemrauschmodell 30 vorherberechnete mittlere Rauschleistung stellt einen Störschwellenwert bereit, um den Grad der in dem aktuellen FFT-Paket vorhandenen Clutter-Leistung zu messen (Schritt 32). Falls beispielsweise der Schätzwert der Gesamtleistung des gefilterten FFT-Pakets die vorherberechnete Rauschleistung für die betreffende Paketgröße um nicht mehr als z.B. 20 dB übersteigt, ist keine wesentliche Clutter-Störung vorhanden. In diesem Fall wird die Wandfilterauswahllogik 36 automatisch die niedrigste Wandfilter-Eckfrequenz aus der Filterkoeffizienten-LUT 22 auswählen.
  • Falls wesentliche Clutter-Störungsleistung in dem FFT-Paket vorhanden ist (Schritt 32), fährt der Algorithmus mit der Berechnung des Mittelwerts und der Varianz der Clutter-Störungsfrequenz über das FFT-Paket hinweg fort (Schritt 34). Es können verschiedene Verfahren zum Schätzen des Mittelwerts eingesetzt werden, z.B. können das erste Moment des FFT-Spektrums des Pakets berechnet und die Nulldurchgänge gezählt werden. Das bevorzugte Verfahren ist jedoch das Standard-Kasai-Estimatorverfahren, das sowohl genau ist als auch rechnerisch effizient ist – weshalb es gewöhnlich auch in Color-Flow-Bildgebungssystemen verwendet wird. Für eine I/Q-Paketgröße von M sind die Kasai-Estimatoren für den Mittelwert und die Varianz der Frequenz der I/Q-Daten definiert durch: E(f) = (1/2π)tan–1(N/D) var(f) = 2/(2πT)2[1 – |R(T)|/R(0)]mit:
    Figure 00110001
    |R(T)| = (N2 + D2)1/2
    Figure 00110002
    und T gleich der Abtastperiode der I/Q-Daten.
  • Der geschätzte Mittelwert und die Varianz der Clutter-Störungsfrequenz werden anschließend an die Filterauswahllogik 36 ausgegeben, die aus der LUT 22 die am besten geeignete Filter-Eckfrequenz für das aktuelle Clutter-Signal auswählt. Beispielsweise kann die Hochpassfilter-Eckfrequenz basierend auf dem Absolutwert der mittleren Frequenz plus, z.B., dem 3fachen der Standardabweichung berechnet werden. Um eine zu starke Schwankung der Wandfilter-Eckfrequenz von einem FFT-Paket zum nächsten zu vermeiden, kann eine (in der Auswahllogik eingesetzte) beliebige Persistenzfunktion auf die vorgeschriebene Wandfilter-Eckfrequenz angewandt werden, indem beispielsweise der Mittelwert der Eckfrequenz-Schätzwerte über mehrere letzte FFT-Pakete gebildet wird. Nachdem die neue optimale Filter-Eckfrequenz ausgewählt ist, verläuft die übrige Verarbeitung wie im Falle einer herkömmlichen Doppler-Wandfilterung (wie sie in 2 gezeigt ist).
  • Während die Erfindung anhand bevorzugter Ausführungsbeispiele beschrieben wurde, ist es dem Fachmann klar, dass vielfältige Änderungen vorgenommen werden können, und dass Elemente davon durch äquivalente Ausführungen substituiert werden können, ohne dass der durch die beigefügten Patentansprüche definierte Schutzumfang der Erfindung berührt ist. Insbesondere sollte dem Fachmann klar sein, dass das in 3 gezeigte Verfahren als Hardware (z.B. als digitaler Signalverarbeitungsprozessor) und/oder als Software verwirklicht werden kann.

Claims (10)

  1. System zur Ultraschallrückstreuungsbildgebung mit: einem Sender (2, 4) zum Aussenden einer Vielzahl von Ultraschallimpulsen in ein ultraschallstreuendes Volumen, wobei die Vielzahl von Impulsen ein Paket bildet, einem Empfänger (4) zur Aufnahme einer Vielzahl von aufeinander folgenden Proben von aus dem ultraschallstreuenden Probenvolumen rückgestreuten Dopplersignalen, einem Anzeigemonitor (18), der Pixelspalten aufweist, einem Prozessor zur Verarbeitung der Dopplersignale zur Erzeugung von Spektralliniendaten, die die Geschwindigkeit repräsentieren, um bei der Dopplersignalakquisition eine Folge von Spektrallinien fortwährend anzuzeigen, wobei die Daten für jede Spektrallinie eine Anzahl von Frequenzabschnitten für unterschiedliche Frequenzintervalle umfasst, wobei jedes Paket Leistungsspektrumsdaten für ein entsprechendes Frequenzintervall enthält, wobei die Leistungsspektrumsdaten in jedem Paket für jede Spektrallinie als eine Skala in einem entsprechenden Pixel einer entsprechenden Pixelspalte auf dem Anzeigemonitor wiedergegeben werden, wobei der Prozessor aufweist: einen Demodulator (6) zur Demodulation der Dopplersignalproben in Phasen- und Quadratur-Komponenten, einen Integrator (8) zur Integration der In-Phase- und der Quadratur-Komponente über ein spezielles Zeitintervall, erste und zweite Wandfilter (10), um eine Clutter-Störung in der integrierten, In-Phase- und der Quadraturkomponente (I/Q) abzuweisen, die stationären oder sich sehr langsam bewegenden Gewebe entsprechen gekennzeichnet durch: Eckfrequenzauswahlmittel (36) zur Auswahl der Eckfrequenzen des ersten und des zweiten Wandfilters (10) in Abhängigkeit von der Clutter-Störungsleistung in dem integrierten In-Phase- und Quadratur-Komponente (I/Q) und eine Spektralanalyseeinrichtung (12) zur Erzeugung von Spektralleistungsdaten mittels Fast-Fourier-Transformation in den wandgefilterten In-Phase- und Quadratur-Komponenten, wobei das Eckfrequenzauswahlmittel aufweist: ein Tiefpassfilter (26) zur Abweisung von Fließfrequenzkomponenten in der integrierten In-Phase- und Quadratur-Komponente, die über einem Clutter-Störungsfrequenzbereich liegen, Mitteln zur Berechnung der Gesamtleistung (28) des Ausgangssignals des Tiefpassfilters, Mittel, die eine Schätzung (30) der mittleren Systemrauschleistung in dem Ausgangssignal des Tiefpassfilters liefern, Mittel zum Vergleich (33) der Gesamtleistung des Tiefpassfilterausgangssignals mit der mittleren Systemrauschleistung in dem Tiefpassfilterausgangssignal und Mittel (36, 34, 22) zur Koeffizientenselektion zur Auswahl von Wandfilterkoeffizienten, die eine Funktion der Ergebnisse des Vergleichs der Gesamtleistung mit der mittleren Systemrauschleistung sind.
  2. System nach Anspruch 1, bei dem das Filterkoeffizientenauswahlmittel aufweist: Mittel zur Berechnung des Mittelwerts und der Varianz (34) der Clutter-Störungsfrequenz des Tiefpassfilterausgangssignals und einen Speicher zur Speicherung (22) von Wandfilterkoeffizienten, die den verschiedenen Eckfrequenzen entsprechen, und eine Auswahllogik (36) zum Abrufen der Wandfilterkoeffizienten aus dem Speicher in Abhängigkeit von dem berechneten Mittelwert und der berechneten Varianz der Clutter-Frequenz, wenn das Vergleichermittel feststellt, dass die Gesamtleistung die mittlere Systemrauschleistung um wenigstens ein vorbestimmtes Maß übersteigt.
  3. System nach Anspruch 2, bei dem das Mittel zur Berechnung des Mittelwerts und der Varianz der Clutter-Frequenz den Standard-Kasai-Estimator nutzt.
  4. System nach Anspruch 2, bei dem die Auswahllogik Wandfilterkoeffizienten auswählt, die den niedrigsten, in dem Speicher gespeicherten Eckfrequenz entsprechen, wenn das Vergleichmittel feststellt, dass die Gesamtleistung die mittlere Systemrauschleistung nicht um wenigstens ein vorbestimmtes Maß überschreitet.
  5. System nach Anspruch 1, bei dem das Eckfrequenzauswahlmittel ein Mittel zur Auswahl der Filterkoeffizienten gemäß einer Persistenzfunktion aufweist, die auf eine Vielzahl von Paketen angewendet wird.
  6. Verfahren zur Anzeige einer Dopplergeschwindigkeits-Zeithüllkurve mit folgenden Schritten: Aussenden einer Vielzahl von Ultraschallimpulsen in ein ultraschallstreuendes Probenvolumen, wobei die Vielzahl der Impulse ein Paket bildet, Akquisition einer Vielzahl von aufeinander folgenden Proben von aus dem ultraschallstreuenden Probenvolumen rückgestreuten Dopplersignalen und Verarbeiten der Dopplersignale zur Erzeugung von Spektralliniendaten, die die Geschwindigkeit repräsentieren, zur fortwährenden Anzeige einer Folge von Spektrallinien während der Dopplersignalakquisition, wobei die Daten für jede Spektrallinie eine Vielzahl von Frequenzabschnitten für unterschiedliche Frequenz intervalle beinhalten, wobei jeder Abschnitt Spektralleistungsdaten für ein entsprechendes Frequenzintervall enthält, wobei die Spektralleistungsdaten jeder Abschnitt für jede Spektrallinie als eine Skale in einem entsprechenden Pixel einer entsprechenden Pixelspalte auf dem Displaymonitor angezeigt werden, wobei der Verarbeitungsschritt durch die folgenden Schritte gekennzeichnet ist: Demodulation der Dopplersignalproben in In-Phase- und Quadratur-Komponenten, Integration der In-Phase-Komponente und der Quadratur-Komponente über ein spezielles Zeitintervall, Wandfilterung der integrierten In-Phase-Komponente und der Quadratur-Komponente, um Clutter-Störungen, die stationärem oder sich sehr langsam bewegenden Gewebe entsprechen, im Wesentlichen abzuweisen, Auswählen von Eckfrequenzen, die in dem Wandfilterungsschritt verwendet werden, in Abhängigkeit von der Clutter-Störungsleistung und den integrierten In-Phase- und Quadratur-Komponenten und Erzeugung von Spektralleistungsdaten durch Fast-Fourier-Transformation der wandgefilterten In-Phase-Komponente und der Quadratur-Komponente, wobei der Schritt der Auswahl der Eckfrequenz die Schritte beinhaltet, dass: die integrierte In-Phase-Komponente und die Quadratur-Komponente tiefpassgefiltert werden, um Flussfrequenzkomponenten im Wesentlichen auszublenden, die oberhalb eines Clutter-Frequenzbereichs liegen, die Gesamtleistung des Ausgangssignals aus dem Tiefpassfilterungsschritt berechnet wird, eine Schätzung der mittleren Systemrauschleistung in dem Tiefpassfilterausgangssignal durchgeführt wird, die Gesamtleistung des Tiefpassfilterausgangssignals mit der mittleren Systemrauschleistung in dem Tiefpassfilterausgangssignal verglichen wird und Bandfilterkoeffizienten ausgewählt werden, die eine Funktion der Ergebnisse des Vergleichs der Gesamtleistung mit der mittleren Systemrauschleistung sind.
  7. Verfahren nach Anspruch 6, bei dem der Schritt der Auswahl der Wandfilterkoeffizienten die Schritte beinhaltet: dass der Mittelwert und die Varianz der Clutter-Frequenz des Tiefpassfilterausgangssignals berechnet werden und dass die Wandfilterkoeffizienten, die den unterschiedlichen Eckfrequenzen entsprechen, vorgespeichert werden und dass die Wandfilterkoeffizienten aus dem Speicher in Abhängigkeit von dem berechneten Mittelwert und der Varianz der Clutter-Frequenz abgerufen werden, wenn die Gesamtleistung die mittlere Systemrauschleistung um wenigstens ein vorbestimmtes Maß überschreitet.
  8. Verfahren nach Anspruch 7, bei dem der Schritt der Berechnung des Mittelwerts und der Varianz der Clutter-Frequenz den Standard-Kasai-Estimator nutzt.
  9. Verfahren nach Anspruch 7, bei dem die Wandfilterkoeffizienten, die den niedrigsten abgespeicherten Eckfrequenzen entsprechen, ausgewählt werden, wenn die Gesamtleistung die mittlere Systemrauschleistung nicht wenigstens um ein vorbestimmtes Maß überschreitet.
  10. Verfahren nach Anspruch 7, bei dem die Wandfilterkoeffizienten gemäß einer Persistenzfunktion ausgewählt werden, die auf eine Vielzahl von Paketen angewandt wird.
DE60026239T 1999-07-09 2000-07-07 Verfahren und Anordnung zur Spektral- Doppler- Bilderzeugung mit adaptivem Zeitbereichwandfilter Expired - Fee Related DE60026239T2 (de)

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
US349586 1999-07-09
US09/349,586 US6296612B1 (en) 1999-07-09 1999-07-09 Method and apparatus for adaptive wall filtering in spectral Doppler ultrasound imaging

Publications (2)

Publication Number Publication Date
DE60026239D1 DE60026239D1 (de) 2006-04-27
DE60026239T2 true DE60026239T2 (de) 2006-11-16

Family

ID=23373061

Family Applications (1)

Application Number Title Priority Date Filing Date
DE60026239T Expired - Fee Related DE60026239T2 (de) 1999-07-09 2000-07-07 Verfahren und Anordnung zur Spektral- Doppler- Bilderzeugung mit adaptivem Zeitbereichwandfilter

Country Status (5)

Country Link
US (2) US6296612B1 (de)
EP (1) EP1067400B1 (de)
JP (1) JP2001137243A (de)
DE (1) DE60026239T2 (de)
IL (1) IL137109A0 (de)

Families Citing this family (83)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6685645B1 (en) 2001-10-20 2004-02-03 Zonare Medical Systems, Inc. Broad-beam imaging
US6733455B2 (en) * 1999-08-20 2004-05-11 Zonare Medical Systems, Inc. System and method for adaptive clutter filtering in ultrasound color flow imaging
GB0030449D0 (en) * 2000-12-13 2001-01-24 Deltex Guernsey Ltd Improvements in or relating to doppler haemodynamic monitors
US6689064B2 (en) 2001-06-22 2004-02-10 Koninklijke Philips Electronics N.V. Ultrasound clutter filter
US6663566B2 (en) * 2002-02-19 2003-12-16 Ge Medical Systems Global Technology Company, Llc Method and apparatus for automatic control of spectral doppler imaging
US6734820B2 (en) * 2002-05-13 2004-05-11 Honeywell International Inc. Methods and apparatus for conversion of radar return data
US7283461B2 (en) * 2002-08-21 2007-10-16 Alcatel Canada Inc. Detection of denial-of-service attacks using frequency domain analysis
US6789427B2 (en) * 2002-09-16 2004-09-14 General Electric Company Phased array ultrasonic inspection method for industrial applications
EP1579244B1 (de) * 2002-12-02 2014-05-21 Koninklijke Philips N.V. Segmentierungswerkzeug zum identifizeirren von flussregionen ineinem abbildungssystem
US6792808B1 (en) 2003-04-30 2004-09-21 General Electric Company Ultrasonic inspection method
EP1636609A1 (de) * 2003-06-10 2006-03-22 Koninklijke Philips Electronics N.V. Benutzeroberfläche für eine 3d farb-ultraschall bildgebende vorrichtung
JP4245428B2 (ja) * 2003-07-29 2009-03-25 東芝医用システムエンジニアリング株式会社 超音波ドプラ診断装置及び画像データ生成方法
US20050054931A1 (en) * 2003-09-09 2005-03-10 Clark David W. Tracking clutter filter for spectral & audio doppler
WO2005036455A2 (en) * 2003-10-08 2005-04-21 Case Western Reserve University Method for displaying spectral trends in complex signals
US7288068B2 (en) * 2003-12-15 2007-10-30 Siemens Medical Solutions Usa, Inc. Automatic optimization for ultrasound medical imaging
US7749166B2 (en) * 2004-04-26 2010-07-06 General Electric Company System and method for filtering in imaging systems
US7118532B2 (en) * 2004-05-06 2006-10-10 General Electric Company Data dependent color wall filters
US7627386B2 (en) 2004-10-07 2009-12-01 Zonaire Medical Systems, Inc. Ultrasound imaging system parameter optimization via fuzzy logic
JPWO2006088094A1 (ja) * 2005-02-17 2008-07-03 松下電器産業株式会社 超音波ドプラ血流計
US8784318B1 (en) 2005-07-22 2014-07-22 Zonare Medical Systems, Inc. Aberration correction using channel data in ultrasound imaging system
US8002705B1 (en) 2005-07-22 2011-08-23 Zonaire Medical Systems, Inc. Continuous transmit focusing method and apparatus for ultrasound imaging system
US8036856B2 (en) * 2005-12-01 2011-10-11 General Electric Company Method and apparatus for automatically adjusting spectral doppler gain
CN101138507B (zh) * 2006-09-05 2010-05-12 深圳迈瑞生物医疗电子股份有限公司 多普勒血流声音信号的处理方法及装置
US20080146922A1 (en) * 2006-10-24 2008-06-19 Zonare Medical Systems, Inc. Control of user interfaces and displays for portable ultrasound unit and docking station
US8013538B2 (en) 2007-01-26 2011-09-06 Integrated Illumination Systems, Inc. TRI-light
JP5247066B2 (ja) * 2007-06-04 2013-07-24 三菱電機株式会社 目標追尾装置
JP5214920B2 (ja) * 2007-07-24 2013-06-19 株式会社東芝 超音波診断装置及び超音波診断装置の音響出力方法
CN101416885B (zh) * 2007-10-22 2012-10-17 深圳迈瑞生物医疗电子股份有限公司 用于超声成像系统的帧平均电路
US20090112096A1 (en) * 2007-10-29 2009-04-30 Aloka Co., Ltd. Methods and apparatus for ultrasound imaging
US9060669B1 (en) 2007-12-20 2015-06-23 Zonare Medical Systems, Inc. System and method for providing variable ultrasound array processing in a post-storage mode
KR101055581B1 (ko) * 2008-03-20 2011-08-23 삼성메디슨 주식회사 초음파 시스템 및 클러터 신호 제거용 필터 설정 방법
US8255487B2 (en) * 2008-05-16 2012-08-28 Integrated Illumination Systems, Inc. Systems and methods for communicating in a lighting network
KR101100551B1 (ko) 2008-12-17 2011-12-29 삼성메디슨 주식회사 클러터 신호를 필터링하는 초음파 시스템 및 방법
KR101120791B1 (ko) * 2009-11-16 2012-03-22 삼성메디슨 주식회사 클러터 필터링을 적응적으로 수행하는 초음파 시스템 및 방법
KR101029407B1 (ko) 2009-12-02 2011-04-14 (주)메디슨 초음파 컬러 도플러 영상 시스템 및 상기 초음파 컬러 도플러 영상 시스템의 클러터 신호 필터링 방법
CN102121991B (zh) * 2010-01-08 2013-01-23 郑州威科姆科技股份有限公司 基于北斗一代卫星信号接收的干扰抑制方法及装置
KR101239809B1 (ko) * 2010-05-07 2013-03-06 삼성메디슨 주식회사 클러터 필터링을 위한 초음파 영상 장치 및 그 방법
KR101232796B1 (ko) * 2010-07-22 2013-02-13 삼성메디슨 주식회사 클러터 필터링을 위한 초음파 영상 장치 및 그 방법
KR101339683B1 (ko) 2012-01-17 2013-12-10 서강대학교산학협력단 초음파 고조파 영상을 위한 적응적 동적 직각 복조 장치 및 방법
US10475138B2 (en) 2015-09-23 2019-11-12 Causam Energy, Inc. Systems and methods for advanced energy network
US10861112B2 (en) 2012-07-31 2020-12-08 Causam Energy, Inc. Systems and methods for advanced energy settlements, network-based messaging, and applications supporting the same on a blockchain platform
US8849715B2 (en) 2012-10-24 2014-09-30 Causam Energy, Inc. System, method, and apparatus for settlement for participation in an electric power grid
US8983669B2 (en) 2012-07-31 2015-03-17 Causam Energy, Inc. System, method, and data packets for messaging for electric power grid elements over a secure internet protocol network
US9513648B2 (en) 2012-07-31 2016-12-06 Causam Energy, Inc. System, method, and apparatus for electric power grid and network management of grid elements
KR20140056597A (ko) 2012-10-30 2014-05-12 삼성전자주식회사 빔 포밍 계수 결정 방법, 빔 포밍 방법 및 초음파 이미징 장치
US10237770B2 (en) 2013-03-15 2019-03-19 DGS Global Systems, Inc. Systems, methods, and devices having databases and automated reports for electronic spectrum management
US10257728B2 (en) 2013-03-15 2019-04-09 DGS Global Systems, Inc. Systems, methods, and devices for electronic spectrum management
US10257727B2 (en) 2013-03-15 2019-04-09 DGS Global Systems, Inc. Systems methods, and devices having databases and automated reports for electronic spectrum management
US8798548B1 (en) 2013-03-15 2014-08-05 DGS Global Systems, Inc. Systems, methods, and devices having databases for electronic spectrum management
US10122479B2 (en) 2017-01-23 2018-11-06 DGS Global Systems, Inc. Systems, methods, and devices for automatic signal detection with temporal feature extraction within a spectrum
US9537586B2 (en) 2013-03-15 2017-01-03 DGS Global Systems, Inc. Systems, methods, and devices for electronic spectrum management with remote access to data in a virtual computing network
US8750156B1 (en) 2013-03-15 2014-06-10 DGS Global Systems, Inc. Systems, methods, and devices for electronic spectrum management for identifying open space
US10299149B2 (en) 2013-03-15 2019-05-21 DGS Global Systems, Inc. Systems, methods, and devices for electronic spectrum management
US10231206B2 (en) 2013-03-15 2019-03-12 DGS Global Systems, Inc. Systems, methods, and devices for electronic spectrum management for identifying signal-emitting devices
US10271233B2 (en) 2013-03-15 2019-04-23 DGS Global Systems, Inc. Systems, methods, and devices for automatic signal detection with temporal feature extraction within a spectrum
US10219163B2 (en) 2013-03-15 2019-02-26 DGS Global Systems, Inc. Systems, methods, and devices for electronic spectrum management
US10244504B2 (en) 2013-03-15 2019-03-26 DGS Global Systems, Inc. Systems, methods, and devices for geolocation with deployable large scale arrays
US8977212B2 (en) 2013-03-15 2015-03-10 DGS Global Systems, Inc. Systems, methods, and devices for electronic spectrum management with remote access to data in a virtual computing network
US11646918B2 (en) 2013-03-15 2023-05-09 Digital Global Systems, Inc. Systems, methods, and devices for electronic spectrum management for identifying open space
US9078162B2 (en) 2013-03-15 2015-07-07 DGS Global Systems, Inc. Systems, methods, and devices for electronic spectrum management
US10257729B2 (en) 2013-03-15 2019-04-09 DGS Global Systems, Inc. Systems, methods, and devices having databases for electronic spectrum management
US8805292B1 (en) 2013-03-15 2014-08-12 DGS Global Systems, Inc. Systems, methods, and devices for electronic spectrum management for identifying signal-emitting devices
US8787836B1 (en) 2013-03-15 2014-07-22 DGS Global Systems, Inc. Systems, methods, and devices having databases and automated reports for electronic spectrum management
US20140336510A1 (en) * 2013-05-08 2014-11-13 Siemens Medical Solutions Usa, Inc. Enhancement in Diagnostic Ultrasound Spectral Doppler Imaging
JP6413616B2 (ja) * 2013-11-06 2018-10-31 コニカミノルタ株式会社 超音波診断装置、超音波診断装置の制御器及び超音波診断装置の制御方法
WO2016153098A1 (ko) * 2015-03-26 2016-09-29 알피니언메디칼시스템 주식회사 신호의 예측 퀄리티 개선 방법 및 장치
CN105738887B (zh) * 2016-01-29 2018-03-06 西安电子科技大学 基于多普勒通道划分的机载雷达杂波功率谱的优化方法
KR102577752B1 (ko) * 2016-02-02 2023-09-12 삼성메디슨 주식회사 대상체의 속도를 출력하는 방법 및 이를 위한 초음파 진단 장치
EP3323350A3 (de) * 2016-11-16 2018-06-20 Huntleigh Technology Limited Ultraschall-doppler-überwachungssystem
US10459020B2 (en) 2017-01-23 2019-10-29 DGS Global Systems, Inc. Systems, methods, and devices for automatic signal detection based on power distribution by frequency over time within a spectrum
US10700794B2 (en) 2017-01-23 2020-06-30 Digital Global Systems, Inc. Systems, methods, and devices for automatic signal detection based on power distribution by frequency over time within an electromagnetic spectrum
US10498951B2 (en) 2017-01-23 2019-12-03 Digital Global Systems, Inc. Systems, methods, and devices for unmanned vehicle detection
US10529241B2 (en) 2017-01-23 2020-01-07 Digital Global Systems, Inc. Unmanned vehicle recognition and threat management
US11076835B2 (en) * 2018-04-30 2021-08-03 Bk Medical, Aps Model-based control of a dynamic range of an ultrasound image
US10943461B2 (en) 2018-08-24 2021-03-09 Digital Global Systems, Inc. Systems, methods, and devices for automatic signal detection based on power distribution by frequency over time
JP6483910B1 (ja) * 2018-11-28 2019-03-13 メロディ・インターナショナル株式会社 超音波検査装置
CN110353727B (zh) * 2019-06-03 2021-11-23 浙江聚康生物工程有限公司 无线掌上彩色成像设备及成像方法
JP2021058583A (ja) * 2019-10-04 2021-04-15 キヤノンメディカルシステムズ株式会社 超音波診断装置、及び検査方法
CN110907930B (zh) * 2019-11-29 2021-09-03 成都纳雷科技有限公司 一种基于角度估计的车载雷达目标检测估计方法及装置
CN111388010B (zh) * 2020-03-26 2022-06-24 深圳开立生物医疗科技股份有限公司 超声多普勒血流成像方法、装置、设备及可读存储介质
CN111388011B (zh) * 2020-03-26 2023-02-21 深圳开立生物医疗科技股份有限公司 超声多普勒血流成像方法、装置、设备及可读存储介质
US20210393242A1 (en) * 2020-06-19 2021-12-23 EchoNous, Inc. Device and methods for motion artifact suppression in auscultation and ultrasound data
CN112039453A (zh) * 2020-08-04 2020-12-04 中电科仪器仪表有限公司 一种低频信号自动电平控制电路

Family Cites Families (12)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP3093823B2 (ja) * 1991-07-02 2000-10-03 フクダ電子株式会社 超音波ドプラ診断装置
US5269308A (en) * 1991-07-25 1993-12-14 Matsushita Electric Industrial Co., Ltd. Ultrasonic doppler imaging apparatus
US5299174A (en) * 1992-04-10 1994-03-29 Diasonics, Inc. Automatic clutter elimination
US5287753A (en) * 1992-05-02 1994-02-22 Advanced Technology Laboratories, Inc. Continuous display of peak and mean blood flow velocities
US5349524A (en) * 1993-01-08 1994-09-20 General Electric Company Color flow imaging system utilizing a time domain adaptive wall filter
US5443071A (en) * 1994-10-11 1995-08-22 Siemens Medical Systems, Inc. Quantitative color flow
US5524629A (en) * 1994-11-02 1996-06-11 General Electric Company Color flow processor having adaptive wall filter
US5544659A (en) * 1994-12-29 1996-08-13 Siemens Medical Systems, Inc. Ultrasonic doppler imager having a reduced hardware adaptive tissue rejection filter arrangement
US5494037A (en) * 1994-12-29 1996-02-27 Siemens Medical Systems, Inc. Ultrasonic doppler imager having a spatially smoothed control signal for an adaptive tissue rejection filter
US5827969A (en) * 1996-06-12 1998-10-27 Medasonics, Inc. Portable hand held doppler fetal heart rate probe with selective power settings
JP3946288B2 (ja) * 1996-10-01 2007-07-18 東芝医用システムエンジニアリング株式会社 超音波カラードプラ診断装置および超音波カラードプライメージングの信号処理方法
US6146331A (en) * 1998-09-30 2000-11-14 Siemens Medical Systems, Inc. Method for improved clutter suppression for ultrasonic color doppler imaging

Also Published As

Publication number Publication date
US6296612B1 (en) 2001-10-02
JP2001137243A (ja) 2001-05-22
EP1067400A3 (de) 2004-05-26
IL137109A0 (en) 2001-06-14
US6390983B1 (en) 2002-05-21
EP1067400A2 (de) 2001-01-10
EP1067400B1 (de) 2006-03-01
DE60026239D1 (de) 2006-04-27

Similar Documents

Publication Publication Date Title
DE60026239T2 (de) Verfahren und Anordnung zur Spektral- Doppler- Bilderzeugung mit adaptivem Zeitbereichwandfilter
DE69937284T2 (de) Automatische Regelung des Geschwindigkeitsgrads und Pulsfolgefrequenz für Doppler- Ultraschall- Spektrogramme
DE3417660C2 (de) Vorrichtung zur Ultraschalluntersuchung eines Kreislaufsystems an lebenden biologischen Strukturen, sowie Verfahren zur Ultraschallabbildung von Parametern eines Kreislaufsystems
DE69937422T2 (de) Ultraschallabbildung mittels kodierter Anregung beim Senden und selektiver Filterung beim Empfang
DE19520920A1 (de) Verfahren zum Bestimmen des Geschwindigkeit-Zeit-Spektrums einer Blutströmung
DE102005034697B9 (de) Kontrastmittelbildgebung mit einer (Kontrast)mittelspezifischen Ultraschalldetektion
DE10058449B4 (de) Verfahren und Vorrichtung zur Visualisierung einer Bewegung bei einer Ultraschallflussabbildung unter Verwendung einer kontinuierlichen Datenerfassung
DE69530686T2 (de) Gerät zur gewebecharakterisierung mittels ultraschall
DE69735927T2 (de) Diagnostik-Bilderzeugung mittels Ultraschall unterschiedlicher Sende- und Empfangsfrequenz
DE3431001C2 (de)
DE69927040T2 (de) Adaptive unterdrückung von abkling-artefakten in der ivus-bildgebung
DE2406630C2 (de) Verfahren und Vorrichtung zur Durchflußgeschwindigkeitsmessung
DE19819893B4 (de) Verfahren und Einrichtung zum Verbessern der Auflösung und Empfindlichkeit bei der Farbströmungs-Ultraschall-Bildgebung
DE3735121C2 (de)
DE10039346B4 (de) Verfahren und Vorrichtung zur dynamischen Rauschverringerung für ein Doppler-Audioausgangssignal
DE19913198A1 (de) Verfahren und Einrichtung zur verbesserten Flußbilderzeugung in B-Modus-Ultraschall
DE2546856A1 (de) Verfahren und vorrichtung zum messen der frequenz bzw. periodendauer eines signals
DE19756730A1 (de) Verfahren, Einrichtung und Anwendungen zur Verknüpfung von Sende-Wellenfunktionen zur Gewinnung einer synthetichen Wellenform in einem Ultraschall-Bildgebungssystem
DE19912089A1 (de) Verfahren und Einrichtung zur Farbfluß-Bildgebung unter Verwendung von Golay-codierter Anregung beim Senden und Pulskomprimierung beim Empfangen
DE60024162T2 (de) Verfahren zur Farbbilderzeugung und Doppler-Anordnung zur Farbbilderzeugung
EP1043602A1 (de) Verfahren zur Messung der Entfernung mindestens eines Ziels
DE102016116658A1 (de) Sparkleartefakt-erkennung im ultraschall-farbfluss
DE19513226C2 (de) Verfahren und Vorrichtung zur Ultraschallabbildung unter Vermeidung des Kordmusterartefakts bei der Parallelflussoperation
DE19615353C2 (de) Verfahren zum aufwandgünstigen Bestimmen einer Impulsantwort eines hochauflösenden bandbegrenzten Radarkanals
DE69937789T2 (de) Verfahren und vorrichtung zur bestimmung des rhythmus bei auskultationssignalen

Legal Events

Date Code Title Description
8364 No opposition during term of opposition
8339 Ceased/non-payment of the annual fee