def test_map_bounds():
locations = [
# North Pole
[90, -180],
[90, 0],
[90, 180],
[90, 179],
# Equator
[0, -180],
[0, 0],
[0, 180],
# South Pole
[-90, -180],
[-90, 0],
[-90, 180]
]
invalid_locations = [
[-90.1, 0],
[90.1, 0],
[0, -180.1],
[0, 180.1],
[0, 360],
]
height_info = HeightInfo()
for location in locations:
assert height_info.get_height(
*location)['altitude_m'] != height_info.NODATA
# out of bounds
for location in invalid_locations:
with pytest.raises(ValueError):
height_info.get_height(*location)
def test_check_for_metadata_get_min_height():
height_info = HeightInfo()
data = height_info.get_min_height(51, 7, 52, 8)
assert data['h_min'] != height_info.NODATA
assert isinstance(data['h_min'], float)
assert isinstance(data['location_min'], list)
assert data['counter_min'] >= 0
assert isinstance(data['source'], str)
assert isinstance(data['attributions'], list)
def test_check_for_metadata_get_height():
height_info = HeightInfo()
for location in [[53.57, 9.98], [52.51, 13.42], [47.94, 8.3], [-41, 172]]:
data = height_info.get_height(*location)
assert data['altitude_m'] != height_info.NODATA
assert isinstance(data['source'], str)
assert isinstance(data['attributions'], list)
assert data['distance_m'] >= 0
assert data['distance_m'] < 654.86
assert isinstance(data['lat_found'], float)
assert isinstance(data['lon_found'], float)
#!/usr/bin/python
# -*- coding: utf-8 -*-
import time
import json
from height_map.terr50 import Terrain50
from height_map.srtm1 import Srtm1
from height_map.dgm200 import Dgm200
from height_map.gebco import Gebco
from height_map.height_info import HeightInfo
dgm200 = Dgm200()
terr50 = Terrain50()
srtm1 = Srtm1()
gebco = Gebco()
height_info = HeightInfo()
def get_height(lat, lon):
srtm1_result = srtm1.get_height(lat, lon)
terr50_result = terr50.get_height(lat, lon)
dgm200_result = dgm200.get_height(lat, lon)
gebco_result = gebco.get_height(lat, lon)
results = {'request': {'lat': lat, 'lon': lon}}
if srtm1_result['altitude_m'] != srtm1.NODATA:
results['SRTM1'] = srtm1_result
if terr50_result['altitude_m'] != terr50.NODATA:
results['TERR50'] = terr50_result
if dgm200_result['altitude_m'] != dgm200.NODATA:
results['DGM200'] = dgm200_result
import logging
import numpy as np
from pydantic import BaseModel, confloat, constr, conlist
from typing import List
from geojson import FeatureCollection, Feature, LineString
import pygeodesy.ellipsoidalVincenty as eV
from simplification.cutil import simplify_coords
from height_map.height_info import HeightInfo
from height_map.cci_land_cover import LandCover
from height_map.timeit import timeit
logger = logging.getLogger(__name__)
hi = HeightInfo()
lc = LandCover()
class Location(BaseModel):
lat: confloat(ge=-90, le=90)
lon: confloat(ge=-180, le=180)
class PositionRequest(BaseModel):
track: List[Location]
distance: confloat(ge=0)
@timeit
def get_track_length(track: List[Location]):
distance = 0
old_location = None
for _location in track:
def test_content_get_height():
height_info = HeightInfo()
# Hambach open pit
assert math.isclose(height_info.get_height(50.91, 6.51)['altitude_m'],
-94.83,
abs_tol=55)
# storage pool in Geeste
assert math.isclose(height_info.get_height(52.588, 7.294)['altitude_m'],
34,
abs_tol=16)
# Black Forest
assert math.isclose(height_info.get_height(47.94, 8.3)['altitude_m'],
927.52,
abs_tol=16)
# farmland in the Emsland region
assert math.isclose(height_info.get_height(52.78, 7.4)['altitude_m'],
32.61,
abs_tol=16)
# high moor in the Emsland region
assert math.isclose(height_info.get_height(52.8, 7.4)['altitude_m'],
24.64,
abs_tol=16)
# lowest location in Germany
assert math.isclose(height_info.get_height(50.925084,
6.529813)['altitude_m'],
-265.5,
abs_tol=16)
# highest location in the bounding box of Germany
assert math.isclose(height_info.get_height(47.422239,
10.986123)['altitude_m'],
2920.32,
abs_tol=16)
# lowest location
assert math.isclose(height_info.get_height(11.366667,
142.5875)['altitude_m'],
-10952,
abs_tol=16)
# highest location
assert math.isclose(height_info.get_height(27.9875, 86.925)['altitude_m'],
8613.2,
abs_tol=16)
# London, River Thames
assert math.isclose(height_info.get_height(51.499, -0.122)['altitude_m'],
6.98,
abs_tol=16)
# London, The Serpentine
assert math.isclose(height_info.get_height(51.5052, -0.1666)['altitude_m'],
20.96,
abs_tol=16)
# London, Round Pond
assert math.isclose(height_info.get_height(51.5058, -0.1834)['altitude_m'],
28.96,
abs_tol=16)
# London, Hyde Park
assert math.isclose(height_info.get_height(51.509, -0.163)['altitude_m'],
28.96,
abs_tol=16)
# London, Kensington Gardens
assert math.isclose(height_info.get_height(51.507, -0.179)['altitude_m'],
25.95,
abs_tol=16)
# North Sea
assert math.isclose(height_info.get_height(53.8, 6.9)['altitude_m'],
-22.04,
abs_tol=16)
# Dead Sea
assert math.isclose(height_info.get_height(31.52, 35.48)['altitude_m'],
-415,
abs_tol=16)
# IJsselmeer
assert math.isclose(height_info.get_height(52.74, 5.42)['altitude_m'],
-4.28,
abs_tol=16)
# Markermeer
assert math.isclose(height_info.get_height(52.54, 5.22)['altitude_m'],
-3.89,
abs_tol=16)
# Naples harbour
assert math.isclose(height_info.get_height(40.836, 14.262)['altitude_m'],
-14.21,
abs_tol=16)
# Lake Constance
assert math.isclose(height_info.get_height(47.56, 9.5)['altitude_m'],
392,
abs_tol=16)
# Salt Lake City
assert math.isclose(height_info.get_height(40.8, -112)['altitude_m'],
1289.59,
abs_tol=16)
# Tokyo
assert math.isclose(height_info.get_height(35.68, 139.77)['altitude_m'],
20.31,
abs_tol=16)
def test_check_bounds_get_min_max_height():
height_info = HeightInfo()
# in bounds
_result = height_info.get_min_max_height(52.5, 5, 53, 5.5)
assert _result['h_max'] != height_info.NODATA
assert _result['h_min'] != height_info.NODATA
# invalid coordinates
with pytest.raises(ValueError):
height_info.get_min_max_height(52.5, 5, 91, 5.5)
with pytest.raises(ValueError):
height_info.get_min_max_height(-91, 5, 53, 5.5)
with pytest.raises(ValueError):
height_info.get_min_max_height(52.5, -181, 53, 5.5)
with pytest.raises(ValueError):
height_info.get_min_max_height(52.5, 5, 53, 181)
# incorrect rectangle
_result = height_info.get_min_max_height(54.886907, 15.570925, 47.240591,
6.093066)
assert _result['h_max'] == height_info.NODATA
assert _result['h_min'] == height_info.NODATA
# highest and lowest location in the bounding box of Germany
_result = height_info.get_min_max_height(47.240591, 6.093066, 54.886907,
15.570925)
assert math.isclose(_result['h_max'], 2920, abs_tol=16)
assert math.isclose(_result['h_min'], -265.5, abs_tol=16)
# highest and lowest location of the world
_result = height_info.get_min_max_height(-90, -180, 90, 180)
assert math.isclose(_result['h_max'], 8613.2, abs_tol=16)
assert math.isclose(_result['h_min'], -10952, abs_tol=16)