def main(args):
if args.outdir=='GT':
marker_files = map(string.strip,os.popen('ls '+args.indir+'*-GT.marker').readlines())
substack_ids = [f.split('/')[-1].split('-')[0] for f in marker_files]
else:
marker_files = map(string.strip,os.popen('ls '+args.outdir+'/*/ms.marker').readlines())
substack_ids = [f.split('/')[-2] for f in marker_files]
C_final=set()
plist = None
hue = 0.0
# FIXME: This is way too slow, should use divide & conquer using a binary tree
for marker_file,substack_id in zip(marker_files,substack_ids):
substack=SubStack(args.indir,substack_id, plist)
plist = substack.plist
print('Merging', marker_file, substack_id)
C=substack.load_markers(marker_file,args.outdir=='GT')
# I had forgotten that markers are shifted by 1 in save_markers()
# that was because Vaa3D coordinates starts from 1 rather than 0
# ==> repair original values
for c in C:
c.x -= 1
c.y -= 1
c.z -= 1
hue = hue + 0.31
if hue > 1:
hue = hue -1
merge(C_final, C, substack, hue, args.verbose)
substack.save_markers(args.outfile, C_final)
def main(args):
if args.outdir == 'GT':
marker_files = map(
string.strip,
os.popen('ls ' + args.indir + '*-GT.marker').readlines())
substack_ids = [f.split('/')[-1].split('-')[0] for f in marker_files]
else:
marker_files = map(
string.strip,
os.popen('ls ' + args.outdir + '/*/ms.marker').readlines())
substack_ids = [f.split('/')[-2] for f in marker_files]
C_final = set()
plist = None
hue = 0.0
# FIXME: This is way too slow, should use divide & conquer using a binary tree
for marker_file, substack_id in zip(marker_files, substack_ids):
substack = SubStack(args.indir, substack_id, plist)
plist = substack.plist
print('Merging', marker_file, substack_id)
C = substack.load_markers(marker_file, args.outdir == 'GT')
# I had forgotten that markers are shifted by 1 in save_markers()
# that was because Vaa3D coordinates starts from 1 rather than 0
# ==> repair original values
for c in C:
c.x -= 1
c.y -= 1
c.z -= 1
hue = hue + 0.31
if hue > 1:
hue = hue - 1
merge(C_final, C, substack, hue, args.verbose)
substack.save_markers(args.outfile, C_final)
