UNICORT
Unified segmentation based correction of R1 brain maps for RF transmit field inhomogeneities (UNICORT)
(Weiskopf et al, 2011) is an algo available in SPM (also SPM BIDS App) to improve SNR of the MP2RAGE acquisition using an additional B1 acquisition a purely data driven approach.
We use this before running freesufer on
- 7T Anatomical Acquisitions
/Volumes/Hera/Projects/corticalmyelin_development/code/corticalmyelin_maturation/image_processing/B1+_transmitfield_correction/unicort
- SPA MR:
/Volumes/Hera/Projects/SPA/pipeline_testing/scripts/mp2rage_processing
Input images should be uncompressed nii files (not .nii.gz)
- "/Volumes/Hera/Projects/SPA/scripts/mri/preproc/spm_unicort.m"
anat = spm_BIDS(bdir,'data', 'modality','anat', 'type','T1w'); % cell of file paths if isempty(anat), error('Cannot find T1-weighted (UNIT1 as T1w) image in', bdir); end % anat = cellfun(@(x)replace(x, '.gz',''), anat, 'Uni', 0) matlabbatch{1}.spm.spatial.preproc.channel.vols = cellstr(anat{1}); % cellstr not needed? matlabbatch{1}.spm.spatial.preproc.channel.biasreg = 0.001; matlabbatch{1}.spm.spatial.preproc.channel.biasfwhm = 60; matlabbatch{1}.spm.spatial.preproc.channel.write = [1, 1]; matlabbatch{1}.spm.spatial.preproc.warp.mrf = 1; matlabbatch{1}.spm.spatial.preproc.warp.cleanup = 1; matlabbatch{1}.spm.spatial.preproc.warp.reg = [0, 0.001, 0.5, 0.05, 0.2]; matlabbatch{1}.spm.spatial.preproc.warp.affreg = 'mni'; matlabbatch{1}.spm.spatial.preproc.warp.fwhm = 0; matlabbatch{1}.spm.spatial.preproc.warp.samp = 3; matlabbatch{1}.spm.spatial.preproc.warp.write = [0, 0];