DE69917708T2 - Blutgefässabbildung mit adaptiver Mittelwertbildung - Google Patents
Blutgefässabbildung mit adaptiver Mittelwertbildung Download PDFInfo
- Publication number
- DE69917708T2 DE69917708T2 DE69917708T DE69917708T DE69917708T2 DE 69917708 T2 DE69917708 T2 DE 69917708T2 DE 69917708 T DE69917708 T DE 69917708T DE 69917708 T DE69917708 T DE 69917708T DE 69917708 T2 DE69917708 T2 DE 69917708T2
- Authority
- DE
- Germany
- Prior art keywords
- images
- image
- correlation
- cross
- roi1
- 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 - Lifetime
Links
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T5/00—Image enhancement or restoration
- G06T5/50—Image enhancement or restoration by the use of more than one image, e.g. averaging, subtraction
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B6/00—Apparatus for radiation diagnosis, e.g. combined with radiation therapy equipment
- A61B6/46—Apparatus for radiation diagnosis, e.g. combined with radiation therapy equipment with special arrangements for interfacing with the operator or the patient
- A61B6/461—Displaying means of special interest
- A61B6/464—Displaying means of special interest involving a plurality of displays
Description
- Die Erfindung bezieht sich auf Bildbearbeitung und insbesondere auf die Abbildung von Strukturen, die eine durch einen Zyklus definierte Bewegung haben.
- Sich bewegende Strukturen sind schwierig abzubilden, insbesondere wenn eine Bildgebungsebene in den Raum gestellt ist und die Struktur sich in die Bildgebungsebene hinein und aus dieser heraus bewegt. Dies ist besonders schwierig, wenn eine zweite periodische Bewegung hinzugefügt ist. Beispielsweise bewirkt die Bildgebung von Strukturen in einem Patienten, der atmet, eine periodische Bewegung von inneren Strukturen. Dies wird weiter kompliziert durch die Schlagbewegung des Herzens, wenn die Struktur auf oder nahe dem Herzen ist. Der Schlag des Herzens ist eine komplexe Kompressions- und Verdrehbewegung, die der Atmungsbewegung überlagert ist.
- Es ist häufig erwünscht, Bilder von sich bewegenden Strukturen zu gewinnen, wie beispielsweise bei der Koronar-Angiographie, wo Bilder von Gefäßen in einem schlagenden Herz gewonnen werden. Für Bildverbesserungszwecke werden mehrere Bilder Bemittelt, um Bildrauschen zu verringern und die Sichtbarmachung zu verbessern.
- MR Koronar-Angiographie ist üblicherweise dadurch ausgeführt worden, dass entweder schräge Einscheiben-Techniken bei angehaltenem Atem oder durch das Atmen gesteuerte 3D Techniken verwendet werden. Das wiederholte Anhalten des Atmens kann für einige Koronar-Patienten nicht ratsam sein, und Navigationstechniken haben bis heute nicht für ein robustes Verfahren gesorgt, das über einem Bereich unterschiedlicher Atmungsmuster bei einer Vielfalt von Patienten arbeitet. Weiterhin ist bei diesen beiden Lösungen Erfolg oder Versagen häufig nicht deut lich für eine gewisse Zeit nach dem Start der Bildgebung. Der Zweck dieser Studie war, eine robuste Technik ohne Anhalten des Atmens für Koronar-MRI auf der Basis einer adaptiven Mittelwertbildung zu entwickeln.
- Gegenwärtig besteht ein Bedürfnis für ein Bildgebungsverfahren, das klarere Bilder von sich bewegenden Strukturen erzeugt als es bisher möglich war.
- In dem Bildgebungsverfahren der Erfindung werden Realzeitbilder von hoher Qualität von einer gewünschten Struktur mit einer periodischen Bewegung von einem Patienten gewonnen.
- Es werden mehrere Bilder eingefangen und für einen Benutzer bildlich dargestellt. Der Benutzer wählt interaktiv ein Bild, das die gewünschte Struktur als ein Referenzbild R enthält.
- Ein Bereich von Interesse (ROI1) wird auf dem Referenzbild R identifiziert, der die gewünschte Struktur umschließt.
- ROI1 wird dann quer- bzw. kreuzkorreliert mit jedem nachfolgenden eingefangenen Bild, um ein Quer- bzw. Kreuzkorrelationsbild zur Folge zu haben. In dem effizientesten Ausführungsbeispiel wird die Kreuzkorrelation durch Multiplikation in dem Fourier- oder k-Raum-Bereich ausgeführt.
- Ein zweiter Bereich von Interesse wird in der Mitte von dem Kreuzkorrelationsbild definiert. Wenn das Kreuzkorrelationsbild einen Spitzenwert in dem Bereich von Interesse zeigt und mehr als eine vorbestimmte Strecke von dem Rand des interessierenden Bereiches weg ist, wird ein "nutzbares" Bild identifiziert.
- Es wird eine Anzahl nutzbarer Bilder gewonnen und sie können verschoben und gemittelt werden oder nach ihren Versetzungen bzw. Offsets sortiert werden, und nur diejenigen mit ähnlichen Versetzungen werden gemittelt.
- Dies erzeugt gemittelte Bilder mit weniger Verschmieren, da es eine Kompensation für Versetzungen gibt. Weiterhin sorgt die Mittelung für ein Bild höherer Qualität.
- In einem alternativen Ausführungsbeispiel wird für eine Rückführung zu der Datenquelle gesorgt, wenn das gegenwärtige Bild ein "nutzbares" Bild ist. Es wird zusätzliche Information gewonnen, um die ursprünglichen Rohdaten zu verbessern, und es entsteht ein Bild mit höherer Auflösung.
- Es ist eine Aufgabe der vorliegenden Erfindung, ein System bereitzustellen, das für Bilder mit einem hohen Signal/Rausch-Verhältnis von Strukturen sorgt, die eine periodische Bewegung haben.
- Es ist eine weitere Aufgabe der Erfindung, Bilder mit einem hohen Signal/Rausch-Verhältnis von Strukturen von Patienten bereitzustellen, ohne dass es notwendig ist, dass sie während der Bildgebung ihren Atem anhalten.
- Ferner ist es Aufgabe der vorliegenden Erfindung, Bilder mit einem hohen Signal/Rausch-Verhältnis von Koronar-Strukturen von Patienten bereitzustellen, während die Patienten normal atmen.
- Die Merkmale der Erfindung, die für neuartig gehalten werden, sind insbesondere in den beigefügten Ansprüchen angegeben. Die Erfindung selbst jedoch kann sowohl bezüglich ihres Aufbaues und des Arbeitsverfahrens zusammen mit weiteren Aufgaben und Vorteilen davon am besten unter Bezugnahme auf die folgende Beschreibung in Verbindung mit den beigefügten Zeichnungen verstanden werden, in denen:
-
1 –4 ungesteuerte Realzeitbilder ohne Anhalten des Atems von einem Querschnitt von einem Patientenherz sind; -
5 –8 auf entsprechende Weise Kreuzkorrelationen von einem interessierenden Bereich 1 (ROI1) mit jeder der1 –4 sind; -
9 ein Querschnittsbild von den Patienten ist, das aus einer Mittelung der Bilder gemäß den1 –3 gewonnen ist; -
10 und11 zusammen ein vereinfachtes Blockdiagramm von einem Ausführungsbeispiel gemäß der Erfindung sind. - Bei der Bildgebung von sich bewegenden Strukturen treten mehrere Probleme auf. Diese betreffen viele unterschiedliche Disziplinen, aber eine der besten Anwendungen würde in der Koronar-Angiographie sein. Sie wird für diese Anwendung erläutert, aber sie kann für jede Bildgebung von Strukturen verwendet werden, die eine periodische Bewegung durchlaufen, wo die Struktur in eine Bildgebungsebene hineintritt und aus dieser heraustritt.
- Gefäße des Herzens können sich aufgrund des Schlagens des Patientenherzens in die Bildgebungsebene hinein und aus dieser heraus bewegen. Das Atmen des Patienten kann auch bewirken, dass sich Koronargefäße in die Bildgebungsebene hinein und aus diesen heraus bewegen, wenn sich der Brustkorb des Patienten bewegt.
- Wenn man Bilder identifizieren sollte, in denen ein gewünschtes Koronargefäß in einer Bildgebungsebene ist, und diese Bilder verwenden wollte, um ein Durchschnittsbild zu erzeugen oder eine Reihe von Bildern, würde die Bildqualität stark verbessert sein. Wenn man weiter die Versetzung (Offset) des Koronargefäßes in jedem Bild wiedergewinnen wollte, könnte man diese Versetzung verwenden, um die Position des Bildes vor der Mittelung mit anderen in der Sequenz zu korrigieren und somit ein Verschmieren aufgrund der Bewegung des Gefäßes in der Bildgebungsebene vermeiden.
- Der Benutzer arbeitet mit dem System interaktiv zusammen, während er ein Realzeit-Bildgebungsfenster betrachtet, um einen Bereich von Interesse (ROI1) um ein gewünschtes Gefäß herum zu definieren. Dieser ROI1 wird benutzt, um einen Kernel zu erzeugen, der mit jedem neuen Rahmen quer- bzw. kreuzkorreliert wird, wenn er gewonnen wird. Um die Kreuzkorrelation zu beschleunigen, wird sie durch Multiplikation im Fourier (k-Raum)-Bereich ausgeführt.
- Rahmen, die mit dem Referenzkernel eng zusammenpassen, erzeugen eine Kreuzkorrelation mit einem relativ großen zentralen Spitzenwert. Diejenigen Rahmen, in denen die Anatomie innerhalb der Ebene verschoben sind, erzeugen eine Versetzung in dem Spitzenwert, und diejenigen Rahmen, wo sich das Gefäß aus der Ebene heraus bewegt hat, haben ein relativ kleines Signal im Mittelbereich. Es wird ein Spitzenschwellenwert als eine Basis verwendet, um Rahmen zurückzuweisen. Bilder können vor der Mittelung verschoben werden, wenn es keine Rotation oder Deformation gibt, und die Versetzung des Spitzenwertes aus der Mitte des Kreuzkorrelationsbildes gibt die erforderliche Translation an, um das ursprüngliche Bild auszurichten. Bilder, die Spitzenwerte haben, die größer als eine vorbestimmte räumliche Versetzung von der Mitte sind, werden ebenfalls zurückgewiesen.
- In den
1 –4 wurden angiographische Herzbilder gewonnen, die nicht mit entweder den Atmungs- oder Herzzyklen gesteuert sind, während der Patient normal atmete. Jedes Bild stellt eine spezielle Zeitperiode während des Herzzyklus dar. - ROI1 wurde um das Gefäß herum gezogen in einem der Bilder, wo das Gefäß deutlich sichtbar war. In
1 ist ROI1 ein weißes Kästchen, das den größten Teil der Koronararterie umschließt, die die helle gefärbte Linie ist. Diese gleiche Arterie kann in den2 und3 gesehen werden, jedoch hat sich in4 das Herz in eine Position bewegt, wo diese Arterie nicht länger in der Bildgebungsebene ist. - Die Quer- bzw. Kreuzkorrelationen für die
1 ,2 ,3 ,4 sind in den5 ,6 ,7 bzw.8 gezeigt.5 ,6 ,7 zeigen deutliche Korrelationsspitzen, die aus der engen Entsprechung zwischen den Basisbildern (1 ,2 ,3 ) und dem Kernel entstehen (5 ist tatsächlich eine Autokorrelation). Das Koronargefäß in3 ist in der Ebene aufgrund der Atmungsbewegung verschoben, und die Korrelationsspitze in7 ist ähnlich nach oben verschoben. Ein Schwellenwert war als ein Akzeptanz-Kriterium gesetzt, und akzeptierte Rahmen aus dem vollen Datensatz wurden Bemittelt. Die Mittelung von1 und die verschobenen2 und3 hatten9 zur Folge, die klarer als jedes Bild für sich selbst ist. -
10 und11 zeigen zusammen ein vereinfachtes Blockdiagramm von dem Koronar-Angiographiesystem gemäß der Erfindung. Ein Benutzer1 arbeitet interaktiv mit einem Benutzer-Interface13 zusammen, das für eine Eingabe in eine Bildgewinnungsvorrichtung3 sorgt, die möglicherweise ein Magnetresonanz (MR)-Scanner oder eine andere Bildgewinnungsvorrichtung sein kann, die Rohdaten im k-Raum gewinnt. Die Bildgewinnungsvorrichtung3 gewinnt Bilddaten, wie beispielsweise Koronar-Angiographiescheiben, und speichert sie in einem Rohpuffer5 , die unbearbeitete MR Daten und noch in dem Fourier-Zeitbereich oder "k-Raum" sind. Diese Rohdaten stellen mehrere MR Bilder dar. Die Bildgewinnungsvorrichtung und der Rohpuffer arbeiten als eine Datenquelle6 zusammen. - Die Rohbilder werden an eine zweidimensionale schnelle Fourier-Transformations (2DFT)-Vorrichtung
7 geliefert, die diese zurück in den Bildraumbereich transformiert und einen Kreispuffer9 mit Bildern füllt. - Der Benutzer
1 arbeitet interaktiv mit dem Benutzer-Interface13 , um eines der Bilder im Kreispuffer9 zu wählen, das interaktiv auf einer Displayvorrichtung11 als ein Referenzbild R dargestellt wird. - Der Benutzer
1 arbeitet dann interaktiv mit dem Benutzer-Interface13 zusammen, um einen Bereich von Interesse (ROI1) auf dem Referenzbild R zu wählen, der eine gewünschte Struktur umschließt, die bildlich dargestellt werden soll, wie beispielsweise ein Koronargefäß. - ROI1 aus dem Interface
13 wird an eine Maskenvorrichtung15 geliefert, die eine Maske M von Nullen auf allen Flächen außerhalb des ROI1 und Einsen innerhalb des ROI1 erzeugt. - Die Maske M wird an eine Multipliziervorrichtung
17 zusammen mit dem gewählten Referenzbild R aus dem Kreispuffer9 geliefert und multipliziert, um ein maskiertes Bild Im zur Folge zu haben, das im wesentlichen Nullen außerhalb des ROI1 und einen Abschnitt des Referenzbildes R innerhalb des ROI1 hat. - Die Maske M wird an eine inverse Fourier-Transformationsvorrichtung
21 geliefert, um eine k-Raum-Darstellung der Maske M zu liefern. - Diese k-Raum-Darstellung der Maske M hat alle seine Koordinaten in Bezug auf die X-Achse in einer Flip- bzw. Klapp-X-Vorrichtung
25 geklappt, und dieses Ergebnis wird an eine Flip- bzw. Klapp-Y-Vorrichtung29 weitergeleitet, die alle diese Koeffizienten in Bezug auf die Y-Achse klappt, um einen Satz von Gewichtungs-Koeffizienten W zur Folge zu haben, die in einem Gewichtungs-Speicher33 gespeichert werden. - Ein ähnlicher Prozess erfolgt für das maskierte Bild Im mit der inversen FFT Vorrichtung
19 , der Flip- bzw. Klapp-X-Vorrichtung23 bzw. der Flip- bzw. Klapp-Y-Vorrichtung27 . Die Ausgangsgröße der Klapp-Y-Vorrichtung27 ist ein Kernelsatz K, der in dem Kernel-Speicher31 gespeichert wird. - Sowohl der Kernelsatz K als auch der Gewichtungssatz W sind im k-Raum. Der Kernelsatz K und der Gewichtungssatz W werden verwendet, um k-Raum-Bilddaten D aus dem Rohpuffer
5 zu multiplizieren. Eine k-Raum-Bildscheibe D aus dem Rohpuffer5 wird an die Multiplizierer35 und37 mit jedem Gewichtungssatz W bzw. Kernelsatz K geliefert. - Die Ergebnisse der Multiplizierer
35 und37 werden an zweidimensionale Vorwärts-Fourier-Transformations (2DFT)-Vorrichtungen39 bzw.41 geliefert, um Bildraumdarstellungen (A1, ϕ1); A2, ϕ2) zu erzeugen, wobei A die Amplitudeninformation und ϕ die Phaseninformation sind. - Die Größenvorrichtungen
43 und45 empfangen die Phasen und Amplituden von ihren entsprechenden Vorwärts-2DFT-Vorrichtungen und erzeugen Bildgrößen |I1| und |I2|, die die Größen des Bildes an jedem Punkt über den zwei Dimensionen des Bildes darstellen. An jedem Ort wird I2 durch I1 dividiert, um Ic zur Folge zu haben, die die Intensität der Korrelation an diesem Ort ist. Dies hat eine zweidimensionale Korrelationsmap bzw. -abbildung zur Folge. - Die Elemente, die durch die gestrichelte Linie umschlossen sind, führen eine Korrelation im k-Raum aus und werden gemeinsam als eine k-Raum-Korrelationsvorrichtung 30 bezeichnet.
- Ein Finde-Max-Vorrichtung
49 ermittelt, ob ein Spitzenwert, der eine lokale maximale Intensität von Ic über einem zweiten Bereich von Interesse (ROI2) ist, in dem zweidimensionalen Kreuzkorrelations-Bildfeld zentriert ist. ROI2 kann vorbestimmt und gespeichert sein oder er kann vom Benutzer über das Benutzer-Interface13 geliefert werden. - Der Ort der größten Bildintensität Ic auf der Korrelationsmap in ROI2 wird von der Mitte des Bildes subtrahiert, um eine Versetzung zu ermitteln. Die Höhe der maximalen Korrelationsbildintensität Ic ist der maximale Korrelationswert R.
- Eine Bildwählvorrichtung
51 empfängt R in der Versetzung von der Finde-Max-Vorrichtung49 und empfängt auch ein k-Raum-Bild von dem Rohpuffer5 . - Ein Bildwähler
51 hat eine vordefinierte Verknüpfung in sich um zu ermitteln, ob es einen Spitzenwert in der Korrelationsmap gibt, und er ermittelt auch, ob der Spitzenwert in einem vorbestimmten Abstand von den Rändern des Korrelationsbildes ist. Er kann auch ermitteln, ob es keinen Spitzenwert gibt, aber die Korrelationsintensitätswerte Ic ständig in Richtung auf den Rand des ROI zunehmen. Dies zeigt an, dass der ROI1 die Koronargefäße des korrelierten Bildes nicht überlappt. Wenn dies der Fall ist, wird das Bild aussortiert, und es wird ein anderes Bild an die Multiplizierer35 und37 geliefert, wobei der Prozess bis zu diesem Punkt wiederholt wird. - Wenn jedoch die Bildwählvorrichtung
51 ermittelt, dass es einen Spitzenwert Ic gibt und dieser mehr als eine minimale Strecke von dem Rand des ROI entfernt ist, bestimmt er, dass dies ein gutes Bild, ein "nutzbares" Bild ist, und eine Vorwärts-2DFT-Vorrichtung52 transformiert es in den Bildraum und liefert ein Bild an eine Router-Vorrichtung53 zusammen mit der Versetzungs-Information. - In einem Ausführungsbeispiel speichert der Router
53 nur Bilder in einem segmentierten Speicher55 , der mehrere Orte hat, die jeweils die Information halten können, die eine Bildscheibe darstellt. - In einem zweiten Ausführungsbeispiel empfängt der Router
53 mehr Bilder und Versetzungen als die Anzahl der Bins bzw. Fächer im Speicher55 . Ein Sortierer52 in dem Router53 speichert dann Bilder mit einem vorbestimmten Versetzungsbereich in jedem der Bins in dem segmentierten Speicher55 . Jedes zusätzliche Bild, das in einem Bin angeordnet ist, wird mit einem ankommenden Mittelwert von allen vorherigen Bildern Bemittelt, die in diesem Bin gespeichert sind. - Dies hat einen "Bemittelten Film" zur Folge, in dem die Herzbewegung im wesentlichen beseitigt ist, und es wird die Koronararterie gesehen, wie sie sich mit der Atmung auf und ab bewegt.
- In einem noch anderen Ausführungsbeispiel kann der Router
53 durch eine Schiebevorrichtung57 verkörpert sein, die mit der Bildwählvorrichtung51 verbunden ist. Sie empfängt die Bilder und ihre Versetzungen und verschiebt jedes Bild um seine Versetzung, um die Bilder zu normieren, und dann liefert sie die verschobenen Bilder an einen Mittelwertbildner59 . Der Mittelwertbildner59 empfängt auch den Inhalt, der sich bereits in einem laufenden Speicher61 befindet, der somit ein laufender Bildmittler ist. Der Mittler mittelt dann diese Bilder und speichert sie als einen neuen laufenden Mittelwert im laufenden Speicher61 . Jedes dieser Ausführungsbeispiele wird unterschiedliche Anwendungen und Resultate haben. - Dieses adaptiv gemittelte Bild zeigt eine markante Verbesserung im Signal/Rausch-Verhältnis (SNR) relativ zu den Realzeit-Bildern. Eine adaptive Mittlung ist für die Koronar-Bildgebung brauchbar.
- In einem noch anderen Ausführungsbeispiel kann die Bildgewinnungsvorrichtung
3 mit der Bildwählvorrichtung51 verbunden sein und auf diese ansprechen. Bilder können in einer normalen oder abgekürzten Weise gewonnen werden, um die Geschwindigkeit der Gewinnung während der normalen Bildgebung zu vergrößern. Wenn jedoch die Bildwählvorrichtung51 ermittelt, dass sich die gewünschte Struktur in dem laufenden Bild D befindet, wie es durch eine hohe Korrelationsspitze in dem Korrelationsbild angezeigt wird, wird ein Signal an die Bildgewinnungsvorrichtung3 gesendet, um zusätzliche Bildinformation zu gewinnen. Dies kann bedeuten, dass die Bildgewinnungsvorrichtung3 kleinere räumliche Frequenzen in dem k-Raum während schneller Gewinnungen abtastet, aber höhere räumliche Frequenzen beim Empfangen des Signals aus der Bildwählvorrichtung51 abtastet. Die zusätzliche Information kann direkt zur Bildwählvorrichtung51 geleitet werden oder zunächst durch den Rohpuffer hindurchtreten und soll die bereits gewonnenen Rohbilddaten D verbessern. Dies gestattet, dass Bilder mit höherer Auflösung gewonnen werden, wenn die gewünschte Struktur sichtbar ist, während zu anderen Zeiten "schnelle" Bilder mit niedrigerer Auflösung gewonnen werden. - Da es einen ersten Schirm
11 und einen zweiten Schirm71 gibt, können wenigstens zwei unterschiedliche Darstellungen gleichzeitig für den Benutzer1 angezeigt werden. Der erste Schirm11 kann ein Realzeitbild des schlagenden Herzens des Patienten zeigen oder kann abgetrennt sein und irgendwelche oder alle Bilder spielen, die in dem Kreispuffer9 gespeichert sind. - Der zweite Schirm kann den "gemittelten Film" spielen, während er diesen mit der Realzeitdarstellung der ersten Schirmvorrichtung
11 vergleicht, oder er kann ein gemitteltes Bild des laufenden Speichers61 haben, der gemäß der Benutzer-Eingabe aktualisiert ist. Es ist möglich, dass die Aktualisierung kontinuierlich fortschreitet. - Der zweite Schirm
71 kann auch alle oder einige der Bilder im Speicher55 zeigen, die kontinuierlich aktualisiert oder statisch sein können. - Resultate
- Das Koronarbild gemäß den
1 –4 wurde unter Verwendung einer spiralförmigen Fluoreszenzpulssequenz auf einem Hochleistungs-Gradientensystem (Signa Horizon EchoSpeed, GE Medical Systems) gewonnen. Der Scanner wurde von einem UltraSparc2 gesteuert, und Bilder wurden zurückübertragen, rekonstruiert und bildlich dargestellt, um eine Realzeitansicht von dem schlagenden Herz zu erzeugen. Die 6-Sprung-Spiralpulssequenz erzeugte eine Matrixgröße von 128×128 mit 5 vollständigen neuen Bildern pro Sekunde und 9 Rekonstruktionen pro Sekunde. Es wurden Schirm-basierte Werkzeuge benutzt, um interaktiv optimale Koronar-Abtastebenen innerhalb 10–20 Sekunden vor der adaptiven Koronar-Bildgebung anzuordnen. - Realzeit-Spiralfluoreszenz mit adaptiver Mittelung zeigt eine Möglichkeit als eine robuste Technik zur Koronar-MRI. Sie erfordert kein Anhalten des Atems, keine separate Navigationssequenz und keine EKG Steuerung. Darüber hinaus gestattet sie eine kontinuierliche Sichtbarmachung des Gefäßes von Interesse vom Start der Bildgebung. Eine interaktive Lokalisierung von Koronar-Abtastebenen kann auf einfache Weise verschachtelt werden mit adaptiver Mittelwertbildung während einer Realzeit-Bildgebung dieser Abtastebenen.
Claims (10)
- Verfahren zum Gewinnen von Bildern hoher Qualität von einer gewünschten Struktur mit einer periodischen Bewegung in einer Reihe von Bildern, enthaltend die Schritte: a) Anzeigen der Reihe von Bildern zum Akzeptieren einer vom Benutzer gelieferten Eingabe, die ein gewähltes Bild identifiziert, das die gewünschte Struktur als ein Referenzbild R enthält, b) Identifizieren eines interessierenden Bereiches ROI1 im Referenzbild R, das die gewünschte Struktur einschließt, c) Kreuz- bzw. Quer-Korrelieren ROI1 mit jedem der anderen Bilder aus der Reihe von Bildern, um eine Reihe von Kreuz- bzw. Quer-Korrelationsbildern zur Folge zu haben, d) Definieren eines zweiten interessierenden Bereiches ROI2, der in jedem Kreuz- bzw. Quer-Korrelationsbild zentriert ist, e) Identifitzieren als "nutzbare" Bilder, die einen Spitzenwert in ROI2 von ihrem Kreuz- bzw. Quer-Korrelationsbild mehr als eine vorbestimmte Strecke entfernt von einem Rand von ROI2, f) Verarbeiten nur der "nutzbaren" Bilder, um eine verarbeitetes Bild höherer Qualität zur Folge zu haben.
- Verfahren zum Gewinnen von Bildern hoher Qualität nach Anspruch 1, wobei der Schritt des Ausführens einer Kreuz- bzw. Quer-Korrelation die Schritte enthält: a) Ausführen einer Fourier-Transformation an dem Referenzbild R in dem Bereich ROI1 in ein k-Raum-Bild, b) Klappen des k-Raum-Bildes in bezug auf seine X-Koordinaten, c) Klappen des k-Raum-Bildes in bezug auf seine Y-Koordinaten, um einen Kernel-Satz K zur Folge zu haben, d) Erzeugen einer Maske M mit Nullen ausserhalb des Bereiches ROI1 und Einsen innerhalb von ROI1, e) inverses Fourier-Tranformieren der Maske M in ein k-Raum-Bild, Klappen des k-Raum-Bildes der Maske M in bezug auf seine X-Koordinaten, f) Klappen des k-Raum-Bildes der Maske M in bezug auf seine Y-Koordinaten, um einen Gewichtungssatz W zur Folge zu haben, g) Multiplizieren des Kernel-Satzes K und des Gewichtungssatzes W jeweils mit einem Rohdatenbild D, um das Produkt K bzw. Produkt W zur Folge zu haben, h) Ausführen einer 2D Vorwärts-Fourier-Transformation am Produkt K und Produkt W, um eine Amplituden- und Phaseninformation (A2, ϕ2) bzw. (A1, ϕ1) zur Folge zu haben, i) Ermitteln der Grössen |I2|, |I1| aus der Information (A2, ϕ2) bzw. (A1, ϕ1), und j) Ermitteln eines Quer-Korrelationsbildes Ic gemäss:
- Verfahren zum Gewinnen von Bildern hoher Qualität nach Anspruch 1, wobei die Quer-Korrelation durch Multiplikation im k-Raum ausgeführt wird, um eine schnelle Quer-Korrelation zur Folge zu haben.
- Verfahren zum Gewinnen von Bildern hoher Qualität nach Anspruch 1, wobei ferner die Schritte vorgesehen sind: a) Berechnen von Versetzungen von der Lage von Quer-Korrelations-Spitzenwerten nutzbarer Bilder von einem Mittelpunkt ihrer Quer-Korrelations-Bilder, b) Verschieben der nutzbaren Bilder um den Versetzungsbetrag, um die Versetzung zu korrigieren, und c) Mitteln der verschobenen Bilder, um bearbeitete Bilder mit verminderten Verzerrungen zur Folge zu haben.
- Verfahren zum Gewinnen von Bildern hoher Qualität nach Anspruch 1, wobei ferner die Schritte vorgesehen sind: a) Berechnen von Versetzungen von der Lage von Quer-Korrelations-Spitzenwerten von nutzbaren Bildern von einem Mittelpunkt ihrer Quer-Korrelations-Bilder, b) Sortieren von Bildern gemäss ihren Versetzungen, c) Rückspielen der Bilder in ihrer sortierten Reihenfolge zum Liefern eines "gemittelten Films", der sich gemäss der Versetzung ändert.
- Verfahren zum Gewinnen von Bildern hoher Qualität nach Anspruch 3, wobei: a) die Bilder Koronar-Angiographie-Bilder von einem Patienten sind, der während der Gewinnung atmet, b) die abgebildete Struktur ein Koronar-Gefäss ist, c) die Versetzungen aufgrund von Atmungsbewegung bestehen und d) der "gemittelte Film" Herzbewegung im wesentlichen beseitigt.
- System zum bildlichen Darstellen von Bildern hoher Qualität, die von einer Reihe von Bildern von einer gewünschten Struktur mit einer periodischen Bewegung abgeleitet sind, enthaltend: a) eine Datenquelle zum Gewinnen einer Anzahl von Rohbildern im Fourier-k-Raum und zum Liefern dieser Bilder an verbundene Elemente, b) eine mit der Datenquelle verbundene interaktive Bildwählvorrichtung zum Wählen und bildlichen Darstellen Benutzer definierter Bilder und die zum Definieren eines ROI1 auf einem Benutzer-selektierten Referenzbild R arbeitet, c) eine Maskenvorrichtung, die mit der interaktiven Bildwählvorrichtung verbunden ist und zum Empfangen von ROI1 und zum Erzeugen einer Maske M arbeitet, die alle Nullen ausserhalb von ROI1 und alle Einsen innerhalb von ROI1 enthält, d) eine k-Raum Korrelationsvorrichtung, die mit der interaktiven Bildwählvorrichtung, der Maskenvorrichtung und den Datenquellen verbunden ist zum Erzeugen eines Kreuzkorrelationsbildes Ic über einem zweiten interessierenden Bereich ROI2 aus dem Referenzbild R, der Maske M und dem Rohdatenbild D, e) eine Finde-Max-Vorrichtung, die mit der k-Raum Querkorrelationsvorrichtung verbunden ist und die Funktion hat zu berechnen, ob ein Spitzenwert existiert, und eine Spitzenhöhe und eine Versetzung des Spitzenwertes von der Mitte des Querkorrelationsbildes Ic für existierende Spitzenwerte, f) eine Bildwählvorrichtung, die mit der Finde-Max-Vorrichtung verbunden ist und aus einem Rohpuffer, der die Funktion hat, die Versetzung und Spitzenwerthöhe, Rohdatenbild D zu empfangen und zu ermitteln, ob die Spitzenwerthöhe grösser als ein vorbestimmter Schwellenwert ist, und ob die Versetzung kleiner als ein vorbestimmter Schwellenwert ist, der ein "nutzbares" Bild definiert, und die Funktion hat, eine 2D Fourier-Tranformation an den "nutzbaren" Bildern auszuführen, g) einen Bildspeicher, der eine Anzahl von Bildern speichern und laden kann, h) einen Router, der mit der Bildwählvorrichtung verbunden ist und die Funktion hat, die 2DFT Darstellungen der nutzbaren Bilder zum empfangen und diese Bilder weiter zu verarbeiten, dann diese bearbeiteten Bilder im Bildspeicher zu speichern und i) eine Anzeigevorrichtung, die mit dem Bildspeicher verbunden ist, zum bildlichen Darstellen der bearbeiteten Bilder, die in dem Bildspeicher gespeichert sind.
- System zum bildlichen Darstellen von Bildern hoher Qualität nach Anspruch 7, wobei der Router einen Sortierer aufweist, der Bilder bearbeitet, indem er nutzbare Bilder sortiert, die gemäss ihrer Versetzung an ihn geliefert worden sind.
- System zum bildlichen Darstellen von Bildern hoher Qualität nach Anspruch 7, wobei der Router enthält: a) eine Schiebevorrichtung zum Empfangen von Bild und Versetzung, die jedes Bild an seiner Versetzung verschiebt, b) einen laufenden Speicher, der Bilddaten speichern kann, die einen laufenden Bildmittelwert enthalten, und c) einen Mittelwertbildner, der mit dem laufenden Speicher und der Schiebevorrichtung verbunden ist, zum Empfangen und Bilden des Mittelwertes des laufenden Mittelwertbildes und des verschobenen Bildes aus der Schiebevorrichtung und zum Speichern des Ergebnisses in dem laufenden Speicher.
- System zum bildlichen Darstellen von Bildern hoher Qualität nach Anspruch 7, wobei: a) eine Bildgewinnungsvorrichtung in wenigstens zwei Modi arbeitet, wobei der erste ein "schnelle" Gewinnung von abgekürzten Bildern ist, und nach Empfang eines Signals in einem zweiten Modus arbeitet, der ein Gewinnungsmodus hoher Auflösung ist, b) die Bildwählvorrichtung mit der Bildgewinnungsvorrichtung verbunden ist und ein Signal an die Bildgewinnungsvorrichtung sendet, wenn ein gegenwärtiges "nutzbares" Bild erfasst wird, wodurch sie veranlasst wird, zusätzliche Information zu gewinnen, die die ursprünglichen Rohdaten D ergänzen.
Applications Claiming Priority (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
US62417 | 1998-04-17 | ||
US09/062,417 US6088488A (en) | 1998-04-17 | 1998-04-17 | Vascular imaging with adaptive averaging |
Publications (2)
Publication Number | Publication Date |
---|---|
DE69917708D1 DE69917708D1 (de) | 2004-07-08 |
DE69917708T2 true DE69917708T2 (de) | 2005-07-21 |
Family
ID=22042356
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
DE69917708T Expired - Lifetime DE69917708T2 (de) | 1998-04-17 | 1999-03-31 | Blutgefässabbildung mit adaptiver Mittelwertbildung |
Country Status (4)
Country | Link |
---|---|
US (1) | US6088488A (de) |
EP (1) | EP0952547B1 (de) |
JP (1) | JP4316725B2 (de) |
DE (1) | DE69917708T2 (de) |
Families Citing this family (43)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US8788020B2 (en) * | 1998-10-23 | 2014-07-22 | Varian Medical Systems, Inc. | Method and system for radiation application |
US6937696B1 (en) | 1998-10-23 | 2005-08-30 | Varian Medical Systems Technologies, Inc. | Method and system for predictive physiological gating |
AU3528600A (en) * | 1999-10-21 | 2001-04-30 | Navlet.Com, Inc. | Context-sensitive switching in a computer network environment |
WO2001067960A2 (en) * | 2000-03-16 | 2001-09-20 | Cedara Software Corp. | System and method for processing an image sequence |
US7209777B2 (en) | 2000-11-30 | 2007-04-24 | General Electric Company | Method and apparatus for automated tracking of non-linear vessel movement using MR imaging |
US6980846B2 (en) | 2001-03-30 | 2005-12-27 | General Electric Company | Robust coronary MR angiography without respiratory navigation |
GB0115714D0 (en) * | 2001-06-27 | 2001-08-22 | Imperial College | Structure determination of macromolecules |
US7127096B2 (en) * | 2001-11-20 | 2006-10-24 | Accuimage Diagnostics Corp. | Method and software for improving coronary calcium scoring consistency |
CA2421352A1 (en) * | 2002-03-08 | 2003-09-08 | Giovanni Battista Mancini (A.K.A. Mancini, G.B. John) | Vessel evaluation methods, apparatus, computer-readable media and signals |
DE10220295A1 (de) * | 2002-05-07 | 2003-11-20 | Philips Intellectual Property | Verfahren zur Verbesserung der Bildqualität |
DE10221643A1 (de) | 2002-05-15 | 2003-12-04 | Siemens Ag | Auswertung zeitlich veränderlicher diagnostischer Vorgänge |
US7149331B1 (en) | 2002-09-03 | 2006-12-12 | Cedara Software Corp. | Methods and software for improving thresholding of coronary calcium scoring |
US6835177B2 (en) * | 2002-11-06 | 2004-12-28 | Sonosite, Inc. | Ultrasonic blood vessel measurement apparatus and method |
US7620221B2 (en) * | 2003-01-29 | 2009-11-17 | Koninklijke Philips Electronics N.V. | System and method for enhancing an object of interest in noisy medical images |
EP1685535B1 (de) | 2003-08-21 | 2014-04-30 | Philips Intellectual Property & Standards GmbH | Einrichtung und verfahren zum kombinieren zweier bilder |
EP1665130A4 (de) | 2003-09-25 | 2009-11-18 | Paieon Inc | System und verfahren zur dreidimensionalen rekonstruktion eines röhrenförmigen organs |
US20050125843A1 (en) * | 2003-11-05 | 2005-06-09 | Okezie Charles E. | Television viewer/studio interactive commentary |
DE102004061507B4 (de) * | 2004-12-21 | 2007-04-12 | Siemens Ag | Verfahren zur Korrektur von Inhomogenitäten in einem Bild sowie bildgebende Vorrichtung dazu |
DE102006026542B3 (de) * | 2006-06-07 | 2008-01-03 | Siemens Ag | Verfahren zur Bilderzeugung mittels magnetischer Kernresonanz |
WO2009153794A1 (en) | 2008-06-19 | 2009-12-23 | Sync-Rx, Ltd. | Stepwise advancement of a medical tool |
US9629571B2 (en) | 2007-03-08 | 2017-04-25 | Sync-Rx, Ltd. | Co-use of endoluminal data and extraluminal imaging |
US9375164B2 (en) | 2007-03-08 | 2016-06-28 | Sync-Rx, Ltd. | Co-use of endoluminal data and extraluminal imaging |
US8542900B2 (en) | 2007-03-08 | 2013-09-24 | Sync-Rx Ltd. | Automatic reduction of interfering elements from an image stream of a moving organ |
US11197651B2 (en) | 2007-03-08 | 2021-12-14 | Sync-Rx, Ltd. | Identification and presentation of device-to-vessel relative motion |
US10716528B2 (en) | 2007-03-08 | 2020-07-21 | Sync-Rx, Ltd. | Automatic display of previously-acquired endoluminal images |
EP2129284A4 (de) | 2007-03-08 | 2012-11-28 | Sync Rx Ltd | Bildgebung und werkzeuge zur verwendung mit beweglichen organen |
US9968256B2 (en) | 2007-03-08 | 2018-05-15 | Sync-Rx Ltd. | Automatic identification of a tool |
US11064964B2 (en) | 2007-03-08 | 2021-07-20 | Sync-Rx, Ltd | Determining a characteristic of a lumen by measuring velocity of a contrast agent |
CN102007422B (zh) | 2008-04-18 | 2015-05-27 | 皇家飞利浦电子股份有限公司 | 磁共振成像中的回顾运动补偿 |
US10667727B2 (en) | 2008-09-05 | 2020-06-02 | Varian Medical Systems, Inc. | Systems and methods for determining a state of a patient |
US9974509B2 (en) | 2008-11-18 | 2018-05-22 | Sync-Rx Ltd. | Image super enhancement |
US9101286B2 (en) | 2008-11-18 | 2015-08-11 | Sync-Rx, Ltd. | Apparatus and methods for determining a dimension of a portion of a stack of endoluminal data points |
US8855744B2 (en) | 2008-11-18 | 2014-10-07 | Sync-Rx, Ltd. | Displaying a device within an endoluminal image stack |
US9095313B2 (en) | 2008-11-18 | 2015-08-04 | Sync-Rx, Ltd. | Accounting for non-uniform longitudinal motion during movement of an endoluminal imaging probe |
US10362962B2 (en) | 2008-11-18 | 2019-07-30 | Synx-Rx, Ltd. | Accounting for skipped imaging locations during movement of an endoluminal imaging probe |
US9144394B2 (en) | 2008-11-18 | 2015-09-29 | Sync-Rx, Ltd. | Apparatus and methods for determining a plurality of local calibration factors for an image |
US11064903B2 (en) | 2008-11-18 | 2021-07-20 | Sync-Rx, Ltd | Apparatus and methods for mapping a sequence of images to a roadmap image |
EP2723231A4 (de) | 2011-06-23 | 2015-02-25 | Sync Rx Ltd | Luminale hintergrundreinigung |
CA2875346A1 (en) | 2012-06-26 | 2014-01-03 | Sync-Rx, Ltd. | Flow-related image processing in luminal organs |
US9594144B2 (en) | 2014-04-23 | 2017-03-14 | General Electric Company | Low-noise magnetic resonance imaging using low harmonic pulse sequences |
GB2542738B (en) * | 2014-07-18 | 2020-07-29 | Hewlett Packard Development Co | Frequency domain range determination for a periodic or quasi-periodic target |
US10557923B2 (en) * | 2015-02-25 | 2020-02-11 | The Government Of The United States Of America, As Represented By The Secretary Of The Navy | Real-time processing and adaptable illumination lidar camera using a spatial light modulator |
FR3036189B1 (fr) * | 2015-05-13 | 2018-07-27 | Universite Paris-Sud | Procede et dispositif d'imagerie par resonance magnetique |
Family Cites Families (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5020113A (en) * | 1988-08-08 | 1991-05-28 | Hughes Aircraft Company | Maskable bilevel correlators |
US5054045A (en) * | 1990-11-14 | 1991-10-01 | Cedars-Sinai Medical Center | Coronary tracking display |
US5293574A (en) * | 1992-10-23 | 1994-03-08 | General Electric Company | Digital x-ray imaging system with automatic tracking |
US5550937A (en) * | 1992-11-23 | 1996-08-27 | Harris Corporation | Mechanism for registering digital images obtained from multiple sensors having diverse image collection geometries |
DE69331982T2 (de) * | 1992-11-27 | 2003-01-23 | Fuji Photo Film Co Ltd | Verfahren zur Lageanpassung von Röntgenbildern |
US5361307A (en) * | 1993-03-25 | 1994-11-01 | General Electric Company | Correlation methods of identifying defects in imaging devices |
US5431161A (en) * | 1993-04-15 | 1995-07-11 | Adac Laboratories | Method and apparatus for information acquistion, processing, and display within a medical camera system |
JP3410843B2 (ja) * | 1994-12-27 | 2003-05-26 | 株式会社東芝 | 超音波診断装置 |
US5577181A (en) * | 1995-06-07 | 1996-11-19 | E-Systems, Inc. | Method for autonomous determination of tie points in imagery |
US5787889A (en) * | 1996-12-18 | 1998-08-04 | University Of Washington | Ultrasound imaging with real time 3D image reconstruction and visualization |
-
1998
- 1998-04-17 US US09/062,417 patent/US6088488A/en not_active Expired - Lifetime
-
1999
- 1999-03-31 DE DE69917708T patent/DE69917708T2/de not_active Expired - Lifetime
- 1999-03-31 EP EP99302574A patent/EP0952547B1/de not_active Expired - Lifetime
- 1999-04-15 JP JP10733299A patent/JP4316725B2/ja not_active Expired - Fee Related
Also Published As
Publication number | Publication date |
---|---|
EP0952547A3 (de) | 1999-12-22 |
JP4316725B2 (ja) | 2009-08-19 |
EP0952547A2 (de) | 1999-10-27 |
EP0952547B1 (de) | 2004-06-02 |
JP2000048185A (ja) | 2000-02-18 |
US6088488A (en) | 2000-07-11 |
DE69917708D1 (de) | 2004-07-08 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
DE69917708T2 (de) | Blutgefässabbildung mit adaptiver Mittelwertbildung | |
DE3426933C2 (de) | Verfahren und Anordnung zum selbsttätigen Korrigieren von Artefakten aufgrund von Deckungsfehlern | |
DE69813087T2 (de) | Verarbeitung von mittels Ultraschall erzeugten und verbesserten intravaskulären Abbildungen und Signalen | |
DE10312018B4 (de) | System und Verfahren zur Ultraschall-Bildverarbeitung | |
DE60028877T2 (de) | Bildverarbeitungsverfahren, -system und -gerät zur bildung eines gesamten bildes in einer lang gestreckten szene | |
DE19851597B4 (de) | Verfahren und Gerät zur Erzeugung verbesserter Bilder von Kranzarterien | |
DE602004003845T2 (de) | Bildverarbeitungsvorrichtung zur Reduktion von Pixelrauschen | |
DE60131337T2 (de) | Diagnostische Bilderzeugung | |
DE19746938A1 (de) | Datenveranschaulichungsverbesserung durch Entfernung von dominierenden Strukturen | |
DE60014338T2 (de) | Verfahren und gerät zur bildverarbeitung | |
DE19613342A1 (de) | Automatisches Bildauswertungsverfahren | |
DE102006051919B4 (de) | Verfahren zum Bereitstellen eines 3D-Röntgenbilddatensatzes des Herzens eines Patienten | |
EP0938063B1 (de) | Verfahren zur zweidimensionalen Abbildung von Strukturen für die medizinische Diagnostik | |
DE3735519C2 (de) | Binäre Rauminterpolation | |
DE19635017A1 (de) | Verfahren und Vorrichtung zur Erhöhung der Bildschärfe | |
DE102006023093A1 (de) | Verfahren und Vorrichtung zur Korrektur eines Bildes | |
EP1302899A2 (de) | Vorrichtung und Verfahren zur Verarbeitung von Digitalbildern | |
DE69629149T2 (de) | Digitales angiographisches Gerät | |
DE102010043975B4 (de) | Verfahren zur Reduktion der verwendeten Strahlendosis im Rahmen einer bildgebenden Röntgenuntersuchung und Computersystem | |
DE10356174A1 (de) | Verfahren und Einrichtung zur Tomosynthese-Bildverbesserung unter Verwendung von Querfilterung | |
DE102007061934A1 (de) | Verfahren zur Qualitätssteigerung von computertomographischen Aufnahmeserien durch Projektionsdatenverarbeitung und CT-System mit Recheneinheit | |
DE69632775T2 (de) | Geschwindigkeitsanpassbare gefilterte angiographie | |
DE10238322A1 (de) | Retrospektive bzw. fenstergesteuerte Filterung von Bildern zur Adaption von Schärfe und Rauschen in der Computer-Tomographie | |
DE69839400T2 (de) | Röntgenstrahl-bildverarbeitung | |
DE102008048045A1 (de) | Verfahren zur Erzeugung von computertomographischen Bilddatensätzen eines Patienten in der Herz-CT bei einer Perfusionskontrolle unter Kontrastmittelapplikation |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
8364 | No opposition during term of opposition | ||
8328 | Change in the person/name/address of the agent |
Representative=s name: ROEGER UND KOLLEGEN, 73728 ESSLINGEN |