楊楊，曾軍，范文亮，雷子喬，余建明

磁共振成像是一種功能強(qiáng)大的診斷、預(yù)后和治療評估工具，但通常只能定性而不是定量的顯示組織之間的對比以用于表征潛在的病理學(xué)信息。近年來，定量MRI技術(shù)的出現(xiàn)克服了這一局限性(如T1-和T2-mapping)，但此類方法一次僅能測量一個參數(shù)，時間效率低下，掃描時間過長，所獲得的圖像需擬合數(shù)學(xué)模型用以評估目標(biāo)參數(shù)，如弛豫時間(T1)或信號衰減時間(T2)[1-2]。后續(xù)的研究提出了幾種方法來縮短采集時間或在單次采集中提供T1和T2的組合測量[3-11]。然而，臨床應(yīng)用提出了更高的要求，希望能既快速又準(zhǔn)確地進(jìn)行定量參數(shù)的測量，磁共振指紋(MRF)技術(shù)的誕生使其成為可能[12-14]。這項技術(shù)旨在實現(xiàn)一次采集即可提供多個參數(shù)的測量，例如T1、T2、相對自旋密度、B0和磁場不均勻性(偏共振頻率)等。 MRF完全改變了定量MRI的執(zhí)行方式，定量化分析組織器官的特性，使醫(yī)師能夠更好地區(qū)分健康和病理組織，在隨訪研究中能更容易客觀地比較不同的檢查結(jié)果[15]，并且較常規(guī)的加權(quán)成像方法能更準(zhǔn)確的反映細(xì)胞水平的潛在變化[16-17]。定量成像對評估具有微觀特征變化的疾病至關(guān)重要，如心臟彌漫性纖維化、鐵或肝臟脂肪沉積[18-20]。

MRF成像技術(shù)原理

近兩年GE公司推出的MAGIC技術(shù)可以通過一次掃描序列獲取幾乎涵蓋臨床應(yīng)用中各種所需的對比度圖像，如T1WI、T2WI、FLAIR圖、STIR圖、相位敏感反轉(zhuǎn)恢復(fù)(phase sensitive inversion recovery,PSIR)圖和雙反轉(zhuǎn)恢復(fù)(Dual IR,DIR)圖等，同時得到T1、T2以及質(zhì)子密度(proton density,PD)的定量信息，是一種全新的MR影像采集方式。實現(xiàn)這一技術(shù)的關(guān)鍵是磁共振弛豫率與質(zhì)子密度定量圖譜的獲取。在磁共振成像中，了解了組織的弛豫信息，即可根據(jù)“Bloch方程”獲得組織在任何一種可能的參數(shù)設(shè)定下所呈現(xiàn)的圖像信息。MRF技術(shù)與Magic技術(shù)兩者雖然在技術(shù)關(guān)鍵點的實現(xiàn)方式上有所不同，但都是一次掃描即可進(jìn)行多參數(shù)測量，并且根據(jù)數(shù)學(xué)模型推導(dǎo)出組織信息。

1.MRF成像序列

