如何绘制好看的大脑皮层图?wb_view +fun.gii格式操作
下面先展示一下效果图,只保留了感兴趣区域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