I have noticed a bug in seeds to targets mode of probtrackx:
The seeds_to_targets files do not have proper NIFTI headers.
When running something like:
probtrackx --seed=<image.nii.gz> --targetmasks=<targetlist.txt>
--dir=<results folder> --samples=<merged> --mask=<nodif_brain_mask>
--stop=<stop mask> --xfm=<xfm> --loopcheck --cthr=0.2 --nsteps=2000
--steplength=0.5 --nsamples=25000 --fibthresh=0.05 --opd --os2t --fibst=1
--modeuler –forcedir
There is output a fdt_paths.nii.gz, which contains all of
the samples that were sent out, and a number of seeds to targets files as
expected. The fdt_paths.nii.gz has an appropriate header that matches the
seed’s header, e.g.:
filename
fdt_paths.nii.gz
sizeof_hdr 348
data_type INT32
dim0
3
dim1
58
dim2
69
dim3
58
dim4
1
dim5
1
dim6
1
dim7
1
vox_units mm
time_units s
datatype 8
nbyper 4
bitpix 32
pixdim0
0.0000000000
pixdim1
3.0000000000
pixdim2
3.0000000000
pixdim3
3.0000000000
pixdim4
1.0000000000
pixdim5
0.0000000000
pixdim6
0.0000000000
pixdim7
0.0000000000
vox_offset 352
cal_max 4.0000
cal_min 0.0000
scl_slope 0.000000
scl_inter 0.000000
phase_dim 0
freq_dim 0
slice_dim 0
slice_name Unknown
slice_code 0
slice_start 0
slice_end 0
slice_duration 0.000000
time_offset 0.000000
intent
Unknown
intent_code 0
intent_name
intent_p1 0.000000
intent_p2 0.000000
intent_p3 0.000000
qform_name Scanner Anat
qform_code 1
qto_xyz:1 3.000000 0.000000
0.000000 -84.500000
qto_xyz:2 0.000000
3.000000 0.000000 -123.500000
qto_xyz:3 0.000000
0.000000 3.000000 -74.500000
qto_xyz:4 0.000000
0.000000 0.000000 1.000000
qform_xorient Left-to-Right
qform_yorient Posterior-to-Anterior
qform_zorient Inferior-to-Superior
sform_name Scanner Anat
sform_code 1
sto_xyz:1 3.000000
0.000000 0.000000 -84.500000
sto_xyz:2 0.000000
3.000000 0.000000 -123.500000
sto_xyz:3 0.000000 0.000000
3.000000 -74.500000
sto_xyz:4 0.000000
0.000000 0.000000 1.000000
sform_xorient Left-to-Right
sform_yorient Posterior-to-Anterior
sform_zorient Inferior-to-Superior
file_type NIFTI-1+
file_code 1
descrip FSL4.0
aux_file
The seeds_to_targets files do not have proper headers, e.g.:
filename
seeds_to_2.nii.gz
sizeof_hdr 348
data_type FLOAT32
dim0
3
dim1
58
dim2
69
dim3
58
dim4
1
dim5
1
dim6
1
dim7
1
vox_units mm
time_units s
datatype 16
nbyper 4
bitpix 32
pixdim0
0.0000000000
pixdim1
1.0000000000
pixdim2
1.0000000000
pixdim3
1.0000000000
pixdim4 1.0000000000
pixdim5
0.0000000000
pixdim6
0.0000000000
pixdim7
0.0000000000
vox_offset 352
cal_max 0.0000
cal_min 0.0000
scl_slope 0.000000
scl_inter 0.000000
phase_dim 0
freq_dim 0
slice_dim 0
slice_name Unknown
slice_code 0
slice_start 0
slice_end 0
slice_duration 0.000000
time_offset 0.000000
intent
Unknown
intent_code 0
intent_name
intent_p1 0.000000
intent_p2 0.000000
intent_p3 0.000000
qform_name Unknown
qform_code 0
qto_xyz:1 1.000000
0.000000 0.000000 0.000000
qto_xyz:2 0.000000
1.000000 0.000000 0.000000
qto_xyz:3 0.000000
0.000000 1.000000 0.000000
qto_xyz:4 0.000000
0.000000 0.000000 1.000000
qform_xorient Left-to-Right
qform_yorient Posterior-to-Anterior
qform_zorient Inferior-to-Superior
sform_name Unknown
sform_code 0
sto_xyz:1 0.000000
0.000000 0.000000 0.000000
sto_xyz:2 0.000000
0.000000 0.000000 0.000000
sto_xyz:3 0.000000
0.000000 0.000000 0.000000
sto_xyz:4 0.000000
0.000000 0.000000 0.000000
sform_xorient Unknown
sform_yorient Unknown
sform_zorient Unknown
file_type NIFTI-1+
file_code 1
descrip FSL4.0
aux_file
Basically, the pixel dimensions are not set, and neither are
the qforms or sforms.
Thanks,
Matt.