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