[prev in list] [next in list] [prev in thread] [next in thread] 

List:       postgis-users
Subject:    Re: [postgis-users] PostGIS - Simple Raster Point Query
From:       Tim-Hinnerk Heuer <th.heuer () gmail ! com>
Date:       2012-12-18 23:42:29
Message-ID: CAApG7tZdEzDqBi2a7EHYBNQeqC9QAudLYL-H8X5rOjP+E=W0tg () mail ! gmail ! com
[Download RAW message or body]

[Attachment #2 (multipart/alternative)]


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@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)
>

[Attachment #5 (text/html)]

<div class="gmail_extra">Thanks Pierre! That works like a charm (with some minor \
changes because of projections and different point values). It is noticeably faster \
now too.</div><div class="gmail_extra"><br></div><div class="gmail_extra"> Thanks, \
again. :-)</div><div class="gmail_extra"><br></div><div \
class="gmail_extra">Cheers,</div><div class="gmail_extra">Tim<br \
clear="all"><div>Tim-Hinnerk Heuer<br><br>Twitter: @geekdenz<br><div>Blog: <a \
href="http://www.thheuer.com/" target="_blank">http://www.thheuer.com</a></div> \
</div><br> <br><br><div class="gmail_quote">On 19 December 2012 10:56, Pierre Racine \
<span dir="ltr">&lt;<a href="mailto:Pierre.Racine@sbf.ulaval.ca" \
target="_blank">Pierre.Racine@sbf.ulaval.ca</a>&gt;</span> wrote:<br><blockquote \
class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc \
solid;padding-left:1ex"> <div id=":1aa">SELECT p, geom, ST_Value(rast,geom) elev<br>
FROM (SELECT p, ST_line_interpolate_point(<br>
                           ST_GeometryFromText(&#39;LINESTRING(-69 46.5, -68 46, -67 \
                46.5, -67 45)&#39;, 4326),<br>
                                                         CAST(p AS numeric)/30) \
geom<br>  FROM generate_series(1,30) p<br>
                 ) foo,<br>
              srtm_tiled_100x100<br>
WHERE ST_Intersects(rast, geom)</div></blockquote></div><br></div>



_______________________________________________
postgis-users mailing list
postgis-users@lists.osgeo.org
http://lists.osgeo.org/cgi-bin/mailman/listinfo/postgis-users


[prev in list] [next in list] [prev in thread] [next in thread]