DE10036207A1 - Verfahren zum Durchführen einer Perfusionsmessung mittels Magnetresonanzbildgebung - Google Patents

Abstract

Ein Verfahren zum Durchführen einer Perfusionsmessung mittels Magnetresonanzbildgebung beinhaltet folgende Merkmale: DOLLAR A - In einer zeitlichen Abfolge werden Bilddatensätze eines in einem Abbildungsvolumen eines Magnetresonanzgeräts positionierten abzubildenden Bereichs eines Untersuchungsobjekts erzeugt, DOLLAR A - während der zeitlichen Abfolge auftretende Lageveränderungen des abzubildenden Bereichs bezüglich des Abbildungsvolumens werden erfasst, und DOLLAR A - entsprechend den erfassten Lageveränderungen erfolgt ein Korrigieren von Auswirkungen der Lageveränderungen auf die Bilddatensätze.

Description

Die Erfindung betrifft ein Verfahren zum Durchführen einer Perfusionsmessung mittels Magnetresonanzbildgebung.

Die Magnetresonanztechnik ist eine bekannte Technik zum Ge­ winnen von Bildern eines Körperinneren eines Untersuchungsob­ jekts. Dabei werden in einem Magnetresonanzgerät einem stati­ schen Grundmagnetfeld schnell geschaltete Gradientenfelder überlagert, die von einem Gradientensystem erzeugt werden. Ferner umfasst das Magnetresonanzgerät ein Hochfrequenzsys­ tem, das zum Auslösen von Magnetresonanzsignalen Hochfre­ quenzsignale in das Untersuchungsobjekt einstrahlt und die erzeugten Magnetresonanzsignale aufnimmt, auf deren Basis Bilddatensätze und Magnetresonanzbilder erstellt werden.

Bei einer Ausführungsform einer funktionellen Magnetresonanz­ bildgebung werden mit einer gleichen Ortskodierung von einem abzubildenden Bereich eines Untersuchungsobjekts in einer zeitlichen Abfolge Bilddatensätze erzeugt. Daran anschließend wird eine retroperspektive Bewegungskorrektur der Bilddaten­ sätze durchgeführt, mit der Unterschiede zwischen den Bildda­ tensätzen, die Folge einer Lageveränderung des abzubildenden Bereichs bezüglich des Geräts während der zeitlichen Abfolge sind, ermittelbar und korrigierbar sind. Ein Verfahren zum Ermitteln von Lageveränderung aus zeitlich aufeinanderfolgend aufgenommenen Bilddatensätzen basiert dabei auf einer Be­ schreibung einer beliebigen Starrkörperbewegung im dreidimen­ sionalen Raum mittels sechs Bewegungsparametern, wobei drei Parameter Translationen und drei Parameter Rotationen kenn­ zeichnen. Vorgenannte Parameter werden beispielsweise in ei­ nem Spaltenvektor notiert. Die Werte aller Voxel oder aus­ gewählter Voxel eines ersten Bilddatensatzes und eines zwei­ ten Bilddatensatzes, der dem ersten zeitlich nachfolgend erzeugt worden ist, werden in einer übereinstimmenden Reihen­ folge in einem ersten Spaltenvektor und einem zweiten Spaltenvektor notiert. Zur Ermittlung einer Lageverände­ rung zwischen den Aufnahmezeitpunkten des ersten und des zweiten Bilddatensatzes, d. h. zur Bestimmung der Bewegungspa­ rameter wird nachfolgende Gleichung, die auf einer Taylor- Entwicklung erster Ordnung basiert, durch ein iteratives Ver­ fahren, beispielsweise ein Gauß-Newtonsches Iterationsverfah­ ren, gelöst:

Zur genaueren Beschreibung wird dazu beispielsweise auf das Buch von R. S. J. Frackowiak et al., "Human Brain Function", Academic Press, 1997, insbesondere Kapitel 3, Seiten 43 bis 58 und auf den Artikel von K. J. Friston et al., "Movement- Related Effects in fMRI Time-Series", Magnetic Resonance in Medicine 35 (1996), Seiten 346 bis 355 verwiesen.

