[postgis-users] PostGIS - Simple Raster Point Query

[Tim-Hinnerk Heuer]

Thanks Pierre! That works like a charm (with some minor changes because of
projections and different point values). It is noticeably faster now too.

Thanks, again. :-)

Cheers,
Tim
Tim-Hinnerk Heuer

Twitter: @geekdenz
Blog: http://www.thheuer.com



On 19 December 2012 10:56, Pierre Racine <Pierre.Racine at sbf.ulaval.ca>wrote:

> SELECT p, geom, ST_Value(rast,geom) elev
> FROM (SELECT p, ST_line_interpolate_point(
>                   ST_GeometryFromText('LINESTRING(-69 46.5, -68 46, -67
> 46.5, -67 45)', 4326),
>                                       CAST(p AS numeric)/30) geom
>             FROM generate_series(1,30) p
>            ) foo,
>          srtm_tiled_100x100
> WHERE ST_Intersects(rast, geom)
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.osgeo.org/pipermail/postgis-users/attachments/20121219/b5c95e9c/attachment.html>