Skip to content

Instantly share code, notes, and snippets.

@andrewxhill

andrewxhill/gist:1447684

Created December 8, 2011 17:18
Show Gist options
  • Star 0 You must be signed in to star a gist
  • Fork 0 You must be signed in to fork a gist
  • Save andrewxhill/1447684 to your computer and use it in GitHub Desktop.
Save andrewxhill/1447684 to your computer and use it in GitHub Desktop.
Download ZIP
Raw
gistfile1.sql
SELECT rid,
encode(
ST_AsPng(
ST_Resample(
ST_ClipEx(
rast,
ST_World2RasterCoordX(
rast,
foo.upper_left),
ST_World2RasterCoordY(
rast,
foo.upper_left),
(ST_World2RasterCoordX(rast,
foo.upper_right)-ST_World2RasterCoordX(rast, foo.upper_left)),
(ST_World2RasterCoordY(rast, foo.upper_left)-ST_World2RasterCoordX(rast, foo.lower_left))),
the_geom_grid)),
'base64') as png
FROM
marso
CROSS JOIN
(SELECT
GeometryFromText('Multipolygon(((-180 0, 0 0, 0 -85.05112877980659, -180 -85.05112877980659, -180 0)))', 4326) as the_geom,
ST_AsRaster(GeometryFromText('Multipolygon(((-180 0, 0 0, 0 -85.05112877980659, -180 -85.05112877980659, -180 0)))', 4326), 256, 256) as the_geom_grid,
GEOMETRYFROMTEXT('POINT(-180 0)',4326) as upper_left,
GEOMETRYFROMTEXT('POINT(-180 -85.05112877980659)',4326) as lower_left,
GEOMETRYFROMTEXT('POINT(0 0)',4326) as upper_right
) foo
WHERE
St_Intersects(rast, foo.the_geom)
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment