DE102010012599A1

DE102010012599A1 - Verfahren zur Erstellung eines Bilddatensatzes mittels einer radialen Abtastung mit Hilfe einer Magnetresonanzanlage sowie entsprechende Magnetresonanzanlage

Info

Publication number: DE102010012599A1
Application number: DE102010012599A
Authority: DE
Inventors: Dr. Block Kai Tobias
Original assignee: Siemens AG
Current assignee: Siemens Healthcare GmbH
Priority date: 2010-03-24
Filing date: 2010-03-24
Publication date: 2011-09-29
Anticipated expiration: 2030-03-25
Also published as: DE102010012599B4; US20110234228A1; US8493067B2

Abstract

Die Erfindung betrifft eine Erstellung eines Bilddatensatzes mittels einer radialen Abtastung eines Rohdatenraums (4) mit Hilfe einer Magnetresonanzanlage (5). Dabei werden folgende Schritte durchgeführt. Durchführen mindestens einer Kalibrierungsmessung für mindestens eine vorbestimmte Speiche (x; y) der radialen Abtastung, wobei eine Gradientenmomentdifferenz zwischen einem angenommenen und einem tatsächlich angelegten Gradientenmoment entlang der mindestens einen vorbestimmten Speiche (x; y) bestimmt wird. Auslesen aller Speichen (x, y, 1) des vorbestimmten Rohdatenraums (4), indem mehrere Magnetfeldgradienten in jeweils eine Raumrichtung geschaltet werden, um jeweils Abtastpunkte einer jeweiligen Speiche (x; y; 1) auszulesen. Korrigieren einer Position jedes Abtastpunkts jeder Speiche (x; y; 1) abhängig von der Gradientenmomentdifferenz, indem die aufgrund der jeweils geschalteten Magnetfeldgradienten angenommene Position des jeweiligen Abtastpunkts um die Gradientenmomentdifferenz verschoben wird.

Description

