DE102008049709A1 - Verfahren zur selektiven Darstellung einer Bewegung der Lunge, Computerprogramm, Bildverarbeitungseinheit und Magnetresonanzgerät - Google Patents
Verfahren zur selektiven Darstellung einer Bewegung der Lunge, Computerprogramm, Bildverarbeitungseinheit und Magnetresonanzgerät Download PDFInfo
- Publication number
- DE102008049709A1 DE102008049709A1 DE102008049709A DE102008049709A DE102008049709A1 DE 102008049709 A1 DE102008049709 A1 DE 102008049709A1 DE 102008049709 A DE102008049709 A DE 102008049709A DE 102008049709 A DE102008049709 A DE 102008049709A DE 102008049709 A1 DE102008049709 A1 DE 102008049709A1
- Authority
- DE
- Germany
- Prior art keywords
- images
- magnetic resonance
- movement
- frequency
- waveform
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
Classifications
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/05—Detecting, measuring or recording for diagnosis by means of electric currents or magnetic fields; Measuring using microwaves or radio waves
- A61B5/055—Detecting, measuring or recording for diagnosis by means of electric currents or magnetic fields; Measuring using microwaves or radio waves involving electronic [EMR] or nuclear [NMR] magnetic resonance, e.g. magnetic resonance imaging
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/72—Signal processing specially adapted for physiological signals or for diagnostic purposes
- A61B5/7235—Details of waveform analysis
- A61B5/7253—Details of waveform analysis characterised by using transforms
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
- G01R33/00—Arrangements or instruments for measuring magnetic variables
- G01R33/20—Arrangements or instruments for measuring magnetic variables involving magnetic resonance
- G01R33/44—Arrangements or instruments for measuring magnetic variables involving magnetic resonance using nuclear magnetic resonance [NMR]
- G01R33/48—NMR imaging systems
- G01R33/54—Signal processing systems, e.g. using pulse sequences ; Generation or control of pulse sequences; Operator console
- G01R33/56—Image enhancement or correction, e.g. subtraction or averaging techniques, e.g. improvement of signal-to-noise ratio and resolution
- G01R33/561—Image enhancement or correction, e.g. subtraction or averaging techniques, e.g. improvement of signal-to-noise ratio and resolution by reduction of the scanning time, i.e. fast acquiring systems, e.g. using echo-planar pulse sequences
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/72—Signal processing specially adapted for physiological signals or for diagnostic purposes
- A61B5/7203—Signal processing specially adapted for physiological signals or for diagnostic purposes for noise prevention, reduction or removal
- A61B5/7207—Signal processing specially adapted for physiological signals or for diagnostic purposes for noise prevention, reduction or removal of noise induced by motion artifacts
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/72—Signal processing specially adapted for physiological signals or for diagnostic purposes
- A61B5/7235—Details of waveform analysis
- A61B5/7253—Details of waveform analysis characterised by using transforms
- A61B5/7257—Details of waveform analysis characterised by using transforms using Fourier transforms
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
- G01R33/00—Arrangements or instruments for measuring magnetic variables
- G01R33/20—Arrangements or instruments for measuring magnetic variables involving magnetic resonance
- G01R33/44—Arrangements or instruments for measuring magnetic variables involving magnetic resonance using nuclear magnetic resonance [NMR]
- G01R33/48—NMR imaging systems
- G01R33/54—Signal processing systems, e.g. using pulse sequences ; Generation or control of pulse sequences; Operator console
- G01R33/56—Image enhancement or correction, e.g. subtraction or averaging techniques, e.g. improvement of signal-to-noise ratio and resolution
- G01R33/5608—Data processing and visualization specially adapted for MR, e.g. for feature analysis and pattern recognition on the basis of measured MR data, segmentation of measured MR data, edge contour detection on the basis of measured MR data, for enhancing measured MR data in terms of signal-to-noise ratio by means of noise filtering or apodization, for enhancing measured MR data in terms of resolution by means for deblurring, windowing, zero filling, or generation of gray-scaled images, colour-coded images or images displaying vectors instead of pixels
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
- G01R33/00—Arrangements or instruments for measuring magnetic variables
- G01R33/20—Arrangements or instruments for measuring magnetic variables involving magnetic resonance
- G01R33/44—Arrangements or instruments for measuring magnetic variables involving magnetic resonance using nuclear magnetic resonance [NMR]
- G01R33/48—NMR imaging systems
- G01R33/54—Signal processing systems, e.g. using pulse sequences ; Generation or control of pulse sequences; Operator console
- G01R33/56—Image enhancement or correction, e.g. subtraction or averaging techniques, e.g. improvement of signal-to-noise ratio and resolution
- G01R33/563—Image enhancement or correction, e.g. subtraction or averaging techniques, e.g. improvement of signal-to-noise ratio and resolution of moving material, e.g. flow contrast angiography
- G01R33/56308—Characterization of motion or flow; Dynamic imaging
Landscapes
- Health & Medical Sciences (AREA)
- Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Engineering & Computer Science (AREA)
- General Health & Medical Sciences (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Molecular Biology (AREA)
- High Energy & Nuclear Physics (AREA)
- Biophysics (AREA)
- Biomedical Technology (AREA)
- Heart & Thoracic Surgery (AREA)
- Medical Informatics (AREA)
- Radiology & Medical Imaging (AREA)
- Surgery (AREA)
- Animal Behavior & Ethology (AREA)
- Pathology (AREA)
- Public Health (AREA)
- Veterinary Medicine (AREA)
- Signal Processing (AREA)
- Condensed Matter Physics & Semiconductors (AREA)
- General Physics & Mathematics (AREA)
- Artificial Intelligence (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Physiology (AREA)
- Psychiatry (AREA)
- Magnetic Resonance Imaging Apparatus (AREA)
Abstract
Es ist eine Aufgabe der vorliegenden Erfindung, die Lungenbildgebung derart zu verbessern, dass dynamische Prozesse selektiv dargestellt werden können. Das erfindungsgemäße Verfahren zur selektiven Darstellung einer Bewegung der Lunge kann die folgenden Schritte aufweisen: Zuerst werden Magnetresonanzbilder (MR-Bilder) der Lunge im zeitlichen Verlauf aufgenommen, d. h. es werden über mehrere Atemzyklen MR-Bilder der Lunge aufgenommen. Die aufgenommenen MR-Bilder werden bezüglich einer Referenzposition registriert, und es wird in ihnen der Signalverlauf über die Zeit bestimmt. Die geschieht vorteilhaft pixelweise, um die Auflösung der MR-Bilder voll ausnutzen zu können. Sodann wird, z. B. mittels einer Fouriertransformation, das Frequenzspektrum von den bestimmten Signalverläufen bestimmmt. Ein bestimmtes Frequenzspektrum wird mit einem Frequenzbandfilter gefiltert, wobei der Frequenzbereich des Frequenzbandfilters der darzustellenden Bewegung angepasst ist. Das gefilterte Frequenzspektrum wird in einen gefilterten Signalverlauf der MR-Bilder rücktransformiert und die durch diese Rücktransformation erhaltenen Magnetresonanzbilder werden mit dem gefilterten Signalverlauf im zeitlichen Verlauf angezeigt. Weiter werden ein Computerprogramm, eine Bildverarbeitungseinheit und ein Magnetresonanzgerät beansprucht, die es ermöglichen, dynamische Prozesse in einer Folge von MR-Bildern im zeitlichen Verlauf selektiv darstellen zu können.
Description
- Die Erfindung betrifft ein Verfahren zur selektiven Darstellung einer Bewegung der Lunge, ein Computerprogramm, eine Bildverarbeitungseinheit und ein Magnetresonanzgerät.
- Die Magnetresonanztechnik (im Folgenden steht „MR” für „Magnetresonanz”) ist eine bekannte Technik, mit der Bilder vom Inneren eines Untersuchungsobjektes erzeugt werden können. Vereinfacht ausgedrückt wird hierzu das Untersuchungsobjekt in einem MR-Gerät in einem vergleichsweise starken statischen, homogenen Grundmagnetfeld (Feldstärken von 0,2 Tesla bis 7 Tesla und mehr) positioniert, so dass sich dessen Kernspins entlang des Grundmagnetfeldes orientieren. Zum Auslösen von Kernspinresonanzen werden Sequenzen hochfrequenter Anregungspulse resonant in das Untersuchungsobjekt eingestrahlt, die ausgelösten Kernspinsignale gemessen und auf deren Basis MR-Bilder rekonstruiert. Die MR-Technik eignet sich insbesondere zur Bildgebung von Weichteilen, da hier ein besonders guter Kontrast erzielt werden kann.
- Auch die Bildgebung der Lunge mittels MR wird seit Jahren vorangetrieben, um vor allem Rückschlüsse auf die Funktion der Lunge zu erhalten. Dabei sind allerdings verschiedene Schwierigkeiten zu beachten. Einerseits sind die Protonendichten im Bereich der mit Luft gefüllten Lunge, und auch die Protonendichte des Lungenparenchyms selbst sehr gering was zu einem schlechten Signal-Rausch-Verhältnis („signal-to-noise ratio”: SNR) führt. Andererseits erschweren starke lokale Feldänderungen, die durch Suszeptibilitätsänderungen entstehen, die MR-Bildgebung. Weiter ist der Bereich der Lunge sowohl durch Atmung als auch durch Herzschlag des zu untersuchenden Patienten Bewegungen unterworfen, die zu Bewegungsartefakten in MR-Bildern führen. Dabei ändert insbesondere die Atembewegung nicht nur die Lage des aufgenommenen Gewebes, sondern auch das Gesamtvolumen der Lunge. Ein weiterer Effekt, der die qualitative Auswertung von Lungendaten erschwert, ist der starke, variierende Blutfluss in der Lunge innerhalb eines Herzzyklus.
- Um diesen Schwierigkeiten zu begegnen, werden MR-Bilder der Lunge beispielsweise unter Verstärkung durch Kontrastmittel, wie z. B. Helium oder Sauerstoff, durchgeführt. Ein Beispiel für Sauerstoff-verstärkte Bildgebung ist z. B. Vu M. Mai et al. „Influence of Oxygen Flow Rate an Signal and T1 Changes in Oxygen-Enhanced Ventilation Imaging" im Journal of Magnetic Resoncance Imaging 16: 37–41, 2002.
- Dank der mittlerweile erreichbaren kurzen Aufnahmezeiten für MR-Bilder, ist es auch möglich, mittels MR-Tomographie Aufnahmen der Lunge ohne starke Bewegungsartefakte zu gewinnen. Dabei muss jedoch oftmals eine Triggerung des Sequenzablaufs mit der äußeren Bewegung, der Atmung bzw. dem Herzschlag, erfolgen, oder die Atmung muss für die Dauer der Aufnahme angehalten werden. Ein Beispiel für eine Herzschlag-getriggerte Bildgebung ist z. B. Knight-Scott J. et al. „Temporal Dynamics of Blood Flow Effects in Half-Fourier Fast Spin Echo 1H Magnetic Resonance Imaging of the Human Lung" im Journal of Magnetic Resonance Imaging 14: 411–418, 2001.
- Bei Grundfeldstärken von größer als 0,7 T ist die T2*-Relaxationszeit suszeptibilitätsbedingt sehr kurz (M. Deimling in Proc. Int. Soc. Mang. Reson. Med. 8 (2000), No. 2202). Daher werden zur Darstellung des Lungenparenchyms Spinecho basierte Sequenzen, wie z. B. single shot HASTE, verwendet, da Spinechosequenzen den Einfluss von Suszeptibilitätsänderungen unterdrücken (siehe z. B. wieder Vu M. Mai et al. im Journal of Magnetic Resoncance Imaging 16: 37–41, 2002). Mit neuzeitlichen Tomographen können sehr kurze Echozeiten erreicht werden, so das auch bei Grundfeldstärken bis etwa 1.5 T Gradientenechosequenzen angewendet werden können (siehe z. B. Marcus J. T. et. al. in Proc. Int. Soc. Mag. Reson. Med. (2007), No. 2777).
- Bei Auswertungen von Bildfolgen der Lunge, z. B. im zeitlichen Verlauf, tritt, wie oben beschrieben, das Problem auf, dass durch Einatmen und Ausatmen sich die Größe der Lunge verändert. Hierdurch kommt es zu einer Verschiebung korrespondierender Bereiche in den verschiedenen MR-Bildern, da die Lunge jeweils eine unterschiedliche Lage in verschiedenen MR-Bildern einnimmt. Es sind jedoch im Stand der Technik Verfahren bekannt, die unterschiedlichen Lungenbildgrößen im Atemzyklus auf eine Referenzgröße fixieren können, und mit denen so eine Registrierung der einzelnen Bilder zueinander möglich ist (H. G. Topf et al. in Proc. Int. Soc. Magn. Reson. Med 11 (2004) No. 671). Damit ist es möglich geworden, Signaländerungen durch Dichteänderungen als Funktion der Zeit zu analysieren. Durch das Ein- und Ausatmen verändert sich das signalgebende Volumen, so dass die Atembewegung zur Signaländerungen führt. Durch die Darstellung der Signaländerungen können beispielsweise Ventilationsdefekte in der Lunge sichtbar gemacht werden. Ein hierbei auftretendes Problem ist aber weiterhin, dass die Signaländerungen nicht nur von der sich ändernden Parenchymdichte abhängen, sondern auch von dem Signal des Blutes beeinflusst werden, das sich im Herzrhythmus in und von den Bronchialgefäßen ausbreitet. Dieses Blutsignal stört die Analyse der Lungensignaländerung.
- Aus der
DE 10 2005 010 093 A1 ist ein Verfahren zur Darstellung von Respirationsmustern durch Trennung der Signalanteile des Blutes von Signalanteilen des Lungenparenchyms in einem Magnetresonanzbild bekannt, das die folgenden Schritte umfasst: – Aufnehmen von Magnetresonanzbildern der Lunge im zeitlichen Verlauf, – Berechnen des Signalverlaufs in den Magnetresonanzbildern über die Zeit, – Bilden der Fouriertransformation des zeitlichen Signalverlaufs, – Extrahieren des zu dem Lungenparenchym gehörigen Fourierspektrums, – Darstellen einer in dem Fourierspektrum enthaltenen Information im Magnetresonanzbild. Dabei wird die gewonnene Infor mation aus dem Spektrum statisch einem anatomischen Abbild der Lunge überlagert. - Ein anderes Beispiel für eine Auswertung von Spektren in MR-Bildern ist beschrieben in Weisskoff, R. M., Baker, J., Belliveau, J., Davis, T. L., Kwong, K. K., Cohen, M. S., & Rosen, B. R., „Power spectrum analysis of functionally-weighted MR data: What's in the noise?", Proceedings of the Society for Magnetic Resonance, 1, 7 (1993).
- Es ist eine Aufgabe der vorliegenden Erfindung, die Lungenbildgebung derart zu verbessern, dass dynamische Prozesse selektiv dargestellt werden können.
- Diese Aufgabe wird erfindungsgemäß durch ein Verfahren gemäß Patentanspruch 1 gelöst. In den abhängigen Ansprüchen sind bevorzugte Ausführungsformen der Erfindung beschrieben.
- Das erfindungsgemäße Verfahren zur selektiven Darstellung einer Bewegung der Lunge kann die folgenden Schritte aufweisen: zuerst werden Magnetresonanzbilder der Lunge im zeitlichen Verlauf aufgenommen, d. h. es werden über mehrere Atemzyklen Magnetresonanzbilder (MR-Bilder) der Lunge aufgenommen. Die aufgenommenen MR-Bilder werden bezüglich einer Referenzposition registriert, und es wird in den aufgenommenen MR-Bildern der Signalverlauf über die Zeit bestimmt. Dies geschieht vorteilhaft pixelweise, um die Auflösung der MR-Bilder voll ausnutzen zu können. Sodann wird, z. B. mittels einer Fouriertransformation, das Frequenzspektrum von den bestimmten Signalverläufen bestimmt. Ein bestimmtes Frequenzspektrum wird mit einem Frequenzbandfilter gefiltert, wobei der Frequenzbereich des Frequenzbandfilters der darzustellenden Bewegung angepasst ist. Das gefilterte Frequenzspektrum wird in einen gefilterten Signalverlauf der MR-Bilder rücktransformiert und die durch diese Rücktransformation erhaltenen Magnetresonanzbilder werden mit dem gefilterten Signalverlauf im zeitlichen Verlauf angezeigt.
- Durch das angegebene Verfahren ist es möglich, selektiv die Dynamik von Lungenparenchym oder Blut eines perfundierten Gewebes ohne Kontrastmittel darzustellen, also ohne externe Agentien wie z. B. Sauerstoff oder Helium. Dabei ist es nicht erforderlich, dass ein Patient den Atem anhält, sondern sie bzw. er kann normal atmen.
- Es ist eine weitere Aufgabe der vorliegenden Erfindung, ein Computerprogramm, eine Bildverarbeitungseinheit und ein Magnetresonanzgerät anzugeben, die es ermöglichen dynamische Prozesse in einer Folge von MR-Bildern im zeitlichen Verlauf selektiv dargestellen zu können.
- Die weitere Aufgabe wird gelöst durch ein Computerprogramm gemäß Anspruch 9, das ein vorbeschriebenes Verfahren auf einer Recheneinheit implementiert, wenn es auf der Recheneinheit ausgeführt wird, durch eine Bildverarbeitungseinheit gemäß Anspruch 10, die zur Durchführung eines vorbeschreibenen Verfahrens ausgebildet ist, und durch ein Magnetresonanzgerät gemäß Anspruch 11 mit einer vorgenannten Bildverarbeitungseinheit. Die oben beschriebenen Vorteile gelten hierbei analog.
- Weitere Vorteile und Einzelheiten der vorliegenden Erfindung ergeben sich aus den im Folgenden beschriebenen Ausführungsbeispielen sowie anhand der Zeichnungen. Die aufgeführten Beispiele stellen keine Beschränkung der Erfindung dar. Es zeigen:
-
1 ein Flussdiagramm zur Verdeutlichung einer beispielhaften Ausführungsform eines Verfahrens zur selektiven Darstellung einer Bewegung der Lunge, -
2 einen beispielhafter Signalverlauf, und -
3A –B Frequenzspektren, die während des Verfahrens auftreten können. -
1 zeigt ein Flussdiagramm eines beispielhaften Verfahrens zur selektiven Darstellung einer Bewegung der Lunge. Dabei werden in einem Schritt101 zunächst mehrere MR-Bilder der Lunge im zeitlichen Verlauf aufgenommen. Beispielsweise wird eine Reihe von N = 100 MR-Bildern oder mehr, zumindest aber N = 20 MR-Bilder, der Lunge im zeitlichen Verlauf einer Bewegung der Lunge aufgenommen. - Die Wahl der bei der Aufnahme dieser MR-Bilder verwendeten Messsequenz hängt üblicherweise von der Wahl des verwendeten Magnetresonanzgeräts ab, hier insbesondere von der Grundfeldstärke. Durch die Anwendung von sehr kurzen Echozeiten TE können auch Gradientenechosequenzen und auch SSFP-(Steady State Free Precession)Sequenzen vorteilhaft verwendet werden, wobei die Messzeit pro MR Bild ebenfalls im Bereich von ca. 0,1–0,5 sec, vorzugsweise 0,3 sec liegen sollte.
- In einem weiteren Schritt
102 werden die aufgenommenen Bilder registriert. Dies geschieht z. B. wie in der bereits zitierten Veröffentlichung von H. G. Topf et al. beschrieben, durch eine elastische Pixeltransformation der Einzelbilder die sich auf eine Referenzgröße beziehen. Damit können gleiche Orte in den einzelnen MR-Bildern mit guter Zuverlässigkeit als gleiche Orte in dem aufgenommenen Untersuchungsobjekt, hier der Lunge, angenommen werden. - In einem weiteren Schritt
103 wird der Signalverlauf21 in den MR-Bildern über die Zeit bestimmt. Dazu wird vorteilhaft pixelweise ein Verlauf einer Intensität des MR-Signals über die Zeit in den MR-Bildern bestimmt. Dabei kann es bereits ausreichend sein, die Signalverläufe der einzelnen Pixel der MR-Bilder nur in einem interessierenden Bereich in jedem MR-Bild, der z. B. mittels einer Segmentierungsmethode bestimmt wurde, zu bestimmen, um Aufwand und Zeit weiter zu verringern. Als interessierender Bereich kommen insbesondere die abgebildete Lunge und/oder die Bronchialgefäße in Frage. Alternativ kann der Signalverlauf auch über Gruppen von Pixeln der MR-Bilder gemittelt werden, wobei dies zu Verlusten in der Auflösung führt. -
2 zeigt beispielhaft einen derartigen Signalverlauf21 eines Pixels i über der Zeit. - Der bestimmte Signalverlauf
21 könnte beispielsweise sinusähnliche Schwingungen enthalten, die sich aufgrund periodischer Bewegungen in dem aufgenommenen Objekt, hier der Lunge ergeben. Typische periodische Bewegungen sind insbesondere Bewegungen durch den Blutfluss und durch die aufgrund der Atmung des Patienten verursachte Dichteänderung des Lungenparenchyms. - Um die in den MR-Bildern enthaltenen Informationen bezüglich der verschiedenen Bewegungen trennen zu können, wird in einem weiteren Schritt
104 das Frequenzspektrum des Signalverlaufs21 bestimmt. Dieser Schritt kann vorteilhaft mittels einer Fourieranalyse durchgeführt werden. Es sind aber auch andere gängige Methoden wie z. B. eine Zerlegung in Sinus- oder sinλ(t)-Funktionen anwendbar. In dem Frequenzspektrum sind einzelne Frequenzanteile des zeitlichen Signalverlaufs21 als Maxima erkennbar. - Ein Beispiel für ein derartiges Frequenzspektrum ist in der später näher beschriebenen
3 gegeben. - Auf ein gewonnenes Frequenzspektrum wird in einem weiteren Schritt
105 ein Frequenzbandfilter35 angewendet, dessen Frequenzbereich der darzustellenden Bewegung angepasst ist. Somit kann der Frequenzanteil der darzustellenden Bewegung aus dem gesamten Frequenzspektrum herausgefiltert werden. Die darzustellende Bewegung kann dabei beispielsweise eine Atembewegung oder eine durch eine Herzaktivität verursachte Bewegung sein. - In einem weiteren Schritt
106 wird das gefilterte Frequenzspektrum34 wieder in einen, dann gefilterten, Signalverlauf rücktransformiert. Dabei wird vorteilhaft die inverse Methode zu der in Schritt104 angewandten Methode zur Bestimmung des Frequenzspektrums angewendet. War der Frequenzbereich des Frequenzbandfilters aus Schritt105 auf die Atembewegung angepasst, entspricht der durch die Rücktransformation erhaltene gefilterte Signalverlauf einer Signaländerung aufgrund der Atembewegung. War der Frequenzbereich des Frequenzbandfilters aus Schritt105 auf die durch die Herzaktivität verursachte Bewegung angepasst, entspricht der durch die Rücktransformation erhaltene gefilterte Signalverlauf einer Signaländerung aufgrund der Herzaktivität. Ein beispielhafter Signalverlauf22 in einem Pixel i der MR-Bilder für eine Signaländerung durch eine Atembewegung und ein beispielhafter Signalverlauf23 in einem Pixel i der MR-Bilder für eine Signaländerung durch Herzaktivität sind in2 dargestellt. - Die so gewonnenen gefilterten Signalverläufe in den MR-Bildern werden in einem letzten Schritt
107 wieder als MR-Bilder' mit dem gefilterten Signalverlauf im zeitlichen Verlauf angezeigt. Auf diese Weise erhält man einfach eine Folge von MR-Bildern' im zeitlichen Verlauf, in denen die Dynamik der durch den Frequenzbandfilter herausgefilterten Bewegung und die Anatomie der Lunge gleichzeitig erkennbar sind. Ein Überlagern der aus dem Frequenzspektrum gewonnenen Informationen mit einem anatomischen Bild kann somit entfallen. Darüber hinaus ist die Darstellung der Bewegung in der Lunge nicht statisch, sondern dynamisch über den zeitlichen Verlauf gegeben, was eine genauere Analyse der Funktion der Lunge ermöglicht. - Das in Schritt
104 bestimmte Frequenzspektrum kann auch nacheinander mit einem ersten Frequenzbandfilter, der z. B. auf eine Atembewegung angepasst ist und mit einem zweiten Frequenzbandfilter, der auf eine durch die Herzaktivität verursachte Bewegung angepasst ist, gefiltert werden, sodass je zwei gefilterte Signalverläufe erzeugt werden. So kann die Dynamik der verschiedenen Bewegungen unabhängig voneinander in den jeweiligen Schritten107 angezeigt werden. Dabei sind die einzelnen Darstellungen jeweils frei von Einflüssen anderer, nicht darzustellender Bewegungen, auf die der entsprechende Frequenzbandfilter nicht angepasst war. - Die in Schritt
105 zu dem gewählten Frequenzbereich zugehörenden frequenzsortierten MR-Bilder können vor der Rücktransformation in Schritt106 durch Bildung z. B. der Quadratsumme der durch die Filterkurve gewichteten MR-Bilder zu einem resultierenden Bild bei der Mittenfrequenz des Frequenzbandfilters zusammengefasst werden. - In
3 oben ist ein typisches Frequenzspektrum mit Maxima31 ,32 und33 dargestellt. Dabei liegt das Maximum31 , das dem Gleichanteil des Spektrums entspricht, bei einer Frequenz ν = 0. Ein weiteres Maximum33 , das dem, durch die Herzfrequenz verursachten Frequenzanteil des Blutes in dem Signalverlauf21 entspricht liegt bei einer höheren Frequenz als ein Maximum32 , das dem Frequenzanteil der Änderung der Dichte des Lungenparenchyms in dem Signalverlauf21 , die sich durch die im Vergleich zur Herzfrequenz niederfrequenten Atmung ergibt, entspricht. Wie durch die gepunktete Darstellung des Maximums31 angedeutet, kann der Gleichanteil vorteilhaft von dem bestimmten Frequenzspektrum subtrahiert werden. - Ein Frequenzbandfilter
35 kann vorteilhaft eine Gaußfunktion mit einem Frequenzbereich36 auf einen einer darzustellenden Bewegung angepassten Frequenzanteil sein. Dabei wird der, der Bewegung angepasste, Frequenzanteil vorteilhaft aus dem Frequenzspektrum, als Frequenz des Maximums der darzustellenden Bewegung bestimmt. Auch der Frequenzbereich des Frequenzbandfilters kann vorteilhaft durch die bestimmten Frequenzanteile in dem Frequenzspektrum des bestimmten Signalverlaufs ermittelt werden. Dies geschieht z. B. vorteilhaft durch eine Bestimmung der Breite des betreffenden Maximums im Frequenzspektrum. Alternativ kann auch ein gewisser Toleranzbereich als Frequenzbereich vorgegeben werden. Dabei schließt der Frequenzbereich des Frequenzbandfilters den Frequenzanteil des Signalverlaufs ein, der der darzustellenden Bewegung entspricht. - Es ist auch möglich, den Frequenzanteil, den der Frequenzbandfilter aus dem Frequenzspektrum filtern soll durch eine externe Messung der darzustellenden Bewegung, z. B. durch Pulsmessung bzw. Atembewegungsmessung, etwa durch einen Atemgürtel, zu bestimmen. Dies ist aber umständlicher und nicht notwendig.
- In
3 ist der Frequenzbandfilter35 als Gaußfunktion über dem Maximum, das den Frequenzanteil einer durch die Herzaktivität verursachten Bewegung entspricht, dargestellt. Das mit diesem Frequenzbandfilter35 gefilterte Spektrum34 enthält somit nur noch ein Maximum34 , das dem Maximum33 , das der durch die Herzaktivität verursachten Bewegung entspricht, entspricht. - Analog hat ein, mit einem Frequenzbandfilter, der über dem der Atembewegung entsprechenden Frequenzanteil liegt, gefiltertes Frequenzspektrum nur noch ein Maximum, das dem Maximum
32 entspricht, das der Atembewegung entspricht. Dies ist der Übersichtlichkeit halber nicht noch einmal explizit dargestellt. - Wie bereits oben erwähnt, wird der Signalverlauf
21 der MR-Bilder vorteilhaft pixelweise bestimmt, und die Analyse des Signalverlaufs21 , Filterung und Rücktransformation in einen gefilterten Signalverlauf22 bzw.23 wird ebenfalls pixelweise vollzogen. Auf diese Weise erhält man MR-Bilder' mit gefiltertem Signalverlauf, die dieselbe Auflösung haben wie die ursprünglich aufgenommenen MR-Bilder, die aber selektiv entweder den Einfluss der Atembewegung oder den Einfluss der durch die Herzaktivität verursachten Bewegung in der Lunge darstellen. - ZITATE ENTHALTEN IN DER BESCHREIBUNG
- Diese Liste der vom Anmelder aufgeführten Dokumente wurde automatisiert erzeugt und ist ausschließlich zur besseren Information des Lesers aufgenommen. Die Liste ist nicht Bestandteil der deutschen Patent- bzw. Gebrauchsmusteranmeldung. Das DPMA übernimmt keinerlei Haftung für etwaige Fehler oder Auslassungen.
- Zitierte Patentliteratur
-
- - DE 102005010093 A1 [0008]
- Zitierte Nicht-Patentliteratur
-
- - Vu M. Mai et al. „Influence of Oxygen Flow Rate an Signal and T1 Changes in Oxygen-Enhanced Ventilation Imaging” im Journal of Magnetic Resoncance Imaging 16: 37–41, 2002 [0004]
- - Knight-Scott J. et al. „Temporal Dynamics of Blood Flow Effects in Half-Fourier Fast Spin Echo 1H Magnetic Resonance Imaging of the Human Lung” im Journal of Magnetic Resonance Imaging 14: 411–418, 2001 [0005]
- - M. Deimling in Proc. Int. Soc. Mang. Reson. Med. 8 (2000), No. 2202 [0006]
- - Vu M. Mai et al. im Journal of Magnetic Resoncance Imaging 16: 37–41, 2002 [0006]
- - Marcus J. T. et. al. in Proc. Int. Soc. Mag. Reson. Med. (2007), No. 2777) [0006]
- - H. G. Topf et al. in Proc. Int. Soc. Magn. Reson. Med 11 (2004) No. 671 [0007]
- - Weisskoff, R. M., Baker, J., Belliveau, J., Davis, T. L., Kwong, K. K., Cohen, M. S., & Rosen, B. R., „Power spectrum analysis of functionally-weighted MR data: What's in the noise?”, Proceedings of the Society for Magnetic Resonance, 1, 7 (1993) [0009]
- - H. G. Topf et al. [0022]
Claims (11)
- Verfahren zur selektiven Darstellung einer Bewegung der Lunge umfassend folgende Schritte: – Aufnehmen einer Reihe von Magnetresonanzbilder der Lunge im zeitlichen Verlauf, – Registrieren der Magnetresonanzbilder bezüglich einer Referenzposition, – Bestimmen des Signalverlaufs in den Magnetresonanzbildern über die Zeit, – Bestimmen eines Frequenzspektrums des Signalverlaufs, – Filtern des Frequenzspektrums mit einem Frequenzbandfilter, dessen Frequenzbereich der darzustellenden Bewegung angepasst ist, – Rücktransformieren des gefilterten Frequenzspektrums in einen gefilterten Signalverlauf der Magnetresonanzbilder, – Anzeigen der durch die Rücktransformation erhaltenen Magnetresonanzbilder mit gefiltertem Signalverlauf im zeitlichen Verlauf.
- Verfahren nach Anspruch 1, wobei der Frequenzbereich des Frequenzbandfilters durch eine Bestimmung von Frequenzanteilen in dem bestimmten Signalverlauf ermittelt wird.
- Verfahren nach Anspruch 2, wobei die Frequenzanteile des Signalverlaufs mittels einer Fourieranalyse des Signalverlaufs ermittelt werden.
- Verfahren nach einem der Ansprüche 2 bis 3, wobei der Frequenzbereich des Frequenzbandfilters einen Frequenzanteil des Signalverlaufs einschließt, der der darzustellenden Bewegung entspricht.
- Verfahren nach Anspruch 4, wobei der Frequenzbandfilter eine Gaussfunktion des Frequenzanteils des Signalverlaufs ist, der der darzustellenden Bewegung entspricht.
- Verfahren nach einem der vorhergehenden Ansprüche, wobei der Signalverlauf über die Zeit in den Magnetresonanzbildern pixelweise bestimmt wird.
- Verfahren nach einem der vorhergehenden Ansprüche, wobei die darzustellende Bewegung eine Atembewegung ist.
- Verfahren nach einem der vorhergehenden Ansprüche, wobei die darzustellende Bewegung eine durch eine Herzaktivität verursachte Bewegung ist.
- Computerprogramm, das ein Verfahren nach einem der Ansprüche 1 bis 8 auf einer Recheneinheit implementiert, wenn es auf der Recheneinheit ausgeführt wird.
- Bildverarbeitungseinheit, die zur Durchführung eines Verfahrens nach einem der Ansprüche 1 bis 8 ausgebildet ist.
- Magnetresonanzgerät mit einer Bildverarbeitungseinheit, die zur Durchführung eines Verfahrens gemäß einem der Ansprüche 1 bis 8 ausgebildet ist.
Priority Applications (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
DE102008049709.6A DE102008049709B4 (de) | 2008-09-30 | 2008-09-30 | Verfahren zur selektiven Darstellung einer Bewegung der Lunge, Computerprogramm, Bildverarbeitungseinheit und Magnetresonanzgerät |
US12/570,156 US8154288B2 (en) | 2008-09-30 | 2009-09-30 | Method, processor, and magnetic resonance apparatus for selective presentation of lung movement |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
DE102008049709.6A DE102008049709B4 (de) | 2008-09-30 | 2008-09-30 | Verfahren zur selektiven Darstellung einer Bewegung der Lunge, Computerprogramm, Bildverarbeitungseinheit und Magnetresonanzgerät |
Publications (2)
Publication Number | Publication Date |
---|---|
DE102008049709A1 true DE102008049709A1 (de) | 2010-04-22 |
DE102008049709B4 DE102008049709B4 (de) | 2019-12-12 |
Family
ID=42034765
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
DE102008049709.6A Active DE102008049709B4 (de) | 2008-09-30 | 2008-09-30 | Verfahren zur selektiven Darstellung einer Bewegung der Lunge, Computerprogramm, Bildverarbeitungseinheit und Magnetresonanzgerät |
Country Status (2)
Country | Link |
---|---|
US (1) | US8154288B2 (de) |
DE (1) | DE102008049709B4 (de) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US9226684B2 (en) | 2011-10-28 | 2016-01-05 | Siemens Aktiengesellschaft | Method for recording magnetic resonance image data when using a ventilation apparatus |
Families Citing this family (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2010138919A2 (en) | 2009-05-28 | 2010-12-02 | Angiodynamics, Inc. | System and method for synchronizing energy delivery to the cardiac rhythm |
US9895189B2 (en) | 2009-06-19 | 2018-02-20 | Angiodynamics, Inc. | Methods of sterilization and treating infection using irreversible electroporation |
GB2474233A (en) | 2009-10-06 | 2011-04-13 | Uk Investments Associates Llc | Cooling pump comprising a detachable head portion |
EP2627274B1 (de) | 2010-10-13 | 2022-12-14 | AngioDynamics, Inc. | System zur elektrischen ableitung des gewebes eines patienten |
US9078665B2 (en) | 2011-09-28 | 2015-07-14 | Angiodynamics, Inc. | Multiple treatment zone ablation probe |
JP6185929B2 (ja) * | 2011-12-15 | 2017-08-23 | コーニンクレッカ フィリップス エヌ ヴェKoninklijke Philips N.V. | Mpiにおけるバックグラウンド信号の除去 |
CN105353420B (zh) * | 2015-11-17 | 2017-08-22 | 中煤科工集团西安研究院有限公司 | 矿井音频电透视数据的一致性校正方法 |
US10905492B2 (en) | 2016-11-17 | 2021-02-02 | Angiodynamics, Inc. | Techniques for irreversible electroporation using a single-pole tine-style internal device communicating with an external surface electrode |
JP6848393B2 (ja) * | 2016-11-25 | 2021-03-24 | コニカミノルタ株式会社 | 動態画像処理装置 |
TW202116253A (zh) * | 2019-07-04 | 2021-05-01 | 新加坡商派拉梅維爾私人有限公司 | 診斷支援程式 |
Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
DE102005010093A1 (de) | 2005-03-04 | 2006-09-14 | Siemens Ag | Verfahren zur Darstellung von Respirationsmustern |
Family Cites Families (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
FR2709656B1 (fr) * | 1993-09-07 | 1995-12-01 | Deemed Int Sa | Installation pour opération de microchirurgie assistée par ordinateur et procédés mis en Óoeuvre par ladite installation. |
JP2004520865A (ja) * | 2000-07-25 | 2004-07-15 | リタ メディカル システムズ インコーポレイテッド | 局在化インピーダンス測定を使用する腫瘍の検出および処置のための装置 |
ATE320767T1 (de) * | 2001-09-28 | 2006-04-15 | Rita Medical Systems Inc | Impedanzgesteuerte vorrichtung zur ablation von gewebe |
US20050033157A1 (en) * | 2003-07-25 | 2005-02-10 | Klein Dean A. | Multi-modality marking material and method |
-
2008
- 2008-09-30 DE DE102008049709.6A patent/DE102008049709B4/de active Active
-
2009
- 2009-09-30 US US12/570,156 patent/US8154288B2/en active Active
Patent Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
DE102005010093A1 (de) | 2005-03-04 | 2006-09-14 | Siemens Ag | Verfahren zur Darstellung von Respirationsmustern |
Non-Patent Citations (11)
Title |
---|
Biswal B., DeYoe E.A., Hyde J.S.: "Reduction of physiological fluctuations in fMRI using digital filters. IN: Magnetic Resonance in Medicine, 1996, Vol. 35, S.107-113 * |
H. G. Topf et al. in Proc. Int. Soc. Magn. Reson. Med 11 (2004) No. 671 |
Knight-Scott J. et al. "Temporal Dynamics of Blood Flow Effects in Half-Fourier Fast Spin Echo 1H Magnetic Resonance Imaging of the Human Lung" im Journal of Magnetic Resonance Imaging 14: 411-418, 2001 |
M. Deimling in Proc. Int. Soc. Mang. Reson. Med. 8 (2000), No. 2202 |
Marcus J. T. et. al. in Proc. Int. Soc. Mag. Reson. Med. (2007), No. 2777) |
Mrsic-Flogel T., Hubener M., Bonhoeffer T.: "Brain Mapping: New Wave Optical Imaging "IN: Current Biology, 2003, Vol. 13, Nr. 19, S. R778-R780 * |
Vu M. Mai et al. "Influence of Oxygen Flow Rate an Signal and T1 Changes in Oxygen-Enhanced Ventilation Imaging" im Journal of Magnetic Resoncance Imaging 16: 37-41, 2002 |
Vu M. Mai et al. im Journal of Magnetic Resoncance Imaging 16: 37-41, 2002 |
Weisskoff R. M., Baker J., Belliveau J., Davis T.L, Kwong K.K., Cohen M.S., Rosen B.R: "Power spectrum analysis of functionally- weighted MR data: Wath's in the noise?", Proceedings of the Society for Magnetic Resonance, 1993, Vol. 1, S.7 * |
Weisskoff R. M., Baker J., Belliveau J., Davis T.L, Kwong K.K., Cohen M.S., Rosen B.R: "Power spectrum analysis of functionally- weighted MR data: Wath's in the noise?", Proceedings of the Society for Magnetic Resonance, 1993, Vol. 1, S.7 Biswal B., DeYoe E.A., Hyde J.S.: "Reduction of physiological fluctuations in fMRI using digital filters. IN: Magnetic Resonance in Medicine, 1996, Vol. 35, S.107-113 Mrsic-Flogel T., Hubener M., Bonhoeffer T.: "Brain Mapping: New Wave Optical Imaging "IN: Current Biology, 2003, Vol. 13, Nr. 19, S. R778-R780 |
Weisskoff, R. M., Baker, J., Belliveau, J., Davis, T. L., Kwong, K. K., Cohen, M. S., & Rosen, B. R., "Power spectrum analysis of functionally-weighted MR data: What's in the noise?", Proceedings of the Society for Magnetic Resonance, 1, 7 (1993) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US9226684B2 (en) | 2011-10-28 | 2016-01-05 | Siemens Aktiengesellschaft | Method for recording magnetic resonance image data when using a ventilation apparatus |
Also Published As
Publication number | Publication date |
---|---|
US8154288B2 (en) | 2012-04-10 |
DE102008049709B4 (de) | 2019-12-12 |
US20100090696A1 (en) | 2010-04-15 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
DE102008049709B4 (de) | Verfahren zur selektiven Darstellung einer Bewegung der Lunge, Computerprogramm, Bildverarbeitungseinheit und Magnetresonanzgerät | |
DE102009055961B4 (de) | Verfahren zur Erfassung von einer Atembewegung eines Untersuchungsobjektes entsprechender Signaldaten mittels Magnetresonanz | |
DE102011007850B4 (de) | Verfahren zur getriggerten Akquisition von Messdaten mittels eines Magnetresonanzgeräts unter Unterdrückung von unerwünschten Signalen, Magnetresonanzgerät, Computerprogramm und elektronisch lesbarer Datenträger | |
DE102012206585B4 (de) | Verfahren zur schnellen ortsaufgelösten Bestimmung eines Magnetresonanz-Relaxationsparameters in einem Untersuchungsgebiet | |
DE102007035176B4 (de) | Verfahren zur Aufzeichnung und Verarbeitung einer Folge von zeitlich aufeinander folgenden Bilddatensätzen sowie Magnet-Resonanz-Gerät | |
DE4428503C2 (de) | Diffusionsgewichtete Bildgebung mit magnetischer Resonanz | |
DE102012215718B4 (de) | Verfahren und Magnetresonanzanlage zur MR-Bildgebung eines vorbestimmten Volumenabschnitts eines lebenden Untersuchungsobjekts mittels Stimulieren des Untersuchungsobjekts | |
DE4432570A1 (de) | Verfahren und Vorrichtung für die Kernresonanzabbildung physiologischer Funktionsinformation | |
DE102011083619B4 (de) | Verfahren zur Erzeugung einer Serie von MR-Bildern zur Überwachung einer Position eines in einem Untersuchungsgebiet befindlichen Interventionsgeräts, Magnetresonanzanlage und elektronisch lesbarer Datenträger | |
DE102014217284A1 (de) | Verfahren zu einer Magnetresonanz-Untersuchung eines Untersuchungsobjekts | |
DE102007011807B3 (de) | Sequenz für die Magnet-Resonanz-Bildgebung und Magnet-Resonanz-Gerät hierzu | |
DE102008014059B4 (de) | Verfahren zur Erstellung eines fettreduzierten ortsaufgelösten Magnetresonanzspektrums und Magnetresonanzanlage hierfür | |
DE102014209351B4 (de) | Magnetresonanz-Spektroskopieverfahren mit kurzer Echozeit, Magnetresonanzanlage und digitales Speichermedium | |
DE102012213551A1 (de) | Verfahren zur bewegungsgemittelten Schwächungskorrektur und Magnetresonanz-Anlage | |
DE102016207641A1 (de) | Parallele Magnetresonanz-Akquisitionstechnik | |
EP3382413A1 (de) | 2d navigatortechnik in der mri | |
DE102007030746B4 (de) | Verfahren zum Erstellen eines angiographischen Bildes und Magnet-Resonanz-Gerät | |
DE102014210218A1 (de) | Verfahren zum Generieren von medizinischen Bilddatensätzen | |
DE102014206724B4 (de) | Dynamische Bildgebung mit variablem Kontrast | |
DE102015200850B4 (de) | Verfahren zur Auswertung von medizinischen Bilddaten | |
DE102008039340A1 (de) | Verfahren zur Aufzeichnung und Darstellung von Kalibrierungsbildern sowie Magnet-Resonanz-Gerät | |
DE102011003606A1 (de) | Verfahren zum Entfernen unerwünschter perfusionsbedingter Signale aus MR-Bildern, Magnetresonanzanlage, Computerprogrammprodukt und elektronisch lesbarer Datenträger | |
EP3290940B1 (de) | Iterative rekonstruktion von quantitativen mr-bildern | |
DE102008039581B4 (de) | Verfahren zur Erstellung einer Magnet-Resonanz-Angiographie und Magnet-Resonanz-Gerät | |
DE10256208B4 (de) | Verfahren zur verbesserten Flussmessung in der Magnetresonanz-Tomographie |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
OP8 | Request for examination as to paragraph 44 patent law | ||
R081 | Change of applicant/patentee |
Owner name: SIEMENS HEALTHCARE GMBH, DE Free format text: FORMER OWNER: SIEMENS AKTIENGESELLSCHAFT, 80333 MUENCHEN, DE |
|
R018 | Grant decision by examination section/examining division | ||
R020 | Patent grant now final |