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
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