-
Feld der Erfindung
-
Diese Patentanmeldung bezieht sich auf Verfahren der Akquisition von Bildern unter Verwendung von Magnetresonanzbildgebung. Insbesondere bezieht sie sich auf automatisierte Verfahren zum Optimieren von Bildakquisitionsparametern, insbesondere bei MRT-Systemen mit Permanentmagneten.
-
Hintergrund der Erfindung
-
Die Magnetresonanzbildgebung (MRT) ist ein Standarddiagnosewerkzeug geworden. Trotz ihrer weitverbreiteten Verwendung ist der eigentliche Betrieb eines MRT-Geräts durch die dem Verfahren inhärente Komplexität aber mehr eine Kunst als eine Wissenschaft geblieben. Die Qualität des während der MRT erhaltenen Bildes hängt in kritischer Weise von der Wahl der Akquisitionsparameter ab, wobei die Optimierung dieser Parameter aber häufig jenseits der Fähigkeiten des durchschnittlichen Bedieners ist. Es gab daher ein erhebliches Bemühen, die Zeit und die Anstrengungen zu reduzieren, die erforderlich sind, um optimale Akquisitionsparameter zu finden, insbesondere durch Automatisierung der Bestimmung dieser Parameter.
-
Beispielsweise offenbart
US-A-4,694,250 ein Verfahren zur Optimierung von Akquisitionsparametern bei dem für gegebene T
1 und T
2 und eine gegebene Protonendichte die Varianz oder Standardabweichung zwischen einem berechneten Bild und dem tatsächlichen Bild als eine Funktion der Scanparameter minimiert wird.
-
US-B-6,781,375 offenbart ein Verfahren zur Optimierung von zumindest einem Scanparameter, bei dem eine Mehrzahl von vorbereitenden Bildern unter Verwendung von verschiedenen Werten eines „Bildqualitätsparameters” erhalten werden. Der Bediener wählt dann das beste Bild aus, und die Scanparameter, die verwendet wurden, um dieses Bild zu erhalten, werden dann für das endgültige Bild verwendet.
-
US-B-7,715,899 offenbart ein Verfahren zur Optimierung von Akquisitionsparametern bei dem ein Vollkörperscan mit niedriger Auflösung durchgeführt wird, ein Bereich von Interesse identifiziert wird und die Akquisitionsparameter für einen nachfolgenden Scan mit hoher Auflösung dann bestimmt werden. Neben dem Erfordernis eines vorläufigen Vollkörperscans offenbart dieses Patent keine Details des Algorithmus, der verwendet wird, um diese Parameter zu finden, oder der Gleichungen, die verwendet werden könnten, um die optimalen Parameter zu finden. Weiterhin wird darin nichts erwähnt über die kritischen Parameter wie das Kontrast-zu-Rauschen-Verhältnis (CNR) oder die Möglichkeit von Voraussagen auf der Grundlage von T
1 oder T
2, gemessen oder abgeschätzt.
-
US-A-2007/0276221 offenbart ein Verfahren zum Erzeugen von MRT-Bildern, das eine Akquise eines Referenzscans umfasst, das Bereitstellen eines Zielwertes eines spezifischen Scanparameters gegenüber dem MRT-Gerät und das Bestimmen eines optimalen Scanparametersatzes entsprechend dem Zielwert dieses spezifischen Scanparameters.
-
Ein automatisiertes Verfahren zum Erhalten von MRT-Bildern, bei dem die Bestimmung von optimalen Akquisitionsparametern nicht auch einen Vollkörperscan oder eine Eingabe des Bedieners erfordert, bleibt daher ein lange bestehendes aber noch unbefriedigtes Bedürfnis.
-
Zusammenfassung der Erfindung
-
Es ist ein Ziel der vorliegenden Erfindung, dieses Bedürfnis zu befriedigen. Insbesondere wird ein Verfahren für eine automatisierte Optimierung von Akquisitionsparametern vorgestellt, um ein Bild eines gewünschten Gewebetyps in Anwesenheit eines zweiten Gewebetyps zu erhalten.
-
Es ist daher ein Ziel der vorliegenden Erfindung, ein Verfahren zum automatischen Bestimmen von optimalen Magnetresonanzbildgebungs-(MRT)-Akquisitionsparametern zum Abbilden einer Probe in einem MRT-Gerät zu offenbaren, die zwei Arten von Gewebe enthält, ein Gewebe A und ein Gewebe B, wobei das besagte Verfahren umfasst: Bestimmen von T
1A, T
2A, T
1B, T
2B, ρ
A und ρ
B, wobei ρ die Dichte von NMR-aktiven Kernen repräsentiert, die untersucht werden; Festlegen von Anfangswerten von T
R und T
E; Bestimmen der Sinalintensitäten S
A und S
B aus der Geleichung S = ρE
1E
2, wobei
und
Berechnen des Kontrast-zu-Rauschen-Verhältnisses für das Gewebe A in Anwesenheit von Gewebe B (CNR
AB) aus der Gleichung
wobei P eine Proportionalitätskonstante ist; und Bestimmen von optimalen Werten von T
R und T
E, die einen Maximalwert von CNR
AB(T
R, T
E) ergeben.
-
Es ist ein zusätzliches Ziel der vorliegenden Erfindung, ein solches Verfahren zu offenbaren, wobei die besagten NMR-aktiven Kerne Protonen sind.
-
In einigen Ausführungsformen der Erfindung umfasst der besagte Schritt des Bestimmens von T1A, T2A, T1B, T2B, ρA und ρB ein Importieren von zumindest einem von T1A, T2A, T1B, T2B, ρA und ρB aus einer Datenbank bekannter Werte. In einigen Ausführungsformen der Erfindung umfasst der besagte Schritt des Bestimmens von T1A, T2A, T1B, T2B, ρA und ρB ein Bestimmen von zumindest einem von T1A, T2A, T1B, T2B, ρA und ρB aus einem vorläufigen MRT-Scan, der an der besagten Probe durchgeführt wurde.
-
Es ist ein weiteres Ziel der vorliegenden Erfindung, das Verfahren wie in irgendeinem des Obigen definiert zu offenbaren, wobei P = 1.
-
Es ist ein weiteres Ziel der vorliegenden Erfindung, das Verfahren wie in irgendeinem des Obigen definiert zu offenbaren, wobei der besagte Schritt des Bestimmens von optimalen Werten von TR und TE umfasst: systematisches Verändern von TR und TE unabhängig voneinander über vorbestimmte Wertebereiche; Speichern von CNRAB(TR, TE) für jedes Paar von Werten TR, TE; und Definieren als optimale Werte TR und TE solche Werte, die den besagten Maximalwert von CNRAB(TR, TE) ergeben. In einigen Ausführungsformen der Erfindung umfasst der besagte Schritt des systematischen Veränderns von TR und TE unabhängig voneinander über vorbestimmte Wertebereiche das Verändern von TE über den Bereich 10–100 ms und das Verändern von TR über den Bereich 0,5–5 s.
-
Es ist ein weiteres Ziel der vorliegenden Erfindung, das Verfahren wie in irgendeinem des Obigen definiert zu offenbaren, wobei der besagte Schritt des Bestimmens von optimalen Werten von TR und TE die Verwendung eines vorprogrammierten Optimierungsalgorithmus umfasst, um den besagten Maximalwert von CNRAB(TR, TE) zu finden. In einigen Ausführungsformen der Erfindung wird der besagte vorprogrammierte Optimierungsalgorithmus ausgewählt aus der Gruppe bestehend aus simulierter Abkühlung („simulated annealing”) Branch-and-Bound-Verfahren („Verzweigung-und-Schranke-Verfahren”) und Monte-Carlo-Auswahlverfahren („Monte-Carlo-Sampling-Methods”).
-
Es ist ein weiteres Ziel der vorliegenden Erfindung, das Verfahren wie in irgendeinem des Obigen definiert zu offenbaren, wobei das Gewebe A und das Gewebe B Gewebearten sind ausgewählt aus der Gruppe bestehend aus grauer Substanz, weißer Substanz und Rückenmarksflüssigkeit.
-
Es ist ein weiteres Ziel der vorliegenden Erfindung, das Verfahren wie in irgendeinem des Obigen definiert zu offenbaren, wobei das Gewebe A ein Tumorgewebe innerhalb eines Organs von Interesse ist und das Gewebe B ein normales Gewebe ist. In einigen Ausführungsformen der Erfindung befinden sich das Gewebe A und das Gewebe B innerhalb eines speziellen Organs.
-
Es ist ein weiteres Ziel der vorliegenden Erfindung, das Verfahren wie in irgendeinem des Obigen definiert zu offenbaren, wobei die Gewebe A und B zwei verschiedene Organe innerhalb eines Gesichtsfeldes des besagten MRT-Geräts sind.
-
Es ist ein weiteres Ziel der vorliegenden Erfindung, das Verfahren wie in irgendeinem des Obigen definiert zu offenbaren, wobei der besagte Schritt des Bestimmens von optimalen Werten von TR und TE umfasst: Verändern von TR und TE innerhalb von Bereichen, die typisch sind für einen Scantyp ausgewählt aus der Gruppe bestehend aus Bereichen typisch für einen T1-gewichteten Scan und Bereichen typisch für einen T2-gewichteten Scan; und Bestimmen des besagten Maximalwertes von CNRAB(TR, TE); wobei das besagte Verfahren automatisch bestimmt, ob ein T1-gewichteter Scan oder ein T2-gewichteter Scan das besagte maximale CNRAB liefert. In einigen Ausführungsformen der Erfindung umfasst der besagte Schritt des Bestimmens von optimalen Werten von TR und TE ein Verändern von TR und TE über zumindest eine Gruppe von Bereichen, die durch eine Gruppe von Randbedingungen begrenzt sind ausgewählt aus der Gruppe bestehend aus TR < 0,75 s, TE < 40 ms und TR > 2 s, TE < 100 ms.
-
Es ist ein weiteres Ziel der vorliegenden Erfindung, das Verfahren wie in irgendeinem des Obigen definiert zu offenbaren, umfassend ein Bestimmen von n zusätzlichen Akquisitionsparametern Pn, n ≥ 1. In einigen Ausführungsformen der Erfindung sind die besagten zusätzlichen Akquisitionsparameter ausgewählt aus der Gruppe bestehend aus Flipwinkel, HF-Pulslänge und HF-Pulsamplitude.
-
Es ist ein weiteres Ziel der vorliegenden Erfindung, das Verfahren wie in irgendeinem des Obigen definiert zu offenbaren, wobei der besagte Schritt des Bestimmens von n zusätzlichen Akquisitionsparametern Pn umfasst: Verändern von jedem der besagten Parameter Pn über einen vorbestimmten Bereich; Bestimmen der besagten optimalen TR und TE für jeden Wert von Pn; und Bestimmen des besagten optimalen Werts von Pn als ein Wert von Pn, der den besagten Maximalwert von CNRAB(TR, TE) ergibt.
-
Es ist ein weiteres Ziel der vorliegenden Erfindung, das Verfahren wie in irgendeinem des Obigen definiert zu offenbaren, wobei der besagte Schritt des Bestimmens von optimalen Werten von n zusätzlichen Akquisitionsparametern Pn ein Bestimmen der optimalen Werte von TR, TE und P1...Pn umfasst, die einen Maximalwert von CNRAB(TR, TE, P1...Pn) ergeben. In einigen Ausführungsformen der Erfindung umfasst der besagte Schritt des Bestimmens der optimalen Werte von TR, TE und P1...Pn, die einen Maximalwert von CNRAB(TR, TE, P1...Pn) ergeben, die Verwendung eines Optimierungsalgorithmus, der das maximale CNRAB über den Funktionsraum von TR, TE, P1...Pn findet. In einigen Ausführungsformen der Erfindung ist der besagte Optimierungsalgorithmus ausgewählt aus der Gruppe bestehend aus simulierter Abkühlung, Branch-and-Bound-Verfahren und Monte-Carlo-Auswahlverfahren.
-
Es ist ein weiteres Ziel der vorliegenden Erfindung, das Verfahren wie in irgendeinem des Obigen definiert zu offenbaren, wobei das besagte Verfahren implementiert ist als Teil der Steuerungs- und/oder Akquisitionssoftware eines MRT-Systems. In einigen bevorzugten Ausführungsformen der Erfindung ist das besagte MRT-System ein MRT-System mit einem Permanentmagneten.
-
Es ist ein weiteres Ziel der vorliegenden Erfindung, ein MRT-System umfassend ein Steuerungs- und/oder Akquisitionsuntersystem zu offenbaren, das dazu programmiert ist, das Verfahren auszuführen, das in irgendeinem des Obigen definiert ist. In einigen bevorzugten Ausführungsformen der Erfindung ist das besagte MRT-System ein MRT-System mit einem Permanentmagneten.
-
Es ist ein weiteres Ziel der vorliegenden Erfindung, ein MRT-System wie in irgendeinem des Obigen definiert zu offenbaren, wobei das besagte Steuerungs- und/oder Akquisitionsuntersystem dazu programmiert ist, das Verfahren wie in irgendeinem des Obigen definiert ohne Eingriff eines Bedieners des besagten Systems auszuführen.
-
Kurze Beschreibung der Figuren
-
Die Erfindung wird nun mit Bezug zu den Figuren beschrieben werden, wobei:
-
1 ein Flussdiagramm der Schritte des hierin offenbarten Verfahrens zum Optimieren des Kontrast-zu-Rauschen-Verhältnisses eines MRT-Bildes als Funktion von TR und TE in einer Ausführungsform zeigt; und
-
2 ein Flussdiagramm der Schritte des hierin offenbarten Verfahrens zum Optimieren des Kontrast-zu-Rauschen-Verhältnisses eines MRT-Bildes als Funktion von TR, TE und zumindest einem anderen Akquisitionsparameter in einer Ausführungsform zeigt.
-
Detaillierte Beschreibung der bevorzugten Ausführungsformen
-
In der folgenden Beschreibung werden verschiedene Aspekte der Erfindung beschrieben werden. Zum Zwecke der Erläuterung werden spezifische Details dargelegt, um ein gründliches Verständnis der Erfindung bereitzustellen. Es wird für den Fachmann ersichtlich sein, dass es andere Ausführungsformen der Erfindung gibt, die in Details abweichen, ohne das Wesen der Erfindung zu betreffen. Daher ist die Erfindung nicht auf das beschränkt, was in den Figuren illustriert und in der Beschreibung beschrieben ist, sondern nur auf das, was in den begleitenden Ansprüchen angegeben ist, wobei der angemessene Umfang der besagten Ansprüche nur durch deren breiteste Interpretation gegeben ist.
-
Wie hierin verwendet beziehen sich die Indices ”A” und ”B” auf einen Parameter, der für die betreffende Substanz spezifisch ist. Beispielsweise ist T1A die T1 der Substanz A, während T1B die T1 der Substanz B ist. Abgesehen vom Zusatz der Indices „A” und „B” werden hierin alle Symbole und Abkürzungen entsprechend der Standard-MRT/NMR-Praxis verwendet.
-
Es ist aus dem Stand der Technik bekannt (siehe
Ting, Y.-L. und Bendel, P., J. Magn. Reson. Imaging 1992, 3, 393–399), dass das CNR zwischen zwei verschiedenen Geweben ausgedrückt werden kann als das Verhältnis zwischen der Differenz in den Signalintensitäten und dem Bildrauschen. Die Signalintensität von Gewebe A (bezeichnet als S
A) ist gegeben durch Gleichung (1):
SA = ρAE1AE2A (1) wobei ρ
A die Dichte der untersuchten NMR-aktiven Kerne (im Allgemeinen werden dies Protonen sein) von Gewebe A ist, und E
1A und E
2A sind gegeben durch Gleichung (2a) bzw. (2b):
-
Diese Gleichungen nehmen an, dass TE << TR, T1.
-
Das CNR zwischen Gewebe A und einem zweiten Gewebe B (CNR
AB) ist gegeben durch Gleichung (3):
wobei P eine Proportionalitätskonstante ist, die auf 1 gesetzt werden kann, und wobei die Division durch die Wurzel von T
R den Ausdruck auf eine konstante Bildgebungszeit normalisiert.
-
In dem hierin offenbarten Verfahren können T1, T2 und ρ experimentell für jeden Gewebetyp aus einem vorläufigen Scan bestimmt werden, oder sie können aus einer Datenbank mit vorher gemessenen Werten importiert werden. Sobald T1, T2 und ρ für jeden Gewebetyp bekannt sind, werden dann TR und TE zur Bildgebung von Gewebe A in Anwesenheit von Gewebe B durch Maximieren von CNRAB als Funktion von TR und TE optimiert. Jeder im Stand der Technik bekannte Optimierungsalgorithmus kann verwendet werden. Als nicht beschränkendes Beispiel kann eine Art Holzhammerverfahren („brute-force”-Herangehensweise) verwendet werden, bei dem eine Abfolge von vorläufigen Scans aufgenommen wird, in denen TR und TE systematisch und unabhängig voneinander über vorbestimmte Bereich von vernünftigen Werten (z. B. 10–100 ms für TE und 0,5–5 s für TR) variiert werden und das Paar von TR- und TE-Werten, das die höchsten CNRAB-Werte ergibt, verwendet wird. Andere Optimierungsverfahren, die TR und TE in vorläufigen Scans variieren, um das Maximum von CNRAB zu finden, wie beispielsweise eine simulierte Abkühlung („simulated Annealing”), Branch-and-Bound-Verfahren und Monte-Carlo-Auswahlverfahren können ebenfalls verwendet werden. Es wird nun auf die 1 Bezug genommen, die ein Flussdiagramm zeigt, das die Schritte des Verfahrens skizziert.
-
Das hierin beschriebene Verfahren kann verwendet werden, um die optimalen TR und TE für die Detektion von jedwedem Gewebetyp A in Anwesenheit von Gewebetyp B zu finden. Nicht beschränkende Beispiele schließen die Maximierung des CNR für graue Substanz gegenüber weißer Substanz (oder umgekehrt) oder von grauer Substanz oder weißer Substanz gegenüber Rückenmarksflüssigkeit in einem Hirnscan ein, die Maximierung des CNR für Tumorgewebe gegenüber normalem Gewebe (in einigen Ausführungsformen sind beide Gewebetypen in einem einzigen Organ von Interesse lokalisiert), die Maximierung des CNR für ein Organ gegenüber einem anderen Organ, wenn sich beide innerhalb des Gesichtsfeldes des MRT-Geräts befinden, usw.
-
Es ist in dem Bereich der Erfindung wobei das hierin offenbarte Verfahren verwendet wird, um automatisch zu bestimmen, ob ein T1-gewichteter Scan oder ein T2-gewichteter Scan den optimalen Kontrast zwischen den beiden Gewebetypen bereitstellt. In dieser Ausführungsform der Erfindung werden TR und TE innerhalb von Bereichen optimiert, die typisch sind entweder für einen T1-gewichteten Scan (z. B. TR < 0,75 s, TE < 40 ms) oder für einen T2-gewichteten Scan (z. B. TR > 2 s, TE < 100 ms), und es wird dann bestimmt, welcher Satz TR oder TE das maximale CNRAB bereitstellt.
-
Es ist auch innerhalb des Bereiches der Erfindung wobei das hierin offenbarte Verfahren verwendet wird, um automatisch andere Akquisitionsparameter P zu optimieren einschließlich, aber nicht beschränkt auf Flipwinkel, HF-Pulslänge und HF-Pulsamplitude. In einer exemplarischen und nicht beschränkenden Ausführungsform wird für jeden zu optimierenden Parameter P eine Schleife zu dem Optimierungsalgorithmus hinzugefügt, in der der betreffende Parameter innerhalb vorbestimmter Grenzen variiert wird, und die optimalen TR und TE werden wie oben beschrieben gefunden. Für jeden Wert von Pn wird das maximale CNRAB(TR, TE) aufgezeichnet. Der Wert von Pn, der zu dem maximalen Wert von CNRAB(TR, TE) führt, wird bei der Akquisition des MRT-Bildes verwendet. Es wird nun auf 2 Bezug genommen, die ein Flussdiagramm zeigt, das diese Ausführungsform des Verfahrens illustriert. In anderen Ausführungsformen des Verfahrens wird CNRAB als seine Funktion von (TR, TE, P1...Pn) behandelt und es wird ein Optimierungsalgorithmus verwendet, der das maximale CNRAB über den Funktionsraum von allen Akquisitionsparametern von Interesse findet. Wie im Fall der Optimierung von CNRAB(TR, TE), was oben offenbart wurde, kann jeder im Stand der Technik bekannte Optimierungsalgorithmus (z. B. simulierte Abkühlung, Branch-and-Bound-Verfahren und Monte-Carlo-Auswahlverfahren) verwendet werden, um den maximalen Wert von CNRAB(TR, TE, P1...Pn) zu finden.
-
In bevorzugten Ausführungsformen des Verfahrens ist es als Teil der Steuer- und Akquisitionssoftware des MRT-Systems implementiert. In anderen Ausführungsformen ist es als eigenständiges Softwarepaket implementiert. Das heißt, dass der Optimierungsalgorithmus automatisch durch das MRT-System ohne irgendeinen Eingriff durch den Bediener des Systems durchgeführt wird.
-
Es ist innerhalb des Bereiches der Erfindung, ein MRT-System zu offenbaren, bei dem das Akquisitionssystem dazu programmiert ist, das hierin offenbarte Optimierungsverfahren als Teil der Bildakquisitionssoftware auszuführen. In bevorzugten Ausführungsformen der Erfindung ist das MRT-Gerät, das ein Steuersystem umfasst, das dazu programmiert ist, das hierin offenbarte Verfahren auszuführen, ein MRT-System mit einem Permanentmagneten.
-
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 Patentliteratur
-
- US 4694250 A [0003]
- US 6781375 B [0004]
- US 7715899 B [0005]
- US 2007/0276221 A [0006]
-
Zitierte Nicht-Patentliteratur
-
- Ting, Y.-L. und Bendel, P., J. Magn. Reson. Imaging 1992, 3, 393–399 [0030]