def test_xy_array_bilinear_window():
    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
integration test
dataset/ETHPy150Open perrygeo/python-rasterstats/tests/test_point.py/test_xy_array_bilinear_window