[postgis-users] Raster + Geometry Intersection

[Ed Linde]

Hi,
I am trying to compute elevations for a road network that I have downloaded
and installed into my pg database using osm2pgsql. I also have loaded SRTM
tif files and they are stored as type RASTER in my table.
I am now trying to compute an intersection between the SRTM rasters and the
OSM road geometries... with no luck! :(
I am following this tutorial -->
http://trac.osgeo.org/postgis/wiki/WKTRasterTutorial01
but I think it might be slightly outdated.

 SELECT o.osm_id,
        ST_Intersection(s.rast, o.way)
 FROM planet_osm_line o,
      srtm_tiled s
 WHERE ST_Intersects(s.rast, o.way);

way type -->   *way geometry(LineString,900913)*

CREATE TABLE public.srtm_tiled
(
  rid integer NOT NULL DEFAULT nextval('srtm_tiled_rid_seq'::regclass),
 * rast raster,*
  CONSTRAINT srtm_tiled_pkey PRIMARY KEY (rid)
)
WITH (
  OIDS=FALSE
);

performing the above query returns no rows. So I am not sure what is
wrong.. or if I am doing this query correctly. Is the SRID set differently
in the raster and geometry?
I also loaded the tif files using the tutorial, but I used the C version
not the python loader because I understand its been deprecated.
Can someone help?


Cheers,
Ed
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.osgeo.org/pipermail/postgis-users/attachments/20120227/3982aeef/attachment.html>