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
_______________________________________________ postgis-users mailing list postgis-users@postgis.refractions.net http://postgis.refractions.net/mailman/listinfo/postgis-users