批量生成种子点

for j=1:246
    matlabbatch{1}.spm.util.imcalc.input = {'D:\brain_map\brainnetome_246_3mm\BN_Atlas_246_3mm.nii,1'};
    matlabbatch{1}.spm.util.imcalc.output = ['seed_',num2str(j)];
    matlabbatch{1}.spm.util.imcalc.outdir = {'D:\brain_map\brainnetome_246_3mm\seed_246'};
    matlabbatch{1}.spm.util.imcalc.expression = ['i1==',num2str(j)];
    matlabbatch{1}.spm.util.imcalc.var = struct('name', {}, 'value', {});
    matlabbatch{1}.spm.util.imcalc.options.dmtx = 0;
    matlabbatch{1}.spm.util.imcalc.options.mask = 0;
    matlabbatch{1}.spm.util.imcalc.options.interp = 1;
    matlabbatch{1}.spm.util.imcalc.options.dtype = 4;
    spm_jobman('run',matlabbatch)
end

%为了批量导入dpabi，生成mat文件
clear
clc
raw_path='D:\brain_map\brainnetome_246_3mm\tian_4\'
% [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)
    %nift_sub=dir(sub_data_path)
   % sub_paname=fullfile(raw_path,sub_num,'surf\',nift_sub(1).name)
       ROICell(i,1)={sub_paname}
end
save('ROI_List_44.mat','ROICell')