def loaddata(fn,slave=True): """" load a dataset from filename fn, use pickled data if slave=True. for each dataset, this has to be called at least once with slave=False to create the pickled version. """ print "Loading dataset `%s'"%fn D = dataset("%s.dat"%fn,crop=False,upsample="zoom",remove_rohr=True,usepickled=slave) return D
def loaddata(fn, slave=True): """" load a dataset from filename fn, use pickled data if slave=True. for each dataset, this has to be called at least once with slave=False to create the pickled version. """ print "Loading dataset `%s'" % fn D = dataset("%s.dat" % fn, crop=False, upsample="zoom", remove_rohr=True, usepickled=slave) return D
color=(.1, .1, 1), dumb=True) viewer.show_points_onefile(basename + "/missing.txt", dscale=0.5, what="wireframe", color=(1, .3, .3), dumb=True, opacity=0.3) if offscreen: mkimg(fig, "wireframe_diff") mlab.clf() if "satoiso" in token: Dsato = dataset(basename + "/sato.dat", upsample=None, crop=False, usepickled=False, medianfilt=False).D viewer.show_iso(Dsato, 4, "Spectral", 0.2) if offscreen: mkimg(fig, "satoiso") mlab.clf() if "satovol" in token: Dsato = dataset(basename + "/sato.dat", dtype=np.float32).D #viewer.show_volume(Dsato, "bone", 0.00001, 0.01) # barley viewer.show_volume(Dsato, "bone", 0.01, 0.04) # lupine if offscreen: mkimg(fig, "satovol") mlab.clf()
mkimg(fig, "wireframe") mlab.clf() if "wireframe_diff" in token: gt_name = basename.split("roots")[0] viewer.show_points( gt_name+"/vertices.txt", gt_name+"/edges.txt", dscale=0.5, what="wireframe", opacity=0.2,color=(0.8,0.3,0.3)) viewer.show_points( basename+"/vertices.txt", basename+"/edges.txt", dscale=0.5*.85, what="wireframe", opacity=0.2,color=(0.1,0.1,0.8)) viewer.show_points_onefile( basename+"/toomuch.txt", dscale=0.5, what="wireframe", color=(.1,.1,1),dumb=True) viewer.show_points_onefile( basename+"/missing.txt", dscale=0.5, what="wireframe", color=(1,.3,.3),dumb=True, opacity=0.3) if offscreen: mkimg(fig, "wireframe_diff") mlab.clf() if "satoiso" in token: Dsato = dataset(basename+"/sato.dat", upsample=None, crop=False,usepickled=False,medianfilt=False).D viewer.show_iso(Dsato, 4, "Spectral", 0.2) if offscreen: mkimg(fig, "satoiso") mlab.clf() if "satovol" in token: Dsato = dataset(basename+"/sato.dat",dtype=np.float32).D #viewer.show_volume(Dsato, "bone", 0.00001, 0.01) # barley viewer.show_volume(Dsato, "bone", 0.01, 0.04) # lupine if offscreen: mkimg(fig, "satovol") mlab.clf() if "dmap" in token: info = WurzelInfo(basename+".dat")