DE102018209584A1 - Magnetresonanz-Fingerprinting-Verfahren - Google Patents

Magnetresonanz-Fingerprinting-Verfahren Download PDF

Info

Publication number
DE102018209584A1
DE102018209584A1 DE102018209584.1A DE102018209584A DE102018209584A1 DE 102018209584 A1 DE102018209584 A1 DE 102018209584A1 DE 102018209584 A DE102018209584 A DE 102018209584A DE 102018209584 A1 DE102018209584 A1 DE 102018209584A1
Authority
DE
Germany
Prior art keywords
parameter
pixel
signal
determined
time series
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Pending
Application number
DE102018209584.1A
Other languages
English (en)
Inventor
Gregor Körzdörfer
Mathias Nittka
Josef Pfeuffer
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Siemens Healthineers Ag De
Original Assignee
Siemens Healthcare GmbH
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Siemens Healthcare GmbH filed Critical Siemens Healthcare GmbH
Priority to DE102018209584.1A priority Critical patent/DE102018209584A1/de
Priority to US16/441,383 priority patent/US10921406B2/en
Priority to CN201910515214.5A priority patent/CN110604569B/zh
Publication of DE102018209584A1 publication Critical patent/DE102018209584A1/de
Pending legal-status Critical Current

Links

Images

Classifications

    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
    • G01R33/00—Arrangements or instruments for measuring magnetic variables
    • G01R33/20—Arrangements or instruments for measuring magnetic variables involving magnetic resonance
    • G01R33/44—Arrangements or instruments for measuring magnetic variables involving magnetic resonance using nuclear magnetic resonance [NMR]
    • G01R33/48—NMR imaging systems
    • G01R33/4828—Resolving the MR signals of different chemical species, e.g. water-fat imaging
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
    • G01R33/00—Arrangements or instruments for measuring magnetic variables
    • G01R33/20—Arrangements or instruments for measuring magnetic variables involving magnetic resonance
    • G01R33/24—Arrangements or instruments for measuring magnetic variables involving magnetic resonance for measuring direction or magnitude of magnetic fields or magnetic flux
    • G01R33/243—Spatial mapping of the polarizing magnetic field
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00—Measuring for diagnostic purposes; Identification of persons
    • A61B5/05—Detecting, measuring or recording for diagnosis by means of electric currents or magnetic fields; Measuring using microwaves or radio waves 
    • A61B5/055—Detecting, measuring or recording for diagnosis by means of electric currents or magnetic fields; Measuring using microwaves or radio waves  involving electronic [EMR] or nuclear [NMR] magnetic resonance, e.g. magnetic resonance imaging
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
    • G01R33/00—Arrangements or instruments for measuring magnetic variables
    • G01R33/20—Arrangements or instruments for measuring magnetic variables involving magnetic resonance
    • G01R33/28—Details of apparatus provided for in groups G01R33/44 - G01R33/64
    • G01R33/38—Systems for generation, homogenisation or stabilisation of the main or gradient magnetic field
    • G01R33/385—Systems for generation, homogenisation or stabilisation of the main or gradient magnetic field using gradient magnetic field coils
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
    • G01R33/00—Arrangements or instruments for measuring magnetic variables
    • G01R33/20—Arrangements or instruments for measuring magnetic variables involving magnetic resonance
    • G01R33/44—Arrangements or instruments for measuring magnetic variables involving magnetic resonance using nuclear magnetic resonance [NMR]
    • G01R33/48—NMR imaging systems
    • G01R33/50—NMR imaging systems based on the determination of relaxation times, e.g. T1 measurement by IR sequences; T2 measurement by multiple-echo sequences
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
    • G01R33/00—Arrangements or instruments for measuring magnetic variables
    • G01R33/20—Arrangements or instruments for measuring magnetic variables involving magnetic resonance
    • G01R33/44—Arrangements or instruments for measuring magnetic variables involving magnetic resonance using nuclear magnetic resonance [NMR]
    • G01R33/48—NMR imaging systems
    • G01R33/54—Signal processing systems, e.g. using pulse sequences ; Generation or control of pulse sequences; Operator console
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
    • G01R33/00—Arrangements or instruments for measuring magnetic variables
    • G01R33/20—Arrangements or instruments for measuring magnetic variables involving magnetic resonance
    • G01R33/44—Arrangements or instruments for measuring magnetic variables involving magnetic resonance using nuclear magnetic resonance [NMR]
    • G01R33/48—NMR imaging systems
    • G01R33/54—Signal processing systems, e.g. using pulse sequences ; Generation or control of pulse sequences; Operator console
    • G01R33/56—Image enhancement or correction, e.g. subtraction or averaging techniques, e.g. improvement of signal-to-noise ratio and resolution
    • G01R33/5608—Data processing and visualization specially adapted for MR, e.g. for feature analysis and pattern recognition on the basis of measured MR data, segmentation of measured MR data, edge contour detection on the basis of measured MR data, for enhancing measured MR data in terms of signal-to-noise ratio by means of noise filtering or apodization, for enhancing measured MR data in terms of resolution by means for deblurring, windowing, zero filling, or generation of gray-scaled images, colour-coded images or images displaying vectors instead of pixels
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
    • G01R33/00—Arrangements or instruments for measuring magnetic variables
    • G01R33/20—Arrangements or instruments for measuring magnetic variables involving magnetic resonance
    • G01R33/44—Arrangements or instruments for measuring magnetic variables involving magnetic resonance using nuclear magnetic resonance [NMR]
    • G01R33/48—NMR imaging systems
    • G01R33/54—Signal processing systems, e.g. using pulse sequences ; Generation or control of pulse sequences; Operator console
    • G01R33/56—Image enhancement or correction, e.g. subtraction or averaging techniques, e.g. improvement of signal-to-noise ratio and resolution
    • G01R33/561—Image enhancement or correction, e.g. subtraction or averaging techniques, e.g. improvement of signal-to-noise ratio and resolution by reduction of the scanning time, i.e. fast acquiring systems, e.g. using echo-planar pulse sequences
    • G01R33/5613—Generating steady state signals, e.g. low flip angle sequences [FLASH]
    • G01R33/5614—Generating steady state signals, e.g. low flip angle sequences [FLASH] using a fully balanced steady-state free precession [bSSFP] pulse sequence, e.g. trueFISP
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T11/00—2D [Two Dimensional] image generation
    • G06T11/003—Reconstruction from projections, e.g. tomography
    • G06T11/005—Specific pre-processing for tomographic reconstruction, e.g. calibration, source positioning, rebinning, scatter correction, retrospective gating
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T11/00—2D [Two Dimensional] image generation
    • G06T11/003—Reconstruction from projections, e.g. tomography
    • G06T11/008—Specific post-processing after tomographic reconstruction, e.g. voxelisation, metal artifact correction
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00—Image analysis
    • G06T7/0002—Inspection of images, e.g. flaw detection
    • G06T7/0012—Biomedical image inspection
    • G06T7/0014—Biomedical image inspection using an image reference approach
    • G06T7/0016—Biomedical image inspection using an image reference approach involving temporal comparison
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
    • G01R33/00—Arrangements or instruments for measuring magnetic variables
    • G01R33/20—Arrangements or instruments for measuring magnetic variables involving magnetic resonance
    • G01R33/24—Arrangements or instruments for measuring magnetic variables involving magnetic resonance for measuring direction or magnitude of magnetic fields or magnetic flux
    • G01R33/246—Spatial mapping of the RF magnetic field B1
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
    • G01R33/00—Arrangements or instruments for measuring magnetic variables
    • G01R33/20—Arrangements or instruments for measuring magnetic variables involving magnetic resonance
    • G01R33/44—Arrangements or instruments for measuring magnetic variables involving magnetic resonance using nuclear magnetic resonance [NMR]
    • G01R33/48—NMR imaging systems
    • G01R33/4818—MR characterised by data acquisition along a specific k-space trajectory or by the temporal order of k-space coverage, e.g. centric or segmented coverage of k-space
    • G01R33/4824—MR characterised by data acquisition along a specific k-space trajectory or by the temporal order of k-space coverage, e.g. centric or segmented coverage of k-space using a non-Cartesian trajectory
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00—Indexing scheme for image analysis or image enhancement
    • G06T2207/10—Image acquisition modality
    • G06T2207/10072—Tomographic images
    • G06T2207/10088—Magnetic resonance imaging [MRI]
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00—Indexing scheme for image analysis or image enhancement
    • G06T2207/20—Special algorithmic details
    • G06T2207/20172—Image enhancement details
    • G06T2207/20182—Noise reduction or smoothing in the temporal domain; Spatio-temporal filtering

Landscapes

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

Abstract

Die Erfindung betrifft ein Magnetresonanz-Fingerprinting-Verfahren zur verbesserten Bestimmung von lokalen Parameterwerten eines Untersuchungsobjektes, bei dem mindestens zwei Signalvergleiche von erfassten Bildpunkt-Zeit-Serien mit Vergleichssignalverläufen zur Bestimmung von Parameterwerten durchgeführt werden, wobei ein weiterer (nachfolgender) Signalvergleich Ergebnisse eines vorhergehenden Signalvergleichs berücksichtigt. Die dadurch mehrstufige Bestimmung von Parameterwerten erlaubt eine Erhöhung der räumlichen Auflösung und der Genauigkeit mit der die Parameterwerte bestimmt werden können.