Im letztgenannten Artikel ist darüber hinaus beschrieben, dass sich selbst bei einem optimalen Zurückdrehen bzw. Zu­ rückverschieben der Bilddatensätze hinsichtlich eines Refe­ renz-Bilddatensatzes nicht alle unerwünschten Signaldifferen­ zen infolge von Bewegungen beheben lassen. Ursächlich dafür ist, dass nach einer Lageveränderung des abzubildenden Be­ reichs auf bestimmte Volumenbereiche des abzubildenden Be­ reichs gegenüber ihrer Ausgangslage bei unveränderter Ortsko­ dierung Gradienten- und Hochfrequenzfelder anders wirken und sich damit Anrege-, Resonanz- und Relaxationseigenschaften der Volumenbereiche ändern. Dadurch wird das Signalverhalten dieser Volumenbereiche nicht nur für einen unmittelbar nach­ folgend aufgenommenen Bilddatensatz, sondern nachhaltig auch für weitere aufzunehmende Bilddatensätze verändert. In dem Artikel von K. J. Friston et al. ist dazu ein Näherungsverfah­ ren vorgeschlagen, mit dem im Nachgang an das Erzeugen von Bilddatensätzen auch letztgenannte, bewegungsverursachte Sig­ naldifferenzen aus Bilddatensätzen ausfilterbar sind.

Bei einem weiteren Verfahren zur bilddatensätzebasierten La­ geveränderungserfassung werden alle oder bestimmte ausgewähl­ te Punkte eines ersten im k-Raum beschriebenen Bilddatensatz und eines zweiten Bilddatensatzes, der dem ersten zeitlich nachfolgend erzeugten worden ist, miteinander verglichen. Dabei beruht das Verfahren darauf, dass sich infolge einer Lageveränderung zwischen den Aufnahmezeitpunkten beider Bild­ datensätze bei einem Vergleich von innerhalb beider Bildda­ tensätzen gleich angeordneten Datenpunkten Translationen und/oder Rotationen des abzubildenden Bereichs in einer Ver­ änderung von Phase und/oder Betrag der Datenpunkte widerspie­ geln. Ausführungsformen vorgenannten Verfahrens sind bei­ spielsweise im Artikel von L. C. Maas et al., "Decoupled Auto­ mated Rotational and Translational Registration for Functio­ nal MRI Time Series Data: The DART Registration Algorithm", Magnetic Resonance in Medicine 37 (1997), Seiten 131 bis 139 sowie im Artikel von Q. Chen et al., "Symmetric Phase-Only Matched Filtering of Fourier-Mellin Transforms for Image Re­ gistration and Recognition", IEEE Transactions on Pattern Analysis and Machine Intelligence, Vol. 16, No. 12 (1994), Seiten 1156 bis 1168 näher beschrieben.

Ein anderer Ansatz zur Vermeidung unerwünschter, bewegungs­ verursachter Unterschiede bei einer funktionellen Magnetreso­ nanzbildgebung besteht darin, die Bilddatensätze nicht im Nachgang retroperspektiv zu korrigieren, sondern eine pro­ spektive Bewegungskorrektur während eines Ablaufs der funkti­ onellen Magnetresonanzbildgebung durchzuführen. Dazu werden von Bilddatensatz zu Bilddatensatz eventuelle Lageveränderun­ gen des abzubildenden Bereichs beispielsweise durch orbitale Navigatorechos erfasst und eine Ortskodierung während des Ablaufs entsprechend angepasst. Dabei ist ein orbitales Navi­ gatorecho ein Magnetresonanzsignal, das durch einen kreisför­ migen k-Raumpfad gekennzeichnet ist und von einer speziellen Navigatorsequenz erzeugt wird. Anhand von orbitalen Naviga­ torechos, die zu unterschiedlichen Zeitpunkten erzeugt wer­ den, sind Lageveränderungen ermittelbar. Dazu wird beispiels­ weise vor jedem Erzeugen eines Bilddatensatzes die Navigator­ sequenz ausgeführt und ein Navigatorecho aufgenommen, welches zur Bewegungskorrektur mit einem Referenz-Navigatorecho ver­ glichen wird. Vorgenanntes ist beispielsweise in dem Artikel von H. A. Ward et al., "Real-Time Prospective Correction of Complex Multiplanar Motion in fMRI", Proc. of ISMRM 7 (1999), Seite 270 näher beschrieben.