Die vorliegende Erfindung betrifft ein Verfahren zur Erstellung eines Bilddatensatzes mittels einer radialen Abtastung eines Rohdatenraums mit Hilfe einer Magnetresonanzanlage, wobei unbeabsichtigte zeitliche Verzögerungen der zu erzeugenden Gradientenfelder korrigiert werden. Darüber hinaus betrifft die vorliegende Erfindung eine entsprechend ausgestaltete Magnetresonanzanlage.
Üblicherweise wird in der Magnetresonanztomographie ein Rohdatenraum (auch als Messraum oder K-Raum bekannt) zeilenweise abgetastet. Allerdings sind auch andere Abtastmuster bekannt, wobei insbesondere die radiale Abtastung entlang von Speichen (d. h. durch das K-Raum-Zentrum verlaufende Geraden) in den letzten Jahren zunehmendes Interesse gewonnen hat. Die radiale Abtastung bietet verschiedene Vorteile, wie beispielsweise eine reduzierte Bewegungssensitivität oder die Möglichkeit, mit ultrakurzen Echozeiten (UTE) abzutasten. Die radiale Abtastung ist ein seit langer Zeit bekanntes Verfahren aus den frühen Tagen der Magnetresonanz-Technik, wobei sie sich allerdings bisher nicht auf breiter Ebene durchsetzen konnte. Ursächlich dafür sind in erster Linie inhärente technische Schwierigkeiten, welche beim Übergang der Abtasttrajektorie entlang von parallelen Linien auf gegenläufige, überlappende Speichen entstehen.
Ein Kernproblem stellen dabei zeitliche Verzögerungen der bei der Abtastung erzeugten Gradientenfelder dar, welche zu einer Abweichung zwischen der angenommenen und der tatsächlich abgetasteten Koordinate der Fouriertransformierten (d. h. im K-Raum) führen. Bei der konventionellen zeilenweisen Abtastung sind diese Abweichungen irrelevant, da hierbei alle Zeilen gleichermaßen bezüglich der Ausleserichtung verschoben werden. Diese Verschiebung führt aufgrund der Verschiebungseigenschaft der Fouriertransformation (Shift-Theorem) zu einer linearen Phasenmodulation des Objektes im Bildraum bei der zeilenweisen Abtastung. Da jedoch üblicherweise nur die Betragswerte (und nicht die Phase) bei der Bildgebung betrachtet werden, bleibt die Phasenmodulation bei zeilenweiser Abtastung ohne eine Auswirkung auf das dargestellte Bild.
Bei der radialen Abtastung führen die Gradientenverzögerungen ebenfalls zu einer linearen Phasenmodulation der abgetasteten Komponenten im Bildraum. Da sich jedoch die Ausleserichtung bei jeder Speiche von der Ausleserichtung einer jeweils beliebig anderen Speiche unterscheidet, ergibt sich jeweils eine unterschiedliche Phasenmodulation der enthaltenen räumlichen Information. Diese unterschiedliche Phasenmodulation führt zu teilweise starken Bildartefakten durch Interferenzeffekte, welche den diagnostischen Wert von mit radialer Abtastung erstellten MR-Bildern erheblich mindern.
Nach dem Stand der Technik existieren bisher keine genauen und schlüssigen Erkenntnisse über die physikalischen Ursachen der Gradientenverzögerungen. Einen großen Einfluss scheint das bauformbedingte Antwortverhalten von Gradientenspulen zu haben, da die beobachteten Gradientenverzögerungen meistens anisotrop sind, d. h. die Gradientenverzögerung der Gradientenspule in x-Richtung unterscheidet sich von der Gradientenverzögerung der Gradientenspule in y-Richtung. Darüber hinaus hängen die Gradientenverzögerungen von der gewählten Auslesegeschwindigkeit oder Bandbreite (”Bandwidth”) ab, was einerseits auf eine Amplitudenabhängigkeit oder auf Verzögerungen aufgrund der Digitalisierungshardware deuten könnte. Schließlich scheint es einen Einfluss von Wirbelstromeffekten sowie der System-Justage (Shim-Einstellungen) zu geben.
Nach dem Stand der Technik sind im Wesentlichen zwei Methoden zur Korrektur von Gradientenverzögerungen bekannt. Bei der ersten Methode werden die tatsächlich erzeugten Gradientenfelder (und somit die im Frequenzraum erzeugte Trajektorie) vermessen, welche dann anschließend für die Zuordnung der Messdaten im Frequenzraum verwendet werden. Für die dabei notwendigen Trajektorienmessungen sind nach dem Stand der Technik die beiden folgenden Druckschriften bezüglich der ersten Methode bekannt, welche für die Trajektorienmessungen teilweise spezielle Sensorhardware einsetzen.
"Simple Correction Method for k-Space Trajectory Deviations in MRI", J. H. Duyn, Y. Yang, J. A. Frank und J. W. van der Veen, JMR Volume 132, Issue 1, Mai 1998, Seiten 150–153.
"Spiral imaging artefact reduction: A comparision of two k-trajectory measurement methods", S. M. Lechner, P. T. Sipilä, F. Wiesinger, A. B. Kerr, M. W. Vogel, JMRI Volume 29, Issue 6, Seiten 1485–1492.
Bei der zweiten Methode wird die zu erwartende zeitliche Verschiebung abgeschätzt und abhängig davon das Moment des Dephasierungsgradienten derart angepasst, dass der tatsächliche Echozeitpunkt (Abtastung der Ursprungsposition des Frequenzraums) mit der angenommenen Echozeit zusammenfällt. Diese Methode ist in "Centering the Proection Reconstruction Trajectory: Reducing Gradient Delay Errors"; D. C. Peters, J. A. Derbyshire, E. R. McVeigh, Magn Reson Med., Juli 2003, 50(1): 1–6 vorgestellt, wobei die Verzögerung für ein festes Messprotokoll mit einer einmaligen Kalibrierungsmessung bestimmt wird.
In "Robust radial imaging with predetermined isotropic gradient delay correction, P. Speier, F. Trautwein, Proc. Intl. Soc. Mag. Reson. Med. 14 (2006) 2379 wird die zu erwartende Verzögerung anhand eines linearen Modells ermittelt, so dass es nicht notwendig ist, eine neue Kalibrierung bei einer Änderung der Messparameter durchzuführen. In der Praxis hat sich jedoch gezeigt, dass die entsprechende Korrektur nicht ausreicht, da die ursächlichen Effekte anscheinend systemabhängig, justageabhängig und teilweise auch patientenabhängig sind. Daher reicht die sich ergebende Bildqualität nicht für einen klinischen Einsatz, insbesondere nicht für morphologische Untersuchungen, aus.
Die vorliegende Erfindung stellt sich daher die Aufgabe, Gradientenverzögerungen bei einer radialen Abtastung besser zu korrigieren, als dies nach dem Stand der Technik möglich ist.
Erfindungsgemäß wird diese Aufgabe durch ein Verfahren zur Erstellung eines Bilddatensatzes mittels einer radialen Abtastung des Rohdatenraums mit Hilfe einer Magnetresonanzanlage nach Anspruch 1, durch eine Magnetresonanzanlage nach Anspruch 12, durch ein Computerprogrammprodukt nach Anspruch 14 und durch einen elektronisch lesbaren Datenträger nach Anspruch 15 gelöst. Die abhängigen Ansprüche definieren bevorzugte und vorteilhafte Ausführungsformen der vorliegenden Erfindung.
Im Rahmen der vorliegenden Erfindung wird ein Verfahren zur Erstellung eines Bilddatensatzes mit einer radialen Abtastung mittels einer Magnetresonanzanlage dargestellt. Dabei umfasst das Verfahren folgende Schritte:
Für eine oder für mehrere vorbestimmte Speichen wird jeweils eine oder jeweils mehrere Kalibrierungsmessungen durchgeführt. Bei jeder Kalibrierungsmessung wird dabei eine Gradientenmomentdifferenz zwischen einem angenommenen und einem tatsächlich angelegten Gradientenmoment entlang der jeweiligen Speiche bestimmt. Da das angenommene bzw. tatsächlich angelegte Gradientenmoment dem angenommenen bzw. tatsächlich abgetasteten Abtastpunkt im Frequenzraum entspricht, kann aus der Gradientenmomentdifferenz jederzeit die aufgrund der Gradientenverzögerung resultierende Pixelverschiebung zwischen dem angenommenen Abtastpunkt und dem tatsächlichen Abtastpunkt im K-Raum oder Frequenzraum bestimmt werden.
Anschließend werden alle Speichen eines vorbestimmten Rohdatenraums ausgelesen, wobei mehrere Magnetfeldgradienten in jeweils einer Raumrichtung (d. h. zwei Magnetfeldgradienten zum Auslesen einer Schicht oder drei Magnetfeldgradienten zum Auslesen eines dreidimensionalen Raumsegments) geschaltet werden, um jeweils Abtastpunkte einer jeweiligen Speiche zu erfassen.
Nach dem Auslesen des entsprechenden Abtastpunkts wird eine Position jedes Abtastpunkts jeder Speiche abhängig von der Gradientenmomentdifferenz korrigiert. Dazu wird die angenommene Position des jeweiligen Abtastpunkts, welche durch die jeweils geschalteten Magnetfeldgradienten definiert ist, um die Gradientenmomentdifferenz verschoben.
Das erfindungsgemäße Verfahren ermöglicht einen praktischen Einsatz einer radialen Abtastung, da die unerwünschten zeitlichen Verzögerungen der erzeugten Gradientenfelder bei der Bestimmung der tatsächlichen Position des jeweiligen Abtastpunkts korrigiert werden.
Zur Bestimmung der Gradientenmomentdifferenz gibt es mehrere Möglichkeiten:

1. Die Gradientenmomentdifferenz wird nur anhand genau einer vorbestimmten Speiche ermittelt. Dazu kann eine beliebige Speiche gewählt werden; allerdings wird die Speiche meist in der Richtung der zwei oder drei Gradientenspulen gewählt werden.
2. Die Gradientenmomentdifferenz wird anhand mehrerer vorbestimmter Speichen ermittelt. In der Regel entspricht die Anzahl der vorbestimmten Speichen bei dieser Möglichkeit der Anzahl der Magnetfeldgradienten oder der Gradientenspulen, wobei sich jeweils eine dieser zwei oder drei Speichen in der Richtung des jeweiligen Magnetfeldgradienten oder in der Richtung der jeweiligen Gradientenspule erstreckt. Es ist allerdings auch möglich, dass die Anzahl der vorbestimmten Speichen die Anzahl der Magnetfeldgradienten oder Gradientenspulen übersteigt. Bei dieser Möglichkeit wird die Gradientenmomentdifferenz als Mittelwert derjenigen Gradientenmomentdifferenzen berechnet, welche für die vorbestimmten Speichen bestimmt worden sind.
3. Die Anzahl der vorbestimmten Speichen, anhand derer die Gradientenmomentdifferenz bestimmt wird, entspricht der Anzahl der Magnetfeldgradienten oder Gradientenspulen. Dabei erstreckt sich jeweils eine dieser zwei oder drei Speichen in der Richtung des jeweiligen Magnetfeldgradienten oder der jeweiligen Gradientenspule. Es wird für jede Speiche eine Gradientenmomentdifferenz bestimmt, so dass sich die Gradientenmomentdifferenz quasi aus zwei oder drei Komponenten zusammensetzt.

