DE102018200900A1 - Artefaktvermeidung bei schnellen 3D Spinechosequenzen - Google Patents

Artefaktvermeidung bei schnellen 3D Spinechosequenzen Download PDF

Info

Publication number
DE102018200900A1
DE102018200900A1 DE102018200900.7A DE102018200900A DE102018200900A1 DE 102018200900 A1 DE102018200900 A1 DE 102018200900A1 DE 102018200900 A DE102018200900 A DE 102018200900A DE 102018200900 A1 DE102018200900 A1 DE 102018200900A1
Authority
DE
Germany
Prior art keywords
raw data
data set
complete
recorded
raw
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
DE102018200900.7A
Other languages
English (en)
Other versions
DE102018200900B4 (de
Inventor
Flavio Carinci
Dominik Paul
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Siemens Healthineers Ag De
Original Assignee
Siemens Healthcare GmbH
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Siemens Healthcare GmbH filed Critical Siemens Healthcare GmbH
Priority to DE102018200900.7A priority Critical patent/DE102018200900B4/de
Priority to US16/251,424 priority patent/US10605879B2/en
Priority to CN201910053157.3A priority patent/CN110068782B/zh
Publication of DE102018200900A1 publication Critical patent/DE102018200900A1/de
Application granted granted Critical
Publication of DE102018200900B4 publication Critical patent/DE102018200900B4/de
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01RMEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
    • G01R33/00Arrangements or instruments for measuring magnetic variables
    • G01R33/20Arrangements or instruments for measuring magnetic variables involving magnetic resonance
    • G01R33/44Arrangements or instruments for measuring magnetic variables involving magnetic resonance using nuclear magnetic resonance [NMR]
    • G01R33/48NMR imaging systems
    • G01R33/54Signal processing systems, e.g. using pulse sequences ; Generation or control of pulse sequences; Operator console
    • G01R33/56Image enhancement or correction, e.g. subtraction or averaging techniques, e.g. improvement of signal-to-noise ratio and resolution
    • G01R33/561Image 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
    • G01R33/5615Echo train techniques involving acquiring plural, differently encoded, echo signals after one RF excitation, e.g. using gradient refocusing in echo planar imaging [EPI], RF refocusing in rapid acquisition with relaxation enhancement [RARE] or using both RF and gradient refocusing in gradient and spin echo imaging [GRASE]
    • G01R33/5618Echo train techniques involving acquiring plural, differently encoded, echo signals after one RF excitation, e.g. using gradient refocusing in echo planar imaging [EPI], RF refocusing in rapid acquisition with relaxation enhancement [RARE] or using both RF and gradient refocusing in gradient and spin echo imaging [GRASE] using both RF and gradient refocusing, e.g. GRASE
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/05Detecting, measuring or recording for diagnosis by means of electric currents or magnetic fields; Measuring using microwaves or radio waves 
    • A61B5/055Detecting, 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
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01RMEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
    • G01R33/00Arrangements or instruments for measuring magnetic variables
    • G01R33/20Arrangements or instruments for measuring magnetic variables involving magnetic resonance
    • G01R33/44Arrangements or instruments for measuring magnetic variables involving magnetic resonance using nuclear magnetic resonance [NMR]
    • G01R33/48NMR imaging systems
    • G01R33/4818MR characterised by data acquisition along a specific k-space trajectory or by the temporal order of k-space coverage, e.g. centric or segmented coverage of k-space
    • G01R33/482MR characterised by data acquisition along a specific k-space trajectory or by the temporal order of k-space coverage, e.g. centric or segmented coverage of k-space using a Cartesian trajectory
    • G01R33/4822MR characterised by data acquisition along a specific k-space trajectory or by the temporal order of k-space coverage, e.g. centric or segmented coverage of k-space using a Cartesian trajectory in three dimensions
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01RMEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
    • G01R33/00Arrangements or instruments for measuring magnetic variables
    • G01R33/20Arrangements or instruments for measuring magnetic variables involving magnetic resonance
    • G01R33/44Arrangements or instruments for measuring magnetic variables involving magnetic resonance using nuclear magnetic resonance [NMR]
    • G01R33/48NMR imaging systems
    • G01R33/54Signal processing systems, e.g. using pulse sequences ; Generation or control of pulse sequences; Operator console
    • G01R33/543Control of the operation of the MR system, e.g. setting of acquisition parameters prior to or during MR data acquisition, dynamic shimming, use of one or more scout images for scan plane prescription
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01RMEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
    • G01R33/00Arrangements or instruments for measuring magnetic variables
    • G01R33/20Arrangements or instruments for measuring magnetic variables involving magnetic resonance
    • G01R33/44Arrangements or instruments for measuring magnetic variables involving magnetic resonance using nuclear magnetic resonance [NMR]
    • G01R33/48NMR imaging systems
    • G01R33/54Signal processing systems, e.g. using pulse sequences ; Generation or control of pulse sequences; Operator console
    • G01R33/56Image enhancement or correction, e.g. subtraction or averaging techniques, e.g. improvement of signal-to-noise ratio and resolution
    • G01R33/561Image 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
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01RMEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
    • G01R33/00Arrangements or instruments for measuring magnetic variables
    • G01R33/20Arrangements or instruments for measuring magnetic variables involving magnetic resonance
    • G01R33/44Arrangements or instruments for measuring magnetic variables involving magnetic resonance using nuclear magnetic resonance [NMR]
    • G01R33/48NMR imaging systems
    • G01R33/54Signal processing systems, e.g. using pulse sequences ; Generation or control of pulse sequences; Operator console
    • G01R33/56Image enhancement or correction, e.g. subtraction or averaging techniques, e.g. improvement of signal-to-noise ratio and resolution
    • G01R33/561Image 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
    • G01R33/5611Parallel magnetic resonance imaging, e.g. sensitivity encoding [SENSE], simultaneous acquisition of spatial harmonics [SMASH], unaliasing by Fourier encoding of the overlaps using the temporal dimension [UNFOLD], k-t-broad-use linear acquisition speed-up technique [k-t-BLAST], k-t-SENSE
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01RMEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
    • G01R33/00Arrangements or instruments for measuring magnetic variables
    • G01R33/20Arrangements or instruments for measuring magnetic variables involving magnetic resonance
    • G01R33/44Arrangements or instruments for measuring magnetic variables involving magnetic resonance using nuclear magnetic resonance [NMR]
    • G01R33/48NMR imaging systems
    • G01R33/54Signal processing systems, e.g. using pulse sequences ; Generation or control of pulse sequences; Operator console
    • G01R33/56Image enhancement or correction, e.g. subtraction or averaging techniques, e.g. improvement of signal-to-noise ratio and resolution
    • G01R33/561Image 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
    • G01R33/5615Echo train techniques involving acquiring plural, differently encoded, echo signals after one RF excitation, e.g. using gradient refocusing in echo planar imaging [EPI], RF refocusing in rapid acquisition with relaxation enhancement [RARE] or using both RF and gradient refocusing in gradient and spin echo imaging [GRASE]
    • G01R33/5617Echo train techniques involving acquiring plural, differently encoded, echo signals after one RF excitation, e.g. using gradient refocusing in echo planar imaging [EPI], RF refocusing in rapid acquisition with relaxation enhancement [RARE] or using both RF and gradient refocusing in gradient and spin echo imaging [GRASE] using RF refocusing, e.g. RARE
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01RMEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
    • G01R33/00Arrangements or instruments for measuring magnetic variables
    • G01R33/20Arrangements or instruments for measuring magnetic variables involving magnetic resonance
    • G01R33/44Arrangements or instruments for measuring magnetic variables involving magnetic resonance using nuclear magnetic resonance [NMR]
    • G01R33/48NMR imaging systems
    • G01R33/54Signal processing systems, e.g. using pulse sequences ; Generation or control of pulse sequences; Operator console
    • G01R33/56Image enhancement or correction, e.g. subtraction or averaging techniques, e.g. improvement of signal-to-noise ratio and resolution
    • G01R33/565Correction of image distortions, e.g. due to magnetic field inhomogeneities

Landscapes

  • Physics & Mathematics (AREA)
  • Health & Medical Sciences (AREA)
  • Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • High Energy & Nuclear Physics (AREA)
  • Engineering & Computer Science (AREA)
  • General Health & Medical Sciences (AREA)
  • Radiology & Medical Imaging (AREA)
  • General Physics & Mathematics (AREA)
  • Condensed Matter Physics & Semiconductors (AREA)
  • Signal Processing (AREA)
  • Surgery (AREA)
  • Molecular Biology (AREA)
  • Pathology (AREA)
  • Animal Behavior & Ethology (AREA)
  • Medical Informatics (AREA)
  • Public Health (AREA)
  • Veterinary Medicine (AREA)
  • Heart & Thoracic Surgery (AREA)
  • Biomedical Technology (AREA)
  • Biophysics (AREA)
  • Magnetic Resonance Imaging Apparatus (AREA)

Abstract

Die Erfindung betrifft ein Verfahren zum Erstellen eines MR 3D Bilddatensatzes, bei dem mit Spinechosequenzen zwei Rohdatensätze 40,50 unterabgetastet aufgenommen werden, wobei die Anregungspulse oder die Refokussierungspulse bei den beiden Rohdatensätzen eine entgegengesetzte Phasenlage haben. Diese beiden Rohdatensätze werden zu einem nicht unterabgetasteten kombinierten 3D Rohdatensatz 60 kombiniert und es wird eine Gewichtungsmatrix berechnet für die Verwendung bei der Berechnung der bei dem ersten Rohdatensatz und dem zweiten Rohdatensatz jeweils nicht aufgenommenen Rohdatenpunkte. Damit wird ein erster und zweiter vollständiger Rohdatensatz berechnet, die dann kombiniert werden.

Description

  • Die vorliegende Erfindung betrifft ein Verfahren zur Erstellung eines 3D MR Bilddatensatzes eines Untersuchungsobjekts und betrifft die zugehörige MR Anlage. Weiterhin sind ein Computerprogrammprodukt mit Programmmitteln und ein elektronisch lesbarer Datenträger vorgesehen.
  • In der MR Bildgebung ist eine schnelle 3D Spinechosequenz bekannt, bei der nach einem HF Anregungspuls ein Refokussierungspulszug mit bis zu mehreren 100 Refokussierungspulsen geschaltet wird mit teilweise reduzierten Refokussierungskippwinkeln. Ursprünglich wurde diese Bildgebungssequenz mit nichtselektiven HF Anregungspulse konzipiert, womit jedoch nur komplette Volumina aufgenommen werden können. Um derartige 3D Multi-Spinechosequenzen in Untersuchungsbereichen wie der Wirbelsäule, Hüfte oder dem Becken nutzen zu können, wurde ein selektiver Betriebsmodus dieser Sequenz eingeführt, bei dem ein selektiver HF Anregungspuls verwendet wird mit einer Folge von nichtselektiven Refokussierungspulsen.
  • Die Nutzung der nichtselektiven Refokussierungspulse nach dem Anregungspuls führt jedoch zu FID(Free Induction Decay)-Signalen innerhalb des Echozugs aus Bereichen außerhalb des selektiv angeregten Volumens. Dies kann zu Artefakten führen, welche die Bildgebung stören und die Befundung erschweren.
  • Es ist daher eine Aufgabe der vorliegenden Erfindung, diese schnellen 3D Spinecho basierten Bildgebungssequenzen so zu verbessern, dass diese Artefakte aufgrund des FID-Signals reduziert werden.
  • Diese Aufgabe wird durch die Merkmale der unabhängigen Ansprüche gelöst. In den abhängigen Ansprüchen sind bevorzugte Ausführungsformen der Erfindung beschrieben.
  • Gemäß einem ersten Aspekt der Erfindung wird ein Verfahren zum Erstellen eines 3D MR Bilddatensatzes eines Untersuchungsobjekts mit mehreren Empfangsspulen bereitgestellt. Das Verfahren weist folgende Schritte auf: Es wird zumindest ein erster HF Anregungspuls in das Untersuchungsobjekt eingestrahlt gefolgt von einem oder mehreren ersten nicht frequenzselektiven HF Refokussierungspulsen zur Erzeugung von zumindest einem ersten Spinnecho. Das oder die mehreren ersten Spinnechos werden in einem ersten Rohdatensatz im dreidimensionalen Rohdatenraum aufgenommen mit den mehreren Empfangsspulen, wobei hierbei der erste Rohdatensatz des Rohdatenraums nur teilweise derart mit Rohdaten gefüllt wird, dass nach dem Nyquist Theorem der erste Rohdatensatz nicht vollständig mit Rohdaten gefüllt ist. Weiterhin wird zumindest ein zweiter HF Anregungspuls in das Untersuchungsobjekt eingestrahlt gefolgt von einem oder mehreren zweiten nichtselektiven HF Refokussierungspulsen zur Erzeugung des oder der mehreren zweiten Spinnechos. Hierbei haben die zweiten HF Refokussierungspulse jeweils eine entgegengesetzte Phasenlage zu den ersten HF Refokussierungspulsen. Alternativ ist es möglich, dass der erste und zweite HF Anregungspuls eine entgegengesetzte Phasenlage haben. Das oder die mehreren zweiten Spinnechos werden in einem zweiten Rohdatensatz im dreidimensionalen Rohdatenraum mit den mehreren Empfangsspulen aufgenommen, wobei hierbei bei dem zweiten Rohdatensatz dieser nur teilweise derart mit Rohdaten gefüllt wird, dass nach dem Nyquist-Theorem der zweite Teilbereich nicht vollständig mit Rohdaten gefüllt ist, wobei der Rohdatenraum jedoch durch den ersten Rohdatensatz und zweiten Rohdatensatz insgesamt vollständig nach dem Nyquist-Theorem mit Rohdaten gefüllt ist. Der erste Rohdatensatz und der zweite Rohdatensatz werden zu einem kombinierten dreidimensionalen Rohdatensatz zusammengefasst, der vollständig nach dem Nyquist-Theorem mit Rohdaten gefüllt ist. Anschließend erfolgt die Berechnung einer Gewichtungsmatrix für eine Parallelbildgebung auf Grundlage des kombinierten 3D Rohdatensatzes für die Verwendung bei der Berechnung der bei dem ersten Rohdatensatz und dem zweiten Rohdatensatz jeweils nicht aufgenommenen Rohdatenpunkte. Weiterhin erfolgt eine Rekonstruktion von dem im ersten Rohdatensatz nicht aufgenommenen Rohdatenpunkten unter Berücksichtigung von der Gewichtungsmatrix und dem im ersten Rohdatensatz aufgenommenen Rohdaten. Diese Berechnungen für die beiden Rohdatensätze kann sowohl im Rohdatenraum als auch im Bildraum erfolgen, hier dann beispielsweise mit der SENSE Technik.
  • Es ist dann möglich, einen ersten vollständigen Rohdatensatz anhand der für den ersten Rohdatensatz rekonstruierten Rohdatenpunkte und der im ersten Rohdatensatz aufgenommenen Rohdaten zu berechnen. Weiterhin werden die im zweiten Rohdatensatz nicht aufgenommenen Rohdatenpunkte unter Berücksichtigung von der Gewichtungsmatrix und den im zweiten Rohdatensatz aufgenommenen Rohdatenpunkte rekonstruiert und es erfolgt die Berechnung eines zweiten vollständigen Rohdatensatzes anhand der für den zweiten Rohdatensatz rekonstruierten Rohdatenpunkte und der im zweiten Rohdatensatz aufgenommenen Rohdatenpunkte. Schließlich erfolgt die Addition des ersten vollständigen Rohdatensatzes und des zweiten vollständigen Rohdatensatzes zur Erstellung des 3D MR Bilddatensatzes.
  • Durch die Aufnahme der beiden Rohdatensätze mit der entgegengesetzten Phasenlage derart, dass sie zusammen den gesamten Rohdatenraum auffüllen, wird erstens die Aufnahmezeit verkürzt, da der gesamte Rohdatenraum nicht zweimal vollständig mit Rohdaten gefüllt werden muss. Zweitens ist es möglich, mit Hilfe des kombinierten 3D Rohdatensatzes die Gewichtungsmatrix zu berechnen, die anschließend für die Berechnung des ersten vollständigen Rohdatensatzes und des zweiten vollständigen Rohdatensatzes notwendig ist. Anschließend ist es möglich, diese beiden vollständigen Rohdatensätze zu kombinieren zur Erstellung des 3D MR Bilddatensatzes, bei dem dann die Artefakte aufgrund der FID-Signale reduziert bzw. vollständig unterdrückt sind. Durch die Verwendung der Refokussierungspulse im zweiten Rohdatensatz mit entgegengesetzter Phasenlage, d.h. mit einer um 180° gedrehten Phasenlage, werden die FID-Signale außerhalb des angeregten Untersuchungsobjekts destruktiv addiert und heben sich auf, so dass nur noch die Spinnechos den wesentlichen Signalanteil liefern. Weiterhin wird die Messzeit verringert, da jeweils der erste Rohdatensatz und der zweite Rohdatensatz unterabgetastet sind.
  • Der erste vollständige Rohdatensatz und der zweite vollständige Rohdatensatz können im Rohdatenraum oder im Bilddatenraum komplex addiert werden.
  • Weiterhin ist es möglich, dass der erste und der zweite Rohdatensatz derart aufgenommen werden, dass sie keine gemeinsamen Rohdatenpunkte haben. Der 3D Rohdatenraum wird vorzugsweise hälftig in dem ersten Rohdatensatz aufgenommen, wobei die andere Hälfte im zweiten Rohdatensatz aufgenommen wird, so dass insgesamt der 3D Rohdatenraum komplett aufgenommen wird, jedoch aber nur einmal.
  • Der erste vollständige Rohdatensatz und der zweite vollständige Rohdatensatz können mithilfe von Rekonstruktionsverfahren der parallelen Bildgebung, beispielsweise des GRAPPA-Verfahrens oder mithilfe des CAIPIRINHA-Verfahrens rekonstruiert werden. Die fehlenden Bereiche in den zwei Rohdatensätzen können somit und mithilfe der Gewichtungsmatrix rekonstruiert werden, wobei die Gewichtungsmatrix auf Grundlage des kombinierten 3D Rohdatensatzes berechnet werden kann.
  • Der erste HF Anregungspuls und der zweite HF Anregungspuls sind vorzugsweise frequenzselektive Anregungspulse, wobei das Verfahren jedoch auch mit nicht-frequenzselektiven HF Anregungspulsen verwendet werden kann.
  • Die Aufnahme des ersten Rohdatensatzes und des zweiten Rohdatensatzes kann getrennt voneinander nacheinander erfolgen oder in dem s.g. „interleaved“-Schema, bei dem Teile des zweiten Rohdatensatzes aufgenommen werden, bevor das Aufnehmen des ersten Rohdatensatzes vollständig abgeschlossen ist.
  • Vorzugsweise weisen der erste Rohdatensatz und der zweite Rohdatensatz jeweils die Hälfte des gesamten Rohdatenraums auf und füllen diese.
  • Weiterhin ist die zugehörige MR Anlage vorgesehen mit mehreren Empfangsspulen, einer HF Einheit und zumindest einer Bildsequenzsteuerung, die auch die mehreren Empfangsspulen und die HF Einheit steuert derart, dass die Aufnahme der Rohdatenraum wie oben erfolgt. Zusätzlich weist die MR Anlage eine Recheneinheit auf, die den 3D MR Bilddatensatz berechnet wie oben beschrieben.
  • Weiterhin ist ein Computerprogrammprodukt vorgesehen mit Programmmitteln und welches direkt in eine Speichereinheit einer MR Anlage ladbar ist, um die oben beschriebenen Schritte und die noch nachfolgend beschriebenen Schritte auszuführen.
  • Ebenso ist ein elektronisch lesbarer Datenträger vorgesehen mit darauf gespeicherten elektronisch lesbaren Steuerinformationen, wobei diese derart ausgestaltet sind, dass sie bei Verwendung des Datenträgers in der Steuereinheit der MR Anlage das oben beschriebene Verfahren durchführen.
  • Die oben beschriebenen Merkmale und die nachfolgend beschriebenen Merkmale können nicht nur in den entsprechenden explizit dargestellten Kombinationen verwendet werden, sondern auch in anderen Kombinationen, sofern es nicht explizit anderweitig erwähnt ist.
  • Die Erfindung wird unter Bezugnahme auf die beiliegenden Zeichnungen näher erläutert.
    • 1 zeigt schematisch eine MR Anlage, mit der erfindungsgemäß 3D Spinechosequenzen in kurzer Aufnahmezeit und ohne FID-Artefakte erzeugt werden können;
    • 2 zeigt schematisch ein Sequenzdiagramm der bei der Aufnahme der MR Daten verwendeten Bildgebungssequenz;
    • 3 zeigt schematisch das Aufnahmemuster, welches zur Erstellung des ersten Rohdatensatzes und des zweiten Rohdatensatzes verwendet werden kann;
    • 4 zeigt schematisch, wie aus der Kombination der beiden Rohdatensätze von 3 der kombinierte 3D Rohdatensatz erstellt wird, welcher zur Berechnung der Gewichtungsmatrix verwendet werden kann;
    • 5 zeigt schematisch die bei Erstellung des 3D MR Datensatzes mit vermindertem FID-Artefakt die verschiedenen Datensätze aufgenommen und welche weiteren Datensätze berechnet werden;
    • 6 zeigt schematisch ein Flussdiagramm mit den Schritten, die durchgeführt werden zur Erstellung des 3D MR Bilddatensatzes mit reduziertem FID-Artefakt.
  • Nachfolgend wird die vorliegende Erfindung anhand bevorzugter Ausführungsformen unter Bezugnahme auf die beiliegenden Zeichnungen näher erläutert. In den Figuren beschreiben gleiche Bezugszeichen gleiche oder ähnliche Elemente. Weiterhin sind die Figuren schematische Darstellungen verschiedener Ausführungsformen, und die in den Figuren dargestellten Elemente sind nicht notwendigerweise maßstabsgetreu dargestellt. Vielmehr sind die in den Figuren dargestellten Elemente derart wiedergegeben, dass ihre Funktion und ihr Zweck für den Fachmann verständlich werden. Die in den Figuren dargestellten Verbindungen zwischen funktionellen Einheiten oder sonstigen Elementen können auch als indirekte Verbindung implementiert werden, wobei eine Verbindung drahtlos oder drahtgebunden sein kann. Funktionelle Einheiten können als Hardware, Software oder als Kombination aus Hardware und Software implementiert werden.
  • Unter Bezugnahme auf 1 wird eine MR Anlage 9 erläutert, mit der, wie nachfolgend erläutert werden wird, 3D Spinecho basierte MR Bilder mit reduziertem FID-Artefakt erzeugt werden können. Die MR Anlage 9 weist einen Magneten 10 zur Erzeugung des Polarisationsfeldes B0 auf, wobei eine auf einer Liege 12 angeordnete Untersuchungsperson in den Magneten gefahren wird, um dort ortskodierte Magnetresonanzsignale aus der Untersuchungsperson 13 mithilfe von mehreren Empfangsspulen 11 aufzunehmen. Die Erfindung findet bei der s.g. parallelen Bildgebung Anwendung, bei der die MR Signale gleichzeitig mit den mehreren Empfangsspulen 11 aufgenommen werden. Durch Einstrahlen von Hochfrequenzpulsen und Schalten von Magnetfeldgradienten kann die durch das Polarisationsfeld B0 erzeugte Magnetisierung aus der Gleichgewichtslage ausgelenkt und ortskodiert werden, und die sich ergebende Magnetisierung wird von den Empfangsspulen 11 detektiert. Wie durch Einstrahlen der HF Pulse und durch Schalten von Magnetfeldgradienten in verschiedenen Kombinationen und Reihenfolgen MR Bilder erzeugt werden können, ist dem Fachmann grundsätzlich bekannt und wird hier nicht näher erläutert.
  • Die MR Anlage weist eine Steuereinheit 20 auf zur Steuerung der MR Anlage 9. Die Steuereinheit 20 weist eine HF Steuereinheit 14 auf für die Steuerung und Generierung der HF Pulse zur Auslenkung der Magnetisierung. Eine Gradientensteuereinheit 15 ist vorgesehen zur Steuerung und Schaltung der notwendigen Magnetfeldgradienten. Eine Bildsequenzsteuerung 16 steuert die Abfolge der Magnetfeldgradienten, der Signaldetektion und der HF Pulse und damit indirekt die Gradientensteuereinheit 15, die Empfangsspulen 11 und die HF Steuereinheit 14. Über eine Eingabeeinheit 17 kann eine Bedienperson die MR Anlage steuern und auf einer Anzeigeeinheit 18 können MR Bilder und sonstige zur Steuerung notwendigen Informationen dargestellt werden. Eine Recheneinheit 19 mit mindestens einer Prozessoreinheit ist vorgesehen zur Steuerung der verschiedenen Einheiten in der Steuereinheit 20. Weiterhin ist eine Speichereinheit 21 vorgesehen, in der beispielsweise Programmmodule bzw. Programmmittel abgespeichert werden können, die, wenn sie von der Recheneinheit 19 bzw. der darin ausgebildeten Prozessoreinheit ausgeführt werden, den Ablauf der MR Anlage steuern können. Wie nachfolgend erläutert wird, sind insbesondere die Bildsequenzsteuerung 16 und die Recheneinheit 19 derart ausgebildet, dass ein 3D Rohdatenraum auf spezielle Weise mit Spinnechos gefüllt wird, um einen 3D MR Bilddatensatz zu erzeugen, der mit einer verringerten Aufnahmezeit gegenüber dem Stand der Technik die Vermeidung von FID-Artefakten vermeidet.
  • In 2 ist schematisch die 3D MR Bildgebungssequenz dargestellt. Nach einem frequenzselektiven HF Anregungspuls 22 erfolgt ein erster HF Refokussierungspuls 23, gefolgt von weiteren nicht selektiven HF Refokussierungspulsen 24-26, die einen geringeren Fokussierungswinkel als 180° aufweisen können zur Reduktion der in die Untersuchungsperson eingestrahlten Energie. In Schichtselektionsrichtung Gz erfolgt gleichzeitig zu dem Anregungspuls die Schaltung des Schichtselektionsgradienten 27 zur Anregung des gewünschten Bereiches im Untersuchungsobjekt. Weiterhin erfolgen die für die Phasenkodierung notwendigen Gradientenschaltungen 28 und 29, ebenso in Phasenkodierrichtung die Gradientenschaltungen 30 und 31. In Ausleserichtung erfolgt während der Signalauslese 32,33 die Schaltung der Auslesegradienten 34-36, wobei das erste Echo nicht ausgelesen werden kann.
  • Unter zusätzlicher Bezugnahme auf 3 wird nun erläutert, wie mit einer derartigen Bildgebungssequenz der Rohdatenraum aufgenommen wird. In einem ersten Schritt wird nur, wie links in 3 zu erkennen ist, die erste Hälfte des Rohdatenraums aufgenommen, wobei mit Bezugszeichen 45 die nicht aufgenommenen Rohdatenlinien dargestellt sind, die in die Zeichenebene hinein in die kx -Richtung verlaufen, da es sich um einen 3D Datensatz handelt. Weiterhin sind die Rohdaten 46 dargestellt, die mit einem ersten Phasenzyklus aufgenommen werden, was bedeutet, dass die Refokussierungspulse 23-26 eine erste Phase relativ zu dem Anregungspuls haben. Beispielsweise können, wenn der HF Anregungspuls in die x-Richtung eingestrahlt wird, alle Refokussierungspulse in die y-Richtung eingestrahlt werden. Es handelt sich somit um ein Abtastschema, das ähnlich zu dem IPAT-Muster ist, jedoch ohne Aufnahme von Referenzlinien. Das führt zu einem ersten Rohdatensatz 40.
  • Weiterhin wird in einem zweiten Schritt ein zweiter Rohdatensatz 50 aufgenommen, der in 3 rechts dargestellt ist, wobei mit 55 wiederum die nicht aufgenommenen Rohdatenlinien dargestellt sind, während mit 56 die Daten dargestellt sind, die bei dem zweiten Rohdatensatz 50 aufgenommen werden, wobei sich die Phase der Refokussierungspulse von der Phase der Refokussierungspulse bei der Aufnahme der Rohdatenpunkte 46 um 180° unterscheidet. Bezugnehmend auf 2 bedeutet dies, dass einmal die Refokussierungspulse mit einer ersten Phasenlage aufgenommen werden, und ein zweites Mal für den zweiten Rohdatensatz 50 mit einer zweiten entgegengesetzten Phasenlage, d.h. um 180° versetzten Phasenlage. Wenn beispielsweise die Refokussierungspulse 23-26 für den ersten Rohdatensatz entlang der +y-Richtung eingestrahlt werden, so ist es möglich, sie bei der Einstrahlung für den zweiten Rohdatensatz entlang der -y-Richtung einzustrahlen. Alternativ ist es möglich, eine entgegengesetzte Phasenlage bei den Anregungspulsen zu haben.
  • Die beiden Rohdatensätze 40, 50 werden nun so aufgenommen, dass letztendlich der 3D Rohdatenraum komplett aufgenommen wird, aber nur einmal während jeder der beiden Rohdatensätze von 3 unterabgetastet ist.
  • Wie in 4 zu erkennen ist, ist es dann möglich, die beiden Rohdatensätze 40 und 50 zu kombinieren, so dass ein kombinierter 3D Rohdatensatz erzeugt wird, wie er in 4 mit Bezugszeichen 60 dargestellt ist. Hierbei werden die aufgenommenen Rohdaten der beiden Rohdatensätze 40 und 50 kombiniert. Bei Verwendung von entgegengesetzter Phasenlage bei den Anregungspulsen werden die beiden Rohdatensätze dann in der Kombination subtrahiert. Aus dem kombinierten 3D Rohdatensatz 60 ist es dann möglich, aus einem zentralen Bereich 61 des kombinierten Rohdatensatzes eine Gewichtungsmatrix zu berechnen, beispielsweise den s.g. GRAPPA-Kernel, mit dem dann die fehlenden Bereiche in den beiden Rohdatensätzen 40 und 50 rekonstruiert werden können, nämlich die Rohdatenpunkte mit Bezugszeichen 45 und 55. Damit werden jeweils ein erster vollständiger Rohdatensatz berechnet anhand der Gewichtungsmatrix bzw. des Kernels und ein zweiter vollständiger Rohdatensatz anhand der Gewichtungsmatrix. Diese beiden vollständigen Rohdatensätze können dann addiert werden, wobei es sich um eine komplexe Addition handelt, um die FID-Artefakte zu vermeiden. Ebenso ist es möglich, dass die beiden Rohdatenräume erst komplex addiert und dann rekonstruiert werden.
  • Da die komplette Aufnahme des Rohdatenraums normalerweise mehrere Echozüge notwendig macht, so müssen die beiden Rohdatensätze nicht zeitlichen nacheinander aufgenommen werden, sie können auch in einem s.g. „interleaved“-Schema aufgenommen werden.
  • Option A (konventionelles Schema):
    R1,+, R2,+, ... Rn/2,+, Rn/2+1,-, Rn/2+2,- ... Rn,-
  • Option B (interleaved Schema):
    R1,+, R2,-, R3,+, R4,- ... Rn-1,+, Rn,-
  • Hierbei ist mit R die Anzahl der Echozüge bezeichnet, und + bzw. - bezeichnet die Phase des Refokussierungspulses, wobei N die Anzahl der Wiederholungen ist.
  • 5 zeigt schematisch noch einmal die bei den verschiedenen Schritten aufgenommenen bzw. berechneten Datensätze. Mit den ersten HF Anregungspulsen und den ersten nicht frequenzselektiven Refokussierungspulsen wird der erste Rohdatensatz 40 aufgenommen, wobei mit den zweiten Spinnechos und den entgegengesetzten Refokussierungspulsen der zweite Rohdatensatz 50 aufgenommen wird. Diese können dann zu einem kombinierten 3D Rohdatensatz 60 kombiniert werden, wobei mithilfe des zentralen Bereichs 61 die Gewichtungsmatrix berechnet werden kann. Mit dieser Gewichtungsmatrix kann dann mithilfe der aufgenommenen Rohdaten des Rohdatensatzes 40 ein erster vollständiger Rohdatensatz 70 berechnet werden, wobei ebenso anhand der im zweiten Rohdatensatz gemessenen Rohdatenpunkte und der Gewichtungsmatrix der zweite vollständige Rohdatensatz 80 berechnet werden kann. Durch Kombination dieser beiden vollständigen Rohdatensätze ist es möglich, den 3D Rohdatensatz 100 zu erzeugen, wobei die Addition eine komplexe Addition ist.
  • Bei dem in 4 verwendeten Schema wurden abwechselnd Rohdatenpunkte bzw. Rohdatenlinien in kz bzw. ky Richtung aufgenommen und nicht aufgenommen. Das Verfahren kann jedoch beliebig in kz oder ky Richtung angewandt werden, d.h. die beiden Rohdatenräume 40 und 50 können beliebig ausgebildet sein, solange sie nicht vollständig mit Rohdaten gefüllt sind und sich die beiden Rohdatensätze 40 und 50 nicht in den Rohdatenpunkten überlappen.
  • Für die Berechnung des ersten vollständigen Rohdatensatzes und des zweiten vollständigen Rohdatensatzes können neben dem GRAPPA-Verfahren auch andere Verfahren der parallelen Bildgebung wie CAIPIRINHA verwendet werden.
  • In 6 werden die Schritte noch einmal zusammengefasst. In einem Schritt S61 erfolgt die Aufnahme des ersten Rohdatensatzes 40 mit der Einstrahlung der ersten HF Anregungspulse und der ersten nicht frequenzselektiven HF Refokussierungspulse. Dieser erste Rohdatensatz ist nicht vollständig mit Rohdaten gefüllt und ist somit unterabgetastet nach dem Nyquist-Theorem. Ebenso wird in Schritt S62 der zweite Rohdatensatz 50 aufgenommen, wie es in Zusammenhang mit 2 und 3 erläutert wurde. Wie oben erläutert, müssen die Schritte S61 und S62 nicht nacheinander durchgeführt werden, sondern können auch parallel durchgeführt werden im s.g. „interleaved“-Verfahren, bei dem Teile des ersten Rohdatensatzes 40 und Teile des zweiten Rohdatensatzes 50 abwechselnd aufgenommen werden. Im Schritt S63 erfolgt dann die Kombination der aufgenommenen Rohdaten zu dem kombinierten 3D Rohdatensatz 60. Anschließend können im Schritt S64 die Gewichtungsmatrix bzw. der Kernel berechnet werden auf Grundlage des kombinierten 3D Rohdatensatzes. Hiermit ist es schließlich möglich, in Schritt S65 jeweils den ersten vollständigen Rohdatensatz 70 und den zweiten vollständigen Rohdatensatz 80 zu berechnen, wobei durch Addition dieser vollständigen Rohdatensätze schließlich der 3D MR Bilddatensatz 100 im Schritt S66 berechnet werden kann. Die mit diesem 3D MR Bilddatensatz erzeugten MR Bilder weisen keine oder nur stark unterdrückte FID-Artefakte auf, wobei insgesamt gegenüber dem Stand der Technik die Messzeit um den Faktor 2 zumindest reduziert werden kann.

Claims (19)

  1. Verfahren zum Erstellen eines MR 3D Bilddatensatzes (100) eines Untersuchungsobjekts (13) mit mehreren Empfangsspulen (11), mit den folgenden Schritten: - Einstrahlen von zumindest einem ersten HF Anregungspuls in das Untersuchungsobjekt, - Einstrahlen von zumindest einem ersten nicht frequenzselektiven HF Refokussierungspuls zur Erzeugung von zumindest einem ersten Spinecho, die - Aufnehmen des zumindest einen ersten Spinechos in einem ersten Rohdatensatz (40) im dreidimensionalen Rohdatenraum mit den mehreren Empfangsspulen (11), wobei hierbei der erste Rohdatensatz des Rohdatenraums nur teilweise derart mit Rohdaten gefüllt wird, dass nach dem Nyquist-Theorem der erste Rohdatensatz nicht vollständig mit Rohdaten gefüllt ist, - Einstrahlen von zumindest einem zweiten HF Anregungspuls in das Untersuchungsobjekt, - Einstrahlen von zumindest einem zweiten nicht selektiven HF Refokussierungspuls zur Erzeugung von zumindest einem zweiten Spinecho, wobei der zumindest eine zweite HF Refokussierungspuls eine entgegengesetzte Phasenlage zu dem zumindest einen ersten HF Refokussierungspuls hat, oder der zumindest eine zweite HF Anregungspuls eine entgegengesetzte Phasenlage zu dem zumindest einen ersten HF Anregungspuls hat, - Aufnehmen des zumindest einen zweiten Spinechos in einem zweiten Rohdatensatzes (50) im dreidimensionalen Rohdatenraum mit den mehreren Empfangsspulen, wobei hierbei bei dem zweiten Rohdatensatz (50) der Rohdatenraum nur teilweise derart mit Rohdaten gefüllt wird, dass nach dem Nyquist-Theorem der zweite Rohdatensatz nicht vollständig mit Rohdaten gefüllt ist, und der Rohdatenraum durch den ersten Rohdatensatz und zweiten Rohdatensatz insgesamt vollständig nach dem Nyquist-Theorem mit Rohdaten gefüllt ist, - Zusammenfassen des ersten Rohdatensatzes (40) und des zweiten Rohdatensatzes (50) zu einem kombinierten 3D Rohdatensatz (60), der vollständig nach dem Nyquist-Theorem mit Rohdaten gefüllt ist, - Berechnen einer Gewichtungsmatrix auf Grundlage des kombinierten 3D Rohdatensatz (60) für die Verwendung bei der Berechnung der bei dem ersten Rohdatensatz und dem zweiten Rohdatensatz jeweils nicht aufgenommenen Rohdatenpunkte, - Rekonstruieren von den im ersten Rohdatensatz (40) nicht aufgenommenen Rohdatenpunkten unter Berücksichtigung von der Gewichtungsmatrix und den im ersten Rohdatensatz (40) aufgenommenen Rohdaten, - Berechnen eines ersten vollständigen Rohdatensatzes (70) anhand der für den ersten Rohdatensatz rekonstruierten Rohdatenpunkte und der im ersten Rohdatensatz (40) aufgenommenen Rohdaten, - Rekonstruieren von den im zweiten Rohdatensatz (50) nicht aufgenommenen Rohdatenpunkten unter Berücksichtigung von der Gewichtungsmatrix und den im zweiten Rohdatensatz aufgenommenen Rohdaten, - Berechnen eines zweiten vollständigen Rohdatensatzes (80) anhand der für den zweiten Rohdatensatz rekonstruierten Rohdatenpunkte und der im zweiten Rohdatensatz aufgenommenen Rohdaten, - Addieren des ersten vollständigen Rohdatensatzes und des zweiten vollständigen Rohdatensatzes zur Erstellung des MR 3D Bilddatensatzes (100).
  2. Verfahren nach Anspruch 1, dadurch gekennzeichnet, dass der erste vollständige Rohdatensatz (70) und der zweite vollständige Rohdatensatz (80) im Rohdatenraum komplex addiert werden.
  3. Verfahren nach Anspruch 1, dadurch gekennzeichnet, dass der erste vollständige Rohdatensatz (70) und der zweite vollständige Rohdatensatz (80) im Bildraum unter Berücksichtigung der Phaseninformation komplex addiert werden.
  4. Verfahren nach einem der vorstehenden Ansprüche, dadurch gekennzeichnet, dass der erste und zweite Rohdatensatz (40, 50) derart aufgenommen werden, dass sie keine gemeinsamen Rohdatenpunkte haben.
  5. Verfahren nach einem der vorstehenden Ansprüche, dadurch gekennzeichnet, dass der erste vollständige Rohdatensatz (70) und der zweite vollständige Rohdatensatz (80) mit Hilfe von einem Rekonstruktionsverfahren der parallelen Bildgebung berechnet werden.
  6. Verfahren nach einem der vorstehenden Ansprüche, dadurch gekennzeichnet, dass der zumindest eine erste HF Anregungspuls und zumindest eine zweite HF Anregungspuls frequenzselektive Anregungspulse sind.
  7. Verfahren nach einem der vorstehenden Ansprüche, dadurch gekennzeichnet, dass erst der erste Rohdatensatz (40) komplett aufgenommen wird, bevor der zweite Rohdatensatz (50) komplett aufgenommen wird.
  8. Verfahren nach einem der Ansprüche 1 bis 6, dadurch gekennzeichnet, dass der zweite Rohdatensatz (50) zumindest teilweise aufgenommen wird, bevor das Aufnehmen des ersten Rohdatensatzes vollständig abgeschlossen ist.
  9. Verfahren nach einem der vorstehenden Ansprüche, dadurch gekennzeichnet, dass der erste Rohdatensatz (40) und der zweite Rohdatensatz (50) so aufgenommen werden, dass sie im Wesentlichen gleich groß sind.
  10. MR Anlage ausgebildet zur Erstellung eines MR 3D Bilddatensatzes eines Untersuchungsobjekts, welche aufweist: - mehrere Empfangsspulen (11), - eine HF Einheit (14) zum Einstrahlen von HF Pulsen, - eine Bildsequenzsteuerung (16), die ausgebildet ist zur Steuerung von zumindest den mehreren Empfangsspulen (11) und der HF Einheit (14) derart, dass die folgenden Schritte ausgeführt werden: - Einstrahlen von zumindest einem ersten HF Anregungspuls in das Untersuchungsobjekt, - Einstrahlen von zumindest einem ersten nicht frequenzselektiven HF Refokussierungspuls zur Erzeugung von zumindest einem ersten Spinecho, - Aufnehmen des zumindest einen ersten Spinechos in einem ersten Rohdatensatz (40) im dreidimensionalen Rohdatenraum mit den mehreren Empfangsspulen (11), wobei hierbei der erste Rohdatensatz (40) des Rohdatenraums nur teilweise derart mit Rohdaten gefüllt wird, dass nach dem Nyquist-Theorem der erste Rohdatensatz nicht vollständig mit Rohdaten gefüllt ist, - Einstrahlen von zumindest einem zweiten HF Anregungspuls in das Untersuchungsobjekt, - Einstrahlen von zumindest einem zweiten nicht frequenzselektiven HF Refokussierungspuls zur Erzeugung von zumindest einem zweiten Spinecho, wobei der zumindest eine zweite HF Refokussierungspuls eine entgegengesetzte Phasenlage zu dem zumindest einen ersten HF Refokussierungspuls hat, oder der zumindest eine zweite HF Anregungspuls eine entgegengesetzte Phasenlage zu dem zumindest einen ersten HF Anregungspuls hat - Aufnehmen des zumindest einen zweiten Spinechos in einem zweiten Rohdatensatzes (50) im dreidimensionalen Rohdatenraum mit den mehreren Empfangsspulen (11), wobei hierbei der zweite Rohdatensatz des Rohdatenraums nur teilweise derart mit Rohdaten gefüllt wird, dass nach dem Nyquist-Theorem der zweite Teilbereich nicht vollständig mit Rohdaten gefüllt ist und der Rohdatenraum durch den ersten Rohdatensatz (40) und zweiten Rohdatensatz (50) insgesamt vollständig nach dem Nyquist-Theorem mit Rohdaten gefüllt ist, - eine Recheneinheit (19) die ausgebildet ist, - den ersten Rohdatensatz (40) und den zweiten Rohdatensatz (50) zu einem kombinierten 3D Rohdatensatz (60) zusammenzufassen, der vollständig nach dem Nyquist-Theorem mit Rohdaten gefüllt ist, - eine Gewichtungsmatrix auf Grundlage des kombinierten 3D Rohdatensatz (60) zu berechnen für die Verwendung bei der Berechnung der bei dem ersten Rohdatensatz und dem zweiten Rohdatensatz jeweils nicht aufgenommenen Rohdatenpunkte, - die im ersten Rohdatensatz nicht aufgenommenen Rohdatenpunkte zu rekonstruieren unter Berücksichtigung von der Gewichtungsmatrix und den im ersten Rohdatensatz aufgenommenen Rohdaten, - einen ersten vollständigen Rohdatensatz (70) anhand der für den ersten Rohdatensatz rekonstruierten Rohdatenpunkte und der im ersten Rohdatensatz (40) aufgenommenen Rohdaten zu berechnen, - die im zweiten Rohdatensatz nicht aufgenommenen Rohdatenpunkte zu rekonstruieren unter Berücksichtigung von der Gewichtungsmatrix und den im zweiten Rohdatensatz (50) aufgenommenen Rohdaten, - einen zweiten vollständigen Rohdatensatz (80) anhand der für den zweiten Rohdatensatz rekonstruierten Rohdatenpunkte und der im zweiten Rohdatensatz aufgenommenen Rohdaten zu berechnen, - den ersten vollständigen Rohdatensatz (70) und den zweiten vollständigen Rohdatensatz (80) zu addieren zur Erstellung des MR 3D Bilddatensatzes (100).
  11. MR Anlage nach Anspruch 10, dadurch gekennzeichnet, dass die Recheneinheit ausgebildet ist, den ersten vollständigen Rohdatensatz (70) und den zweiten vollständigen Rohdatensatz (80) im Rohdatenraum komplex zu addieren.
  12. MR Anlage nach Anspruch 10 oder 11, dadurch gekennzeichnet, dass die Recheneinheit (19) ausgebildet ist, den ersten vollständigen Rohdatensatz (70) und den zweiten vollständigen Rohdatensatz (80) im Bildraum unter Berücksichtigung der Phaseninformation komplex zu addieren.
  13. MR Anlage nach einem der Ansprüche 10 bis 12, dadurch gekennzeichnet, dass die Bildsequenzsteuerung (16) ausgebildet ist, den ersten und zweiten Rohdatensatz (40, 50) jeweils derart aufzunehmen, dass sie keine gemeinsamen Rohdatenpunkte haben.
  14. MR Anlage nach einem der Ansprüche 10 bis 13, dadurch gekennzeichnet, dass die Recheneinheit (19) ausgebildet ist, den ersten vollständigen Rohdatensatz (70) und den zweiten vollständigen Rohdatensatz (80) mit Hilfe eines GRAPPA-Verfahrens zu berechnen.
  15. MR Anlage nach einem der Ansprüche 10 bis 14, dadurch gekennzeichnet, dass die HF Einheit (14) ausgebildet ist, den zumindest einen ersten HF Anregungspuls und den zumindest einen zweiten HF Anregungspuls als frequenzselektive Anregungspulse auszubilden.
  16. MR Anlage nach einem der Ansprüche 10 bis 15, dadurch gekennzeichnet, dass die Bildsequenzsteuerung (16) ausgebildet ist, erst den ersten Rohdatensatz (40) komplett aufzunehmen, bevor der zweite Rohdatensatz (50) komplett aufgenommen wird.
  17. MR Anlage nach einem der Ansprüche 10 bis 15, dadurch gekennzeichnet, dass die Bildsequenzsteuerung (16) ausgebildet ist, den zweiten Rohdatensatz (50) zumindest teilweise aufzunehmen, bevor das Aufnehmen des ersten Rohdatensatzes (40) vollständig abgeschlossen ist.
  18. Computerprogrammprodukt, welches Programmmittel umfasst und direkt in eine Speichereinheit einer Steuereinheit einer MR Anlage ladbar ist, um alle Schritte des Verfahrens nach einem der Ansprüche 1-9 auszuführen, wenn die Programmmittel in der Steuereinheit ausgeführt werden.
  19. Elektronisch lesbarer Datenträger mit darauf gespeicherten elektronisch lesbaren Steuerinformationen, welche derart ausgestaltet sind, dass sie bei Verwendung des Datenträgers in einer Steuereinheit einer MR Anlage das Verfahren nach einem der Ansprüche 1-9 durchführen.
DE102018200900.7A 2018-01-22 2018-01-22 Verfahren und Magnetresonanzanlage zur Artefaktvermeidung bei schnellen 3D Spinechosequenzen Active DE102018200900B4 (de)

Priority Applications (3)

Application Number Priority Date Filing Date Title
DE102018200900.7A DE102018200900B4 (de) 2018-01-22 2018-01-22 Verfahren und Magnetresonanzanlage zur Artefaktvermeidung bei schnellen 3D Spinechosequenzen
US16/251,424 US10605879B2 (en) 2018-01-22 2019-01-18 Magnetic resonance method and apparatus for artifact prevention in fast 3D spin echo sequences
CN201910053157.3A CN110068782B (zh) 2018-01-22 2019-01-21 快速三维自旋回波序列中的伪影避免

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
DE102018200900.7A DE102018200900B4 (de) 2018-01-22 2018-01-22 Verfahren und Magnetresonanzanlage zur Artefaktvermeidung bei schnellen 3D Spinechosequenzen

Publications (2)

Publication Number Publication Date
DE102018200900A1 true DE102018200900A1 (de) 2019-07-25
DE102018200900B4 DE102018200900B4 (de) 2023-02-02

Family

ID=67145203

Family Applications (1)

Application Number Title Priority Date Filing Date
DE102018200900.7A Active DE102018200900B4 (de) 2018-01-22 2018-01-22 Verfahren und Magnetresonanzanlage zur Artefaktvermeidung bei schnellen 3D Spinechosequenzen

Country Status (3)

Country Link
US (1) US10605879B2 (de)
CN (1) CN110068782B (de)
DE (1) DE102018200900B4 (de)

Family Cites Families (11)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
DE102004017852B4 (de) * 2004-04-13 2008-11-27 Siemens Ag Bewegungskorrigiertes Multi-Shot-Verfahren zur diffusionsgewichteten Bildgebung in der Magnetresonanztomographie
DE102006042998B4 (de) * 2006-09-13 2008-07-03 Siemens Ag Messsequenz für die dreidimensionale MR-Bildgebung sowie MR-Gerät
US7898254B2 (en) * 2008-02-19 2011-03-01 Advanced Mri Technologies, Llc Arterial spin labeled, segmented, interleaved 3D GRASE MRI
US8981776B2 (en) * 2011-04-22 2015-03-17 The General Hospital Corporation Method for magnetic resonance imaging with controlled aliasing
DE102011082010B4 (de) * 2011-09-01 2013-03-14 Siemens Aktiengesellschaft Rephasierung von Spinsystemen in einer ersten und zweiten Schicht in Schicht-Multiplexing-Messsequenzen zur Magnetresonanzbildgebung
US9678190B2 (en) * 2012-04-06 2017-06-13 General Electric Company System and method for generating MR phase contrast images near metal
DE102012216353B3 (de) * 2012-09-14 2014-02-27 Siemens Aktiengesellschaft Ansteuerung eines Magnetresonanzsystems zur selektiven Anregung von Teilbereichen
DE102013205930B4 (de) * 2013-04-04 2014-11-20 Siemens Aktiengesellschaft Bestimmung einer Resonanzfrequenzabweichung bei räumlich verzerrter Schichtanregung
CN104714199B (zh) * 2013-12-17 2018-04-24 西门子(深圳)磁共振有限公司 一种磁共振成像方法和装置
WO2017207700A1 (en) 2016-06-02 2017-12-07 Koninklijke Philips N.V. Dixon-type water/fat separation mr imaging
DE102016212632A1 (de) * 2016-07-12 2018-01-18 Siemens Healthcare Gmbh Reduzierung von Artefakten in der Magnetresonanztechnik

Also Published As

Publication number Publication date
CN110068782A (zh) 2019-07-30
US20190227139A1 (en) 2019-07-25
DE102018200900B4 (de) 2023-02-02
US10605879B2 (en) 2020-03-31
CN110068782B (zh) 2021-04-09

Similar Documents

Publication Publication Date Title
DE102011080254B4 (de) Bewegungskorrigierte Multi-Shot-Verfahren zur diffusionsgewichteten MR-Bildgebung mit zusätzlichen Referenz-Rohdaten und entsprechende Vorrichtungen
DE102007015040B3 (de) Verfahren zur Beschleunigung der Datenaufnahme bei dynamischen Messungen der Magnetresonanz (MRT) mittels paralleler Bildgebung
DE102009025753A1 (de) System und Verfahren zur Nutzung der Parallel Imaging Technik mit der Compressed Sensing Technik
DE102010032080B4 (de) Getriggerte Magnetresonanzbildgebung auf der Grundlage einer partiellen parallelen Akquisition (PPA)
DE102007044463A1 (de) Verfahren zur Bestimmung der räumlichen Verteilung von Magnetresonanzsignalen durch mehrdimensionale HF-Anregungspulse
DE102014206395B4 (de) Aufnahme und iterative Rekonstruktion einer Parameterkarte eines Zielbereichs
DE112015001951T5 (de) System und Verfahren zur Magnetresonanz-Bildgebung mit reduziertem Sichtfeld
DE60122818T2 (de) Selbstkorrektur von Bewegungsartefakten in MR-Projektionsbildern
DE102011005084B3 (de) Reduktion von Artefakten in der diffusionsgewichteten Bildgebung
DE102013220012B4 (de) Bildrekonstruktion mit mehreren virtuellen Spulen
DE102006055933A1 (de) Verfahren zur Ermittlung einer Bewegung bei der Aufzeichnung von MR-Messdaten und Magnet-Resonanz-Gerät hierzu
EP3557276A1 (de) Verfahren und steuereinrichtung zur erzeugung von magnetresonanzaufnahmen innerhalb und ausserhalb des b0-homogenitätsvolumens
EP3376246A1 (de) Beschleunigtes erzeugen einer serie von magnetresonanzbildern mit simultaner magnetresonanz-mehrschichtbildgebung
DE102009003566B4 (de) Erfassung von MRI-Daten für Impulssequenzen mit mehreren Phasenkodierungsrichtungen und periodischer Signalmodulation
DE102016224175A1 (de) EPI mit Zickzack-Trajektorie und paralleler Bildgebung
DE112019000571T5 (de) MR-Bildgebung mithilfe einer Stack-Of-Stars-Erfassung mit intrinsischer Bewegungskorrektur
DE102014209803B4 (de) 4D-Geschwindigkeitsprofil-Bildaufnahme mit einer Magnetresonanzanlage
DE102020209382A1 (de) Verfahren zur Aufnahme von Messdaten mittels einer Magnetresonanzanlage mit einer Korrektur der verwendeten k-Raumtrajektorien
DE102010041659B4 (de) Erzeugung eines optimierten MR-Bildes eines Untersuchungsobjekts durch Einstrahlen einer mindestens zwei HF-Pulse umfassenden Pulsfolge
DE102010003552B4 (de) Verfahren zur Homogenisierung der Auflösung in Magnetresonanz-Tomographie-Aufnahmen unter Verwendung von nicht-linearen Kodierfeldern
DE102014200562B4 (de) MRI-Sequenz mit einer Anzahl von schichtselektiven Inversionspulsen zur Präparation der Magnetisierung sowie einer hierzu disjunkten Anzahl von Sättigungspulsen
DE102017203082B4 (de) Parallele Bildgebung mit verbessertem Referenzdatensatz zur Bestimmung der Gewichtungsmatrix
DE102015209838B4 (de) Verfahren zu einer Korrektur von Magnetresonanz-Messdaten
DE102018200900A1 (de) Artefaktvermeidung bei schnellen 3D Spinechosequenzen
DE102013220010B4 (de) Bildkonstruktion mit mehreren virtuellen Spulen

Legal Events

Date Code Title Description
R012 Request for examination validly filed
R016 Response to examination communication
R018 Grant decision by examination section/examining division
R020 Patent grant now final
R081 Change of applicant/patentee

Owner name: SIEMENS HEALTHINEERS AG, DE

Free format text: FORMER OWNER: SIEMENS HEALTHCARE GMBH, MUENCHEN, DE