SPM批量将功能项配准到结构像和MNI

最新推荐文章于 2024-11-11 11:07:54 发布

原创最新推荐文章于 2024-11-11 11:07:54 发布 · 1.8k 阅读

13 ·

CC 4.0 BY-SA版权

文章标签：

#matlab

本文介绍了一个使用SPM12进行正电子发射断层成像(PET)数据标准化处理的MATLAB脚本示例。该脚本通过读取指定路径下的文件，实现了PET图像与模板之间的配准及标准化过程。

petmni.m

global sub filelist

pathname='C:\Users\86191\Desktop\AD数据整理';
filelist=dir(pathname);
filelist=filelist(~ismember({filelist.name},{'.','..'}));

for i=1:length(filelist)
    sub=filelist(i).name;


% List of open inputs
nrun = 1; % enter the number of runs here
jobfile = {'C:\Users\86191\Desktop\2.0ruijinprocess\petmni_job.m'};
jobs = repmat(jobfile, 1, nrun);
inputs = cell(0, nrun);
for crun = 1:nrun
end
spm('defaults', 'PET');
spm_jobman('run', jobs, inputs{:});

end

petmni_job.m

global sub
inputpath=['C:\Users\86191\Desktop\AD数据整理\',sub];

%load raw data
pathname1=[inputpath '/'];
name_t1=dir([pathname1,'*FDGcoT1.nii']);
imgfile_t1=[pathname1 name_t1.name];

name_FDG=dir([pathname1,'*_oFDG.nii']);
imgfile_FDG=[pathname1 name_FDG.name];


%-----------------------------------------------------------------------
% Job saved on 09-Nov-2020 22:10:48 by cfg_util (rev $Rev: 7345 $)
% spm SPM - SPM12 (7487)
% cfg_basicio BasicIO - Unknown
%-----------------------------------------------------------------------
matlabbatch{1}.spm.spatial.normalise.estwrite.subj.vol = {imgfile_t1};
matlabbatch{1}.spm.spatial.normalise.estwrite.subj.resample = {imgfile_FDG};
matlabbatch{1}.spm.spatial.normalise.estwrite.eoptions.biasreg = 0.0001;
matlabbatch{1}.spm.spatial.normalise.estwrite.eoptions.biasfwhm = 60;
matlabbatch{1}.spm.spatial.normalise.estwrite.eoptions.tpm = {'C:\Program Files\MATLAB\R2018a\toolbox\spm12\tpm\TPM.nii'};
matlabbatch{1}.spm.spatial.normalise.estwrite.eoptions.affreg = 'eastern';
matlabbatch{1}.spm.spatial.normalise.estwrite.eoptions.reg = [0 0.001 0.5 0.05 0.2];
matlabbatch{1}.spm.spatial.normalise.estwrite.eoptions.fwhm = 0;
matlabbatch{1}.spm.spatial.normalise.estwrite.eoptions.samp = 3;
matlabbatch{1}.spm.spatial.normalise.estwrite.woptions.bb = [NaN NaN NaN
                                                             NaN NaN NaN];
matlabbatch{1}.spm.spatial.normalise.estwrite.woptions.vox = [2 2 2];
matlabbatch{1}.spm.spatial.normalise.estwrite.woptions.interp = 4;
matlabbatch{1}.spm.spatial.normalise.estwrite.woptions.prefix = 'w';