DE19954925A1

DE19954925A1 - Verfahren zum Korrigieren von Feldinhomogenitäten höherer Ordnung in einer Apparatur der magnetischen Resonanz

Info

Publication number: DE19954925A1
Application number: DE19954925A
Authority: DE
Inventors: Arno Nauerth
Original assignee: Bruker Medizintechnik GmbH
Current assignee: Bruker Biospin MRI GmbH
Priority date: 1999-11-16
Filing date: 1999-11-16
Publication date: 2001-05-23
Anticipated expiration: 2019-11-17
Also published as: GB2357852B; GB0027995D0; US6392412B1; GB2357852A; DE19954925C2

Abstract

Die Erfindung betrifft ein Verfahren zum Korrigieren von Feldinhomogenitäten höherer Ordnung in einer Apparatur der magnetischen Resonanz. Nach einem Anregungspuls und dem Einwirken eines Phasengradienten in einer vorgegebenen Richtung p wird nach zwei fest vorgegebenen Zeiten t¶d1¶ bzw. t¶d2¶ ein Datenpunkt genommen. Dies wird für systematisch veränderte Stärken des Phasengradienten und gleichbleibenden t¶d1¶ bzw. t¶d2¶ wiederholt und daraus zwei Sätze des Signalverlaufs in p-Richtung gewonnen. Die beiden Phasenverläufe der beiden Fourier-transformierten Datensätze in p-Richtung werden subtrahiert. Das Profil der Differenzphase in p-Richtung entspricht einem Feldfehlerprofil. Das Verfahren wird für mehrere p-Richtungen wiederholt und daraus Koeffizienten höherer Ordnung einer Reihenentwicklung der Feldinhomogenitäten nach Kugelflächenfunktionen gewonnen.

Description

Die vorliegende Erfindung steht in Zusammenhang mit der zeitgleich eingereichten Erfindung "Verfahren zum Korrigieren linearer Feldinho mogenitäten in einer Apparatur der magnetischen Resonanz" (Anwalts zeichen P7439), auf die vollinhaltlich bezug genommen wird. Die bei den Verfahren können sich in vorteilhafter Weise ergänzen. Insbeson dere kann das zitierte Verfahren dem vorliegenden vorangestellt wer den, um lineare Feldinhomogenitäten bereits vorab zu beseitigen.

Die vorliegende Erfindung betrifft ein Verfahren zum Bestimmen und Korrigieren von Magnetfeldinhomogenitäten höherer Ordnung eines nahezu homogenen Magnetfelds B₀ in Untersuchungsvolumen einer Apparatur der magnetischen Resonanz wobei in einer im Untersu chungsvolumen angeordnete Probe durch einen Hochfrequenzpuls magnetische Resonanz angeregt wird, mindestens ein zusätzliches, lineares, magnetisches Gradientenfeld angelegt und ein magnetisches Resonanzsignal gemessen wird.

Ein solches Verfahren ist beispielsweise für eine Apparatur der Kern spintomographie bekannt aus US 5,345178.

Bei dem bekannten Verfahren wirkt nach einer Hochfrequenzanregung auf eine Probe eine Gradienten- oder Spin-EchoSequenz gemessen. Das NMR-Signal wird daraufhin Fourier-transformiert und in einem vor gegebenen Bereich ein Phasenprofil ermittelt. Dies wird für eine Viel zahl von Projektionsrichtungen wiederholt und die gewonnenen Pha senkurven werden mit einer Fit-Methode analysiert und daraus Koeffizi enten einer Reihenentwicklung des Magnetfeldverlaufs nach Kugelflä chenfunktionen bestimmt. Daraus werden wiederum Korrekturströme für einzelne, den Koeffizienten zugeordnete Shim- bzw. Gradienten spulen berechnet.

