[postgis-devel] [PostGIS] #318: 900913->4326 Transform

[PostGIS]

#318: 900913->4326 Transform
----------------------+-----------------------------------------------------
  Reporter:  pramsey  |       Owner:  mcayland     
      Type:  defect   |      Status:  new          
  Priority:  medium   |   Milestone:  PostGIS 1.4.1
 Component:  postgis  |     Version:  trunk        
Resolution:           |    Keywords:               
----------------------+-----------------------------------------------------
Comment (by arthurnederlof):

 Hallo,
 a more trivial example (osm form open street map, worldcountries has
 geometries in 4326 srid format:
 this works:
 select world.gid, osm.admin_level from planet_osm_polygon as
 osm,worldcountries as world  where (osm.admin_level='4')  AND
 (ST_Contains(ST_Transform(world.the_geom,900913),ST_Centroid(osm.way)))
 limit 10

 this one does not work:
 select world.gid, osm.admin_level from planet_osm_polygon as
 osm,worldcountries as world  where (osm.admin_level='4')  AND
 (ST_Contains(world.the_geom,ST_Centroid(ST_Transform(osm.way,4326))))
 limit 10

 ERROR:  transform: couldn't project point (5.1755e+06 6.1799e+06 0):
 failed to load NAD27-83 correction file (-38)
 HINT:  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.

-- 
Ticket URL: <http://trac.osgeo.org/postgis/ticket/318#comment:11>
PostGIS <http://trac.osgeo.org/postgis/>
The PostGIS Trac is used for bug, enhancement & task tracking, a user and developer wiki, and a view into the subversion code repository of PostGIS project.