MRI-批处理-roi提取

为了实现可以全部批量处理，编写了一个循环，个人笔记，不一定适用，欢迎探讨

clear
clc
raw_path='D:\2023_12_NSSI_YLL\2024_01_15\results_nsa_nsi_hc\smri\cat_nsa43_nsi28_nsi38\'%原始数据路径
[date,name,biao]=xlsread('D:\2023_12_NSSI_YLL\2024_01_15\results_nsa_nsi_hc\smri\nsa43_nsi28_hc38.xlsx',1);%表格在此处无用
cell_data=dir([raw_path,'*_*'])
for i=1:length(cell_data)%提取数据所在路径及名称，并且全部放在元胞数组中
    sub_num=cell_data(i).name
    sub_data_path=fullfile(raw_path,sub_num,'label\','*catROI*.xml')
    nift_sub=dir(sub_data_path)
    sub_paname=fullfile(raw_path,sub_num,'label\',nift_sub(2).name)
    sub_aim(i,1)={sub_paname}
end
%cat12中提取代码
% save('gmv_roi_get_file.mat', 'sub_aim');
matlabbatch{1}.spm.tools.cat.tools.calcroi.roi_xml = sub_aim;%将输入数据的xml替换为已获得的元胞数组
matlabbatch{1}.spm.tools.cat.tools.calcroi.point = '.';
matlabbatch{1}.spm.tools.cat.tools.calcroi.outdir = {'D:\2023_12_NSSI_YLL\2024_01_15\results_nsa_nsi_hc\smri\gmv\male_female_nsa43_nsi28_hc38'};
matlabbatch{1}.spm.tools.cat.tools.calcroi.calcroi_name = 'sbm_roi';
spm_jobman('run',matlabbatch)