I did a "yum install proj-nad" on my Fedora c 15 box, and the projection error is gone. Now I have a new one! :(<br><br><br>SELECT o.osm_id,<br> ST_Intersection(ST_Transform(o.way,4326), s.rast)<br> FROM planet_osm_line o,<br>
srtm_tiled_4326 s<br> WHERE ST_Intersects(ST_Transform(o.way,4326), s.rast)<br> and o.highway = 'primary';<br><br>I keep seeing the following messages below! <br><br>NOTICE: Ring Self-intersection at or near point 8.50333 55.5392<br>
CONTEXT: PL/pgSQL function "st_intersection" line 10 at RETURN QUERY<br>NOTICE: Ring Self-intersection at or near point 8.87833 55.5408<br>CONTEXT: PL/pgSQL function "st_intersection" line 10 at RETURN QUERY<br>
NOTICE: Ring Self-intersection at or near point 8.63 55.5192<br>CONTEXT: PL/pgSQL function "st_intersection" line 10 at RETURN QUERY<br>NOTICE: Ring Self-intersection at or near point 8.685 55.5083<br><br>Can someone please explain what these messages mean? Am I doing something wrong?<br>
<br>Cheers,<br>Ed<br><br><br> <br><br><div class="gmail_quote">On Tue, Feb 28, 2012 at 9:21 AM, Ed Linde <span dir="ltr"><<a href="mailto:edolinde@gmail.com">edolinde@gmail.com</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
I have loaded my SRTM tif files as srid = 4326. And use a ST_Transform on the road geometry (As suggested).<br>Running Query: <br><br>SELECT o.osm_id,<br> ST_Intersection(ST_Transform(o.way,4326), s.rast)<br> FROM planet_osm_line o,<br>
srtm_tiled_4326 s<br> WHERE ST_Intersects(ST_Transform(o.way,4326), s.rast);<br><br>I get this error:<br><br><i>NOTICE: PostGIS was unable to transform the point because either no grid shift files were found, or the point does <br>
not lie within the range for which the grid shift is defined. Refer to the ST_Transform() section of the PostGIS manual <br>for details on how to configure PostGIS to alter this behaviour.<br><br>ERROR: transform: couldn't project point (977705 5.98119e+06 0): failed to load NAD27-83 correction file (-38)<br>
********** Error **********<br>ERROR: transform: couldn't project point (977705 5.98119e+06 0): failed to load NAD27-83 correction file (-38)<br>SQL state: XX000<br><br></i>I looked up the postgis manual on ST_Transform and it suggests doing ---> <br>
<pre>UPDATE spatial_ref_sys SET proj4text = '+proj=longlat +ellps=clrk66 +nadgrids=@conus,@alaska,@ntv2_0.gsb,@ntv1_can.dat,null +no_defs' WHERE srid = 4326;</pre>Is this the correct way to add a grid? Or is there some better way? Some other threads on this mailing list suggest copying some files over <br>
and placing them in another folder.. not sure what is the better way. Can someone please suggest?<br><br>I am on --><br>-------------------------------------------------------------------------------------------------------------------------------------------------------<br>
<b> POSTGIS="2.0.0alpha7SVN" GEOS="3.3.0-CAPI-1.7.0" PROJ="Rel. 4.7.1, 23 September 2009" GDAL="GDAL 1.7.3, released 2010/11/10" LIBXML="2.7.8" USE_STATS<br></b>(1 row)<br>
<br>Cheers,<br>Ed<div class="HOEnZb"><div class="h5"><br><br><br><br><br><div class="gmail_quote">On Mon, Feb 27, 2012 at 4:05 PM, Pierre Racine <span dir="ltr"><<a href="mailto:Pierre.Racine@sbf.ulaval.ca" target="_blank">Pierre.Racine@sbf.ulaval.ca</a>></span> wrote:<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div>> ok, just to get this right.. :) , you are suggesting I keep my OSM as is (id: 900913).<br>
> Load in the tif rasters as srid (4326) and then "reproject" the OSM road geometry<br>
> to srid = 4326 right? So it matches the srid of my raster! I ask because you said<br>
> to reproject to 900913 in your mail.<br>
<br>
</div>Correct. We're on the right way.<br>
<br>
You can create a new raster table (after loading them properly with 4326) with your rasters reprojected to 900913 but as I said you might hit our bug:<br>
<br>
<a href="http://trac.osgeo.org/postgis/ticket/1618" target="_blank">http://trac.osgeo.org/postgis/ticket/1618</a><br>
<div><div><br>
Pierre<br>
_______________________________________________<br>
postgis-users mailing list<br>
<a href="mailto:postgis-users@postgis.refractions.net" target="_blank">postgis-users@postgis.refractions.net</a><br>
<a href="http://postgis.refractions.net/mailman/listinfo/postgis-users" target="_blank">http://postgis.refractions.net/mailman/listinfo/postgis-users</a><br>
</div></div></blockquote></div><br>
</div></div></blockquote></div><br>