Bei einem weiteren Verfahren werden Lageveränderungen des abzubildenden Bereichs optisch erfasst, indem beispielsweise am abzubildenden Bereich optische Reflektoren angebracht sind, die von einem optischen Erfassungssystem hinsichtlich ihrer Position überwacht werden. Näheres ist dazu beispiels­ weise in dem Artikel von H. Eviatar et al., "Real Time Head Motion Correction for Functional MRI", Proc. of ISMRM 7 (1999), Seite 269 erläutert. Desweiteren wird dazu auf die US 5,828,770 und die US 5,923,417 hingewiesen.

Bei einem weiteren Verfahren zur prospektiven Bewegungskor­ rektur werden die im bereits eingangs zitierten Buch von R. S. J. Frackowiak und Artikel von K. J. Friston beschriebenen Verfahren zum Ermitteln von Lageveränderungen aus zeitlich aufeinanderfolgend aufgenommenen Bilddatensätzen eingesetzt. Näheres dazu ist im Artikel von S. Thesen et al., "Prospec­ tive Acquisition Correction for Head Motion with Image-based Tracking for Real-Time fMRI", Proc. of ISMRM 8 (2000), Seite 56 beschrieben.

Bei einer Perfusionsmessung mittels Magnetresonanztechnik werden von einem selben abzubildenden Bereich eines Untersu­ chungsobjekts, beispielsweise einem Hirn eines Patienten, in einer möglichst schnellen zeitlichen Abfolge Volumendatensät­ ze aufgenommen. Dies geschieht unabhängig davon, ob mit oder ohne einer Kontrastmittelgabe gearbeitet wird. Aus einer zeitlichen Änderung eines Werts eines Voxels, das innerhalb der aufgenommenen Volumendatensätze gleich positioniert ist, ist eine Aussage über eine lokale Perfusion gewinnbar. Findet während der Aufnahme der Volumendatensätze bei einer gleichen Ortskodierung eine Lageveränderung des abzubildenden Bereichs statt, so führt dies zu einer Translation und/oder Rotation der einzelnen Volumendatensätze gegeneinander. Dadurch ent­ stehen systematische Fehler in vorgenannten Voxelwert- Zeitreihen, die wiederum zu einer verfälschten Aussage über die lokale Perfusion führen.

Aufgabe der vorliegenden Erfindung ist es daher, ein verbes­ sertes Verfahren zum Durchführen einer Perfusionsmessung mit­ tels Magnetresonanzbildgebung zu schaffen, das unter anderem vorgenannten Nachteil des Standes der Technik bei Perfusions­ messungen verhindert.

Die Aufgabe wird erfindungsgemäß durch ein Verfahren zum Durchführen einer Perfusionsmessung mittels Magnetresonanz­ bildgebung gelöst, das folgende Merkmale beinhaltet:

• - In einer zeitlichen Abfolge werden Bilddatensätze eines in einem Abbildungsvolumen eines Magnetresonanzgeräts positio­ nierten abzubildenden Bereichs eines Untersuchungsobjekts erzeugt,
• - während der zeitlichen Abfolge auftretende Lageveränderun­ gen des abzubildenden Bereichs bezüglich des Abbildungsvo­ lumens werden erfasst, und
• - entsprechend den erfassten Lageveränderungen erfolgt ein Korrigieren von Auswirkungen der Lageveränderungen auf die Bilddatensätze.

Dadurch sind entsprechende Voxelwert-Zeitreihen der korri­ gierten Bilddatensätze frei von einem systematischen Fehler, so dass eine unverfälschte Aussage über eine lokale Perfusion gewinnbar ist.

In einer vorteilhaften Ausgestaltung erfolgt das Korrigieren zeitlich nach einem abgeschlossenen Erzeugen aller Bilddatensätze durch ein Korrigieren der Bilddatensätze. Dabei werden Bilddatensätze, für die hinsichtlich eines vorgebbaren Refe­ renz-Bilddatensatzes eine Lageveränderung erfasst wurde, ent­ sprechend der erfassten Lageveränderung zurückgedreht und/oder zurückverschoben.

In einer anderen vorteilhaften Ausgestaltung erfolgt das Kor­ rigieren während der zeitlichen Abfolge durch ein Anpassen einer Ortskodierung des Magnetresonanzgeräts von Bilddaten­ satz zu Bilddatensatz entsprechend den ermittelten Lageverän­ derungen. Dadurch ist eine retroperspektive Bewegungskorrek­ tur, verbunden mit einem Zurückdrehen und/oder Zurückver­ schieben von Bilddatensätzen, überflüssig.

