Exemplo n.º 1
0
def test_xy_array_bilinear_window():
    """ integration test
    """
    x, y = (245309, 1000064)

    with rasterio.open(raster) as src:
        win, unitxy = point_window_unitxy(x, y, affine)
        arr = src.read(1, window=win)

    val = bilinear(arr, *unitxy)
    assert round(val) == 74
Exemplo n.º 2
0
def test_bilinear():
    import numpy as np
    arr = np.array([[1.0, 2.0],
                    [3.0, 4.0]])

    assert bilinear(arr, 0, 0) == 3.0
    assert bilinear(arr, 1, 0) == 4.0
    assert bilinear(arr, 1, 1) == 2.0
    assert bilinear(arr, 0, 1) == 1.0
    assert bilinear(arr, 0.5, 0.5) == arr.mean()
    assert bilinear(arr, 0.95, 0.95) < 4.0
    assert bilinear(arr, 0.05, 0.95) > 1.0