DE102009050662B4 - Radiale MR-Datenerfassung in einem dreidimensionalen K-Raum mit einer Anordnung der Speichen gemäß einer Spiral-Phyllotaxis - Google Patents

Radiale MR-Datenerfassung in einem dreidimensionalen K-Raum mit einer Anordnung der Speichen gemäß einer Spiral-Phyllotaxis Download PDF

Info

Publication number
DE102009050662B4
DE102009050662B4 DE102009050662A DE102009050662A DE102009050662B4 DE 102009050662 B4 DE102009050662 B4 DE 102009050662B4 DE 102009050662 A DE102009050662 A DE 102009050662A DE 102009050662 A DE102009050662 A DE 102009050662A DE 102009050662 B4 DE102009050662 B4 DE 102009050662B4
Authority
DE
Germany
Prior art keywords
spokes
angle
points
spoke
space
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.)
Active
Application number
DE102009050662A
Other languages
English (en)
Other versions
DE102009050662A1 (de
Inventor
Davide Piccini
Michael Zenge
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 AG
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 AG filed Critical Siemens AG
Priority to DE102009050662A priority Critical patent/DE102009050662B4/de
Priority to US12/912,096 priority patent/US8648594B2/en
Publication of DE102009050662A1 publication Critical patent/DE102009050662A1/de
Application granted granted Critical
Publication of DE102009050662B4 publication Critical patent/DE102009050662B4/de
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01RMEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
    • G01R33/00Arrangements or instruments for measuring magnetic variables
    • G01R33/20Arrangements or instruments for measuring magnetic variables involving magnetic resonance
    • G01R33/44Arrangements or instruments for measuring magnetic variables involving magnetic resonance using nuclear magnetic resonance [NMR]
    • G01R33/48NMR imaging systems
    • G01R33/54Signal processing systems, e.g. using pulse sequences ; Generation or control of pulse sequences; Operator console
    • G01R33/56Image enhancement or correction, e.g. subtraction or averaging techniques, e.g. improvement of signal-to-noise ratio and resolution
    • G01R33/561Image 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
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01RMEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
    • G01R33/00Arrangements or instruments for measuring magnetic variables
    • G01R33/20Arrangements or instruments for measuring magnetic variables involving magnetic resonance
    • G01R33/44Arrangements or instruments for measuring magnetic variables involving magnetic resonance using nuclear magnetic resonance [NMR]
    • G01R33/48NMR imaging systems
    • G01R33/4818MR 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/4824MR 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
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01RMEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
    • G01R33/00Arrangements or instruments for measuring magnetic variables
    • G01R33/20Arrangements or instruments for measuring magnetic variables involving magnetic resonance
    • G01R33/44Arrangements or instruments for measuring magnetic variables involving magnetic resonance using nuclear magnetic resonance [NMR]
    • G01R33/48NMR imaging systems
    • G01R33/4818MR 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/4824MR 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
    • G01R33/4826MR 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 in three dimensions

Landscapes

  • Physics & Mathematics (AREA)
  • High Energy & Nuclear Physics (AREA)
  • Condensed Matter Physics & Semiconductors (AREA)
  • General Physics & Mathematics (AREA)
  • Health & Medical Sciences (AREA)
  • General Health & Medical Sciences (AREA)
  • Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
  • Radiology & Medical Imaging (AREA)
  • Engineering & Computer Science (AREA)
  • Signal Processing (AREA)
  • Magnetic Resonance Imaging Apparatus (AREA)

Abstract

Verfahren zur radialen Datenerfassung in einem dreidimensionalen K-Raum bei einer MR-Messung mittels einer Magnetresonanzanlage (5), wobei Daten im K-Raum entlang von geradlinigen Speichen erfasst werden, wobei jede der Speichen durch einen Punkt (1) auf einer Kugeloberfläche (4) und den Mittelpunkt dieser Kugel definiert ist, und wobei dieser Mittelpunkt dem Zentrum des zu erfassenden K-Raums entspricht, dadurch gekennzeichnet, dass die Punkte (1) auf der Kugeloberfläche (4) derart angeordnet werden, dass die Verteilung der Punkte (1) der Spiral-Phyllotaxis gehorcht.

Description

  • Die vorliegende Erfindung betrifft ein Verfahren und eine Vorrichtung zur radialen Datenerfassung in einem dreidimensionalen K-Raum bei einer MR-Messung für eine Magnetresonanzanlage. Darüber hinaus betrifft die vorliegende Erfindung eine Magnetresonanzanlage mit der vorab beschriebenen Vorrichtung.
  • Die US 6,281,681 B1 offenbart eine zweidimensionale Anordnung von K-Raum-Punkten längs einer Fibonacci-Spirale.
  • Die US 7,245,124 B1 beschreibt das Erfassen von MR-Daten mit einer Reihe von sphärisch verlaufenden Trajektorien.
  • Die US 4,651,096 beschreibt ein Verfahren zur Erstellung von MR-Bildern, wobei zeitabhängige Magnetfeldgradienten eingesetzt werden. Ein Integral dieser Gradienten über der Zeit weist dabei einen spiralförmigen Verlauf auf.
  • In den letzten Jahren wurde die radiale dreidimensionale Datenerfassung im K-Raum auf dem Gebiet der Magnetresonanztomographie aus mehreren Gründen populär. Erstens ist die Datenerfassung sehr robust gegenüber einer Bewegung (z. B. der Bewegung eines Patienten bei einer MR-Untersuchung), da bei den gängigen radialen dreidimensionalen Verfahren jeder Auslesevorgang (Vorgang, bei welchem in einem Schritt mehrere Punkte des K-Raums erfasst werden) durch das Zentrum des K-Raums verläuft. Darüber hinaus ermöglicht die radiale dreidimensionale Datenerfassung sehr einfach eine radiale Unterabtastung, wodurch die Abtastdichte (erfasste Punkte pro Volumeneinheit) an den Rändern des K-Raums relativ schwach ist, während sie in der Umgebung des Zentrums und im Zentrum selbst relativ groß ist.
  • Die radiale dreidimensionale Datenerfassung weist allerdings die folgenden zwei Problemfelder auf:
    • 1. Es ist schwierig, eine möglichst gleichförmige Verteilung der im K-Raum erfassten Datenpunkte zu erzielen.
    • 2. Die meisten Verfahren zur radialen dreidimensionalen Datenerfassung weisen ein verschachteltes Vorgehen auf, wobei der K-Raum durch Erfassungsschritte mehrfach abgetastet wird. Dabei wird der K-Raum von jedem Erfassungsschritt grob abgetastet, wobei jeder Erfassungsschritt im Wesentlichen andere Punkte des K-Raums abtastet oder erfasst als die jeweils anderen Erfassungsschritte. Die radiale dreidimensionale Datenerfassung des K-Raums setzt sich dabei aus den Ergebnissen aller verschachtelten Erfassungsschritte zusammen. Das Problem bei dieser verschachtelten Vorgehensweise besteht darin, eine Reihenfolge, in welcher die Punkte des K-Raums für einen der Erfassungsschritte erfasst werden, derart zu wählen, dass dies mit einer möglichst gleichmäßigen Gradientenänderung des Magnetfelds einhergeht, um die von Magnetfeldänderungen hervorgerufenen Wirbelstromeffekte zu minimieren.
  • Gemäß einem weithin akzeptierten Ansatz zur radialen Datenerfassung in einem dreidimensionalen K-Raum setzt sich die Datenerfassung aus mehreren Auslesevorgängen zusammen, wobei pro Auslesevorgang Punkte entlang einer Speiche (d. h. einer geradlinigen Strecke) erfasst werden, wobei diese Speiche durch einen Punkt auf einer Kugel und den Mittelpunkt dieser Kugel definiert ist. Mit anderen Worten verläuft jede Speiche, auf welcher die von dem entsprechenden Auslesevorgang erfassten Punkte in dem K-Raum erfasst werden, durch diesen Mittelpunkt, welcher sich im Zentrum des K-Raums befindet, und durch den entsprechenden Punkt auf der Kugel. Eine Speiche unterscheidet sich somit durch den entsprechenden Punkt auf der Kugel von den anderen Speichen, da jede Speiche durch den Mittelpunkt verläuft.
  • Die Punkte auf der Kugel, welche jeweils eine der Speichen definieren, liegen auf einer Trajektorie, welche die Form einer dreidimensionalen archimedischen Spirale aufweist. Je mehr Punkte im K-Raum abgetastet werden, umso mehr Speichen existieren und umso mehr Windungen weist die archimedische Spirale auf der Kugel auf, wodurch sich der Abstand von benachbarten Windungen der archimedischen Spirale verringert. Bei einem verschachtelten Vorgehen wird für jeden Erfassungsschritt nur jede m-te Speiche abgetastet, wenn m der Anzahl der Erfassungsschritte entspricht. Mit anderen Worten werden beispielsweise für den ersten Erfassungsschritt die Punkte auf den Speichen 1, (m + 1), (2·m + 1), (3·m + 1), usw. erfasst, während für den k-ten Erfassungsschritt die Punkte auf den Speichen k, (m + k), (2·m + k), (3·m + k), usw. erfasst werden (k ≤ m).
  • Bei diesem Vorgehen führt eine starke Verschachtelung, d. h. m >> 1, zu einer großen Gradientenänderung des Magnetfelds beim Übergang von einer Speiche zur nächsten, was nachteiliger Weise zu Artefakten, welche durch die Auswirkungen des Wirbelstroms auftreten.
  • In 1 sind für eine radiale Datenerfassung in einem dreidimensionalen K-Raum, wobei die Daten mittels 1600 Speichen erfasst werden, die diese Speichen definierenden Punkte 1 auf einer Kugel 4 dargestellt. Dabei zeigt die 1a die Kugel 4 von oben und die 1b zeigt die Kugel 4 schräg von oben. Die dargestellten Punkte 1 liegen auf einer dreidimensionalen archimedischen Spirale, welche die Trajektorie bildet, mittels welcher die Reihenfolge bestimmt wird, in welcher die einzelnen Speichen abgearbeitet werden.
  • Für den Fall, dass die gesamte Datenerfassung aus 20 Erfassungsschritten besteht, stellt die 1c die Trajektorie 15 für einen dieser Erfassungsschritte dar. Die geradlinigen Teilstücke repräsentieren dabei den Übergang von einer Speiche zur nächsten. Da diese geradlinigen Teilstücke zumindest teilweise relativ groß sind, führt dies zu einer relativ großen Veränderung bei den zur Datenerfassung notwendigen Gradienten des Magnetfelds, was wiederum starke Wirbelströme induziert, welche schließlich zu Artefakten bei der Bildgebung in Abhängigkeit der Datenerfassung führt.
  • Der vorab beschriebene Ansatz zur radialen Datenerfassung in einem dreidimensionalen K-Raum weist allerdings den Vorteil einer äußerst gleichmäßigen Abtastung auf, wodurch eine Kompensation der Abtastdichte stark erleichtert wird. Dabei wird unter der Kompensation der Abtastdichte der Vorgang verstanden, bei welchem die bereits oben angemerkte unterschiedliche Abtastdichte (groß im Zentrum und gering am Rand des K-Raums) für die Bildgebung derart kompensiert wird, dass die Dichte der aus den erfassten Daten ermittelten Bildpunkte im gesamten abgetasteten Volumen möglichst gleichmäßig ist.
  • In „Temporal Stability of Adaptive 3D Radial MRI Using Multidimensional Golden Means”, R. W. Chan und andere, Magnetic Resonance in Medicine 61, Seiten 354–363, 2009, ist ein weiterer Ansatz zur radialen Datenerfassung in einem dreidimensionalen K-Raum beschrieben. Das vorgestellte Verfahren führt zumindest im Zeitbereich eine gleichmäßige Abtastung eines zweidimensionalen normalisierten Raums aus, wobei das Konzept des goldenen Schnitts zusammen mit einer modifizierten Fibonacci-Reihe eingesetzt wird, um die Punkte 1 auf einer Kugel 4 zu verteilen, wobei diese Punkte 1 jeweils eine Speiche zur radialen Datenerfassung definieren. Dieses Abtastmuster wird auf eine Kugeloberfläche übertragen, indem die Koordinaten der planaren Abtastpunkte mit dem Polarwinkel und dem Azimutwinkel innerhalb des dreidimensionalen Raumes gekoppelt werden. Ein Ergebnis ist in 2 dargestellt, wobei 2a eine Ansicht der Kugel 4 von oben und 2b eine Ansicht der Kugel 4 von schräg oben darstellt.
  • Dieser Ansatz zielt insbesondere auf eine möglichst gleichmäßige Abtastung im zeitlichen Verlauf ab. Dennoch kommt man mit einer Voronoi-Analyse zu dem Ergebnis, dass die Abtastverteilung innerhalb eines Zeitrahmens nicht besonders gleichmäßig ist. Folglich erfordert dieser Ansatz eine komplizierte Kompensation der Abtastdichte. Schließlich werden durch diesen Ansatz Wirbelströme nicht optimal vermieden.
  • Zusammenfassend ist festzuhalten, dass die Ansätze zur radialen Datenerfassung in einem dreidimensionalen K-Raum nach dem Stand der Technik entweder hinsichtlich einer gleichmäßigen Abtastdichte – räumlich oder zeitlich – oder hinsichtlich einer geeigneten Kompensation von Wirbelströmen (hinsichtlich einer Vermeidung von zu großen Wirbelströmen) optimiert sind. Dabei ist in den meisten Fällen eine komplizierte Kompensation der Abtastdichte erforderlich.
  • Daher ist es die Aufgabe der vorliegenden Erfindung, eine radiale Datenerfassung in einem dreidimensionalen K-Raum bereitzustellen, bei welcher sowohl eine gleichmäßige räumliche Abtastdichte als auch eine Vermeidung von zu großen Wirbelströmen gewährleistet ist und bei welcher nur eine einfache Kompensation der Abtastdichte erforderlich ist.
  • Erfindungsgemäß wird diese Aufgabe durch ein Verfahren zur radialen Datenerfassung in einem dreidimensionalen K-Raum nach Anspruch 1, durch eine Vorrichtung für eine Magnetresonanzanlage zur radialen Datenerfassung in einem dreidimensionalen K-Raum nach Anspruch 13, durch eine Magnetresonanzanlage nach Anspruch 24, durch ein Computerprogrammprodukt nach Anspruch 25 und durch einen elektronisch lesbaren Datenträger nach Anspruch 26 gelöst. Die abhängigen Ansprüche definieren bevorzugte und vorteilhafte Ausführungsformen der vorliegenden Erfindung.
  • Im Rahmen der vorliegenden Erfindung wird ein Verfahren zur radialen Datenerfassung in einem dreidimensionalen K-Raum bei einer MR-Messung für eine Magnetresonanzanlage bereitgestellt. Dabei werden Daten in dem K-Raum entlang von geradlinigen Speichen erfasst. Jede dieser Speichen wird durch einen Punkt auf einer Kugel und den Mittelpunkt dieser Kugel definiert (d. h. der Punkt auf der Kugel und der Mittelpunkt liegen auf dieser Speiche), wobei der Mittelpunkt der Kugel dem Zentrum des K-Raums entspricht. Die Punkte werden dabei auf der Kugel bzw. Kugeloberfläche derart angeordnet, dass die Punkte oder die Verteilung der Punkte der Spiral-Phyllotaxis, insbesondere der Fibonacci-(Spiral-)Phyllotaxis, entsprechen.
  • Unter Phyllotaxis wird dabei die Lehre von der Blattstellung bzw. dem Blattstand bei Pflanzen verstanden. Unter der Blattstellung versteht man dabei die Anordnung von blattähnlichen Organen bei einer Pflanze. Blattähnliche Organe umfassen dabei Blätter und Fruchtblätter einer Pflanze, wie z. B. die Schuppen eines Nadelbaumzapfens, die Fruchtspelzen einer Ananas und die Samen in einer Sonnenblumenblüte. Die Anordnung der blattähnlichen Organe gemäß der Phyllotaxis ist dabei häufig auch spiralförmig, wobei der Aufbau dieser Spirale bestimmten Regeln gehorcht. Diese Regeln sind gemäß der Spiral-Phyllotaxis definiert, wozu auf „A Collision-based Model of Spiral Phyllotaxis”, D. R. Fowler und andere, Proceedings of SIGGRAPH 1992, in Computer Graphics, July 1992, ACM SIGGRAPH, New York, Seiten 361–368, verwiesen wird.
  • Bei der Spiral-Phyllotaxis ist ein Divergenzwinkel, welchen ein beliebiger Punkt zu seinem Vorgänger-Punkt aufweist, immer konstant. Der Divergenzwinkel ist dabei bezüglich eines Mittelpunktes definiert, um welchen sich eine Spirale, welche von den Punkten gebildet wird, windet. Anders ausgedrückt beschreibt der Divergenzwinkel einen Winkel, welcher in einem Dreieck, welches aus dem Punkt, seinem Vorgängerpunkt und dem Mittelpunkt aufgebaut ist, am Mittelpunkt auftritt.
  • Bei der vorliegenden Erfindung liegen die Punkte auf einer Kugeloberfläche. Bestimmt man die Lage eines beliebigen dieser Punkte mittels Kugelkoordinaten (Koordinatenursprung = Mittelpunkt der Kugel), dann gilt für den Fall, dass die Punkte auf der Kugel gemäß der Spiral-Phyllotaxis angeordnet sind, dass ein Differenzwinkel zwischen dem Azimutwinkel eines beliebigen dieser Punkte und dem Azimutwinkel seines Vorgänger-Punktes immer denselben Wert aufweist bzw. konstant ist. Dieser Differenzwinkel kann auch als der Divergenzwinkel angesehen werden.
  • Die Fibonacci-Phyllotaxis ist ein Spezialfall der Spiral-Phyllotaxis, wobei der Divergenzwinkel bei der Fibonacci-Phyllotaxis einem der goldenen Winkel (siehe Gleichungen (1) und (2) unten) entspricht.
  • In diesem Zusammenhang sei zum besseren Verständnis insbesondere der Fibonacci-Phyllotaxis auch auf „A Better Way to Construct the Sunflower Head, H. Vogel, Mathematical Biosiences 44: Seiten 179–189, 1979, verwiesen.
  • Indem die Punkte zur Definition der Speichen auf der Kugel bzw. Kugeloberfläche entsprechend der Spiral-Phyllotaxis, insbesondere entsprechend der Fibonacci-Phyllotaxis, angeordnet werden, ist die Verteilung der Abtastpunkte gleichmäßiger als es bei dem Ansatz der Fall ist, welcher in „Temporal Stability of Adaptive 3D Radial MRI Using Multidimensional Golden Means”, R. W. Chan und andere, Magnetic Resonance in Medicine 61, Seiten 354–363, 2009, beschrieben ist. Gegenüber dem Ansatz, bei welchem die Punkte auf der Kugel in Form einer archimedischen Spirale angeordnet sind, ergibt sich vorteilhafterweise eine Verringerung der Wirbelstromeffekte.
  • Um die Punkte zur Definition der Speichen auf der Kugel entsprechend der Spiral-Pyhllotaxis, insbesondere der Fibonacci-Phyllotaxis, anzuordnen, können die Punkte entweder mittels des Voronoi-Diagramms oder mittels des „Contact Pressure Modells” (siehe „Contact Pressure models for spiral Phyllotaxis and their computer simulation”, H. Hellwig und andere, Journal of Theoretical Biology 240 (2006), Seiten 489–500) angeordnet werden.
  • Ein Voronoi-Diagramm bezeichnet eine Zerlegung des Raumes in Regionen, die durch eine vorgegebene Menge an Punkten des Raumes, welche als Zentren bezeichnet werden, bestimmt werden. Jede Region wird durch genau ein Zentrum bestimmt und umfasst alle Punkte des Raumes, die in Bezug zur euklidischen Metrik näher an dem Zentrum der Region liegen als an jedem anderen Zentrum. Derartige Regionen werden auch als Voronoi-Regionen bezeichnet. Aus allen Punkten, die mehr als ein nächstgelegenes Zentrum besitzen und somit die Grenzen der Regionen bilden, entsteht das Voronoi-Diagramm.
  • Gemäß einer bevorzugten erfindungsgemäßen Ausführungsform werden N Speichen zur Datenerfassung des dreidimensionalen K-Raums abgearbeitet, d. h. die Anzahl der abgearbeiteten Speichen beträgt N. Jeder Punkt auf der Kugel, durch welchen eine dieser Speichen definiert ist, wird vom Mittelpunkt der Kugel durch einen Polarwinkel und einen Azimutwinkel definiert. Mit anderen Worten beschreibt der Polarwinkel unter der Annahme eines kartesischen Koordinatensystems mit dem Ursprung gleich dem Mittelpunkt der Kugel denjenigen Winkel zwischen der positiven x-Achse dieses Koordinatensystems und der durch den entsprechenden Punkt definierten Speiche während der Azimutwinkel den Winkel zwischen der positiven z-Achse dieses Koordinatensystems und dieser Speiche definiert. Dabei ist der Polarwinkel der Speiche n oder n-ten Speiche abhängig von (z. B. proportional zu) der Wurzel aus n bzw. √ n oder abhängig von (z. B. proportional zu) der Wurzel aus der Anzahl der Speichen (N) bzw. √ N . Es ist auch möglich, dass der Polarwinkel sowohl abhängig von (z. B. proportional zu) √ n als auch abhängig von (z. B. proportional zu) √ N ist.
  • Gemäß einer weiteren erfindungsgemäßen Ausführungsform ist der Azimutwinkel der n-ten Speiche abhängig von (z. B. proportional zu) n oder abhängig von (z. B. proportional zu) dem goldenen Winkel. Darüber hinaus ist es möglich, dass der Azimutwinkel sowohl abhängig von (z. B. proportional zu) n als auch abhängig von (z. B. proportional zu) einem der goldenen Winkel ist. Dabei wird bei der vorliegenden Erfindung unter einem goldenen Winkel Ψ entweder der so genannte kleine goldene Winkel Ψ1 oder der so genannte große goldene Winkel Ψ2 verstanden. Der kleine goldene Winkel Ψ1 ist gemäß der folgenden Gleichung (1) definiert:
    Figure 00090001
  • Der große goldene Winkel Ψ2 ist gemäß der folgenden Gleichung (2) definiert:
    Figure 00090002
  • Gemäß einer bevorzugten erfindungsgemäßen Ausführungsform wird der Polarwinkel θn der Speiche n bzw. der n-ten Speiche entsprechend der folgenden Gleichung (3) ermittelt:
    Figure 00090003
  • Gemäß einer weiteren bevorzugten erfindungsgemäßen Ausführungsform wird der Azimutwinkel Φn der Speiche n bzw. der n-ten Speiche entsprechend der folgenden Gleichung (4) ermittelt: Φn = mod(n × Ψ, 2 × π) (4)
  • Mit anderen Worten existieren erfindungsgemäß zwei unterschiedliche Definitionen des Azimutwinkels Φn. Bei der ersten Definition wird der Azimutwinkel Φn mit dem kleinen goldenen Winkel definiert (d. h. Ψ = Ψ1) und bei der zweiten Definition wird der Azimutwinkel mit dem großen goldenen Winkel definiert (d. h. Ψ = Ψ2). Dabei wird insbesondere bei der Definition der Punkte auf der Kugeloberfläche und damit der Speichen entweder nur mit der ersten Definition des Azimutwinkels Φn oder nur mit zweiten Definition des Azimutwinkels Φn gearbeitet.
  • Mit Hilfe der Gleichungen (3) und (4) kann sowohl der Polarwinkel als auch der Azimutwinkel in einer geschlossenen Form nur abhängig von der Gesamtzahl der Speichen N und abhängig von dem jeweiligen Index n der jeweils zu bestimmenden Speiche berechnet werden, wodurch dann die jeweilige Speiche definiert ist. Dadurch wird eine zeitaufwändige Iteration oder Optimierung zur Bestimmung der Verteilung der Punkte auf der Kugeloberfläche vorteilhafterweise vermieden.
  • Erfindungsgemäß ist es auch möglich, dass in Gleichung (4) anstelle des kleinen goldenen Winkels Ψ1 bzw. anstelle des großen goldenen Winkels ein kleiner Winkel Ψ1' bzw. ein großer Winkel Ψ2' eingesetzt wird, welcher im Wesentlichen dem kleinen goldenen Winkel bzw. dem großen goldenen Winkel entspricht, indem der kleine Winkel Ψ1' der folgenden Gleichung (5) und der große Winkel Ψ2' der folgenden Gleichung (6) gehorcht: 0,9·Ψ1 ≤ Ψ1' ≤ 1,1·Ψ1 (5) 0,9·Ψ2 ≤ Ψ2' ≤ 1,1·Ψ2 (6)
  • Gemäß einer erfindungsgemäßen Ausführungsform werden alle Speichen, welche zur radialen Datenerfassung abzuarbeiten sind, in M Mengen unterteilt. Entsprechend einer vorab bestimmten Reihenfolge der Speichen werden die Speichen derart einer der M Mengen zugeteilt, dass die k-te Speiche, die (k + M)-te Speiche, die (k + 2·M)-te Speiche, ... zu der k-ten Menge gehört. Bei der Datenerfassung werden erst die Daten in dem K-Raum für alle Speichen einer dieser Mengen erfasst, bevor die Daten aller Speichen für die nächste Menge erfasst werden. Dabei wird die Anzahl M der Mengen vorteilhafterweise derart gewählt, dass es eine Fibonacci-Zahl ist. Da die Fibonacci-Zahlen auch die Zahl 1 umfassen, ist es demnach auch möglich, dass alle Speichen zu derselben Menge gehören.
  • Wenn die Anzahl M einer Fibonacci-Zahl entspricht, weisen entsprechend der vorbestimmten Reihenfolge zwei aufeinanderfolgende Speichen derselben Menge vorteilhafterweise einen geringen Abstand auf. Dadurch ist die Änderung zwischen einem Magnetfeldgradienten zur Erfassung der einen dieser beiden Speichen und einem Magnetfeldgradienten zur Erfassung der anderen dieser beiden Speichen vorteilhafterweise gering, so dass nur geringe Wirbelstromeffekte auftreten, so dass auf eine Kompensation dieser Wirbelstromeffekte mittels eines komplizierten Algorithmus verzichtet werden kann.
  • Natürlich ist es auch möglich, dass die Reihenfolge, in welcher die Speichen einer bestimmten Menge (einer der M Mengen) zur Datenerfassung abgearbeitet werden, derart bestimmt wird, dass der Abstand zwischen direkt aufeinanderfolgend abgearbeiteten Speichen minimal ist. Beispielsweise könnte mittels eines Optimierungsverfahrens diejenige Reihenfolge ermittelt werden, in welcher eine Summe der Abstände zwischen zwei aufeinanderfolgenden abgearbeiteten Speichen ein Minimum aufweist. Dabei könnte der Abstand zwischen zwei Speichen durch die Strecke definiert sein, welche die beiden Punkte, über die die jeweilige Speiche definiert ist, auf der Kugeloberfläche verbindet.
  • Gemäß einem erfindungsgemäßen Beispiel könnten alle N Speichen mit Hilfe der Gleichungen (2) und (3) definiert werden, wobei die Reihenfolge durch den Index n festgelegt ist. Mit anderen Worten folgt gemäß dieser vorbestimmten Reihenfolge die Speiche, welche über den Polarwinkel θn+1 und den Azimutwinkel Φn+1 definiert ist, direkt derjenigen Speiche, welche über den Polarwinkel θn und den Azimutwinkel Φn definiert ist. Diese derart definierten Speichen werden nun in M verschiedene Mengen aufgeteilt, so dass in der Menge 1 die Speiche mit dem Index 1 (definiert durch den Polarwinkel θ1 und den Azimutwinkel Φ1), die Speiche mit dem Index (1 + M) (definiert durch den Polarwinkel θ1+M und den Azimutwinkel Φ1+M), die Speiche mit dem Index (1 + 2·M) (definiert durch den Polarwinkel θ1+2·M und den Azimutwinkel Φ1+2·M), die Speiche mit dem Index (1+3·M) (definiert durch den Polarwinkel θ1+3·M und den Azimutwinkel Φ1+3·M), usw. vorhanden ist. In der k-ten Menge befinden sich demnach die Speiche mit dem Index k (definiert durch den Polarwinkel θk und den Azimutwinkel Φk), die Speiche mit dem Index (k + M) (definiert durch den Polarwinkel θk+M und den Azimutwinkel Φk+M), die Speiche mit dem Index (k+2·M) (definiert durch den Polarwinkel θk+2·M und den Azimutwinkel Φk+2·M), die Speiche mit dem Index (k + 3·M) (definiert durch den Polarwinkel θk+3·M und den Azimutwinkel Φk+3·M), usw..
  • Wenn die Anzahl M einer Fibonacci-Zahl entspricht, werden die Speichen einer beliebigen der M Mengen in einer Reihenfolge abgearbeitet, welche durch den Index gegeben ist. Mit anderen Worten wird für die k-te Menge zuerst die Speiche mit dem Index k, anschließend die Speiche mit dem Index (k + M), dann die Speiche mit dem Index (k + 2·M), dann die Speiche mit dem Index (k + 3·M) usw. abgearbeitet. Aufgrund der Tatsache, dass die Speichen mit Hilfe der Gleichungen (2) und (3) definiert sind, ist dann gewährleistet, dass aufgrund der geringen Gradientenänderungen zwischen zwei aufeinanderfolgenden Speichen nahezu keine Wirbelstromeffekte auftreten.
  • Wenn dagegen die Anzahl M nicht einer Fibonacci-Zahl entspricht, existieren in der Regel große Gradientenänderungen zwischen zwei aufeinanderfolgenden Speichen, so dass es ratsam ist, die Reihenfolge, in welcher die Speichen einer Menge abgearbeitet werden, nicht gleich derjenigen Reihenfolge zu definieren, welche durch den Index gegeben ist. In diesem Fall kann die Reihenfolge derart gewählt werden, dass die Abstände zwischen zwei aufeinander folgenden Speichen jeweils möglichst klein ist, was natürlich eine aufwendige Optimierung bedeuten kann.
  • Im Rahmen der vorliegenden Erfindung wird auch eine Vorrichtung für eine Magnetresonanzanlage zur radialen Datenerfassung in einem dreidimensionalen K-Raum bei einer MR-Messung bereitgestellt. Dabei umfasst die Vorrichtung eine Ansteuereinheit zur Ansteuerung der Magnetresonanzanlage, eine Empfangsvorrichtung zum Empfang von von der Magnetresonanzanlage aufgenommenen MR-Daten und eine Auswertevorrichtung zur Auswertung der MR-Daten. Die Vorrichtung ist derart ausgestaltet, dass die Vorrichtung die Magnetresonanzanlage mittels der Ansteuervorrichtung derart ansteuert, dass Daten in dem K-Raum entlang von geradlinig verlaufenden Speichen erfasst werden. Dabei ist jede Speiche durch einen Punkt auf einer Kugel und den Mittelpunkt der Kugel definiert, wobei der Mittelpunkt der Kugel dem Zentrum des K-Raums entspricht. Die Vorrichtung ordnet die Punkte auf der Kugeloberfläche derart an, dass die Verteilung der Punkte der Spiral-Phyllotaxis, noch besser der Fibonacci-Phyllotaxis, gehorcht.
  • Die Vorteile der erfindungsgemäßen Vorrichtung entsprechen im Wesentlichen den Vorteilen des erfindungsgemäßen Verfahrens, welche vorab im Detail ausgeführt sind, weshalb hier auf eine Wiederholung verzichtet wird.
  • Darüber hinaus offenbart die vorliegende Erfindung eine Magnetresonanzanlage, welche die vorab beschriebene erfindungsgemäße Vorrichtung umfasst.
  • Des Weiteren beschreibt die vorliegende Erfindung ein Computerprogrammprodukt, insbesondere eine Software, welche man in einen Speicher einer programmierbaren Steuerung bzw. einer Recheneinheit einer Magnetresonanzanlage laden kann. Mit diesem Computerprogrammprodukt können alle oder verschiedene vorab beschriebene Ausführungsformen des erfindungsgemäßen Verfahrens ausgeführt werden, wenn das Computerprogrammprodukt in der Steuerung läuft. Dabei benötigt das Computerprogrammprodukt eventuell Programmmittel, z. B. Bibliotheken und Hilfsfunktionen, um die entsprechenden Ausführungsformen der Verfahren zu realisieren. Mit anderen Worten soll mit dem auf das Computerprogrammprodukt gerichteten Anspruch insbesondere eine Software unter Schutz gestellt werden, mit welcher eine der oben beschriebenen Ausführungsformen des erfindungsgemäßen Verfahrens ausgeführt werden kann. Dabei kann es sich bei der Software um einen Quellcode (z. B. in C++), der noch compiliert und gebunden oder der nur interpretiert werden muss, oder um einen ausführbaren Softwarecode handeln, der zur Ausführung nur noch in die entsprechende Recheneinheit zu laden ist.
  • Schließlich offenbart die vorliegende Erfindung einen elektronisch lesbaren Datenträger, z. B. eine DVD, ein Magnetband oder einen USB-Stick, auf welchem elektronisch lesbare Steuerinformationen, insbesondere Software (vgl. oben), gespeichert ist. Wenn diese Steuerinformationen (Software) von dem Datenträger gelesen und in eine Steuerung bzw. Recheneinheit einer Magnetresonanzanlage gespeichert werden, können alle erfindungsgemäßen Ausführungsformen des vorab beschriebenen Verfahrens durchgeführt werden.
  • Die vorliegende Erfindung ist insbesondere zur radialen Datenerfassung in einem dreidimensionalen K-Raum bei einer MR-Messung (Magnetresonanz-Messung) geeignet. Selbstverständlich ist die vorliegende Erfindung nicht auf diesen bevorzugten Anwendungsbereich beschränkt, da die vorliegende Erfindung beispielsweise auch bei anderen bildgebenden Verfahren (z. B. Röntgen-Computertomographie, Positronen-Emissionstomographie) eingesetzt werden könnte.
  • Im Folgenden wird die vorliegende Erfindung anhand bevorzugter erfindungsgemäßer Ausführungsformen im Detail erläutert.
  • 1 stellt eine erste Ausführungsform zur radialen Datenerfassung in einem dreidimensionalen K-Raum nach dem Stand der Technik dar.
  • 2 stellt eine zweite Ausführungsform zur radialen Datenerfassung in einem dreidimensionalen K-Raum nach dem Stand der Technik dar.
  • 3a und 3b stellen eine Verteilung von Punkten zur Definition von Speichen auf einer Kugeloberfläche gemäß einer erfindungsgemäßen Ausführungsform dar.
  • 3c stellt eine Trajektorie gemäß einer nicht erfindungsgemäßen Abarbeitung der entsprechend der 3a und 3b definierten Speichen dar.
  • 3d stellt eine Trajektorie gemäß einer erfindungsgemäßen Abarbeitung der entsprechend der 3a und 3b definierten Speichen dar.
  • In 4 sind Ergebnisse von verschiedenen Phantom-Experimenten dargestellt, welche die Auswirkungen des Wirbelstroms bei einer dreidimensionalen True-FISP-Sequenz für unterschiedliche Ausführungsformen einer radialen Datenerfassung in einem dreidimensionalen K-Raum zeigen.
  • 5 stellt eine erfindungsgemäße Magnetresonanzanlage mit einer erfindungsgemäßen Vorrichtung dar.
  • In 3a und 3b sind Punkte 1 auf einer Kugeloberfläche 4 dargestellt, welche mit den vorab beschriebenen Gleichungen (2) und (3) berechnet worden sind. Dabei wurde N mit 1600 gewählt, so dass 1600 Punkte 1 auf der Kugeloberfläche 4 vorhanden sind, welche 1600 Speichen definieren. Die 3a zeigt die Kugel 4 von oben, während die 3b die Kugel 4 schräg von oben darstellt. Aus den 3a und 3b ist die gleichmäßige Verteilung der Punkte 1 auf der Kugeloberfläche 4 ersichtlich. Durch den Laufindex n in den Gleichungen (2) und (3) ist auch eine implizite Reihenfolge der 1600 Punkte 1 definiert.
  • Wenn die 1600 Speichen in 20 Mengen oder Gruppen aufgeteilt werden, wobei der k-ten Menge die k-te Speiche, die (k + 20)-te Speiche, die (k + 2·20)-te Speiche, die (k + 3·20)-te Speiche, usw. zugeteilt wird, ergibt sich für eine dieser Mengen die in 3c dargestellte Trajektorie 15. Dabei verbindet die in 3c dargestellte Trajektorie 15 die Punkte 1 der entsprechenden Menge gemäß der oben erwähnten impliziten Reihenfolge. Man erkennt, dass die geradlinigen Abschnitte der Trajektorie 15, deren Länge ein Maß einer Änderung der Magnetfeldgradienten darstellt, relativ lang sind, was nachteiligerweise zu deutlichen Wirbelstromeffekten führt. Der Grund für diese nachteiligen Wirbelstromeffekte liegt in der Wahl der Anzahl der Mengen, da die gewählte Anzahl 20 keine Fibonacci-Zahl ist.
  • Wenn die 1600 Speichen dagegen in 21 Mengen aufgeteilt werden, wobei der k-ten dieser 20 Mengen die k-te Speiche, die (k + 21)-te Speiche, die (k + 2·21)-te Speiche, die (k + 3·21)-te Speiche, usw. zugeteilt wird, ergibt sich für eine dieser Mengen die in 3d dargestellte Trajektorie 15. Man erkennt, dass die geradlinigen Abschnitte dieser Trajektorie 15 in 3d wesentlich kürzer sind als bei der Trajektorie 15 in 3c. Daher führt eine Abarbeitung der Speichen, welche durch die Punkte 1 definiert sind, gemäß der in 3d dargestellten Trajektorie 15 zu wesentlich geringeren Wirbelstromeffekten als es bei der in 3c dargestellten Trajektorie 15 der Fall ist. Der Grund ist einfach, dass es sich bei der gewählten Anzahl 21 um eine Fibonacci-Zahl handelt.
  • In 4 sind für verschiedene Ausführungsformen einer radialen Datenerfassung in einem dreidimensionalen K-Raum Ergebnisse für Phantom-Experimente dargestellt.
  • In 4a ist das Ergebnis für 1588 Speichen dargestellt, welche mittels Punkten auf der Kugeloberfläche definiert sind, welche entsprechend einer archimedischen Spirale auf der Kugeloberfläche angeordnet sind. Die Speichen wurden dabei nicht verschachtelt erfasst, was bedeutet, dass die 1588 Speichen zur Datenerfassung nicht in unterschiedliche Mengen aufgeteilt wurden. Man erkennt, dass in 4a keine Artefakte zu erkennen sind, da zwei zur Datenerfassung aufeinanderfolgende Speichen aufgrund der nicht vorhandenen Verschachtelung jeweils einen äußerst geringen Abstand aufweisen, so dass nahezu keine Wirbelstromeffekte auftreten.
  • Wenn die 1588 Speichen dagegen in 377 Mengen unterteilt werden und die Speichen einer dieser Mengen in der von der archimedischen Spirale vorgegebenen Reihenfolge zur Datenerfassung abgearbeitet werden, ergibt sich das in 4b dargestellte Bild. Man erkennt in 4b starke Artefakte, welche durch große Wirbelstromeffekte hervorgerufen werden. Dabei ist anzumerken, dass es sich bei der Anzahl der Mengen von 377 um eine Fibonacci-Zahl handelt.
  • Für das in 4c dargestellte Ergebnis wurden die Punkte auf der Kugeloberfläche, welche die Speichen definieren, dagegen erfindungsgemäß angeordnet, so dass eine Verteilung der Punkte der Fibonacci-Phyllotaxis gehorcht. Die Speichen wurden ebenfalls in 377 Mengen aufgeteilt, und die radiale Datenerfassung erfolgte, indem die Speichen Menge für Menge (für alle 377 Mengen) abgearbeitet wurden. Das daraus resultierende Ergebnis, welches in 4c dargestellt ist, zeigt, dass nahezu kaum wahrnehmbare Artefakte aufgrund von Wirbelstromeffekten aufgetreten sind. Es zeigt sich also, dass die Wahl einer Fibonacci-Zahl für die Anzahl der Mengen nur dann einen Vorteil bringt, wenn die Punkte auf der Kugeloberfläche, welche die Speichen zur radialen Datenerfassung definieren, eine der Fibonacci-Phyllotaxis gehorchende Verteilung aufweisen, wie es erfindungsgemäß gefordert ist.
  • In 5 ist eine erfindungsgemäße Magnetresonanzanlage 5 schematisch dargestellt. Die Magnetresonanzanlage 5 umfasst im Wesentlichen einen Tomograph 3, mit welchem das für die MR-Untersuchung notwendige Magnetfeld in einem Messraum 4 erzeugt wird, einen Tisch 2, eine Steuereinrichtung 6, mit welcher der Tomograph 3 gesteuert wird und MR-Daten von dem Tomograph 3 erfasst werden, und ein an die Steuereinrichtung 6 angeschlossenes Terminal 7.
  • Die Steuereinrichtung 6 umfasst ihrerseits eine Ansteuereinheit 11, eine Empfangsvorrichtung 12 und eine Auswertevorrichtung 13. Während einer MR-Untersuchung werden MR-Daten mittels des Tomograph 3 von der Empfangsvorrichtung 12 erfasst, wobei der Tomograph 3 und der Tisch 2 von der Ansteuereinheit 11 derart angesteuert werden, dass MR-Daten in einem Messvolumen, welches sich im Körperinneren eines auf dem Tisch 2 liegenden Patienten O befindet, erfasst werden.
  • Die Auswertevorrichtung 13 bereitet dann die MR-Daten derart auf, dass sie auf einem Bildschirm 8 des Terminals 7 grafisch dargestellt werden können und dass erfindungsgemäß erstellte Bilder angezeigt werden. Neben der grafischen Darstellung der MR-Daten kann mit dem Terminal 7, welches neben dem Bildschirm 8 eine Tastatur 9 und eine Maus 10 umfasst, von einem Anwender z. B. ein zu vermessender dreidimensionaler Volumenabschnitt vorgegeben werden und weitere Parameter zur Durchführung des erfindungsgemäßen Verfahrens bestimmt werden. Über das Terminal 7 kann auch die Software für die Steuereinrichtung 6 in die Steuereinrichtung 6, insbesondere in die Auswertevorrichtung 13, geladen werden. Diese Software der Steuereinrichtung 6 kann dabei auch das erfindungsgemäße Verfahren umfassen. Es ist dabei auch möglich, dass ein erfindungsgemäßes Verfahren in einer Software enthalten ist, welche in dem Terminal 7 abläuft. Unabhängig davon, in welcher Software das erfindungsgemäße Verfahren enthalten ist, kann die Software auf einer DVD 14 gespeichert sein, so dass diese Software dann von dem Terminal 7 von der DVD 14 gelesen und entweder in die Steuereinrichtung 6 oder in eine Recheneinheit des Terminals 7 selbst kopiert werden kann.

Claims (26)

  1. Verfahren zur radialen Datenerfassung in einem dreidimensionalen K-Raum bei einer MR-Messung mittels einer Magnetresonanzanlage (5), wobei Daten im K-Raum entlang von geradlinigen Speichen erfasst werden, wobei jede der Speichen durch einen Punkt (1) auf einer Kugeloberfläche (4) und den Mittelpunkt dieser Kugel definiert ist, und wobei dieser Mittelpunkt dem Zentrum des zu erfassenden K-Raums entspricht, dadurch gekennzeichnet, dass die Punkte (1) auf der Kugeloberfläche (4) derart angeordnet werden, dass die Verteilung der Punkte (1) der Spiral-Phyllotaxis gehorcht.
  2. Verfahren nach Anspruch 1, dadurch gekennzeichnet, dass die Punkte (1) auf der Kugeloberfläche (4) derart angeordnet werden, dass die Verteilung der Punkte (1) der Fibonacci-Phyllotaxis gehorcht.
  3. Verfahren nach Anspruch 1 oder 2, dadurch gekennzeichnet, dass die Punkte (1) auf der Kugeloberfläche (4) mittels des Voronoi-Diagramms angeordnet werden.
  4. Verfahren nach einem der vorhergehenden Ansprüche, dadurch gekennzeichnet, dass die Punkte (1) auf der Kugeloberfläche (4) mittels des „Contact Pressure Modells” angeordnet werden.
  5. Verfahren nach einem der vorhergehenden Ansprüche, dadurch gekennzeichnet, dass die Anzahl von Speichen, welche zur Datenerfassung abgearbeitet werden, gleich N ist, dass jeder Punkt (1) auf der Kugeloberfläche (4) durch einen Polarwinkel und einen Azimutwinkel definiert ist, dass der Polarwinkel einer n-ten Speiche abhängig von der Wurzel aus n und/oder abhängig von der Wurzel aus der Anzahl N der Speichen gebildet wird, und/oder dass der Azimutwinkel der n-ten Speiche abhängig von n und/oder abhängig von dem goldenen Winkel gebildet wird, wobei der goldene Winkel (Ψ) entweder dem kleinen goldenen Winkel (Ψ1) oder dem großen goldenen Winkel (Ψ2) entspricht.
  6. Verfahren nach einem der Ansprüche 1–4, dadurch gekennzeichnet, dass die Anzahl von Speichen, welche zur Datenerfassung abgearbeitet werden, gleich N ist, dass jeder Punkt (1) auf der Kugeloberfläche (4) durch einen Polarwinkel und einen Azimutwinkel definiert ist, und dass der Polarwinkel der n-ten Speiche proportional zu der Wurzel aus n und/oder proportional zu dem Kehrwert der Wurzel aus der Anzahl N der Speichen gebildet wird.
  7. Verfahren nach Anspruch 5 oder 6, dadurch gekennzeichnet, dass der Polarwinkel θn der n-ten Speiche wie folgt ermittelt wird:
    Figure 00200001
  8. Verfahren nach einem der Ansprüche 5–7, dadurch gekennzeichnet, dass der Azimutwinkel der n-ten Speiche proportional zu n und/oder proportional zu dem goldenen Winkel ist.
  9. Verfahren nach einem der Ansprüche 5–7, dadurch gekennzeichnet, dass der Azimutwinkel Φn der n-ten Speiche abhängig von dem goldenen Winkel Ψ wie folgt ermittelt wird: Φn = mod(n × Ψ, 2 × π).
  10. Verfahren nach einem der vorhergehenden Ansprüche, dadurch gekennzeichnet, dass die Speichen in eine Anzahl M von Mengen aufgeteilt werden, dass einer k-ten dieser Mengen gemäß einer vorab bestimmten Reihenfolge der Speichen jeweils die k-te, (k + M)-te, (k + 2·M)-te usw. Speiche zugeordnet wird, und dass erst die Daten in dem K-Raum entlang aller Speichen einer jeweiligen der Mengen erfasst werden, bevor die Daten in dem K-Raum entlang der Speichen einer anderen der Mengen erfasst werden.
  11. Verfahren nach Anspruch 10, dadurch gekennzeichnet, dass die Anzahl M entsprechend einer Fibonacci-Zahl gewählt wird, und dass eine weitere Reihenfolge (15), in welcher die Speichen einer der M Mengen zur Datenerfassung abgearbeitet werden, derart gewählt wird, dass diese weitere Reihenfolge der vorab bestimmten Reihenfolge entspricht.
  12. Verfahren nach Anspruch 10, dadurch gekennzeichnet, dass eine weitere Reihenfolge (15), in welcher die Speichen einer der M Mengen zur Datenerfassung abgearbeitet werden, derart gewählt wird, dass der Abstand zwischen direkt aufeinander folgenden Speichen möglichst gering ist.
  13. Vorrichtung für eine Magnetresonanzanlage (5) zur radialen Datenerfassung in einem dreidimensionalen K-Raum bei einer MR-Messung, wobei die Vorrichtung (6) eine Ansteuereinheit (11) zur Ansteuerung der Magnetresonanzanlage (5), eine Empfangsvorrichtung (12) zum Empfang von von der Magnetresonanzanlage (5) aufgenommenen MR-Daten und eine Auswertevorrichtung (13) zur Auswertung der MR-Daten umfasst, dass die Vorrichtung (6) derart ausgestaltet ist, dass sie die Magnetresonanzanlage (5) über die Ansteuervorrichtung (11) derart ansteuert, dass Daten in dem K-Raum entlang von geradlinigen Speichen erfasst werden, wobei jede der Speichen durch einen Punkt (1) auf einer Kugeloberfläche (4) und den Mittelpunkt der Kugel definiert ist, und wobei der Mittelpunkt dem Zentrum des K-Raums entspricht, dadurch gekennzeichnet, dass Vorrichtung (6) derart ausgestaltet ist, dass sie die Punkte (1) auf der Kugeloberfläche (4) derart anordnet, dass eine Verteilung der Punkte (1) der Spiral-Phyllotaxis gehorcht.
  14. Vorrichtung nach Anspruch 13, dadurch gekennzeichnet, dass die Vorrichtung (6) derart ausgestaltet ist, dass sie die Punkte (1) auf der Kugeloberfläche (4) derart anordnet, dass die Verteilung der Punkte (1) der Fibonacci-Phyllotaxis gehorcht.
  15. Vorrichtung nach Anspruch 13 oder 14, dadurch gekennzeichnet, dass die Vorrichtung (6) derart ausgestaltet ist, dass sie die Punkte (1) auf der Kugeloberfläche (4) mittels des Voronoi-Diagramms anordnet.
  16. Vorrichtung nach einem der Ansprüche 13–15, dadurch gekennzeichnet, dass die Vorrichtung (6) derart ausgestaltet ist, dass sie die Punkte (1) auf der Kugeloberfläche (4) mittels des „Contact Pressure Modells” anordnet.
  17. Vorrichtung nach einem der Ansprüche 13–16, dadurch gekennzeichnet, dass eine Anzahl von Speichen, welche zur Datenerfassung abgearbeitet werden, gleich N ist, dass jeder Punkt (1) auf der Kugeloberfläche (4) durch einen Polarwinkel und einen Azimutwinkel definiert ist, und dass die Vorrichtung (6) derart ausgestaltet ist, dass sie den Polarwinkel einer n-ten Speiche abhängig von der Wurzel aus n und/oder abhängig von der Wurzel aus der Anzahl N der Speichen bildet, und/oder dass die Vorrichtung (6) den Azimutwinkel der n-ten Speiche abhängig von n und/oder abhängig von dem goldenen Winkel bildet, wobei der goldene Winkel (Ψ) entweder dem kleinen goldenen Winkel (Ψ1) oder dem großen goldenen Winkel (Ψ2) entspricht.
  18. Vorrichtung nach einem der Ansprüche 13–16, dadurch gekennzeichnet, dass eine Anzahl von Speichen, welche zur Datenerfassung abgearbeitet werden, gleich N ist, dass jeder Punkt (1) auf der Kugeloberfläche (4) durch einen Polarwinkel und einen Azimutwinkel definiert ist, und dass die Vorrichtung (6) derart ausgestaltet ist, dass sie den Polarwinkel der n-ten Speiche proportional zu der Wurzel aus n und/oder proportional zu dem Kehrwert der Wurzel aus der Anzahl N der Speichen bildet.
  19. Vorrichtung nach Anspruch 17 oder 18, dadurch gekennzeichnet, dass die Vorrichtung (6) derart ausgestaltet ist, dass sie den Polarwinkel θn der n-ten Speiche wie folgt ermittelt:
    Figure 00230001
  20. Vorrichtung nach einem der Ansprüche 17–19, dadurch gekennzeichnet, dass die Vorrichtung (6) derart ausgestaltet ist, dass sie den Azimutwinkel der n-ten Speiche proportional zu n und/oder proportional zu dem goldenen Winkel bildet.
  21. Vorrichtung nach einem der Ansprüche 17–20, dadurch gekennzeichnet, dass die Vorrichtung (6) derart ausgestaltet ist, dass sie den Azimutwinkel Φn der n-ten Speiche abhängig von dem goldenen Winkel Ψ wie folgt ermittelt: Φn = mod(n × Ψ, 2 × π).
  22. Vorrichtung nach einem der Ansprüche 13–21, dadurch gekennzeichnet, dass die Vorrichtung (6) derart ausgestaltet ist, dass sie die Speichen in eine Anzahl M von Mengen derart aufteilt, dass sie einer k-ten dieser Mengen gemäß der gewählten Reihenfolge jeweils die k-te, (k + M)-te, (k + 2·M)-te usw. Speiche zuordnet, und dass sie ferner erst die Daten im K-Raum entlang aller Speichen einer jeweiligen der Mengen erfasst, bevor die Vorrichtung (6) die Daten in dem K-Raum entlang der Speichen einer anderen der Mengen erfasst.
  23. Vorrichtung nach Anspruch 22, dadurch gekennzeichnet, dass die Vorrichtung (6) derart ausgestaltet ist, dass sie die Anzahl M entsprechend einer Fibonacci-Zahl wählt und dass sie eine weitere Reihenfolge (15), in welcher sie die Speichen einer der M Mengen zur Datenerfassung abarbeitet, derart wählt, dass diese weitere Reihenfolge der vorab bestimmten Reihenfolge entspricht.
  24. Magnetresonanzanlage mit einer Vorrichtung (6) nach einem der Ansprüche 13–23.
  25. Computerprogrammprodukt, welches direkt in einen Speicher einer programmierbaren Steuereinrichtung (6) einer Magnetresonanzanlage (5) ladbar ist, mit Programm-Mitteln, um alle Schritte eines Verfahrens nach einem der Ansprüche 1–12 auszuführen, wenn das Programm in der Steuereinrichtung (6) der Magnetresonanzanlage (5) ausgeführt wird.
  26. Elektronisch lesbarer Datenträger mit darauf gespeicherten elektronisch lesbaren Steuerinformationen, welche derart ausgestaltet sind, dass sie bei Verwendung des Datenträgers (14) in einer Steuereinrichtung (6) einer Magnetresonanzanlage (5) ein Verfahren nach einem der Ansprüche 1–12 durchführen.
DE102009050662A 2009-10-26 2009-10-26 Radiale MR-Datenerfassung in einem dreidimensionalen K-Raum mit einer Anordnung der Speichen gemäß einer Spiral-Phyllotaxis Active DE102009050662B4 (de)

Priority Applications (2)

Application Number Priority Date Filing Date Title
DE102009050662A DE102009050662B4 (de) 2009-10-26 2009-10-26 Radiale MR-Datenerfassung in einem dreidimensionalen K-Raum mit einer Anordnung der Speichen gemäß einer Spiral-Phyllotaxis
US12/912,096 US8648594B2 (en) 2009-10-26 2010-10-26 Method and device for uniform radial data acquisition in three-dimensional K-space in an MR measurement for a magnetic resonance system

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
DE102009050662A DE102009050662B4 (de) 2009-10-26 2009-10-26 Radiale MR-Datenerfassung in einem dreidimensionalen K-Raum mit einer Anordnung der Speichen gemäß einer Spiral-Phyllotaxis

Publications (2)

Publication Number Publication Date
DE102009050662A1 DE102009050662A1 (de) 2011-04-28
DE102009050662B4 true DE102009050662B4 (de) 2012-09-27

Family

ID=43796820

Family Applications (1)

Application Number Title Priority Date Filing Date
DE102009050662A Active DE102009050662B4 (de) 2009-10-26 2009-10-26 Radiale MR-Datenerfassung in einem dreidimensionalen K-Raum mit einer Anordnung der Speichen gemäß einer Spiral-Phyllotaxis

Country Status (2)

Country Link
US (1) US8648594B2 (de)
DE (1) DE102009050662B4 (de)

Families Citing this family (16)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
DE102009050662B4 (de) * 2009-10-26 2012-09-27 Siemens Aktiengesellschaft Radiale MR-Datenerfassung in einem dreidimensionalen K-Raum mit einer Anordnung der Speichen gemäß einer Spiral-Phyllotaxis
DE102011080591A1 (de) 2011-08-08 2013-02-14 Siemens Aktiengesellschaft Verfahren zur Verarbeitung medizintechnischer Bilddaten
DE102011086401B4 (de) * 2011-11-15 2013-08-08 Siemens Aktiengesellschaft MR-Datenerfassung in einem dreidimensionalen K-Raum mit einer Anordnung von Datenzeilen gemäß einer Spiral-Phyllotaxis
CN104125876B (zh) 2011-12-31 2018-07-20 圣戈班磨料磨具有限公司 具有开口的不均匀分布的研磨制品
US10429469B2 (en) * 2012-04-11 2019-10-01 Dignity Health System and method for magnetic resonance imaging using three-dimensional, distributed, non-Cartesian sampling trajectories
DE102012222149B4 (de) 2012-12-04 2014-09-18 Siemens Aktiengesellschaft Verfahren zur Erzeugung von Magnetresonanzmessdaten
DE102014203068B4 (de) * 2014-02-20 2015-11-26 Siemens Aktiengesellschaft Pseudo-zufällige Erfassung von MR-Daten eines zweidimensionalen Volumenabschnitts
DE102014207236A1 (de) * 2014-04-15 2015-10-15 Siemens Aktiengesellschaft Verfahren und Steuerungseinrichtung zum Erstellen von Magnetresonanzaufnahmen
JP6072723B2 (ja) * 2014-04-21 2017-02-01 株式会社日立製作所 磁気共鳴イメージング装置、及び画像撮像方法
JP6362950B2 (ja) * 2014-07-24 2018-07-25 ジーイー・メディカル・システムズ・グローバル・テクノロジー・カンパニー・エルエルシー 磁気共鳴装置
WO2016069602A1 (en) 2014-10-31 2016-05-06 The Trustees Of The University Of Pennsylvania A method and device for magnetic resonance imaging data acquisition guided by physiologic feedback
US10281542B2 (en) * 2015-09-30 2019-05-07 General Electric Company Magnetic resonance imaging system and method
CN106308799B (zh) * 2016-11-04 2019-04-23 深圳先进技术研究院 用于三维动态磁共振成像的球形k空间采集方法和装置
CN107703467A (zh) * 2017-08-21 2018-02-16 深圳先进技术研究院 三维动态磁共振成像的采集方法、装置、设备及存储介质
EP3764117A1 (de) * 2019-07-09 2021-01-13 Koninklijke Philips N.V. Verbesserte dreidimensionale radiale mr-bildgebung
US11802926B2 (en) * 2020-11-12 2023-10-31 Oxford University Innovation Limited Methods and systems for optimizing sampling efficiency of dynamically sampled magnetic resonance imaging (MRI)

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4651096A (en) * 1985-02-19 1987-03-17 Resonex Imaging method for nuclear magnetic resonance utilizing balanced gradient fields
US6281681B1 (en) * 1999-01-28 2001-08-28 General Electric Company Magnetic resonance imaging with interleaved Fibonacci spiral scanning
US7245124B2 (en) * 2005-04-12 2007-07-17 Mayo Foundation For Medical Education And Research Under-sampled 3D MRI using a shells k-space sampling trajectory

Family Cites Families (14)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5349294A (en) * 1993-05-12 1994-09-20 Picker International Inc. Two and three-dimensionally selective RF pulses for magnetic resonance imaging
US20030135105A1 (en) * 2000-04-26 2003-07-17 Jack Clifford R. Alignment of multiple MR images using navigator signals
US7127092B2 (en) * 2000-04-26 2006-10-24 Mayo Foundation For Medical Education And Research Reduction of motion artifact in NMR images using spherical navigator signals
DE10119660B4 (de) * 2001-04-20 2006-01-05 Siemens Ag Verfahren zur schnellen Gewinnung eines Magnetresonanzbildes
US7109707B2 (en) * 2003-04-08 2006-09-19 University Of Queensland Computer readable magnetic resonance method calculating on-line gradients during imaging
EP2013634A2 (de) * 2006-04-21 2009-01-14 Koninklijke Philips Electronics N.V. Magnetresonanzvorrichtung und -verfahren
US8095202B2 (en) * 2006-08-15 2012-01-10 Mayo Foundation For Medical Education And Research Accelerated shells trajectory MRI acquisition
DE102006042998B4 (de) * 2006-09-13 2008-07-03 Siemens Ag Messsequenz für die dreidimensionale MR-Bildgebung sowie MR-Gerät
DE102006058316B4 (de) * 2006-12-11 2010-10-14 Siemens Ag Verfahren zur Aufnahme von Bilddaten einer Gefäßwand und Magnet-Resonanz-Gerät hierzu
US8214013B2 (en) * 2007-01-26 2012-07-03 Mayo Foundation For Medical Education And Research Accelerated shells trajectory MRI acquisition
WO2008134891A1 (en) * 2007-05-03 2008-11-13 National Research Counsil Of Canada Rf based spatially selective excitation in mri
US7548062B1 (en) * 2008-03-26 2009-06-16 The Board Of Trustees Of The Leland Stanford Junior University 3D cones gradient waveform generation
DE102009050662B4 (de) * 2009-10-26 2012-09-27 Siemens Aktiengesellschaft Radiale MR-Datenerfassung in einem dreidimensionalen K-Raum mit einer Anordnung der Speichen gemäß einer Spiral-Phyllotaxis
DE102011086369B4 (de) * 2011-11-15 2022-08-11 Siemens Healthcare Gmbh Verfahren zur Erstellung von MR-Angiographiebildern und entsprechende Magnetresonanzanlage

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4651096A (en) * 1985-02-19 1987-03-17 Resonex Imaging method for nuclear magnetic resonance utilizing balanced gradient fields
US6281681B1 (en) * 1999-01-28 2001-08-28 General Electric Company Magnetic resonance imaging with interleaved Fibonacci spiral scanning
US7245124B2 (en) * 2005-04-12 2007-07-17 Mayo Foundation For Medical Education And Research Under-sampled 3D MRI using a shells k-space sampling trajectory

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
R.W. Chan et al.: Temporal stability of adaptive 3D radial MRI using multidimensional golden means. In: Magn. Reson. Med., 61, 2009, S. 354-363. *

Also Published As

Publication number Publication date
DE102009050662A1 (de) 2011-04-28
US20110095762A1 (en) 2011-04-28
US8648594B2 (en) 2014-02-11

Similar Documents

Publication Publication Date Title
DE102009050662B4 (de) Radiale MR-Datenerfassung in einem dreidimensionalen K-Raum mit einer Anordnung der Speichen gemäß einer Spiral-Phyllotaxis
DE102004026616B4 (de) Verfahren zum Messen eines Untersuchungsbereichs mit einem Magnetresonanzgerät
DE69811049T2 (de) Elektronisches bildverarbeitungsgerät zur detektion von dimensionnellen änderungen
DE102012212402B3 (de) Ermittlung einer MR-Sequenz unter Berücksichtigung der Energieverteilungsfunktion im k-Raum
DE102017222359B4 (de) Automatische Bestimmung von Korrekturfaktoren für eine Magnetresonanzanlage
EP3120765B1 (de) Mpi-verfahren
DE102011088828B4 (de) Erstellung eines MR-Bildes eines Untersuchungsobjekts unter Verwendung einer für ein Empfangsspulenelement erstellten Maske
EP3485457B1 (de) System, insbesondere magnetresonanzsystem, zum erzeugen von bildern
DE102009014054A1 (de) Verfahren und Vorrichtung zur Steuerung eines Ablaufs einer MR-Messung bei einer Magnetresonanzanlage
DE102015205694B3 (de) MR-Sättigung unter Berücksichtigung der abzubildenden anatomischen Strukturen
EP1882232B1 (de) Verfahren und vorrichtung zum bestimmen von materialgrenzen eines prüfobjekts
EP3649614B1 (de) Verfahren zur bestimmung von unsicherheiten in messdaten aus einer messung eines objekts
WO2018206225A1 (de) Computerimplementiertes verfahren zur vermessung eines objekts aus einer digitalen darstellung des objekts
DE102009033606B4 (de) Bestimmung der Stärke eines magnetischen Störfelds eines Partikels
DE102013205868B4 (de) Verfahren zum Zuordnen von K-Raum-Zeilen zu Echozügen und Verfahren zur Erfassung von MR-Daten sowie entsprechend ausgestaltete Magnetresonanzanlagen
DE102011086401A1 (de) MR-Datenerfassung in einem dreidimensionalen K-Raum mit einer Anordnung von Datenzeilen gemäß einer Spiral-Phyllotaxis
DE19529636A1 (de) Verfahren zur MR-Bildgebung aus mehreren Schichten und Anordnung zur Durchführung des Verfahrens
DE102009055122B4 (de) Verfahren und Vorrichtung zur phasensensitiven Flussmessung mittels einer Magnetresonanzanlage
EP3622477B1 (de) Computerimplementiertes verfahren zur bestimmung einer lokalen abweichung einer geometrie eines objektes von einer soll-geometrie des objektes
DE102016202669A1 (de) Spiralförmiges MR-Erfassungsmuster
DE102016209885B4 (de) Automatische Charakterisierung von Lebergewebe
DE19900737A1 (de) Verfahren zur Korrektur der Meßergebnisse eines Koordinatenmeßgerätes und Koordinatenmeßgerät
DE102007010047A1 (de) Verfahren und Steuereinrichtung zur Positionsbestimmung von Empfangsvorrichtungen für eine Magnetresonanzanlage
DE102015215778B4 (de) Verfahren zu einer Darstellung von Magnetresonanz-Daten
DE102012211799A1 (de) Verfahren zur rechnergestützten Verarbeitung von mit einer Mehrzahl von Elektroden gemessenen Aktionspotentialen des menschlichen oder tierischen Körpers

Legal Events

Date Code Title Description
OP8 Request for examination as to paragraph 44 patent law
R016 Response to examination communication
R018 Grant decision by examination section/examining division
R020 Patent grant now final

Effective date: 20121228

R081 Change of applicant/patentee

Owner name: SIEMENS HEALTHCARE GMBH, DE

Free format text: FORMER OWNER: SIEMENS AKTIENGESELLSCHAFT, 80333 MUENCHEN, DE

R081 Change of applicant/patentee

Owner name: SIEMENS HEALTHINEERS AG, DE

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