init_connectome_tract.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)
algs = [("tens", "dtidet", "lmax2"), ("detr", "csddet", "lmax8"), ("prob", "csdprob", "lmax8")]

mapping = {}
mapping["stn"] = {
    "input": "/N/dc2/projects/lifebid/HCP/Brent/cogs610/reps_data/stn_{}_{}_{}_rep{}.mat",
    "output": root + "O3D_STN/derivatives/connectome_tract/sub-{}/dwi/sub-{}_dwi_var-{}life_run-{}"
}

mapping["hcp3t"] = {
    "input": "/N/dc2/projects/lifebid/HCP/Brent/cogs610/reps_data/hcp_{}_{}_{}_rep{}.mat",
    "output": root + "O3D_HCP3T/derivatives/connectome_tract/sub-{}/dwi/sub-{}_dwi_var-{}life_run-{}"
}

mapping["hcp7t"] = {
    "input": "/N/dc2/projects/lifebid/HCP/Brent/cogs610/7t_reps_data/7T_{}_{}_{}_rep{}.mat",
    "output": root + "O3D_HCP7T/derivatives/connectome_tract/sub-{}/dwi/sub-{}_dwi_var-{}life_run-{}"
}

for a in algs:
    inalg, outalg, lmax = a
    in_str = mapping[dataset]["input"].format(name, inalg, lmax, rep_str)
    out_str = mapping[dataset]["output"].format(subject, subject, outalg, rep_str)
    copier.copy(in_str, out_str, action = "mat2fiber", dummy = False)