Pulmonale strukturelle MRT mit freiatmender, selbstgesteuerter Ultrakurz-Echozeit-Bildgebung

Zusammenfassung

Es wird ein Protokoll zur Erstellung hochauflösender struktureller Bilder der Lunge mit Hilfe der Ultrakurzechozeit (UTE) Magnetresonanztomographie (MRT) beschrieben. Dieses Protokoll ermöglicht die Aufnahme von Bildern mit einer einfachen MRT-Pulssequenz während der freien Atmung.

Zusammenfassung

Eine qualitativ hochwertige MRT der Lunge wird durch eine geringe Gewebedichte, eine schnelle MRT-Signalrelaxation sowie Atem- und Herzbewegungen erschwert. Aus diesen Gründen wird die strukturelle Bildgebung der Lunge fast ausschließlich mittels Computertomographie (CT) durchgeführt. Die CT-Bildgebung liefert jedoch ionisierende Strahlung und ist daher für bestimmte gefährdete Bevölkerungsgruppen (z. B. Kinder) oder für Forschungsanwendungen weniger gut geeignet. Als Alternative stößt die MRT mit ultrakurzen Echozeiten (UTE) auf Interesse. Diese Technik kann während der freien Atmung im Verlauf eines ~5-10-minütigen Scans durchgeführt werden. Informationen über Atembewegungen werden zusammen mit Bildern kodiert; Diese Informationen können verwendet werden, um Bilder selbst zu steuern. Durch das Self-Gating entfällt somit die Notwendigkeit einer fortschrittlichen MRT-Pulssequenzprogrammierung oder der Verwendung eines Atembalgs, was die Bildaufnahme vereinfacht. In diesem Protokoll werden einfache, robuste und recheneffiziente Akquisitions- und Rekonstruktionsmethoden für die Aufnahme einer qualitativ hochwertigen UTE-MRT der Lunge vorgestellt. Dieses Protokoll wurde für die Verwendung auf einem 3T-MRT-Scanner entwickelt, aber die gleichen Prinzipien können bei geringerer Magnetfeldstärke implementiert werden. Das Protokoll enthält empfohlene Parametereinstellungen für die radiale 3D-UTE-Bilderfassung sowie Anweisungen für die selbstgesteuerte Bildrekonstruktion, um Bilder in verschiedenen Atemphasen zu erzeugen. Durch die Implementierung dieses Protokolls können Benutzer hochauflösende UTE-Bilder der Lunge mit minimalen bis minimalen bis gar keinen Bewegungsartefakten erstellen. Diese Bilder können zur Beurteilung der Lungenstruktur verwendet werden, die für Forschungszwecke bei einer Vielzahl von Lungenerkrankungen eingesetzt werden kann.

Einleitung

Die hochauflösende Bildgebung der Lungenstruktur ist bei vielen Lungenerkrankungen ein wesentlicher Bestandteil der Diagnostik. In der Regel erfolgt dies mit Hilfe der Computertomographie (CT), die sich hervorragend eignet, um hochauflösende Bilder der Lunge^{zu erstellen 1}. Die CT-Bildgebung liefert jedoch eine nicht unerhebliche Dosis ionisierender Strahlung, was sie für die regelmäßige Wiederholungsbildgebung, die Bildgebung in mehreren verschiedenen Atemphasen oder die Bildgebung bestimmter Bevölkerungsgruppen (z. B. Kinder) ungeeignet macht. Die Magnetresonanztomographie (MRT) birgt nicht das gleiche Risiko ionisierender Strahlung und ist daher für solche bildgebenden Aufgaben zugänglich. Es ist jedoch schwierig, die Lunge mittels MRT abzubilden, da die Gewebedichte gering ist, die Atem- und Herzbewegung sowie die sehr schnelle Signalrelaxation ^2,3,4.