Description

  • Die Erfindung betrifft ein Magnetresonanz-Fingerprinting-Verfahren zur verbesserten Bestimmung von lokalen Parameterwerten eines Untersuchungsobjektes.
  • Die Magnetresonanz-Technik (im Folgenden steht die Abkürzung MR für Magnetresonanz) ist eine bekannte Technik, mit der Bilder vom Inneren eines Untersuchungsobjektes erzeugt werden können. Vereinfacht ausgedrückt wird hierzu das Untersuchungsobjekt in einem Magnetresonanzgerät in einem vergleichsweise starken statischen, homogenen Grundmagnetfeld, auch B0-Feld genannt, mit Feldstärken von 0,2 Tesla bis 7 Tesla und mehr positioniert, so dass sich dessen Kernspins entlang des Grundmagnetfeldes orientieren. Zum Auslösen von Kernspinresonanzen werden hochfrequente Anregungspulse (RF-Pulse) in das Untersuchungsobjekt eingestrahlt, die ausgelösten Kernspinresonanzen als sogenannte k-Raumdaten gemessen und auf deren Basis MR-Bilder rekonstruiert oder Spektroskopiedaten ermittelt. Zur Ortskodierung der Messdaten werden dem Grundmagnetfeld schnell geschaltete magnetische Gradientenfelder überlagert, die die Trajektorien festlegen, entlang derer die Messdaten im k-Raum ausgelesen werden. Die aufgezeichneten Messdaten werden digitalisiert und als komplexe Zahlenwerte in einer k-Raum-Matrix abgelegt. Aus der mit Werten belegten k-Raum-Matrix ist z.B. mittels einer mehrdimensionalen Fourier-Transformation ein zugehöriges MR-Bild rekonstruierbar. Eine hierzu verwendete auf bestimmte Art und Weise geordnete Folge von einzustrahlenden RF-Pulsen, zu schaltenden Gradienten und Auslesevorgängen wird als Sequenz bezeichnet.
  • Es sind verschiedene Sequenztypen bekannt, die unterschiedlich sensitiv auf die in einem untersuchten Untersuchungsobjekt enthaltenen Stoffe beschreibende Parameter (z.B. die longitudinale Longitudinalrelaxation T1, die Transversalrelaxation T2 und die Protonendichte) sind. Die aus mit einem bestimmten Sequenztypen aufgenommenen Messdaten rekonstruierten MR-Bilder zeigen entsprechend der Sensitivitäten des verwendeten Sequenztyps gewichtete Abbildungen des Untersuchungsobjektes.
  • Eine Magnetresonanz-Bildgebung mittels einer Magnetresonanzanlage kann dazu dienen, eine Anwesenheit und/oder eine Verteilung eines Stoffs, welcher sich in einem Untersuchungsobjekt befindet, zu bestimmen. Der Stoff kann dabei beispielsweise ein, möglicherweise pathologisches, Gewebe des Untersuchungsobjekts, ein Kontrastmittel, eine Markierungssubstanz oder ein Stoffwechselprodukt sein.
  • Informationen über die vorliegenden Stoffe können dabei auf vielfältige Weise aus den aufgenommenen Messdaten gewonnen werden. Eine relativ einfache Informationsquelle sind z.B. aus den Messdaten rekonstruierte Bilddaten. Es gibt jedoch auch komplexere Verfahren, die, z.B. aus Bildpunkt-Zeit-Serien von aus sukzessive gemessenen Messdatensätzen rekonstruierten Bilddaten, Informationen über das untersuchte Untersuchungsobjekt ermitteln.
  • Mit Hilfe von quantitativen MR-Bildgebungstechniken lassen sich absolute Eigenschaften des gemessenen Objekts bestimmen, am Menschen z.B. die gewebespezifische T1- und T2-Relaxation. Im Gegensatz dazu erzeugen die in der klinischen Routine meist verwendeten konventionellen Sequenzen lediglich eine relative Signalintensität unterschiedlicher Gewebetypen (sogenannte Wichtungen), so dass die diagnostische Interpretation in hohem Maße der subjektiven Einschätzung des Radiologen unterliegt. Quantitative Techniken bieten somit den offensichtlichen Vorteil einer objektiven Vergleichbarkeit, finden aber aufgrund langer Messzeiten derzeit in der Routine kaum Verwendung.
  • Neuere quantitative Messverfahren wie Magnetresonanz-Fingerprinting-Verfahren (MRF-Verfahren), könnten den genannten Nachteil der langen Messzeiten auf ein akzeptables Maß senken. Bei MRF-Verfahren werden Signalverläufe von aus zeitlich nacheinander mit verschiedenen Aufnahmeparametern aufgenommenen Messdaten rekonstruierten Bilddaten für jeden Bildpunkt oder zumindest interessierende Bildpunkte der Bilddaten als Bildpunkt-Zeit-Serien betrachtet, wobei der Signalverlauf einer Bildpunkt-Zeit-Serie als „Fingerabdruck“ (engl. „fingerprint“) der in dem Bildpunkt Ort des Untersuchungsobjektes während der Messung vorliegenden Parameter betrachtet wird. Diese Signalverläufe werden mittels Mustererkennungsverfahren („pattern recognition“) mit Signalverläufen einer zuvor ermittelten Datenbank von für bestimmte Stoffe charakteristischen Signalverläufen (dem sogenannten „Dictionary“) verglichen, um die in den aus den Messdaten rekonstruierten Bilddaten repräsentierten Stoffe bzw. die räumliche Verteilung von gewebespezifischen Parametern (wie die Transversalrelaxation T2 oder die Longitudinalrelaxation T1; sogenannte T1- und T2-Karten) in dem abgebildeten Untersuchungsobjekt zu ermitteln. Die in einem derartigen Dictionary enthaltenen Signalverläufe können hierbei auch durch Simulationen erstellt worden sein.
  • Das Prinzip dieser Methode ist somit, gemessene Signalverläufe mit einer Vielzahl im Voraus bekannter Signalverläufe zu vergleichen. Dabei können Signalverläufe für verschiedene Kombinationen aus T1- und T2-Relaxationszeiten sowie auch anderen Parametern für das Dictionary ermittelt worden sein. Die Parameterwerte, z.B. T1- und T2-Zeiten, eines Bildpunktes (Pixels/Voxels) im Bild werden dann insbesondere bestimmt, indem der gemessene Signalverlauf mit allen oder einem Teil der simulierten Signalverläufe verglichen wird. Dieser Vorgang wird als „Matching“ bezeichnet. Derjenige Signalverlauf des Dictionarys, der dem gemessenen Signalverlauf am ähnlichsten ist, determiniert die Parameter, z.B. Relaxationsparameter T1 und T2, des jeweiligen Bildpunktes.
  • Grundsätzlich können dabei neben den bereits genannten gewebespezifischen Parametern eines untersuchten Objektes auch messungsspezifische Parameter, wie z.B. die Feldstärken der applizierten Magnetfelder oder auch die lokale Verteilung der Stärke eines eingestrahlten Hochfrequenzfeldes B1+, ermittelt werden, da mittels MR-Techniken aufgenommene Signale von den in einem untersuchten Objekt vorliegenden gewebespezifischen Parametern sowie von messungsspezifischen Parametern, die die während der Messung vorliegenden Bedingungen beschreiben, abhängen können. Die verwendeten Aufnahmeparameter sind hierbei derart gewählt, dass die aufgenommenen Messdaten eine Abhängigkeit von den gewünschten zu bestimmenden Parametern zeigen. Beispielsweise können Sequenztypen für das MRF-Verfahren eingesetzt werden, die sensitiv auf die gewünschten Parameter sind. Durch die Abhängigkeiten und die Variation der Aufnahmeparameter und deren Berücksichtigung in den Vergleichssignalverläufen sind die gewünschten Parameter aus derartig aufgenommenen Bildpunkt-Zeit-Serien bestimmbar.
  • Ein gängiges MRF-Verfahren ist beispielsweise in dem Artikel von Ma et al., „Magnetic Resonance Fingerprinting", Nature, 495: S. 187-192 (2013) beschrieben. Dort wird eine TrueFISP („true fast imaging with steady-state free precession“) basierte Sequenz in Kombination mit einer spiralförmigen k-Raumabtastung benutzt. Auf Grund der Abhängigkeit eines mittels einer TrueFISP-Sequenz aufgenommenen Signals vom statischen Magnetfeld B0 ist dieser Parameter ebenfalls Teil des das Signal beschreibenden Signalmodells und kann, wenn die Bildpunkt-Zeit-Serie unter entsprechender Variation der Aufnahmeparameter erstellt wurde, grundsätzlich durch einen MRF-„Dictionary-Match“ neben T1 und T2 in jedem Pixel identifiziert werden. Die Stärke eines aufgenommenen Signals hängt bei TrueFISP-Sequenzen maßgeblich vom lokalen Grundmagnetfeld B0 ab, welches innerhalb eines Untersuchungsobjektes räumlich variiert. Durch Variation der Wiederholzeit TR kann die B0-Abhängigkeit des Signals verändert werden. Allerdings wird bei MRF-Verfahren (vgl. den bereits genannten Artikel von Ma et al.) eine relativ lange Auslesedauer für die spiralförmige k-Raumabtastung benutzt, die möglichen Dauern einer Wiederholzeit TR nach unten limitiert, sodass eine vergleichsweise hohe untere Grenze für die Wiederholzeit TR existiert. Daher kann es zu Artefakten in MRF-Messungen kommen, die sogenannten Banding-Artefakten in konventionellen TrueFISP Sequenzen sehr ähnlich sind (vgl. den bereits genannten Artikel von Ma et al.). Des Weiteren ist eine erreichbare räumliche Auflösung, mit der der messungsspezifische Parameter Grundmagnetfeld B0 mit einer solchen MRF-Technik bestimmt werden kann, wegen der beschriebenen hohen unteren Grenze für die Wiederholzeit TR ebenfalls limitiert. Diese Limitierung der räumlichen Auflösung führt dazu, dass die auf diese Weise bestimmten räumlichen Verteilungen des Grundmagnetfeldes B0 (B0-Karten) für manche Anwendungen, z.B. Korrekturen von aus den lokalen Abweichungen des Grundmagnetfeldes von einem globalen Wert resultierenden Artefakten, nicht ausreichend hoch aufgelöst sind.
  • Für MRF-Verfahren kann grundsätzlich jede Echotechnik in Kombination mit jeglichen Verfahren zur k-Raumabtastung (z.B. kartesisch, spiralförmig, radial) benutzt werden. Nichtkartesische k-Raumabtastungen, wie insbesondere radiale oder spiralförmige k-Raumabtastungen, decken den zentralen Bereich des k-Raums inhärent stärker ab als periphere Bereiche des k-Raums, was für Messungen mit einer Unterabtastung (gemäß dem Theorem von Nyquist) des k-Raums vorteilhaft ist. Ortskodierungen mit spiralförmigen Trajektorien bei der k-Raumabtastung, wie sie z.B. in dem genannten Artikel von Ma et al. verwendet werden, und wie ein Beispiel in der 1 gezeigt ist, werden wegen ihrer hohen Effizienz, insbesondere ihrer Eigenschaft eine starke Unterabtastung des k-Raums bezüglich des Theorems von Nyquist zu erlauben, die besonders kurze Messzeiten ermöglicht, besonders häufig für MRF-Verfahren eingesetzt. Allerdings sind derartige spiralförmige Ortskodierungen besonders empfindlich auf Off-Resonanzen, d.h. lokale Störungen des statischen Magnetfeldes B0, wodurch sich in den aus den entlang spiralförmig ausgelesenen Messdaten rekonstruierten Bildern unscharfe Bereiche ergeben können, sogenannte Blurring-Artefakte.
  • Der Erfindung liegt die Aufgabe zugrunde, eine verbesserte Bestimmung von Parameterwerten mittels MRF-Verfahren zu ermöglichen.
  • Die Aufgabe wird gelöst durch ein Verfahren zur Bestimmung von Parameterwerten in Bildpunkten eines Untersuchungsobjektes mittels einer Magnetresonanz-Fingerprinting(MRF)-Technik gemäß Anspruch 1, eine Magnetresonanzanlage gemäß Anspruch 12, ein Computerprogramm gemäß Anspruch 13, sowie einen elektronisch lesbaren Datenträger gemäß Anspruch 14.
  • Ein erfindungsgemäßes Verfahren zur Bestimmung von Parameterwerten in Bildpunkten eines Untersuchungsobjektes mittels einer Magnetresonanz-Fingerprinting(MRF)-Technik umfasst die Schritte:
    • - Erfassen mindestens einer Bildpunkt-Zeit-Serie mit Hilfe eines MRF-Aufnahmeverfahrens, aus denen je ein Wert zumindest eines in dem jeweiligen Bildpunkt dargestellten Ort des Untersuchungsobjekts vorliegenden Parameters bestimmt werden soll,
    • - Durchführen eines Signalvergleichs mindestens eines Abschnitts des jeweiligen Signalverlaufs der erfassten Bildpunkt-Zeit-Serien mit einem entsprechenden Abschnitt von Vergleichssignalverläufen zur Bestimmung der jeweiligen Werte des mindestens einen Parameters,
    • - Durchführen eines weiteren Signalvergleichs mindestens eines Abschnittes des jeweiligen Signalverlaufs unter Berücksichtigung der Ergebnisse des vorhergehenden Signalvergleichs mit einem entsprechenden Abschnitt der Vergleichssignalverläufe zur verbesserten Bestimmung der jeweiligen Werte des mindestens einen Parameters,
    • - Ausgeben der für den jeweiligen Bildpunkt verbessert bestimmten jeweiligen Werte der zu bestimmenden Parameter.
  • Durch das erfindungsgemäße Durchführen eines weiteren Signalvergleichs unter Berücksichtigung eines vorhergehenden Signalvergleichs und die damit einhergehende schrittweise bzw. mehrstufige Bestimmung von Parameterwerten können die zu bestimmenden Parameter mit einer höheren räumlichen Auflösung und mit einer höheren Genauigkeit bestimmt werden.
  • Hierbei kann auf Grundlage eines für einen Bildpunkt mit einem erfindungsgemäßen Signalvergleich bestimmten Wertes eines Parameters die dem Bildpunkt entsprechende erfasste Bildpunkt-Zeit-Serie korrigiert werden, sodass eine korrigierte Bildpunkt-Zeit-Serie erstellt wird.
  • Insbesondere kann bereits durch den ersten Signalvergleich mindestens eines Abschnitts des Signalverlaufs der erfassten Bildpunkt-Zeit-Serie mit einem entsprechenden Abschnitt der Vergleichssignalverläufe („dictionary“) ein Wert eines Parameters, insbesondere eines messungsspezifischen Parameters wie das Grundmagnetfeld B0, bestimmt werden, und die erfasste Bildpunkt-Zeit-Serie auf Grundlage des in dem ersten Signalvergleich bestimmten Wertes des Parameters korrigiert werden, sodass eine korrigierte Bildpunkt-Zeit-Serie erstellt wird.
  • Korrekturen von Bilddaten auf Grundlage von insbesondere messungsspezifischen Parametern sind bekannt. Beispielsweise werden Einflüsse von Inhomogenitäten der verwendeten Magnetfelder und/oder Hochfrequenzfelder häufig retrospektiv nach Ermittlung der lokalen Verteilung der jeweiligen Magnet- oder Hochfrequenzfelder korrigiert. Solche Korrekturverfahren sind auch hier anwendbar und können Artefakte in den Bildpunkt-Zeit-Serien beseitigen. Somit können sie die Genauigkeit und Artefaktfreiheit der einem weiteren Signalvergleich zugrunde liegenden korrigierten Bildpunkt-Zeit-Serien erhöhen, was sich positiv auf die Genauigkeit der bestimmten Werte der Parameter auswirkt.
  • Wurden in einem vorhergehenden Signalvergleich Werte mehrerer verschiedener Parameter bestimmt, kann die Bildpunkt-Zeit-Serie, z.B. sukzessive, auf Grundlage mehrerer der bestimmten Werte korrigiert werden.
  • Die Ergebnisse eines vorhergehenden Signalvergleichs können bei einem weiteren Signalvergleich berücksichtigt werden, indem der weitere Signalvergleich eine derartig auf Grundlage von in einem vorhergehenden Signalvergleich bestimmten Parameterwerten korrigierte Bildpunkt-Zeit-Serie mit den Vergleichssignalverläufen vergleicht. Es kann also zur Berücksichtigung der Ergebnisse des vorhergehenden Signalvergleichs bei einem weiteren Signalvergleich der weitere Signalvergleich mindestens einen Abschnitt des Signalverlaufs der korrigierten Bildpunkt-Zeit-Serie mit einem entsprechenden Abschnitt der Vergleichssignalverläufe zur verbesserten Bestimmung der Werte der zu bestimmenden Parameter vergleichen.
  • Durch die Verwendung von korrigierten Bildpunkt-Zeit-Serien bei dem weiteren Signalvergleich kann eine höhere Genauigkeit der verbessert bestimmten Parameterwerte erreicht werden.
  • Eine korrigierte Bildpunkt-Zeit-Serie kann auf Grundlage eines durch Signalvergleich der korrigierten Bildpunkt-Zeit-Serie mit den Vergleichssignalverläufen verbessert bestimmten Wertes des Parameters erneut korrigiert werden, und erneut der Wert des Parameters durch Vergleichen der erneut korrigierten Bildpunkt-Zeit-Serie mit den Vergleichssignalverläufen verbessert bestimmt werden. Somit kann eine solche Korrektur einer Bildpunkt-Zeit-Serie auch mehrfach, z.B. iterativ, durchgeführt werden.
  • Durch die erneute Korrektur kann eine höhere Genauigkeit sowohl der korrigierten Bildpunkt-Zeit-Serie als auch der verbessert bestimmten Parameterwerte erreicht werden.
  • Eine derartige erneute Korrektur kann solange, z.B. iterativ, durchgeführt werden, bis ein Unterschied von erneut verbessert bestimmten Werten des Parameters zu den entsprechenden vor der erneuten Korrektur bestimmten Werten des Parameters kleiner als ein vorgegebener Schwellwert ist.
  • Der Vergleich einer Differenz der in verschiedenen Signalvergleichen bestimmten Werte eines zu bestimmenden Parameters und Vergleich dieser Differenz mit einem vorgegebenen Schwellwert gewährleisten ein Erreichen einer durch den Schwellwert definierten Genauigkeit des jeweiligen Wertes.
  • Es ist weiterhin möglich, dass Werte von mindestens zwei verschiedenen Parametern bestimmt werden, und zur Berücksichtigung der Ergebnisse des vorhergehenden Signalvergleichs bei Signalvergleichen zur verbesserten Bestimmung der Werte mindestens eines zu bestimmenden Parameters, der bzw. die Werte mindestens eines anderen als des bzw. der aktuell verbessert zu bestimmenden Parameter der zu bestimmenden Parameter auf eine vorgegebene maximale Schwankung gegenüber dem Ergebnis des vorhergehenden Signalvergleichs begrenzt, insbesondere konstant gehalten, werden.
  • Durch ein derartiges Begrenzen einer erlaubten Schwankung von Werten für einige Parameter während einer verbesserten Bestimmung eines Wertes mindestens eines anderen Parameters kann der Rechenaufwand bei der verbesserten Bestimmung gering gehalten werden und der Signalvergleich für den verbessert zu bestimmenden Parameterwert optimiert durchgeführt werden, sodass die Genauigkeit der verbessert bestimmten Parameterwerte erhöht ist.
  • Hierbei können solange Signalvergleiche zur verbesserten Bestimmung eines zu bestimmenden Parameters, durchgeführt werden, bis die Werte aller insgesamt zu bestimmender Parameter mindestens einmal bei der verbesserten Bestimmung bestimmt wurden oder bis alle Werte zu bestimmender Parameter mit ausreichender Genauigkeit, z.B. gemäß einer oben beschriebenen Überprüfung mit einem vorgegebenen Schwellwert, bestimmt wurden.
  • Eine erfindungsgemäße Magnetresonanzanlage umfasst eine Magneteinheit, eine Gradienteneinheit, eine Hochfrequenzeinheit und eine zur Durchführung eines erfindungsgemäßen Verfahrens ausgebildete Steuereinrichtung mit einer Parameterwertbestimmungseinheit.
  • Ein erfindungsgemäßes Computerprogramm implementiert ein erfindungsgemäßes Verfahren auf einer Steuereinrichtung, wenn es auf der Steuereinrichtung ausgeführt wird.
  • Das Computerprogramm kann hierbei auch in Form eines Computerprogrammprodukts vorliegen, welches direkt in einen Speicher einer Steuereinrichtung ladbar ist, mit Programmcode-Mitteln, um ein erfindungsgemäßes Verfahren auszuführen, wenn das Computerprogrammprodukt in der Recheneinheit des Rechensystems ausgeführt wird.
  • Ein erfindungsgemäßer elektronisch lesbarer Datenträger umfasst darauf gespeicherte elektronisch lesbare Steuerinformationen, welche zumindest ein erfindungsgemäßes Computerprogramm umfassen und derart ausgestaltet sind, dass sie bei Verwendung des Datenträgers in einer Steuereinrichtung einer Magnetresonanzanlage ein erfindungsgemäßes Verfahren durchführen.
  • Die in Bezug auf das Verfahren angegebenen Vorteile und Ausführungen gelten analog auch für die Magnetresonanzanlage, das Computerprogrammprodukt und den elektronisch lesbaren Datenträger.
  • Weitere Vorteile und Einzelheiten der vorliegenden Erfindung ergeben sich aus den im Folgenden beschriebenen Ausführungsbeispielen sowie anhand der Zeichnungen. Die aufgeführten Beispiele stellen keine Beschränkung der Erfindung dar. Es zeigen:
    • 1 eine beispielhafte spiralförmige k-Raum-Trajektorie,
    • 2 ein schematisches Ablaufdiagramm eines erfindungsgemäßen Verfahrens,
    • 3 ein mögliches Aufnahmeschema für eine Bildpunkt-Zeit-Serie,
    • 4 Darstellungen möglicher räumlicher Verteilungen eines bestimmten Parameters zur Veranschaulichung eines Ausgleichsprozesses,
    • 5 eine schematisch dargestellte erfindungsgemäße Magnetresonanzanlage.
  • 1 zeigt gestichelt dargestellt eine beispielhafte spiralförmige k-Raum-Trajektorie entlang derer nach einem RF-Anregungspuls ein Satz an Messdaten ausgelesen wird, aus dem ein Bilddatensatz rekonstruiert wird. Zur Erzeugung einer derartigen k-Raumtrajektorie werden nach dem Einstrahlen eines RF-Anregungspulses in das Untersuchungsobjekt Gradientenfelder in zwei senkrecht zueinander stehenden Phasenkodierrichtungen (hier x- und y-Richtung) derart auf grundsätzlich bekannte Art und Weise geschaltet, dass die spiralförmige k-Raumtrajektorie abgetastet wird. Durch wiederholte Aufnahme derartiger Messdatensätze und Rekonstruktion zugehöriger Bilddatensätze entsteht eine Serie von Bilddatensätzen, in denen ein Signalverlauf in jeweils einem gleichen der Bildpunkte der Bilddatensätze in MRF-Verfahren als Bildpunkt-Zeit-Serie betrachtet und mit Vergleichssignalverläufen verglichen werden kann, um in dem betrachteten Bildpunkt in dem Untersuchungsobjekt zum Zeitpunkt der Messung vorliegende Parameter zu bestimmen.
  • 2 ist ein schematisches Ablaufdiagramm eines erfindungsgemäßen Verfahrens zur Bestimmung von Parameterwerten in Bildpunkten eines Untersuchungsobjektes mittels einer MRF-Technik.
  • Dabei wird mindestens eine Bildpunkt-Zeit-Serie BZS mit Hilfe eines MRF-Aufnahmeverfahrens, aus denen je ein Wert P zumindest eines in dem jeweiligen Bildpunkt dargestellten Ortes des Untersuchungsobjekts vorliegenden Parameters bestimmt werden soll, erfasst (Block 101). Genauer werden typischerweise bei MRF-Verfahren Serien von Bildern in einer gewissen zeitlichen Folge unter Variation von Aufnahmeparametern aufgenommen. Jeder den Bildern der Serie von Bildern gemeinsame Bildpunkt einer solchen Serie von Bildern kann als Bildpunkt-Zeit-Serie BZS betrachtet werden. Es können alle Bildpunkte der Serie von Bildern als jeweilige Bildpunkt-Zeit-Serien BZS betrachtet werden oder auch nur interessierende Bildpunkte, beispielsweise diejenigen, welche einen Untersuchungsbereich (ROI, „region of interest“) darstellen. Ein bei der Erfassung der Bildpunkt-Zeit-Serien verwendetes MRF-Aufnahmeverfahren legt die bestimmbaren Parameter fest. Vorteilhafterweise wird ein MRF-Aufnahmeverfahren verwendet, das die Bestimmung mehrerer, insbesondere gewebespezifischer und auch messungsspezifischer, Parameter erlaubt. Ein Beispiel für ein geeignetes MRF-Aufnahmeschema ist weiter unten mit Bezug auf 3 beschrieben. Generell sind MRF-Aufnahmetechniken bevorzugt, bei denen zumindest die Parameter B0, T1 und T2, gerne auch noch weitere Parameter wie B1+, durch einen Signalvergleich mittels eines MRF-Matching-Verfahrens bestimmt werden können.
  • Ein Signalvergleich wird durchgeführt (Block 103), der den jeweiligen Signalverlauf zumindest eines Abschnitts einer erfassten Bildpunkt-Zeit-Serie mit einem entsprechenden Abschnitt von Vergleichssignalverläufen D vergleicht, um einen jeweiligen Wert P des mindestens einen zu bestimmenden Parameters in dem jeweiligen Bildpunkt zu bestimmen. Die Bestimmung der Werte P der zu bestimmenden Parameter kann dabei mit Hilfe eines bei MRF-Verfahren üblichen „Matching“-Verfahrens mit in einem Dictionary enthaltenen Vergleichssignalverläufen D erfolgen.
  • Ein weiterer Signalverlauf (Block 105) wird durchgeführt, der mindestens einen Abschnitt eines jeweiligen Signalverlaufs unter Berücksichtigung der Ergebnisse des vorhergehenden Signalvergleichs 103 oder, wenn bereits zumindest einmal erfolgt, eines Signalverlaufs 105 mit einem entsprechenden Abschnitt der Vergleichssignalverläufe D vergleicht, um einen jeweiligen Wert P' des mindestens einen Parameters in dem jeweiligen Bildpunkt verbessert zu bestimmen.
  • Die verbessert bestimmten jeweiligen Werte P' der zu bestimmenden Parameter können gespeichert und/oder an einer Ausgabe ausgegeben werden, z.B., insbesondere in Form von Parameterkarten, angezeigt werden (Block 107).
  • Somit werden die Werte P' von zu bestimmenden Parametern schrittweise in zumindest zwei Signalvergleichen 103 und 105 verbessert bestimmt.
  • Bei den verwendeten Signalverläufen kann bei mindestens einem Signalvergleich der Signalverlauf einer gesamten Bildpunkt-Zeit-Serie BZS mit dem gesamten Verlauf der Vergleichssignalverläufe D verglichen werden.
  • Bei den verwendeten Signalverläufen kann auch bei mindestens einem Signalvergleich der Signalverlauf in Abschnitten einer gesamten zu vergleichenden Bildpunkt-Zeit-Serie BZS, die mit einem auf den zu bestimmenden Parameter sensitiven Sequenztyp erfasst wurden, mit dem entsprechenden Abschnitten der Vergleichssignalverläufe D verglichen werden. Dabei bleiben zwar diejenigen Signale der nicht verglichenen Abschnitte der Signalverläufe unberücksichtigt, es kann jedoch erreicht werden, dass mögliche, in diesen unberücksichtigten Abschnitten enthaltene, einen zu bestimmenden Parameter möglicherweise verfälschende Artefakte keinen oder zumindest weniger Einfluss auf die bestimmten Werte der zu bestimmenden Parameter haben, insbesondere wenn die gewählten Abschnitte nicht mit den zu vermeidenden Artefakten behaftet sind.
  • Es kann durch einen Signalvergleich 103 (oder auch 105) mindestens eines Abschnitts des Signalverlaufs einer erfassten Bildpunkt-Zeit-Serie BZS mit einem entsprechenden Abschnitt der Vergleichssignalverläufe D ein Wert P (oder auch P') eines Parameters bestimmt werden, der mit Hilfe von aus zumindest Teilen der erfassten Bildpunkt-Zeit-Serie BZS erhaltenen Ausgleichsdaten ausgeglichen wird (Block 109), sodass ein ausgeglichener Wert P* des Parameters bestimmt wird, der den zuvor bestimmten Wert P (oder P') des Parameters ersetzen kann.
  • Insbesondere kann, wie weiter unten mit Bezug auf 4 beschrieben wird, durch den Signalvergleich für jeden Bildpunkt ein Wert eines lokal vorliegenden Grundmagnetfeldes B0 bestimmt werden, und aus den Signalen der erfassten Bildpunkt-Zeit-Serie Ausgleichsdaten erhalten werden, die z.B. in Form von grob aufgelösten Parameterwerten, insbesondere Grundmagnetfeldwerten, ermittelt werden können, und mit Hilfe derer als Ausgleich Artefakte in den Parameterwerten entfernt oder zumindest reduziert werden, wodurch die somit ausgeglichenen Parameterwerte artefaktfrei oder zumindest weniger mit Artefakten behaftet sind und ein im Gegensatz zu dem Wertebereich der noch nicht ausgeglichenen Parameter vergrößerter Wertebereich der ausgeglichenen Parameterwerte erzielt werden kann.
  • Auf Grundlage von mit einem erfindungsgemäßen Signalvergleich 103 oder 105 bestimmten Werten P oder P' eines Parameters in einem Bildpunkt kann die zu diesem Bildpunkt erfasste Bildpunkt-Zeit-Serie korrigiert werden (Block 111), sodass eine korrigierte Bildpunkt-Zeit-Serie BZS' erstellt wird.
  • Insbesondere kann durch den ersten Signalvergleich 103 mindestens eines Abschnitts des Signalverlaufs einer erfassten Bildpunkt-Zeit-Serie BZS mit einem entsprechenden Abschnitt der Vergleichssignalverläufe D ein Wert P eines messungsspezifischen Parameters wie z.B. das Grundmagnetfeld B0, bestimmt werden. Auf Grundlage des bestimmten Wertes P des Parameters kann die erfasste Bildpunkt-Zeit-Serie BZS zu einer korrigierten Bildpunkt-Zeit-Serie BZS' korrigiert werden.
  • Somit ist es möglich die für das MRF-Verfahren erfassten Bildpunkt-Zeit-Serien BZS mittels aus den Bildpunkt-Zeit-Serien selbst abgeleiteten Parameterwerten zu korrigieren.
  • Werden die Bildpunkt-Zeit-Serien beispielsweise gemäß einer spiralförmigen Abtastung des k-Raums erfasst, können die Bildpunkt-Zeit-Serien wie oben beschrieben durch Off-Resonanzen mit Blurring-Artefakten behaftet sein. Diese können mit dem beschriebenen Verfahren auf besonders effiziente Weise durch lokale BO-Werte, die durch einen Signalvergleich aus den Bildpunkt-Zeit-Serien selbst bestimmt wurden, korrigiert werden.
  • Hierzu sind bereits Off-Resonanz-Korrekturmethoden (auch „Deblurring“-Methoden), insbesondere für spiralförmige k-Raumabtastungen, bekannt. Diese erforderten in der Regel bisher eine separate Messung einer BO-Karte. Mit dem hier beschriebenen Verfahren kann eine derartige separate Messung einer BO-Karte entfallen und direkt die aus dem Signalvergleich bestimmten Werte für das Grundmagnetfeld B0 verwendet werden. Anhand der gewonnenen Information über die räumliche Verteilung des Grundmagnetfeldes B0 können die mit der spiralförmigen k-Raumabtastung gemessenen Daten derart manipuliert werden, dass das genannte Blurring (Verschmieren) nicht mehr auftritt. Gängig ist beispielsweise eine frequenzsegmentierte Methode, bei der eine bisher separat gemessene B0-Karte nach n Off-Resonanzen segmentiert wird. Die mit der spiralförmigen k-Raumabtastung gemessenen Daten werden unter Einbeziehung verschiedener Off-Resonanzen rekonstruiert, d.h. n Bilder für n Off-Resonanzen werden rekonstruiert. Anhand der derart segmentierten BO-Karte wird für jeden Bildpunkt ausgewählt, welches der n Bilder den korrekt rekonstruierten Bildpunkt enthält. Es können auch andere bekannte Korrekturverfahren eingesetzt werden, die eine Korrektur von Bilddaten auf Basis von Parameterkarten durchführen, wobei die benötigten Parameterkarten im Rahmen des hier beschriebenen Verfahrens nicht separat bestimmt, sondern direkt aus den Bildpunkt-Zeit-Serien durch mindestens einen Signalvergleich mit Vergleichssignalverläufen bestimmt werden können.
  • Zur Berücksichtigung der Ergebnisse eines vorhergehenden Signalvergleichs 103 oder 105 bei einem weiteren Signalvergleich 105 kann der weitere Signalvergleich mindestens einen Abschnitt des Signalverlaufs einer derart korrigierten Bildpunkt-Zeit-Serie BZS' mit einem entsprechenden Abschnitt der Vergleichssignalverläufe D zur verbesserten Bestimmung der Werte P' der zu bestimmenden Parameter verglichen werden.
  • Durch die Korrektur 111 der Bildpunkt-Zeit-Serie BZS und erneute Bestimmung der Werte P' der zu bestimmenden Parameter kann die räumliche Auflösung und auch die Genauigkeit der verbessert bestimmten Werte P' erhöht werden.
  • Eine derartige Korrektur 111 kann mehrmals durchgeführt werden, sodass eine korrigierte Bildpunkt-Zeit-Serie BZS' auf Grundlage eines durch Signalvergleich der korrigierten Bildpunkt-Zeit-Serie BZS' mit den Vergleichssignalverläufen D bestimmten Wertes P' des Parameters erneut korrigiert wird, und erneut der Wert P' des Parameter durch Vergleichen der erneut korrigierten Bildpunkt-Zeit-Serie BZS' mit den Vergleichssignalverläufen D bestimmt wird.
  • Hierzu kann eine Abfrage 115 beispielsweise prüfen, ob eine gewünschte Anzahl an Korrekturen durchgeführt wurde. Ist dies der Fall (Abfrage 115, y), können die nach der letzten Korrektur der Bildpunkt-Zeit-Serie BZS bzw. BZS' bestimmten Werte P' der zu bestimmenden Parameter gespeichert und/oder ausgegeben werden (Block 107). Ist dies nicht der Fall (Abfrage 115, n), kann eine erneute Korrektur der Bildpunkt-Zeit-Serie BZS oder BZS` durchgeführt werden.
  • Es ist auch denkbar, solange erneute Korrekturen (iterativ) durchzuführten, bis ein Unterschied von erneut bestimmten Werten P' des Parameters zu den entsprechenden vor der erneuten Korrektur bestimmten Werten P oder P' des Parameters kleiner als ein vorgegebener Schwellwert S ist.
  • Hierzu kann die Abfrage 115 zusätzlich oder alternativ die Differenz der nach den beiden letzten Korrekturen bestimmten Werte P,P` bilden und mit dem vorgegebenen Schwellwert S vergleichen. Ist die genannte Bedingung erfüllt (Abfrage 115, y), werden die nach der letzten Korrektur der Bildpunkt-Zeit-Serie BZS bzw. BZS' bestimmten Werte P' der zu bestimmenden Parameter gespeichert und/oder ausgegeben (Block 107). Ist sie nicht erfüllt (Abfrage 115, n), kann eine erneute Korrektur der Bildpunkt-Zeit-Serie BZS oder BZS' durchgeführt werden.
  • Somit erlaubt das beschriebene Verfahren eine, auch iterativ mögliche, Korrektur und damit Vermeidung von Artefakten, z.B. von Off-Resonanz-Artefakten („Deblurring“), in für eine Bestimmung von Parameterwerten in MRF-Verfahren verwendeten Bildpunkt-Zeit-Serien und somit ebenfalls eine Korrektur und damit Vermeidung derartiger Artefakte in den bestimmten Parameterwerten.
  • Sollen die lokalen Werte P' von mindestens zwei verschiedenen Parametern bestimmt werden, können die Werte der insgesamt zu bestimmenden Parameter in zwei Gruppen G1, G2 eingeteilt werden, von denen nur Werte der Gruppe G1 verbessert bestimmt werden, wohingegen die Werte der anderen Gruppe G2 derart beschränkt werden, dass sie nur gemäß einer vorgegebenen Schwankung SW von dem in einem vorhergehenden Signalvergleich bestimmten Wert abweichen dürfen (Block 113). Es kann somit zur Berücksichtigung der Ergebnisse eines vorhergehenden Signalvergleichs bei Signalvergleichen zur verbesserten Bestimmung der Werte P' mindestens eines zu bestimmenden Parameters G1, der bzw. die Werte P, P' mindestens eines anderen (G2) als des bzw. der aktuell verbessert zu bestimmenden Parameter G1 auf eine vorgegebene maximale Schwankung SW gegenüber dem Ergebnis des vorhergehenden Signalvergleichs begrenzt, insbesondere konstant gehalten, werden (Block 113). Somit können in einem weiteren Signalvergleich 105 nur Werte von verbessert zu bestimmenden Parametern G1 verbessert bestimmt werden, während die Werte P2 der übrigen Parameter allenfalls gering von einem in einem vorhergehenden Signalvergleich bestimmten Wert abweichen dürfen. Auf diese Weise kann eine dedizierte verbesserte Bestimmung der Werte P' einzelner oder weniger Parameter G1 erreicht werden.
  • Es können solange Signalvergleiche zur verbesserten Bestimmung eines Wertes P' eines zu bestimmenden Parameters, durchgeführt werden, bis alle Werte P' der zu bestimmenden Parameter mindestens einmal bei der verbesserten Bestimmung bestimmt wurden.
  • Hierzu kann die Abfrage 115 zusätzlich oder alternativ überwachen, ob zu allen zu bestimmenden Parameter Werte P' verbessert bestimmt wurden, und ggf. ob alle dediziert verbessert zu bestimmenden Werte von Parametern G1 bestimmt wurden. Ist jede gewünschte der genannten Bedingungen erfüllt (Abfrage 115, y) werden verbessert bestimmten Werte P' der zu bestimmenden Parameter gespeichert und/oder ausgegeben (Block 107). Ist dies nicht der Fall (Abfrage 115, n) kann ein weiterer Signalvergleich zur verbesserten Bestimmung eines noch nicht verbessert bestimmten Parameters durchgeführt werden.
  • 3 zeigt schematisch ein Aufnahmeverfahren zur Aufnahme von mindestens einer Bildpunkt-Zeit-Serie, wie sie für das erfindungsgemäße Verfahren eingesetzt werden kann. Das dargestellte Beispiel zeigt ein Aufnahmeverfahren zur Erfassung von Bildpunkt-Zeit-Serien, bei welchem drei verschiedene Sequenztypen mit Eigenschaften von Sequenztypen aus der Gruppe der Sequenztypen TrueFISP („True Fast Imaging with Steady-state free Precession“), FISP („Fast Imaging with Steady-state free Precession“) und FLASH („Fast Low-Angle Shot“) eingesetzt werden. Die Eigenschaften der Sequenztypen betreffen insbesondere die Sensitivität der jeweiligen Sequenztypen auf Variation von gewebespezifischen und/oder messungsspezifischen Parametern. Beispielsweise ist eine FISP-Sequenz wenig sensitiv auf Variationen des Grundmagnetfeldes B0, eine trueFISP-Sequenz hingegen stärker sensitiv auf Variationen des Grundmagnetfeldes B0. FLASH-Sequenzen und FISP-Sequenzen sind sensitiv auf lokale Variationen eines eingestrahlten Hochfrequenzfeldes B1+.
  • In dem in 3 dargestellten Beispiel sind auf der Achse 26 die jeweiligen Nummern von aufgenommenen Bilddatensätzen in der Zeitserie und auf der Achse 27 unterschiedliche Größen aufgetragen. Als erste Größe ist der Flipwinkel in ° von 0° beim Ursprung bis 90° beim Achsenpunkt 28 aufgetragen. Die Achse 26 läuft in dem dargestellten Beispiel vom Bilddatensatz 1 bis zum Bilddatensatz 3000.
  • Die 3000 Bilddatensätze sind auf zwölf Abschnitte 29, 30, 31, 32, 33, 34, 35, 36, 37, 38 und 39 verteilt.
  • Im ersten Abschnitt 29 ist über die Kurve 40 für zweihundert Bilddatensätze der Flipwinkel aufgetragen, der bei der Aufnahme verwendet wurde, wobei für die Aufnahme in Abschnitt 29 eine FISP-Sequenz eingesetzt werden kann. Wie zu 1 beschrieben wird nach dem Anlegen eines RF-Anregungspulses mit einem bestimmten Flipwinkel ein kompletter Bilddatensatz aufgenommen und dann der nächste RF-Anregungspuls mit dem nächsten Flipwinkel angelegt und ein weiterer Bilddatensatz aufgenommen. 3 zeigt in Abschnitt 29 eine Flipwinkelverteilung, die einer sin2-Halbkurve entspricht. Der maximale Flipwinkel kann z.B. 24° betragen und es können konstante Phasen verwendet werden.
  • Für den hundertsten Bilddatensatz ist rein exemplarisch eine Linie 41 eingetragen. Der korrespondierende Flipwinkel ist der maximale Flipwinkel der Kurve 40.
  • Im zweiten Abschnitt 30 werden im dargestellten Beispiel vierhundert Bilddatensätze mit einem anderen Sequenztypen, beispielsweise mit einer TrueFISP-Sequenz akquiriert. Dabei werden Flipwinkel gemäß der Kurven 42 und 43 eingesetzt. Bei der Kurve 42 reichen diese bis 45° und bei der Kurve 43 bis 72°.
  • Auch für den Abschnitt 30 ist rein exemplarisch beim Flipwinkel für den vierhundertsten Bilddatensatz eine Linie 44 eingezeichnet. Hier beträgt der Flipwinkel 1°.
  • Eine Besonderheit stellt im Abschnitt 30 die Verwendung zweier unterschiedlicher Phasenzyklen dar. Beim Durchlaufen der Flipwinkel der Kurve 42 wird ein 00-Phasenzyklus bzw. kein Phasenzyklus verwendet und beim Durchlaufen der Kurve 43 ein 180°-Phasenzyklus. Ein 00-Phasenzyklus bezeichnet eine feststehende Phase.
  • Im folgenden Abschnitt 31 sind in der Kurve 45 die Flipwinkel zur Aufnahme von vierhundertfünfzig Bilddatensätzen mit einem weiteren Sequenztyp, beispielsweise einer FLASH-Sequenz angegeben. Diese sind kleiner als in der FISP- oder TrueFISP-Sequenz und laufen bis 6°. Auch ihre Verteilung ist eine sin2-Verteilung.
  • Zusätzlich zur Variation der Flipwinkel wird beim wiederholten Durchlaufen der FLASH-Sequenz ein Phasenzyklus zur Realisierung eines RF-Spoiling angelegt. Dabei wird wie beschrieben die Phase um Vielfache von 117° erhöht.
  • Die Abfolge der verschiedenen verwendeten Sequenztypen in den Abschnitten 29, 30 und 31 bilden zusammen einen Block 45. Dieser wird in 3 insgesamt dreimal verwendet. Dabei wird allein auf den Typ der Sequenz aber nicht auf die Anzahl der Bilddatensätze oder die Flipwinkelkurve abgestellt.
  • Im Abschnitt 32 werden 200 Bilddatensätze wieder mit dem ersten Sequenztyp des Blocks 45, also beispielsweise mit einer FISP-Sequenz, aufgenommen. Die Phase ist wie in Abschnitt 29 konstant, aber der maximale Flipwinkel beträgt 45°. Die verwendeten Flipwinkel liegen auf der Kurve 46.
  • In Abschnitt 33 folgen 200 Bilddatensätze, die mit dem zweiten Sequenztyp des Blocks 45, also beispielsweise mit einer TrueFISP-Sequenz, zu akquirieren sind. Hier wird ein 90°-Phasenzyklus eingesetzt, der maximale Flipwinkel liegt bei 50°. Die Flipwinkel sind auf der Kurve 47 aufgetragen.
  • Die nächsten ca. 450 Bilddatensätze in Abschnitt 34 sind, wie in Abschnitt 31, mit dem dritten Sequenztyp des Blocks 45, also beispielsweise mit einer FLASH-Sequenz, aufzunehmen. Die Kurve 48 zeigt eine sin2-Verteilung mit einem Maximalwert von 14°.
  • Kurve 49 in Abschnitt 35 läuft bis 72° und zeigt die Flipwinkel des Hochfrequenzimpulses 19 beim drittmaligen Verwenden des ersten Sequenztyps des Blocks 45, also beispielsweise mit einer FISP-Sequenz. Auch in diesem Durchlauf ist die Phase konstant.
  • Beim Akquirieren weiterer zweihundert Bilddatensätze mit dem zweiten Sequenztyp des Blocks 45, also beispielsweise mit einer TrueFISP-Sequenz, wird ein 270°-Phasenzyklus eingesetzt. Die Flipwinkel, die in der Kurve 50 in Abschnitt 36 aufgetragen sind, laufen bis 65°.
  • Die nächsten ca. 450 Bilddatensätze in Abschnitt 37 werden mit dem dritten Sequenztyp des Blocks 45, also beispielsweise mit der FLASH-Sequenz, aufgenommen. Die Kurve 51 stellt einen Flipwinkelverlauf bis maximal 20° dar, wiederum sin2-verteilt.
  • Im letzten Abschnitt 38 liegen zwei Kurven 52 und 53 zur Aufnahme von Bilddatensätzen mit dem ersten Sequenztyp des Blocks 45, also beispielsweise mit einer FISP-Sequenz. Diese stellen wiederum Flipwinkelverläufe dar. Wie bereits in den vorhergehenden Abschnitten wird beispielsweise bei der Messdatenaufnahme mit einer FISP-Sequenz eine konstante Phase verwendet.
  • Das dargestellte Beispiel erlaubt durch seine Ausgestaltung mit Abschnitten verschiedener Sequenztypen mit unterschiedlichen Sensitivitäten eine Bestimmung von mehreren gewebespezifischen und messungsspezifischen Parametern in einem Matching-Schritt, insbesondere können die Parameter T1, T2 sowie B0 und B1+ bestimmt werden. Das dargestellte Beispiel ist jedoch nicht einschränkend zu lesen. Das Erfassen einer Bildpunkt-Zeit-Serie kann grundsätzlich auch mit Hilfe von nur einem Sequenztyp erfolgen. Die Verwendung mehrerer verschiedener Sequenztypen mit unterschiedlichen Eigenschaften, insbesondere mit Bezug auf die jeweiligen Sensitivitäten auf gewebespezifische und/oder messungsspezifische Parameter, erhöht jedoch die Anzahl der aus einer derartig aufgenommenen Bildpunkt-Zeit-Serie bestimmbaren Parameter und/oder die Qualität der bestimmten Werte der Parameter entsprechend der Eigenschaften der verwendeten Sequenztypen.
  • Auch die dargestellten Anzahlen an mit einem Sequenztyp aufgenommenen Bilddatensätzen sowie die gezeigten Flipwinkelverläufe sind nur als Beispiel zu sehen.
  • In 4 sind zur Veranschaulichung eines Ausgleichsprozesses als Beispiel räumliche Verteilungen von bestimmten Grundmagnetfeldwerten B0 in einem Untersuchungsobjekt dargestellt. Während der Aufnahme der Messdaten aus denen die Bilddaten rekonstruiert wurden, deren Signalverläufe pro Bildpunkt als Bildpunkt-Zeit-Serien zur Bestimmung zumindest der lokalen Werte des Grundmagnetfeldes B0 mit Vergleichssignalverläufen verglichen werden sollen, wurde zur Verdeutlichung des Effektes ein linearer Verlauf des Grundfeldmagnetes B0 „von oben nach unten“ in dem Untersuchungsobjekt erzeugt.
  • Die links dargestellte räumlich aufgelöste Verteilung von Grundmagnetfeldwerten (B0-Karte) in einem Untersuchungsobjekt entspricht den mit einem Signalvergleich von Bildpunkt-Zeit-Serien mit Vergleichssignalverläufen bildpunktweise bestimmten Grundmagnetfeldwerten, wobei die für den Signalvergleich verwendeten Vergleichssignalverläufe, Vergleichssignalverläufe eines Dictionarys sind, dessen Vergleichssignalverläufe einen kleineren B0-Wertebereich abdecken als BO-Werte im Untersuchungsobjekt auftreten. Die links dargestellte BO-Karte deckt somit z.B. BO-Werte in einem durch die Aufnahmeart der den Bildpunkt-Zeit-Serien zugrundeliegenden Messdaten und die entsprechenden Vergleichssignalverläufe (den Dictionary) vorgegebenen Bereich von -40 bis +40 in beliebigen Einheiten ab.
  • Die in der Mitte dargestellte Differenz-B0-Karte entspricht einem Differenzbild der links dargestellten BO-Karte mit einer räumlich grob aufgelösten BO-Karte, die aus bestimmten, B0-sensitiven Abschnitten der Bildpunkt-Zeit-Serie abgeleiteten B0-Werten besteht, und somit BO-Werte in infinitesimal kleinen Schritten abdeckt. Die grob aufgelöste BO-Karte kann beispielsweise aus Abschnitten der Bildpunkt-Zeit-Serie, in denen die Messdaten mit einem Sequenztyp mit Eigenschaften einer FISP-Sequenz aufgenommen wurden, durch Bestimmung der Phasendifferenzen der Signale in den genannten Abschnitten der Bildpunkt-Zeit-Serie, die zu unterschiedlichen Echozeiten aufgenommen wurden, ggf. unter Verwendung einer Glättungsoperation, bestimmt werden.
  • Rundet man die Differenz-B0-Karte der genannten räumlich grob aufgelösten BO-Karte und der links dargestellten durch Signalvergleich erstellten BO-Karte auf Vielfache von 1/TR, ergeben sich Segmente, die Vielfachen von 1/TR entsprechen (in dem gezeigten Beispiel kann man fünf derartige Segmente erkennen). Zieht man das Ergebnis dieser Rundung von der artefaktbehafteten BO-Karte (links) ab, erhält man eine ausgeglichene, entfaltete BO-Karte wie sie rechts dargestellt ist. Die rechts dargestellte BO-Karte gibt den während der Messung künstlich erzeugten linearen Verlauf des Grundmagnetfeldes B0 sogar durch die beschriebene Verarbeitung mit der groben B0-Karte in einem größeren Wertebereich, in dem dargestellten Beispiel in einem Bereich von -200 bis +200 in denselben Einheiten wie in der links dargestellten BO-Karte, gut wieder.
  • Es wird hierbei ausgenutzt, dass aus den Signalen der verwendeten Bildpunkt-Zeit-Serien bereits ohne einen Signalvergleich Ausgleichsdaten, z.B. durch zumindest grob aufgelöste Bestimmung lokaler Werte eines zu bestimmenden Parameters, insbesondere eines messungsspezifischen Parameters, wie dem Grundmagnetfeld B0, erhalten werden können, welche Ausgleichsdaten zum Ausgleichen, z.B. entfalten von mittels Signalvergleich bestimmten Parameterwerten genutzt werden können. Obwohl die Auflösung einer solchen grob aufgelösten Parameterkarte nicht direkt als Ergebnis einer lokalen Verteilung der Parameterwerte ausreichend ist, können die Parameterwerte der grob aufgelösten Parameterkarte als Ausgleichsdaten genutzt werden, um höher aufgelöste Parameterkarten, die mittels Signalvergleich ermittelt wurden, auszugleichen, beispielsweise wie in dem oben beschriebenen Beispiel zu entfalten.
  • Somit können in der hohen Auflösung der erstellten Bildpunkt-Zeit-Serien ausgeglichene Werte der zu den Bildpunkt-Zeit-Serien bestimmten Parameter bestimmt werden, die zuvor bestimmte Werte des jeweiligen Parameters verbessern können und so Artefakte ausgleichen, wobei sogar der Bereich, der bei der Bestimmung der ausgeglichenen Parameterwerte erreichbaren Werte im Vergleich zu den bei einem reinen Signalvergleich erreichbaren Wertebereich vergrößert werden kann.
  • Mit einem MRF-Aufnahmeverfahren, das die Bestimmung der Parameter T1, T2, B0 und B1+ erlaubt, könnte ein beispielhafter möglicher Ablauf eines erfindungsgemäßen Verfahrens zur Bestimmung der lokalen Parameterwerte der Parameter T1, T2, B0 und B1 mit einer mehrstufigen Bestimmung der lokalen Parameterwerte wie folgt aussehen.
  • Es können nach einem ersten Signalvergleich der erfassten Bildpunkt-Zeit-Serien mit Vergleichssignalverläufen, bei dem insbesondere der gesamte Signalverlauf der Bildpunkt-Zeit-Serien mit dem jeweils gesamten Verlauf der Vergleichssignalverläufe verglichen wurde, lokale Werte (Karten) für alle Parameter T1, T2, B0 und B1+ bestimmt werden.
  • Eine so bestimmte BO-Karte kann mit einfaltungsartefaktähnlichen Artefakten behaftet sein und zunächst, wie oben mit Bezug auf 4 beschrieben, ausgeglichen, insbesondere entfaltet („unwrapped“) werden, sodass man eine ausgeglichene BO-Karte erhält.
  • Auf Grundlage der mittels des ersten Signalvergleichs bestimmten und ggf. ausgeglichenen BO-Karte kann zunächst, wie oben eine Korrektur, insbesondere ein Deblurring, der erfassten Bildpunkt-Zeit-Serien erfolgen, wodurch mögliche Blurring-Artefakte entfernt oder zumindest reduziert und entsprechend korrigierte Bildpunkt-Zeit-Serien erhalten werden.
  • In einem zweiten Signalvergleich können die so korrigierten Bildpunkt-Zeit-Serien mit den Vergleichssignalverläufen verglichen werden, um eine verbesserte Bestimmung der Parameter T1, T2, B0 und B1+ zu erreichen. Hierbei kann wieder der vollständige Signalverlauf der korrigierten Bildpunkt-Zeit-Serien mit dem gesamten Verlauf der Vergleichssignalverläufe verglichen werden.
  • Um beispielsweise die bestimmten Werte des Parameters B1+ weiter zu verbessern, können in einem nächsten Signalvergleich nur diejenigen Abschnitte der (korrigierten) Bildpunkt-Zeit-Serien mit entsprechenden Abschnitten der Vergleichssignalverläufe verglichen werden, die mit einer FLASH(-artigen) Sequenz aufgenommen wurden, wobei nur der Wert des Parameters B1+ bestimmt wird, die Werte der anderen Parameter T1, T2, B0 jedoch auf eine vorgegebene maximale Schwankung gegenüber den entsprechenden vorhergehend bestimmten Werten begrenzt werden, wobei insbesondere der Wert des Parameters B0 hierbei auch konstant gehalten werden kann. Auf diese Weise können möglicherweise die bestimmten Werte von B1+ verfälschende Flussartefakte eliminiert oder zumindest verringert werden, da die bei diesem Signalvergleich verwendeten Abschnitte der Signalverläufe nicht oder zumindest weniger sensitiv auf einen möglichen Fluss von Spins sind. Somit erhält man eine verbesserte B1+-Karte, die weiterhin mit einem Glättungsverfahren, das z.B. bildgeführte Filter umfasst, geglättet werden kann.
  • In einem weiteren Signalvergleich kann nochmals der gesamte Signalverlauf der korrigierten Bildpunkt-Zeit-Serien mit den gesamten Vergleichssignalverläufen verglichen werden, um die Werte der Parameter T1, T2 und B0 erneut und verbessert zu bestimmen, wobei die Werte der im vorhergehenden Signalvergleich bereits verbessert bestimmten Werte des Parameters B1+ konstant gehalten werden können.
  • Auf diese Weise erhält man schrittweise verbessert bestimmte Werte für die zu bestimmenden Parameter, wobei eine besonders hohe Genauigkeit erreicht und die bestimmten Werte verfälschende Artefakte eliminiert oder verringert werden können. Jeder dieser Schritte kann insbesondere auch iterativ durchgeführt werden, um eine höhere Genauigkeit zu erreichen.
  • 5 stellt schematisch eine erfindungsgemäße Magnetresonanzanlage 1 dar. Diese umfasst eine Magneteinheit 3 zur Erzeugung des Grundmagnetfeldes, eine Gradienteneinheit 5 zur Erzeugung der Gradientenfelder, eine Hochfrequenzeinheit 7 zur Einstrahlung und zum Empfang von Hochfrequenzsignalen und eine zur Durchführung eines erfindungsgemäßen Verfahrens ausgebildete Steuereinrichtung 9. In der 5 sind diese Teileinheiten der Magnetresonanzanlage 1 nur grob schematisch dargestellt. Insbesondere kann die Hochfrequenzeinheit 7 aus mehreren Untereinheiten, beispielsweise aus mehreren Spulen wie den schematisch gezeigten Spulen 7.1 und 7.2 oder mehr Spulen bestehen, die entweder nur zum Senden von Hochfrequenzsignalen oder nur zum Empfangen der ausgelösten Hochfrequenzsignale oder für beides ausgestaltet sein können.
  • Zur Untersuchung eines Untersuchungsobjektes U, beispielsweise eines Patienten oder auch eines Phantoms, kann dieses auf einer Liege L in die Magnetresonanzanlage 1 in deren Messvolumen eingebracht werden. Die Schicht S stellt ein exemplarisches Zielvolumen des Untersuchungsobjekts dar, aus dem Messdaten aufgenommen werden sollen.
  • Die Steuereinrichtung 9 dient der Steuerung der Magnetresonanzanlage und kann insbesondere die Gradienteneinheit 5 mittels einer Gradientensteuerung 5' und die Hochfrequenzeinheit 7 mittels einer Hochfrequenz-Sende-/Empfangs-Steuerung 7' steuern. Die Hochfrequenzeinheit 7 kann hierbei mehrere Kanäle umfassen, auf denen Signale gesendet oder empfangen werden können.
  • Die Hochfrequenzeinheit 7 ist zusammen mit ihrer Hochfrequenz-Sende-/Empfangs-Steuerung 7' für die Erzeugung und das Einstrahlen (Senden) eines Hochfrequenz-Wechselfeldes zur Manipulation der Spins in einem zu manipulierenden Bereich (beispielsweise in zu messenden Schichten S) des Untersuchungsobjekts U zuständig. Dabei wird die Mittenfrequenz des, auch als B1-Feld bezeichneten, Hochfrequenz-Wechselfeldes möglichst so eingestellt, dass sie nahe der Resonanzfrequenz der zu manipulierenden Spins liegt. Abweichungen von der Mittenfrequenz von der Resonanzfrequenz werden als Off-Resonanz bezeichnet. Zur Erzeugung des B1-Feldes werden in der Hochfrequenzeinheit 7 mittels der Hochfrequenz-sende/empfangs-Steuerung 7' gesteuerte Ströme an den HF-Spulen angelegt.
  • Weiterhin umfasst die Steuereinrichtung 9 eine Parameterwertbestimmungseinheit 15, mit welcher erfindungsgemäße Signalvergleiche zur Bestimmung von Parameterwerten durchgeführt werden können. Die Steuereinrichtung 9 ist insgesamt dazu ausgebildet, ein erfindungsgemäßes Verfahren durchzuführen.
  • Eine von der Steuereinrichtung 9 umfasste Recheneinheit 13 ist dazu ausgebildet alle für die nötigen Messungen und Bestimmungen nötigen Rechenoperationen auszuführen. Hierzu benötigte oder hierbei ermittelte Zwischenergebnisse und Ergebnisse können in einer Speichereinheit S der Steuereinrichtung 9 gespeichert werden. Die dargestellten Einheiten sind hierbei nicht unbedingt als physikalisch getrennte Einheiten zu verstehen, sondern stellen lediglich eine Untergliederung in Sinneinheiten dar, die aber auch z.B. in weniger oder auch in nur einer einzigen physikalischen Einheit realisiert sein können.
  • Über eine Ein-/Ausgabeeinrichtung E/A der Magnetresonanzanlage 1 können, z.B. durch einen Nutzer, Steuerbefehle an die Magnetresonanzanlage geleitet werden und/oder Ergebnisse der Steuereinrichtung 9 wie z.B. Bilddaten angezeigt werden.
  • Ein hierin beschriebenes Verfahren kann auch in Form eines Computerprogrammprodukts vorliegen, welches ein Programm umfasst und das beschriebene Verfahren auf einer Steuereinrichtung 9 implementiert, wenn es auf der Steuereinrichtung 9 ausgeführt wird. Ebenso kann ein elektronisch lesbarer Datenträger 26 mit darauf gespeicherten elektronisch lesbaren Steuerinformationen vorliegen, welche zumindest ein solches eben beschriebenes Computerprogrammprodukt umfassen und derart ausgestaltet sind, dass sie bei Verwendung des Datenträgers 26 in einer Steuereinrichtung 9 einer Magnetresonanzanlage 1 das beschriebene Verfahren durchführen.
  • ZITATE ENTHALTEN IN DER BESCHREIBUNG
  • Diese Liste der vom Anmelder aufgeführten Dokumente wurde automatisiert erzeugt und ist ausschließlich zur besseren Information des Lesers aufgenommen. Die Liste ist nicht Bestandteil der deutschen Patent- bzw. Gebrauchsmusteranmeldung. Das DPMA übernimmt keinerlei Haftung für etwaige Fehler oder Auslassungen.
  • Zitierte Nicht-Patentliteratur
    • Ma et al., „Magnetic Resonance Fingerprinting“, Nature, 495: S. 187-192 (2013) [0010]

