第一期|腦功能MRI數據處理手冊系列教程

2020-12-04 行上行下

作者按

我希望我的投稿能給所有從事腦功能研究的同學、同行、老師帶來幫助,能讓他們在山重水複疑無路時有所借鑑。本文章是原創文章,如需轉載,請聯繫原作者。

作者:範嘉晨

指導老師:劉加成

單位:東南大學 江蘇省分子影像與功能影像重點實驗室

轉載自:我愛腦科學網

Part1

數據處理手冊介紹

本手冊主要詳細講解如何進行頭顱磁共振數據的影像處理及相關步驟的原理。主要處理任務態磁共振(Stroop、N-back、facialrecognition)、靜息態磁共振、彌散張量成像、3D結構像。

本手冊使用以下軟體進行圖像處理(建議下載最新版):MRIConvert、MRIcron(自帶dcm2nii)、MATLAB(R2015b)、SPM8、VBM8、xjview96、FSL、FreeSurfer。

以上軟體除了MATLAB是付費軟體外(東南大學有購買正版版權,可以利用東大學生一卡通號進行註冊激活),其餘皆可在官網下載。

需注意FSL、FreeSurfer目前主要在非Windows系統下使用(Windows下安裝子系統Ubuntu或安裝虛擬機也可以運行,但是不穩定,具體參考FSL官網安裝步驟)。

Part2

相關軟體的介紹

1.圖像格式轉換軟體:MRIConvert、MRIcron

1.1 下面講一下MRIConvert使用方法,MRIcron類似,首先點擊Input下Add files或者Add folders,添加文件或文件夾;也可以在DICOM series toconvert 右側直接點擊Add files或者Addfolders。

在Output下選擇Directory選擇保存目錄即保存位置,options可以設置保存文件的許多屬性(比如名字);也可以不點擊Output,直接在首頁點擊Options或Directory。注意兩白色框之間的文件格式選擇(比如勾選FSL NifTI),最後點擊Convert all(或者Convert-Convert All)。

1.2 下面講一下MRIcron的用法,官網下載mricron,解壓後雙擊裡面的dcm2niigui.exe,打開其圖形界面。可以通過File下不同轉換方式(如DICOM to NIfTI)添加目標文件,Help下Preference可以選擇輸出文件的屬性(如命名),Output Format可以選擇輸出文件的格式。當然也可以在linux的終端使用dcm2nii的命令來轉換圖像,具體命令如下:

2. 統計參數圖SPM

SPM是倫敦大學學院基於MATLAB開發的功能磁共振成像分析軟體包(當然其也可以處理PET,EEG等),可以從官網免費下載,但是要注意必須在MATLAB中設置路徑裡添加SPM這個軟體的文件夾,然後在MATLAB裡輸入spm(小寫),點擊fMRI後才能打開SPM這個軟體。

3.可視化工具xjview

官網下載此軟體後,同樣需加入MATLAB設置路徑(set path)裡,終端輸入xjview即可打開。

4.基於體素的形態學分析包VBM

此軟體可以用來分析大腦結構相(灰質、白質、腦脊液),官網下載此軟體後,需將其複製到SPM文件夾下toolbox中,隨SPM添加到MATLAB的設置路徑裡。此時在SPM8界面中Toolbox下拉即可出現VBM選項。注意VBM版本必須與SPM版本一致才能兼容。

5. FSL

FMRIB創建於牛津大學,可以分析靜息態、任務態、動脈自旋標記成像、結構像、彌散張量成像等影像數據。但是FSL目前在Windows系統下不是很兼容,建議在UNIX平臺下使用。本文主要在centos 7上面使用FSL。

請注意:首先一定要正確安裝FSL,在bash裡設置正確的環境變量,否則在centos7終端輸入fsl,不會喚起fsl圖形界面(這個問題筆者遇到過很多次,在windows虛擬機上安裝時、在windows的子系統Ubuntu上安裝時、多次在centos7上安裝時都遇到過,網上也有許多人遇到這個問題,所以一定要重視);其次,利用FSL官網提供的測試數據(非練習數據)自動運行一遍,檢測安裝正確與否。

