@@ -10,6 +10,7 @@ DEFAULT_HCP_PIPELINES_DIR="${HOME}/pipeline_tools/HCPpipelines"
10
10
DEFAULT_FSL_DIR=" /export/HCP/fsl-6.0.1b0"
11
11
DEFAULT_FREESURFER_DIR=" /export/freesurfer-6.0"
12
12
DEFAULT_WORKBENCH_DIR=" /export/HCP/workbench-v1.3.2"
13
+ DEFAULT_FIELDMAP_NUMBER=1
13
14
14
15
inform ()
15
16
{
@@ -32,7 +33,8 @@ get_options()
32
33
g_fsl_dir=" ${DEFAULT_FSL_DIR} "
33
34
g_freesurfer_dir=" ${DEFAULT_FREESURFER_DIR} "
34
35
g_workbench_dir=" ${DEFAULT_WORKBENCH_DIR} "
35
-
36
+ g_fieldmap_number=" ${DEFAULT_FIELDMAP_NUMBER} "
37
+
36
38
# parse arguments
37
39
local num_args=${# arguments[@]}
38
40
local argument
@@ -78,6 +80,10 @@ get_options()
78
80
g_workbench_dir=${argument/* =/ " " }
79
81
index=$(( index + 1 ))
80
82
;;
83
+ --fieldmap-number=* )
84
+ g_fieldmap_number=${argument/* =/ " " }
85
+ index=$(( index + 1 ))
86
+ ;;
81
87
* )
82
88
inform " Unrecognized option: ${argument} "
83
89
exit 1
@@ -146,6 +152,13 @@ get_options()
146
152
inform " Workbench dir: ${g_workbench_dir} "
147
153
fi
148
154
155
+ if [ -z " ${g_fieldmap_number} " ]; then
156
+ inform " --fieldmap-number= required"
157
+ error_count=$(( error_count + 1 ))
158
+ else
159
+ inform " Fieldmap number: ${g_fieldmap_number} "
160
+ fi
161
+
149
162
if [ ${error_count} -gt 0 ]; then
150
163
inform " ABORTING"
151
164
exit 1
@@ -233,8 +246,8 @@ ${g_run_dir}/StructuralPreprocessing.SINGULARITY_PROCESS \\
233
246
--fnirtconfig=T1_2_MNI152_2mm.cnf \\
234
247
--gdcoeffs=Prisma_3T_coeff_AS82.grad \\
235
248
--topupconfig=b02b0.cnf \\
236
- --se-phase-pos=${g_session} _SpinEchoFieldMap1_PA .nii.gz \\
237
- --se-phase-neg=${g_session} _SpinEchoFieldMap1_AP .nii.gz \\
249
+ --se-phase-pos=${g_session} _SpinEchoFieldMap ${g_fieldmap_number} _PA .nii.gz \\
250
+ --se-phase-neg=${g_session} _SpinEchoFieldMap ${g_fieldmap_number} _AP .nii.gz \\
238
251
239
252
EOF
240
253
# --processing-phase=PreFreeSurfer
0 commit comments