WO2022180067A1

WO2022180067A1 - Mr imaging with reduction of zipper artifacts

Info

Publication number: WO2022180067A1
Application number: PCT/EP2022/054468
Authority: WO
Inventors: Enrico Pannicke
Original assignee: Otto-Von-Guericke-Universität Magdeburg
Priority date: 2021-02-24
Filing date: 2022-02-23
Publication date: 2022-09-01
Also published as: DE102021104344A1

Abstract

The invention relates to a method for MR imaging of an object (10) which is located in a main magnetic field inside an examination volume of an MR device (1). The aim of the invention is to provide a method for MR imaging which reduces zipper artifacts (21, 51). The aim is achieved in that the method proposed by the invention comprises the following steps: - applying an imaging sequence to the object (10), the imaging sequence comprising a series of HF pulses and connected magnetic field gradients, in order to generate frequency- and phase-coded MR signals, - recording the MR signals from the object (10) in a specified k-space region (41), - detecting interference signal components (31) on the basis of a spectral analysis of the MR signals, - synthesising the detected interference signal components, - correcting the recorded MR signals by subtracting the synthesised interference signal components (32) from the recorded MR signals, and - reconstructing an MR image from the corrected MR signals. The invention further relates to an MR device (1) for carrying out the method, and to a computer program to be run on an MR device (1).

Description

MR-Bildqebunq mit Reduktion von Zipper-Artefakten Die Erfindung betrifft das Gebiet der Magnetresonanztomographie (MR). Sie betrifft ein Verfahren zur MR-Bildgebung eines Objekts, das sich in einem Hauptmagnetfeld innerhalb eines Untersuchungsvolumens eines MR-Geräts befindet. Die Erfindung bezieht sich auch auf ein MR-Gerät und auf ein Computerprogramm, das auf einem MR-Gerät ausgeführt wird. Bildgebende MR-Verfahren, die die Wechselwirkung zwischen Magnetfeldern und Kernspins nutzen, um zwei- oder dreidimensionale Bilder zu erzeugen, sind heute vor allem in der medizinischen Diagnostik weit verbreitet, weil sie zur Darstellung von Weichteilgewebe anderen bildgebenden Verfahren in vielerlei Hinsicht überlegen sind, keine ionisierende Strahlung benötigen und in der Regel nicht invasiv sind. MR image generation with reduction of zipper artefacts The invention relates to the field of magnetic resonance tomography (MR). It relates to a method for MR imaging of an object located in a main magnetic field within an examination volume of an MR device. The invention also relates to an MR device and to a computer program that runs on an MR device. Imaging MR methods, which use the interaction between magnetic fields and nuclear spins to generate two- or three-dimensional images, are widespread today, especially in medical diagnostics, because they are superior to other imaging methods for displaying soft tissue in many respects, none require ionizing radiation and are usually non-invasive.

Bei der MR-Methode im Allgemeinen wird das Objekt, d.h. der Körper des zu untersuchenden Patienten in einem starken, gleichmäßigen Magnetfeld Bo angeordnet, dessen Richtung gleichzeitig eine Achse (normalerweise die z- Achse) des Koordinatensystems definiert, das der Messung zugrunde liegt. Das Magnetfeld Bo erzeugt in Abhängigkeit von der Magnetfeldstärke unterschiedliche Energieniveaus für die einzelnen Kernspins, die durch Anlegen eines elektromagnetischen Wechselfeldes (HF-Feld) definierter Frequenz (sog. Larmor- Frequenz oder MR-Resonanzfrequenz) angeregt werden können (Spinresonanz). Aus makroskopischer Sicht ergibt sich aus der Verteilung der einzelnen Kernspins eine Gesamtmagnetisierung, die durch Anlegen eines elektromagnetischen Impulses geeigneter Frequenz (HF-Puls) senkrecht zur z-Achse aus dem Gleichgewichtszustand ausgelenkt werden kann, so dass die Magnetisierung eine Präzessionsbewegung um die z-Achse ausführt. Die Präzessionsbewegung der Magnetisierung beschreibt eine Oberfläche eines Kegels, dessen Öffnungswinkel als Flip-Winkel bezeichnet wird. Die Größe des Flip-Winkels ist abhängig von der Stärke und der Dauer des angelegten elektromagnetischen Impulses. Im Falle eines sogenannten 90°-Pulses werden die Spins von der z- Achse in die Querebene abgelenkt (Flip-Winkel = 90°). In the MR method in general, the object, ie the body of the patient to be examined, is placed in a strong, uniform magnetic field Bo, the direction of which also defines an axis (usually the z-axis) of the coordinate system on which the measurement is based. Depending on the magnetic field strength, the magnetic field Bo generates different energy levels for the individual nuclear spins, which can be excited (spin resonance) by applying an alternating electromagnetic field (HF field) with a defined frequency (the so-called Larmor frequency or MR resonance frequency). From a macroscopic point of view, the distribution of the individual nuclear spins results in a total magnetization that can be deflected from the equilibrium state perpendicular to the z-axis by applying an electromagnetic pulse of a suitable frequency (HF pulse), so that the magnetization performs a precession movement around the z-axis. The precession movement of the magnetization describes a surface of a cone, the opening angle of which is called the flip angle. The size of the flip angle depends on the strength and duration of the applied electromagnetic pulse. In the case of a so-called 90° pulse, the spins are deflected from the z-axis into the transverse plane (flip angle = 90°).