Das bekannte Verfahren hat u. a. die Nachteile, daß die chemische Ver schiebung zwischen Wasser- und Fettanteilen des Signals eine uner wünschte geometrische Verschiebungen aufprägt und bei Proben mit kurzen T₁/T*₂-Zeiten nur schlecht oder gar nicht anwendbar ist.

Es besteht daher der Bedarf nach einem schnellen, unkomplizierten, direkten Shimverfahren, das unempfindlich gegen Einflüsse der chemi schen Verschiebungen gemacht werden kann und auch bei Proben mit kurzen Relaxationszeiten erfolgreich angewendet werden kann.

Die Aufgabe wird durch ein eingangs genanntes Verfahren gelöst mit den folgenden Schritten:

A) Auf die Probe wird ein erster Anregungshochfrequenzpuls eingestrahlt;
B) Ein erster Phasengradient G_ix1 in einer vorgegebenen Richtung x wird angelegt;
C) Zu einem festen Zeitpunkt t_dx1 nach dem ersten Anre gungshochfrequenzpuls wird ein Wert S_ix1 des von der Probe ausgehenden magnetischen Resonanzsignals ge messen, digitalisiert und abgespeichert;
D) Auf die Probe wird ein zweiter Anregungshochfrequenz puls eingestrahlt;
E) Ein zweiter Phasengradient G_ix2 in der vorgegebenen Richtung x wird angelegt;
F) Zu einem festen Zeitpunkt t_dx2 nach dem zweiten Anre gungshochfrequenzpuls wird ein Wert S_ix2 des von der Probe ausgehenden magnetischen Resonanzsignals ge messen, digitalisiert und abgespeichert, wobei G_ix2 und t_dx2 so gewählt sind, daß die integrale
gleich sind;
G) Die Schritte A) bis F) werden mehrmals mit systematisch geänderter Stärke des Phasengradienten G_ix wiederholt;
H) Die Werte der gemessenen Resonanzsignale S_ix1 werden in Abhängigkeit von der zugehörigen Gradientenstärke G_ix zu einem Quasi-Spin-Echo-Datensatz S_x1 zusammenge faßt;
I) Die Werte der gemessenen Resonanzsignale S_ix2 werden in Abhängigkeit von der zugehörigen Gradientenstärke G_ix zu einem Quasi-Spin-Echo-Datensatz S_x2 zusammenge faßt;
J) Der Datensatz S_x1 wird Fourier-transformiert und ggf. pha senkorrigiert, so daß die Phasen ϕ_ix1 der phasenkorrigier ten Resonanzsignale S'_ix1 alle im wesentlichen denselben Wert haben;
K) Der Datensatz S_x2 wird ebenfalls Fourier-transformiert und mit denselben Parametern phasenkorrigiert wie der Da tensatz S_x1 in Schritt J), wodurch die Differenz der Phasen ϕ'_ix1-ϕ'_ix2 der phasenkorrigierten Phasen der Signale S'_ix1 und S'_ix2 ein Maß für ein Profil der Magneffeldinhomoge nität entlang der Richtung x darstellen;
L) In folgenden Messungen der magnetischen Resonanz in der Apparatur wird ein Korrekturmagnetfeld B(x) zur Ho mogenisierung des Magnetfeldes B₀ im Untersuchungs volumen angelegt, das die in Schritt K) ermittelte Magnet feldinhomogenität kompensiert.