In einer vorteilhaften Ausgestaltung wird dazu eine Shim- Einstellung des Magnetresonanzgeräts, umfassend Shim-Ströme eines aktiven Shim-Systems und Offset-Ströme eines Gradien­ tenspulensystems des Magnetresonanzgeräts, mit angepasst.

In einer vorteilhaften Ausgestaltung werden die Lageverände­ rungen optisch erfasst. Dazu sind Verfahren und Vorrichtungen entsprechend denen bekannt, die in der Beschreibungseinlei­ tung aufgeführt sind.

In einer vorteilhaften Ausgestaltung werden die Lageverände­ rungen durch orbitale Navigatorechos erfasst.

In einer vorteilhaften Ausgestaltung werden die Lageverände­ rungen aus zeitlich aufeinanderfolgend erzeugten Bilddaten­ sätzen ermittelt. Dies ist besonders vorteilhaft, weil dabei keine zusätzlichen Vorrichtungen zum Magnetresonanzgerät, wie bei den optischen Lageveränderungs-Erfassungsverfahren, benö­ tigt werden und auch keine zusätzlichen Pulssequenzen, wie bei einer Lageveränderungserfassung durch orbitale Navigator­ echos notwendig sind. Andererseits werden an Verfahren zur bilddatensätzebasierten Lageveränderungserfassung bei Perfu­ sionsmessungen infolge der starken zeitlichen Kontrastschwankungen von Bilddatensatz zu Bilddatensatz besondere Anforde­ rungen, insbesondere hinsichtlich einer Stabilität der Ver­ fahren gestellt. Dabei sind insbesondere die beiden nachfol­ gend als vorteilhafte Ausgestaltungen beschriebene Verfahren hinsichtlich der Kontrastschwankungen besonders robust und stabil. Ein Übertragen der aus der funktionellen Magnetreso­ nanzbildgebung bekannten bilddatensätzebasierten Verfahren zur Lageveränderungserfassung auf Perfusionsmessungen ist dabei nicht naheliegend, weil diese Verfahren im Allgemeinen keine oder nur sehr kleine Kontrastschwankungen zwischen zwei Bilddatensätzen tolerieren. Dahingegen steht ein Ermitteln von Kontrastschwankungen zwischen zwei Bilddatensätzen bei der Perfusionsmessung gerade im Vordergrund. Des weiteren scheitert die Lageveränderungserfassung bei einer Vielzahl der bilddatensätzebasierten Verfahren der funktionellen Mag­ netresonanzbildgebung bei Lageveränderungen, die größer weni­ ge Grad und/oder wenige Millimeter sind. Bei der funktionel­ len Magnetresonanzbildgebung stellen kleine Lageveränderungen von wenigen Grad und/oder bis ca. 1 mm die Hauptproblematik dar, wohingegen derartige kleine Lageveränderungen bei der Perfusionsmessung verglichen mit den wesentlich stärkeren Kontrastschwankungen nicht so kritisch sind. Bei der Perfusi­ onsmessung ist gerade ein Ermitteln vergleichsweise großer Lageveränderungen interessant, dies gilt beispielsweise ins­ besondere für eine Perfusionsmessung an einem Schlaganfalls­ patienten, bei dem zustandsbedingt mit größeren Bewegungen während der Messung zu rechnen ist.

In einer vorteilhaften Ausgestaltung wird das bereits ein­ gangs erläuterte Verfahren zum Ermitteln von Lageveränderun­ gen aus Bilddatensätzen auf Basis einer Taylor-Entwicklung erster Ordnung, insbesondere in Verbindung mit einem Gauß- Newtonschen Iterationsverfahren eingesetzt. Dieses aus der funktionellen Magnetresonanzbildgebung bekannte Verfahren ist bereits nach überraschend geringen Modifikationen in stabiler und robuster Weise auch bei Perfusionsmessungen einsetzbar. Dabei zählen zu den Modifikationen, dass bei der Perfusionsmessung gegenüber der funktionellen Magnetresonanzbildgebung bei Bilddatensätzen vergleichbarer Größe eine deutlich höhere Anzahl an Werten je Bilddatensatz ausgewählt wird und dass als ein Referenzbilddatensatz nicht ein in einer zeitlichen Mitte der Abfolge aufgenommener Bilddatensatz, sondern einer der zeitlich ersten Bilddatensätze herangezogen wird. Des weiteren werden die partiellen Ableitungen der Jacobischen Funktionalmatrix für die ausgewählten Werte eines Bilddaten­ satzes nach den Bewegungsparametern in Form einfacher Diffe­ renzquotienten mit einer linearen Interpolation ermittelt. Vorgenanntes Vorgehen weist in vorteilhafter Weise als eine intrinsische Eigenschaft auf, dass sich für alle Werte von Bilddatensatz zu Bilddatensatz nach einer Bewegungskorrektur ein "weicher" Übergang ergibt.

