医学影像批量配准，matlab 代码
配准后的文件与配准前的文件保存在统一目录下，命名以‘r’开头。
%-----------------------------------------------------------------------
% Job saved on 21-Nov-2019 15:35:09 by cfg_util (rev $Rev: 7345 $)
% spm SPM - SPM12 (7487)
% cfg_basicio BasicIO - Unknown
% 用于批量配准
%-----------------------------------------------------------------------
clear
spm('defaults', 'fmri');
spm_jobman('initcfg');
subjectsdir = {'E:\nii_path'};
% 这里是data文件夹的绝对位置
path = 'E:\nii_path';
subjects = dir(path);
subjects = {subjects.name};
% 单个或多个被试的文件夹
subjects(:,[1,2]) = [];
funcdir = fullfile('T1C');
% reference的文件夹，向哪个序列配准，哪个是reference
anatdir = fullfile('T1');
% source的文件夹，被配准的序列
nsubj = length(subjects);
% 配准的病例数
for csubj = 1:nsubj
subjdir = {spm_select('CPath', subjects{csubj}, subjectsdir)};
% source file
adir = spm_select('CPath', anatdir, subjdir);
% 选择当前source文件夹
anatfile = spm_select('FPList', adir, '/*.nii');
% 该文件夹中.nii结尾的文件
if size(anatfile, 1) == 0
sprintf('No anat file found for %s',subjects{csubj})
continue;
end
% reference file
fdir = {spm_select('CPath', funcdir, subjdir)};
% 选择当前reference文件夹
ffiles = spm_select('FPList', fdir, '/*.nii');
% 该文件夹中.nii结尾的文件
if size(ffiles,1) == 0
sprintf('No source file found for %s', subjects{csubj})
continue;
end
% Coregister
n_nii = size(anatfile,1);
% nii文件个数，被配准的序列可能不止一个nii
for n = 1:n_nii
clear matlabbatch
% 以下spm12默认设置
matlabbatch{1}.spm.spatial.coreg.estwrite.ref = {ffiles};
matlabbatch{1}.spm.spatial.coreg.estwrite.source = {anatfile(n,:)};
matlabbatch{1}.spm.spatial.coreg.estwrite.other = {''};
matlabbatch{1}.spm.spatial.coreg.estwrite.roptions.interp = 4;
matlabbatch{1}.spm.spatial.coreg.estwrite.roptions.wrap = [0 0 0];
matlabbatch{1}.spm.spatial.coreg.estwrite.roptions.mask = 0;
matlabbatch{1}.spm.spatial.coreg.estwrite.roptions.prefix = 'r';
% 运行
spm_jobman('run', matlabbatch);
end
【matlab的spm,医学图像的批量配准（Matlab，SPM）,医学影像,matlabSPM】end