Eine MRT-Technik, die in der Lage ist, diese Herausforderungen zu mildern, ist die MRT ^4,5,6 mit ultrakurzer Echozeit (UTE). Bei der UTE-MRT wird das MRT-Signal unmittelbar nach der Signalanregung abgetastet, wodurch die Auswirkungen einer schnellen Signalrelaxation reduziert werden. Darüber hinaus wird bei dieser Technik der k-Raum von der Mitte nach außen abgetastet, was zu einer signifikanten Überabtastung im Zentrum des k-Raums führt. Dieses Oversampling in der Mitte des k-Raums macht dieses bildgebende Verfahren robust gegenüber Bewegungen. Zusätzlich zu dieser inhärenten Robustheit gegenüber Bewegung kodiert die wiederholte Abtastung des Zentrums des k-Raums Informationen über die Atembewegung, was das Selbstgating von Bildern ermöglicht ^7,8,9. Dieses Self-Gating kann verwendet werden, um Bilder in einer Vielzahl von Atemphasen zu erzeugen. Da der Mensch den größten Teil der Atemphase bei der Exspiration verbringt, ist es üblich, ein Bild für die Endexspiration zu erstellen, da in dieser Phase die meisten Bilddaten erfasst werden.

Es gibt eine Vielzahl von Strategien für das respiratorische Self-Gating in der Lungen-MRT. Die erste Unterscheidung, die getroffen werden muss, ist bildbasiert vs. k-Space-basiertes Gating¹⁰ (Abbildung 1). Beim bildbasierten Gating wird ein Satz von Bildern mit hoher zeitlicher Auflösung erzeugt, indem kleine zeitliche Teilmengen der Bilddaten rekonstruiert werden. Anschließend wird die Position des Zwerchfells in diesen Bildern verwendet, um die Atmungsphase für einen gegebenen Satz von Bildprojektionen^{zu identifizieren 10, 11}. Beim k-Raum-basierten Gating werden Daten aus dem Zentrum des k-Raums ("k₀") untersucht ^8,9,12. Die Signalintensität des Bildes ist in k₀ kodiert, und somit variiert die Intensität des k_0-Punktes mit der Atmung. Projektionen können so basierend auf der Intensität von k₀ in verschiedene Atmungsphasen eingeteilt werden. Sowohl beim bildbasierten als auch beim k-Raum-basierten Gating werden Projektionen mit gleichrespiratorischen Phasen für die Bildrekonstruktion gruppiert. Es wurde vermutet, dass das bildbasierte Gating eine verbesserte Genauigkeit bei der Abschätzung der Atmungsphase bietet und dadurch Bilder mit reduzierter Unschärfe liefert^10,13.

Abbildung 1: Bildbasierte und k-Raum-basierte Self-Gating-Techniken. (A) Beim bildbasierten Gating werden Bilder mit geringer räumlicher Auflösung und hoher zeitlicher Auflösung, die das Zwerchfell zeigen, aus zeitlichen Teilmengen der Gesamtdaten erzeugt. Mit Hilfe einer Linie über dem Zwerchfell kann die Atembewegung visualisiert und für die Bildrekonstruktion in einen Bin-Bereich aufgenommen werden. (B) Beim k-Raum-basierten Gating wird der erste Punkt auf einer zentrierten k-Raumprojektion ("k₀") verwendet, um die Atembewegung zu visualisieren. Nach der Glättung von k₀ sind Signalintensitätsunterschiede basierend auf dem Atemzyklus deutlich sichtbar und können zur Identifizierung verschiedener Atmungsphasen verwendet werden. Bitte klicken Sie hier, um eine größere Version dieser Abbildung anzuzeigen.

Sowohl bild- als auch k-raum-basiertes Gating können entweder mit hartem Gating oder weichem Gating durchgeführt werden^11,14. Beim Hard Gating werden nur die Vorsprünge rekonstruiert, die der gewünschten Atmungsphase entsprechen. Diese Verwerfung unerwünschter Projektionen kann jedoch zu einem verringerten Signal-Rausch-Verhältnis (SNR) von Bildern und zu vermehrten Undersampling-Artefakten führen. Diese unerwünschten Effekte können durch die Verwendung von Soft Gating gemildert werden. Beim Soft-Gating werden alle Projektionen für die Bildrekonstruktion verwendet, Projektionen aus einer unerwünschten Atmungsphase werden jedoch so gewichtet, dass sie einen geringeren Einfluss auf das endgültige Bild haben. Auf diese Weise können Bilder mit minimalen Artefakten und hohem SNR rekonstruiert werden, während die Auswirkungen der Atembewegung weiterhin unterdrückt werden.