Nach Beendigung des FIF-Pulses relaxiert die Magnetisierung wieder in den ursprünglichen Gleichgewichtszustand, indem die Magnetisierung in z-Richtung mit einer ersten Zeitkonstante Ti (Spin-Gitter- oder longitudinale Relaxationszeit) wieder aufgebaut wird und die Magnetisierung in der Richtung senkrecht zur z- Richtung mit einer zweiten Zeitkonstante T2 (Spin-Spin- oder transversale Relaxationszeit) relaxiert. Die Variation der Magnetisierung kann mittels HF- Antennen als MR-Signal detektiert werden, wobei die FIF-Antennen innerhalb des Untersuchungsvolumens des MR-Gerätes so angeordnet und orientiert sind, dass die Variation der Magnetisierung in der Richtung senkrecht zur z-Achse gemessen wird. Der Zerfall der transversalen Magnetisierung geht nach Anlegen z.B. eines 90°-Pulses mit einem Übergang der Kernspins (induziert durch lokale Magnetfeld-Inhomogenitäten) von einem geordneten Zustand mit gleicher Phase in einen Zustand einher, in dem alle Phasenwinkel gleichmäßig verteilt sind (Dephasierung). Die Dephasierung kann durch einen HF-Refokussierungspuls (z.B. einen 180°-Puls) kompensiert werden. Dadurch wird in den HF-Antennen ein Echosignal erzeugt. After termination of the FIF pulse, the magnetization relaxes back to the original equilibrium state by rebuilding the magnetization in the z-direction with a first time constant Ti (spin lattice or longitudinal relaxation time) and the magnetization in the direction perpendicular to the z-direction with a second time constant T2 (spin-spin or transverse relaxation time). The variation in magnetization can be detected as an MR signal using HF antennas, with the FIF antennas being arranged and oriented within the examination volume of the MR device in such a way that the variation in magnetization is measured in the direction perpendicular to the z-axis. After the application of a 90° pulse, for example, the decay of the transverse magnetization is accompanied by a transition of the nuclear spins (induced by local magnetic field inhomogeneities) from an ordered state with the same phase to a state in which all phase angles are evenly distributed (dephasing). The dephasing can be compensated by an RF refocusing pulse (e.g. a 180° pulse). This generates an echo signal in the HF antennas.

Um eine Ortsauflösung im Körper zu realisieren, werden geschaltete Magnetfeldgradienten, die sich entlang der drei Hauptachsen erstrecken, dem gleichförmigen Magnetfeld Bo überlagert, was zu einer linearen Ortsabhängigkeit der MR-Resonanzfrequenz führt. Das in den HF-Antennen aufgenommene MR- Signal enthält dann Komponenten mit unterschiedlichen Frequenzen, die verschiedenen Orten im Körper zugeordnet werden können. Die über die HF- Antennen erfassten Signaldaten entsprechen dem Ortsfrequenzbereich (k-Raum genannt) und werden als k-Raumdaten bezeichnet. Die k-Raum-Daten umfassen in der Regel mehrere Datenzeilen, mit Frequenzkodierung entlang jeder Zeile, wobei die Datenzeilen mit jeweils unterschiedlicher Phasenkodierung erfasst werden. Jede digitale Datenzeile wird durch das Aufzeichnen einer Anzahl von Abtastwerten des empfangenen MR-Signals erhalten. Ein Satz solcher k-Raum- Daten wird in ein MR-Bild umgewandelt, z.B. mittels Fourier-Transformation. In order to achieve spatial resolution in the body, switched magnetic field gradients, which extend along the three main axes, are superimposed on the uniform magnetic field Bo, which leads to a linear spatial dependence of the MR resonance frequency. The MR signal recorded in the HF antennas then contains components with different frequencies that can be assigned to different locations in the body. The signal data acquired via the HF antennas correspond to the spatial frequency range (called k-space) and are referred to as k-space data. The k-space data typically comprises multiple data lines, with frequency encoding along each line, the data lines being acquired with different phase encoding in each case. Each line of digital data is created by recording a number of Obtain samples of the received MR signal. A set of such k-space data is converted into an MR image, for example by means of Fourier transformation.

Im rekonstruierten MR-Bild zeigen sich gelegentlich sogenannte Zipper-Artefakte in Form von reißverschlussartigen Bändern oder Linien von Störsignalen, die an einer bestimmten Position in der Frequenzkodierungsrichtung, parallel zur Phasenkodierungsrichtung durch das Bild verlaufen. Zipper-Artefakte sind zumeist auf HF-Signale von Störquellen im Bereich der MR-Resonanzfrequenz zurückzuführen, die entweder innerhalb der HF-Abschirmung des verwendeten MR-Gerätes (z.B. von in der Nähe des MR-Gerätes befindlichen elektronischen Geräten, wie etwa Anästhesie-Überwachungsgeräten) emittiert werden oder außerhalb davon und trotz HF-Abschirmung bis zur Empfangselektronik des MR- Gerätes durchdringen. The reconstructed MR image occasionally shows so-called zipper artifacts in the form of zipper-like bands or lines of interference signals that run through the image at a specific position in the frequency-encoding direction, parallel to the phase-encoding direction. Zipper artifacts are mostly due to HF signals from sources of interference in the range of the MR resonance frequency, which are either within the HF shielding of the MR device used (e.g. from electronic devices located in the vicinity of the MR device, such as anesthesia monitors ) are emitted or penetrate outside of it and despite HF shielding to the receiving electronics of the MR device.

Es ist Aufgabe der Erfindung, ein Verfahren für die MR-Bildgebung bereitzustellen, das Zipper-Artefakte reduziert. It is the object of the invention to provide a method for MR imaging that reduces zipper artifacts.

Diese Aufgabe löst die Erfindung durch ein Verfahren zur MR-Bildgebung eines Objekts, das in einem Hauptmagnetfeld innerhalb eines Untersuchungsvolumens eines MR-Gerätes befindlich ist, wobei das Verfahren die folgenden Schritte umfasst: The invention achieves this object by a method for MR imaging of an object that is located in a main magnetic field within an examination volume of an MR device, the method comprising the following steps:

Anwendung einer Bildgebungssequenz auf das Objekt, wobei die Bildgebungssequenz eine Abfolge von HF-Pulsen und geschalteten Magnetfeldgradienten umfasst, um frequenz- und phasenkodierte MR-Signale zu erzeugen, applying an imaging sequence to the subject, the imaging sequence comprising a train of RF pulses and switched magnetic field gradients to generate frequency and phase encoded MR signals,

Erfassung der MR-Signale von dem Objekt in einem vorgegebenen k-Raum Bereich, acquisition of the MR signals from the object in a predetermined k-space region,

Detektion von Störsignalanteilen auf Basis einer spektralen Analyse der MR-Signale, Detection of interference signal components based on a spectral analysis of the MR signals,

Synthese der detektierten Störsignalanteile, Synthesis of the detected interference signal components,

Bereinigung der erfassten MR-Signale durch Subtraktion der synthetisierten Störsignalanteile von den erfassten MR-Signalen, undcleaning up the detected MR signals by subtracting the synthesized noise signal components from the detected MR signals, and

Rekonstruktion eines MR-Bildes aus den bereinigten MR-Signalen Zipper-Artefakte der oben beschriebenen Art bilden ein Band oder eine Linie in Phasenkodierungsrichtung. Das Artefakt entsteht durch Überlagerung des eigentlichen Nutzsignals mit dem Störsignal bei einer bestimmten Frequenz fo. Reconstruction of an MR image from the cleaned MR signals Zipper artifacts of the type described above form a band or line in the phase encode direction. The artifact is caused by the interference signal being superimposed on the actual useful signal at a specific frequency fo.

