Artykuł metodologiczny

Ocena odkształcenia mięśnia sercowego w rezonansie magnetycznym z wykorzystaniem śledzenia cech (feature tracking)

14.3K wyświetleń

DOI:

10.3791/2356

12 lutego 2011

W tym artykule

Informacja o sprostowaniu

Important: There has been an erratum issued for this article. View Erratum Notice

Podsumowanie

Dokładna i praktyczna metoda pomiaru parametrów takich jak odkształcenie (strain) w tkance mięśnia sercowego ma dużą wartość kliniczną, ponieważ wykazano, że odkształcenie jest czulszym i wcześniejszym markerem dysfunkcji kurczliwości niż często stosowany parametr EF.

Streszczenie

Cel: Dokładna i praktyczna metoda pomiaru parametrów, takich jak odkształcenie (strain) w tkance mięśnia sercowego, ma ogromną wartość kliniczną, ponieważ wykazano, że odkształcenie jest czulszym i wcześniejszym markerem dysfunkcji kurczliwej niż często stosowany parametr EF. Obecne technologie CMR są czasochłonne i trudne do wdrożenia w praktyce klinicznej. Feature tracking to technologia, która może prowadzić do większej automatyzacji i powtarzalności ilościowej analizy obrazów medycznych przy mniejszym nakładzie czasu niż w przypadku porównywalnych metod.

Metody: Automatyczne lub manualne wprowadzenie danych w pojedynczej fazie służy jako inicjalizacja, od której system rozpoczyna śledzenie przemieszczenia poszczególnych wzorców reprezentujących struktury anatomiczne w czasie. Specyfiką tej metody jest to, że obrazy nie wymagają żadnej wcześniejszej manipulacji, takiej jak np. tagowanie obrazów CMR.

Wyniki: Metoda ta bardzo dobrze nadaje się do śledzenia tkanki mięśniowej, co pozwala na ilościową analizę mięśnia sercowego, a także przepływu krwi.

Wnioski: Ta nowa metoda oferuje solidną i oszczędzającą czas procedurę ilościowego określania parametrów przemieszczenia, prędkości i odkształcenia tkanki mięśnia sercowego oraz krwi w standardowych sekwencjach obrazowania CMR. Dzięki temu może ona zostać wdrożona w praktyce klinicznej.

Protokół

1. Wstęp

Automatyczne wykrywanie krawędzi jest fundamentalną kwestią w analizie obrazów. W obrazowaniu serca możliwość automatycznego wykrycia granicy wsierdzia w obrazowaniu lewej komory pozwoliłaby na obiektywny pomiar objętości komór oraz odkształcenia mięśnia sercowego (strain). Zostało to osiągnięte w echokardiografii dzięki technice śledzenia plam (speckle tracking). Opracowanie niezawodnych metod automatycznego wykrywania granic jest trudnym zadaniem, dla którego nie znaleziono jeszcze powszechnie uznanego i niezawodnego rozwiązania w rezonansie magnetycznym serca (CMR). W rzeczywistości w praktyce klinicznej granice są albo rysowane ręcznie przez operatora, albo wykrywane przez oprogramowanie jako interfejs między mięśniem sercowym a jamą 1,2. W niniejszym artykule przedstawiamy inne podejście, w którym granice nie są „wykrywane”, lecz „śledzone”, tzn. obserwowane w czasie, zaczynając od jednego wiarygodnego, istniejącego obrysu chwilowego, który zazwyczaj – choć niekoniecznie – jest rysowany ręcznie przez doświadczonego operatora na pojedynczej klatce obrazu. Poszczególne punkty tworzące taki pierwszy wiarygodny obrys są śledzone w czasie poprzez wyszukiwanie tych samych cech, które znajdują się w sąsiedztwie danego punktu w kolejnych klatkach. Śledzone cechy mogą być granicą jama-tkanka lub elementami anatomicznymi, które różnią się w obrębie tkanki. Są one wyszukiwane za pomocą metod największej wiarygodności w dwóch obszarach zainteresowania pomiędzy dwiema klatkami.