Nach jeder Anregung wird nicht ein vollständiges Signalecho bzw. FID gemessen sondern es wird bei zwei fest vorgegebenen Zeiten nach der Anregung jeweils ein Meßpunkt aufgenommen. Die Zeiten t_d1 und t_d2 sind immer dieselben, daher auch die Dephasierung aufgrund der In homogenitäten des Magnetfelds R₀. Zusätzlich wird durch Anlegen ei nes Phasengradienten im Intervall zwischen t = 0 und t = t_d1 bzw. zwi schen t = 0 und t = t_d2 jeweils eine zusätzliche, definiert einstellbare De phasierung aufgeprägt. Zu den Zeitpunkten der Datenaufnahme über lagern sich die beiden Effekte. Die Gradientenstärken sind so gewählt, daß es zu jedem Meßwert bei t_d1 einen entsprechenden bei t_d2 gibt, der dieselbe Dephasierung aufgrund des angelegten Phasengradienten aufweist, nicht aber dieselbe aufgrund der Feldinhomogenität von B₀. Die Auswertung nach den oben angegebenen Schritten ergibt ein Pha sendifferenzprofil in Richtung des angelegten Phasengradienten, das, bis auf T₂-Relaxationseffekte, allein auf der Wirkung der B₀- Inhomogenität in dieser Richtung beruht. Die in Schritt J) angesproche ne Phasenkorrektur ist optional, wichtig ist nur, daß auf beide Datensät ze die gleiche Phasenkorrektur angewendet wird.

Die Anzahl der notwendigen Anregungen kann verringert werden, wenn die Zeiten t_d1 und t_d2 und die Gradientenstärken in einem rationalen Verhältnis stehen. Man kann dann zumindest für einige Gradientenstär ken des Phasengradienten sowohl bei t_d1 als auch bei t_d1 einen Meß punkt nehmen, wobei der bei t_d2 einer Messung bei t_d1 mit größerer Stärke des Gradienten zugeordnet wird, so daß für dieses Paar die Be dingung aus Schritt F) eingehalten wird.

Der Einfluß der chemischen Verschiebung ist dadurch eliminiert, daß ein "Quasi-Spin-Echo" in Phasenrichtung und nicht in Leserichtung er zeugt und ausgewertet wird, d. h. die Meßpunkte liegen jeweils beim selben relativen Zeitpunkt (t_d1, t_d2) nach der Anregung und unterschei den sich jeweils nicht bezüglich der Dephasierung durch die chemische Verschiebung. Es empfiehlt sich, die Zeiten t_d1, t_d2 so zu wählen, daß zumindest näherungsweise Fett- und Wasseranteile der Signale im B₀- Feld in Phase sind.

Die Zeit t_d kann weitgehend frei gewählt werden, insbesondere kann sie für Proben mit kurzen Relaxationszeiten T₁/T*₂ auch sehr kurz sein.

Das zunächst für eine Richtung beschriebene Verfahren läßt sich zwanglos auf mehr Richtungen übertragen. Eine Kompensation in zwei Dimensionen bietet sich insbesondere bei der Untersuchung von Scheiben aus einem Objekt an. Dabei kann das Feld bei Mehrfach scheibenuntersuchungen für jede Scheibe separat homogenisiert wer den.

Selbstverständlich kann das Verfahren entsprechend für dreidimensio nale Volumina erweitert werden, indem es analog für eine weitere Di mension durchgeführt wird. Damit kann das Feld im gesamten Proben volumen homogenisiert werden oder in Verbindung mit volumenselekti ven Maßnahmen auch für ausgewählte Teilvolumina, ggf. viele unter schiedliche innerhalb eines Objekts.

Besonders bevorzugt wird das Verfahren analog für eine Vielzahl vor gegebener Richtungen durchgeführt wird, die so ausgewählt werden, daß sich daraus nach an sich bekannten Verfahren die wesentlichen Entwicklungskoeffizienten einer Reihenentwicklung des Magnetfeldes B₀ nach Kugelflächenfunktionen bestimmen lassen, insbesondere die jenigen zweiter Ordnung. Ein Shimverfahren mit dem Namen "Fastmap" ist beispielsweise in dem Artikel von R. Gruetter in Magnetic Re sonance in Medicine 29: 804-811 (1993) beschrieben. In einer bevor zugten Variante wird zur Bestimmung der Koeffizienten bis einschließ lich zweiter Ordnung entlang von nur sechs Projektionsrichtungen ge messen. Diese Richtungen können für das vorliegende, erfindungsge mäße Verfahren übernommen werden.