Claims (14)

  1. Verfahren zur Bestimmung von Parameterwerten in Bildpunkten eines Untersuchungsobjektes mittels einer Magnetresonanz-Fingerprinting(MRF)-Technik, umfassend die Schritte: - Erfassen mindestens einer Bildpunkt-Zeit-Serie (BZS) mit Hilfe eines MRF-Aufnahmeverfahrens, aus denen je ein Wert (P') zumindest eines in dem jeweiligen Bildpunkt dargestellten Ort des Untersuchungsobjekts vorliegenden Parameters bestimmt werden soll, - Durchführen eines Signalvergleichs (103) mindestens eines Abschnitts des jeweiligen Signalverlaufs der erfassten Bildpunkt-Zeit-Serien (BZS) mit einem entsprechenden Abschnitt von Vergleichssignalverläufen (D) zur Bestimmung der jeweiligen Werte (P) des mindestens einen Parameters, - Durchführen eines weiteren Signalvergleichs (105) mindestens eines Abschnittes des jeweiligen Signalverlaufs unter Berücksichtigung der Ergebnisse (P) des vorhergehenden Signalvergleichs mit einem entsprechenden Abschnitt der Vergleichssignalverläufe (D) zur verbesserten Bestimmung der jeweiligen Werte (P') des mindestens einen Parameters, - Ausgeben der für den jeweiligen Bildpunkt verbessert bestimmten, jeweiligen Werte (P') der zu bestimmenden Parameter.
  2. Verfahren nach Anspruch 1 , wobei bei mindestens einem Signalvergleich der Signalverlauf in Abschnitten einer gesamten zu vergleichenden Bildpunkt-Zeit-Serie (BZS), die mit einem auf den zu bestimmenden Parameter sensitiven Sequenztyp erfasst wurden, mit den entsprechenden Abschnitten der Vergleichssignalverläufe (D) verglichen werden.
  3. Verfahren nach einem der vorhergehenden Ansprüche, wobei durch einen Signalvergleich mindestens eines Abschnitts des Signalverlaufs der erfassten Bildpunkt-Zeit-Serie (BZS) mit einem entsprechenden Abschnitt der Vergleichssignalverläufe (D) ein Wert (P, P') eines Parameters, insbesondere des Grundmagnetfeldes B0, bestimmt wird, und der bestimmte Wert des Parameters mit Hilfe von aus zumindest Teilen der erfassten Bildpunkt-Zeit-Serie erhaltenen Ausgleichsdaten ausgeglichen wird, sodass ein ausgeglichener Wert (P*) des Parameters bestimmt wird, der den zuvor bestimmten Wert (P, P') des Parameters ersetzt.
  4. Verfahren nach einem der vorhergehenden Ansprüche, wobei auf Grundlage eines für einen Bildpunkt mit einem erfindungsgemäßen Signalvergleich bestimmten Wertes (P, P') eines Parameters die Bildpunkt-Zeit-Serie (BZS) des Bildpunktes korrigiert wird, sodass eine korrigierte Bildpunkt-Zeit-Serie (BZS') erstellt wird.
  5. Verfahren nach Anspruch 4, wobei zur Berücksichtigung der Ergebnisse eines vorhergehenden Signalvergleichs bei einem weiteren Signalvergleich der weitere Signalvergleich mindestens einen Abschnitt des Signalverlaufs einer korrigierten Bildpunkt-Zeit-Serie (BZS') mit einem entsprechenden Abschnitt der Vergleichssignalverläufe (D) zur verbesserten Bestimmung der Werte (P') der zu bestimmenden Parameter vergleicht.
  6. Verfahren nach einem der Ansprüche 4 oder 5, wobei eine korrigierte Bildpunkt-Zeit-Serie (BZS') auf Grundlage eines durch Signalvergleich der korrigierten Bildpunkt-Zeit-Serie (BZS') mit den Vergleichssignalverläufen (D) verbessert bestimmten Wertes (P') des Parameters erneut korrigiert wird, und erneut der Wert (P') des Parameter durch Vergleichen der erneut korrigierten Bildpunkt-Zeit-Serie (BZS') mit den Vergleichssignalverläufen (D) verbessert bestimmt wird.
  7. Verfahren nach Anspruch 6, wobei eine erneute Korrektur solange (iterativ) durchgeführt wird, bis ein Unterschied von erneut verbessert bestimmten Werten (P') des Parameters zu den entsprechenden vor der erneuten Korrektur bestimmten Werten (P, P`) des Parameters kleiner als ein vorgegebener Schwellwert (S) ist.
  8. Verfahren nach einem der vorhergehenden Ansprüche, wobei Werte (P) von mindestens zwei verschiedenen Parametern bestimmt werden sollen, und zur Berücksichtigung der Ergebnisse des vorhergehenden Signalvergleichs bei Signalvergleichen zur verbesserten Bestimmung der Werte (P') mindestens eines zu bestimmenden Parameters (G1), der bzw. die Werte (P, P') mindestens eines anderen (G2) als des bzw. der aktuell verbessert zu bestimmenden Parameter (G1) auf eine vorgegebene maximale Schwankung (SW) gegenüber dem Ergebnis des vorhergehenden Signalvergleichs begrenzt, insbesondere konstant gehalten, wird.
  9. Verfahren nach Anspruch 8, wobei solange Signalvergleiche zur verbesserten Bestimmung eines Wertes (P') eines zu bestimmenden Parameters, durchgeführt werden, bis alle Werte der zu bestimmenden Parameter mindestens einmal bei der verbesserten Bestimmung bestimmt wurden.
  10. Verfahren nach einem der vorhergehenden Ansprüche, wobei bei der Erfassung der Bildpunkt-Zeit-Serie (BZS) mindestens ein Sequenztyp mit Eigenschaften, insbesondere Sensitivitäten, von Sequenztypen aus der Gruppe der Sequenztypen TrueFISP („True Fast Imaging with Steady-state free Precession“), FISP („Fast Imaging with Steady-state free Precession“) und FLASH („Fast Low-Angle Shot“), insbesondere alle genannten Sequenztypen, eingesetzt wird.
  11. Verfahren nach einem der vorhergehenden Ansprüche, wobei die Erfassung der Bildpunkt-Zeit-Serie (BZS) entlang einer nicht-kartesischen, insbesondere entlang einer spiralförmigen, Abtastung des k-Raums erfolgt.
  12. Magnetresonanzanlage (1) umfassend, eine Magneteinheit (3), eine Gradienteneinheit (5), eine Hochfrequenzeinheit (7) und eine Steuereinrichtung (9) mit einer Hochfrequenz-Sende-/Empfangs-Steuerung (7') und einer Parameterwertbestimmungseinheit (15), wobei die Steuereinrichtung (9) dazu ausgebildet ist, ein Verfahren nach einem der Ansprüche 1 bis 11 auf der Magnetresonanzanlage (1) auszuführen.
  13. Computerprogramm, welches direkt in einen Speicher einer Steuereinrichtung (9) einer Magnetresonanzanlage (1) ladbar ist, mit Programm-Mitteln, um die Schritte des Verfahrens nach einem der Ansprüche 1 bis 11 auszuführen, wenn das Programm in der Steuereinrichtung (9) der Magnetresonanzanlage (1) ausgeführt wird.
  14. Elektronisch lesbarer Datenträger mit darauf gespeicherten elektronisch lesbaren Steuerinformationen, welche zumindest ein Computerprogramm nach Anspruch 13 umfassen und derart ausgestaltet sind, dass sie bei Verwendung des Datenträgers in einer Steuereinrichtung (9) einer Magnetresonanzanlage (1) ein Verfahren nach einem der Ansprüche 1 bis 11 durchführen.
DE102018209584.1A 2018-06-14 2018-06-14 Magnetresonanz-Fingerprinting-Verfahren Pending DE102018209584A1 (de)

Priority Applications (3)

Application Number Priority Date Filing Date Title
DE102018209584.1A DE102018209584A1 (de) 2018-06-14 2018-06-14 Magnetresonanz-Fingerprinting-Verfahren
US16/441,383 US10921406B2 (en) 2018-06-14 2019-06-14 Magnetic resonance fingerprinting method and apparatus
CN201910515214.5A CN110604569B (zh) 2018-06-14 2019-06-14 确定图像元中的参数值的方法、磁共振设备和数据载体

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
DE102018209584.1A DE102018209584A1 (de) 2018-06-14 2018-06-14 Magnetresonanz-Fingerprinting-Verfahren

Publications (1)

Publication Number Publication Date
DE102018209584A1 true DE102018209584A1 (de) 2019-12-19

Family

ID=68724680

Family Applications (1)

Application Number Title Priority Date Filing Date
DE102018209584.1A Pending DE102018209584A1 (de) 2018-06-14 2018-06-14 Magnetresonanz-Fingerprinting-Verfahren

Country Status (3)

Country Link
US (1) US10921406B2 (de)
CN (1) CN110604569B (de)
DE (1) DE102018209584A1 (de)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110623666A (zh) * 2018-06-21 2019-12-31 西门子医疗有限公司 磁共振指纹方法

Families Citing this family (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP3629048A1 (de) * 2018-09-27 2020-04-01 Siemens Healthcare GmbH Niederfeld-magnetresonanz-fingerprinting
JP7270520B2 (ja) * 2019-10-09 2023-05-10 キヤノンメディカルシステムズ株式会社 磁気共鳴イメージング装置および画像処理装置
EP3828581B8 (de) * 2019-11-28 2024-02-21 Siemens Healthineers AG Verfahren zur erzeugung wenigstens eines bilddatensatzes und eines referenzbilddatensatzes, computerprogrammprodukt sowie magnetresonanzanlage

Family Cites Families (36)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP2233941A1 (de) * 2009-03-25 2010-09-29 Koninklijke Philips Electronics N.V. MR-Bildgebung mit parallelem Signalerhalt
US9008396B2 (en) * 2011-04-29 2015-04-14 Magnetic Resonance Innovations, Inc. Tissue similarity mapping
US9097781B2 (en) * 2012-04-12 2015-08-04 Mark Griswold Nuclear magnetic resonance (NMR) fingerprinting with parallel transmission
DE102012206585B4 (de) * 2012-04-20 2013-12-12 Siemens Aktiengesellschaft Verfahren zur schnellen ortsaufgelösten Bestimmung eines Magnetresonanz-Relaxationsparameters in einem Untersuchungsgebiet
GB201304728D0 (en) * 2013-03-15 2013-05-01 Isis Innovation Medical imaging
US10209335B2 (en) * 2014-04-21 2019-02-19 Case Western Reserve University Nuclear magnetic resonance (NMR) fingerprinting with singular value decomposition (SVD) compression
US10261154B2 (en) * 2014-04-21 2019-04-16 Case Western Reserve University Nuclear magnetic resonance (NMR) fingerprinting tissue classification and image segmentation
US10281547B2 (en) * 2014-04-21 2019-05-07 Case Western Reserve University Diffusion-weighted double-echo magnetic resonance fingerprinting (MRF)
DE102014224656A1 (de) * 2014-12-02 2016-06-02 Siemens Aktiengesellschaft Verfahren und Vorrichtung zum Segmentieren eines medizinischen Untersuchungsobjekts mit quantitativen MR-Bildgebungsmethoden
CN106157287B (zh) * 2015-04-03 2019-09-06 东芝医疗系统株式会社 数据处理装置、数据处理方法和医学图像设备
US10209334B2 (en) * 2015-04-08 2019-02-19 The General Hospital Corporation System and method for controlling partial volume effects in magnetic resonance fingerprinting
US10761171B2 (en) * 2015-06-22 2020-09-01 Case Western Reserve University Systems and methods for free-breathing three-dimensional magnetic resonance fingerprinting
DE102015215778B4 (de) * 2015-08-19 2023-08-03 Siemens Healthcare Gmbh Verfahren zu einer Darstellung von Magnetresonanz-Daten
US10175321B2 (en) * 2015-10-23 2019-01-08 Siemens Healthcare Gmbh Iterative magnetic resonance fingerprinting reconstruction
CN108351395B (zh) * 2015-10-27 2021-02-26 皇家飞利浦有限公司 来自磁共振图像的虚拟ct图像
US9633455B1 (en) * 2015-12-10 2017-04-25 Siemens Healthcare Gmbh Dictionary-free MR parameter inference for fingerprinting reconstruction
JP6687383B2 (ja) * 2015-12-21 2020-04-22 キヤノンメディカルシステムズ株式会社 磁気共鳴イメージング装置
US10671939B2 (en) * 2016-04-22 2020-06-02 New York University System, method and computer-accessible medium for learning an optimized variational network for medical image reconstruction
US10634748B2 (en) * 2016-06-22 2020-04-28 Comsats Institute Of Information Technology Graphical processing unit (GPU) implementation of magnetic resonance fingerprinting (MRF)
WO2018006058A1 (en) * 2016-07-01 2018-01-04 Cubisme, Inc. System and method for forming a super-resolution biomarker map image
DE102016214775A1 (de) * 2016-08-09 2018-02-15 Siemens Healthcare Gmbh Bestimmung einer Eigenschaft eines Organs
DE102016217675A1 (de) * 2016-09-15 2018-03-15 Siemens Healthcare Gmbh Verbesserte Erzeugung von Bildpunkt-Zeit-Serien eines Untersuchungsobjektes mittels Magnetresonanztechnik
WO2018106713A1 (en) * 2016-12-06 2018-06-14 Darmiyan, Inc. Methods and systems for identifying brain disorders
WO2018108643A1 (en) * 2016-12-15 2018-06-21 Koninklijke Philips N.V. Multi-state magnetic resonance fingerprinting
US10705168B2 (en) * 2017-01-17 2020-07-07 Case Western Reserve University System and method for low rank approximation of high resolution MRF through dictionary fitting
US20180231626A1 (en) * 2017-02-10 2018-08-16 Case Western Reserve University Systems and methods for magnetic resonance fingerprinting for quantitative breast imaging
US10552955B2 (en) * 2017-04-04 2020-02-04 Yale University Imaging acceleration methods for MRI parameter mapping
US10459055B2 (en) * 2017-04-07 2019-10-29 Case Western Reserve University System and method for reduced field of view MR fingerprinting for parametric mapping
WO2018187760A1 (en) * 2017-04-07 2018-10-11 Case Western Reserve University System and method for dynamic multiple contrast enhanced, magnetic resonance fingerprinting (dmce-mrf)
US11232853B2 (en) * 2017-04-21 2022-01-25 Cubisme, Inc. System and method for creating, querying, and displaying a MIBA master file
US11335464B2 (en) * 2018-01-12 2022-05-17 Siemens Medical Solutions Usa, Inc. Integrated precision medicine by combining quantitative imaging techniques with quantitative genomics for improved decision making
CN110346741A (zh) * 2018-04-06 2019-10-18 凯斯西储大学 使用多个脉冲序列类型的磁共振指纹的系统和方法
EP3572825A3 (de) * 2018-05-24 2019-12-11 Siemens Healthcare GmbH Wörterbuchkompression für magnetresonanzfingerprinting
EP3579010B1 (de) * 2018-06-07 2021-07-28 Siemens Healthcare GmbH Verbesserte bestimmung von quantitativen gewebeparametern
US11062488B2 (en) * 2018-06-28 2021-07-13 Siemens Healthcare Gmbh Magnetic resonance fingerprinting based on bloch manifold and spatial regularizations in parallel
DE102018221695A1 (de) * 2018-12-13 2020-06-18 Siemens Healthcare Gmbh Verfahren zur quantitativen Magnetresonanzbildgebung, Magnetresonanzeinrichtung, Computerprogramm und elektronisch lesbarer Datenträger

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
Ma et al., „Magnetic Resonance Fingerprinting", Nature, 495: S. 187-192 (2013)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110623666A (zh) * 2018-06-21 2019-12-31 西门子医疗有限公司 磁共振指纹方法
CN110623666B (zh) * 2018-06-21 2023-07-21 西门子医疗有限公司 磁共振指纹方法

Also Published As

Publication number Publication date
CN110604569A (zh) 2019-12-24
US20190383892A1 (en) 2019-12-19
US10921406B2 (en) 2021-02-16
CN110604569B (zh) 2023-07-18

Similar Documents

Publication Publication Date Title
EP3296757B1 (de) Artefaktarme magnetresonanz-fingerprinting-messung
EP3078978B1 (de) Verfahren zur magnetresonanz-bildgebung
DE102011006851B4 (de) Gruppierung von diffusionsgewichteten MR-Bildern nach den verwendeten Diffusionsgradienten zur Erstellung einer Diffusionsinformation
DE102016218955B4 (de) Optimiertes Magnetresonanz-Schichtmultiplexing-Verfahren
DE102011083406B4 (de) Verfahren zur Auswahl eines Unterabtastungsschemas für eine MR-Bildgebung, Verfahren zur Magnetresonanz-Bildgebung und Magnetresonanzanlage
EP3413072A1 (de) Mrt-schicht-multiplexing-verfahren
DE102014206395B4 (de) Aufnahme und iterative Rekonstruktion einer Parameterkarte eines Zielbereichs
DE102018209584A1 (de) Magnetresonanz-Fingerprinting-Verfahren
DE102014200006A1 (de) Rekonstruktion von Magnetresonanz-Rohdaten
DE102014203068B4 (de) Pseudo-zufällige Erfassung von MR-Daten eines zweidimensionalen Volumenabschnitts
DE102019206827A1 (de) Magnetresonanz-Fingerprinting-Verfahren
DE102016207314A1 (de) Verfahren zur Anzeige von quantitativen Magnetresonanz-Bilddaten
DE102014214844A1 (de) Verfahren und Magnetresonanzanlage zum Erfassen von MR-Daten
EP3457159A1 (de) Verbessertes magnetresonanz-fingerprinting-verfahren
DE102014222964A1 (de) Dynamische Verzeichnungskorrektur von Magnetresonanz-Bildaufnahmen
DE102020212250B4 (de) Verfahren zur Ermittlung einer Point-Spread-Funktion (PSF) für eine Rekonstruktion von Bilddaten aus mittels einer Mag-netresonanzanlage aufgenommenen Messdaten
DE102014222495A1 (de) Korrektur von Magnetresonanz-Bildaufnahmen mit Hilfe konvertierter Feldkarten
DE102020209382A1 (de) Verfahren zur Aufnahme von Messdaten mittels einer Magnetresonanzanlage mit einer Korrektur der verwendeten k-Raumtrajektorien
DE102016200889A1 (de) Rekonstruktion von Bilddaten
EP3550318A1 (de) Verfahren zur erzeugung einer b0-karte mittels eines magnetresonanz-fingerabdruckverfahrens
EP3828581B1 (de) Verfahren zur erzeugung wenigstens eines bilddatensatzes und eines referenzbilddatensatzes, computerprogrammprodukt sowie magnetresonanzanlage
DE102021203269A1 (de) Verfahren zur Reduzierung von Artefakten in mittels Magnetresonanz aufgenommenen Bilddatensätzen
EP3629047A1 (de) Magnetresonanz-fingerprinting-verfahren basierend auf ähnlichkeitswerten
EP3629048A1 (de) Niederfeld-magnetresonanz-fingerprinting
DE102013223191A1 (de) Erfassen von MR-Daten einer Schicht eines Volumenabschnitts mit einer variablen Anzahl von Phasenkodierschritten

Legal Events

Date Code Title Description
R012 Request for examination validly filed
R016 Response to examination communication
R081 Change of applicant/patentee

Owner name: SIEMENS HEALTHINEERS AG, DE

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