最初的定量MR成像方法需要進(jìn)行多次掃描獲得所需數(shù)據(jù)，不斷重復(fù)相同的采集模式，如射頻激發(fā)角(翻轉(zhuǎn)角)、TR和梯度模式，直到獲得傅里葉域(也被稱為k空間)中的所有所需數(shù)據(jù)，然后使用傅立葉變換重建每個圖像，并且對每個體素進(jìn)行非線性擬合處理。對于MRF成像，翻轉(zhuǎn)角、TR和采集軌跡在整個采集過程中以偽隨機(jī)的方式變化，為每個組織生成不相干的信號，從而提供用于識別組織信息的獨特“指紋”。MRF序列是基于平衡穩(wěn)態(tài)自由進(jìn)動(balance steady state free precession,bSSFP)或真穩(wěn)態(tài)自由進(jìn)動快速成像(true fast imaging with steady state precession,TrueFISP)序列，它對T1、T2和偏共振頻率敏感，并且該序列產(chǎn)生的穩(wěn)態(tài)信號已經(jīng)被我們熟知[21]。其中翻轉(zhuǎn)角以正弦方式變化從而平滑地改變磁化的瞬態(tài)，變化范圍在0～60°和0～30°之間交替，產(chǎn)生具有250個時間點(圖像)的周期。在該信號的基礎(chǔ)上，添加隨機(jī)變化以引起具有相似參數(shù)的組織時間演變的差異。在每半個周期之后，將50個翻轉(zhuǎn)角設(shè)置為0以允許信號恢復(fù)。TR的變化基于Perlin噪聲[22]，范圍為9.34～12.00 ms。上述只是參數(shù)如何隨機(jī)變化的1個例子，其它的隨機(jī)模式已經(jīng)過測試，證實MRF技術(shù)并不局限于一組特定的參數(shù)設(shè)置[23]。在采集序列開始時發(fā)射反轉(zhuǎn)恢復(fù)脈沖，以增強(qiáng)組織之間T1的差異。對于每個TR，重建嚴(yán)重欠采樣的圖像。得到的原始圖像中的每個體素都包含一個簽名指紋，稍后將用于匹配識別。獲取的圖像總數(shù)(也稱為時間點)可能因采集而異，范圍為1000～2500，取決于圖像分辨率、欠采樣率和使用的匹配方法等。采集軌跡在大多數(shù)情況下使用可變密度螺旋軌跡[24](圖1b)，因其具有最小時間梯度和零力矩補(bǔ)償。例如，有一種填充軌跡的方法已成功應(yīng)用于128×128的矩陣，需要1個插入來完全采樣k空間的中心和48個交錯插入以完全采樣k空間的外部區(qū)域。在256×256矩陣的情況下，需要24個插入以完全采樣內(nèi)部區(qū)域和48個插入以完全采樣外部區(qū)域。在每個TR內(nèi)，獲取一個插入并用于重建圖像。隨后TR中的插入與前一個相比旋轉(zhuǎn)7.5°(≈2π/48)。隨著MRI技術(shù)的發(fā)展，其它類型的序列也可以實現(xiàn)MRF技術(shù)。例如，MRF技術(shù)應(yīng)用于穩(wěn)態(tài)進(jìn)動序列FISP以避免出現(xiàn)在寬視場掃描或高場強(qiáng)掃描儀中出現(xiàn)的條帶偽影。FISP序列仍然對T1和T2敏感，但對偏共振頻率不敏感。由于每個TR內(nèi)的不平衡梯度導(dǎo)致產(chǎn)生的信號是體素內(nèi)的自旋總和，使得序列免受條帶偽影的影響。不平衡梯度導(dǎo)致FISP序列與TrueFISP相比具有更短的瞬態(tài)，由于這個原因，翻轉(zhuǎn)角需要不同的偽隨機(jī)變化以保持信號和欠采樣圖像之間的不相干性并且識別產(chǎn)生的指紋：翻轉(zhuǎn)角基于正弦變化，在達(dá)到最大值的半周期內(nèi)進(jìn)行隨機(jī)變化，范圍為5°～90°。TR的變化基于Perlin噪聲模式，其范圍為11.5～14.0 ms。

2.字典生成

字典可以被視為MRF框架的核心。它是包含從采集數(shù)據(jù)中可以觀察到的所有可能的信號演變的數(shù)據(jù)庫，并且可以識別每個體素內(nèi)的組織。與法醫(yī)指紋識別過程一樣，MRF僅在數(shù)據(jù)庫足夠大時才有效。在MRF中，使用在采集期間模擬自旋行為的算法在計算機(jī)上生成字典，從而預(yù)測實際的信號演變。在基于TrueFISP的采集情況下，給定一組感興趣的組織參數(shù)，采用Bloch方程模擬采集序列對自旋的各種影響[26]。因此，MRF可以檢索的信息與模擬的物理效應(yīng)相關(guān)。在研究初期，MRF中僅包括對T1、T2和偏共振的模擬，但現(xiàn)在可以模擬和提取更多的組織特征，如部分體積、擴(kuò)散和灌注[13,27]。字典的一個關(guān)鍵方面是它的大小：為了確保識別采集中存在的任何可能的組織參數(shù)，需要模擬T1、T2和偏共振頻率的廣泛組合。例如一種參數(shù)范圍的標(biāo)準(zhǔn)TrueFISP字典總共有363,624種可能的組合，在使用基于C++腳本的標(biāo)準(zhǔn)臺式計算機(jī)計算1000個時間點的這種字典需要大約2.5min、2.5GB內(nèi)存。字典容量和分辨率的進(jìn)一步增加將提高所獲得圖譜的準(zhǔn)確性，代價是重建時間的增加和對存儲器有更高的要求。