In einer vorteilhaften Ausgestaltung beinhaltet ein Verfahren zum Ermitteln der Lageveränderung aus zeitlich aufeinander­ folgend erzeugten Bilddatensätzen folgende Merkmale:

• - Wenigstens für ausgewählte Werte eines ersten Bilddatensat­ zes wird eine erste Häufigkeitsverteilung n(x) gebildet,
• - für entsprechend den ausgewählten Werten des ersten Bildda­ tensatzes ausgewählte Werte eines zweiten Bilddatensatzes, der dem ersten zeitlich nachfolgend erzeugt worden ist, wird eine zweite Häufigkeitsverteilung n(y) gebildet,
• - aus in den Bilddatensätzen entsprechend positionierten Wer­ ten werden Wertepaare gebildet, mit denen eine weitere Häu­ figkeitsverteilung n(x, y) gebildet wird,
• - aus der ersten, zweiten und weiteren Häufigkeitsverteilung wird eine Mutual Information
  gebildet,
• - einer der Bilddatensätze wird mit Parametern versehen, so dass entsprechend einer beliebigen Lageveränderung des ab­ zubildenden Bereichs im dreidimensionalen Raum der Bildda­ tensatz anpassbar ist, und
• - durch ein Optimierungsverfahren werden die Parameter so bestimmt, dass die Mutual Information minimal wird.

Vorgenanntes Verfahren ist selbst bei stärksten Kontrast­ schwankungen noch stabil. Dabei ist mit sechs Parametern die beliebige Lageveränderung im dreidimensionalen Raum be­ schreibbar, wobei mit drei der Parameter eine beliebige Translationsbewegung und mit den übrigen drei Parametern eine beliebige Rotationsbewegung beschreibbar ist. Für das Opti­ mierungsverfahren sind unter anderem das Downhill-Simplex-, das Powell's-, das Konjugierte-Gradienten- und/oder das Vari­ abel-Metrische-Verfahren einsetzbar, die beispielsweise im Buch von W. H. Press et al. "Numerical Recipes in C. The art of scientific computing", Cambridge Univ. Pr., 1992, Seiten 408 bis 430, beschrieben sind. Des weiteren wird für eine genauere Erläuterung hinsichtlich der Mutual Information auf die Promotionsarbeit von P. A. Viola "Alignment by Maximiza­ tion of Mutual Information", AI-TR1548 Massachusetts Insti­ tute of Technology, Dept. of Electrical Engineering and Com­ puter Science, Juni 1995, hingewiesen.

Weitere vorteilhafte Ausgestaltungen sind in den Unteransprü­ chen 10 bis 13 beschrieben. Weitere Vorteile, Merkmale und Einzelheiten ergeben sich aus dem im folgenden beschriebenen Ausführungsbeispiel anhand der Zeichnung.

Die Figur zeigt als ein Ausführungsbeispiel der Erfindung ein Flussdiagramm für einen Abschnitt einer Perfusionsmessung mittels Magnetresonanztechnik mit einer prospektiven Bewe­ gungskorrektur. Dazu werden in einer zeitlichen Abfolge von einem selben abzubildenden Bereich eines Untersuchungsobjekts Bilddatensätze aufgenommen. Dabei wird in einem ersten Schritt 10 ein Referenz-Bilddatensatz des abzubildenden Be­ reichs aufgenommen.