Durch die Kombination der UTE-MRT-Aufnahme mit dem Self-Gating nach der Aufnahme können qualitativ hochwertige Bilder erzeugt werden, die zwar nicht der CT entsprechen, aber einen Kontrast und eine Auflösung aufweisen, die sich der der CT-Bildgebung 6,15,16,17,18,19 annähern. Darin wird ein einfaches Protokoll zum Sammeln und Rekonstruieren von UTE-MRT-Bildern bereitgestellt, um qualitativ hochwertige Bilder der Lungenstruktur zu erzeugen.

Dieses Protokoll wurde hauptsächlich für 3T-MRT-Scanner geschrieben. 3T ist die gebräuchlichste Feldstärke, die für die Forschungs-MRT verwendet wird. Niedrigere Magnetfeldstärken wie 1,5T oder das kürzlich verfügbare 0,55 T²⁰ können für eine verbesserte Bildqualität und Signalintensität in der Lunge sorgen, da die Signalrelaxation in der Lunge bei diesen Feldstärken langsamer ist.

Obwohl alle Anstrengungen unternommen wurden, um Klarheit und Einfachheit in diesem Protokoll und dem bereitgestellten Bildrekonstruktionscode zu schaffen, wird das Protokoll wahrscheinlich einen dedizierten MRT-Physiker (oder einen ähnlichen MRT-Experten) erfordern, um eine geeignete UTE-MRT-Sequenz auf dem MRT-Scanner zu erstellen. Die MRT-Sequenz sollte eine nicht-kartesische 3D-Kodierungsstrategie mit zentrierten k-Raum-Trajektorien implementieren. Beispiele hierfür sind radiale 3D- oder 3D-Spiralsequenzen (z. B. "FLORET")^21,22 Bildsequenzen. Wichtig ist, dass die Reihenfolge der Projektionen eine gute zeitliche Stabilität aufweist: Über eine gegebene Teilmenge der Zeit sollten die Projektionen den gesamten Bereich des k-Raums²³ abdecken. Beispiele für Projektionsordnungsstrategien mit guter zeitlicher Stabilität sind die goldenen Mittel oder die Halton-randomisierte archimedische Spirale. Wenn eine Projektionsordnung mit geringer zeitlicher Stabilität verwendet wird, werden beim Self-Gating nach der Erfassung große Bereiche des k-Raums weggelassen, was zu Bildartefakten führt. Schließlich sollte die Sequenz in der Lage sein, eine Echozeit (TE) von <100 μs zu erreichen. Die T₂*-Relaxationszeit in der Lunge bei 3T beträgt <1 ms²⁴, daher ist die Verwendung einer sehr kurzen TE für die Erzeugung qualitativ hochwertiger Bilder unerlässlich.

Protokoll

Alle Bildgebungen am Menschen wurden mit Genehmigung des KUMC IRB durchgeführt. Von allen Teilnehmern wurde eine schriftliche Einverständniserklärung eingeholt. Die Bilder in dieser Studie wurden im Rahmen eines generischen technischen Entwicklungsprotokolls aufgenommen, und die Ein-/Ausschlusskriterien waren bewusst weit gefasst. Einschlusskriterien: Alter ≥ 18 Jahre. Ausschlusskriterien: MRT kontraindiziert auf der Grundlage der Antworten auf den MRT-Screening-Fragebogen und Schwangerschaft. Das Zubehör und die Ausrüstung, die für diese Studie verwendet werden, sind in der Materialtabelle aufgeführt.

1. UTE-Bilderfassung

1. Bereiten Sie die Bildgebungssequenz vor. Bereiten Sie die Bildgebungssequenz einmal vor und verwenden Sie dieselbe Sequenz für alle Teilnehmer.
   1. Stellen Sie die Parameter gemäß Tabelle 1 ein.
   2. Platzieren Sie ein MRT-Phantom in der Mitte des MRT und führen Sie die Bildgebungssequenz aus.
      HINWEIS: Da diese Sequenz eine schnelle Gradientenleistung und viele HF-Impulse erfordert, ist es wichtig zu überprüfen, ob die Protokolleinrichtung vor dem Test an einem Menschen ausgeführt werden kann.
