Gist c32a8597373bd00c10029d4f21f2400e

100%

Gist c32a8597373bd00c10029d4f21f2400e

✓ Published0🌍 Public

Last edited Jun 28, 2021

Created on Jun 28, 2021

This example demonstrates how to convert a GDAL VRT (Virtual Dataset) into a single raster file using the `ogr2vrt` utility. It shows a command-line workflow where a shell script first installs the `gdal-utils` package via `pip`, then locates the installed samples directory to make the utility scripts executable, and finally runs `ogr2vrt.py` to perform the conversion. The code relies on Python’s `osgeo_utils` module and standard shell commands to manage the GDAL sample tools.

AI-generated description

Similar vizzes

Loading thumbnail…

VRT to merge vector data

This example demonstrates how to merge multiple vector datasets into a single unified dataset using the Virtual Dataset (VRT) format, specifically for French IRIS geographic zones. It shows the workflow of combining shapefiles from different territories and projections into one WGS 84 layer without creating an intermediate file, then converting the result to GeoPackage and PostGIS formats. The code uses the `ogrmerge.py` GDAL utility to create the VRT, `ogr2ogr` for format conversions, and psql commands for database setup.

TThomasG77

69% 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

68% match

Loading thumbnail…

VRT recipe for CSV

This example demonstrates how to convert a CSV file containing building bounding box coordinates into a polygon vector layer using a VRT (Virtual Dataset) recipe. The visualization shows a map of building footprints as rectangles derived from the xmin, xmax, ymin, and ymax columns in the CSV. It works by using an OGRVRTDataSource with SQLite dialect in the `demo-sql-csv.vrt` file, where the `ST_PolygonFromText` function constructs spatial geometries from the coordinate data. The technique is presented as a lightweight alternative to processing GeoJSON files, and can be viewed in QGIS or inspected via GDAL command-line tools like `ogrinfo`.

TThomasG77

65% match

Loading thumbnail…

Gist 56db21b335e25668ff86c2e8359bf3d1

This example visualizes tree inventory data from a CSV file, mapping each tree’s location as a point based on longitude and latitude coordinates. The data includes attributes such as species, planting date, height, and diameter. The implementation uses a VRT file to define the CSV’s geometry, with the QGIS XML configuration specifying the point layer’s symbology and coordinate system. The visualization is rendered through QGIS, leveraging the OGR provider to read the spatial data and display it as a map layer with a single symbol renderer.

TThomasG77

61% match

Loading thumbnail…

Download data from WFS using GDAL Python bindings

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.

TThomasG77

61% match

Loading thumbnail…

Reprojecting CSV with ogr2ogr

This example demonstrates how to convert a CSV file containing projected coordinates (EPSG:2263) into unprojected longitude and latitude values (EPSG:4326). The conversion relies on a virtual dataset defined in a VRT file, which specifies the source CSV, its geometry type, and the source coordinate system. The command-line utility `ogr2ogr` from the GDAL library performs the reprojection, outputting a new CSV with the transformed coordinates. The VRT uses the `PointFromColumns` geometry encoding to interpret the X and Y columns as spatial points.

Mmbostock

61% 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

60% match

Loading thumbnail…

Gist 61afda568887fcc603e1f9ef963925ec

This example demonstrates a QGIS Processing model that exports spatial data to a PostgreSQL/PostGIS database. The model chains two algorithms: `qgis:postgisexecuteandloadsql` to run SQL and load the result, followed by `qgis:importintopostgis` to write the output to a table named "out2" in the public schema. It uses the PostgreSQL connection "EPSG demo," sets UTF-8 encoding, and creates a geometry column. The model shows a reusable workflow for automating database ingestion directly from the QGIS Processing framework.

TThomasG77

60% match