Hi,<br>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. <br>I am now trying to compute an intersection between the SRTM rasters and the OSM road geometries... with no luck! :(<br>
I am following this tutorial --> <a href="http://trac.osgeo.org/postgis/wiki/WKTRasterTutorial01">http://trac.osgeo.org/postgis/wiki/WKTRasterTutorial01</a><br>but I think it might be slightly outdated.<br><br> SELECT o.osm_id, <br>
        ST_Intersection(s.rast, o.way)<br> FROM planet_osm_line o, <br>      srtm_tiled s<br> WHERE ST_Intersects(s.rast, o.way);<br><br>way type -->   <b>way geometry(LineString,900913)</b><br><br>CREATE TABLE public.srtm_tiled<br>
(<br>  rid integer NOT NULL DEFAULT nextval('srtm_tiled_rid_seq'::regclass),<br> <b> rast raster,</b><br>  CONSTRAINT srtm_tiled_pkey PRIMARY KEY (rid)<br>)<br>WITH (<br>  OIDS=FALSE<br>);<br><br>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? <br>
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.<br>Can someone help?<br><br><br>Cheers,<br>Ed<br><br>