Jak Funkcjonują Skanery MRI — Magnetyczny Rezonans Nuklearny i Obrazowanie Medyczne

Cząsteczka protonu i chwilowe momenty magnetyczne

Każdy proton (jądro wodoru) ma wewnętrzny moment kątowy nazywany spinem. Niezależnie od tego, że spin jest cechą kwantową bez odpowiednika klasycznego, zachowuje się macroscopycznie jak mały magnes: każdy proton nosi moment magnetyczny μ położony wzdłuż osi swego spinu.

μ = γ · ℏ · I γ (iloraz gyromagnetyczny dla ¹H) = 2,675 × 10⁸ rad s⁻¹ T⁻¹ ℏ = stała Plancka zmniejszona = 1,055 × 10⁻³⁴ J·s I = licznik kwantowy spinu = ½ dla ¹H Ciało ludzkie stanowi około 60% wody (H₂O), co sprawia, że wodor jest daleko najpopularniejszym nukleusem aktywnym w NMR. Zatem MRI kliniczne przede wszystkim obrazuje gęstość protonów i właściwości rozluźniania tkanek.

μ = γ · ℏ · I
γ (gyromagnetic ratio for ¹H) = 2.675 × 10⁸ rad s⁻¹ T⁻¹
ℏ = reduced Planck constant = 1.055 × 10⁻³⁴ J·s
I = spin quantum number = ½ for ¹H

Aliniacja w silnym statycznym polu (B₀)

Zewnętrznie w polu magnetycznym cząsteczki protonowe mają spin, który wskazuje w losowych kierunkach i się nawzajem wywalam. W wnętrzu boju skanera MRI generowane jest silne statyczne pole B₀ (typowo 1,5 T lub 3 T — co do 50,000 razy mocniejsze od pola Ziemi). Mechanika kwantowa pozwala na tylko dwie stanowiska energetyczne dla jądra z spinem ½ w polu B₀, oddzielone różnica energii ΔE = γℏB₀. Na ciele ludzkiego temperatury istnieje mała przewaga Boltzmanna o około 7 spinów na milion w niższej energetycznej (paralelnej) stanowisku, i ta niewielka nieporządkowanie tworzy netto magnetyzację makroskopową M₀ skierowaną w kierunku B₀ — proporcjonalną do siły pola, co jest jednym z powodów, dla których skaner 3 T daje ostrojsze obrazy niż skaner 1,5 T.

Impuls optymalny i rezonans Larmorowski

Aby uzyskać sygnał do detekcji, skaner ustawia netową magnetyzację M₀ pod kątem do osi B₀ za pomocą krótkiego impulsu optymalnego frekwencji radio (RF). Kluczem jest rezonans: częstotliwość impulsu musi dokładnie pasować do częstotliwości Larmorowskiej ω₀.

ω₀ = γ · B₀ (równanie Larmorowskie) Dla ¹H przy 1,5 T: ω₀ ≈ 64 MHz Dla ¹H przy 3,0 T: ω₀ ≈ 128 MHz Dla ¹H przy 7,0 T: ω₀ ≈ 298 MHz (zakres FM radiowych) Impuls o kącie 90° ustawia M₀ całkowicie w płaszczyźnie prostopadłej do B₀. Magnetyzacja następnie precesuje wokół B₀ z częstotliwością ω₀, przekraczając przeciwnik odbioru i spowodowując sinusoidalny napięcie — sygnał Free Induction Decay (FID). Równania Blocha opisują pełny wektor magnetyzacji M(t) pod wpływem kombinacji pola statycznego B₀, pola optymalnego frekwencji radio B₁ oraz relaksacji.

ω₀ = γ · B₀  (Larmor equation)
For ¹H at 1.5 T: ω₀ ≈ 64 MHz
For ¹H at 3.0 T: ω₀ ≈ 128 MHz
For ¹H at 7.0 T: ω₀ ≈ 298 MHz (FM radio band)

Relaksacja T1 i T2 — źródło kontrastu

Po wyłączeniu impulsu RF magnetyzacja powraca do równowagi poprzez dwa niezależne procesy, a stałe czasowe charakterystyczne są głównym źródłem kontrastu tkankowego w MRI. Relaksacja longitudinalna (T1) to ekspansyjny powrót Mz do M₀, gdy spiny przekazują energię do otaczającego lattisa molekularnego: Mz(t) = M₀·[1 − e^(−t/T1)]. Relaksacja transwersalna (T2) to zanikanie Mxy, gdy spiny odchylają się od lokalnych zmian pola magnetycznego: Mxy(t) = Mxy(0)·e^(−t/T2). T2 zawsze jest ≤ T1.

Typowe wartości przy 1.5 T: T1 T2 Biała makula ~780 ms ~90 ms Czarne makula ~920 ms ~100 ms Tłuszcz ~260 ms ~80 ms Woda (CSF) ~2400 ms ~1800 ms Lipid zwracają się w taki sposób, że efektywnie przekazują energię, co daje krótki T1 i jasne obrazy na obrazach ocharakteryzowanych T1. Woda wolna ma długi T1 i T2; objawy choroby często wydzielają się jako nadjasne na obrazach ocharakteryzowanych T2, ponieważ procesy patologiczne zwiększają zawartość wody tkankowej. Wybierając TR (czas powtórki) i TE (czas eco), radiologowie podkreślają charakterystykę T1 (krótki TR/TE — tłuszcz jasny, woda ciemna), T2 (długi TR/TE — woda jasna, patologia widoczna), koncentrację protonów lub FLAIR (T2 z CSF supresjonowanym, używane do wykrywania plakatów MS blisko komor).