Je nach Ausstattung der Apparatur mit Shim- und Gradientenspulen und nach Art der Aufgabenstellung wird das erfindungsgemäße Verfah ren vorzugsweise analog für eine Vielzahl vorgegebener Richtungen durchgeführt, die so ausgewählt werden, daß sich daraus nach an sich bekannten Verfahren die Korrekturströme durch Shimspulen oder Gra dientenspulen vorgegebener Korrekturfeldprofile bestimmen lassen. Diese Feldprofile können - müssen aber nicht - mit den Profilen über einstimmen, die durch die entsprechenden reinen Kugelflächenfunktio nen vorgegeben sind.

Bevorzugt handelt es sich bei dem Verfahren um ein Verfahren der kernmagnetischen Resonanz, insbesondere ist es in ein Verfahren der bildgebenden magnetischen Resonanz integriert. Bildgebende Kernre sonanzapparaturen, aber auch immer mehr solche der hochauflösen den Kernresonanzspektroskopie, sind standardmäßig mit Gradienten spulen und Shimsystemen ausgestattet. Insbesondere in den bildge benden Apparaturen, kann es ohne Hardwareprobleme in die beste hende Software integriert werden, insbesondere kann es der eigentli chen Messung vorgeschaltet werden, aber auch in diese vollständig integriert sein, so daß das Feld während eines Untersuchungspro gramms neu homogenisiert wird. Wie weiter oben bereits ausgeführt, kann die Homogenisierung einzeln für getrennte Untersuchungsteilbe reiche erfolgen.

In einer Ausführungsform der Erfindung wird das Verfahren in zwei Di mensionen durchgeführt und es geht ihm, ggf. inklusive Schritt A) ein Scheibenselektionsschritt voran, der im Untersuchungsobjekt eine Scheibe senkrecht zu einer Richtung selektiert.

Alternativ wird ihm ggf. inklusive Schritt A) ein Volumenselektionsschritt vorangehen, der im Untersuchungsobjekt ausgewähltes Volumen se lektiert.

Bei Mehrfachvolumenexperimenten bzw. Mehrfachscheibenexperi menten werden mehrere Sätze von Korrekturfeldprofilen ermittelt, die unterschiedlichen Volumina bzw. Scheiben zugeordnet sind.

Aus den ermittelten Profilen können Korrekturströme zum Einspeisen in Shimspulen der Apparatur ermittelt werden.

Alternativ oder auch zusätzlich können aus den ermittelten Profilen Offsetströme zum Einspeisen in Gradientenspulen der Apparatur er mittelt werden. Dadurch lassen sich diese Ströme besonders einfach auch während eines Pulsprogramms verändern. Anstelle von Shim spulen können grundsätzlich auch andere Mittel zur Feldhomogenisie rung eingesetzt werden, z. B. ferromagnetische oder permanentmagne tische Shimplättchen, deren Anzahl und/oder Position zu bestimmen ist. In einer Ausführungsform umfaßt die Probe biologisches Gewebe. Die obengenannten Vorteile der Erfindung wirken sich besonders positiv aus, bei nicht homogenen Proben mit Fett- und Wasseranteil, deren Position ggf. nicht exakt bekannt ist. Dies trifft auf biologische Proben wie Versuchstier oder auch menschliche Patienten zu. Eine nicht exakt zentrierte Probe erzeugt bei dem erfindungsgemäßen Verfahren zu nächst bei der Bestimmung von Koeffizienten einer bestimmten Ord nung Meßfehler niedrigerer Ordnungen, die ihrerseits aber wieder durch weitere Iterationen, die gerade diese Ordnung erfassen, eliminiert wer den können.

Weitere Vorteile ergeben sich aus der Zeichnung und der Beschrei bung. Ebenso können die vorstehend genannten und die noch weiter aufgeführten Merkmale erfindungsgemäß einzeln für sich und zu meh reren in beliebigen Kombinationen Verwendung finden. Die gezeigten und beschriebenen Ausführungsformen sind nicht als abschließende Aufzählung zu verstehen, sondern haben vielmehr beispielhaften Cha rakter für die Schilderung der Erfindung.

