Dear Ben,
thank you for looking into this problem.
I used "coregister: reslice" with nearest neighbour to reslice to the VBM output.
In detail this is the jobfile:
%-----------------------------------------------------------------------
% Job configuration created by cfg_util (rev $Rev: 4252 $)
%-----------------------------------------------------------------------
matlabbatch{1}.spm.spatial.coreg.write.ref = {'C:\...\spmT_0002.img,1'};
matlabbatch{1}.spm.spatial.coreg.write.source = {'C:\...\ROI_MNI_V4_Cingulum_Ant_LR.nii,1'
'C:\...\WFU_AAL_ACC_LR.nii,1'
};
matlabbatch{1}.spm.spatial.coreg.write.roptions.interp = 0;
matlabbatch{1}.spm.spatial.coreg.write.roptions.wrap = [0 0 0];
matlabbatch{1}.spm.spatial.coreg.write.roptions.mask = 0;
matlabbatch{1}.spm.spatial.coreg.write.roptions.prefix = 'r_nn_';
Then I compared the two resliced masks:
M_AAL=spm_vol('r_nn_ROI_MNI_V4_Cingulum_Ant_LR.nii')
Vol_AAL=spm_read_vols(M_AAL);
M_WFU=spm_vol('r_nn_WFU_AAL_ACC_LR.nii')
Vol_WFU=spm_read_vols(M_WFU);
sum(Vol_AAL(:))
ans: 6328
sum(Vol_WFU(:))
ans: 6439
and to check the overlap:
sum_d_AAL_WFU=sum(abs(d_AAL_WFU(:)))
ans: 379
Best wishes,
Katharina
|