S(n, p ) = S_MR(n, p ) + S_RF Dabei ist S_MR(n,p ) das Nutzsignal, d.h. das durch die Bildgebungssequenz erzeugte eigentliche MR-Signal, S_RF ist das Störsignal, S(n,p ) ist das erfasste, digitale Signal im k-Raum. n und p geben die Koordinaten im k-Raum in Frequenz- bzw. Phasenkodierungsrichtung in Form der Indizes der entsprechenden digitalen Datenpunkte an. Das Störsignal kann dabei z.B. wie folgt beschrieben werden:

S(n, p ) = S _MR (n, p ) + S _RF where S _MR (n,p ) is the useful signal, ie the actual MR signal generated by the imaging sequence, S _RF is the interference signal, S(n, p ) is the acquired digital signal in k-space. n and p indicate the coordinates in k-space in the direction of frequency and phase encoding, respectively, in the form of the indices of the corresponding digital data points. The interference signal can be described as follows, for example:

Dabei ist A die Amplitude des Störsignals, t_n gibt den zu dem Datenpunkt bei der k-Raum-Koordinate n in Frequenzkodierungsrichtung gehörigen Abtastzeitpunkt an. Die Phase cp des Störsignals ist bei der Erfassung der MR-Signale zufällig und damit im Allgemeinen von Phasenkodierungsschritt zu Phasenkodierungsschritt unterschiedlich. Auch die Amplitude kann von Phasenkodierungsschritt zu Phasenkodierungsschritt unterschiedlich ausfallen, wenn das Störsignal zeitlich variiert. Durch spektrale Analyse der erfassten MR-Signale werden erfindungsgemäß die Störsignalanteile S_RF detektiert, und zwar individuell für jeden Phasenkodierungsschritt, d.h. einzeln für jede Datenzeile der erfassten MR- Signale. Auf Basis des obigen Modells des Störsignals werden die Parameter Amplitude, Frequenz und Phase des Störsignals ermittelt. Daraus wird das Störsignal für jede Datenzeile synthetisiert, und die erfassten MR-Signale werden durch Subtraktion des synthetischen Störsignals entsprechend bereinigt. Von Imperfektionen bei der Detektion und des zugrunde gelegten Modells des Störsignals abgesehen ist das so bereinigte MR-Signal von dem Störsignal befreit und enthält nur noch das Nutzsignal. Aus letzterem wird schließlich das MR-Bild rekonstruiert, das entsprechend keine oder zumindest nur noch deutlich reduzierte Zipper-Artefakte aufweist. Bei einer bevorzugten Ausgestaltung des Verfahrens ist die spektrale Analyse auf einen oder mehrere periphere Teilbereiche des vorgegebenen k-Raum Bereichs beschränkt. Das Nutzsignal, d.h. das eigentliche MR-Signal ist hauptsächlich im Zentrum des k-Raums lokalisiert. Dort ist die wesentliche Bildinformation kodiert. In der Peripherie des k-Raums, d.h. außerhalb des k-Raum-Zentrums hat das Nutzsignal nur geringe Intensität, dort befindet sich weniger relevante Bildinformation. Das Störsignal ist demgegenüber in der k-Raum-Peripherie und im k-Raum-Zentrum gleichermaßen vorhanden. Dies nutzt die Erfindung aus, indem gezielt nur die k-Raum-Peripherie verwendet wird, um die Störsignalanteile zu detektieren. Dort kann bei der Detektion der Störsignalanteile das Nutzsignal vernachlässigt werden. Die auf dieser Basis dann synthetisierten Störsignalanteile werden allerdings verwendet, um die erfassten MR-Signale im gesamten vorgegebenen k-Raum-Bereich, d.h. sowohl im Zentrum als auch in der Peripherie zu bereinigen. In this case, A is the amplitude of the interference signal, t _n indicates the sampling time associated with the data point at the k-space coordinate n in the frequency coding direction. The phase cp of the interference signal is random when the MR signals are acquired and thus generally differs from phase encoding step to phase encoding step. The amplitude can also vary from phase encoding step to phase encoding step if the interference signal varies over time. According to the invention, the interference signal components S _RF are detected by spectral analysis of the acquired MR signals, specifically individually for each phase coding step, ie individually for each data line of the acquired MR signals. Based on the above model of the interference signal, the parameters amplitude, frequency and phase of the interference signal are determined. From this, the interference signal is synthesized for each data line, and the acquired MR signals are cleaned accordingly by subtracting the synthetic interference signal. Apart from imperfections in the detection and the underlying model of the interference signal, the MR signal cleaned up in this way is freed from the interference signal and only contains the useful signal. Finally, the MR image is reconstructed from the latter, which correspondingly has no or at least only significantly reduced zipper artefacts. In a preferred embodiment of the method, the spectral analysis is limited to one or more peripheral sub-areas of the specified k-space area. The useful signal, ie the actual MR signal, is mainly localized in the center of k-space. The essential image information is encoded there. In the periphery of k-space, ie outside the center of k-space, the useful signal has only a low intensity, there is less relevant image information there. In contrast, the interference signal is equally present in the k-space periphery and in the k-space center. The invention exploits this in that only the k-space periphery is used in a targeted manner in order to detect the interference signal components. There, the useful signal can be neglected when detecting the interference signal components. However, the interference signal components then synthesized on this basis are used in order to clean up the acquired MR signals in the entire specified k-space region, ie both in the center and in the periphery.

Bei einer bevorzugten Ausgestaltung erfolgt die spektrale Analyse durch Fourier- Transformation der erfassten MR-Signale in der Frequenzkodierungsrichtung. Es eignet sich die bekannte Technik der diskreten Fourier-Transformation, wie sie auch für die Bildrekonstruktion verwendet wird. In a preferred embodiment, the spectral analysis is performed by Fourier transformation of the detected MR signals in the frequency encoding direction. The well-known technique of discrete Fourier transformation is suitable, as is also used for image reconstruction.