Während die Möglichkeiten 1 und 2 nur eine isotrope Korrektur der Gradientenverzögerungen ermöglicht, ermöglicht die Möglichkeit 3 auch eine anisotrope Korrektur. Natürlich ist es bei allen vorab beschriebenen Möglichkeiten möglich, dass die Gradientenmomentdifferenz für die jeweilige vorbestimmte Speiche mehrfach bestimmt wird, um dann beispielsweise mit einem Mittelwert dieser Gradientenmomentdifferenzen als der Gradientenmomentdifferenz für diese Speiche zu arbeiten.
Bei der im Folgenden beschriebenen erfindungsgemäßen Ausführungsform erfolgt eine isotrope Korrektur, d. h. die Gradientenmomentdifferenz wird gemäß der vorab beschriebenen Möglichkeiten 1 oder 2 bestimmt. Dabei sei angenommen, dass die angenommene Position eines Abtastpunkts durch eine x-Koordinate k'_x und durch eine y-Koordinate k'_y und die korrigierte Position durch eine x-Koordinate k_x und durch eine y-Koordinate k_y definiert wird. Dann berechnet sich abhängig von der Gradientenmomentdifferenz GMD und dem Winkel φ, welcher zwischen der jeweiligen Speiche, auf welchem der Abtastpunkt liegt, und der x-Achse vorhanden ist, die x-Koordinate k_x der korrigierten Position aus der x-Koordinate k'_x der angenommenen Position gemäß Gleichung (1) und die y-Koordinate k_y der korrigierten Position aus der y-Koordinate k'_y der angenommenen Position gemäß Gleichung (2). k_x = k'_x – GMD × cos(φ) (1) k_y = k'_y – GMD × sin(φ) (2)
Mit anderen Worten wird bei dieser Ausführungsform die angenommene Position des jeweiligen Abtastpunkts um die Gradientenmomentdifferenz entlang der Speiche, auf welcher sich der Abtastpunkt befindet, verschoben.
Über die folgende Gleichung (3) kann die Gradientenmomentdifferenz GMD in eine Pixelverschiebung PV umgerechnet werden.
Dabei steht N für die Anzahl der Abtastpunkte entlang einer Speiche und k'_max für das maximale angenommene Gradientenmoment. Die Verzögerung wird demnach in Form einer Pixel- oder Sampleverschiebung bestimmt. Eine Pixelverschiebung PV von 2 bedeutet z. B., dass das Echozentrum 2 Sampleabstände neben der angenommenen Echoposition liegt.
Die Pixelverschiebung oder Echoverschiebung kann bei einer Bildrekonstruktion verwendet werden, um Fehler bei der Positionsbestimmung der Abtastpunkte auszugleichen. Ähnlich wie bei den Gleichungen (1) und (2) erfordern die radialen Trajektorien eine spezielle Rekonstruktionsmethode, welche auch als ”Gridding” bekannt ist, wobei die Messdaten von den jeweiligen Speichenpositionen auf ein reguläres Gitter interpoliert werden. Dafür wird die Position der einzelnen Abtastpunkte im Frequenzraum gemäß der folgenden Gleichungen (4) und (5) berechnet. k_x = (–k_max + i × Δk) × cos(φ) (4) k_y = (–k_max + i × Δk) × sin(φ) (5)
Dabei steht k_max für ein um die Gradientenmomentdifferenz korrigiertes maximales Gradientenmoment, i steht für den jeweiligen Index der Abtastposition (d. h. i läuft von 0 bis N) und Δk ist die Differenz zwischen dem Gradientenmoment eines Abtastpunkts und dem Gradientenmoment desjenigen Abtastpunkts, welcher ihm direkt benachbart ist. Demnach lässt sich Δk durch folgende Gleichung (6) berechnen:
Um die Gradientenverzögerung bei den Gleichungen (4) und (5) zu berücksichtigen, wird durch die folgende Gleichung (7) quasi die Außenposition ±k_max der jeweiligen Speiche um die ermittelte Verzögerung oder Pixelverschiebung PV verschoben.
Werden mit den Gleichungen (4) bis (7) die Positionen der Abtastpunkte im K-Raum bestimmt, werden die an den entsprechenden Abtastpunkten erfassten Messdaten Gitterpositionen zugeordnet, welche näherungsweise mit den tatsächlichen Frequenzraum-Positionen zum Abtastzeitpunkt übereinstimmen, d. h. Gradientenverzögerungen werden dabei berücksichtigt oder korrigiert.
Die vorab beschriebenen Ausführungsformen, mit welchen die korrigierten Positionen der Abtastpunkte bestimmt werden, gelten für eine zweidimensionale radiale Abtastung. Selbstverständlich ist die vorliegende Erfindung auch für eine dreidimensionale radiale Abtastung gültig. In diesem Fall werden drei Gleichungen zur Ermittlung der Koordinaten k_x, k_y, k_z benötigt. Die Gleichung (7) behält bei einer isotropen Korrektur ihre Gültigkeit.
Im Folgenden wird eine erfindungsgemäße Ausführungsform zur Bestimmung der Gradientenmomentdifferenz entlang einer der vorbestimmten Speichen beschrieben. Dazu wird ein erstes Signal S₀(t) entlang dieser Speiche in einer vorbestimmten Richtung (es gibt nur zwei Möglichkeiten für diese Richtung, nämlich die Hinrichtung und die Rückrichtung entlang der Speiche) über eine vorbestimmte Strecke dieser Speiche ausgelesen. Darüber hinaus wird ein zweites Signal S₁₈₀(t) entlang dieser Speiche in einer der vorbestimmten Richtung entgegengesetzten Richtung (also in der anderen Richtung) über dieselbe Strecke ausgelesen.
Legt man die x-Achse entlang dieser Speiche, so entspricht das erste Signal S₀(t) einer Messung mit dem Winkel 0° und das zweite Signal S₁₈₀(t) einer Messung mit dem Winkel 180°. Wenn keine Gradientenverzögerung vorliegt, ist das erste Signal S₀(t) lediglich bezüglich der Reihenfolge gegenüber dem zweiten Signal S₁₈₀(t) vertauscht. Mit anderen Worten sind beide Signale identisch, wenn man die Sortierung der aufgenommenen Werte des zweiten Signals S₁₈₀(t) umdreht. Wird nun ein weiteres zweites Signal S'₁₈₀(t) gemäß der folgenden Gleichung (8) definiert, bei welchem die Reihenfolge der Samples gegenüber dem zweiten Signal S₁₈₀(t) vertauscht ist, dann ist dieses weitere zweite Signal S'₁₈₀(t) gleich dem ersten Signal S₀(t), wenn keine Gradientenverzögerung vorhanden ist. S'₁₈₀(t) = S₁₈₀(N – 1 – t) (8)
Mit N gleich der Anzahl der Samples (d. h. t = 0, 1, ..., N – 1) auf der entsprechenden Speiche bzw. in diesem Fall auf der x-Achse.
Liegt allerdings eine Gradientenverzögerung vor, so existiert die doppelte Pixelverschiebung 2 × PV zwischen den beiden Signalen S₀(t) und S'₁₈₀(t), so dass folgende Gleichung (9) gilt. S₀(t) = S'₁₈₀(t + 2 × PV) (9)
Mit anderen Worten ergibt sich, da sowohl das erste Signal S₀(t) als auch das weitere zweite Signal S'₁₈₀(t) verzögert sind, die Echoverzögerung oder Pixelverschiebung PV aus dem halben Pixelversatz zwischen dem ersten Signal S₀(t) und dem weiteren zweiten Signal S'₁₈₀(t).
Zur Bestimmung der Pixelverschiebung PV aus den Kalibrierungsmessungen wird eine Kreuzkorrelationsanalyse durchgeführt, wie es im Folgenden beschrieben wird.
Um das erfindungsgemäße Verfahren robust gegenüber Phasenmodulationen auszugestalten, wird nach der Umsortierung zunächst der Betrag der komplexen Signale S₀(t) und S'₁₈₀(t) berechnet. Anschließend wird die Fouriertransformierte der Beträge für beide Signale berechnet. Durch die bereits vorab angesprochene Verschiebungseigenschaft der Fouriertransformation wird der Pixelversatz von S'₁₈₀(t) in eine Phasenmodulation gegenüber der Fouriertransformierten von S₀(t) übersetzt. Diese Phasenmodulation wird ermittelt, indem die Fouriertransformierte des Betrages von S₀(t) mit der komplexkonjugierten Fouriertransformierten des Betrages von S'₁₈₀(t) multipliziert wird, wodurch sich eine Funktion g(x) gemäß Gleichung (10) ergibt. g(x) = FT(|S₀(t)|) × Conj(FT(|S'₁₈₀(t)|)) (10)
Abhängig von dieser Funktion g(x) kann nun die Gradientenmomentdifferenz bezüglich dieser Speiche bestimmt werden.
Gemäß einem ersten erfindungsgemäßen Ansatz wird dazu die inverse Fouriertransformation von g(x) bestimmt, wodurch man die Kreuzkorrelationsfunktion von S₀(t) und S'₁₈₀(t) erhält. Ein Abstand zwischen der Position des Maximums dieser Kreuzkorrelationsfunktion invFT(g(x)) und dem k-Raum-Zentrum entspricht der zweifachen Pixelverschiebung 2 × PV.
Gemäß einem zweiten erfindungsgemäßen Ansatz wird die Steigung des Phasenverlaufs von g(x) bestimmt. Bei einer einfachen (d. h. gleichartigen Verzögerung entlang der Speiche) ist dieser Phasenverlauf rein linear. Da jedoch der Messraum in der Regel größer als das Messobjekt ist, ist die Phase von g(x) nur innerhalb des Messobjekts bestimmt. Um diesem Sachverhalt Rechnung zu tragen, wird in einem ersten Schritt das Betragsmaximum von g(x) ermittelt. Ausgehend von diesem Maximum wird ein Intervall derart bestimmt, dass der Betrag von g(x) in diesem Intervall, in welchem sich auch das Maximum befindet, niemals kleiner als ein vorbestimmter Betragsschwellenwert ist. Dabei ist dieser Betragsschwellenwert einen vorbestimmten Prozentsatz (z. B. 5%) kleiner als das Betragsmaximum. Innerhalb dieses Intervalls wird anschließend eine lineare Funktion beispielsweise mit einer linearen Ausgleichsrechnung an den Phasenverlauf von g(x) angepasst (gefittet). Aus der Steigung S der linearen Funktion bzw. Geraden ergibt. sich dann die Pixelverschiebung PV gemäß der folgenden Gleichung (11).
Mit Hilfe der Gleichung (3) lässt sich daraus gemäß der folgenden Gleichung (12) die Gradientenmomentdifferenz GMD entlang der Speiche bestimmen.
Um eine anisotrope Korrektur der Verzögerung durchzuführen, werden bei einer zweidimensionalen radialen Abtastung eine Kalibrierungsmessung entlang der x-Achse zur Ermittlung der Gradientenmomentdifferenz GMD_x entlang der x-Achse und eine Kalibrierungsmessung entlang der y-Achse zur Ermittlung der Gradientenmomentdifferenz GMD_y entlang der y-Achse durchgeführt. Bei einer dreidimensionalen radialen Abtastung müsste zusätzlich noch eine Kalibrierungsmessung entlang der z-Achse zur Ermittlung der Gradientenmomentdifferenz GMD_z entlang der z-Achse durchgeführt werden. Über die Gleichung (3) kann die entsprechende Gradientenmomentdifferenz GMD_x, GMD_y bzw. GMD_z einfach in die entsprechende Pixelverschiebung PV_x, PV_y bzw. PV_z umgerechnet werden. Zur anisotropen Korrektur wird nun die Gleichung (7) durch die folgende Gleichung (13) ersetzt.
Über die Gleichungen (4) und (5) kann dann unter Berücksichtigung der Gleichung (13) die anisotrop korrigierte Position des entsprechenden Abtastpunkts bestimmt werden.
Es ist aber auch möglich, aus der angenommenen Position (k'_x, k'_y) des jeweiligen Abtastpunkts abhängig von den Gradientenmomentdifferenzen GMD_x und GMD_y entlang der x-Achse und der y-Achse die korrigierte Position (k_x, k_y) des Abtastpunkts gemäß der folgenden Gleichungen (14) und (15) zu bestimmen.
Eine anisotrope Gradientenverzögerung kann dazu führen, dass Speichen, welche durch eine Überlagerung der X- und Y-Gradienten (und gegebenenfalls mit einer Überlagerung des Z-Gradienten im dreidimensionalen Fall) erzeugt werden, gegenüber dem K-Raum-Zentrum versetzt sind. Diese Verschiebung gegenüber dem K-Raum-Zentrum lässt sich durch die vorab mit den Gleichungen (13) bis (15) beschriebenen Ausführungsformen korrigieren.
Die Durchführung der Kalibrierungsmessungen zur Bestimmung der Gradientenmomentdifferenzen bzw. Pixelverschiebungen kann in einer Vorbereitungsphase der eigentlichen MR-Messung durchgeführt werden. Insbesondere bei gradientenechobasierten Messsequenzen, welche eine kurze Repetitionszeit verwenden (z. B. 3D FLASH oder 2D FLASH), ist es üblicherweise erforderlich, dass vor der Aufnahme der eigentlichen Messdaten so genannte ”Prep-Shots” durchgeführt werden, damit die Magnetisierung innerhalb des Messobjekts in einen Gleichgewichtszustand gebracht wird. Diese vorbereitenden Messungen sind insbesondere bei radialen Sequenzen wichtig, da alle Speichen gleich wichtige Informationen (Messdaten) zum Gesamtergebnis beitragen, so dass auch die über die ersten Speichen erfassten Informationen für das Gesamtergebnis wichtig sind. Werden keine vorbereitenden Messungen (Prep-Shots) durchgeführt, so ist die Signalstärke bei den Messungen der ersten Speichen durch die noch nicht vorhandene Absättigung der Magnetisierung im Messobjekt zu hoch, was zu so genannten Verschmierungsartefakten führt.
Bei beispielsweise 200 Prep-Shots können 100 Gradientenmomentdifferenzen entlang einer Speiche berechnet werden. Da bei einer zweidimensionalen radialen Abtastung im Wesentlichen nur die Speichen entlang der x-Gradientenspule und entlang der y-Gradientenspule interessant sind, kann also die Gradientenmomentdifferenz entlang der x-Achse und die Gradientenmomentdifferenz entlang der y-Achse jeweils 50-Mal bestimmt werden. Die endgültige Gradientenmomentdifferenz entlang der x-Achse und die endgültige Gradientenmomentdifferenz entlang der y-Achse werden dann über eine Mittelwertbildung berechnet. Dazu besteht zum einen die Möglichkeit die erfassten Signale (z. B. S₀(t), S₁₈₀(t)) zu mitteln und anhand dieser gemittelten Signale dann die endgültige Gradientenmomentdifferenz zu bestimmen. Die andere Möglichkeit besteht darin, aus jeweils zwei entsprechenden Signalen (z. B. S₀(t), S₁₈₀(t)) die jeweilige Gradientenmomentdifferenz zu bestimmen und dann aus diesen z. B. 50 Gradientenmomentdifferenzen die endgültige Gradientenmomentdifferenz über Mittelwertbildung zu bestimmen.
Gemäß einer erfindungsgemäßen Ausführungsform kann eine beliebige vorher beschriebene Ausprägung des erfindungsgemäßen Verfahrens mit dem vorher beschriebenen Verfahren zur Vorkompensation der Gradientenverzögerung durch Anpassung des Moments des Vordephasierungsgradienten (siehe D. C. Peters und andere oder P. Speier und andere) kombiniert werden. Die vorliegende Erfindung korrigiert dann insbesondere diejenigen Verzögerungen der Gradientenfelder, welche trotz des angepassten Vordephasierungsgradienten noch vorhanden sind.
Im Rahmen der vorliegenden Erfindung wird auch eine Magnetresonanzanlage zur Erstellung eines Bilddatensatzes mittels einer radialen Abtastung bereitgestellt. Dabei umfasst die Magnetresonanzanlage eine Ansteuereinheit zur Ansteuerung eines Tomographen der Magnetresonanzanlage, eine Empfangsvorrichtung zum Empfang von von dem Tomographen aufgenommenen Signalen und eine Auswertevorrichtung zur Auswertung der Signale und zur Erstellung des Bilddatensatzes. Die Magnetresonanzanlage ist derart ausgestaltet, dass die Magnetresonanzanlage eine oder mehrere Kalibrierungsmessungen für eine oder mehrere vorbestimmte Speichen durchführt, um eine Gradientenmomentdifferenz entlang der einen oder der mehreren Speichen zu bestimmen. Die Magnetresonanzanlage liest alle Speichen des Rohdatenraums aus, indem die Magnetresonanzanlage mehrere Magnetfeldgradienten in jeweils eine Raumrichtung schaltet, um jeweils Abtastpunkte einer dieser Speichen auszulesen. Die Magnetresonanzanlage korrigiert eine Position jedes Abtastpunkts jeder Speiche abhängig von der Gradientenmomentdifferenz, indem die Magnetresonanzanlage eine aufgrund der jeweils geschalteten Magnetfeldgradienten angenommene Position des jeweiligen Abtastpunkts um die Gradientenmomentdifferenz verschiebt.
Die Vorteile der erfindungsgemäßen Magnetresonanzanlage entsprechen dabei im Wesentlichen den Vorteilen des erfindungsgemäßen Verfahrens, welche vorab im Detail ausgeführt worden sind, so dass hier auf eine Wiederholung verzichtet wird.
Des Weiteren beschreibt die vorliegende Erfindung ein Computerprogrammprodukt, insbesondere ein Computerprogramm oder eine Software, welche man in einen Speicher einer programmierbaren Steuerung bzw. einer Recheneinheit einer Magnetresonanzanlage laden kann. Mit diesem Computerprogrammprodukt können alle oder verschiedene vorab beschriebene Ausführungsformen des erfindungsgemäßen Verfahrens ausgeführt werden, wenn das Computerprogrammprodukt in der Steuerung läuft. Dabei benötigt das Computerprogrammprodukt eventuell Programmmittel, z. B. Bibliotheken und Hilfsfunktionen, um die entsprechenden Ausführungsformen der Verfahren zu realisieren. Mit anderen Worten soll mit dem auf das Computerprogrammprodukt gerichteten Anspruch insbesondere ein Computerprogramm oder eine Software unter Schutz gestellt werden, mit welcher eine der oben beschriebenen Ausführungsformen des erfindungsgemäßen Verfahrens ausgeführt werden kann bzw. welche diese Ausführungsform ausführt. Dabei kann es sich bei der Software um einen Quellcode (z. B. C++), der noch compiliert und gebunden oder der nur interpretiert werden muss, oder um einen ausführbaren Softwarecode handeln, der zur Ausführung nur noch in die entsprechende Recheneinheit zu laden ist.
Schließlich offenbart die vorliegende Erfindung einen elektronisch lesbaren Datenträger, z. B. eine DVD, ein Magnetband oder einen USB-Stick, auf welchem elektronisch lesbare Steuerinformationen, insbesondere Software (vgl. oben), gespeichert ist. Wenn diese Steuerinformationen (Software) von dem Datenträger gelesen und in eine Steuerung bzw. Recheneinheit einer Magnetresonanzanlage gespeichert werden, können alle erfindungsgemäßen Ausführungsformen des vorab beschriebenen Verfahrens durchgeführt werden.
Die vorliegende Erfindung ist vorteilhafterweise adaptiv, benötigt keine Benutzerinteraktion und erfordert keine vorherigen Kalibrierungsschritte, so dass die Messzeit nicht verlängert wird. Die vorliegende Erfindung korrigiert vorteilhafterweise gerätespezifische und patientenabhängige Verzögerungen, welche nach dem Stand der Technik ein deutliches Problem für den routinemäßigen Einsatz von radialen Messtechniken darstellen.
Die vorliegende Erfindung ist insbesondere für MR-Verfahren, welche mit einer radialen Abtastung arbeiten, geeignet. Darunter fallen auch MR-Verfahren, bei welchen eine Speiche in mehreren Schritten, d. h. mit mehreren HF-Anregungen abgetastet wird. Selbstverständlich ist die vorliegende Erfindung nicht auf diese bevorzugten Anwendungsbereiche eingeschränkt, da die vorliegende Erfindung beispielsweise auch bei einer Einzelpunkt-Bildgebung eingesetzt werden könnte.
Im Folgenden wird die vorliegende Erfindung anhand bevorzugter erfindungsgemäßer Ausführungsformen im Detail mit Bezug zu den Figuren beschrieben.
In 1 sind eine radiale Abtastung und eine zeilenweise Abtastung dargestellt.
In 2 ist schematisch eine erfindungsgemäße Magnetresonanzanlage dargestellt.
In 3 ist ein Flussplan eines erfindungsgemäßen Verfahrens dargestellt.
In 1a ist eine radiale Abtastung und in 1b eine konventionelle zeilenweise Abtastung dargestellt. Bei der gewöhnlicherweise eingesetzten zeilenweisen Abtastung wird ein Messraum 4 oder eine Schicht des K-Raums Zeile für Zeile abgetastet, wobei pro Zeile mehrere Punkte erfasst werden. Dadurch wird die abgetastete Schicht sowohl in horizontaler als auch in vertikaler Richtung gleichmäßig abgetastet. D. h. Abtastpunkte haben in horizontaler Richtung und in vertikaler Richtung meist denselben Abstand zu ihrem Nachbarn.
Dagegen wird bei der in 1a dargestellten zweidimensionalen radialen Abtastung der Messraum 4 oder die abgetastete Schicht mittels Speichen 1 abgetastet, welche wie die Speichen eines Rades durch das Zentrum Z der Schicht verlaufen. Dabei werden pro Speiche 1 mehrere Punkte erfasst, welche normalerweise äquidistant auf der Speiche 1 angeordnet sind.
Daher werden bei einer radialen Abtastung im Bereich um das Zentrums Z auf derselben Fläche mehr Abtastpunkte erfasst als in den Randbereichen des Messraums 4.
Zur Orientierung sind in 1a eine Speiche x als x-Achse und eine dazu senkrecht stehende Speiche y als y-Achse gekennzeichnet. Diese Bezeichnung ist willkürlich, d. h. es können zwei beliebige aufeinander senkrecht stehende Speichen 1 als x-Achse und y-Achse gekennzeichnet werden. Ein Winkel φ bezeichnet denjenigen Winkel, welchen die entsprechende Speiche 1 mit der x-Achse einschließt.
Bei einer dreidimensionalen radialen Abtastung würden die Speichen 1 ebenfalls durch dasselbe Zentrum verlaufen. In diesem Fall liegen die Speichen allerdings nicht in einer Ebene, sondern sind derart angeordnet, dass ihre Enden auf einer Art Kugeloberfläche liegen, wobei das Zentrum der Mittelpunkt der zugehörigen Kugel ist.
In 2 ist eine erfindungsgemäße Magnetresonanzanlage 5 schematisch dargestellt. Die Magnetresonanzanlage 5 umfasst im Wesentlichen einen Tomograph 3, mit welchem das für die MR-Untersuchung notwendige Magnetfeld in einem Messraum 4 erzeugt wird, einen Tisch oder Liegenbrett 2, eine Steuereinrichtung 6, mit welcher der Tomograph 3 gesteuert wird und MR-Daten von dem Tomograph 3 erfasst werden, und ein an die Steuereinrichtung 6 angeschlossenes Terminal 7.
Die Steuereinrichtung 6 umfasst ihrerseits eine Ansteuereinheit 11, eine Empfangsvorrichtung 12 und eine Auswertevorrichtung 13. Während der Erstellung eines Bilddatensatzes werden MR-Daten mittels des Tomograph 3 von der Empfangsvorrichtung 12 erfasst, wobei der Tomograph 3 und der Tisch 2 von der Ansteuereinheit 11 derart angesteuert werden, dass MR-Daten in einem Messvolumen, welches sich im Körperinneren eines auf dem Tisch 2 liegenden Patienten O befindet, erfasst werden.
Die Auswertevorrichtung 13 bereitet dann die MR-Daten derart auf, dass sie auf einem Bildschirm 8 des Terminals 7 grafisch dargestellt werden können und dass erfindungsgemäß erstellte Bilder, bei welchen Verzögerungen der erzeugten Gradientenfelder korrigiert sind, angezeigt werden. Neben der grafischen Darstellung der MR-Daten kann mit dem Terminal 7, welches neben dem Bildschirm 8 eine Tastatur 9 und eine Maus 10 umfasst, von einem Anwender z. B. ein zu vermessender dreidimensionaler Volumenabschnitt oder eine quasi zweidimensionale Schicht vorgegeben werden und weitere Parameter zur Durchführung des erfindungsgemäßen Verfahrens bestimmt werden. Über das Terminal 7 kann auch die Software für die Steuereinrichtung 6 in die Steuereinrichtung 6 geladen werden. Diese Software der Steuereinrichtung 6 kann dabei auch das erfindungsgemäße Verfahren umfassen. Es ist dabei auch möglich, dass ein erfindungsgemäßes Verfahren in einer Software enthalten ist, welche in dem Terminal 7 abläuft. Unabhängig davon, in welcher Software das erfindungsgemäße Verfahren enthalten ist, kann die Software auf einer DVD 14 gespeichert sein, so dass diese Software dann von dem Terminal 7 von der DVD 14 gelesen und entweder in die Steuereinrichtung 6 oder in eine Recheneinheit des Terminals 7 selbst kopiert werden kann.
In 3 ist ein Flussplan eines erfindungsgemäßen Verfahrens abgebildet.
Im Schritt S1 werden so genannte Prep-Shots durchgeführt, um die Magnetisierung innerhalb des Messobjekts in einen Gleichgewichtszustand zu bringen. Bei diesen Prep-Shots werden für die auf der x-Achse und auf der y-Achse liegenden Speichen mehrfach die Gradientenmomentdifferenzen zwischen einem angenommenen und einem tatsächlich angelegten Gradientenmoment z. B. mittels des mit den Gleichungen (8) bis (11) beschriebenen erfindungsgemäßen Ansatzes bestimmt.
Im Schritt S2 werden mehrere Abtastpunkte entlang einer der Speichen 1 abgetastet. Anschließend wird die Position jedes ausgelesenen Abtastpunkts mit Hilfe der in Schritt S1 bestimmten Gradientenmomentdifferenz(en) korrigiert. Diese Korrektur erfolgt insbesondere mit einem mit den Gleichungen (1) bis (7) und (13) bis (15) beschriebenen erfindungsgemäßen Ansatz.
Im Schritt S4 wird überprüft, ob bereits alle Abtastpunkte im Messraum erfasst worden sind, d. h. ob alle Abtastpunkte aller Speichen ausgelesen worden sind. Wenn dies nicht der Fall ist (nein bei Schritt S4), kehrt das Verfahren zu dem Schritt S2 zurück, in welchem noch fehlende Abtastpunkte einer bereits abgetasteten Speiche oder Abtastpunkte einer weiteren noch nicht abgetasteten Speiche 1 ausgelesen werden. Wenn bereits alle Abtastpunkte erfasst worden sind (ja bei Schritt S4), endet das Verfahren.
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 Nicht-Patentliteratur