Die Erfindung ist in der Zeichnung dargestellt und wird anhand von Ausführungsbeispielen näher erläutert. Es zeigen:

Fig. 1: schematisch HF-Anregungen, Phasengradienten und Meßsi gnale für eine vorgegebene Richtung p.

Fig. 1 zeigt in der ersten Zeile schematisch den zeitlichen Verlauf eines angeregten NMR-Signals nach dem Zeitpunkt eines Hochfrequenzan regungspulses (t = 0), der hier jeweils durch einen vertikalen, schwar zen Balken symbolisiert ist. Die sich jeweils anschließenden vier Kurven sollen die von der Stärke des jeweiligen Phasengradienten abhängigen Verläufe der Signaleinhüllenden nach der Anregung darstellen. Sie sind noch jeweils durch die überlagerten, zu bestimmenden Feldinhomoge nitäten des Magnetfelds B₀ leicht unterschiedlich modifiziert (nicht dar gestellt). Zu den Zeitpunkten t_d1 bzw. t_d1 erfolgt jeweils die Messung. In Fig. 1 sind die Zeiten t_d1 und t_d2 voneinander verschieden jedoch jeweils für alle unterschiedlichen Stärken der in Richtung p angelegten Pha sengradienten G_p1 bzw. G_p2 (siehe zweite Zeile) gleich. Für eine geän derte Richtung p des Phasengradienten können die Zeiten t_d1 bzw. t_d2 ungeändert bleiben, dies ist jedoch nicht zwingend.

In der zweiten Zeile der Fig. 1 ist jeweils der zeitliche Verlauf der Pha sengradienten G_p1 bzw. G_p2 dargestellt, deren Richtungen p überein stimmen. Im Beispiel wird der Phasengradient erst nach dem HF-Puls eingeschaltet und steht während der Datenaufnahme an. Dies ist nicht zwingend. Wichtig ist nur, daß das Integral der Gradientenstärke im Intervall zwischen HF-Puls (t = 0) und Datenaufnahme t = t_d1 bzw. t = t_d2 systematisch variiert wird und daß für jede einmal eingestellte Richtung p t_d1 und t_d2 jeweils konstant bleibt. Wichtig ist auch, daß jeweils für ein Meßpaar mit Detektion bei t_d1 und Detektion bei t_d2 das zeitliche Integral über die Gradientenstärke von t = 0 bis t = t_d1 bzw. von t = 0 bis t = t_d2 gleich ist, d. h. die nach der Anregung aufgrund des Phasengradienten G_p1 bzw. G_p2 erzeugte Dephasierung übereinstimmt (nicht jedoch die aufgrund anderer Feldinhomogenitäten, da t_d1 und t_d2 verschieden sind). Im einfachsten Fall wirkt der Phasengradient im gesamten Inter vall von t = 0 bis t = t_d1 bzw. bis t = t_d2 mit konstanter Stärke G_p1 bzw. G_p2. Dann gilt für ein Meßpaar G_p1*t_d1 = G_p2* t_d2.

