GeoData DareData Blog

reading_raster_data_and_getting_values.py

import rasterio
from rasterio.plot import show
 
url = "zip+file:data/mdt.zip!mdt.tif"
lisbon_elevation = rasterio.open(url)

# Plot the raster data to get a sense of it
show(lisbon_elevation, cmap="terrain")

# Get the elevation from the raster data
house_data_gdf["elevation"] = (
    house_data_gdf
    .to_crs(lisbon_elevation.crs)
    .geometry
    .apply(
        lambda x: next(lisbon_elevation.sample([(x.x, x.y)]))[0]
    )
)