”Simple Correction Method for k-Space Trajectory Deviations in MRI”, J. H. Duyn, Y. Yang, J. A. Frank und J. W. van der Veen, JMR Volume 132, Issue 1, Mai 1998, Seiten 150–153 [0007]
”Spiral imaging artefact reduction: A comparision of two k-trajectory measurement methods”, S. M. Lechner, P. T. Sipilä, F. Wiesinger, A. B. Kerr, M. W. Vogel, JMRI Volume 29, Issue 6, Seiten 1485–1492 [0008]
”Centering the Proection Reconstruction Trajectory: Reducing Gradient Delay Errors”; D. C. Peters, J. A. Derbyshire, E. R. McVeigh, Magn Reson Med., Juli 2003, 50(1): 1–6 [0009]
”Robust radial imaging with predetermined isotropic gradient delay correction, P. Speier, F. Trautwein, Proc. Intl. Soc. Mag. Reson. Med. 14 (2006) 2379 [0010]
D. C. Peters [0045]
P. Speier [0045]

Claims

Verfahren zur Erstellung eines Bilddatensatzes mittels einer radialen Abtastung eines Rohdatenraumes (4) mit Hilfe einer Magnetresonanzanlage (5), wobei das Verfahren folgende Schritte umfasst: Durchführen mindestens einer Kalibrierungsmessung für mindestens eine vorbestimmte Speiche (x; y) der radialen Abtastung, wobei eine Gradientenmomentdifferenz zwischen einem angenommenen und einem tatsächlich angelegten Gradientenmoment entlang der mindestens einen vorbestimmten Speiche (x; y; 1) bestimmt wird, Auslesen aller Speichen (x, y, 1) des vorbestimmten Rohdatenraums (4), indem mehrere Magnetfeldgradienten in jeweils einer Raumrichtung geschaltet werden, um jeweils Abtastpunkte einer jeweiligen Speiche (x; y; 1) auszulesen, und Korrigieren einer Position jedes Abtastpunkts jeder Speiche (x; y; 1) abhängig von der Gradientenmomentdifferenz, indem die aufgrund der jeweils geschalteten Magnetfeldgradienten angenommene Position des jeweiligen Abtastpunkts um die Gradientenmomentdifferenz verschoben wird.
Verfahren nach Anspruch 1, dadurch gekennzeichnet, dass die mindestens eine vorbestimmte Speiche (x; y) genau eine vorbestimmte Speiche ist.
Verfahren nach Anspruch 1, dadurch gekennzeichnet, dass die mindestens eine vorbestimmte Speiche mehrere vorbestimmte Speichen (x, y) umfasst, dass eine Anzahl der vorbestimmten Speichen (x, y) gleich einer Anzahl der Magnetfeldgradienten ist, und dass sich die vorbestimmten Speichen (x, y) jeweils in der Richtung der Magnetfeldgradienten erstrecken.
Verfahren nach Anspruch 3, dadurch gekennzeichnet, dass pro vorbestimmte Speiche (x; y) eine Speichen-Gradientenmomentdifferenz bestimmt wird, und dass die Gradientenmomentdifferenz als Mittelwert der Speichen-Gradientenmomentdifferenzen bestimmt wird.
Verfahren nach einem der vorherigen Ansprüche, dadurch gekennzeichnet, dass die angenommene Position durch eine x-Koordinate k'_x und eine y-Koordinate k'_y definiert wird, dass die korrigierte Position durch eine x-Koordinate k_x und eine y-Koordinate k_y definiert wird, dass die x-Koordinate k_x der korrigierten Position wie folgt bestimmt wird k_x = k'_x – GMD × cos(φ), und dass die y-Koordinate k_y der korrigierten Position wie folgt bestimmt wird k_y = k'_y – GMD × sin(φ), wobei GMD die Gradientenmomentdifferenz und φ ein Winkel ist, welcher zwischen der jeweiligen Speiche (x; y; 1) und der x-Achse ausgebildet ist.
Verfahren nach einem der vorherigen Ansprüche, dadurch gekennzeichnet, dass zur Bestimmung der Gradientenmomentdifferenz entlang einer der vorbestimmten Speichen ein erstes Signal S₀(t) entlang dieser Speiche (x; y) in einer vorbestimmten Richtung über eine vorbestimmte Strecke dieser Speiche (x; y) ausgelesen wird, dass ein zweites Signal S₁₈₀(t) entlang dieser Speiche (x; y) in einer der vorbestimmten Richtung entgegengesetzten Richtung über die vorbestimmte Strecke dieser Speiche (x; y) ausgelesen wird, dass aus dem zweiten Signal S₁₈₀(t) ein weiteres zweites Signal S'₁₈₀(t) entlang dieser Speiche (x; y) bestimmt wird, indem die Reihenfolge der Abtastwerte des weiteren zweiten Signals S'₁₈₀(t) gegenüber dem zweiten Signal S₁₈₀(t) vertauscht wird, dass gemäß folgender Gleichung eine Funktion g(x) bestimmt wird, indem die Fouriertransformierte des Betrags des ersten Signals S₀(t) mit der komplexkonjugierten Fouriertransformierten des Betrags des weiteren zweiten Signals S'₁₈₀(t) multipliziert wird: g(x) = FT(|S₀(t)|) × Conj(FT(|S'₁₈₀(t)|)), und dass abhängig von g(x) die Gradientenmomentdifferenz dieser Speiche bestimmt wird.
Verfahren nach Anspruch 6, dadurch gekennzeichnet, dass ein Ort eines Maximums einer inversen Fouriertransformierten von g(x) bestimmt wird, und dass die Gradientenmomentdifferenz als Abstand zwischen dem Ort des Maximums und dem Zentrum des Rohdatenraums (4) bestimmt wird.
Verfahren nach Anspruch 6, dadurch gekennzeichnet, dass ein Betragsmaximum von g(x) bestimmt wird, dass ein Betragsschwellenwert bestimmt wird, welcher einen vorbestimmten Prozentsatz geringer als das Betragsmaximum ist, dass ein zusammenhängender Bereich von g(x) bestimmt wird, in welchem das Betragsmaximum liegt und in welchem der Betrag von g(x) nicht kleiner als der Betragsschwellenwert liegt, dass in diesem Bereich eine Gerade an den Phasenverlauf von g(x) angepasst wird, dass eine Steigung S der Gerade bestimmt wird, und dass die Gradientenmomentdifferenz GMD wie folgt bestimmt wird:
wobei k_max das maximale Gradientenmoment ist.
Verfahren nach einem der vorherigen Ansprüche, dadurch gekennzeichnet, dass eine Kalibrierungsmessung entlang der x-Achse (x) durchgeführt wird, wobei eine Gradientenmomentdifferenz GMD_x zwischen einem angenommenen und einem tatsächlich abgetasteten Abtastpunkt entlang der x-Achse (x) bestimmt wird, dass eine Kalibrierungsmessung entlang der y-Achse (y) durchgeführt wird, wobei eine Gradientenmomentdifferenz GMD_y zwischen einem angenommenen und einem tatsächlich abgetasteten Abtastpunkt entlang der y-Achse (y) bestimmt wird, dass die angenommene Position durch eine x-Koordinate k'_x und eine y-Koordinate k'_y definiert wird, dass die korrigierte Position durch eine x-Koordinate k_x und eine y-Koordinate k_y definiert wird, dass die x-Koordinate k_x der korrigierten Position gemäß der folgenden Gleichung bestimmt wird
dass die y-Koordinate k_y der korrigierten Position gemäß der folgenden Gleichung bestimmt wird
wobei φ ein Winkel ist, welcher zwischen der jeweiligen Speiche (x; y; 1) und der x-Achse ausgebildet ist.
Verfahren nach einem der vorherigen Ansprüche, dadurch gekennzeichnet, dass vor dem Auslesen aller Speichen (x; y; 1) des Rohdatenraums (4) vorbereitende Messungen durchgeführt werden, um die Magnetisierung von Gradientenspulen der Magnetresonanzanlage (5), mit welchen die Magnetfeldgradienten erzeugt werden, in einen Gleichgewichtszustand zu bringen, und dass die mindestens eine Kalibrierungsmessung im Rahmen der vorbereitenden Messungen durchgeführt wird.
Verfahren nach einem der vorherigen Ansprüche, dadurch gekennzeichnet, dass eine Gradientenmomentdifferenz zwischen dem angenommenen und dem tatsächlich angelegten Gradientenmoment abgeschätzt wird, und dass ein Vordephasierungsgradient bei der mindestens einen Kalibrierungsmessung geschaltet wird, welcher derart an die geschätzte Gradientenmomentdifferenz angepasst ist, dass die geschätzte Gradientenmomentdifferenz durch den Vordephasierungsgradient im Wesentlichen ausgeregelt wird.
Magnetresonanzanlage zur Erstellung eines Bilddatensatzes mittels einer radialen Abtastung eines Rohdatenraums (4), wobei die Magnetresonanzanlage (5) eine Ansteuereinheit (11) zur Ansteuerung eines Tomographen (3), eine Empfangsvorrichtung (12) zum Empfang von von dem Tomographen (3) aufgenommenen Signalen und eine Auswertevorrichtung (13) zur Auswertung der Signale und Erstellung des Bilddatensatzes umfasst, wobei die Magnetresonanzanlage (5) derart ausgestaltet ist, dass die Magnetresonanzanlage (5) mindestens eine Kalibrierungsmessung für mindestens eine vorbestimmte Speiche (x; y) der radialen Abtastung durchführt, um eine Gradientenmomentdifferenz zwischen einem angenommenen und einem tatsächlich von dem Tomographen (3) angelegten Gradientenmoment entlang der mindestens einen Speiche (x; y) zu bestimmen, dass die Magnetresonanzanlage (5) alle Speichen (x, y, 1) des vorbestimmten Rohdatenraums (4) ausliest, indem die Magnetresonanzanlage (5) mehrere Magnetfeldgradienten in jeweils einer Raumrichtung schaltet, um jeweils Abtastpunkte einer jeweiligen Speiche (x; y; 1) auszulesen, und dass die Magnetresonanzanlage (5) eine Position jedes Abtastpunkts jeder Speiche (x; y; 1) abhängig von der Gradientenmomentdifferenz korrigiert, indem die Magnetresonanzanlage (5) eine aufgrund der jeweils geschalteten Magnetfeldgradienten angenommene Position des jeweiligen Abtastpunkts um die Gradientenmomentdifferenz verschiebt.
Magnetresonanzanlage nach Anspruch 12, dadurch gekennzeichnet, dass die Magnetresonanzanlage (5) zur Durchführung des Verfahrens nach einem der Ansprüche 1–11 ausgestaltet ist.
Computerprogrammprodukt, welches ein Programm umfasst und direkt in einen Speicher einer programmierbaren Steuereinrichtung (6) einer Magnetresonanzanlage (5) ladbar ist, mit Programm-Mitteln, um alle Schritte des Verfahrens nach einem der Ansprüche 1–11 auszuführen, wenn das Programm in der Steuereinrichtung (6) der Magnetresonanzanlage (5) ausgeführt wird.
Elektronisch lesbarer Datenträger mit darauf gespeicherten elektronisch lesbaren Steuerinformationen, welche derart ausgestaltet sind, dass sie bei Verwendung des Datenträgers (14) in einer Steuereinrichtung (6) einer Magnetresonanzanlage (5) das Verfahren nach einem der Ansprüche 1–11 durchführen.