與上述方法相比，F(xiàn)ISP采集序列的模擬則不同。由于FISP采集需要模擬不同頻率的多個等色線，然后將它們組合在一起，通過Bloch方程的模擬過程是耗時的。一種高效的模擬是擴(kuò)展相位圖(extended phase graph，EPG)[28]，其中受序列影響的自旋系統(tǒng)可以表示為離散的相位狀態(tài)集，非常適合模擬由不平衡梯度導(dǎo)致強(qiáng)烈失相位的自旋信號的演化過程。上文提到，與TrueFISP采集相比，F(xiàn)ISP序列對偏共振效應(yīng)不太敏感，因此相應(yīng)的字典僅包括T1和T2作為目標(biāo)參數(shù)。在標(biāo)準(zhǔn)臺式計算機(jī)上18,838個字典條目可以在大約8分鐘內(nèi)計算完成、1.2GB內(nèi)存。

3.匹配

在數(shù)據(jù)采集之后，將每個體素的指紋標(biāo)準(zhǔn)化為規(guī)范單位并與所有標(biāo)準(zhǔn)化后的字典條目進(jìn)行比較以識別給定體素中的組織，通過獲取體素信號和每個模擬指紋信號之間的內(nèi)積，返回最高值的條目被認(rèn)為是最能代表組織特性的條目，并且相應(yīng)的T1、T2和偏共振值被分配給該體素。相對自旋密度(M0)圖被計算為獲取的指紋和模擬的指紋之間的縮放因子。內(nèi)積法已被證明是一種合理且可靠的操作方法，即便在由于欠采樣或甚至存在有限量的運動偽影而導(dǎo)致低SNR的情況下也能正確地分辨組織。

由于不同信號演變之間的不連貫性，該方法還具有區(qū)分單個體素內(nèi)存在的不同組織成分(部分體積效應(yīng))的潛力。含有不同組織體素的指紋(S)可以看作不同組分(D)的加權(quán)(S=Dw)。如果不同的組分是先驗已知的，前一個等式的適當(dāng)逆解D-1S=w，將為每個體素提供每個不同組織的權(quán)重[29-30]。對于MR框架的臨床實用性而言，在合理的時間內(nèi)執(zhí)行此操作至關(guān)重要，而使用內(nèi)積的直接匹配是準(zhǔn)確的。

通過在時間維度或參數(shù)組合維度壓縮字典可以潛在地加速匹配，從而減少需要執(zhí)行的比較總數(shù)。對于TrueFISP字典，應(yīng)用奇異值分解(singular value decomposition,SVD)在時間維度上壓縮字典，可以將匹配時間減少3.4倍，F(xiàn)ISP字典最多為4.8倍?；赟VD的字典壓縮方法的估計參數(shù)精度減少率不到2%，字典被投影到由SVD獲得的前25～200個奇異向量所跨越的較低維度的子空間。將獲取的指紋投影到相同的子空間上，并且使用投影的信號和壓縮的字典來執(zhí)行匹配。該框架減少了計算的數(shù)量，盡管在子空間上添加了數(shù)據(jù)投影的操作，但是仍然減少了最終的計算時間。