2. Bereiten Sie den Teilnehmer auf die MRT vor. Verwenden Sie ein MRT-Sicherheitsscreening nach institutionellem Standard, um sicherzustellen, dass der Teilnehmer sicher in die MRT eintreten kann.
3. Positionieren Sie den Teilnehmer auf der MRT-Liege und legen Sie eine Brustspirale über den Oberkörper des Teilnehmers. Positionieren Sie die Spule nahe am Kinn des Teilnehmers, um eine vollständige Abdeckung der Lungenspitzen zu gewährleisten.
4. Bewegen Sie den Teilnehmer in den MRT-Scanner. Platzieren Sie den Positionierungsorientierungspunkt direkt unter dem Brustbein des Teilnehmers.
5. Führen Sie einen Lokalisatorscan durch, um sicherzustellen, dass sich die Lungen des Teilnehmers innerhalb des Sichtfelds für den UTE-Scan befinden. Verschieben Sie nicht die Geometrie des UTE-Scans. Wenn sich die Lungen des Teilnehmers nicht im Sichtfeld befinden, bewegen Sie den Teilnehmer und nehmen Sie zusätzliche Lokalisationsscans auf, bis sich die Lungen vollständig im Sichtfeld befinden.
6. Führen Sie die UTE-Sequenz aus. Während dieser Sequenz kann der Teilnehmer normal atmen.
7. Exportieren Sie die Rohdaten aus dem Scanner. Abhängig von der verwendeten Bildsequenz kann der Scanner Bilder auf dem Scanner rekonstruieren oder nicht. Für die vorgeschlagene retrospektive Gating-Rekonstruktion werden Rohdaten benötigt, um zu bestimmen, ob Bilder auf dem Scanner erzeugt werden oder nicht. Beachten Sie, dass die Rohdaten groß sein werden (>10 GB).
8. Exportieren oder berechnen Sie k-Raum-Trajektorien (d. h. die Position jedes Rohdatenpunkts im k-Raum).
   HINWEIS: Bei einigen Bildgebungssequenzen können K-Raum-Trajektorien zusammen mit den Rohdaten auf dem MRT-Scanner gespeichert und direkt exportiert werden. Für andere Bildgebungssequenzen müssen die k-Raum-Trajektorien auf der Grundlage von Bildgebungsparametern berechnet werden.

Parameter	Generische empfohlene Einstellungen	Hierin implementierte Einstellungen
Bildgebende Sequenz	Nicht-kartesische 3D-Trajektorien mit zentrierten k-Raum-Trajektorien	3D-Radial mit Goldener Mittel-Projektion Reihenfolge
Sichtfeld	400 x 400 x 400 mm3	400 x 400 x 400 mm3
Matrix-Größe	Wie gewünscht für die Zielauflösung	320 x 320 x 320 (1,25 mm isotrope Auflösung)
Bandbreite	Bei Bedarf für die Auslesedauer < 1,0 ms	888 Hz/Pixel
TE	< 0,1 ms	0,07 ms
TR	Minimum (Ziel 3 – 4 ms)	3,5 ms
Flip-Winkel	Ungefähr 5°	4.8°
Anzahl der Projektionen	Mindestens 100.000	1,35,386
Dauer des Bildes	Mindestens 5 Minuten	7 Minuten, 54 Sek.

Tabelle 1: Empfohlene Einstellungen für die UTE-Bildgebung. Es werden generische empfohlene Einstellungen bereitgestellt, die zur Einrichtung des Protokolls verwendet werden können. Spezifische empfohlene Einstellungen, die für die Daten verwendet wurden, werden ebenfalls bereitgestellt, wie als repräsentative Ergebnisse dargestellt. Die Parameterspezifikationen sind bei allen Anbietern generisch, mit Ausnahme der Bandbreite. Einige große MRT-Anbieter geben die Bandbreite als Hz/Pixel an. Andere große MRT-Anbieter geben absolute Bandbreite an. Die empfohlene Bandbreite (888 Hz/Pixel) entspricht einer absoluten Bandbreite von 284.160 Hz.

2. UTE-Bildrekonstruktion mittels bildbasiertem respiratorischem Soft-Gating

HINWEIS: MATLAB-Code zum Ausführen der folgenden Schritte finden Sie unter https://github.com/pniedbalski3/UTE_Reconstruction.

