Example #1
0
def graph_rgb(view_xy):
    view = gview.app.view_manager.get_active_view_window()
    layers = view.viewarea.list_layers()
    to_plot = []
    for layer in layers:
        lclass = ev_mosaic_viewer.layer_class(layer)
        raster = layer.get_data(0)
        date = get_dataset_date(raster.get_dataset())
        if date is None:
            continue
        
        pixel_xy = layer.view_to_pixel(view_xy[0], view_xy[1])

        data_point = [date]
        
        for comp in range(3):
            raster = layer.get_data(comp)
            data_point.append(raster.get_sample(pixel_xy[0], pixel_xy[1]))
        
        print raster.get_dataset().GetDescription(), data_point

        to_plot.append(data_point)

    to_plot = Numeric.array(to_plot)
    
    return gvplot.plot(to_plot,
                       xaxis = 'Days',
                       yaxis = 'Value',
                       ymin = 6500.0,
                       ymax = 11000.0,
                       terminal='openev',
                       multiplot=True)
Example #2
0
def graph_rgb(view_xy):
    view = gview.app.view_manager.get_active_view_window()
    layers = view.viewarea.list_layers()
    to_plot = []
    for layer in layers:
        lclass = ev_mosaic_viewer.layer_class(layer)
        raster = layer.get_data(0)
        date = get_dataset_date(raster.get_dataset())
        if date is None:
            continue

        pixel_xy = layer.view_to_pixel(view_xy[0], view_xy[1])

        data_point = [date]

        for comp in range(3):
            raster = layer.get_data(comp)
            data_point.append(raster.get_sample(pixel_xy[0], pixel_xy[1]))

        print raster.get_dataset().GetDescription(), data_point

        to_plot.append(data_point)

    to_plot = Numeric.array(to_plot)

    return gvplot.plot(to_plot,
                       xaxis='Days',
                       yaxis='Value',
                       ymin=6500.0,
                       ymax=11000.0,
                       terminal='openev',
                       multiplot=True)
Example #3
0
def graph_darkness(view_xy):
    view = gview.app.view_manager.get_active_view_window()
    layers = view.viewarea.list_layers()
    to_plot = []
    for layer in layers:
        lclass = ev_mosaic_viewer.layer_class(layer)
        if lclass != ev_mosaic_viewer.LYR_LANDSAT8:
            continue
        
        raster = layer.get_data(0)

        # ignore rasters that aren't source scenes with a date.
        date = get_dataset_date(raster.get_dataset())
        if date is None:
            continue
        
        pixel_xy = layer.view_to_pixel(view_xy[0], view_xy[1])

        data_point = [date]

        for xoff in range(-2,4,2):
            for yoff in range(-2,4,2):
                rgb_sum = 0
                for comp in range(3):
                    raster = layer.get_data(comp)
                    rgb_sum += raster.get_sample(pixel_xy[0] + xoff,
                                                 pixel_xy[1] + yoff)

                rgb_sum *= 0.3333

                to_plot.append(rgb_sum)

    to_plot.sort()
    
    return gvplot.plot(to_plot,
                       xaxis = 'Sample',
                       yaxis = 'Brightness',
                       ymin = 6500.0,
                       ymax = 11000.0,
                       terminal='openev')
Example #4
0
def graph_darkness(view_xy):
    view = gview.app.view_manager.get_active_view_window()
    layers = view.viewarea.list_layers()
    to_plot = []
    for layer in layers:
        lclass = ev_mosaic_viewer.layer_class(layer)
        if lclass != ev_mosaic_viewer.LYR_LANDSAT8:
            continue

        raster = layer.get_data(0)

        # ignore rasters that aren't source scenes with a date.
        date = get_dataset_date(raster.get_dataset())
        if date is None:
            continue

        pixel_xy = layer.view_to_pixel(view_xy[0], view_xy[1])

        data_point = [date]

        for xoff in range(-2, 4, 2):
            for yoff in range(-2, 4, 2):
                rgb_sum = 0
                for comp in range(3):
                    raster = layer.get_data(comp)
                    rgb_sum += raster.get_sample(pixel_xy[0] + xoff,
                                                 pixel_xy[1] + yoff)

                rgb_sum *= 0.3333

                to_plot.append(rgb_sum)

    to_plot.sort()

    return gvplot.plot(to_plot,
                       xaxis='Sample',
                       yaxis='Brightness',
                       ymin=6500.0,
                       ymax=11000.0,
                       terminal='openev')