Da, wie oben erwähnt, das Störsignal in keiner festen Phasenbeziehung zur MR- Detektion steht, wird zweckmäßig jedes der unterschiedlich phasenkodierten MR- Signale, d.h. jede Datenzeile einzeln den Schritten der Detektion und der Synthese von Störsignalanteilen und der entsprechenden Bereinigung des erfassten MR-Signals unterzogen. Häufig wird zur Verbesserung des Signal-zu- Rausch-Verhältnisses eine Signalakkumulation durchgeführt, d.h. jedes der frequenz- und phasenkodierten MR-Signale wird zum Zwecke der Erhöhung des Signal-zu-Rausch-Verhältnisses zwei- oder mehrfach erfasst, und die mehrfach erfassten MR-Signale werden addiert. Auch in diesem Falle hat das Störsignal bei der Erfassung jedes einzelnen MR-Signals eine zufällige Phase, so dass jedes einzelne der mehrfach erfassten MR-Signale vor der Addition den Schritten der Detektion und der Synthese von Störsignalanteilen und der entsprechenden Bereinigung des erfassten MR-Signals unterzogen werden sollte. Bei einer weiteren bevorzugten Ausgestaltung erfolgt das Erfassen der MR- Signale durch Empfang mittels einer HF-Antenne und durch Digitalisierung der von der Antenne empfangenen Signale, wobei das von der Antenne empfangene Signal mit zumindest der doppelten MR-Resonanzfrequenz abgetastet wird. Es erfolgt mit anderen Worten nicht, wie im Stand der Technik üblich, zunächst eine Mischung der über die HF-Antenne empfangenen MR-Signale mit der MR- Resonanzfrequenz mit anschließender Tiefpassfilterung und Digitalisierung im Audiofrequenzbereich. Stattdessen werden die MR-Signale ohne Mischung direkt im Frequenzband der MR-Resonanzfrequenz detektiert. Dies erfordert nach dem Nyquist-Shannon-Theorem, dass die Abtastfrequenz zumindest gleich der doppelten MR-Resonanzfrequenz ist. Die digitalen MR-Signaldaten liegen dann in Frequenzkodierungsrichtung im HF-Bereich vor, was vorteilhaft für die erfindungsgemäße Detektion der Störsignalanteile ist. Bei dieser Vorgehensweise lassen sich Nutzsignal und Störsignal bei der spektralen Analyse der erfassten MR-Signale optimal voneinander unterscheiden. Since, as mentioned above, the interference signal is not in a fixed phase relationship to the MR detection, each of the differently phase-coded MR signals, ie each data line, is individually subjected to the steps of detection and synthesis of interference signal components and the corresponding cleanup of the detected MR signal subjected. Signal accumulation is often carried out to improve the signal-to-noise ratio, ie each of the frequency- and phase-coded MR signals is recorded twice or more in order to increase the signal-to-noise ratio, and the multiple recorded MR -Signals are added. In this case, too, the interference signal has a random phase during the acquisition of each individual MR signal, so that each of the multiple acquired MR signals before the addition of the steps of detection and synthesis of interference signal components and the corresponding cleaning of the acquired MR signal should be subjected to. In a further preferred embodiment, the MR signals are detected by receiving them using an HF antenna and by digitizing the signals received by the antenna, the signal received by the antenna being sampled at at least twice the MR resonance frequency. In other words, the MR signals received via the HF antenna are not initially mixed with the MR resonant frequency with subsequent low-pass filtering and digitization in the audio frequency range, as is customary in the prior art. Instead, the MR signals are detected directly in the frequency band of the MR resonance frequency without being mixed. According to the Nyquist-Shannon theorem, this requires that the sampling frequency is at least twice the MR resonance frequency. The digital MR signal data are then present in the HF range in the frequency coding direction, which is advantageous for the inventive detection of the interference signal components. With this procedure, useful signal and interference signal can be optimally distinguished from one another in the spectral analysis of the detected MR signals.

Wie oben angesprochen, kann die Synthese der Störsignalanteile im einfachsten Fall auf einem Modell basieren, das jedem Störsignal eine durch Amplitude, Frequenz und Phase definierte harmonische Oszillation zuordnet. Andere, komplexere Modelle sind denkbar, falls sich im konkreten Anwendungsfall zeigt, dass die Unterdrückung der Zipper-Artefakte auf Basis des einfachen Modells unzureichend ist. As mentioned above, in the simplest case, the synthesis of the interference signal components can be based on a model that assigns a harmonic oscillation defined by amplitude, frequency and phase to each interference signal. Other, more complex models are conceivable if the specific application shows that the suppression of zipper artifacts based on the simple model is insufficient.

Das bisher beschriebene erfindungsgemäße Verfahren kann mittels eines MR- Gerätes durchgeführt werden, das mindestens eine Hauptmagnetspule zur Erzeugung eines im Wesentlichen gleichförmigen, statischen Hauptmagnetfeldes Bo innerhalb eines Untersuchungsvolumens, eine Anzahl von Gradientenspulen zur Erzeugung geschalteter Magnetfeldgradienten in unterschiedlichen Raumrichtungen innerhalb des Untersuchungsvolumens aufweist, mindestens eine HF-Antenne zum Erzeugen von HF-Pulsen innerhalb des Untersuchungsvolumens und/oder zum Empfangen von MR-Signalen von einem im Untersuchungsvolumen positionierten Objektes (Körper eines Patienten), eine Steuereinheit zum Steuern der zeitlichen Abfolge von HF-Pulsen und geschalteten Magnetfeldgradienten und eine Rekonstruktionseinheit zum Rekonstruieren von MR-Bildern aus den empfangenen MR-Signalen. Das Verfahren der Erfindung kann durch eine entsprechende Programmierung der Rekonstruktionseinheit und/oder der Steuereinheit des MR-Gerätes realisiert werden. The method according to the invention described so far can be carried out using an MR device which has at least one main magnet coil for generating a substantially uniform, static main magnetic field Bo within an examination volume, a number of gradient coils for generating switched magnetic field gradients in different spatial directions within the examination volume, at least one HF antenna for generating HF pulses within the examination volume and/or for receiving MR signals from an object positioned in the examination volume (body of a patient), a control unit for controlling the time sequence of HF pulses and switched magnetic field gradients and a reconstruction unit for reconstructing MR images from the received MR signals. That The method of the invention can be implemented by appropriate programming of the reconstruction unit and/or the control unit of the MR device.