減少匹配計算時間的另一種方法是減少參數(shù)組合維度。其中一種快速匹配算法中具有強(qiáng)相關(guān)性的字典條目被組合在一起，生成最能代表該組織的新信號[31]。首先，所獲取的指紋與每個組的代表信號匹配，僅考慮返回最高相關(guān)性的組然后匹配；然后，找到指紋與參數(shù)分配的剩余字典條目之間的最佳擬合即可。

4.欠采樣和運動

MRF獲得的參數(shù)圖是模式識別算法的結(jié)果，與傳統(tǒng)的重建技術(shù)相反。這使得MRF對各種圖像偽影更加穩(wěn)定。通過FA、TR和填充軌跡的隨機(jī)變化使得這項技術(shù)得到加強(qiáng)，不僅可以區(qū)分不同組織的指紋，而且還能增加指紋之間的不相干性。只要噪聲或欠采樣偽影與信號不相干，匹配就可以識別基礎(chǔ)信號演變，即使在低信噪比或加速條件下也是如此。另外，即使使用模糊或部分指紋的情況下，MRF也夠在有運動狀態(tài)影響的情況下提供參數(shù)圖而沒有任何殘余運動偽影。

MRF在臨床研究中的應(yīng)用

通過MRF研究獲得的T1、T2值與常規(guī)T1-和T2-mapping以及文獻(xiàn)中公布的方法得出的結(jié)果顯示出良好的相關(guān)性[13-14,32-33]。在臨床前和動物研究中也已證明了MRF的效用[34-35]。MRF的臨床應(yīng)用主要集中在大腦和前列腺上，其在腹部、肌肉骨骼和心臟應(yīng)用中也顯示出了潛力。

1.腦弛豫測量

關(guān)于MRF描述中顯示灰質(zhì)和白質(zhì)的T1和T2值具有良好的文獻(xiàn)相關(guān)性，而CSF的T2值存在輕微的不匹配[12]。Badve等[36]使用相同的MRF框架評估了56例健康志愿者中不同腦區(qū)的T1和T2測量值。 在這項研究中，每層的采集時間為31 s，分辨率為1.2 mm×1.2 mm×5.0 mm。與常規(guī)mapping相比，MRF衍生的弛豫測量同時提供了T1和T2值的快速估計，并證明了區(qū)域差異以及與年齡和性別的相關(guān)性。

MRF還被用于表征和區(qū)分各種腦腫瘤。在一項對31例腦腫瘤患者的研究中，MRF衍生的T1和T2值用于鑒別膠質(zhì)瘤和轉(zhuǎn)移灶。與轉(zhuǎn)移灶相比，低級別膠質(zhì)瘤的平均T2值顯著不同[(172±53) vs. (105±27) ms，P=0.004)[37]。最近正在進(jìn)行的關(guān)于腦腫瘤的MRF研究是針對使用容積3D MRF序列來改善病變可視化、評估治療反應(yīng)以及識別放射性壞死和腫瘤復(fù)發(fā)[38]。

2.前列腺

在已發(fā)表的關(guān)于MRF應(yīng)用于前列腺的研究中，評估了使用MRF-自由進(jìn)動穩(wěn)態(tài)序列(steady state free precession,SSFP)組合和ADC映射檢查的可行性，用于將癌癥與前列腺外周區(qū)域的正常前列腺組織分離。將MRF-SSFP序列添加到常規(guī)前列腺MRI掃描方案中，并使用常規(guī)T2WI作為參考，在MRF生成的T1和T2圖以及DWI上繪制ROI。MRF采集時間為每層39 s，分辨率為1 mm×1 mm×5 mm[39]。癌癥的平均T1、T2和ADC值顯著低于正常周邊區(qū)域[(1628±344) vs. (2247±450)ms，(73±27) vs. (169±61)ms，(0.773×10-3±0.331×10-3) vs. (1.711×10-3±0.269×10-3)mm2/s；P均<0.0001]，三者聯(lián)合應(yīng)用獲得的AUC達(dá)0.99。重要的是，ADC和T2聯(lián)合應(yīng)用鑒別高中級別與低級別前列腺癌的AUC達(dá)0.83。該研究表明，MRF具有識別高級別前列腺癌的潛力。最近正在進(jìn)行的關(guān)于前列腺MRF的應(yīng)用包括對T1和T2值與腫瘤侵襲性的關(guān)系進(jìn)行相關(guān)性分析，并指導(dǎo)腫瘤的靶向活檢。

