Przestrzeń k i rekonstrukcja

Obraz, zanim zacznie wyglądać jak obraz

W CT dane wejściowe rekonstrukcji można naturalnie opisać jako projekcje: pomiary wzdłuż dróg przez pacjenta są porządkowane i przekształcane w obraz. MRI podąża inną drogą. Podczas odczytu gradienty pola magnetycznego sprawiają, że faza i częstotliwość odbieranego sygnału zależą od położenia przestrzennego. Skaner próbkuje sygnał złożony, zamiast mierzyć niezależnie każdy piksel obrazu. Próbki te układa się w przestrzeni k, która jest reprezentacją informacji przestrzennej obiektu w dziedzinie częstotliwości.

Dla wystarczająco spróbkowanej akwizycji kartezjańskiej odwrotna transformacja Fouriera przekształca zespolone dane przestrzeni k w obraz. Współrzędne w przestrzeni k nie są położeniami w ciele pacjenta. Opisują częstotliwości przestrzenne, czyli szybkość zmian sygnału w przestrzeni. Dlatego przestrzeń k może wydawać się oderwana od wyświetlanej anatomii, choć stanowi bezpośrednie dane wejściowe rekonstrukcji.

Częstotliwość przestrzenna daje użyteczny model myślowy

Najprostsza intuicja traktuje obraz jako sumę fal. Próbki w pobliżu środka przestrzeni k reprezentują niższe częstotliwości przestrzenne i zwykle odpowiadają za szerokie zmiany intensywności, ogólny kontrast i zgrubny kształt. Próbki położone dalej od środka reprezentują wyższe częstotliwości przestrzenne i odpowiadają za ostre granice oraz drobne szczegóły. Jest to użyteczne pierwsze przybliżenie, a nie ostry podział anatomiczny: struktury są reprezentowane w pewnym zakresie częstotliwości, a dokładny efekt zależy od wzorca próbkowania i rekonstrukcji.

W siatce kartezjańskiej zakres spróbkowanej przestrzeni k, czyli najwyższa osiągnięta częstotliwość przestrzenna, jest związany z nominalną rozdzielczością. Odstęp między próbkami, równoważnie gęstość próbkowania w danym zakresie, jest związany z polem widzenia. Zależności te wynikają z teorii próbkowania, lecz trajektorie niekartezjańskie, ważenie gęstością, częściowe próbkowanie i założenia rekonstrukcji sprawiają, że proste hasło „środek/obrzeża” jest niepełne. Pole widzenia i rozdzielczość są właściwościami całej akwizycji i rekonstrukcji, a nie jednego odizolowanego piksela.

Schemat syntetyczny, a nie akwizycja MRI: w przestrzeni k nie ma uniwersalnej ostrej granicy między środkiem a obrzeżami, a rzeczywiste trajektorie próbkowania i maski niedostatecznego próbkowania są dyskretne i zróżnicowane. Akceleracja (obrazowanie równoległe, compressed sensing, rekonstrukcja DL) obejmuje znacznie więcej niż pominięcie obrzeży przestrzeni k. Żaden przedstawiony element nie jest ilościowym pomiarem MRI.

Rysunek konkretyzuje intuicję. Zachowanie wyłącznie centralnego obszaru niskich częstotliwości może ocalić znaczną część szerokiego kontrastu, jednocześnie rozmywając krawędzie i małe struktury. Zachowanie wyłącznie obszaru zewnętrznego nie tworzy kompletnego obrazu o wysokiej rozdzielczości, ponieważ wkład niskich częstotliwości jest również potrzebny do odtworzenia szerokich zmian sygnału. Znaczenie mają zarówno faza, jak i moduł próbek zespolonych; wyświetlany obraz modułu przestrzeni k jest wizualizacją, a nie danymi, które otrzymuje odwrotna transformacja.

Próbkowanie i rekonstrukcja należą do toru pomiarowego

Akceleracja zmienia sposób próbkowania przestrzeni k, aby skrócić badanie. Obrazowanie równoległe wykorzystuje informacje o czułości wielu cewek odbiorczych do rekonstrukcji niedostatecznie spróbkowanej akwizycji. Compressed sensing narzuca założenia takie jak rzadkość lub strukturalna zmienność. Rekonstrukcja wyuczona wykorzystuje informacje poznane na przykładach, aby odwzorować dane niedostatecznie spróbkowane lub w inny sposób zdegradowane na obraz. Metody te nie są równoważne i żadnej z nich nie należy sprowadzać do stwierdzenia, że skaner po prostu odrzucił krawędzie obrazu.

Każda rekonstrukcja uzupełnia obraz przy określonych założeniach o sygnale i akwizycji. Założenia te mogą zmieniać pozorny szum, wygładzenie, ostrość krawędzi, teksturę i widoczność małych struktur. Metoda wyuczona może utworzyć wiarygodnie wyglądający szczegół, którego spróbkowane dane niezależnie nie potwierdzają, albo usunąć szczegół uznawany przez jej założenia aprioryczne za szum. Jest to powód, aby oceniać metodę w zamierzonym zadaniu, a nie twierdzenie, że każda rekonstrukcja przyspieszona lub wyuczona zachowuje się w ten sam sposób.

Konsekwencja ilościowa jest bezpośrednia. Pomiar wykonany na obrazie zależy od sekwencji, wzorca próbkowania, algorytmu rekonstrukcji i ustawień, które utworzyły ten obraz. Obraz zrekonstruowany metodą wyuczoną nie staje się automatycznie wymienny z obrazem zrekonstruowanym konwencjonalnie, nawet jeśli oba mają te same wymiary macierzy i przedstawiają tę samą anatomię. Warunki rekonstrukcji i próbkowania należy zapisywać wraz z obrazem, tak jak rozdział o CT zapisuje kernel rekonstrukcyjny.

Gdy tworzenie obrazu rozumie się jako łańcuch obejmujący akwizycję i rekonstrukcję, następnym pytaniem jest zachowanie modelu po zmianie jednego z ogniw. Prowadzi to do zmiany domeny wywołanej akwizycją.

Źródła

• Stanford Medicine Body MRI Research Group, MRI Physics Education Materials — publiczne wizualne wprowadzenie do tworzenia obrazu MRI, częstotliwości przestrzennej, przestrzeni k, rozdzielczości i pola widzenia.
• Ridgway JP. Cardiovascular magnetic resonance physics for clinicians: part I. J Cardiovasc Magn Reson. 2010;12:71. doi:10.1186/1532-429X-12-71. (Pełny tekst)
• Larson P. Principles of MRI: The MRI Signal Equation. Educational resource, 2025. doi:10.5281/zenodo.5547020. (Dokumentacja techniczna)
• Moratal D, Valles-Luch A, Marti-Bonmati L, Brummer ME. k-Space tutorial: an MRI educational tool for a better understanding of k-space. Biomed Imaging Interv J. 2008;4:e15. doi:10.2349/biij.4.1.e15. (Pełny tekst)
• Heckel R, Jacob M, Chaudhari A, Perlman O, Shimron E. Deep learning for accelerated and robust MRI reconstruction. Magn Reson Mater Phy. 2024;37:335-368. doi:10.1007/s10334-024-01173-8. (Pełny tekst)