[postgis-users] Query using ST_transform fails

[Torsten Mohr]

Hello,

thanks for your help.

Regarding your questions:

No, i don't use www.pgrpms.org, my system is openSuse 11.2, i installed their 
postgresql, which is version 8.4.4.

Yes, i installed Proj4 from source.  I have to say i'm not sure any more if i
installed Proj-datumgrid in the /nad subdirectory of Proj 4 before
configuring, but in /usr/local/share/proj i got:

-rw-r--r-- 1 root root    694 10. Jan 2010  GL27
-rw-r--r-- 1 root root   6385 10. Jan 2010  nad.lst
-rw-r--r-- 1 root root  19501 10. Jan 2010  nad27
-rw-r--r-- 1 root root  16559 10. Jan 2010  nad83
-rw-r--r-- 1 root root   7043 10. Jan 2010  world
-rw-r--r-- 1 root root    261 10. Jan 2010  proj_def.dat
-rw-r--r-- 1 root root 551012 10. Jan 2010  epsg
-rw-r--r-- 1 root root 453436 10. Jan 2010  esri
-rw-r--r-- 1 root root  76843 10. Jan 2010  esri.extra
-rw-r--r-- 1 root root   3702 10. Jan 2010  other.extra
-rw-r--r-- 1 root root  77820 10. Jan 2010  IGNF


Would it help to re-install Proj 4 ?

Would you expect problems doing so as i installed PostGis after Proj 4 ?


Best regards,
Torsten.


Am Dienstag, 2. November 2010 22:34:05 schrieb Mike Toews:
> It appears the grid shift files are missing. These are not necessarily
> installed by default.
> 
> Are you using http://www.pgrpms.org/ ? Is proj-nad installed? Did you
> install proj.4 from source? If so, did you get the ZIP files and put
> them in the right place before configuring?
> 
> -Mike
> 
> On 2 November 2010 14:09, Torsten Mohr <tmohr at s.netic.de> wrote:
> > Hello,
> >
> > (i accidentially posted this on the openstreetmap mailing list, but it
> > belongs here:)
> >
> >
> > I once got a hint on this mailing list to use a query like this to get
> > the lat/lon of the world capitals:
> >
> > A)
> > select st_X(wayLL), st_Y(wayLL), name from (select
> > ST_AsText(ST_Transform(way,4326)) as wayLL, name from planet_osm_point
> > where capital='yes') as foo limit 5;
> >
> > B)
> > Based on that hint i used this query:
> > select st_X(st_transform(way,4326)), st_Y(st_transform(way,4326)), name
> > from planet_osm_point where place='city' and capital='yes';
> >
> > That query worked fine and i did not change my system since then (that
> > somehow can't be true).  I now get errors for both queries:
> >
> > FEHLER:  transform: couldn't project point (653103 6.63036e+06 0): failed
> > to load NAD27-83 correction file (-38)
> > TIP:  PostGIS was unable to transform the point because either no grid
> > shift files were found, or the point does not lie within the range for
> > which the grid shift is defined. Refer to the ST_Transform() section of
> > the PostGIS manual for details on how to configure PostGIS to alter this
> > behaviour.
> >
> >
> > Could it be that due to an RPM update of PostgreSQL some scripts need to
> > be reinstalled?  I can still generate maps using mapnik.
> >
> >
> > What do i need to do to make those queries work again?
> >
> >
> > In a second step i'd like to extract the polygons that make up the
> > borders of a country as lat/lon values.  I wonder if this is possible by
> > just applying the SQL query above to the way of a polygon?  Or is there
> > another way to get the data i want?
> >
> >
> > Thanks for any hints,
> > Torsten.
> > _______________________________________________
> > postgis-users mailing list
> > postgis-users at postgis.refractions.net
> > http://postgis.refractions.net/mailman/listinfo/postgis-users
> 
> _______________________________________________
> postgis-users mailing list
> postgis-users at postgis.refractions.net
> http://postgis.refractions.net/mailman/listinfo/postgis-users
>