Typical values at 1.5 T:
              T1        T2
White matter  ~780 ms   ~90 ms
Grey matter   ~920 ms   ~100 ms
Fat           ~260 ms   ~80 ms
Water (CSF)   ~2400 ms  ~1800 ms

Obwody gradientowe, przestrzeń k i transformata Fouriera

Po impulsie RF wszystkie protony generują sygnał przy tej samej częstotliwości Larmorowej, więc skaner nadal nie może określić źródła tego sygnału. Trzy zestawy obwodów gradientowych (Gx, Gy, Gz) dodają małe, liniowe zmiany do B₀, które kodują położenie przestrzenne w częstotliwości i fazie: obwód Gz wybiera(slice) wybraną warstwę, obwód Gx koduje jedno osiowe półprzecięcie w czasie odczytu przez częstotliwość, a obwód Gy koduje drugie półprzecięcie przez fazę. Raw dane są przechowywane w 2-D macierzy nazywanej przestrzenią k, gdzie każda para [kx, ky] reprezentuje częstotliwość przestrzenną obrazu.

S(kx, ky) = ∫∫ ρ(x,y) · e^(−i2π(kx·x + ky·y)) dx dy → S(kx, ky) jest transformata Fouriera 2-D gęstości protonów ρ(x,y) Rekonstrukcja obrazu: ρ(x,y) = 2D-IFFT{ S(kx, ky) } Środek przestrzeni k zawiera niskie częstotliwości przestrzenne — ogólne jasność i kontrast — a okoliczność zawiera wysokie częstotliwości — krawędzie i szczegółowe szczegóły. Stąd powstają szybsze, ale mniej ostre obrazy przy zapisywaniu tylko części przestrzeni k. Technika Echo Planar Imaging (EPI) wypełnia całą przestrzeń k po jednym impulecie RF, szybko oscylując Gx i Gy w wzorcu zatrzaskiwania, co trwa około 50 ms na warstwę — wystarczająco szybko dla MRI funkcjonalnego i obrazowania serca.

S(kx, ky) = ∫∫ ρ(x,y) · e^(−i2π(kx·x + ky·y)) dx dy
→ S(kx, ky) is the 2-D Fourier transform of proton density ρ(x,y)
Image reconstruction: ρ(x,y) = 2D-IFFT{ S(kx, ky) }

Agenty kontrastowe, fMRI i obrazowanie dyfuzji

Gadolinium (Gd³⁺) jest silnie paramagnetyczny i znacznie skróca T1 sąsiednich protonów; wtrącanie go intraweenousowo spowoduje jego akumulację tam, gdzie zrównoważonej barierki krwiono-mozkowej nie ma, tworząc jasne punkty na obrazach o wysokiej wagowosci T1. W fMRI deoksy-hemoglobina jest paramagnetyczna i skróca T2*, podczas gdy oksy-hemoglobina jest diamagnetyczna i tego nie robi — obszary aktywnego mózgu mają wyższą oksydację krwi, dłuższy T2* oraz jaśniejszy sygnał na sekwencjach z gradientami, czyli BOLD kontrast, który stanowi podstawę funkcjonalnego MRI. Obrazowanie Tensora Dyfuzji (DTI) stosuje silne pociągi gradacyjne w wielu kierunkach, aby pomiar dyfuzji wodnej anizotropicznej w osiadłych białych chrząstek nerwowych, co pozwala odkryć ścieżki neuronów i wykryć wcześniejszy szkodliwy wpływ na axony w ataku krwiopuklowym i MS.

Bezpieczeństwo i siły pola

Niezależnie od CT lub PET, MRI nie korzysta z promieniowania ionizującego. Przepisy bezpieczeństwa ograniczają depozycję energii RF (SAR) do 4 W/kg. Glowne zagrożenia to projekcyjny efekt na magnetyczne obiekty, ciepło w implantach spowodowane prądami eddy, dźwięk akustyczny do 130 dB podczas przełączania zwojów gradientowych oraz claustrofobia (rozwiązana przez projektory otwartych i szerokoborewnych). Siły pola kliniczne wynoszą od 0,5–1,0 T dla scannerów otwartych, przechodzą przez 1,5 T – standard złoty od 1985 roku, aż do systemów o 3 T, które podwójnie zwiększają czułość sygnału na neuro i kardiozawodach. Systemy o 7 T (zatwierdzone przez FDA w 2017 roku) są używane do obrazowania korzy otwartych w skali submilimetrycznej, a magnet Iseult o siłce pola 11,7 T – pierwsze obrazy ludzkie opublikowano w marcu 2024 roku. Szybkie przełączanie zwojów gradientowych jest ograniczone przez podnoszenie potencjału nerwowego (PNS), co wpływa na wydajność EPI i sekwencji rozproszenia w systemach o wysokich siłach pola.

Często zadawane pytania

Wypróbuj na żywo

Wszystko powyżej działa bezpośrednio w Twojej przeglądarce — otwórz MRI / NMR Physics — Bloch Equations i zmieniaj parametry podczas działania. Nic nie jest instalowane ani przesyłane na serwer, cały model działa w jednej karcie.