Пример #1
0
def test_graph_from_image(visual=False):
    im = imread("segmentation.inr.gz")
    graph = graph_from_image(im)
    if visual:
        Viewer.display(
            graph2pglscene(graph, graph.vertex_property('barycenter'),
                           graph.vertex_property('L1')))
Пример #2
0
def test_graph_from_simple_image(visual=False):
    im = np.array([[1, 2, 7, 7, 1, 1], [1, 6, 5, 7, 3, 3], [2, 2, 1, 7, 3, 3],
                   [1, 1, 1, 4, 1, 1]])
    borders = SpatialImageAnalysis(im).cells_in_image_margins()
    print '**', SpatialImageAnalysis(im).boundingbox(real=False)
    #~ print borders
    #~ print np.unique(im)
    graph = graph_from_image(im)

    #~ print list(graph.vertices())
    #~ print map(graph.edge_vertices,graph.edges())
    #~ print list(graph.vertex_property_names())
    #~ print list(graph.edge_property_names())
    for propname in graph.vertex_property_names():
        print propname
        print graph.vertex_property(propname)

    #~ for propname in graph.edge_property_names():
    #~ print propname
    #~ print [(i,graph.edge_vertices(i),j) for i,j in graph.edge_property(propname).iteritems()]
    if visual:
        Viewer.display(
            graph2pglscene(graph, graph.vertex_property('barycenter'),
                           graph.vertex_property('border')))
def test_graph_from_simple_image(visual = False):
    im = np.array([[1, 2, 7, 7, 1, 1],
                  [1, 6, 5, 7, 3, 3],
                  [2, 2, 1, 7, 3, 3],
                  [1, 1, 1, 4, 1, 1]])
    borders = SpatialImageAnalysis(im).cells_in_image_margins()
    print '**', SpatialImageAnalysis(im).boundingbox(real=False)
    #~ print borders
    #~ print np.unique(im)
    graph = graph_from_image(im)
    
    #~ print list(graph.vertices())
    #~ print map(graph.edge_vertices,graph.edges())
    #~ print list(graph.vertex_property_names())
    #~ print list(graph.edge_property_names())
    for propname in graph.vertex_property_names():
        print propname
        print graph.vertex_property(propname)

    #~ for propname in graph.edge_property_names():
        #~ print propname
        #~ print [(i,graph.edge_vertices(i),j) for i,j in graph.edge_property(propname).iteritems()]
    if visual :
        Viewer.display(graph2pglscene(graph,graph.vertex_property('barycenter'),graph.vertex_property('border')))
def test_graph_from_image(visual = False):
    im =  imread("segmentation.inr.gz")
    graph = graph_from_image(im)
    if visual :
        Viewer.display(graph2pglscene(graph,graph.vertex_property('barycenter'),graph.vertex_property('L1')))
Пример #5
0
 def plot(self):
     from openalea.plantgl.gui import Viewer
     Viewer.display(self.toScene())