Skip to content

Instantly share code, notes, and snippets.

@lrq3000
Forked from MarcCote/trk2tck.py
Created December 1, 2018 23:55
Show Gist options
  • Star 0 You must be signed in to star a gist
  • Fork 0 You must be signed in to fork a gist
  • Save lrq3000/fee2b6ab26f5cbfdf5d5f6917f62a884 to your computer and use it in GitHub Desktop.
Save lrq3000/fee2b6ab26f5cbfdf5d5f6917f62a884 to your computer and use it in GitHub Desktop.
Script that converts TRK to TCK using https://github.com/MarcCote/nibabel/tree/streamlines_tck
import os
import argparse
import nibabel as nib
def build_argparser():
DESCRIPTION = "Convert tractograms (TRK -> TCK)."
p = argparse.ArgumentParser(description=DESCRIPTION)
p.add_argument('tractograms', metavar='bundle', nargs="+", help='list of tractograms.')
p.add_argument('-f', '--force', action="store_true", help='overwrite existing output files.')
return p
def main():
parser = build_argparser()
args = parser.parse_args()
for tractogram in args.tractograms:
if nib.streamlines.detect_format(tractogram) is not nib.streamlines.TrkFile:
print("Skipping non TRK file: '{}'".format(tractogram))
continue
output_filename = tractogram[:-4] + '.tck'
if os.path.isfile(output_filename) and not args.force:
print("Skipping existing file: '{}'. Use -f to overwrite.".format(output_filename))
continue
trk = nib.streamlines.load(tractogram)
nib.streamlines.save(trk.tractogram, output_filename)
if __name__ == '__main__':
main()
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment