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

Submitted by yubing on Mon, 12/23/2013 - 10:19

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

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

  

Submitted by admin on Tue, 12/24/2013 - 11:29

Check the following codes, it is easy to get each voxel's time course inside a ROI mask.
 http://restfmri.net/pub/rest_20090422/doc/rest_20090422/fc.html

0146         if ~IsDefinedROITimeCourse,% I need retrieving the ROI averaged time course manualy
0147             maskROI =find(maskROI);
0148             % [rangeX, rangeY, rangeZ] = ind2sub(size(maskROI), find(maskROI));
0149             % theTimeCourses =zeros(length(unique(rangeX)), length(unique(rangeY)), length(unique(rangeZ)));
0150             for t=1:sampleLength,
0151                 theTimePoint = squeeze(AllVolume(:,:,:, t));
0152                 theTimePoint = theTimePoint(maskROI);
0153                 AROITimeCourse(t) =mean(theTimePoint);
0154             end    %The Averaged Time Course within the ROI now comes out! 20070903
0155             %Make sure the ROI averaged time course is an col vector
0156             AROITimeCourse =reshape(AROITimeCourse, sampleLength,1);
0157         end
0158         %Save the ROI averaged time course to disk for further study
0159         [pathstr, name, ext, versn] = fileparts(AResultFilename);
0160         theROITimeCourseLogfile =[fullfile(pathstr,['ROI_', name]), '.txt'];
0161         save(theROITimeCourseLogfile, 'AROITimeCourse', '-ASCII', '-DOUBLE','-TABS')