Resting-statefMRI:

DataProcessing嚴超贛Ph.D.Candidateycg.

?2008SKLCNL嚴超贛的靜息態數據處理DataProcessingOutline

預處理

ALFF,ReHo,FC計算統計結果呈現文獻管理嚴超贛的靜息態數據處理DataProcessing預處理

數據整理

DICOM->NIFTI

去前10個時間點

SliceTimingRealignNormalizeSmooth

去線性漂移

Filter:0.01-0.08嚴超贛的靜息態數據處理DataProcessing數據整理被試信息整理嚴超贛的靜息態數據處理DataProcessing數據整理被試信息整理嚴超贛的靜息態數據處理DataProcessing數據整理被試信息整理原始數據整理靜息功能數據結構數據DTI數據嚴超贛的靜息態數據處理DataProcessing數據整理被試信息整理原始數據整理靜息功能數據結構數據DTI數據嚴超贛的靜息態數據處理DataProcessing預處理

數據整理

DICOM->NIFTI

去前10個時間點

SliceTimingRealignNormalizeSmooth

去線性漂移

Filter:0.01-0.08嚴超贛的靜息態數據處理DataProcessingDICOM->NIFTIMRIcroN’sdcm2niiguiSPM5’sDICOMImport嚴超贛的靜息態數據處理DataProcessing預處理

數據整理

DICOM->NIFTI

去前10個時間點

SliceTimingRealignNormalizeSmooth

去線性漂移

Filter:0.01-0.08嚴超贛的靜息態數據處理DataProcessing預處理

數據整理

DICOM->NIFTI

去前10個時間點

SliceTimingRealignNormalizeSmooth

去線性漂移

Filter:0.01-0.08嚴超贛的靜息態數據處理DataProcessingSliceTimingWhy?嚴超贛的靜息態數據處理DataProcessingSliceTimingWhy?Huetteletal.,2004嚴超贛的靜息態數據處理DataProcessingSliceTiming1:2:33,2:2:323322-(2/33)33嚴超贛的靜息態數據處理DataProcessingSliceTiming嚴超贛的靜息態數據處理DataProcessing預處理

數據整理

DICOM->NIFTI

去前10個時間點

SliceTimingRealignNormalizeSmooth

去線性漂移

Filter:0.01-0.08嚴超贛的靜息態數據處理DataProcessingRealignWhy?嚴超贛的靜息態數據處理DataProcessingRealign嚴超贛的靜息態數據處理DataProcessingRealign檢查頭動?b=load(‘rp_name.txt’);c=max(abs(b));c(4:6)=c(4:6)*180/pi;嚴超贛的靜息態數據處理DataProcessing預處理

數據整理

DICOM->NIFTI

去前10個時間點

SliceTimingRealignNormalizeSmooth

去線性漂移

Filter:0.01-0.08嚴超贛的靜息態數據處理DataProcessingNormalizeWhy?Huetteletal.,2004嚴超贛的靜息態數據處理DataProcessingmean_name.imgr*.imgEPI.nii-90-126-72;9090108333嚴超贛的靜息態數據處理DataProcessingNormalize嚴超贛的靜息態數據處理DataProcessing預處理

數據整理

DICOM->NIFTI

去前10個時間點

SliceTimingRealignNormalizeSmooth

去線性漂移

Filter:0.01-0.08嚴超贛的靜息態數據處理DataProcessingSmoothWhy?減少被試間的空間結構差異；有利統計檢驗的有效性;嚴超贛的靜息態數據處理DataProcessingw*.img平滑核Smooth嚴超贛的靜息態數據處理DataProcessing預處理

數據整理

DICOM->NIFTI

去前10個時間點

SliceTimingRealignNormalizeSmooth

去線性漂移

Filter:0.01-0.08嚴超贛的靜息態數據處理DataProcessing去線性漂移嚴超贛的靜息態數據處理DataProcessing預處理

數據整理

DICOM->NIFTI

去前10個時間點

SliceTimingRealignNormalizeSmooth

去線性漂移

Filter:0.01-0.08嚴超贛的靜息態數據處理DataProcessing濾波Why?關心低頻的活動嚴超贛的靜息態數據處理DataProcessing濾波嚴超贛的靜息態數據處理DataProcessingOutline

預處理

ALFF,ReHo,FC計算統計結果呈現文獻管理嚴超贛的靜息態數據處理DataProcessingALFFZangetal.,2007嚴超贛的靜息態數據處理DataProcessingALFF嚴超贛的靜息態數據處理DataProcessingReHo注意：預處理時不做平滑，計算完ReHo后再進行平滑。Zangetal.,2004嚴超贛的靜息態數據處理DataProcessingReHo嚴超贛的靜息態數據處理DataProcessing功能連接提取協變量

頭動參數:rp_name.txt