Das Verfahren der Erfindung kann bei den meisten MR-Geräten, die derzeit im klinischen Einsatz sind, mit Vorteil durchgeführt werden. Dazu ist es lediglich notwendig, ein Computerprogramm zu verwenden, mit dem das MR-Gerät so gesteuert wird, dass es die oben erläuterten Verfahrensschritte der Erfindung ausführt. Das Computerprogramm kann sich entweder auf einem Datenträger oder in einem Datennetz befinden, so dass es zur Installation in die Steuereinheit des MR-Gerätes heruntergeladen werden kann. The method of the invention can be advantageously performed on most MR devices currently in clinical use. All that is necessary for this is to use a computer program with which the MR device is controlled in such a way that it executes the method steps of the invention explained above. The computer program can be located either on a data medium or in a data network so that it can be downloaded to the control unit of the MR device for installation.

Ausführungsbeispiele der Erfindung werden nachfolgend anhand der Zeichnungen näher erläutert. Es zeigen: Exemplary embodiments of the invention are explained in more detail below with reference to the drawings. Show it:

Figur 1 ein MR-Gerät zur Durchführung des Verfahrens der Erfindung, FIG. 1 an MR device for carrying out the method of the invention,

Figur 2 ein MR-Bild mit Zipper-Artefakt, Figur 3 eine Illustration der spektralen Analyse des empfangenen MR-Signals gemäß der Erfindung, FIG. 2 shows an MR image with a zipper artifact, FIG. 3 shows an illustration of the spectral analysis of the received MR signal according to the invention,

Figur 4 ein Diagramm des k-Raums mit Aufteilung in verschiedene Teilbereiche und FIG. 4 shows a diagram of k-space divided into different partial areas and

Figur 5 zwei MR-Bilder ohne (links) und mit (rechts) Figure 5 two MR images without (left) and with (right)

Artefakt-Reduktion gemäß der Erfindung. Artifact reduction according to the invention.

Unter Bezugnahme auf Figur 1 ist ein MR-Gerät 1 als Blockschaltbild dargestellt. Das Gerät besteht aus supraleitenden oder resistiven Flauptmagnetspulen 2, so dass entlang einer z-Achse durch ein Untersuchungsvolumen ein im Wesentlichen gleichmäßiges, zeitlich konstantes Flauptmagnetfeld Bo erzeugt wird. Das Gerät umfasst ferner einen Satz von Shimspulen 2' (1., 2. und - gegebenenfalls - 3. Ordnung), wobei der Stromfluss durch die einzelnen Shimspulen des Satzes 2' steuerbar ist, um Bo-Abweichungen innerhalb des Untersuchungsvolumens zu minimieren. With reference to FIG. 1, an MR device 1 is shown as a block diagram. The device consists of superconducting or resistive primary magnet coils 2, so that a substantially uniform, time-constant primary magnetic field Bo is generated along a z-axis through an examination volume. The device also includes a set of shim coils 2 '(1st, 2nd and - if necessary - 3rd order), wherein the current flow through the individual Shim coils of set 2' is controllable to minimize Bo-deviations within the examination volume.

Ein System zur Erzeugung und Manipulation von Magnetresonanz wendet eine Reihe von HF-Pulsen und geschalteten Magnetfeldgradienten an, um magnetische Kernspins zu invertieren oder anzuregen, Magnetresonanz zu induzieren, Magnetresonanz zu refokussieren, Magnetresonanz zu manipulieren, die Magnetresonanz räumlich und anderweitig zu kodieren, Spins zu sättigen und ähnliches, um MR-Bildgebung durchzuführen. A magnetic resonance generation and manipulation system applies a series of RF pulses and switched magnetic field gradients to invert or excite nuclear magnetic spins, induce magnetic resonance, refocus magnetic resonance, manipulate magnetic resonance, spatially and otherwise encode magnetic resonance, spins saturate and the like to perform MR imaging.

Genauer gesagt, ein Gradientenimpulsverstärker 3 legt Stromimpulse an ausgewählte Ganzkörper-Gradientenspulen 4, 5 und 6 entlang der x-, y- und z- Achse des Untersuchungsvolumens an. Ein HF-Sender 7 sendet FIF-Pulse oder Pulspakete über einen Sende-/Empfangsumschalter 8 und eine Ganzkörper-FIF- Antenne 9 in das Untersuchungsvolumen. Eine typische MR- Bildgebungssequenz besteht aus einem Paket von FIF-Pulsen kurzer Dauer, die zusammen mit geschalteten Magnetfeldgradienten eine gezielte Manipulation der kernmagnetischen Resonanz bewirken. Die FIF-Pulse werden verwendet, um die Resonanz zu sättigen, die Resonanz anzuregen, die Magnetisierung umzukehren, die Resonanz zu refokussieren oder die Resonanz zu manipulieren und einen Teil (z.B. eine Schicht) eines im Untersuchungsvolumen positionierten Körpers 10 eines Patienten auszuwählen. Die MR-Signale werden auch von der FIF-Antenne 9 aufgenommen. More specifically, a gradient pulse amplifier 3 applies current pulses to selected whole body gradient coils 4, 5 and 6 along the x, y and z axes of the examination volume. An HF transmitter 7 transmits FIF pulses or pulse packets via a transmit/receive switch 8 and a whole-body FIF antenna 9 into the examination volume. A typical MR imaging sequence consists of a packet of short-duration FIF pulses which, together with switched magnetic field gradients, result in a targeted manipulation of nuclear magnetic resonance. The FIF pulses are used to saturate the resonance, excite the resonance, reverse the magnetization, refocus or manipulate the resonance, and select a portion (e.g., a slice) of a patient's body 10 positioned in the examination volume. The MR signals are also picked up by the FIF antenna 9 .

Für die Erzeugung von MR-Bildern von begrenzten Regionen des Körpers 10 wird ein Satz lokaler Array-FIF-Antennen 11 , 12, 13 angrenzend an die für die Bildgebung ausgewählte Region platziert. Die Array-Antennen 11 , 12, 13 können zum Empfang von MR-Signalen verwendet werden, die durch FIF-Übertragungen über die Körper-FIF-Antenne 9 induziert werden. For the generation of MR images of limited regions of the body 10, a set of local array FIF antennas 11, 12, 13 is placed adjacent to the region selected for imaging. The array antennas 11 , 12 , 13 can be used to receive MR signals induced by FIF transmissions via the body FIF antenna 9 .