3.腹部

采用MRF進(jìn)行腹部成像具有獨特的挑戰(zhàn)性，并且需要修改原始的MRF序列。其中包括開發(fā)快速成像序列，實現(xiàn)在一次屏氣中采集所需的數(shù)據(jù)，提高空間分辨率，覆蓋腹部內(nèi)的各種組織，并補(bǔ)償B0和B1不均勻性。目前已經(jīng)提出了幾種使用MRF框架對腹部進(jìn)行成像的方法[40]。陳等[14]采用了一種方法，通過基于Bloch-Seigert方法的單獨掃描測量B1變化，并將相應(yīng)的值納入字典中，以減少B1不均勻性對組織特性的影響。另外，針對B0不均勻性的魯棒性，使用基于SSFP的MRF采集序列。該序列在無癥狀受試者的腹部MRI檢查中得到了驗證，獲得的各種腹部器官的T1和T2值與文獻(xiàn)一致[14]。該序列的臨床應(yīng)用也在6例肝臟轉(zhuǎn)移性腺癌患者(20處病變)中進(jìn)行了測試。轉(zhuǎn)移瘤的平均T1和T2值分別為(1673±331)和(43±13) ms，與周圍肝實質(zhì)[分別為(840±113)和(28±3) ms，P<0.0001和P<0.01]和對照組肝實質(zhì)[(745±65)和(31±6) ms，P<0.0001和P=0.021]顯著不同。此次采集的總掃描時間(包括B1-mapping的時間)為每層19 s，分辨率為1.9 mm×1.9 mm×5.0 mm。MRF框架減少了B1不均勻性對估計組織特性的影響。

4.肌肉骨骼

MRI是肌肉骨骼成像中的重要一環(huán)，能夠提供卓越的軟組織對比，并區(qū)分肌肉骨骼成像中觀察到的各種組織，即骨髓、肌肉、脂肪、軟骨、肌腱和軟組織。在存在整形外科植入物的情況下成像是MRI的一個具有挑戰(zhàn)性的方面，因為圖像受到RF場變化的影響，并且植入物導(dǎo)致某些區(qū)域信號丟失?？寺1不均勻性的MRF方法，如Cloos等[40]提出的方法在存在整形外科植入物的情況下成功獲得了沒有偽影的組織圖。生成這些圖譜的分辨率為1.4 mm×1.4 mm×5.0 mm，18個層面，生成時間約為8 min。

5.心臟

心肌組織特性映射已經(jīng)成為心臟成像的焦點，因為它可以比傳統(tǒng)的心臟成像更早地檢測出病理學(xué)改變。Hamilton等[41]描述的心臟MRF使用16個心動周期的屏氣獲取，同時估計T1、T2值和M0值，分辨率為1.6 mm×1.6 mm×8.0 mm。在健康志愿者中比較MRF與傳統(tǒng)的mapping方法獲得的T1和T2值，顯示出良好的一致性。進(jìn)一步的技術(shù)發(fā)展涉及減少完成整個心臟范圍容積采集的掃描時間并優(yōu)化所獲得的M0值。

6.微血管結(jié)構(gòu)

