下面先展示一下效果图,只保留了感兴趣区域1,10,15的ROI.
下面以我们熟悉的Glasser 360 模板为例,
只展示感兴趣区域1,10,15的区域,先看一下代码:
clc;
clear all
%让1,10,15为1,其余为0
out = repelem(0, 180);
out(1) = 1;
out(10) = 1;
out(15) = 1;
csvwrite('parcelNums.csv', out')
mmpL = gifti('mmpL.func.gii'); % load the left side gifti MMP atlas
mmpR = gifti('mmpR.func.gii'); % and the right side MMP atlas
newvals = csvread('parcelNums.csv'); % 180 integers; new value for each parcel
Lout = mmpL; % output gifti
Lout.cdata(:,1) = repelem(0, size(mmpL.cdata,1)); % replace the values with zeros
Rout