Lokalne przemieszczenie z klatki na klatkę jest równoważne ocenie lokalnej prędkości (stosunku przemieszczenia do odstępu czasu). Automatyczna ocena prędkości w danym punkcie jest określana na podstawie porównania przemieszczenia danych obrazowych wokół tego punktu w dwóch kolejnych klatkach. Metody te, w kilku różnych sformułowaniach, były stosowane w wielu dziedzinach badawczych. W zaawansowanej analizie obrazu zalicza się je do ogólnej kategorii znanej jako przepływ optyczny (Optical Flow) 3,4. W obrazowaniu echograficznym, gdy takie prędkości są wykorzystywane do śledzenia ruchu fizjologicznego, określa się je zazwyczaj jako śledzenie plamek (Speckle Tracking) 5, 6, ale znajdują one zastosowanie również w każdej innej metodzie obrazowania, np. w CMR, gdzie metody te nazywa się śledzeniem cech (feature tracking) lub śledzeniem granic (border tracking).

2. Materiały i metody

Metoda śledzenia cech

Brzeg wsierdziowy lub epikardialny obrazu cine 2D CMR jest ręcznie obrysowywany na jednej dowolnej klatce (patrz ryc. 1). Można również obrysować cechy śródmięśniowe. Taki brzeg definiuje się następnie jako sekwencję N punktów, określonych przez ich pary współrzędnych (Proszę podać tekst źródłowy do tłumaczenia.i,yiz i=1...NŚledzenie krawędzi odbywa się poprzez śledzenie każdego pojedynczego punktu; proces ten opiera się na algorytmie hierarchicznym w wielu skalach oraz na połączeniu technik śledzenia 1D, które gwarantują wyższą dokładność, i śledzenia 2D, które jest niezbędne do prawidłowego wykrywania cech rozciągniętych przestrzennie w 2D.

Aby w pierwszej kolejności uchwycić duże przemieszczenie geometryczne brzegu, śledzenie przeprowadza się w kierunku prostopadłym do samego brzegu, gdzie granica między jamą a tkanką jest najlepiej rozpoznawalna. Śledzenie w tym kierunku wykonuje się przy użyciu metody cięć transmuralnych w następujący sposób (patrz ryc. 3). Rysuje się linię przecinającą ścianę, przechodzącą przez dany punkt i prostopadłą do niej. Piksele pobrane wzdłuż linii transmuralnej są rozmieszczone w kolumnach, gdzie każda kolumna odpowiada jednej klatce sekwencji obrazów. W ten sposób ewolucja wzdłuż cięcia transmuralnego może zostać przedstawiona dla wszystkich chwil jednocześnie w reprezentacji dwuwymiarowej, w której jedną osią jest odległość wzdłuż linii, a drugą oś czas (patrz ryc. 2). Reprezentacja ta jest podobna do tzw. trybu M w echokardiografii, natomiast w CMR odpowiada ona bardziej funkcji „scout”. Aby poprawić jakość analizy, w przypadku obrazów słabej jakości z niskim stosunkiem sygnału do szumu, reprezentacja czasoprzestrzenna jest budowana z wykorzystaniem linii cięcia transmuralnego o grubości 5 pikseli. Następnie śledzenie brzegu jest wykonywane na tak uzyskanej mapie czasoprzestrzennej.

W drugim kroku, aby uwzględnić dwuwymiarowe przemieszczenie krawędzi, dla każdego punktu niezależnie przeprowadza się standardowe śledzenie 2D (w oparciu o przepływ optyczny) w ruchomym oknie o rozmiarze MxM, które jest zawsze wyśrodkowane na wcześniej oszacowanym punkcie krawędzi. Śledzenie 2D wykonuje się w dwóch etapach, w których połowa pierwszej estymacji jest wykorzystywana do wyśrodkowania ruchomych okien w drugim przejściu śledzenia. Następnie rozmiar okna jest redukowany z 32 do 16 w dwóch dodatkowych przejściach.

Aby zwiększyć dokładność śledzenia ruchu wzdłuż krawędzi, który służy do szacowania rotacji i torsji, śledzenie 1D jest wykonywane na obrazach czasoprzestrzennych zbudowanych z grubych przekrojów „równoległych” do zakrzywionej krawędzi (patrz rycina 3). W każdym punkcie, niezależnie, piksele pobrane wzdłuż ruchomej krawędzi, wycentrowane w punktach tej krawędzi, są układane w kolumnach, gdzie każda kolumna odpowiada jednej klatce z sekwencji obrazów. Aby poprawić jakość analizy i jak najlepiej uchwycić cechy krawędzi, linię przedłużono o 5 pikseli w głąb tkanki (podwsierdzie). Śledzenie krawędzi jest następnie wykonywane na obrazie czasoprzestrzennym zgodnie z opisaną powyżej procedurą. Aby zapewnić spójność przestrzenną śledzonej krawędzi, do przesunięć obliczonych w sąsiednich punktach na każdym etapie zastosowano medianowy filtr 3-punktowy oraz gaussowski filtr 3-punktowy (o wagach 0,25, 0,5, 0,25).

Śledzenie wzdłuż dwuwymiarowego obrazu czasoprzestrzennego

W tej sekcji opisano procedurę śledzenia krawędzi wzdłuż jednego kierunku na obrazie dwuwymiarowym (zbliżonym do trybu M), rozpoczynając od znanej pozycji w danym momencie.

Oś X definiuje kierunek poziomy, a oś y kierunek pionowy. Kolumny są oznaczone jako xi, i=1...M, gdzie M to liczba kolumn na obrazie. Śledzenie polega na wyznaczeniu dyskretnej sekwencji liczb rzeczywistych yi=y(xi), zaczynając od znanego punktu yk odpowiadającego kolumnom xk.

Przesunięcie z punktu znanego yk do punktu yk+1 jest szacowane poprzez ocenę korelacji wzajemnej między całą kolumną w xk a całą kolumną w xk+1. Funkcja korelacji wzajemnej wykaże maksimum, którego pozycja wskazuje wartość przesunięcia pionowego niezbędnego do zmaksymalizowania podobieństwa między dwiema kolumnami; zatem yk+1 jest szacowane poprzez dodanie takiego przesunięcia do yk. Procedura ta jest powtarzana dla wszystkich par sąsiednich kolumn, a wynikiem jest oszacowanie całej krawędzi yi, i=1...M. Korelacja wzajemna jest tutaj obliczana przy użyciu algorytmu szybkiej transformaty Fouriera (Fast Fourier Transform) w celu skrócenia czasu obliczeń.

Pierwsza estymata yi jest następnie doprecyzowana iteracyjnie. Aby osiągnąć ten cel, wyodrębnia się podzbiór obrazu, pobierając kilka punktów powyżej i poniżej poprzedniej estymaty yi, a następnie generuje się nowy obraz, którego centrum odpowiada sekwencji yi, i wykorzystuje go do śledzenia korekcji. Proces doprecyzowania jest powtarzany do momentu, aż nie zostanie stwierdzona żadna korekta.

Ulepszony i bardziej naturalny wynik uzyskuje się następnie za pomocą końcowej procedury „snake” [5], aby w obrazie czasoprzestrzennym śledzić poziom jasności obrazu przechodzący przez punkt stały yk. Cały proces wykorzystuje okresowość czasową, aby zapewnić okresowy wynik i uniknąć efektu dryftu.

Ograniczenia techniczne śledzenia cech

Technika śledzenia krawędzi, podobnie jak każda metoda śledzenia plamek (speckle tracking), opiera się na kwantyfikacji zmian jasności pikseli między poszczególnymi klatkami. Narzuca to dolną granicę prędkości, wynikającą z konieczności zaobserwowania plamki, która w jednej klatce znajduje się w jednym pikselu, a w następnej przesuwa się do sąsiedniego piksela. Granica ta jest zatem określona wzorem:
Równanie prędkości \( \varepsilon_v = k \frac{\Delta x}{\Delta t} \), wzór dla analizy dynamicznej.
Równanie 1
gdzie Δx to rozmiar piksela, a Δt to odstęp czasu między dwiema klatkami. Współczynnik k zależy od jakości algorytmu śledzenia oraz jego zdolności do oceny dynamicznych zmian w skali subpikselowej. Oznacza to, że prędkości znacznie przekraczające tę granicę są szacowane z dużą dokładnością, natomiast dokładność ta maleje, gdy wartości prędkości zbliżają się do tej granicy lub spadają poniżej niej.

Ograniczenie to oznacza również, że zwiększenie częstotliwości próbkowania klatek (zmniejszenie Δt) z jednej strony pozwala na łatwiejszą ocenę dużych prędkości i ich gwałtownych zmian (np. podczas faz izowolumetrycznych). Z drugiej strony, zwiększenie częstotliwości klatek (zmniejszenie Δt) podnosi tę granicę i wiąże się z obniżeniem dokładności oceny niższych prędkości, o ile nie towarzyszy temu podobny wzrost rozdzielczości przestrzennej (zmniejszenie Δx).

Przygotowanie obrazu fantoma

Przygotowano serię sztucznych, wygenerowanych komputerowo pętli, aby umożliwić testowanie procedury analizy obrazu w prostych i w pełni kontrolowanych warunkach. W tym celu przygotowano fantom w projekcji krótkiej osi idealnej lewej komory zgodnie z poniższym opisem.

Granice wsierdzia i osierdzia są reprezentowane przez dwa koncentryczne okręgi o promieniach odpowiednio R0(t) i R1(t). Obraz przygotowano poprzez nadanie pierścieniowi, który reprezentuje tkankę pomiędzy obiema granicami, jednolitego szarego koloru na czarnym tle. Następnie zastosowano liniowy filtr top-hat o rozmiarze 8x8, aby uniknąć niefizycznych nieciągłości.

Ruch epikardium wyrażono w [mm] jako R0(t)=10+5cos(2πt/T), gdzie T to okres bicia serca przyjęty jako T=1s. Teoretyczna kinematyka endokardium jest stała wzdłuż krawędzi i zależy wyłącznie od czasu; prędkość jest jedynie radialna i dana wzorem V0(t)= dR0/dt=-π sin(2πt/T), w [cm/s]. Odkształcenie procentowe, obliczone w stosunku do długości krawędzi w czasie zero, wynosi St0(t)=100x(R0(t)-R0(0))/R0(0)=100(cos(2πt/T) -1)/3, a szybkość odkształcenia wynika z równania (1) jako SR0(t)=10 V0/R0, w [s-1]. Przyjęto, że epikardium porusza się zgodnie ze stałą grubością, R1(t)= R0(t)+5mm, lub pozostaje w spoczynku R1(t)= R0(0)+5mm.

Każdy obraz jest kwadratem o rozmiarze 48mm, wycentrowanym na pierścieniu tkanki, i posiada rozdzielczość NxN. Przykładowe obrazy przedstawiono na rysunku 4, płyty a i b; profile czasowe odkształcenia i szybkości odkształcenia pokazano na rysunku 4, płyty c i d. Pętle są przygotowywane poprzez zmianę rozdzielczości N, częstotliwości klatek FR oraz rodzaju ruchu wsierdziowego.
Metoda śledzenia wsiercdzia jest stosowana do takich obrazów poprzez wyznaczenie na pierwszej klatce liczby Np punktów równomiernie rozmieszczonych wzdłuż okrągłego wsierdzia.

3. Reprezentatywne wyniki

Badanie na fantomie

W niniejszej pracy przeanalizowano zastosowanie metody analizy obrazu do obrazów fantomów generowanych komputerowo. Globalną miarę końcowego błędu obliczono za pomocą pierwiastka średniokwadratowego procentowej różnicy. Pierwiastek średniokwadratowy, błąd średni oraz maksymalny błąd odkształcenia wsierdzia zdefiniowano jako
Wzory matematyczne do szacowania błędów, przedstawiające równania dla analizy rezydualnej.
Równanie 2
gdzie St0(t) jest wartością dokładną, St(t) jest wartością obliczoną poprzez analizę obrazu, a sumowanie obejmuje wszystkie klatki NF=FRxT. Tą samą definicję zastosowano dla promienia, prędkości i szybkości odkształcenia. Śledzenie jest w zasadzie niezależne od pozycji wzdłuż wsierdzia, a różnice między poszczególnymi punktami są znacznie mniejsze niż 1%.

Wyniki podsumowano w Tabeli I dla 15 fantomów o różnej rozdzielczości przestrzennej, częstotliwości klatek oraz typie ruchu granicy epikardium; przedstawiono również efekt zmiany liczby punktów wykorzystanych do śledzenia granicy endokardium.

Błędy są we wszystkich przypadkach bardzo małe dla wartości całkowych (promień i odkształcenie) i nieco większe dla wartości różniczkowych (prędkość i szybkość odkształcenia), które są powiązane z pochodną tych pierwszych. Było to oczekiwane, ponieważ operator pochodnej potęguje błędy. Jakość wyników pogarsza się wraz ze zmniejszeniem rozdzielczości; w rzeczywistości dokładność jest powiązana z rozmiarem piksela, który reprezentuje (w szerokim sensie) minimalne przemieszczenie możliwe do odczytania z jednej klatki do drugiej. Rozdzielczość czasowa nie wpływa znacząco na wyniki, dopóki liczba klatek na sekundę jest wystarczająca; przy bardzo wysokiej częstotliwości klatek wyniki nie ulegają poprawie, ponieważ przemieszczenia między poszczególnymi klatkami stają się mniejsze niż rozmiar piksela. Pokazuje to, że zwiększenie liczby klatek na sekundę ma niewielką lub żadną użyteczność, gdy nie towarzyszy mu zwiększenie rozdzielczości przestrzennej.

Jednak rozważony tutaj prosty ruch sinusoidalny nie wymaga ekstremalnej rozdzielczości czasowej. Podobnie, zastosowanie zaledwie 8 punktów jest wystarczające, aby odwzorować prosty, kolisty kształt wsierdzia. Wyniki dotyczące wsierdzia nie są w znacznym stopniu zależne od rodzaju ruchu, jakiego poddawany jest epikardium. Potwierdziliśmy również, że wyniki nie są istotnie wpływowane przez przyjęte filtrowanie obrazu.

Wizualna prezentacja wyników przedstawiona jest na Rysunku 4, gdzie obliczona granica endokardium w dwóch chwilach została nałożona na obrazy fantoma (tablice a i b). Odkształcenie oraz szybkość odkształcania są przedstawione na (tablicach c i d) dla przypadku nr 1 oraz przypadku nr 8 o niższej rozdzielczości. Odkształcenie i szybkość odkształcania w przypadku nr 1 (kwadraty) wykazują doskonałą zgodność z wartością teoretyczną, przy średnim błędzie wynoszącym odpowiednio 0,6% i 3%. Zgodność jest nieco mniejsza w przypadku nr 8, w którym rozdzielczość obrazu jest dwukrotnie mniejsza, a błędy dla odkształcenia i szybkości odkształcania wynoszą odpowiednio 0,9% i 4,5%.

Walidacja kliniczna 1.

Porównano obwodowy odkształcenie mięśnia sercowego (εcc) w całym przekroju środkowej części lewej komory (LV) przy użyciu technik Harmonic Phase Imaging (HARP) oraz FT u 191 pacjentów z dystrofią mięśniową Duchenne'a, podzielonych według wieku i stopnia niewydolności serca, oraz u 42 osób kontrolnych dopasowanych pod względem wieku. Przeprowadzono retrospektywną analizę off-line dopasowanych przekrojów z obrazowania tagowanego oraz SSFP. Dla całej populacji badanej (n=233) średnie wartości FT εcc (-13,3 ± 3,8 %) były silnie skorelowane z HARP εcc (-13,6±3,4 %), przy współczynniku korelacji Pearsona wynoszącym 0,899. Średnie wartości εcc u pacjentów z DMD wyznaczone za pomocą HARP (-12,52 ± 2,69 %) i FT (-12,16 ± 3,12 %) nie różniły się istotnie (p=NS). Podobnie, średnie wartości εcc osób kontrolnych wyznaczone za pomocą HARP (-18,85 ± 1,86) i FT (-18,81 ± 1,83) nie różniły się istotnie (p=NS). Wyciągnięto wniosek, że ocena εcc oparta na metodzie FT wykazuje wysoką korelację z εcc pochodzącym z obrazów tagowanych w dużej populacji pacjentów z DMD z szerokim zakresem niewydolności serca.

Fantom nrCzęstotliwość klatekRozdzielczość. NRuchy osierdziaNpεRεRmaksymalnyεVεVmaksymalnyεSRεSRmaksymalneεStεStmaksymalne
132401z endo160.571.103.988.262.956.550.591.34
232401brak ruchu160.270.603.628.532.905.990.751.29
364401z endo160.350.621.754.923.0810.470.340.86
464401brak ruchu160.390.731.785.482.386.460.320.63
5128401z endo160.220.401.573.662.738.390.350.89
616401z endo160.410.9212.1018.277.1215.170.881.93
732201z endo160.430.823.847.034.217.790.831.52
832201brak ruchu160.490.904.008.544.5410.840.931.64
932101z endo162.583.875.7017.449.2218.274.446.77
1032401z endo640.320.733.988.093.118.860.481.10
1132401brak ruchu640.330.533.567.832.787.440.700.98
1232401z endo320.441.003.998.162.996.820.631.49
1332401z endo80.200.413.646.843.099.100.430.74
1432401brak ruchu80.220.413.367.422.765.930.581.01
15* 32401z endo161.352.425.549.687.0216.582.624.46

Tabela 1. Analiza fantomowa śledzenia brzegów wsierdziowych: pierwiastka średniokwadratowego oraz maksymalne błędy procentowe [%] obliczono dla głównych wielkości w odniesieniu do różnych parametrów fantomów. Parametry zaznaczone pogrubieniem wskazują zmiany względem Fantomu nr 1. Rozważono zależność od częstotliwości klatek, rozdzielczości oraz liczby śledzonych punktów. Wpływ rodzaju ruchu osierdzia rozważono dla dwóch przypadków granicznych: gdy brzeg osierdzia się nie porusza (brak ruchu) lub porusza się wraz z wsierdziem (brak pogrubienia). Ostatni fantom (*) został zbudowany bez filtrowania podstawowych obrazów skokowych z gwałtowną zmianą jasności w jednym pikselu. Błędy przekraczające 10% zaznaczono pogrubieniem.

Obrazowanie rezonansem magnetycznym (MRI) z zaznaczonymi obszarami, analiza obrazowania medycznego, narzędzie diagnostyczne.
Rysunek 1. Obraz cMR lewej komory w przekroju podłużnym (lewy obraz) i w przekroju poprzecznym (prawy obraz), z naniesionym obrysem granicy wsierdzia.

Obrazowanie ultradźwiękowe; tryb jasności; progresja klatek; wykres; analiza ilościowa; ruch mięśnia sercowego.
Rycina 2. Reprezentacja czasoprzestrzenna sekwencji obrazów, w której przestrzeń przebiega wzdłuż przekroju transmuralnego. Przekrój transmuralny wykonano zgodnie z punktem początkowym na rycinie 3. Przedstawiono ewolucję czasową automatycznie śledzonego punktu początkowego.

Analiza MRI serca, schemat: cięcia transmuralne a cięcia równoległe do zakrzywionej krawędzi, ocena tkanki serca.
Rysunek 3. Obraz lewej komory w widoku osi podłużnej, z cięciami transmuralnymi oraz cięciami równoległymi do zakrzywionej krawędzi.

Schematy modeli równowagi statycznej z wykresami fal sinusowych dla wyników analizy przestrzennej.
Rysunek 4. Badanie z użyciem fantomów. Dwa obrazy (przypadek nr 2) przy maksymalnym rozszerzeniu (płyta a) i skurczu (płyta b), obliczone punkty granicy wsiekoliczne są nałożone na siebie. Odkształcenie (płyta c) oraz szybkość odkształcenia (płyta d) obliczone dla dwóch różnych fantomów (przypadki nr 1 i nr 8) przedstawiono w porównaniu z wartościami efektywnymi.

Ocena funkcji serca; porównawcze przebiegi fal i parametry; wykres i analiza danych.
Rycina 5. Przykłady globalnego odkształcenia obwodowego (krzywa czarna) i odkształcenia segmentowego (krzywe kolorowe) u zdrowych pacjentów (a). Przykład globalnego odkształcenia obwodowego (krzywa czarna) i segmentowego odkształcenia obwodowego (krzywe kolorowe) u pacjentów z obniżoną funkcją lewej komory i blokiem lewej odnogi pęczka Hisa (b). Zwrócić uwagę na różny czas wystąpienia szczytowego odkształcenia obwodowego, co wskazuje na obecność dyssynchronii lewej komory

Dyskusja

Opisujemy nową metodę „feature tracking” (śledzenie cech) do pomiaru parametrów deformacji i przemieszczenia w tkance mięśnia sercowego na podstawie obrazów SSFP. Metoda ta została niedawno zwalidowana w odniesieniu do HARP w ocenie szczytowego globalnego odkształcenia obwodowego w dużej populacji pacjentów z dystrofią mięśniową Duchenne'a (1). Stwierdzono doskonałą korelację między tymi dwiema technikami. Technologia feature tracking jest stosunkowo prosta w użyciu i nie wymaga żadnych dodatkowych sekwencji obrazowania, takich jak tagging, SENSE czy DENSE. Oprócz klinicznie zwalidowanej oceny szczytowego globalnego Ecc, można uzyskać również regionalne dane o odkształceniu; FT wydaje się być solidną techniką zdolną do szacowania odkształcenia podłużnego i promieniowego LV, a także odkształcenia podłużnego prawej komory.

Konieczne jest przeprowadzenie dalszych badań, aby dokładnie określić ograniczenia tej techniki w zakresie rozdzielczości czasowej i przestrzennej obrazu. FT umożliwia porównania z innymi modalnościami, takimi jak echokardiografia, ponieważ podstawowa zasada działania jest taka sama.

Podsumowanie

Śledzenie cech (feature tracking) ma potencjał do wdrożenia w praktyce klinicznej, ponieważ eliminuje potrzebę czasochłonnej analizy oraz dodatkowych procedur akwizycji danych. FT może być wykorzystywane do oceny globalnego i segmentowego odkształcenia (strain) lewej i prawej komory oraz dyssynchronii lewej komory. Trwają dalsze badania kliniczne mające na celu ocenę użyteczności FT w kardiologii klinicznej.

Oświadczenia

Rolf Baumann jest pracownikiem firmy TomTec Imaging Systems GmbH, Niemcy

Materiały

Lista materiałów użytych w tym artykule
NazwaFirmaNumer katalogowyKomentarze
Oprogramowanie do śledzenia cech Diogenes®TomTec Imaging Systems, Fulda, Niemcy

Bibliografia

  1. Hor, K. N., Gottliebson, W. M., Carson, C., Wash, E., Cnota, J., Fleck, R., Wansapura, J., Klimeczek, P., Al-Khalidi, H. R., Chung, E. S., Benson, D. W., Mazur, W. Magnetic Resonance Feature Tracking: Comparison of Feature Tracking Method for Strain Calculation with Harmonic Phase Imaging Analysis . JACC Cardiovasc Imaging. 3 (2), 152-154 (2010).
  2. Maret, E., Todt, T., Brudin, L., Nylander, E., Swahn, E., Ohlsson, J. L., Engvall, J. E. Functional measurement based on feature tracking of cine magnetic resonance images identify left ventricular segments with myocardial scar. Cardiovascular Ultrasound. 7, 53-53 (2009).
  3. Singh, A. Optic Flow Computation: A Unified Perspective. , IEEE Comput. Soc. Press. Piscataway, NJ. (1992).
  4. Barron, J. L., Fleet, D. J., Beauchemin, S. Performance of optical flow techniques. International Journal of Computer Vision. 12, 43-77 (1994).
  5. Bohs, L. N., Geiman, B. J., Anderson, M. E., Gebhart, S. C., Trahey, G. E. Speckle tracking for multi-dimensional flow estimation. Ultrasonics. 38, 369-375 (2000).
  6. Malpica, N., Santos, A., Zuluaga, M. A., Ledesma, M. J., Pérez, E., Garcia-Fernandez, M. A., Desco, M., M, Tracking of Regions-of-Interest in myocardial contrast echocardiography. Ultrasound in Med. & Biol. 30, 303-309 (2004).

Przedruki i uprawnienia

Sprostowanie


Formal Correction: Erratum: Magnetic Resonance Derived Myocardial Strain Assessment Using Feature Tracking
Posted by JoVE Editors on 11/10/2014. Citeable Link.

A correction was made to Magnetic Resonance Derived Myocardial Strain Assessment Using Feature Tracking. An author's given name was updated.

The author's name was updated from:

Woodrow Benson

to:

D. Woodrow Benson

Tagi

Rezonans magnetyczny sercaodkszta cenie obwodoweanaliza odkszta ceniasekwencja w osi kr tkiejgranica endokardiumledzenie przemieszczepomiar pr dko ciwalidacja na fantomie