全腦信號、腦脊液、白質信號嚴超贛的靜息態數據處理DataProcessing提取協變量嚴超贛的靜息態數據處理DataProcessing提取協變量嚴超贛的靜息態數據處理DataProcessing提取協變量嚴超贛的靜息態數據處理DataProcessing提取協變量嚴超贛的靜息態數據處理DataProcessing提取協變量嚴超贛的靜息態數據處理DataProcessing提取協變量嚴超贛的靜息態數據處理DataProcessing提取協變量嚴超贛的靜息態數據處理DataProcessing功能連接提取協變量

頭動參數:rp_name.txt

全腦信號、腦脊液、白質信號協變量組合：RPCov=load('rp_name.txt');BCWCov=load('ROI_FCMap_name.txt');Cov=[RPCov,BCWCov];save('Cov.txt','Cov','-ASCII','-DOUBLE','-TABS');嚴超贛的靜息態數據處理DataProcessing功能連接Voxel-wiseROI-wiser=0.36嚴超贛的靜息態數據處理DataProcessingVoxel-wise嚴超贛的靜息態數據處理DataProcessingVoxel-wise嚴超贛的靜息態數據處理DataProcessingVoxel-wise嚴超贛的靜息態數據處理DataProcessingVoxel-wise嚴超贛的靜息態數據處理DataProcessingVoxel-wise嚴超贛的靜息態數據處理DataProcessingVoxel-wise嚴超贛的靜息態數據處理DataProcessingVoxel-wise1by1Cov.txt嚴超贛的靜息態數據處理DataProcessingROI-wise嚴超贛的靜息態數據處理DataProcessingROI-wise嚴超贛的靜息態數據處理DataProcessingROI-wise1by1Cov.txt嚴超贛的靜息態數據處理DataProcessingOutline

預處理

ALFF,ReHo,FC計算統計結果呈現文獻管理嚴超贛的靜息態數據處理DataProcessing統計單樣本T檢驗雙樣本T檢驗嚴超贛的靜息態數據處理DataProcessing單樣本T檢驗統計準備

ALFF:mALFFReHo:smReHoFC:zFC嚴超贛的靜息態數據處理DataProcessing單樣本T檢驗mALFF-1,smReHo-1m圖像1by1i1-1嚴超贛的靜息態數據處理DataProcessing單樣本T檢驗m*-1圖像;zFC圖像存放SPM.mat的目錄嚴超贛的靜息態數據處理DataProcessing單樣本T檢驗上一步生成的spm.mat嚴超贛的靜息態數據處理DataProcessing單樣本T檢驗上一步生成的spm.mat1嚴超贛的靜息態數據處理DataProcessing單樣本T檢驗HeightthresholdT=2.687659{p<0.05(FDR)}Extentthresholdk=10voxelsspmT_0001.imgPD_1T.img準備制作雙樣本T檢驗的MASK嚴超贛的靜息態數據處理DataProcessing雙樣本T檢驗統計準備

ALFF：mALFFReHo:smReHoFC:zFCMask嚴超贛的靜息態數據處理DataProcessing雙樣本T檢驗MaskMask=PD_1T.img+CON_1T.img>0PD_1T.imgCON_1T.imgMask_2T.imgi1+i2>0嚴超贛的靜息態數據處理DataProcessing雙樣本T檢驗兩組的m*圖像;zFC圖像存放SPM.mat的目錄Mask_2T.img嚴超贛的靜息態數據處理DataProcessing雙樣本T檢驗上一步生成的spm.mat嚴超贛的靜息態數據處理DataProcessing雙樣本T檢驗上一步生成的spm.mat1-1嚴超贛的靜息態數據處理DataProcessing雙樣本T檢驗HeightthresholdT=Inf{p<0.05(FDR)}Extentthresholdk=0voxelsspmT_0001.img嚴超贛的靜息態數據處理DataProcessingOutline

預處理

ALFF,ReHo,FC計算統計結果呈現文獻管理嚴超贛的靜息態數據處理DataProcessing結果呈現SPM5xjviewrest_sliceviewerMRIcroN嚴超贛的靜息態數據處理DataProcessing結果呈現—SPM5嚴超贛的靜息態數據處理DataProcessing結果呈現—SPM5嚴超贛的靜息態數據處理DataProcessing結果呈現—xjview嚴超贛的靜息態數據處理DataProcessing結果呈現—xjview嚴超贛的靜息態數據處理DataProcessing結果呈現—rest_sliceviewer嚴超贛的靜息態數據處理DataProcessing結果呈現—MRIcroN嚴超贛的靜息態數據處理DataProcessing結果呈現—MRIcroN嚴超贛的靜息態數據處理DataProcessingOutline

預處理

ALFF,ReHo,FC計算統計結果呈現文獻管理嚴超贛的靜息態數據處理DataProcessing文獻管理嚴超贛的靜息態數據處理DataProcessing文獻管理嚴超贛的靜息態數據處理DataProcessing文獻管理嚴超贛的靜息態數據