5 posts
Make a 3D terrain map in Python with PyVista. Turn an Environment Agency LIDAR DEM into a StructuredGrid, drape a Sentinel-2 image on it, compare vertical exaggeration and lighting, drape a GPX route over Striding Edge, and render off-screen to PNG and a rotating MP4 on Windows, macOS or Linux.
Delineate a catchment and its stream network from a free DEM in Python with pysheds: fill pits and depressions, D8 flow direction, flow accumulation, snapping the pour point, catchment polygon and area, Strahler stream order and drainage density, all checked against the official NRFA boundary for the River Derwent at Yorkshire Bridge in the Peak District.
Stream a Copernicus 30 m DEM of the Edale valley, reproject it to British National Grid, scale it to a 150 mm print and write a watertight binary STL using only NumPy, with checks that prove the mesh is printable.
Read only the pixels you need from a cloud-hosted Copernicus DEM, reproject to British National Grid, derive hillshade, slope and aspect with NumPy, compute zonal statistics and draw an elevation profile — all in one script.
From a digital elevation model to a physical 3D print: preparing terrain in QGIS and exporting it for slicing.