Load GTFS as GPKG with GDAL 3.7

100%

Load GTFS as GPKG with GDAL 3.7

✓ Published0🌍 Public

Last edited Jun 28, 2023

Created on Jun 28, 2023

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.

AI-generated description

Similar vizzes

Loading thumbnail…

Gist c32a8597373bd00c10029d4f21f2400e

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.

TThomasG77

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

66% 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…

Multiple approaches to solve question at https://lists.osgeo.org/pipermail/gdal-dev/2022-June/055949.html

This example demonstrates three equivalent methods for converting a CSV file with latitude and longitude columns into a GeoPackage with point geometry using GDAL’s `gdal.VectorTranslate`, GeoPandas’ `read_file` and `to_file`, and Fiona’s manual schema definition and feature writing. The CSV contains a single data point (61, -150), and each approach specifies the coordinate column names via `X_POSSIBLE_NAMES` and `Y_POSSIBLE_NAMES` options. The code shows how to set the output CRS to EPSG:4326 and assign a layer name (`layername`) in the resulting GeoPackage.

TThomasG77

63% match

Loading thumbnail…

mapshaper-info-recreated-via-gdal

This example recreates the command-line tool mapshaper-info using the gdal-async library to inspect a GeoJSON TopoJSON file containing German administrative regions. It downloads the dataset from GitHub, then prints layer metadata including geometry type, record count, bounding box, and coordinate reference system, followed by a formatted table of attribute names and their first values. The script uses the got library for streaming downloads and the pipeline function from stream/promises to handle the file writing asynchronously.

TThomasG77

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

62% match

Loading thumbnail…

Play with gdal.vectorTranslate in Node.js

This example demonstrates a complete geospatial data processing pipeline using the Node.js `gdal-async` library. It downloads a Natural Earth countries dataset, extracts the ZIP archive, and converts the shapefile to GeoPackage format with `gdal.vectorTranslate`, then transforms the GeoPackage into GeoJSON. The script uses the `got` library for HTTP streaming with progress reporting and `adm-zip` for extraction. This showcases how GDAL's vector translation capabilities can be invoked programmatically in a Node.js environment for format conversion between common geospatial file types.

TThomasG77

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