1. Importieren Sie Daten und k-Raum-Trajektorien in MATLAB. Der Code zum Importieren von MRT-Rohdaten ist für alle großen MRT-Anbieter verfügbar.
2. Verwerfen Sie die ersten 1000 Projektionen, um sicherzustellen, dass sich die Daten in stationärer Magnetisierung befinden.
   HINWEIS: Wenn die verwendete Bildgebungssequenz Dummy-Scans vor der Datenerfassung umfasst, kann dieser Schritt übersprungen werden.
3. Rekonstruieren Sie ein Bild mit niedriger Auflösung mit einer sehr kleinen Teilmenge von Daten.
   1. Rekonstruieren Sie das Bild mit einer ungleichmäßigen schnellen Fourier-Transformation auf eine Matrixgröße von 96 x 96 x 96.
   2. Verwenden Sie ca. 200 Projektionen, was Daten von 0,6 s bis 0,8 s entspricht.
   3. Rekonstruieren und speichern Sie Bilder von allen Spulenelementen sowie ein endgültiges, spiralkombiniertes Bild.
4. Wählen Sie im resultierenden spiralkombinierten Bild eine koronale Schicht aus, die das Zwerchfell deutlich zeigt.
   HINWEIS: Der bereitgestellte Code fordert den Benutzer auf, eine Scheibe auszuwählen, die die Membran enthält.
5. Nachdem Sie diese Scheibe ausgewählt haben, sehen Sie sich die einzelnen Spulenbilder für diese Scheibe an und wählen Sie ein oder zwei Spulenelemente aus, die die Membran am besten darstellen (Abbildung 2).
   HINWEIS: Der bereitgestellte Code fordert den Benutzer auf, Spulenelemente auszuwählen.
6. Rekonstruieren Sie Bilder mit einem gleitenden Fenster , um Bilder mit einer zeitlichen Auflösung von ~0,5 s zu generieren (Abbildung 2).
   1. Rekonstruieren Sie nur die Daten aus den in Schritt 2.4 ausgewählten Spulenelementen.
      HINWEIS: Obwohl alle Spulenelemente rekonstruiert werden können, werden nur die Elemente benötigt, die der Membran am nächsten sind, um die Membran zum Zwecke der Atemwegsselbstkontrolle zu visualisieren. Indem nur die Spulenelemente rekonstruiert werden, die dem Diagramm am nächsten sind, werden die Rekonstruktionszeit und der Rechenaufwand drastisch reduziert.
   2. Verwenden Sie die ersten 200 Projektionen, um ein Bild mit einer ungleichmäßigen schnellen Fourier-Transformation zu rekonstruieren (Abbildung 2). Bewahren Sie nur die Scheibe auf, die die Membran zeigt (wie in Schritt 2.4 beschrieben).
      HINWEIS: Letztendlich werden bis zu 1500 Bilder generiert; Es wird nur eine 2D-Schicht benötigt, um die Membranposition zu visualisieren, und das Speichern von 3D-Bildern für jede der Schiebefensterstufen wäre unerschwinglich.
   3. Verschiebung um 100 Projektionen (d.h. das erste Bild wird mit den Projektionen 1-200 rekonstruiert. Die zweite wird mit den Projektionen 101 - 300) rekonstruiert und rekonstruiert ein zusätzliches Bild, wobei die in Schritt 2.4 ausgewählte Schicht gespeichert wird.
   4. Fahren Sie fort, bis alle Projektionen zum Generieren von Bildern verwendet wurden.
