Dear Gian Marco,
Can't speak for CAT12 but you can do something approaching with SPM
functions. It's not a clean way to do it but it looks visually similar
to your attachment. I just added the support of FreeSurfer annotation
files in @gifti so you would have to use FreeSurfer's read_annotation.m
in the meantime.
Best regards,
Guillaume.
M = gifti('lh.inflated');
P = gifti('lh.Schaefer2018_400Parcels_17Networks_order.annot');
P = P.cdata(:,2);
%[~,P]=read_annotation('lh.Schaefer2018_400Parcels_17Networks_order.annot');
T = randn(size(M.vertices,1),1);
T = spm_mesh_smooth(M,T,1000);
for i=unique(P)'
T(P==i) = mean(T(P==i));
end
H = spm_mesh_render('Disp',M);
H = spm_mesh_render('Overlay',H,T);
hold on % tight
for i=unique(P)'
C = spm_mesh_isoline(M,P==i,0.5);
for j=1:numel(C)
plot3(C(j).xdata,C(j).ydata,C(j).zdata,...
'-','Color',.2*[1 1 1],'LineWidth',2);
end
end
On 26/02/2021 10:10, Gian Marco Duma wrote:
> Dear everyone, I have a question about plotting the value of a ROI with CAT12. I computed the segmentation with CAT12, I extraced the ROI value, and I did some analysis with a inlab script. My results is basically a vaule for each ROI, I'm using Schaefer 100 surface ROIs. I would like to plot the value of each ROI on the surface, to obtain something like the image that I attached. I think that it could be possible using CAT12 but I am not sure. Can you suggest me some possible solutions to do it?
> Thanks
>
>
>
--
Guillaume Flandin, PhD
Wellcome Centre for Human Neuroimaging
UCL Queen Square Institute of Neurology
London WC1N 3BG
|