else:
aparc = bg.copy_file(aparc, outdirs["anatomydir"][0])
aparcnii = bg.convert(aparc, aparcnii)
bg.remove_files([aparc])
# make pons
refdir, _ = outdirs["refdir"]
brainstem = os.path.join(refdir, "brainstem.nii.gz")
if os.path.isfile(brainstem):
logging.warning("brainstem %s exists, skipping" % (brainstem))
else:
# copy aseg+aparc to refdir
try:
caparcnii = bg.copy_file(aparcnii, refdir)
bg.make_brainstem(caparcnii)
brainstem = bg.unzip_file(brainstem)
bg.remove_files([caparcnii.replace(".gz", "")])
except:
logging.warning("Check %s" % (brainstem))
rawtracer, exists = outdirs["rawtracer"]
rawtracer_base, _ = os.path.split(rawtracer)
os.system("rm -rf %s" % rawtracer_base)
tracerdir, _ = outdirs["tracerdir"]
newname = "%s_%s" % (subid, tracer)
copied_ecats = bg.copy_files(ecats, tracerdir)
bg.convertallecat(copied_ecats, newname)
logging.info("ecats converted for %s " % (subid))
logging.error('%s NO APARC ASEG: %s' % (subid, aparc))
else:
aparc = bg.copy_file(aparc, outdirs['anatomydir'][0])
aparcnii = bg.convert(aparc, aparcnii)
bg.remove_files([aparc])
# make cerebellum
refdir, _ = outdirs['refdir']
cerebellum = os.path.join(refdir, 'grey_cerebellum.nii')
if os.path.isfile(cerebellum):
logging.warning('cerebellum %s exists, skipping' % (cerebellum))
else:
# copy aseg+aparc to refdir
try:
caparcnii = bg.copy_file(aparcnii, refdir)
bg.make_cerebellum_nibabel(caparcnii)
bg.remove_files([caparcnii])
except:
logging.warning('Fail: unable to make %s' % (cerebellum))
rawtracer, exists = outdirs['rawtracer']
os.system('rm %s' % rawtracer)
tracerdir, _ = outdirs['tracerdir']
newname = '%s_%s' % (subid, tracer)
copied_ecats = bg.copy_files(ecats, tracerdir)
bg.convertallecat(copied_ecats, newname)
logging.info('ecats converted for %s ' % (subid))
bg.remove_files([aparc])
# make cerebellum
refdir, _ = outdirs['refdir']
cerebellum = os.path.join(refdir, 'grey_cerebellum.nii')
if os.path.isfile(cerebellum):
logging.warning('cerebellum %s exists, skipping' % (cerebellum))
continue
else:
# copy aseg+aparc to refdir
try:
caparcnii = bg.copy_file(aparcnii, refdir)
bg.make_cerebellum_nibabel(caparcnii)
bg.remove_files([caparcnii])
except:
logging.warning('Fail: unable to make %s' % (cerebellum))
rawtracer, exists = outdirs['rawtracer']
tracerdir, _ = outdirs['tracerdir']
if exists:
logging.info('RAW %s data exists,'\
'remove %s to rerun'%(subid,rawtracer))
continue
else:
# Copy PET data, convert to nifti
newraw = bg.copy_dir(petdir, rawtracer)
ecats = bg.copy_dir(rawtracer, tracerdir, pattern='*.v')
newname = '%s_%s' % (subid, tracer)
bg.convertallecat(ecats, newname)
logging.info('ecats converted for %s ' % (subid))
# make brainstem, find pons
refdir,_ = outdirs['refdir']
pons = os.path.join(refdir, 'pons_tu.nii')
if os.path.isfile(pons):
logging.warning('pons %s exists, skipping'%(pons))
continue
else:
# copy aseg+aparc to refdir
try:
caparcnii = bg.copy_file(aparcnii, refdir)
bg.make_brainstem(caparcnii)
bg.remove_files([caparcnii])
except:
logging.warning('Fail: unable to make %s'%(pons))
rawtracer, exists = outdirs['rawtracer']
tracerdir, _ = outdirs['tracerdir']
if exists:
logging.info('RAW %s data exists,'\
'remove %s to rerun'%(subid,rawtracer))
continue
else:
# Copy PET data, convert to nifti
newraw = bg.copy_dir(petdir, rawtracer)
ecats = bg.copy_dir(rawtracer, tracerdir, pattern='*.v')
newname = '%s_%s' % (subid, tracer)
bg.convertallecat(ecats, newname)
logging.info('ecats converted for %s ' % (subid))