In einem sich anschließenden Schritt 20 wird wiederum ein Bilddatensatz des abzubildenden Bereichs aufgenommen. Vor einem Aufnehmen eines weiteren Bilddatensatzes in einer Wie­ derholung des Schritts 20 wird der Bilddatensatz bezüglich dem Referenz-Bilddatensatz in einem Schritt 30 auf eventuelle Lageveränderungen des abzubildenden Bereichs hin untersucht und in dem Fall, dass eine Lageveränderung ermittelt wird, wird vor einer Aufnahme des weiteren Bilddatensatzes in einem Schritt 40 eine Ortskodierung und eine Shim-Einstellung des Gradienten- und Shim-Systems entsprechend der ermittelten Lageveränderung nachgeführt. Dadurch liegen an einem Ende der Perfusionsmessung alle Bilddatensätze bereits bewegungskorri­ giert vor, so dass diese zum Bilden entsprechender Perfusi­ onsdarstellungen direkt voneinander subtrahierbar sind.

Im folgenden wird nun das bilddatensätzebasierte Ermitteln einer Lageveränderung des Schritts 30 näher betrachtet. Dabei wird der aufgenommene Bilddatensatz in einem Schritt 31 zu­ nächst mit einem Parametervektor versehen, mit dem der Bilddatensatz an eine beliebige Lageveränderung anpassbar ist. Der Parametervektor umfasst dabei sechs Parameter, mit denen die beliebige Lageveränderung in dreidimensionalem Raum beschreibbar ist, wobei mit drei der Parameter eine be­ liebige Translationsbewegung und mit den übrigen drei Parame­ tern eine beliebige Rotationsbewegung beschrieben wird.

In einem Schritt 32 wird eine Häufigkeitsverteilung n(x) für den Referenz-Bilddatensatz gebildet. Ausgehend von dreidimen­ sionalen Bilddatensätzen, beispielsweise bestehend aus 128 × 128 × 128 Voxeln, wobei jedem Voxel ein Grauwert auf einer Skala der natürlichen Zahlen von 0 bis 100 zugeordnet ist, wird dabei über den einzelnen Grauwerten eine Anzahl von Vo­ xeln des Referenz-Bilddatensatzes aufgetragen, die den ent­ sprechenden Grauwert aufweisen.

In einem Schritt 34 wird entsprechend dem Schritt 32 eine Häufigkeitsverteilung n(y()) für den mit dem Parametervek­ tor an beliebige Lageveränderungen anpassbaren Bilddaten­ satz durchgeführt. Dabei ergibt sich eine Abhängigkeit der Häufigkeitsverteilung n(y()) von dem Parametervektor .

In einem Schritt 33 wird schließlich eine dritte Häufigkeits­ verteilung n(x, y()) für Wertepaare des Referenz-Bilddaten­ satzes und des anpassbaren Bilddatensatzes gebildet. Dabei entstehen die Wertepaare durch einen ersten Grauwert eines Voxels des Referenz-Bilddatensatzes und einen zweiten Grau­ wert für ein Voxel des anpassbaren Bilddatensatzes, das hin­ sichtlich einer Anordnung innerhalb des Bilddatensatzes gleich einer Anordnung des Voxels des Referenz-Bilddaten­ satzes ist.

In einem letzten Unterschritt 35 des Schritts 30 wird schließlich der Parametervektor und damit die Lageverände­ rung ermittelt. Dies geschieht dadurch, dass die Mutual In­ formation MI minimiert wird. Dabei ist die Mutual Information MI wie folgt definiert:

Dabei wird zum Minimieren auf bekannte Optimierungsverfahren zurückgegriffen. Das im Schritt 30 beschriebene Verfahren zum bilddatensätzebasierten Ermitteln von Lageveränderungen ist für Perfusionsmessungen in besonderer Weise geeignet, weil auch bei großen Kontrastunterschieden zwischen Bilddatensät­ zen, wie sie bei Perfusionsmessungen auftreten, eine Lagever­ änderung des abzubildenden Bereichs zwischen Aufnahmezeit­ punkten der Bilddatensätze sicher detektierbar ist.

Claims (13)

1. Verfahren zum Durchführen einer Perfusionsmessung mittels Magnetresonanzbildgebung, beinhaltend folgende Merkmale:

• - In einer zeitlichen Abfolge werden Bilddatensätze eines in einem Abbildungsvolumen eines Magnetresonanzgeräts positio­ nierten abzubildenden Bereichs eines Untersuchungsobjekts erzeugt,
• - während der zeitlichen Abfolge auftretende Lageveränderun­ gen des abzubildenden Bereichs bezüglich des Abbildungsvo­ lumens werden erfasst, und
• - entsprechend den erfassten Lageveränderungen erfolgt ein Korrigieren von Auswirkungen der Lageveränderungen auf die Bilddatensätze.