Die resultierenden MR-Signale werden von der Körper-FIF-Antenne 9 und/oder von den Array-FIF-Antennen 11 , 12, 13 aufgenommen und von einem Empfänger 14, vorzugsweise mit Vorverstärker (nicht abgebildet), erfasst. Der Empfänger 14 ist über den Sende-/Empfangsumschalter 8 an die FIF-Antennen 9, 11 , 12 und 13 angeschlossen. Ein Hostrechner 15 steuert die Shimspulen 2' sowie den Gradientenimpulsverstärker 3 und den Sender 7 zur Erzeugung der Bildgebungssequenz gemäß der Erfindung. Für die gewählte Sequenz empfängt der Empfänger 14 nach jedem HF-Anregungspuls eine einzelne oder mehrere MR-Datenzeilen in schneller Folge. Ein Datenerfassungssystem 16 führt eine Analog-Digital-Wandlung der empfangenen Signale durch und wandelt jede MR- Datenzeile in ein für die Weiterverarbeitung geeignetes digitales Format um. In modernen MR-Geräten ist das Datenerfassungssystem 16 ein separater Computer, der auf die Erfassung von Rohbilddaten spezialisiert ist. Letztendlich werden die digitalen Rohbilddaten durch einen Rekonstruktionsprozessor 17, der eine Fourier-Transformation oder andere geeignete Rekonstruktionsalgorithmen anwendet, in eine Bilddarstellung rekonstruiert. Das resultierende MR-Bild kann einen planaren Schnitt durch den Patienten, eine Anordnung paralleler planarer Schnitte, ein dreidimensionales Volumen oder ähnliches darstellen. Das Bild wird dann in einem Bildspeicher gespeichert, wo auf das Bild zugegriffen werden kann, um Schnitte, Projektionen oder andere Teile der Bilddarstellung in ein geeignetes Format für die Visualisierung umzuwandeln, z.B. über einen Videomonitor 18, der eine Darstellung des MR-Bildes liefert. Der Hostrechner 15 und der Rekonstruktionsprozessor 17 sind so programmiert, dass sie das Verfahren der Erfindung wie oben und im Folgenden beschrieben ausführen. The resulting MR signals are picked up by the body FIF antenna 9 and/or by the array FIF antennas 11, 12, 13 and are detected by a receiver 14, preferably with a preamplifier (not shown). The receiver 14 is connected to the FIF antennas 9, 11, 12 and 13 via the transmit/receive switch 8. A host computer 15 controls the shim coils 2' as well as the gradient pulse amplifier 3 and the transmitter 7 to generate the imaging sequence according to the invention. For the selected sequence, the receiver 14 receives a single or multiple MR data lines in rapid succession after each RF excitation pulse. A data acquisition system 16 performs analog to digital conversion of the received signals and converts each line of MR data into a digital format suitable for further processing. In modern MR devices, the data acquisition system 16 is a separate computer that specializes in acquiring raw image data. Ultimately, the raw digital image data is reconstructed into an image representation by a reconstruction processor 17 applying a Fourier transform or other suitable reconstruction algorithm. The resulting MR image can represent a planar slice through the patient, an array of parallel planar slices, a three-dimensional volume, or the like. The image is then stored in an image store where the image can be accessed to convert slices, projections or other parts of the image representation into a suitable format for visualization, eg via a video monitor 18 which provides a display of the MR image. The host computer 15 and the reconstruction processor 17 are programmed to carry out the method of the invention as described above and below.

Figur 2 zeigt ein konventionell erfasstes und rekonstruiertes MR-Bild eines Phantoms. Die Frequenzkodierungsrichtung FE und die Phasenkodierungs- richtung PE sind durch Pfeile angedeutet. In dem MR-Bild ist ein Zipper-Artefakt 21 in Form einer reißverschlussartigen Linie zu erkennen. Die Störung befindet sich an einer bestimmten Position in der Frequenzkodierungsrichtung und läuft parallel zur Phasenkodierungsrichtung, d.h. bei dem dargestellten Beispiel in vertikaler Richtung durch das MR-Bild. Ursache des Zipper-Artefakts 21 ist ein HF-Signal einer Störquelle im Bereich der MR-Resonanzfrequenz, d.h. im Bereich der Empfangsbandbreite des MR-Gerätes 1. Das Artefakt 21 entsteht durch Überlagerung des eigentlichen Nutzsignals mit dem Störsignal bei einer bestimmten Frequenz fo, die zu der Bildposition des Artefakts 21 in der Frequenzkodierungsrichtung korrespondiert. Die Phase des Störsignals ist bei der Erfassung der MR-Signale zufällig und damit im Allgemeinen von Phasenkodierungsschritt zu Phasenkodierungsschritt unterschiedlich. Auch die Amplitude kann von Phasenkodierungsschritt zu Phasenkodierungsschritt unterschiedlich ausfallen, wenn das Störsignal zeitlich variiert. Daraus resultiert das reißverschlussartige Aussehen des Artefakts 21. FIG. 2 shows a conventionally recorded and reconstructed MR image of a phantom. The frequency encoding direction FE and the phase encoding direction PE are indicated by arrows. A zipper artifact 21 in the form of a zipper-like line can be seen in the MR image. The disturbance is at a specific position in the frequency encoding direction and runs parallel to the phase encoding direction, ie in the illustrated example in the vertical direction through the MR image. The cause of the zipper artifact 21 is an HF signal from an interference source in the range of the MR resonance frequency, ie in the range of the reception bandwidth of the MR device 1. The artefact 21 is created by superimposing the actual useful signal with the interference signal at a specific frequency fo, which corresponds to the image position of the artefact 21 in the frequency coding direction. The phase of the interference signal is random when the MR signals are acquired and thus generally differs from phase encoding step to phase encoding step. The amplitude can also vary from phase encoding step to phase encoding step if the interference signal varies over time. This results in the zipper-like appearance of Artifact 21.

