Skip to main content
SimilarExport ZIP
100%

Download data from WFS using GDAL Python bindings

✓ Published0🌍 Public
TThomasG77
Last edited Sep 24, 2024
Created on Sep 23, 2024

This example shows how to download vector data from a WFS (Web Feature Service) endpoint using GDAL’s Python bindings. It demonstrates bulk download of a full layer, plus filtered extracts by bounding box via the `-spat` option, and outputs to GeoPackage, CSV, and GeoJSON formats. The code explicitly configures WFS paging with `OGR_WFS_PAGING_ALLOWED` and `OGR_WFS_PAGE_SIZE`, and uses `gdal.OpenEx` and `gdal.VectorTranslate` to fetch data from the Géoportail WFS service. A separate call requests GeoJSON output directly through the service’s `outputFormat=json` parameter.

AI-generated description

Similar vizzes

Loading thumbnail…

OGC WFS consumption URLs from a service

This example demonstrates how to enumerate all downloadable data layers and their available output formats from an OGC Web Feature Service (WFS). Using the `@camptocamp/ogc-client` JavaScript library for the WFS endpoint, it retrieves the service metadata, layer types, and output formats, then constructs a complete GetFeature URL for every layer-format combination. The code also shows a Python equivalent using `OWSLib` to achieve the same enumeration, highlighting the interoperability of the two approaches. The visualization likely presents these generated URLs as a list or table, connecting users directly to data downloads.

TThomasG77
71% match
Loading thumbnail…

Gist 2bd25e856943977444bac9a5fd8c6408

This example shows a QGIS workflow for automatically populating the application with WMS, WMTS, and WFS service connections from the French IGN geoservices portal. The script scrapes the portal’s documentation page to find a CSV file containing service metadata, parses it using the `csv` module, and deduplicates entries via Python dictionaries. It then defines functions to add or remove connections in QGIS’s settings using `QgsSettings`, configuring parameters like URLs, versions, and axis orientation. The code demonstrates a practical approach to batch-importing OGC web services into a desktop GIS environment.

TThomasG77
68% match
Loading thumbnail…

Load GTFS as GPKG with GDAL 3.7

This example demonstrates converting an IDFM GTFS transit feed into a GeoPackage database using GDAL’s new GTFS driver. The Bash script first checks that GDAL 3.7 is installed, then uses `ogrinfo` with `-json` to list all available layers. For each layer, it runs `ogr2ogr` with the GeoPackage format, appending layers incrementally with the `-update` flag. The source data is loaded directly from a URL via `/vsicurl/`, and the script allows customizing input and output names through environment variables.

TThomasG77
66% match
Loading thumbnail…

Copie données WFS depuis le géoportail de l'urbanisme

This example automates the bulk download of France’s urban planning WFS data from the Géoportail de l’Urbanisme, specifically extracting the GPKG files. It shows how to list available remote directories via `lftp` over SFTP, identify the latest dated export, and mirror the chosen folder with parallel transfers. The workflow uses shell commands with `timeout`, `ssh`, and `lftp` to handle authentication and file retrieval. It relies on public credentials and a documented PDF manual, making the process reproducible for users needing recent, structured land-use datasets.

TThomasG77
65% match
Loading thumbnail…

GDAL VSICURL tests to remote BDTOPO Geopackage France

This example tests GDAL’s VSICURL virtual file system by streaming a remote GeoPackage of French administrative data, specifically focusing on mairie (town hall) geometries. It demonstrates two ogr2ogr commands that query the BDTOPO dataset using SQL joins to extract mairie footprints as multipolygons and their centroids as points, converting them to GeoJSON with EPSG:4326. The visualization shows how both polygon boundaries and point locations are derived from the same remote source, highlighting the ability to perform spatial SQL operations directly on cloud-hosted vector data without local downloads.

TThomasG77
65% match
Loading thumbnail…

Log related to answer https://gis.stackexchange.com/questions/446578/how-does-qgis-extract-features-from-the-xml-response-of-an-inspire-wfs-server/446580

This example shows a debug log from a GDAL/OGR-based workflow that extracts cadastral parcel features from an INSPIRE WFS server and writes them to a GeoPackage. The log details sequential GDAL API calls, including `GDALOpen`, `GDALDriver::Create`, and repeated HTTP `Fetch` requests, revealing how the WFS driver negotiates capabilities, disables multi-layer definitions, and paginates through `GetFeature` requests. It also demonstrates the internal use of an Expat XML parser for GML reading and a temporary in-memory file (`/vsimem/tempwfs_0x21f0b20/file.gml`) to process the downloaded data. The output shows the explicit `STARTINDEX` and `COUNT` parameters used to fetch batches of up to 10,000 features.

TThomasG77
62% match
Loading thumbnail…

Gist 8dafcac4cab52f5242b00c077430c08d

This example demonstrates how to select cadastral parcels from the IGN Cadastre API based on a drawn polygon’s extent. It shows the process of converting a GeoJSON polygon into an encoded URL parameter, then sending a GET request via Bash with `jq` or Python with `urllib.parse` and `urllib.request`. The API returns parcel geometries, which are saved as a GeoJSON file. The code also warns against using overly complex or overly large geometries due to URL length limits and server load.

TThomasG77
61% match
Loading thumbnail…

Gist 80f99a186a44e6fa1930f2a08618d498

This visualization displays a QGIS project containing world map layers—countries, states, provinces, and disputed borders—from the bundled `world_map.gpkg` GeoPackage, alongside JSON-derived point layers for contributors, hackfests, and user groups. Using the QGIS Python API, the script loads each vector layer via `QgsVectorLayer` with OGR drivers, adds them to the project through `QgsProject.instance().addMapLayers`, and enables a retro canvas style by setting `iface.mapCanvas().setProperty("retro", True)` before refreshing. The example demonstrates procedural map composition from built-in geodata files.

TThomasG77
61% match