6.FreeSurfer

FreeSurfer是一個解剖磁共振分析軟體包,它本身不是功能磁共振成像分析軟體包。因為它提供了用最少的人力投入自動生成皮質表面模型和解剖分割的手段,所以它對於功能磁共振成像數據分析越來越有用。它可以基於表面配準的方法對齊被試間的數據,該方法在某些情況下比標準的基於體積配準的被試間對齊更精確;它還可以把使用FSL或SPM獲得的統計結果投射到重建的皮質表面上,並允許基於表面的組間統計分析。

請注意:一定要正確安裝FresSurfer(本人在windows/Ubuntu/centos上都安裝成功過),並正確設置其環境變量,否則打不開其圖形界面;其次一定要按照官網安裝步驟裡提到的example運行一遍,檢測所有環節有無紕漏。

7.主要fMRI分析軟體包

8.處理流程概述

Part3

圖像格式的介紹

本實驗主要涉及DICOM、NIfTI格式。在從PHILIPS Ingenia 3.0T 磁共振機子上拷貝數據時注意:

1.必須要拷貝所有圖像的DICOM格式,尤其是DTI數據拷貝DICOM時,split data選項絕對不要勾(防止DTI數據轉換成NIfTI時b0與b1000的圖像分開)。

2.除了所有圖像的DICOM格式需要外,非DICOM格式如NIfTI格式及FSL NIfTI格式的數據也需要。(雖然可以使用圖像轉換軟體將DICOM格式轉換為NIfTI格式,但是部分機型磁共振可能不兼容,影響轉換結果)。