Gemäß der Erfindung erfolgt durch spektrale Analyse, z. B. durch diskrete Fourier- Transformation, der einzelnen erfassten Datenzeilen die Detektion der Störsignalanteile, die das Zipper-Artefakt 21 verursachen, d.h. individuell für jeden Phasenkodierungsschritt. Dies illustriert die Figur 3. Die Figur 3 zeigt das diskrete Spektrum (Leistungsspektrum) einer Datenzeile des empfangenen MR-Signals. Dabei ist spektrale Analyse auf einen peripheren Teilbereich B des vorgegebenen k-Raum Bereichs beschränkt, wie in Figur 4 schematisch dargestellt. Die Figur 4 zeigt ein Diagramm des k-Raums mit Phasenkodierungsrichtung kpE und Frequenzkodierungsrichtung kFE. Der von der Bildgebungssequenz vorgegebene k-Raum Bereich ist durch das umgrenzende Rechteck 41 angedeutet. Das Nutzsignal, d.h. das eigentliche MR-Signal ist hauptsächlich im Zentrum des k- Raums lokalisiert, d.h. im Teilbereich A. Dort ist die wesentlichen Bildinformation kodiert. In der Peripherie des k-Raums, d.h. außerhalb des k-Raum-Zentrums (bei kpE=kFE=0, angedeutet durch die sich kreuzenden gestrichelten Linien) hat das Nutzsignal nur geringe Intensität, dort befindet sich weniger relevante Bildinformation. Das Störsignal ist demgegenüber in der k-Raum-Peripherie B und im k-Raum-Zentrum A gleichermaßen vorhanden. Dies nutzt die Erfindung aus, indem gezielt nur die k-Raum-Peripherie, d.h. die Teilbereiche B verwendet werden, um die Störsignalanteile zu detektieren. Die Figur 3 zeigt entsprechend das aus den Teilbereichen B berechnete Spektrum einer Datenzeile des empfangenen MR-Signals. In diesen Teilbereichen kann bei der Detektion der Störsignalanteile das Nutzsignal vernachlässigt werden. Das Störsignal zeigt sich in dem dargestellten Spektrum als einzelne Spektrallinie 31, die vom Signaluntergrund separiert werden kann. Das diskrete Spektrum führt zu Ungenauigkeiten bei der Bestimmung der Frequenz des Störsignals. Aus diesem Grund wird bei dem dargestellten Ausführungsbeispiel eine Gauß-Funktion 32 als Modellfunktion an die Spektrallinie angepasst (z.B. durch die Methode der kleinsten Abweichungsquadrate). Diese Prozedur ergibt gleichzeitig die Amplitude des Störsignals. Die Phase des Störsignals kann mit Kenntnis dieser Parameter dann in einem separaten Schritt z.B. anhand des komplexen Spektrums der analysierten Datenzeile (nicht dargestellt) ermittelt werden. Das Störsignal wird dann auf dieser Basis synthetisiert, und die erfassten MR-Signale werden durch Subtraktion des synthetischen Störsignals entsprechend bereinigt. Die synthetischen Störsignalanteile werden verwendet, um die erfassten MR- Signale im gesamten k-Raum-Bereich, d.h. sowohl in den Teilbereichen B als auch im Teilbereich A zu bereinigen, d.h. von den Störsignalanteilen zu befreien. Das so bereinigte MR-Signal enthält (theoretisch) nur noch das Nutzsignal. Die Zipper-Artefakte sind in dem aus dem bereinigten MR-Signal rekonstruierten MR- Bild entsprechend nicht mehr vorhanden oder zumindest schwächer ausgeprägt. According to the invention, by spectral analysis, e.g. B. by discrete Fourier transformation of the individual recorded data lines, the detection of the interference signal components that cause the zipper artefact 21, ie individually for each phase encoding step. This is illustrated in FIG. 3. FIG. 3 shows the discrete spectrum (power spectrum) of a data line of the received MR signal. In this case, the spectral analysis is restricted to a peripheral partial area B of the specified k-space area, as shown schematically in FIG. FIG. 4 shows a diagram of k-space with phase encoding direction kpE and frequency encoding direction kFE. The k-space region specified by the imaging sequence is indicated by the bounding rectangle 41 . The useful signal, ie the actual MR signal, is mainly localized in the center of k-space, ie in partial area A. The essential image information is encoded there. In the periphery of k-space, ie outside the k-space center (at kpE=kFE=0, indicated by the crossing dashed lines), the useful signal has only a low intensity, there is less relevant image information there. In contrast, the interference signal is equally present in the k-space periphery B and in the k-space center A. The invention makes use of this in that only the k-space periphery, ie the partial areas B, are used in a targeted manner in order to detect the interference signal components. FIG. 3 correspondingly shows the spectrum of a data line of the received MR signal calculated from the partial regions B. In these sub-areas, the useful signal can be neglected when detecting the interference signal components. The interference signal appears in the spectrum shown as a single spectral line 31, which can be separated from the signal background. The discrete spectrum leads to inaccuracies in determining the frequency of the interference signal. For this reason, in the illustrated embodiment, a Gaussian function 32 as Model function fitted to the spectral line (e.g. using the least squares method). This procedure also gives the amplitude of the interference signal. Knowing these parameters, the phase of the interference signal can then be determined in a separate step, for example using the complex spectrum of the analyzed data line (not shown). The noise signal is then synthesized on this basis and the acquired MR signals are cleaned up accordingly by subtracting the synthetic noise signal. The synthetic interference signal components are used in order to clean up the acquired MR signals in the entire k-space region, ie both in the partial areas B and in the partial area A, ie to free them from the interference signal components. The MR signal cleaned up in this way contains (theoretically) only the useful signal. The zipper artefacts are no longer present in the MR image reconstructed from the cleared MR signal, or are at least less pronounced.

Dies illustriert die Figur 5. Die Figur 5 zeigt links ein MR-Bild mit Zipper-Artefakten 51 , das konventionell aufgenommen und rekonstruiert wurde. Im rechten Bild, bei dem das Verfahren der Erfindung angewendet wurde, sind die Zipper-Artefakte 51 deutlich schwächer ausgeprägt. Dass die Zipper-Artefakte 51 im rechten Bild immer noch schwach erkennbar sind, kann auf Imperfektionen bei der Detektion und auf das zugrunde gelegten Modell der Störsignalanteile zurückzuführen sein, das die realen Störsignale möglicherweise nicht vollständig korrekt wiedergibt. This is illustrated in FIG. 5. On the left, FIG. 5 shows an MR image with zipper artefacts 51 that was conventionally recorded and reconstructed. In the image on the right, in which the method of the invention was applied, the zipper artefacts 51 are clearly less pronounced. The fact that the zipper artefacts 51 are still faintly discernible in the image on the right can be attributed to imperfections in the detection and to the underlying model of the interference signal components, which may not reproduce the real interference signals completely correctly.

Claims

Patentansprüche patent claims

1. Verfahren zur MR-Bildgebung eines Objekts (10), das in einem Hauptmagnetfeld innerhalb eines Untersuchungsvolumens eines MR-Gerätes (1 ) befindlich ist, wobei das Verfahren die folgenden Schritte umfasst: 1. A method for MR imaging of an object (10) that is located in a main magnetic field within an examination volume of an MR device (1), the method comprising the following steps:

