关于REST的时间序列提取功能

Submitted by yubing on

    REST 提供了很好很强大的ROI时间序列提取(Extract ROI Signal )功能,可以非常方便的提取ROI的数值,我以前的研究从中受益很多。

    最近的研究,要求不仅仅提取出一个ROI的平均时间序列,还要提取这个ROI中每一个Voxel的时间序列,不知需要如何将该程序修改或重写?多谢指教!

Qing Zhao

Mon, 12/23/2013 - 04:12

个人觉得没必要去改动这个函数,不如用rest_readfile读入ROI和全部功能像的数据,直接从矩阵上取ROI内的时间序列就行了。
比如
[ROI,v,header]=rest_readfile('ROI.nii')
[data,v,header]=rest_readfile('4drestdata.nii')
for i=1:n %n :timepoints
X=data(:,:,:,i)
ROItimecourse(:,i)=X(ROI==1);
end;
即可



 
 

<