init_dissection_afq_clean_csddet_trk.py

import sys
import utils.copy_template as copier

dataset, subject, root = copier.parseCommandLine(sys.argv)

name = copier.subjectNameFromNumber(subject)
repititions = ['0' + `i` for i in range(1,10)] + ["10"]
rep_str = copier.arrToInterpolateString(repititions)
anatomy = copier.getAnatomy(root, dataset, subject)

lifebid_root = "/N/dc2/projects/lifebid/"

mapping = {}
mapping["stn"] = {
    "mat_in": lifebid_root + "code/ccaiafa/Caiafa_Pestilli_paper2015/Results/ETC_Dec2015/Single_TC/fe_structure_{}_96dirs_b2000_1p5iso_STC_run01_500000_SD_STREAM_lmax8_connNUM{}_TRACTS.mat",
    "mat_out": root + "O3D_{}/derivatives/dissection_afq_csddet_trk/sub-{}/dwi/sub-{}_dwi_var-csddetlife_run-{}_var-afq"
}

mapping["hcp3t"] = {
    "mat_in": lifebid_root + "code/ccaiafa/Caiafa_Pestilli_paper2015/Results/ETC_Dec2015/Single_TC/fe_structure_{}_STC_run01_500000_SD_STREAM_lmax8_connNUM{}_TRACTS.mat",
    "mat_out": root + "O3D_{}/derivatives/dissection_afq_csddet_trk/sub-{}/dwi/sub-{}_dwi_var-csddetlife_run-{}_var-afq"
}

# mapping["hcp7t"] = {
#     "mat_in": lifebid_root + "code/ccaiafa/Caiafa_Pestilli_paper2015/paper_datasets/HCP7T/sub-{}/fe_structures/fe_structure_{}_STC_run01_SD_STREAM_lmax8_connNUM{}.mat",
#     "anatomy": lifebid_root + "code/ccaiafa/Caiafa_Pestilli_paper2015/paper_datasets/HCP7T/sub-{}/dwi/data_b2000.nii.gz",
#     "mat_out": root + "O3D_HCP7T/derivatives/tracking_csddet_trk/sub-{}/dwi/sub-{}_dwi_var-csddetlife_run-{}_tract.trk"
# }

if dataset == "hcp7t":
    print "Cannot run for hcp7t yet, these files don't exist"
    sys.exit()


mat_in = mapping[dataset]["mat_in"].format(name, rep_str)
mat_out = mapping[dataset]["mat_out"].format(dataset.upper(), subject, subject, rep_str)
copier.copy(mat_in, mat_out, anatomy = anatomy, action = "clean_afq2trk", dummy = False)