2. Verfahren nach Anspruch 1, wobei das Korrigieren zeit­ lich nach einem abgeschlossenen Erzeugen aller Bilddatensätze durch ein Korrigieren der Bilddatensätze erfolgt.

3. Verfahren nach Anspruch 1, wobei das Korrigieren wäh­ rend der zeitlichen Abfolge durch ein Anpassen einer Ortsko­ dierung des Magnetresonanzgeräts von Bilddatensatz zu Bildda­ tensatz entsprechend den ermittelten Lageveränderungen er­ folgt.

4. Verfahren nach Anspruch 3, wobei eine Shim-Einstellung des Magnetresonanzgeräts, umfassend Shim-Ströme eines aktiven Shim-Systems und Offset-Ströme eines Gradientenspulensystems des Magnetresonanzgeräts, mit angepasst wird.

5. Verfahren nach einem der Ansprüche 1 bis 4, wobei die Lageveränderungen optisch erfasst werden.

6. Verfahren nach einem der Ansprüche 1 bis 5, wobei die Lageveränderungen durch orbitale Navigatorechos erfasst wer­ den.

7. Verfahren nach einem der Ansprüche 1 bis 6, wobei die Lageveränderungen aus zeitlich aufeinanderfolgend erzeugten Bilddatensätzen ermittelt werden.

8. Verfahren nach Anspruch 7, beinhaltend folgende Merkmale:

• - Wenigstens ausgewählte Werte eines ersten Bilddatensatzes werden in einem ersten Vektor notiert,
• - entsprechend den ausgewählten Werten des ersten Bilddaten­ satzes werden Werte eines zweiten Bilddatensatzes, der dem ersten zeitlich nachfolgend erzeugt worden ist, ausgewählt und in einem zweiten Vektor notiert,
• - in einem dritten Vektor werden sechs Parameter notiert, mit denen eine beliebige Lageveränderung im dreidimensionalen Raum beschreibbar ist,
• - gemäß einer Taylor-Entwicklung erster Ordnung wird eine Gleichung gebildet, bei der eine Differenz zwischen dem zweiten und dem ersten Vektor gleich einem Produkt einer Jacobischen Funktionalmatrix mit dem dritten Vektor gesetzt wird, wobei die Jacobische Funktionalmatrix je Zeile par­ tielle Ableitungen des entsprechenden Wertes des ersten Vektors nach den sechs Parametern umfasst, und
• - die Gleichung zur Bestimmung der sechs Parameter wird durch ein Iterationsverfahren gelöst.

9. Verfahren nach Anspruch 7, beinhaltend folgende Merkmale:

• - Wenigstens für ausgewählte Werte eines ersten Bilddatensat­ zes wird eine erste Häufigkeitsverteilung n(x) gebildet,
• - für entsprechend den ausgewählten Werten des ersten Bildda­ tensatzes ausgewählte Werte eines zweiten Bilddatensatzes, der dem ersten zeitlich nachfolgend erzeugt worden ist, wird eine zweite Häufigkeitsverteilung n(y) gebildet,
• - aus in den Bilddatensätzen entsprechend positionierten Wer­ ten werden Wertepaare gebildet, mit denen eine weitere Häu­ figkeitsverteilung n(x, y) gebildet wird,
• - aus der ersten, zweiten und weiteren Häufigkeitsverteilung wird eine Mutual Information
  gebildet,
• - einer der Bilddatensätze wird mit Parametern versehen, so dass entsprechend einer beliebigen Lageveränderung des ab­ zubildenden Bereichs im dreidimensionalen Raum der Bildda­ tensatz anpassbar ist, und
• - durch ein Optimierungsverfahren werden die Parameter so bestimmt, dass die Mutual Information minimal wird.

10. Verfahren nach einem der Ansprüche 1 bis 9, wobei die Bilddatensätze dreidimensionale Bilddatensätze sind.

11. Verfahren nach einem der Ansprüche 1 bis 10, wobei die Bilddatensätze mit einer Schnellbildtechnik erzeugt werden.

12. Verfahren nach Anspruch 11, wobei die Schnellbildtech­ nik ein Echoplanarverfahren umfasst.

13. Verfahren nach einem der Ansprüche 1 bis 12, wobei dem Untersuchungsobjekt zum Erzeugen von Bilddatensätzen ein Kon­ trastmittel verabreicht wird.