In der dritten Zeile ist für die verwendete Richtung p des Phasengra dienten (G_p1, G_p2) der Verlauf der Signalphase ϕ₁ bzw. ϕ₂ in Phasen richtung p aufgetragen, d. h. als Funktion der den Phasengradienten (G_p1 bzw. G_p2) entsprechenden Raumrichtung D_p. In Spalte 1 (t_d1, G_p1, ϕ₁wurde dabei zur besseren Veranschaulichung das Fourier transformierte Signal phasenkorrigiert, d. h. die Phase ist jetzt konstant. In Spalte 2 (t_d2, G_p2, ϕ₂) wurde dieselbe Phasenkorrektur wie in Spalte 1 verwendet. In Spalte 3 ist schließlich die Differenz der beiden Phasen verläufe dargestellt. Sie ist ein direktes Maß für das durch die Feldin homogenitäten des Grundmagnetfelds B₀ erzeugte Feldfehlerprofil ent lang der gemessenen Projektion.

Das Verfahren nach Fig. 1 wird nun für eine weitere Richtung p erneut durchgeführt und ein Phasendifferenz- und damit Feldfehlerprofil auch in dieser Richtung gewonnen. Die Richtungen p können nach an sich bekannten Kriterien, z. B. wie in dem als "Fastmap" bekannt geworde nen Verfahren beschrieben, ausgewählt und daraus Koeffizienten einer Reihenentwicklung des Feldes und/oder Korrekturströme durch Shim- und/oder Gradientenspulen gewonnen werden.

Es versteht sich, daß das Verfahren mit vielen in der magnetischen Re sonanz üblichen Pulssequenzen kombiniert bzw. in diese integriert werden kann.

Claims

1. Verfahren zum Bestimmen und Korrigieren von Magnetfeldinho mogenitäten höherer Ordnung eines nahezu homogenen Ma gnetfelds B₀ in Untersuchungsvolumen einer Apparatur der ma gnetischen Resonanz wobei in einer im Untersuchungsvolumen angeordnete Probe durch einen Hochfrequenzpuls magnetische Resonanz angeregt wird, mindestens ein zusätzliches, lineares magnetisches Gradientenfeld angelegt und ein magnetisches Resonanzsignal gemessen wird, gekennzeichnet durch folgende Schritte:

A) Auf die Probe wird ein erster Anregungshochfrequenzpuls eingestrahlt;
B) Ein erster Phasengradient G_ix1 in einer vorgegebenen Richtung x wird angelegt;
C) Zu einem festen Zeitpunkt t_dx1 nach dem ersten Anre gungshochfrequenzpuls wird ein Wert S_ix1 des von der Probe ausgehenden magnetischen Resonanzsignals ge messen, digitalisiert und abgespeichert;
D) Auf die Probe wird ein zweiter Anregungshochfrequenz puls eingestrahlt;
E) Ein zweiter Phasengradient G_ix2 in der vorgegebenen Richtung x wird angelegt;
F) Zu einem festen Zeitpunkt t_dx2 nach dem zweiten Anre gungshochfrequenzpuls wird ein Wert S_ix2 des von der Probe ausgehenden magnetischen Resonanzsignals ge messen, digitalisiert und abgespeichert, wobei G_ix2 und t_dx2 so gewählt sind, daß die Integrale
gleich sind;
G) Die Schritte A) bis F) werden mehrmals mit systematisch geänderter Stärke des Phasengradienten G_ix wiederholt;
H) Die Werte der gemessenen Resonanzsignale S_ix1 werden in Abhängigkeit von der zugehörigen Gradientenstärke G_ix zu einem Quasi-Spin-Echo-Datensatz S_x1 zusammenge faßt;
I) Die Werte der gemessenen Resonanzsignale Six2 werden in Abhängigkeit von der zugehörigen Gradientenstärke G_ix zu einem Quasi-Spin-Echo-Datensatz S_x2 zusammenge faßt;
J) Der Datensatz S_x1 wird Fourier-transformiert und ggf. pha senkorrigiert, so daß die Phasen ϕ_ix1 der phasenkorrigier ten Resonanzsignale S'_ix1 alle im wesentlichen denselben Wert haben;
K) Der Datensatz S_X2 wird ebenfalls Fourier-transformiert und mit denselben Parametern phasenkorrigiert wie der Da tensatz S_x1 in Schritt J), wodurch die Differenz der Phasen ϕ'_ix1-ϕ'_ix2der phasenkorrigierten Phasen der Signale S'_ix1 und S'_ix2 ein Maß für ein Profil der Magnetfeldinhomoge nität entlang der Richtung x darstellen;
L) In folgenden Messungen der magnetischen Resonanz in der Apparatur wird ein Korrekturmagnetfeld B(x) zur Ho mogenisierung des Magnetfeldes B₀ im Untersuchungs volumen angelegt, das die in Schritt K) ermittelte Magnet feldinhomogenität kompensiert.

