if not os.path.exists(work_dir):
os.makedirs(work_dir)
# go into working dir
os.chdir(work_dir)
lesion_mask_dil = os.path.join(work_dir,
'lesion_mask_mni_dilated.nii.gz')
print lesion_mask_dil
###### Step 1: dilating lesion mask #############################
dilate = DilateImage()
dilate.inputs.in_file = lesion_mask
dilate.inputs.operation = 'modal'
dilate.inputs.out_file = lesion_mask_dil
#print dilate.cmdline
dilate.run()
##### Step 2: get intra-subject gm mask #########################
rest_list = []
rest_str = []
i = 0
print "intra-subject rest mask ..."
for name in glob.glob(data_dir + subject_id + '*' +
'/preprocessed/func/connectivity/' +
'rest_mask_mni3.nii.gz' ):
print name
rest_list.append(name)
i += 1
if i != 1:
work_dir = os.path.join(out_dir, subject_id)
if not os.path.exists(work_dir):
os.makedirs(work_dir)
# go into working dir
os.chdir(work_dir)
lesion_mask_dil = os.path.join(work_dir, 'lesion_mask_mni_dilated.nii.gz')
print lesion_mask_dil
###### Step 1: dilating lesion mask #############################
dilate = DilateImage()
dilate.inputs.in_file = lesion_mask
dilate.inputs.operation = 'modal'
dilate.inputs.out_file = lesion_mask_dil
#print dilate.cmdline
dilate.run()
##### Step 2: get intra-subject gm mask #########################
rest_list = []
rest_str = []
i = 0
print "intra-subject rest mask ..."
for name in glob.glob(data_dir + subject_id + '*' +
'/preprocessed/func/connectivity/' +
'rest_mask_mni3.nii.gz'):
print name
rest_list.append(name)
i += 1
if i != 1:
