基于彌散張量成像技術(shù)的腦白質(zhì)纖維跟蹤技術(shù)研究進(jìn)展

基于彌散張量成像技術(shù)的腦白質(zhì)纖維跟蹤技術(shù)研究進(jìn)展

人類大腦是自然界最復(fù)雜的系統(tǒng)之一，白色腦纖維是構(gòu)成不同腦區(qū)域信息交換的物質(zhì)基礎(chǔ)。傳統(tǒng)的解剖染色方法在研究腦白質(zhì)纖維時(shí),因其有創(chuàng)性而只能局限于動(dòng)物和人類尸體,而不適合活體人腦的研究。隨著核磁共振技術(shù)的興起,尤其是彌散張量成像(diffusiontensorimaging,DTI)技術(shù)和腦白質(zhì)纖維束跟蹤(whitemattertractography,WMT)技術(shù)的出現(xiàn),使得活體人腦的白質(zhì)纖維解剖結(jié)構(gòu)能夠無創(chuàng)地三維重建出來DTI成像的基本原理彌散(diffusion)是指分子的隨機(jī)移動(dòng),即布朗運(yùn)動(dòng)。DTI是一種利用組織中水分子彌散的各向異性(anisotropy)引起的磁共振信號(hào)變化,來探測(cè)組織微觀結(jié)構(gòu)的成像方法,主要用于研究中樞神經(jīng)系統(tǒng)的解剖神經(jīng)束和白質(zhì)纖維解剖結(jié)構(gòu)的顯示。實(shí)驗(yàn)表明:水分子在有髓纖維中的擴(kuò)散表現(xiàn)出較高的各向異性,尤其在平行于腦白質(zhì)纖維的方向上彌散最快。因此,可以通過張量場(chǎng)中各向異性彌散的方向信息來跟蹤腦白質(zhì)纖維的走行方向,獲取白質(zhì)纖維束的微觀解剖結(jié)構(gòu)信息,從而提取和重建白質(zhì)纖維束描述DTI圖像特征的常用參數(shù)包括:反映體內(nèi)水分子向各個(gè)方向彌散平均值的表觀彌散系數(shù)(apparentdiffusioncoefficient,ADC);相對(duì)各向異性(relativeanisotropy,RA)和部分各向異性(fractionalanisotropy,FA)。其中,RA和FA是最常用的描述彌散各向異性的參數(shù)。實(shí)驗(yàn)表明:水分子彌散越明顯,ADC值越高;并且RA和FA的值具有一致性,代表著水分子在彌散主向量軸上的運(yùn)動(dòng)強(qiáng)度。1999年Mori等首次提出WMT算法,之后隨著計(jì)算機(jī)技術(shù)的發(fā)展和臨床應(yīng)用的需求,涌現(xiàn)出各種WMT算法。隨著研究的深入,纖維束跟蹤算法從最初的單張量模型(single-tensormodel)發(fā)展到多張量模型(multi-tensormodel),出現(xiàn)了確定性算法(deterministicfibertracking)和概率算法(probabilisticfibertracking)兩類比較典型的算法。近幾年,隨著數(shù)據(jù)采集新方法和圖像處理新理論的出現(xiàn),基于高角分辨率彌散成像(highangularresolutiondiffusionimaging,HARDI)數(shù)據(jù)的跟蹤算法和基于壓縮感知(compressedsensingorcompressedsampling,CS)理論的跟蹤算法成為該領(lǐng)域新的研究熱點(diǎn)。腦白質(zhì)纖維束跟蹤算法的基本原理腦白質(zhì)纖維跟蹤算法是基于擴(kuò)散線(streamlinepropagation)的方法,即線性方法,其基本原理是:利用局部張量信息進(jìn)行纖維跟蹤,即從一個(gè)被稱為“種子點(diǎn)”的體素為起始點(diǎn)開始跟蹤,計(jì)算該種子點(diǎn)前進(jìn)方向的向量,沿該向量方向跟蹤一段距離之后,再以軌跡上新的體素作為起點(diǎn)繼續(xù)跟蹤,反復(fù)迭代上述過程,最后將這些點(diǎn)連接起來就可以代表白質(zhì)纖維的走行比較有代表性的線性追蹤算法有:流線性跟蹤算法(streamlinetracking,STT)確定性纖維束追蹤算法以上述線性算法為基本原理的確定性算法是比較成熟的一類WMT算法。該類算法通常利用體素所含的張量信息,根據(jù)設(shè)定的追蹤策略來連接相鄰體素,算法的實(shí)現(xiàn)過程一般包括選取種子點(diǎn)、確定纖維追蹤終止條件和纖維選擇策略等3個(gè)步驟。根據(jù)采用的數(shù)學(xué)模型不同可分為單張量模型(single-tensormodel)的算法和多張量模型(multi-tensormodel)的算法。單張量模型的確定性算法單張量模型認(rèn)為每個(gè)體素僅穿過一條纖維,可利用求解張量場(chǎng)主特征向量的特征值來跟蹤白質(zhì)纖維束的走行STT算法:由Mori等TEND算法:該算法與線性算法略有不同,它不以當(dāng)前種子點(diǎn)的主特征向量方向?yàn)樽粉櫡较?而是以主特征向量發(fā)生偏轉(zhuǎn)的方向?yàn)闇?zhǔn)則進(jìn)行追蹤。這種算法重建出的纖維束長且平滑,無角度突變,在各向異性程度小的區(qū)域具有一定的優(yōu)勢(shì)。但是,該算法對(duì)所有的彌散張量都不加區(qū)別地采取偏轉(zhuǎn)處理,在各向異性程度大的區(qū)域跟蹤纖維束時(shí)容易出現(xiàn)偏差;而且該算法要用整個(gè)彌散張量進(jìn)行計(jì)算,數(shù)據(jù)量較大,算法的準(zhǔn)確性還有待進(jìn)一步的臨床驗(yàn)證FACT算法(fiberassignmentbycontinuoustracking):該算法認(rèn)為纖維束在神經(jīng)系統(tǒng)中的軌跡是空間三維曲線,曲線上每點(diǎn)的方向?yàn)樵擖c(diǎn)處正切向量,該向量對(duì)應(yīng)于該點(diǎn)的主特征向量。在跟蹤時(shí),FACT算法以人為設(shè)置的感興趣區(qū)ROI(regionofinterested)內(nèi)的體素為種子點(diǎn),以種子點(diǎn)的主特征向量方向?yàn)槔w維束走形方向,并向相反兩個(gè)方向延長擴(kuò)展纖維束,直至到達(dá)種子點(diǎn)的邊界為止,然后反復(fù)這樣迭代。其中,FA值、偏轉(zhuǎn)角度、延長長度是該算法設(shè)置追蹤終止條件的3個(gè)重要參數(shù)。計(jì)算簡單、跟蹤高效是該算法的特點(diǎn),同時(shí)還可以結(jié)合BruteForce方法和選取多個(gè)感興趣區(qū)域來減少追蹤的錯(cuò)誤。快速、簡單以及能在各向異性較強(qiáng)區(qū)域較好地跟蹤纖維束走行方向是FACT算法的優(yōu)點(diǎn)。但是,如果體素有多條纖維通過,由于無法辨別各條纖維的走形方向,造成該算法在追蹤此類神經(jīng)纖維束時(shí)會(huì)產(chǎn)生一定的誤差GTRACT算法(guidedtensorrestoreanatomicalconnectivitytractography):該算法的核心思想是在追蹤迭代過程中,用搜索的解剖信息來引導(dǎo)追蹤,纖維路徑上每一點(diǎn)方向僅由當(dāng)前體素或鄰近體素決定。由于GTRACT算法對(duì)圖像噪聲不敏感,所以可以得到較好效果,而且該算法還能夠部分解決復(fù)雜區(qū)域的纖維束交叉問題單張量模型算法的數(shù)據(jù)模型相對(duì)簡單,容易實(shí)現(xiàn),但是追蹤效果受到追蹤方向、ROI選擇、終止準(zhǔn)則的設(shè)定等諸多因素的影響。為提高算法的準(zhǔn)確性,很多研究從這些因素入手對(duì)算法進(jìn)行優(yōu)化,如:用貝葉斯張量歸一化(Bayesiantensorregularization)算法來優(yōu)化追蹤方向,用自動(dòng)投影的算法精確ROI的選擇,結(jié)合先驗(yàn)解剖知識(shí)設(shè)置終止準(zhǔn)則等多張量模型的確定性算法DTI成像的分辨率通常為毫米級(jí),而神經(jīng)纖維直徑的實(shí)際大小是微米級(jí),因此單個(gè)體素本身含有的很多纖維束的信息被忽略了。多張量模型以體素的單張量信息為基本點(diǎn),將其分解為多個(gè)特征向量來跟蹤單個(gè)體素內(nèi)的多條腦白質(zhì)纖維束走行方向雙張量模型(two-tensormodel)的確定性算法:該算法認(rèn)為單個(gè)體素(多為扁平體素)含有兩條交叉纖維束的信息,以預(yù)先設(shè)定的幾何限制為終止條件,為跟蹤兩條交叉纖維各自的走行,采用最小二乘法將扁平體素代表的單張量分解為雙張量,其中每個(gè)張量分別表示一條纖維束的走行方向,從而模擬出交叉纖維束的走行方向高階張量模型(highorder-tensormodel)的確定性算法:在雙張量模型的基礎(chǔ)上,有些研究者利用HARDI數(shù)據(jù),基于任意階次的正定高階張量理論概率纖維束追蹤算法與確定性算法不同,概率纖維束跟蹤算法將單個(gè)張量看作是各向異性與各向同性模型的合成形式,通過引入相關(guān)的概率(決策)算法完成張量求解。概率纖維束追蹤算法以兩個(gè)或更多特定區(qū)域之間相互連接的可能性概率分布為著眼點(diǎn),用概率分布函數(shù)(probabilitydistributionfunction,PDF)確定纖維束追蹤的主方向,并對(duì)體素之間的相關(guān)性進(jìn)行量化,使用張量模型來估計(jì)每個(gè)體素最大可能的纖維束走形方向概率纖維追蹤算法的優(yōu)點(diǎn)是:可以減少部分容積效應(yīng)與噪聲的影響,提高追蹤的精度;既能顯示大纖維束,又能顯示小纖維束,還可進(jìn)行纖維束分叉的計(jì)算。但是,概率纖維追蹤算法的計(jì)算量比確定性追蹤算法大、耗時(shí)長,所以達(dá)不到臨床應(yīng)用的要求,而且數(shù)據(jù)獲取的方式對(duì)交叉纖維束追蹤結(jié)果也有一定影響。下面介紹一種研究得比較多的概率跟蹤算法———全局優(yōu)化算法(globaloptimizationalgorithm)。由于噪音和方向的不確定性受部分容積效應(yīng)的影響,按彌散張量的主特征向量進(jìn)行白質(zhì)纖維束跟蹤時(shí)易出現(xiàn)錯(cuò)誤累積纖維束追蹤算法研究的新動(dòng)向基于HARDI數(shù)據(jù)的纖維束追蹤算法臨床MRI數(shù)據(jù)的分辨率通常為3~15mm與基于DTI數(shù)據(jù)的纖維束追蹤技術(shù)相比,HARDI技術(shù)可選多個(gè)方向,并采用現(xiàn)有的各種算法進(jìn)行追蹤,能夠較準(zhǔn)確地解決纖維束交叉的問題基于CS理論的纖維束追蹤算法目前比較成熟的WMT算法(如確定性追蹤算法、概率追蹤算法等)由于受到物理設(shè)備和病理特征等客觀因素和算法本身局限性的制約,在處理時(shí)間和精度上不能達(dá)到臨床應(yīng)用的要求CS理論是由美國科學(xué)家于2006年提出的一種新的信息獲取指導(dǎo)理論CS理論的特點(diǎn)非常符合醫(yī)學(xué)成像領(lǐng)域的應(yīng)用需求,在MRI成像方面有很好的應(yīng)用前景,這主要基于MRI成像本身的特點(diǎn):首先,MRI成像原理決定了其在頻率空間中對(duì)圖像進(jìn)行采樣,而不是直接在像素的級(jí)別上進(jìn)行采樣;其次,MRI圖像本身具有可壓縮性,能夠在變換域中稀疏編碼(如小波變換、傅立葉變換)。因此,根據(jù)MRI成像原理,利用探測(cè)到的自由感應(yīng)衰減信號(hào)與體素灰度圖像之間的二維傅立葉變換關(guān)系(即體素灰度圖像在二維傅立葉空間中所具有的稀疏性),在CS理論指導(dǎo)下,隨機(jī)稀疏采樣獲取信號(hào),然后根據(jù)圖像的先驗(yàn)知識(shí)重構(gòu)出完整圖像HARDI技術(shù)能對(duì)諸如交叉、分叉、扇型等復(fù)雜結(jié)構(gòu)的纖維束較準(zhǔn)確成像,但要獲取多個(gè)梯度的數(shù)據(jù),代價(jià)較高,而CS理論采樣的稀疏性恰好可以彌補(bǔ)HARDI方法的這一缺陷。國外有研究者將HARDI與CS結(jié)合起來進(jìn)行了初步嘗試:在采集少量3T-MRI數(shù)據(jù)的前提下,以HARDI+CS為技術(shù)路線,用多方向的確定性纖維束追蹤算法(multidirectionaldeterministicfibertractography)三維重建出白質(zhì)纖維束的復(fù)雜結(jié)構(gòu)WMT的主要應(yīng)用DTI是當(dāng)前神經(jīng)解剖和生理研究領(lǐng)域唯一的非侵入性醫(yī)學(xué)成像手段,而各類WMT算法在生理研究和疾病的診斷與治療中發(fā)揮了巨大作用,是傳統(tǒng)MRI方法的有力補(bǔ)充,其應(yīng)用主要包括生理研究和臨床應(yīng)用兩個(gè)方面。生理研究21世紀(jì)以來,人腦結(jié)構(gòu)連接網(wǎng)絡(luò)的研究一直是神經(jīng)科學(xué)中的重要課題臨床應(yīng)用DTI成像技術(shù)在評(píng)價(jià)神經(jīng)疾病的病理改變以及組織結(jié)構(gòu)與功能間的關(guān)系等方面具有重要的臨床價(jià)值,尤其是在腦白質(zhì)纖維束的追蹤觀察對(duì)腦部病變的診斷和治療方面。由于許多腦部病變都會(huì)引起ADC和FA值發(fā)生改變,因此利用腦白質(zhì)纖維的追蹤技術(shù)可以觀察纖維束與腦部腫瘤的位置關(guān)系,分析腫瘤浸潤腦內(nèi)纖維束的情況,從而確定白質(zhì)纖維束的受損程度此外,WMT技術(shù)在研究腦神經(jīng)的生長、發(fā)育、退化等方面也極具價(jià)值總結(jié)與展望WMT是目前研究腦白質(zhì)纖維束唯一的非侵入無創(chuàng)性方法,使許多領(lǐng)域的研究有了質(zhì)的飛躍,例如:一些主要的傳導(dǎo)束三維重建的結(jié)果與傳統(tǒng)神經(jīng)解剖有較高的一致性。臨床上導(dǎo)航信息與纖維束跟蹤技術(shù)的結(jié)合已成為數(shù)字醫(yī)學(xué)科研單位的研究重點(diǎn),其目的是在術(shù)中準(zhǔn)確定位白質(zhì)纖維束與腫瘤之間的位置關(guān)系。本文對(duì)幾類比較有代表性的WMT算法和最近幾年出現(xiàn)新的纖維束算法進(jìn)行了歸納總結(jié)。但是,至今還沒有一種通用于各類情況的普適跟蹤算法,每一種算法都有自己的特點(diǎn)、優(yōu)勢(shì)和局限性,在實(shí)踐中要結(jié)合實(shí)際情況選擇使用。纖維束追蹤技術(shù)從出現(xiàn)到應(yīng)用于特定的臨床領(lǐng)域,為滿足精度要求經(jīng)歷了各種嘗試,目前確定性追蹤算法被廣泛使用