7. Wählen Sie im ersten Bild des Schiebefensters eine Linie über der Blende aus. Stellen Sie sicher, dass die Leitung lang genug ist, um sich um 5-10 Voxel in die Lunge und um 5-10 Voxel in das Zwerchfell zu erstrecken.
8. Visualisieren Sie die Atembewegung, indem Sie diesen Atemnavigator für alle Projektionen anzeigen.
9. Bestimmen Sie die Position des Zwerchfells für alle Atemnavigationsgeräte. Es gibt eine Vielzahl von Möglichkeiten, dies zu tun, aber eine einfache Methode besteht darin, die Methode²⁵ von Otsu zu verwenden, um die dunklere Seite (Lunge) von der helleren Seite (Zwerchfell) zu trennen.
10. Verwenden Sie die Position des Zwerchfells, um Projektionen als zu einem bestimmten Atemwegsbehälter gehörend zu kennzeichnen. Wenn ein gegebener Atemnavigator das Zwerchfell an "Position 1" anzeigt, dann gehören die 200 Projektionen, die zur Erzeugung des Bildes für diesen Navigator verwendet wurden, zu "Bin 1".
    HINWEIS: Da Bilder mit einem Schiebefenster mit einer Überlappung von 100 Projektionen erstellt wurden, können einige Projektionen als zu mehreren Klassen gehörend gekennzeichnet sein. Die grobe räumliche Auflösung von Schiebefensterbildern führt zu insgesamt ~4-6 Bins, die die gesamte Bandbreite der Inspiration bis zum Ablauf abdecken.
11. Wählen Sie den Abschnitt aus, der rekonstruiert werden soll, indem Sie bestimmen, welcher Abschnitt die größte Anzahl von Projektionen aufweist, die dem Endablauf entsprechen sollten.
    1. Alternativ können Sie Bilder für die gewünschten Atemphasen auf der Grundlage einer visuellen Inspektion des Atemnavigators rekonstruieren.
12. Generieren von Gewichten für Soft-Gating¹⁴.
    1. Verwenden Sie einen Exponentialfilter , um den Projektionen innerhalb des primären Behälters eine Gewichtung von 1 und den Vorsprüngen in verschiedenen Atemwegsbehältern eine stark reduzierte Gewichtung zu verleihen.