MRF框架也已擴(kuò)展到表征微血管特性，稱為MR血管指紋圖譜(MRvF)。與磁共振腦氧代謝成像不同[42]，Christen等[43]描述了一種新的MRF方法來測量微血管特性，如腦血容量(CBV)、平均血管半徑(R)和血氧飽和度(SO2)。采集序列為基于自由感應(yīng)衰減和自旋回波序列的梯度回波序列。在注射鐵基對比劑(ferumoxytol)之前和之后的2分鐘進(jìn)行采集。在對比劑注射之前和之后使用不同的數(shù)學(xué)模型模擬二項式，并且使用對比劑注射前、后的對比度信號演變的比率作為指紋。通過使用最小二乘法來確定觀察到的指紋與字典之間的最佳擬合，以生成基于體素的CBV、R和SO2圖。使用1.6 mm×1.6 mm×1.5 mm的分辨率進(jìn)行鐵對比劑注射前、后的MRI掃描，掃描時間為4 min。初步結(jié)果顯示灰質(zhì)和白質(zhì)在CBV圖上具有較高的對比度，表明灰質(zhì)和白質(zhì)量存在血流灌注差異，而R和SO2圖上上述表現(xiàn)更均勻。其它相關(guān)研究進(jìn)一步建立在此基礎(chǔ)上并在大鼠研究中應(yīng)用MRvF以獲得更真實的血管圖或評估中風(fēng)和腦腫瘤中微血管特性的差異[44]。 Lemasson等[44]測試了115只大鼠的MRvF，分為3種模型，即腦腫瘤、中風(fēng)和健康動物，同時獲得微血管結(jié)構(gòu)(即血容量分?jǐn)?shù)，血管直徑)和功能(血氧)圖譜。將這些發(fā)現(xiàn)與常規(guī)MRI方法和組織病理學(xué)結(jié)果進(jìn)行比較以驗證。 MRvF可以在腫瘤和中風(fēng)模型有效地區(qū)分具有不同行為的健康和病理性腦組織。這種微觀特性的體內(nèi)評估，如果有效地轉(zhuǎn)化在臨床應(yīng)用中，有望改善對腦疾病的診斷和管理。

7.血流灌注分析

ASL是一種前景廣闊的不使用對比劑的MR灌注成像技術(shù)。由于ASL信號受到多個血流動力學(xué)參數(shù)的影響，Su等[45]最近探討了基于MRF框架的ASL技術(shù)的可行性，一次掃描可獲得多個血流動力學(xué)參數(shù)。采用EPI讀數(shù)進(jìn)行MRF-ASL采集，全腦覆蓋掃描時間為3 min，分辨率為2.8 mm×2.8 mm×10.0 mm。根據(jù)所使用的數(shù)學(xué)模型提取不同的血液動力學(xué)參數(shù)并與字典值進(jìn)行比較。該序列在志愿者中使用高碳酸血癥激發(fā)驗證，顯示出CO2分壓增好導(dǎo)致的血流灌注增加，以及在Moya-Moya疾病患者中在具有狹窄的頸內(nèi)動脈的一側(cè)進(jìn)行對比劑推注后對比劑到達(dá)時間延遲。

小結(jié)與展望

磁共振指紋識別是MRI的一種新技術(shù)，其脈沖序列設(shè)計不旨在獲取圖像，而是直接測量組織特性。在MRF中，每個不同的組織產(chǎn)生獨特的信號演變，并將其與一組理論信號演變相匹配，以測量單次采集中的若干組織特性，這些特性可以協(xié)同地提供所有信息以幫助診斷，預(yù)后和治療評估。MRF的模式識別特性使得采集對欠采樣和運動等偽影具有魯棒性，這對于將多參數(shù)mapping技術(shù)成功運用到臨床環(huán)境中是至關(guān)重要的。未來，MRF可能會在成像生物標(biāo)志物和新型療法的開發(fā)中發(fā)揮作用。序列設(shè)計、效率、采集速度、后處理步驟和圖像質(zhì)量的改進(jìn)是連續(xù)不斷的過程。MRF的臨床驗證要求建立跨學(xué)科的生物醫(yī)學(xué)工程和臨床團(tuán)隊，這些團(tuán)隊可以同時開發(fā)新的假設(shè)并在各種臨床應(yīng)用和器官系統(tǒng)中測試該技術(shù)。這可以幫助MRF成為有效的臨床工具，并促進(jìn)開發(fā)快速的“一次掃描，多參數(shù)定量”的MR成像方法。