def test_utm(self, patch_config):
lats = [0.2, 0.8, 0.6]
lons = [10.2, 10.8, 11.5]
dataset = config.load_datasets()[SRTM_DATASET_NAME]
z = backend.get_elevation(lats, lons, dataset)
dataset_utm = config.load_datasets()[SRTM_UTM_DATASET_NAME]
z_utm = backend.get_elevation(lats, lons, dataset_utm)
assert np.allclose(z, z_utm)
def test_dataset_crs_format_equivalence(self, patch_config):
lats = [43.597009, 45.534601]
lons = [1.455697, 10.698698]
dataset_epsg = config.load_datasets()[EU_DEM_DATASET_NAME]
z_epsg = backend.get_elevation(lats, lons, dataset_epsg)
dataset_wkt = config.load_datasets()[EU_DEM_NO_EPSG_DATASET_NAME]
z_wkt = backend.get_elevation(lats, lons, dataset_wkt)
assert np.allclose(z_epsg, z_wkt)
def test_multi_dataset_mask(self, patch_config):
lats = [47.625765, 0.1, 70]
lons = [9.418759, 10.5, 150]
datasets = [
config.load_datasets()[ETOPO1_RESAMPLED_DATASET_NAME],
config.load_datasets()[EU_DEM_DATASET_NAME],
config.load_datasets()[SRTM_DATASET_NAME],
]
z, dataset_names = backend.get_elevation(lats, lons, datasets)
assert all(np.isfinite(z))
assert dataset_names == [ETOPO1_RESAMPLED_DATASET_NAME] * len(lats)
def _load_datasets():
"""Load datasets defined in config
Returns:
Dict of {dataset_name: config.Dataset object} items.
"""
return config.load_datasets()
def test_srtm_tiles(self, patch_config):
lats = [0.1, 0.9]
lons = [10.5, 11.5]
dataset = config.load_datasets()[SRTM_DATASET_NAME]
z = backend.get_elevation(lats, lons, dataset)
assert all(z)
assert all(np.isfinite(z))
def test_single_dataset(self, patch_config):
lats = [0.1, 0.9]
lons = [10.5, 11.5]
dataset = config.load_datasets()[SRTM_DATASET_NAME]
z, names = backend.get_elevation(lats, lons, [dataset])
assert all(z)
assert all(np.isfinite(z))
assert names == [SRTM_DATASET_NAME] * len(lats)
def test_single_file_dataset(self):
lats = [0.1, -9]
lons = [-50.5, 12.11]
interpolation = "cubic"
dataset = config.load_datasets()[ETOPO1_DATASET_NAME]
elevations_by_dataset = backend.get_elevation(lats, lons, dataset,
interpolation)
elevations_by_path = backend._get_elevation_from_path(
lats, lons, ETOPO1_GEOTIFF_PATH, interpolation)
assert elevations_by_dataset == elevations_by_path
def _load_datasets_memcache():
return config.load_datasets()
def test(self):
config.load_datasets()
def test_alternate_tiled_dataset(self, patch_config):
lats = [47.625765]
lons = [9.418759]
dataset = config.load_datasets()[EU_DEM_DATASET_NAME]
z = backend.get_elevation(lats, lons, dataset)
assert np.isfinite(z)
def test_out_of_srtm_bounds_utm(self, patch_config):
lats = [70]
lons = [10.5]
dataset = config.load_datasets()[SRTM_DATASET_NAME]
z = backend.get_elevation(lats, lons, dataset)[0]
assert z is None
def test_oob(self, patch_config):
lats = [1.5, -0.5, 0.5, 0.5]
lons = [10.5, 11.5, 9.5, 12.5]
dataset = config.load_datasets()[SRTM_DATASET_NAME]
z = backend.get_elevation(lats, lons, dataset)
assert all([x is None for x in z])