13. Verwenden Sie die Berkely Advanced Reconstruction Toolbox (BART; https://mrirecon.github.io/bart/)^26,27, um ein hochauflösendes Bild am gewünschten Atemwegsbehälter zu rekonstruieren.
    HINWEIS: BART ist eine frei verfügbare Toolbox für die MRT-Bildrekonstruktion.
    1. Berechnen Sie Dichtekompensationsgewichte mithilfe der iterativen Dichtekombination.
    2. Skalieren Sie die Dichteausgleichsgewichte mit den Soft-Gating-Gewichten.
    3. Skalendaten basierend auf Dichtekompensation und Soft-Gating-Gewichten
    4. Führen Sie eine grundlegende ungleichmäßige schnelle Fourier-Transformation (NUFFT) durch, um die Spulenkombination zu erleichtern.
    5. Konvertieren Sie das NUFFT-Bild in einen gerasterten k-Raum, der für die Spulenkombination verwendet werden soll.
    6. Generieren Sie eine Spulenkombinationsmatrix und verwenden Sie sie, um Spulen sowohl für die Rohdaten als auch für den gerasterten k-Raum zu kombinieren.
    7. Schätzen Sie die Empfindlichkeit der Spule.
    8. Führen Sie eine parallele Imaging-Rekonstruktion des komprimierten Senses mit der gewichteten Dichtekompensation, kombinierten Spulendaten und Spulenempfindlichkeitskarten durch.
14. Speichern Sie das endgültige Bild. Das NIFTI-Format ist einfach zu implementieren. Wenn das Bild auf ein PACs-System hochgeladen werden soll, kann ein DICOM-Format erforderlich sein.

Abbildung 2: Bildbasiertes Self-Gating. (1) Identifizieren Sie anhand eines niedrig aufgelösten Bildes, das aus einer kleinen Anzahl von Projektionen rekonstruiert wurde (zur Recheneffizienz), eine koronale Schnitte, die das Zwerchfell deutlich zeigt. (2) Wählen Sie durch Untersuchen von Bildern von einzelnen Spulenelementen die Spulenelemente aus, die der Membran am nächsten sind. (3) Durchführen einer Schiebefensterrekonstruktion nur der Spulenelemente, die der Membran am nächsten sind (für die rechnerische Effizienz). Bilder können aus Teilmengen von 200 Projektionen (entsprechend ~0,8 s) erzeugt werden; Durch überlappende Projektionen kann eine pseudo-zeitliche Auflösung von ~0,5 s in Bildern erreicht werden. (4) Identifizieren einer Linie, die senkrecht zum Zwerchfell verläuft und als Atemnavigator verwendet werden soll. (5) Die Visualisierung der Bilddaten auf dieser Linie zeigt die Atembewegung, die zur Ablage von Bildern verwendet werden kann. Bitte klicken Sie hier, um eine größere Version dieser Abbildung anzuzeigen.

3. UTE-Bildrekonstruktion mit k-Space-basiertem respiratorischem Soft-Gating

1. Die Schritte 2.1 bis 2.4 sind so durchzuführen, dass das Spulenelement, das der Membran am nächsten liegt, identifiziert werden kann.
2. Generieren Sie eine k 0-Zeitreihenverfolgung, indem Sie den absoluten Wert des ersten Punktes in der Projektion für alle Projektionen für das ausgewählte Windungselement verwenden. Dies ermöglicht die Visualisierung einer Atemwellenform.
3. Normalisieren Sie in Schritten von 5000 Projektionen k₀ durch die mittlere Signalintensität derselben k₀ Punkte²⁸. Dies verringert die Drift der Signalintensität im Laufe der Zeit und bietet eine verbesserte Möglichkeit, Projektionen quantitativ zu klassifizieren.
4. Beschriften Sie jeden k_0-Punkt als während der Inspiration oder des Ablaufs auftretend.
   1. Glätten Sie die k 0-Zeitreihe, und nehmen Sie die Ableitung, um die Steigung für jeden Punkt auf der Gating-Kurve zu bewerten.
   2. Beschriften Sie Inspirationspunkte basierend auf dem Zeichen der Steigung. Eine positive Steigung entspricht der Exspiration, während eine negative Steigung der Inspiration entspricht.
5. Bin-Projektionen basierend auf der Signalintensität. Da die Atemtiefe variabel sein kann, basieren Bin-Projektionen auf der Signalamplitude und nicht auf der Position in der Atemphase.
   HINWEIS: Eine einfache und schnelle Methode, um dies zu erreichen, ist die Implementierung von k-Means-Clustering, um verschiedene Signalintensitätsstufen zu identifizieren.
6. Identifizieren Sie für Bins, die zwischen Endinspiration und Endablauf liegen, Projektionen, die während der Inspiration und des Ablaufs basierend auf Schritt 3.4 auftreten.
7. Schließen Sie die Bildrekonstruktion mit den Schritten in Schritt 2.10 bis Schritt 2.13 ab.
8. Falls gewünscht, rekonstruieren Sie Bilder für alle Atemwegsbehälter und nicht nur am Ende des Ablaufs.

Ergebnisse

Repräsentative Ergebnisse (Abbildung 3) wurden mit den in Tabelle 1 gezeigten Einstellungen generiert. Die verwendete Bildgebungsdauer liefert qualitativ hochwertige Bilder, die von den meisten Teilnehmern toleriert werden können.

Abbildung 3: Repr?...

Diskussion

Bei der UTE-Bildgebung der Lunge können viele Variationen sowohl der Erfassung als auch der Rekonstruktion verwendet werden, um Bilder der Lunge zu erstellen. Dieses Protokoll konzentriert sich auf die einfache Implementierung und die Recheneffizienz. Die Bildgebung mit radialer 3D-UTE ist relativ einfach, da Bildgebungssequenzen in der Regel von den großen MRT-Anbietern erhältlich sind. MATLAB-basierte Tools werden für die Datenverarbeitung und das Self-Gating bereitgestellt. Da die...

Materialien

Name	Company	Catalog Number	Comments
Chest MRI Coil	Siemens, GE, Philips,, Other Clinical MRI Imaging Coil Vendor	N/A	A 26 - 32 channel Chest coil should be used
High Performance Workstation	HP, Apple, or other Computer Hardware company	N/A	A computer with a minimum of 64 GB of Memory is needed for image reconstruction
Matlab	Mathworks	R2016A or newer	A Matlab license is needed to run the provided computer code
MRI Phantom	Siemens, GE, Philips, or Other MRI Phantom Vendor	N/A	Any Phantom can be used to test the MRI sequence prior to its use in human subjects.
MRI Scanner	Siemens, GE, Philips, or Other Clinical MRI Scanner Vendor	N/A	The protocol was developed on a 3T scanner, but 1.5T or 0.55T would also work with minimal adaptation