Skip to main content
100%

A temporary fix to raster2pgsql.py to work with current PostGIS HEAD

✓ Published0🌍 Public
Aandrewxhill
Last edited Dec 17, 2011
Created on Dec 17, 2011

A Python script loads GDAL raster datasets and converts them into hex-encoded WKB streams for insertion into a PostGIS raster table. It uses the `osgeo.gdal` and `osgeo.osr` bindings to read geospatial data, with `optparse` handling command-line flags for SRID assignment, band extraction, block-size tiling, and overview table creation. The output is a SQL-ready binary representation suitable for PostGIS’s current raster format.

AI-generated description

Similar vizzes

Loading thumbnail…

Rasters to polygons for CartoDB

A PostGIS workflow converts nightly satellite raster data into vector polygons for display in CartoDB. The SQL first uses `ST_MapAlgebraExpr` to reclassify raster band values into discrete categories, then `ST_DumpAsPolygons` extracts each distinct value region as a geometry, wrapping results with `ST_MakeValid` and `ST_Multi`. The resulting polygons, along with their original raster values, are inserted into a new `light_polygons` table with a `band` column, enabling efficient cartographic rendering of classified light intensity zones.

Aandrewxhill
62% match
Loading thumbnail…

Gist 1353427

This example shows how to sum raster pixel values from a PostGIS database that fall within a specified multi-polygon region, using SQL queries to perform spatial intersection tests. The code defines a multi-polygon geometry and iterates over a 10×10 grid of pixel offsets, translating each pixel’s envelope and checking for intersection with the polygon. It uses PostGIS functions like `ST_Value`, `ST_Intersects`, `ST_Translate`, and `ST_MakeEnvelope` to calculate the total sum and count of intersecting raster cells, demonstrating a spatial aggregation workflow.

Aandrewxhill
62% match
Loading thumbnail…

Gist 1353003

This example demonstrates a spatial SQL query that computes raster statistics for a custom polygon region. It uses PostGIS functions to sum pixel values, count pixels, and calculate the average value across a 10x10 grid of sub-pixels within each raster tile. The query filters raster tiles intersecting the target polygon and then intersects translated grid cells against the same polygon. The code relies on `ST_Value`, `ST_Intersects`, `ST_Translate`, `ST_MakeEnvelope`, and `generate_series` to perform the analysis, extracting data from the `carbon_query_test_s` table.

Aandrewxhill
62% match
Loading thumbnail…

Gist 1447684

This visualization shows a clipped portion of a Mars elevation raster dataset, rendering it as a base64-encoded PNG image for a defined geographic region. The SQL code uses PostGIS functions including `ST_AsPng`, `ST_Resample`, and `ST_ClipEx` to extract and resample the raster, while `ST_World2RasterCoordX` and `ST_World2RasterCoordY` convert real-world coordinates to raster grid positions. A multipolygon defines the area of interest, which is also converted to a raster grid via `ST_AsRaster` for alignment. The query filters rows with `St_Intersects` to return only relevant raster tiles, demonstrating spatial clipping and coordinate transformation for remote-sensing data delivery.

Aandrewxhill
61% match
Loading thumbnail…

GDAL/OGR

This example demonstrates how to convert a File Geodatabase (file.gdb) into a PostgreSQL database using the GDAL/OGR command-line tool. The code executes `ogr2ogr` with a PostgreSQL connection string specifying the database name and user, directly transferring spatial and attribute data between formats. It highlights the library’s role in migrating geospatial datasets without additional scripting. The command handles the conversion in a single step, preserving the original file’s structure. No visualization or rendering occurs; the focus is purely on the data interchange process via GDAL’s OGR vector driver.

Aapburnes
59% match
Loading thumbnail…

Gist 6114786

This example implements a private torque-style time-series visualization for geospatial data using PostgreSQL and PostGIS. It defines custom SQL functions that aggregate point features into web-Mercator grid cells, grouping counts by time steps and returning arrays of values and dates for rendering. The code uses `ST_SnapToGrid`, `CDB_XYZ_Resolution`, and `CDB_XYZ_Extent` to tile data at a given zoom level, with two variants: a standard response for instantaneous counts and a cumulative response that includes expiration dates for each cell. A metadata function returns the spatial bounds and time range from the source table, supporting both timestamped and numeric columns.

Aandrewxhill
58% match
Loading thumbnail…

Python script to covert GeoTiff .CLR file to a Carto style for CartoDB

This script converts a GeoTiff .CLR color palette file into a CartoCSS style for CartoDB. It reads each line of the .CLR file, extracts a raster cell value and its corresponding RGB color, and generates a rule that assigns that polygon-fill color to matching values in a specified column. The resulting CartoDB.carto.css file uses a default orange fill with white outlines for unmatched values. The code relies on Python’s sys.argv for command-line arguments and standard file I/O, with no external visualization libraries.

Aandrewxhill
57% match
Loading thumbnail…

command-line python interface for manipulating data on CartoDB

This example demonstrates a command-line Python tool for managing and manipulating data stored on CartoDB, allowing users to upload files, overwrite or append rows, add columns, drop tables, and export datasets. It first shows the workflow of interacting with CartoDB’s SQL and Import APIs, then explains the underlying implementation. The script uses the `requests` library for HTTP calls, `urllib2` for polling import status, and `argparse` for command-line options, translating between CartoDB’s data types and PostgreSQL types during column alterations and insertions.

Aandrewxhill
57% match