def createAndRun(sessions, environment, experiment, pipeline, cluster):
from baw_exp import OpenSubjectDatabase
from utilities.misc import add_dict
master_config = {}
for configDict in [environment, experiment, pipeline, cluster]:
master_config = add_dict(master_config, configDict)
database = OpenSubjectDatabase(experiment['cachedir'], ['all'], environment['prefix'], experiment['dbfile'])
database.open_connection()
try:
all_sessions = database.getAllSessions()
if not set(sessions) <= set(all_sessions) and 'all' not in sessions:
missing = set(sessions) - set(all_sessions)
assert len(missing) == 0, "Requested sessions are missing from the database: {0}".format(missing)
elif 'all' in sessions:
sessions = set(all_sessions)
else:
sessions = set(sessions)
for session in sessions:
_dict = {}
_dict['session'] = session
_dict['project'] = database.getProjFromSession(session)
_dict['subject'] = database.getSubjFromSession(session)
_dict['T1s'] = database.getFilenamesByScantype(session, ['T1-15', 'T1-30'])
_dict['T2s'] = database.getFilenamesByScantype(session, ['T2-15', 'T2-30'])
_dict['PDs'] = database.getFilenamesByScantype(session, ['PD-15', 'PD-30'])
_dict['FLs'] = database.getFilenamesByScantype(session, ['FL-15', 'FL-30'])
_dict['OTs'] = database.getFilenamesByScantype(session, ['OTHER-15', 'OTHER-30'])
workflow = create_singleSession(_dict, master_config, 'Linear', 'singleSession_{0}_{1}'.format(_dict['subject'], _dict['session']))
workflow.run(plugin='SGEGraph', plugin_args=master_config['plugin_args'])
except:
raise
finally:
try:
database.close_connection()
except:
pass
def createAndRun(sessions, environment, experiment, pipeline, cluster, useSentinal, dryRun):
from baw_exp import OpenSubjectDatabase
from utilities.misc import add_dict
from workflows.utils import run_workflow
master_config = {}
for configDict in [environment, experiment, pipeline, cluster]:
master_config = add_dict(master_config, configDict)
database = OpenSubjectDatabase(experiment['cachedir'], ['all'], environment['prefix'], experiment['dbfile'])
database.open_connection()
try:
all_sessions = database.getAllSessions()
if not set(sessions) <= set(all_sessions) and 'all' not in sessions:
missing = set(sessions) - set(all_sessions)
assert len(missing) == 0, "Requested sessions are missing from the database: {0}\n\n{1}".format(missing,all_sessions)
elif 'all' in sessions:
sessions = set(all_sessions)
else:
sessions = set(sessions)
print("!=" * 40)
print("Doing sessions {0}".format(sessions))
print("!=" * 40)
for session in sessions:
_dict = {}
subject = database.getSubjFromSession(session)
_dict['session'] = session
_dict['project'] = database.getProjFromSession(session)
_dict['subject'] = subject
_dict['T1s'] = database.getFilenamesByScantype(session, ['T1-15', 'T1-30'])
_dict['T2s'] = database.getFilenamesByScantype(session, ['T2-15', 'T2-30'])
_dict['PDs'] = database.getFilenamesByScantype(session, ['PD-15', 'PD-30'])
_dict['FLs'] = database.getFilenamesByScantype(session, ['FL-15', 'FL-30'])
_dict['OTs'] = database.getFilenamesByScantype(session, ['OTHER-15', 'OTHER-30'])
sentinal_file_basedir = os.path.join(
master_config['resultdir'],
_dict['project'],
_dict['subject'],
_dict['session']
)
sentinal_file_list = list()
sentinal_file_list.append(os.path.join(sentinal_file_basedir))
if 'denoise' in master_config['components']:
# # NO SENTINAL FILE
pass
# # Use t1 average sentinal file if specified.
if 'landmark' in master_config['components']:
sentinal_file_list.append(os.path.join(
sentinal_file_basedir,
"ACPCAlign",
"landmarkInitializer_atlas_to_subject_transform.h5"
))
if 'tissue_classify' in master_config['components']:
for tc_file in ["complete_brainlabels_seg.nii.gz", "t1_average_BRAINSABC.nii.gz"]:
sentinal_file_list.append(os.path.join(
sentinal_file_basedir,
"TissueClassify",
tc_file
))
if 'warp_atlas_to_subject' in master_config['components']:
sentinal_file_list.append(os.path.join(
sentinal_file_basedir,
"WarpedAtlas2Subject",
"rho.nii.gz"
))
sentinal_file_list.append(os.path.join(
sentinal_file_basedir,
"WarpedAtlas2Subject",
"left_hemisphere_wm.nii.gz"
))
if 'malf_2012_neuro' in master_config['components']:
sentinal_file_list.append(os.path.join(
sentinal_file_basedir,
"TissueClassify",
"MALF_HDAtlas20_2015_fs_standard_label.nii.gz"
))
sentinal_file_list.append(os.path.join(
sentinal_file_basedir,
"TissueClassify",
"MALF_HDAtlas20_2015_label.nii.gz"
))
sentinal_file_list.append(os.path.join(
sentinal_file_basedir,
"TissueClassify",
"MALF_HDAtlas20_2015_lobar_label.nii.gz"
))
sentinal_file_list.append(os.path.join(
sentinal_file_basedir,
"TissueClassify",
"MALF_HDAtlas20_2015_CSFVBInjected_label.nii.gz"
))
if 'malf_2015_wholebrain' in master_config['components']:
pass
if master_config['workflow_phase'] == 'atlas-based-reference':
atlasDirectory = os.path.join(master_config['atlascache'], 'spatialImages', 'rho.nii.gz')
else:
atlasDirectory = os.path.join(master_config['previousresult'], subject, 'Atlas', 'AVG_rho.nii.gz')
if os.path.exists(atlasDirectory):
print("LOOKING FOR DIRECTORY {0}".format(atlasDirectory))
else:
print("MISSING REQUIRED ATLAS INPUT {0}".format(atlasDirectory))
print("SKIPPING: {0} prerequisites missing".format(session))
continue
## Use different sentinal file if segmentation specified.
from workflows.baseline import DetermineIfSegmentationShouldBeDone
do_BRAINSCut_Segmentation = DetermineIfSegmentationShouldBeDone(master_config)
if do_BRAINSCut_Segmentation:
sentinal_file_list.append(os.path.join(
sentinal_file_basedir,
"CleanedDenoisedRFSegmentations",
"allLabels_seg.nii.gz"
))
def allPathsExists(list_of_paths):
is_missing = False
for ff in list_of_paths:
if not os.path.exists(ff):
is_missing = True
print("MISSING: {0}".format(ff))
return not is_missing
if useSentinal and allPathsExists(sentinal_file_list):
print("SKIPPING: {0} exists".format(sentinal_file_list))
else:
print("PROCESSING INCOMPLETE: at least 1 required file does not exists")
if dryRun == False:
workflow = _create_singleSession(_dict, master_config, 'Linear',
'singleSession_{0}_{1}'.format(_dict['subject'], _dict['session']))
print("Starting session {0}".format(session))
# HACK Hard-coded to SGEGraph, but --wfrun is ignored completely
run_workflow(workflow, plugin=master_config['plugin_name'],
plugin_args=master_config['plugin_args'])
else:
print("EXITING WITHOUT WORK DUE TO dryRun flag")
except:
raise
finally:
try:
database.close_connection()
except:
pass
def createAndRun(sessions, environment, experiment, pipeline, cluster,
useSentinal, dryRun):
from baw_exp import OpenSubjectDatabase
from utilities.misc import add_dict
from workflows.utils import run_workflow
master_config = {}
for configDict in [environment, experiment, pipeline, cluster]:
master_config = add_dict(master_config, configDict)
database = OpenSubjectDatabase(experiment['cachedir'], ['all'],
environment['prefix'], experiment['dbfile'])
database.open_connection()
try:
all_sessions = database.getAllSessions()
if not set(sessions) <= set(all_sessions) and 'all' not in sessions:
missing = set(sessions) - set(all_sessions)
assert len(
missing
) == 0, "Requested sessions are missing from the database: {0}\n\n{1}".format(
missing, all_sessions)
elif 'all' in sessions:
sessions = set(all_sessions)
else:
sessions = set(sessions)
print("!=" * 40)
print("Doing sessions {0}".format(sessions))
print("!=" * 40)
for session in sessions:
_dict = {}
subject = database.getSubjFromSession(session)
_dict['session'] = session
_dict['project'] = database.getProjFromSession(session)
_dict['subject'] = subject
_dict['T1s'] = database.getFilenamesByScantype(
session, ['T1-15', 'T1-30'])
_dict['T2s'] = database.getFilenamesByScantype(
session, ['T2-15', 'T2-30'])
_dict['PDs'] = database.getFilenamesByScantype(
session, ['PD-15', 'PD-30'])
_dict['FLs'] = database.getFilenamesByScantype(
session, ['FL-15', 'FL-30'])
_dict['OTs'] = database.getFilenamesByScantype(
session, ['OTHER-15', 'OTHER-30'])
sentinal_file_basedir = os.path.join(master_config['resultdir'],
_dict['project'],
_dict['subject'],
_dict['session'])
sentinal_file_list = list()
sentinal_file_list.append(os.path.join(sentinal_file_basedir))
if 'denoise' in master_config['components']:
# # NO SENTINAL FILE
pass
# # Use t1 average sentinal file if specified.
if 'landmark' in master_config['components']:
sentinal_file_list.append(
os.path.join(
sentinal_file_basedir, "ACPCAlign",
"landmarkInitializer_atlas_to_subject_transform.h5"))
if 'tissue_classify' in master_config['components']:
for tc_file in [
"complete_brainlabels_seg.nii.gz",
"t1_average_BRAINSABC.nii.gz"
]:
sentinal_file_list.append(
os.path.join(sentinal_file_basedir, "TissueClassify",
tc_file))
if 'warp_atlas_to_subject' in master_config['components']:
sentinal_file_list.append(
os.path.join(sentinal_file_basedir, "WarpedAtlas2Subject",
"rho.nii.gz"))
sentinal_file_list.append(
os.path.join(sentinal_file_basedir, "WarpedAtlas2Subject",
"left_hemisphere_wm.nii.gz"))
if 'malf_2012_neuro' in master_config['components']:
sentinal_file_list.append(
os.path.join(
sentinal_file_basedir, "TissueClassify",
"MALF_HDAtlas20_2015_fs_standard_label.nii.gz"))
sentinal_file_list.append(
os.path.join(sentinal_file_basedir, "TissueClassify",
"MALF_HDAtlas20_2015_label.nii.gz"))
sentinal_file_list.append(
os.path.join(sentinal_file_basedir, "TissueClassify",
"MALF_HDAtlas20_2015_lobar_label.nii.gz"))
sentinal_file_list.append(
os.path.join(
sentinal_file_basedir, "TissueClassify",
"MALF_HDAtlas20_2015_CSFVBInjected_label.nii.gz"))
if 'malf_2015_wholebrain' in master_config['components']:
pass
if master_config['workflow_phase'] == 'atlas-based-reference':
atlasDirectory = os.path.join(master_config['atlascache'],
'spatialImages', 'rho.nii.gz')
else:
atlasDirectory = os.path.join(master_config['previousresult'],
subject, 'Atlas',
'AVG_rho.nii.gz')
if os.path.exists(atlasDirectory):
print("LOOKING FOR DIRECTORY {0}".format(atlasDirectory))
else:
print(
"MISSING REQUIRED ATLAS INPUT {0}".format(atlasDirectory))
print("SKIPPING: {0} prerequisites missing".format(session))
continue
## Use different sentinal file if segmentation specified.
from workflows.baseline import DetermineIfSegmentationShouldBeDone
do_BRAINSCut_Segmentation = DetermineIfSegmentationShouldBeDone(
master_config)
if do_BRAINSCut_Segmentation:
sentinal_file_list.append(
os.path.join(sentinal_file_basedir,
"CleanedDenoisedRFSegmentations",
"allLabels_seg.nii.gz"))
def allPathsExists(list_of_paths):
is_missing = False
for ff in list_of_paths:
if not os.path.exists(ff):
is_missing = True
print("MISSING: {0}".format(ff))
return not is_missing
if useSentinal and allPathsExists(sentinal_file_list):
print("SKIPPING: {0} exists".format(sentinal_file_list))
else:
print(
"PROCESSING INCOMPLETE: at least 1 required file does not exists"
)
if dryRun == False:
workflow = _create_singleSession(
_dict, master_config, 'Linear',
'singleSession_{0}_{1}'.format(_dict['subject'],
_dict['session']))
print("Starting session {0}".format(session))
# HACK Hard-coded to SGEGraph, but --wfrun is ignored completely
run_workflow(workflow,
plugin=master_config['plugin_name'],
plugin_args=master_config['plugin_args'])
else:
print("EXITING WITHOUT WORK DUE TO dryRun flag")
except:
raise
finally:
try:
database.close_connection()
except:
pass
def createAndRun(sessions, environment, experiment, pipeline, cluster,
useSentinal, dryRun):
from baw_exp import OpenSubjectDatabase
from utilities.misc import add_dict
from collections import OrderedDict
import sys
from workflows.utils import run_workflow
master_config = {}
for configDict in [environment, experiment, pipeline, cluster]:
master_config = add_dict(master_config, configDict)
database = OpenSubjectDatabase(experiment['cachedir'], ['all'],
environment['prefix'], experiment['dbfile'])
database.open_connection()
try:
all_sessions = database.getAllSessions()
if not set(sessions) <= set(all_sessions) and 'all' not in sessions:
missing = set(sessions) - set(all_sessions)
assert len(
missing
) == 0, "Requested sessions are missing from the database: {0}\n\n{1}".format(
missing, all_sessions)
elif 'all' in sessions:
sessions = set(all_sessions)
else:
sessions = set(sessions)
print("!=" * 40)
print("Doing sessions {0}".format(sessions))
print("!=" * 40)
for session in sessions:
_dict = OrderedDict()
t1_list = database.getFilenamesByScantype(session,
['T1-15', 'T1-30'])
if len(t1_list) == 0:
print(
"ERROR: Skipping session {0} for subject {1} due to missing T1's"
.format(session, subject))
print("REMOVE OR FIX BEFORE CONTINUING")
continue
subject = database.getSubjFromSession(session)
_dict['session'] = session
_dict['project'] = database.getProjFromSession(session)
_dict['subject'] = subject
_dict['T1s'] = t1_list
_dict['T2s'] = database.getFilenamesByScantype(
session, ['T2-15', 'T2-30'])
_dict['BadT2'] = False
if _dict['T2s'] == database.getFilenamesByScantype(
session, ['T2-15']):
print("This T2 is not going to be used for JointFusion")
print("This T2 is not going to be used for JointFusion")
print("This T2 is not going to be used for JointFusion")
print("This T2 is not going to be used for JointFusion")
print(_dict['T2s'])
_dict['BadT2'] = True
_dict['PDs'] = database.getFilenamesByScantype(
session, ['PD-15', 'PD-30'])
_dict['FLs'] = database.getFilenamesByScantype(
session, ['FL-15', 'FL-30'])
_dict['EMSP'] = database.getFilenamesByScantype(session, ['EMSP'])
_dict['OTHERs'] = database.getFilenamesByScantype(
session, ['OTHER-15', 'OTHER-30'])
sentinal_file_basedir = os.path.join(master_config['resultdir'],
_dict['project'],
_dict['subject'],
_dict['session'])
sentinal_file_list = list()
sentinal_file_list.append(os.path.join(sentinal_file_basedir))
if 'denoise' in master_config['components']:
# # NO SENTINAL FILE
pass
# # Use t1 average sentinal file if specified.
if 'landmark' in master_config['components']:
sentinal_file_list.append(
os.path.join(
sentinal_file_basedir, "ACPCAlign",
"landmarkInitializer_atlas_to_subject_transform.h5"))
if 'tissue_classify' in master_config['components']:
for tc_file in [
"complete_brainlabels_seg.nii.gz",
"t1_average_BRAINSABC.nii.gz"
]:
sentinal_file_list.append(
os.path.join(sentinal_file_basedir, "TissueClassify",
tc_file))
if 'warp_atlas_to_subject' in master_config['components']:
warp_atlas_file_list = [
"hncma_atlas.nii.gz", "l_accumben_ProbabilityMap.nii.gz",
"l_caudate_ProbabilityMap.nii.gz",
"l_globus_ProbabilityMap.nii.gz",
"l_hippocampus_ProbabilityMap.nii.gz",
"l_putamen_ProbabilityMap.nii.gz",
"l_thalamus_ProbabilityMap.nii.gz",
"left_hemisphere_wm.nii.gz", "phi.nii.gz",
"r_accumben_ProbabilityMap.nii.gz",
"r_caudate_ProbabilityMap.nii.gz",
"r_globus_ProbabilityMap.nii.gz",
"r_hippocampus_ProbabilityMap.nii.gz",
"r_putamen_ProbabilityMap.nii.gz",
"r_thalamus_ProbabilityMap.nii.gz", "rho.nii.gz",
"right_hemisphere_wm.nii.gz",
"template_WMPM2_labels.nii.gz",
"template_headregion.nii.gz",
"template_leftHemisphere.nii.gz",
"template_nac_labels.nii.gz",
"template_rightHemisphere.nii.gz",
"template_ventricles.nii.gz", "theta.nii.gz"
]
for ff in warp_atlas_file_list:
sentinal_file_list.append(
os.path.join(sentinal_file_basedir,
"WarpedAtlas2Subject", ff))
if 'jointfusion_2015_wholebrain' in master_config['components']:
sentinal_file_list.append(
os.path.join(
sentinal_file_basedir, "TissueClassify",
"JointFusion_HDAtlas20_2015_lobar_label.nii.gz"))
sentinal_file_list.append(
os.path.join(sentinal_file_basedir, "TissueClassify",
"lobeVolumes_JSON.json"))
if master_config['workflow_phase'] == 'atlas-based-reference':
atlasDirectory = os.path.join(master_config['atlascache'],
'spatialImages', 'rho.nii.gz')
sentinal_file_list.append(atlasDirectory)
else:
atlasDirectory = os.path.join(master_config['previousresult'],
subject, 'Atlas',
'AVG_rho.nii.gz')
sentinal_file_list.append(atlasDirectory)
sentinal_file_list.append(
os.path.join(master_config['previousresult'], subject,
'Atlas', 'AVG_template_headregion.nii.gz'))
if os.path.exists(atlasDirectory):
print("LOOKING FOR DIRECTORY {0}".format(atlasDirectory))
else:
print(
"MISSING REQUIRED ATLAS INPUT {0}".format(atlasDirectory))
print("SKIPPING: {0} prerequisites missing".format(session))
continue
## Use different sentinal file if segmentation specified.
from workflows.baseline import DetermineIfSegmentationShouldBeDone
do_BRAINSCut_Segmentation = DetermineIfSegmentationShouldBeDone(
master_config)
if do_BRAINSCut_Segmentation:
sentinal_file_list.append(
os.path.join(sentinal_file_basedir,
"CleanedDenoisedRFSegmentations",
"allLabels_seg.nii.gz"))
def allPathsExists(list_of_paths):
is_missing = False
for ff in list_of_paths:
if not os.path.exists(ff):
is_missing = True
print("MISSING: {0}".format(ff))
return not is_missing
if useSentinal and allPathsExists(sentinal_file_list):
print("SKIPPING: {0} exists".format(sentinal_file_list))
else:
print(
"PROCESSING INCOMPLETE: at least 1 required file does not exists"
)
if dryRun == False:
workflow = _create_singleSession(
_dict, master_config, 'Linear',
'singleSession_{0}_{1}'.format(_dict['subject'],
_dict['session']))
print("Starting session {0}".format(session))
# HACK Hard-coded to SGEGraph, but --wfrun is ignored completely
run_workflow(workflow,
plugin=master_config['plugin_name'],
plugin_args=master_config['plugin_args'])
else:
print("EXITING WITHOUT WORK DUE TO dryRun flag")
except:
raise
finally:
try:
database.close_connection()
except:
pass
def createAndRun(sessions, environment, experiment, pipeline, cluster, useSentinal, dryRun):
from baw_exp import OpenSubjectDatabase
from utilities.misc import add_dict
from collections import OrderedDict
import sys
from workflows.utils import run_workflow
master_config = {}
for configDict in [environment, experiment, pipeline, cluster]:
master_config = add_dict(master_config, configDict)
database = OpenSubjectDatabase(experiment['cachedir'], ['all'], environment['prefix'], experiment['dbfile'])
database.open_connection()
try:
all_sessions = database.getAllSessions()
if not set(sessions) <= set(all_sessions) and 'all' not in sessions:
missing = set(sessions) - set(all_sessions)
assert len(missing) == 0, "Requested sessions are missing from the database: {0}\n\n{1}".format(missing,
all_sessions)
elif 'all' in sessions:
sessions = set(all_sessions)
else:
sessions = set(sessions)
print("!=" * 40)
print("Doing sessions {0}".format(sessions))
print("!=" * 40)
for session in sessions:
_dict = OrderedDict()
t1_list = database.getFilenamesByScantype(session, ['T1-15', 'T1-30'])
if len(t1_list) == 0:
print("ERROR: Skipping session {0} for subject {1} due to missing T1's".format(session, subject))
print("REMOVE OR FIX BEFORE CONTINUING")
continue
subject = database.getSubjFromSession(session)
_dict['session'] = session
_dict['project'] = database.getProjFromSession(session)
_dict['subject'] = subject
_dict['T1s'] = t1_list
_dict['T2s'] = database.getFilenamesByScantype(session, ['T2-15', 'T2-30'])
_dict['BadT2'] = False
if _dict['T2s'] == database.getFilenamesByScantype(session, ['T2-15']):
print("This T2 is not going to be used for JointFusion")
print("This T2 is not going to be used for JointFusion")
print("This T2 is not going to be used for JointFusion")
print("This T2 is not going to be used for JointFusion")
print(_dict['T2s'])
_dict['BadT2'] = True
_dict['PDs'] = database.getFilenamesByScantype(session, ['PD-15', 'PD-30'])
_dict['FLs'] = database.getFilenamesByScantype(session, ['FL-15', 'FL-30'])
_dict['EMSP'] = database.getFilenamesByScantype(session, ['EMSP'])
_dict['OTHERs'] = database.getFilenamesByScantype(session, ['OTHER-15', 'OTHER-30'])
sentinal_file_basedir = os.path.join(
master_config['resultdir'],
_dict['project'],
_dict['subject'],
_dict['session']
)
sentinal_file_list = list()
sentinal_file_list.append(os.path.join(sentinal_file_basedir))
if 'denoise' in master_config['components']:
# # NO SENTINAL FILE
pass
# # Use t1 average sentinal file if specified.
if 'landmark' in master_config['components']:
sentinal_file_list.append(os.path.join(
sentinal_file_basedir,
"ACPCAlign",
"landmarkInitializer_atlas_to_subject_transform.h5"
))
if 'tissue_classify' in master_config['components']:
for tc_file in ["complete_brainlabels_seg.nii.gz", "t1_average_BRAINSABC.nii.gz"]:
sentinal_file_list.append(os.path.join(
sentinal_file_basedir,
"TissueClassify",
tc_file
))
if 'warp_atlas_to_subject' in master_config['components']:
warp_atlas_file_list = [
"hncma_atlas.nii.gz",
"l_accumben_ProbabilityMap.nii.gz",
"l_caudate_ProbabilityMap.nii.gz",
"l_globus_ProbabilityMap.nii.gz",
"l_hippocampus_ProbabilityMap.nii.gz",
"l_putamen_ProbabilityMap.nii.gz",
"l_thalamus_ProbabilityMap.nii.gz",
"left_hemisphere_wm.nii.gz",
"phi.nii.gz",
"r_accumben_ProbabilityMap.nii.gz",
"r_caudate_ProbabilityMap.nii.gz",
"r_globus_ProbabilityMap.nii.gz",
"r_hippocampus_ProbabilityMap.nii.gz",
"r_putamen_ProbabilityMap.nii.gz",
"r_thalamus_ProbabilityMap.nii.gz",
"rho.nii.gz",
"right_hemisphere_wm.nii.gz",
"template_WMPM2_labels.nii.gz",
"template_headregion.nii.gz",
"template_leftHemisphere.nii.gz",
"template_nac_labels.nii.gz",
"template_rightHemisphere.nii.gz",
"template_ventricles.nii.gz",
"theta.nii.gz"
]
for ff in warp_atlas_file_list:
sentinal_file_list.append(os.path.join(
sentinal_file_basedir,
"WarpedAtlas2Subject",
ff
))
if 'jointfusion_2015_wholebrain' in master_config['components']:
sentinal_file_list.append(os.path.join(
sentinal_file_basedir,
"TissueClassify",
"JointFusion_HDAtlas20_2015_lobar_label.nii.gz"
))
sentinal_file_list.append(os.path.join(
sentinal_file_basedir,
"TissueClassify",
"lobeVolumes_JSON.json"
))
if master_config['workflow_phase'] == 'atlas-based-reference':
atlasDirectory = os.path.join(master_config['atlascache'], 'spatialImages', 'rho.nii.gz')
sentinal_file_list.append(atlasDirectory)
else:
atlasDirectory = os.path.join(master_config['previousresult'], subject, 'Atlas', 'AVG_rho.nii.gz')
sentinal_file_list.append(atlasDirectory)
sentinal_file_list.append(
os.path.join(master_config['previousresult'], subject, 'Atlas', 'AVG_template_headregion.nii.gz'))
if os.path.exists(atlasDirectory):
print("LOOKING FOR DIRECTORY {0}".format(atlasDirectory))
else:
print("MISSING REQUIRED ATLAS INPUT {0}".format(atlasDirectory))
print("SKIPPING: {0} prerequisites missing".format(session))
continue
## Use different sentinal file if segmentation specified.
from workflows.baseline import DetermineIfSegmentationShouldBeDone
do_BRAINSCut_Segmentation = DetermineIfSegmentationShouldBeDone(master_config)
if do_BRAINSCut_Segmentation:
sentinal_file_list.append(os.path.join(
sentinal_file_basedir,
"CleanedDenoisedRFSegmentations",
"allLabels_seg.nii.gz"
))
def allPathsExists(list_of_paths):
is_missing = False
for ff in list_of_paths:
if not os.path.exists(ff):
is_missing = True
print("MISSING: {0}".format(ff))
return not is_missing
if useSentinal and allPathsExists(sentinal_file_list):
print("SKIPPING: {0} exists".format(sentinal_file_list))
else:
print("PROCESSING INCOMPLETE: at least 1 required file does not exists")
if dryRun == False:
workflow = _create_singleSession(_dict, master_config, 'Linear',
'singleSession_{0}_{1}'.format(_dict['subject'], _dict['session']))
print("Starting session {0}".format(session))
# HACK Hard-coded to SGEGraph, but --wfrun is ignored completely
run_workflow(workflow, plugin=master_config['plugin_name'],
plugin_args=master_config['plugin_args'])
else:
print("EXITING WITHOUT WORK DUE TO dryRun flag")
except:
raise
finally:
try:
database.close_connection()
except:
pass