% Batch file of CAT12 segmentation for SPM12 standalone installation % modified parameters for ENIGMA % %_______________________________________________________________________ % $Id: cat_standalone_segment_enigma.m 1842 2021-06-01 14:41:58Z gaser $ % first undefined data field, that will be dynamically replaced by cat_standalone.sh matlabbatch{1}.spm.tools.cat.estwrite.data = ''; % Entry for choosing TPM % Remove comments and edit entry if you would like to change the parameter. % Otherwise the default value from cat_defaults.m is used. % Or use 1st parameter field, that will be dynamically replaced by cat_standalone.sh %matlabbatch{1}.spm.tools.cat.estwrite.opts.tpm = ''; % Entry for choosing shooting template % Remove comments and edit entry if you would like to change the parameter. % Otherwise the default value from cat_defaults.m is used. % Or use 2nd parameter field, that will be dynamically replaced by cat_standalone.sh %matlabbatch{1}.spm.tools.cat.estwrite.extopts.registration.regmethod.shooting.shootingtpm = ''; % Strength of Shooting registration: 0 - Dartel, eps (fast), 0.5 (default) to 1 (accurate) optimized Shooting, 4 - default Shooting; default 0.5 %matlabbatch{1}.spm.tools.cat.estwrite.extopts.registration.regmethod.shooting.regstr = 0.5; % voxel size for normalized data (EXPERIMENTAL: inf - use Tempate values) matlabbatch{1}.spm.tools.cat.estwrite.extopts.registration.vox = 1.5; % additional bounding box matlabbatch{1}.spm.tools.cat.estwrite.extopts.registration.bb = 12; % Affine regularisation (SPM12 default = mni) - '';'mni';'eastern';'subj';'none';'rigid' matlabbatch{1}.spm.tools.cat.estwrite.opts.affreg = 'mni'; % Strength of the bias correction that controls the biasreg and biasfwhm parameter (CAT only!) % 0 - use SPM parameter; eps - ultralight, 0.25 - light, 0.5 - medium, 0.75 - strong, and 1 - heavy corrections % job.opts.biasreg = min( 10 , max( 0 , 10^-(job.opts.biasstr*2 + 2) )); % job.opts.biasfwhm = min( inf , max( 30 , 30 + 60*job.opts.biasstr )); matlabbatch{1}.spm.tools.cat.estwrite.opts.biasstr = 0.5; % Affine PreProcessing (APP) with rough bias correction and brain extraction for special anatomies (nonhuman/neonates) % 0 - none; 1070 - default; [1 - light; 2 - full; 1144 - update of 1070, 5 - animal (no affreg)] matlabbatch{1}.spm.tools.cat.estwrite.extopts.segmentation.APP = 1070; % Strength of the local adaptation: 0 to 1; default 0.5 matlabbatch{1}.spm.tools.cat.estwrite.extopts.segmentation.LASstr = 0.5; % Strength of the noise correction: 0 to 1; 0 - no filter, -Inf - auto, 1 - full, 2 - ISARNLM (else SANLM), default -Inf matlabbatch{1}.spm.tools.cat.estwrite.extopts.segmentation.NCstr = -Inf; % Strength of skull-stripping: 0 - SPM approach; eps to 1 - gcut; 2 - new APRG approach; -1 - no skull-stripping (already skull-stripped); default = 2 matlabbatch{1}.spm.tools.cat.estwrite.extopts.segmentation.gcutstr = 2; % Strength of the cleanup process: 0 to 1; default 0.5 matlabbatch{1}.spm.tools.cat.estwrite.extopts.segmentation.cleanupstr = 0.5; % resolution handling: 'native','fixed','best', 'optimal' matlabbatch{1}.spm.tools.cat.estwrite.extopts.segmentation.restypes.optimal = [1 0.3]; % use center-of-mass approach for estimating origin matlabbatch{1}.spm.tools.cat.estwrite.extopts.segmentation.setCOM = 1; % modify affine scaling matlabbatch{1}.spm.tools.cat.estwrite.extopts.segmentation.affmod = 0; % use k-means AMAP approach or SPM segmentation for initial segmentation matlabbatch{1}.spm.tools.cat.estwrite.extopts.segmentation.spm_kamap = 0; % Correction of WM hyperintensities: 0 - no correction, 1 - only for Dartel/Shooting % 2 - also correct segmentation (to WM), 3 - handle as separate class; default 1 matlabbatch{1}.spm.tools.cat.estwrite.extopts.segmentation.WMHC = 2; % Stroke lesion correction (SLC): 0 - no correction, 1 - handling of manual lesion that have to be set to zero! % 2 - automatic lesion detection (in development) matlabbatch{1}.spm.tools.cat.estwrite.extopts.segmentation.SLC = 0; % surface options matlabbatch{1}.spm.tools.cat.estwrite.extopts.surface.pbtres = 0.5; matlabbatch{1}.spm.tools.cat.estwrite.extopts.surface.pbtmethod = 'pbt2x'; matlabbatch{1}.spm.tools.cat.estwrite.extopts.surface.reduce_mesh = 1; matlabbatch{1}.spm.tools.cat.estwrite.extopts.surface.scale_cortex = 0.7; matlabbatch{1}.spm.tools.cat.estwrite.extopts.surface.add_parahipp = 0.1; matlabbatch{1}.spm.tools.cat.estwrite.extopts.surface.close_parahipp = 0; matlabbatch{1}.spm.tools.cat.estwrite.extopts.surface.SRP = 22; % set this to 1 for skipping preprocessing if already processed data exist matlabbatch{1}.spm.tools.cat.estwrite.extopts.admin.lazy = 0; % catch errors: 0 - stop with error (default); 1 - catch preprocessing errors (requires MATLAB 2008 or higher); matlabbatch{1}.spm.tools.cat.estwrite.extopts.admin.ignoreErrors = 1; % verbose output: 1 - default; 2 - details; 3 - write debugging files matlabbatch{1}.spm.tools.cat.estwrite.extopts.admin.verb = 2; % display and print out pdf-file of results: 0 - off, 2 - volume only, 2 - volume and surface (default) matlabbatch{1}.spm.tools.cat.estwrite.extopts.admin.print = 2; % surface and thickness creation: 0 - no (default), 1 - lh+rh, 2 - lh+rh+cerebellum, % 3 - lh, 4 - rh, 5 - lh+rh (fast, no registration, only for quick quality check and not for analysis), % 6 - lh+rh+cerebellum (fast, no registration, only for quick quality check and not for analysis) % 9 - thickness only (for ROI analysis, experimental!) % +10 to estimate WM and CSF width/depth/thickness (experimental!) matlabbatch{1}.spm.tools.cat.estwrite.output.surface = 1; % BIDS output matlabbatch{1}.spm.tools.cat.estwrite.output.BIDS.BIDSno = 1; % define here volume atlases matlabbatch{1}.spm.tools.cat.estwrite.output.ROImenu.atlases.neuromorphometrics = 1; matlabbatch{1}.spm.tools.cat.estwrite.output.ROImenu.atlases.lpba40 = 1; matlabbatch{1}.spm.tools.cat.estwrite.output.ROImenu.atlases.cobra = 1; matlabbatch{1}.spm.tools.cat.estwrite.output.ROImenu.atlases.hammers = 1; matlabbatch{1}.spm.tools.cat.estwrite.output.ROImenu.atlases.ibsr = 1; matlabbatch{1}.spm.tools.cat.estwrite.output.ROImenu.atlases.aal3 = 1; matlabbatch{1}.spm.tools.cat.estwrite.output.ROImenu.atlases.mori = 1; matlabbatch{1}.spm.tools.cat.estwrite.output.ROImenu.atlases.thalamus = 1; matlabbatch{1}.spm.tools.cat.estwrite.output.ROImenu.atlases.anatomy3 = 1; matlabbatch{1}.spm.tools.cat.estwrite.output.ROImenu.atlases.julichbrain = 1; matlabbatch{1}.spm.tools.cat.estwrite.output.ROImenu.atlases.Schaefer2018_100Parcels_17Networks_order = 1; matlabbatch{1}.spm.tools.cat.estwrite.output.ROImenu.atlases.Schaefer2018_200Parcels_17Networks_order = 1; matlabbatch{1}.spm.tools.cat.estwrite.output.ROImenu.atlases.Schaefer2018_400Parcels_17Networks_order = 1; matlabbatch{1}.spm.tools.cat.estwrite.output.ROImenu.atlases.Schaefer2018_600Parcels_17Networks_order = 1; matlabbatch{1}.spm.tools.cat.estwrite.output.ROImenu.atlases.ownatlas = {'code/Tian_Subcortex_S2_7T.nii'}; % Writing options (see cat_defaults for the description of parameters) % native 0/1 (none/yes) % warped 0/1 (none/yes) % mod 0/1/2/3 (none/affine+nonlinear/nonlinear only/both) % dartel 0/1/2/3 (none/rigid/affine/both) % GM/WM/CSF/WMH matlabbatch{1}.spm.tools.cat.estwrite.output.GM.native = 1; matlabbatch{1}.spm.tools.cat.estwrite.output.GM.warped = 1; matlabbatch{1}.spm.tools.cat.estwrite.output.GM.mod = 3; matlabbatch{1}.spm.tools.cat.estwrite.output.GM.dartel = 2; matlabbatch{1}.spm.tools.cat.estwrite.output.WM.native = 1; matlabbatch{1}.spm.tools.cat.estwrite.output.WM.warped = 1; matlabbatch{1}.spm.tools.cat.estwrite.output.WM.mod = 3; matlabbatch{1}.spm.tools.cat.estwrite.output.WM.dartel = 2; matlabbatch{1}.spm.tools.cat.estwrite.output.CSF.native = 1; matlabbatch{1}.spm.tools.cat.estwrite.output.CSF.warped = 1; matlabbatch{1}.spm.tools.cat.estwrite.output.CSF.mod = 3; matlabbatch{1}.spm.tools.cat.estwrite.output.CSF.dartel = 2; matlabbatch{1}.spm.tools.cat.estwrite.output.WMH.native = 0; matlabbatch{1}.spm.tools.cat.estwrite.output.WMH.warped = 0; matlabbatch{1}.spm.tools.cat.estwrite.output.WMH.mod = 0; matlabbatch{1}.spm.tools.cat.estwrite.output.WMH.dartel = 0; % stroke lesion tissue maps (only for opt.extopts.SLC>0) - in development matlabbatch{1}.spm.tools.cat.estwrite.output.SL.native = 0; matlabbatch{1}.spm.tools.cat.estwrite.output.SL.warped = 0; matlabbatch{1}.spm.tools.cat.estwrite.output.SL.mod = 0; matlabbatch{1}.spm.tools.cat.estwrite.output.SL.dartel = 0; % Tissue classes 4-6 to create own TPMs matlabbatch{1}.spm.tools.cat.estwrite.output.TPMC.native = 0; matlabbatch{1}.spm.tools.cat.estwrite.output.TPMC.warped = 0; matlabbatch{1}.spm.tools.cat.estwrite.output.TPMC.mod = 0; matlabbatch{1}.spm.tools.cat.estwrite.output.TPMC.dartel = 0; % atlas maps (for evaluation) matlabbatch{1}.spm.tools.cat.estwrite.output.atlas.native = 1; % label % background=0, CSF=1, GM=2, WM=3, WMH=4 (if opt.extopts.WMHC==3), SL=1.5 (if opt.extopts.SLC>0) matlabbatch{1}.spm.tools.cat.estwrite.output.label.native = 1; matlabbatch{1}.spm.tools.cat.estwrite.output.label.warped = 1; matlabbatch{1}.spm.tools.cat.estwrite.output.label.dartel = 0; % bias and noise corrected, global intensity normalized matlabbatch{1}.spm.tools.cat.estwrite.output.bias.native = 1; matlabbatch{1}.spm.tools.cat.estwrite.output.bias.warped = 1; matlabbatch{1}.spm.tools.cat.estwrite.output.bias.dartel = 0; % bias and noise corrected, (locally - if LAS>0) intensity normalized matlabbatch{1}.spm.tools.cat.estwrite.output.las.native = 1; matlabbatch{1}.spm.tools.cat.estwrite.output.las.warped = 1; matlabbatch{1}.spm.tools.cat.estwrite.output.las.dartel = 0; % jacobian determinant 0/1 (none/yes) matlabbatch{1}.spm.tools.cat.estwrite.output.jacobianwarped = 1; % deformations, order is [forward inverse] matlabbatch{1}.spm.tools.cat.estwrite.output.warps = [1 1]; % registration matrix matlabbatch{1}.spm.tools.cat.estwrite.output.rmat = 1; % deface native intensity normalized images in native space matlabbatch{2}.spm.util.deface.images(1) = cfg_dep('CAT12: Segmentation: Native Bias Corr. Image', substruct('.','val', '{}',{1}, '.','val', '{}',{1}, '.','val', '{}',{1}, '.','val', '{}',{1}), substruct('()',{1}, '.','biascorr', '()',{':'})); matlabbatch{2}.spm.util.deface.images(2) = cfg_dep('CAT12: Segmentation: Native LAS Bias Corr. Image', substruct('.','val', '{}',{1}, '.','val', '{}',{1}, '.','val', '{}',{1}, '.','val', '{}',{1}), substruct('()',{1}, '.','ibiascorr', '()',{':'})); % and remove images that are not defaced matlabbatch{3}.cfg_basicio.file_dir.file_ops.file_move.files(1) = cfg_dep('CAT12: Segmentation: Native Bias Corr. Image', substruct('.','val', '{}',{1}, '.','val', '{}',{1}, '.','val', '{}',{1}, '.','val', '{}',{1}), substruct('()',{1}, '.','biascorr', '()',{':'})); matlabbatch{3}.cfg_basicio.file_dir.file_ops.file_move.files(2) = cfg_dep('CAT12: Segmentation: Native LAS Bias Corr. Image', substruct('.','val', '{}',{1}, '.','val', '{}',{1}, '.','val', '{}',{1}, '.','val', '{}',{1}), substruct('()',{1}, '.','ibiascorr', '()',{':'})); matlabbatch{3}.cfg_basicio.file_dir.file_ops.file_move.action.delete = false;