3.數據記得傳PACS以及多方式備份,拷貝數據時一定要核對完整圖片數量(如靜息態數據採集240個時間點,每個時間點有33層,則圖片總計240*33=7920張;任務態數據採集265個時間點,每個時間點有33層,則圖片總計265*33=8745張。

1.DICOM

DICOM保留著在標頭中最大數量的元數據,包括掃描儀和圖像採集的基本信息及被試的相關信息。雖然DICOM時MRI掃描儀輸出數據的標準格式,但在數據分析之前,我們通常需要把DICOM轉換為其他格式。主要原因是,DICOM數據集是很難處理的,因為每個部分被存儲為一個單獨的文件。這可能會很快產生大量的阻礙文件系統運行和減慢分析速度的小文件。

2.Analyze

Analyze將每個數據集存儲在一組兩個文件裡。一個擴展名為.img的數據文件,包含了圖像的二進位數據。另一個擴展名為.hdr的頭文件,包含了圖像的元數據。其主要的缺陷是頭文件中僅有相當有限的圖像元數據。

3.NIfTI

NIfTI格式其中一個最重要的特點是它可以呈現三維像素指數和MRI掃描儀的空間位置關係。如果使用得當,這有助於確保人們總是可以準確地確定方向(如確定是神經學方向還是放射學方向)。

NIfTI圖像標準的文件擴展名是.nii,其中包括標題和圖像數據。它也可以利用單獨的圖像(.img)和頭文件(.hdr)代表NIfTI圖像。單個文件.nii格式的一個方便的特點是文件可以被標準的壓縮軟體(如gzip)壓縮,並且一些軟體包如FSL可以直接讀寫壓縮的nii文件(擴展名為.nii.gz)。

今天就主要講解到到這裡了,分享完畢,希望有所幫助。下期我們繼續分享,歡迎持續關注喲~

作者聯繫郵箱:jia_chen_fan@163.com

相關焦點

  • 行上行下|2020年腦科學技術乾貨合集
    ERPLAB教程 | 畫漂亮的腦電圖14. EMSE軟體使用手冊15. 乾貨 | 腦電-超掃描數據分析16. 實現EEG功能連接的可視化17. 人類大腦連接的網絡可視化工具18. 使用Cartool工具包分析EEG源成像?19. HyPyP包:超掃描可視化python包20.
  • The R-fMRI Course 功能核磁共振中文視頻課程
    抑鬱症腦成像大數據項目介紹                 三、靜息態功能磁共振成像在心理學和腦疾病的應用                 四、靜息態功能磁共振成像數據處理原理與DPARSF使用                 五、基於皮層的腦成像數據分析與DPABISurf
  • 大師兄的10本醫學影像技術類書籍 (含MRI/PET-MRI/腦科學等)
    50本醫療大數據/人工智慧•經典書籍免費贈送!1、醫用影像設備(CT/MR/DSA)成像原理與臨床應用(全國醫用設備使用人員業務能力考評教材)2、醫學影像學 (人衛版 第9版)7、fMRI數據分析手冊 8、功能神經影像學
  • 降維算法讓腦神經大數據處理變簡單
    絕大部分感覺、認知和運動發起功能都依賴腦中大量神經元相互作用。
  • 腦功能成像的新方法------功能性近紅外光譜技術(fNIRS)
    今天我要給大家介紹一種新興的腦功能成像技術:功能性近紅外光譜技術 a functional near-infrared spectroscopy(fNIRS)。 功能性近紅外光譜技術是近年來新興的一種非侵入式腦功能成像技術。fNIRS進行腦功能成像的原理與fMRI相似,即大腦神經活動會導致局部的血液動力學變化。
  • 腦與神經類開放資料庫匯總
    下面是根據文章中關於腦與神經類資料庫 部分的內容整理1000 Functional Connectomes Project最為著名的腦與神經資料庫之一http://fcon_1000.projects.nitrc.org/fcpClassic/FcpTable.html ADNI提供與阿爾茨海默疾病相關的各種數據
  • 可預測兒童腦創傷功能結局的MRI生物標誌物
    如何可靠預測這類患者腦創傷後的功能結局是目前迫切需要解決的問題。MRI在腦創傷患者中廣泛使用,並且具有識別病理生理學關鍵變化的巨大潛力,這些變化可用於預測患者的功能結局。由於豬的腦解剖結構、生理學和發育方面與人類的相似性,其已成為人們關注的腦創傷模型動物。
  • C#初學者教程系列4:C 數據類型示例,int、double、string、var
    本文是C#初學者簡單教程,這是第3篇。感謝觀看,記得關注我,後續還有更多教程文章,謝謝。如無特別說明,本系列教程使用的Visual Studio版本為VS2019,.Net Framework版本為4.7.2,.Net Core版本為3.0,.Net Standard為2.1。
  • 「學科交叉的魅力」名師系列講座第三講:饒毅演繹腦的學科交叉研究
    2020年5月24日晚,北京大學前沿交叉學科研究院「學科交叉的魅力」名師系列講座第三講舉行。講座主題為「腦:物理、化學、生物、心理認知的交叉研究」,由北京大學講席教授、北大麥戈文研究所創始所長、北京腦科學中心創始主任、首都醫科大學校長饒毅主講,北京大學博雅教授、研究生院副院長、前沿交叉學科研究院副院長陳鵬主持。講座在學習強國、北京大學央視頻、抖音、快手、微視、B站、百度平臺以及知識分子平臺播出,在線收看人數超過150萬。
  • 這篇文章,助你搞定XRD數據分析,玩轉Jade!
    它具有強大的粉末XRD圖處理和分析功能,尤其注重定量分析和鑑定物相,可實現非常規少量物相的準確鑑定和定量分析。我們本次推薦的Jade視頻有兩套,分別是Jade開發者周榮生博士的Jadex系列視頻、黃繼武老師的Jade分析和精修視頻。
  • 南方CASS測繪軟體,視頻講解+word版教程,限時打包帶走
    今天小嘉和大家分享的是南方CASS測繪軟體系列教程,包括了視頻講解和word版教程,手把手教你快速入門CASS測繪軟體,視頻講解全面細緻,word版教程圖文並茂淺顯易懂,一起來學習吧。【文末領取方式!!】
  • 助力腦醫學科研和臨床應用,慧腦雲獲A輪融資5000萬元
    慧腦雲CEO周亮表示:本輪融資將用於擴大公司規模,加大對多模態腦影像數據分析計算平臺的銷售,採用私有雲、混合雲和公有雲的形式進行市場布局;進一步完善腦結構定量分析平臺、腦磁圖分析平臺、腦小血管病、AD早期診斷、癲癇病灶精準定位產品研發;建立全國性的銷售和技術服務網絡。慧腦雲公司致力於探索腦科學、助力腦醫學的科研和臨床應用。
  • 色彩系列教程:實際運用
    首頁 > 教程 > 關鍵詞 > 色彩最新資訊 > 正文 色彩系列教程:實際運用
  • Keras入門系列教程:兩分鐘構建你的第一個神經網絡模型
    本教程由深度學習中文社區(Studydl.com)持續發布與更新, 本系列其餘教程地址見文章末尾.這是一個用於構建和訓練模型的高級API,其中包括對TensorFlow特定功能的一流支持,例如急切執行, tf.data pipeline和estimators。 tf.keras使TensorFlow更易於使用而不犧牲靈活性和性能。
  • TEM制樣﹑操作﹑數據分析﹑教程書籍資料大全(附下載)
    3.JEM2100培訓手冊1.TEM圖像分析處理軟體DigitalMicreograph,含使用手冊2.DM軟體教程(一)Digitalmicrograph也可以做顆粒尺寸分析(二)DM中文說明書(三)Multivariate Histogram Analysis Users Guide rev11.TEM透射電鏡粒徑分布統計方法2.透射電鏡圖處理.doc3.常見晶體的標準電鏡衍射花樣standard TEM SAED.pdf4.TEM透射電鏡衍射斑點標定深入淺出5.
  • 一種基於人工智慧的腦圖像處理方法
    雖然其中一些方法取得了很好的效果,但在分析高度複雜的腦成像數據時,尤其是從密集信號區域分析時,大多數方法的效果並不理想。由於意識到現有的腦細胞檢測和計數技術的局限性,UZH的團隊開發了一種新的全自動全腦圖像處理技術,稱為DeNeRD。
  • 大數據分析Python NumPy庫使用教程
    在處理數據時,有一個功能強大的庫可以極大地提高代碼的效率,尤其是在處理大型數據集時:NumPy。大數據分析Python NumPy庫使用教程https://www.aaa-cg.com.cn/data/2424.html這為我們的數據工程師之路的下一課程提供了很好的入門:處理Pandas中的大數據集。 在這兩門課程結束時,您將能夠使用Python技能以及NumPy和Pandas的新知識來處理和處理龐大的數據集,這要比普通Python高效得多。
  • python數據科學系列:numpy入門詳細教程
    數據科學基礎庫主要是三劍客:numpy,pandas以及matplotlib,每個庫都集成了大量的方法接口,配合使用功能強大。平時雖然一直在用,也看過很多教程,但紙上得來終覺淺,還是需要自己系統梳理總結才能印象深刻。本篇先從numpy開始,對numpy常用的方法進行思維導圖式梳理,多數方法僅拉單列表,部分接口輔以解釋說明及代碼案例。最後分享了個人關於axis和廣播機制的理解。
  • 西門子1500高速脈衝採集功能和應用及數據的處理
    兩路高速脈衝採集通道,在一般項目應用中有點少 在工業領域數據採集和處理中,經常會遇到如流量、轉速、扭矩等高速脈衝信號形式的傳感器。這種類型的傳感器精確度高,抗幹擾能力強,得到技術人員的青睞。 舉例說明西門子S7-1500高速脈衝採集功能和應用,包括數據的處理。