BG
在很很多时候,在对FMRI建立网络的时候,很多时候都是需要预先定义好ROI然后计算功能连接,所以在提取ROI的时候,就需要对模版图像进行处理。
Workflow
- 1 利用Dpabi的image reslice 功能,对模版进行重新采样。
- 2 利用MASK对大脑进行分区时间序列的提取,FMRI是4D数据,笨法子就一层一层按ROI进行提取,然后在合并,但是这种效率比较读,相对而言,进行Reshape就比较快,只需要将前三维展开成一维的 ,最后成2维的数据即可进行提取,具体的matlab代码如下。以power 264的alta为例子。
fmridata=y_ReadAll(combineName);##读取FMRI data
[a,b,c,d]=size(fmridata);
if(d==490)##统一时间序列长度
fmridataTran=reshape(fmridata,[],size(fmridata,4));
slide264=zeros(264,490);
for j=1:264
% j=1;
roi1=find(atlasdata==j);
slide264(j,:)=mean(fmridataTran(roi1,:));
end