Getting Elevation Info From the Austrian Elevation Service

Binder

Homepage of the service: https://maegger.github.io/getAustrianElevation.html (Copyright (c) 2017, Manfred Egger)

Related QGIS plugin: https://github.com/maegger/AustrianElevation

Elevation data source: CC BY 3.0 AT http://geoland.at/

In [ ]:
import hvplot.pandas
from geopy.geocoders import Nominatim
from utils.dataaccess import get_elevation
from utils.converting import location_to_gdf
In [ ]:
address = "Stephansdom, Wien"
locator = Nominatim(user_agent="OGD.AT-Lab")
location = locator.geocode(address)
print(location.address)
print("Latitude = {}, Longitude = {}".format(location.latitude, location.longitude))
gdf = location_to_gdf(location, address)

Before we can query the elevation, we need to reproject the coordinates to EPSG:3857

In [ ]:
gdf = gdf.to_crs('epsg:3857')
gdf
In [ ]:
gdf.loc[0, 'elevation'] = get_elevation(gdf.iloc[0].geometry)
gdf.to_crs('epsg:4326')
In [ ]: