Dear Thomas,
this is not supported. CAT12. You could try to tweak the code in cat_surf_surf2roi.m that extract values inside a surface ROI. Here, functions such as min,max or mean are used to summarize the values inside the ROI and you could try to change the code and define FN{ai} as empty and also change the resulting variable and output:
for roii=1:numel(catROI{si}.(rinfo.dataname).ids)
switch catROI{si}.(rinfo.dataname).names{roii}(1)
case 'l', catROI{si}.(rinfo.dataname).data.(fieldname)(roii,1) = ...
eval(sprintf('%s%s(lCS.cdata(lrdata==catROI{si}.(rinfo.dataname).ids(roii)))',nanfunc,FN{ai}));
case 'r', catROI{si}.(rinfo.dataname).data.(fieldname)(roii,1) = ...
eval(sprintf('%s%s(rCS.cdata(rrdata==catROI{si}.(rinfo.dataname).ids(roii)))',nanfunc,FN{ai}));
case 'b', catROI{si}.(rinfo.dataname).data.(fieldname)(roii,1) = ...
eval(sprintf(['%s%s(lCS.cdata(lrdata==catROI{si}.(rinfo.dataname).ids(roii))) + ' ...
'%s%s(rCS.cdata(rrdata==catROI{si}.(rinfo.dataname).ids(roii)))'],nanfunc,FN{ai},nanfunc,FN{ai}));
otherwise, catROI{si}.(rinfo.dataname).data.(fieldname)(roii,1) = nan;
end
end
Sorry, for having no easier suggestion.
Best,
Christian
On Tue, 21 Sep 2021 16:24:43 +0100, Thomas Stephan <[log in to unmask]> wrote:
>Dear Dr. Gaser,
>
>inspired by https://www.nature.com/articles/s41598-020-62832-z we performed a surface based analysis of an existing fMRI data set. We can now look at second level results on the surface using cat_surf_results, overlay e.g. the HCP parcellation, and use 'atlas labeling' to get the overlap of the results with the HCP labels.
>We would now like to extract the results-values of all surface elements inside specific atlas labels, for later comparison between hemispheres. Is this possible with the tools in CAT12, or will I need to modify the code in cat_surf_results.m - select_atlas(atlas) to tackle this challenge?
>
>Best regards, and many thanks for this great toolbox!
>Thomas Stephan
|