2. Verfahren nach Anspruch 1, dadurch gekennzeichnet, daß es analog für mindestens eine weitere vorgegebene Richtung y durchgeführt wird, die vorzugsweise senkrecht auf der Richtung x steht.

3. Verfahren nach Anspruch 1, dadurch gekennzeichnet, daß es analog für eine Vielzahl vorgegebener Richtungen durchgeführt wird, die so ausgewählt werden, daß sich daraus nach an sich bekannten Verfahren die wesentlichen Entwicklungskoeffizienten einer Reihenentwicklung des Magnetfeldes B₀ nach Kugelflä chenfunktionen, insbesondere der zweiter Ordnung, bestimmen lassen.

4. Verfahren nach einem der vorangehenden Ansprüche, dadurch gekennzeichnet, daß es analog für eine Vielzahl vorgegebener Richtungen durchgeführt wird, die so ausgewählt werden, daß sich daraus nach an sich bekannten Verfahren die Korrektur ströme durch Shim- bzw. Gradientenspulen oder die Parameter anderer Homogenisierungsmittel mit jeweils vorgegebenen Kor rekturfeldprofilen bestimmen lassen.

5. Verfahren nach einem der vorangehenden Ansprüche, dadurch gekennzeichnet, daß es sich um ein Verfahren der kernmagneti schen Resonanz handelt.

6. Verfahren nach einem der vorangehenden Ansprüche, dadurch gekennzeichnet, daß es in ein Verfahren der bildgebenden ma gnetischen Resonanz integriert ist.

7. Verfahren nach Anspruch 6, dadurch gekennzeichnet, daß das Verfahren innerhalb einer Ebene (z. B. x, y) durchgeführt wird und ihm, ggf. inklusive Schritt A) ein Scheibenselektionsschritt voran geht, der im Untersuchungsobjekt eine Scheibe (z. B. x, y) senk recht zu einer Richtung (z. B. z) selektiert.

8. Verfahren nach Anspruch 6, dadurch gekennzeichnet, daß ihm ggf. inklusive Schritt A) ein Volumenselektionsschritt vorangeht, der im Untersuchungsobjekt ausgewähltes Volumen selektiert.

9. Verfahren nach Anspruch 8, dadurch gekennzeichnet, daß für Mehrfachvolumenexperimente mehrere Sätze von Korrekturgra dientenfeldern ermittelt werden, die unterschiedlichen Volumina zugeordnet sind.

10. Verfahren nach einem der vorangehenden Ansprüche, dadurch gekennzeichnet, daß aus den ermittelten Korrekturgradienten magnetfeldern Korrekturströme zum Einspeisen in Shimspulen der Apparatur ermittelt werden.

11. Verfahren nach einem der vorangehenden Ansprüche, dadurch gekennzeichnet, daß aus den ermittelten Korrekturgradienten magnetfeldern Offsetströme zum Einspeisen in Gradientenspu len der Apparatur ermittelt werden.

12. Verfahren nach einem der vorangehenden Ansprüche, dadurch gekennzeichnet, daß aus den ermittelten Korrekturgradienten magnetfeldern Anzahl und/der Positionen von ferromagnetischen oder permanentmagnetischen Shimplättchen in der Apparatur ermittelt werden.

13. Verfahren nach einem der vorangehenden Ansprüche, dadurch gekennzeichnet, daß die Probe biologisches Gewebe umfaßt.

14. Verfahren nach einem der vorangehenden Ansprüche, dadurch gekennzeichnet, daß die Probe sehr kurze Relaxationszeiten T₁ und/oder T*₂ aufweist.