Exemplo n.º 1
0
def test_read_ds_data():
    expected = 151.0
    gdal.AllRegister()
    data_source = gdal.Open(DS_FILENAME, GA_ReadOnly)
    actual = geods.read_ds_data(data_source, 529, 477)

    assert abs(expected - actual) <= EPSILON
Exemplo n.º 2
0
def test_read_ds_data_np():
    expected = np.array([195, 182, 176, 177, 175, 160, 136, 130, 109, 113])
    gdal.AllRegister()
    data_source = gdal.Open(DS_FILENAME, GA_ReadOnly)
    actual = geods.read_ds_data(data_source, np.linspace(300, 400, 10, dtype=int),
                                np.linspace(400, 300, 10, dtype=int))

    for exp, act in zip(expected, actual):
        assert abs(exp - act) <= EPSILON