|
30 | 30 |
|
31 | 31 | specifically the 2mm versions of:
|
32 | 32 |
|
33 |
| -- `Joint Fusion Atlas <http://mindboggle.info/data/atlases/jointfusion/OASIS-TRT-20_jointfusion_DKT31_CMA_labels_in_MNI152_2mm.nii.gz>`_ |
| 33 | +- `Joint Fusion Atlas <http://mindboggle.info/data/atlases/jointfusion/OASIS-TRT-20_jointfusion_DKT31_CMA_labels_in_MNI152_2mm_v2.nii.gz>`_ |
34 | 34 | - `MNI template <http://mindboggle.info/data/templates/ants/OASIS-30_Atropos_template_in_MNI152_2mm.nii.gz>`_
|
35 | 35 |
|
36 |
| -The 2mm version was generated with:: |
37 |
| -
|
38 |
| - >>> from nipype.interfaces import freesurfer as fs |
39 |
| - >>> rs = fs.Resample() |
40 |
| - >>> rs.inputs.in_file = 'OASIS-TRT-20_jointfusion_DKT31_CMA_labels_in_MNI152.nii.gz' |
41 |
| - >>> rs.inputs.resampled_file = 'OASIS-TRT-20_jointfusion_DKT31_CMA_labels_in_MNI152_2mm.nii.gz' |
42 |
| - >>> rs.inputs.voxel_size = (2., 2., 2.) |
43 |
| - >>> rs.inputs.args = '-rt nearest -ns 1' |
44 |
| - >>> res = rs.run() |
45 |
| -
|
46 | 36 | """
|
47 | 37 |
|
48 | 38 | import os
|
|
60 | 50 |
|
61 | 51 | mlab.MatlabCommand.set_default_matlab_cmd("matlab -nodisplay")
|
62 | 52 | # If SPM is not in your MATLAB path you should add it here
|
63 |
| -mlab.MatlabCommand.set_default_paths('/cm/shared/openmind/spm/spm12b/spm12b_r5918/') |
| 53 | +# mlab.MatlabCommand.set_default_paths('/software/matlab/spm12') |
64 | 54 |
|
65 | 55 | from nipype.algorithms.rapidart import ArtifactDetect
|
66 | 56 | from nipype.algorithms.misc import TSNR
|
@@ -609,10 +599,10 @@ def create_workflow(files,
|
609 | 599 | """
|
610 | 600 |
|
611 | 601 | art = Node(interface=ArtifactDetect(), name="art")
|
612 |
| - art.inputs.use_differences = [True, False] |
| 602 | + art.inputs.use_differences = [True, True] |
613 | 603 | art.inputs.use_norm = True
|
614 | 604 | art.inputs.norm_threshold = norm_threshold
|
615 |
| - art.inputs.zintensity_threshold = 3 |
| 605 | + art.inputs.zintensity_threshold = 9 |
616 | 606 | art.inputs.mask_type = 'spm_global'
|
617 | 607 | art.inputs.parameter_source = 'SPM'
|
618 | 608 |
|
@@ -915,12 +905,11 @@ def get_names(files, suffix):
|
915 | 905 | fsdir = '/software/temp/brainconnect/fsdata'
|
916 | 906 | files = sorted(glob(os.path.abspath('%s/E?/func/rest.nii' % subj_id)))
|
917 | 907 | anat_file = glob(os.path.abspath('%s/EO/anat/anat.nii' % subj_id))[0]
|
918 |
| - target_file = \ |
919 |
| - os.path.abspath(('OASIS-TRT-20_jointfusion_DKT31_CMA_labels_in_MNI152_' |
920 |
| - '2mm_v2.nii.gz')) |
| 908 | + target_file = os.path.abspath('OASIS-30_Atropos_template_in_MNI152_2mm.nii.gz') |
921 | 909 | wf = create_workflow(files, anat_file, target_file, subj_id,
|
922 | 910 | 2.0, 33, vol_fwhm=6.0, surf_fwhm=15.0,
|
923 | 911 | lowpass_freq=0.1, highpass_freq=0.01,
|
| 912 | + norm_threshold=2., |
924 | 913 | subjects_dir=fsdir,
|
925 | 914 | target_subject=['fsaverage5'],
|
926 | 915 | sink_directory=os.getcwd(),
|
|
0 commit comments