GeoPandas, rasterio, xarray, STAC and OSMnx workflows you can run end-to-end.
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.
Benchmark point-in-polygon joins, buffers, overlays and dissolves on 200,000 points with GeoPandas and Shapely 2, see why spatial indexes give a 100x+ speed-up, and compare Shapefile, GeoPackage, FlatGeobuf and GeoParquet.
Turn GeoDataFrames into shareable Leaflet web maps: a one-line choropleth with .explore(), then a hand-built Folium map with custom styling, tooltips, clustered markers, switchable basemaps and plugins — saved as a single HTML file.
Download Sheffield's walkable street network from OpenStreetMap with OSMnx, find the shortest walk across the city, draw 5–20 minute isochrones and map how far every street is from the nearest supermarket.
Search seven summers of Sentinel-2 imagery with pystac-client, stream only the bands you need from Cloud Optimised GeoTIFFs, mask clouds, compute NDWI and MNDWI and measure how much of Lake Mead dried up between 2019 and 2025.
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.
Load Natural Earth vector data with GeoPandas, reproject to an equal-area CRS, calculate population density, join cities to countries, dissolve to continents and draw publication-ready choropleths.