Exemplo n.º 1
0
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
Exemplo n.º 2
0
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
Exemplo n.º 3
0
                                   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()
Exemplo n.º 4
0
        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")