Anwendung einer Bildgebungssequenz auf das Objekt (10), wobei die Bildgebungssequenz eine Abfolge von HF-Pulsen und geschalteten Magnetfeldgradienten umfasst, um frequenz- und phasenkodierte MR-Signale zu erzeugen, Application of an imaging sequence to the object (10), the imaging sequence comprising a succession of RF pulses and switched magnetic field gradients in order to generate frequency- and phase-coded MR signals,

Erfassung der MR-Signale von dem Objekt (10) in einem vorgegebenen k-Raum Bereich (41), Acquisition of the MR signals from the object (10) in a predetermined k-space region (41),

Detektion von Störsignalanteilen (31) auf Basis einer spektralen Analyse der MR-Signale, - Synthese der detektierten Störsignalanteile, Detection of interference signal components (31) based on a spectral analysis of the MR signals, - synthesis of the detected interference signal components,

Bereinigung der erfassten MR-Signale durch Subtraktion der synthetisierten Störsignalanteile (32) von den erfassten MR-Signalen, undCleaning up the detected MR signals by subtracting the synthesized noise signal components (32) from the detected MR signals, and

Rekonstruktion eines MR-Bildes aus den bereinigten MR-Signalen. Reconstruction of an MR image from the cleaned MR signals.

2. Verfahren nach Anspruch 1 , wobei die spektrale Analyse auf einen oder mehrere periphere Teilbereiche (B) des vorgegebenen k-Raum Bereichs2. The method according to claim 1, wherein the spectral analysis on one or more peripheral sub-areas (B) of the predetermined k-space area

(41) beschränkt wird. (41) is restricted.

3. Verfahren nach Anspruch 1 oder 2, wobei die spektrale Analyse durch Fourier-Transformation der erfassten MR-Signale in der Frequenz kodierungsrichtung (FE) erfolgt. 3. The method according to claim 1 or 2, wherein the spectral analysis is carried out by Fourier transformation of the detected MR signals in the frequency encoding direction (FE).

4. Verfahren nach einem der Ansprüche 1 bis 3, wobei jedes der unterschiedlich phasenkodierten MR-Signale einzeln den Schritten der Detektion und der Synthese von Störsignalanteilen und der entsprechenden Bereinigung des erfassten MR-Signals unterzogen wird. 5. Verfahren nach einem der Ansprüche 1 bis 4, wobei jedes der frequenz- und phasenkodierten MR-Signale zum Zwecke der Erhöhung des Signal-zu-Rausch-Verhältnisses zwei- oder mehrfach erfasst wird und die mehrfach erfassten MR-Signale addiert werden. 4. The method as claimed in one of claims 1 to 3, in which each of the differently phase-coded MR signals is individually subjected to the steps of detecting and synthesizing interference signal components and correspondingly cleaning up the detected MR signal. 5. The method according to any one of claims 1 to 4, wherein each of the frequency- and phase-coded MR signals is detected twice or more for the purpose of increasing the signal-to-noise ratio and the multiply detected MR signals are added.

6. Verfahren nach Anspruch 5, wobei jedes einzelne der mehrfach erfassten MR-Signale vor der Addition den Schritten der Detektion und der6. The method of claim 5, wherein each of the multiply detected MR signals before the addition of the steps of detection and the

Synthese von Störsignalanteilen und der entsprechenden Bereinigung des erfassten MR-Signals unterzogen wird. Synthesis of noise signal components and the corresponding cleanup of the detected MR signal is subjected.

7. Verfahren nach einem der Ansprüche 1 bis 6, wobei das Erfassen der MR-Signale durch Empfang mittels einer HF-Antenne (9, 11, 12, 13) und durch Digitalisierung der von der HF-Antenne (9, 11, 12, 13) empfangenen MR- Signale erfolgt, wobei das von der Antenne empfangene MR-Signal mit zumindest der doppelten MR-Resonanzfrequenz abgetastet wird. 7. The method according to any one of claims 1 to 6, wherein the acquisition of the MR signals by reception by means of an HF antenna (9, 11, 12, 13) and by digitizing the signals from the HF antenna (9, 11, 12, 13) received MR signals, the MR signal received by the antenna being sampled at at least twice the MR resonance frequency.

8. Verfahren nach einem der Ansprüche 1 bis 7, wobei die Synthese der Störsignalanteile auf einem Modell basiert, das jedem Störsignal eine durch Amplitude, Frequenz und Phase definierte harmonische Oszillation zuordnet. 8. The method as claimed in one of claims 1 to 7, in which the synthesis of the interference signal components is based on a model which assigns a harmonic oscillation defined by amplitude, frequency and phase to each interference signal.

9. MR-Gerät mit mindestens einer Hauptmagnetspule (2) zur Erzeugung eines statischen Hauptmagnetfeldes Bo innerhalb eines Untersuchungsvolumens, einer Anzahl von Gradientenspulen (4, 5, 6) zur Erzeugung geschalteter Magnetfeldgradienten in unterschiedlichen räumlichen Richtungen innerhalb des Untersuchungsvolumens, mindestens einer HF- Antenne (9) zur Erzeugung von HF-Pulsen innerhalb des9. MR device with at least one main magnetic coil (2) for generating a static main magnetic field Bo within an examination volume, a number of gradient coils (4, 5, 6) for generating switched magnetic field gradients in different spatial directions within the examination volume, at least one HF antenna (9) to generate RF pulses within the

Untersuchungsvolumens und/oder zum Empfang von MR-Signalen von einem im Untersuchungsvolumen befindlichen Objekt (10), einer Steuereinheit (15) zum Steuern der zeitlichen Abfolge der HF-Pulse und der geschalteten Magnetfeldgradienten, und einer Rekonstruktionseinheit (17) zum Rekonstruieren von MR-Bildern aus den MR-Signalen, wobei das MR-Gerät (1) zur Ausführung der Schritte des Verfahrens nach einem der Ansprüche 1 bis 8 eingerichtet ist. examination volume and/or for receiving MR signals from an object (10) located in the examination volume, a control unit (15) for controlling the time sequence of the RF pulses and the switched magnetic field gradients, and a reconstruction unit (17) for reconstruction of MR images from the MR signals, wherein the MR device (1) is set up to carry out the steps of the method according to one of Claims 1 to 8.

10. Computerprogramm, das auf einem MR-Gerät (1) ausführbar ist, wobei das Computerprogramm Befehle zur Ausführung des Verfahrens nach einem der Ansprüche 1 bis 8 umfasst. 10. Computer program that can be executed on an MR device (1), the computer program comprising instructions for executing the method according to any one of claims 1 to 8.