[postgis-users] ST_Transform problem (bad results)

[Javier Sanchez]

Hello

I'm trying to transform point coordinates from geographic latlon ED50
(EPSG:4230) to WGS84 (EPSG:4326) in Postgis, but the result is incorrect
(input and ouput coordinates are the same, although they must be diferent):

select ST_AsText(ST_Transform(ST_SetSRID(ST_MakePoint(10,20), 4230), 4326));
  st_astext   
--------------
 POINT(10 20)

The same problem with cs2cs from proj4 in a shell

echo 10 20 | cs2cs +init=epsg:4230 +to +init=epsg:4326
10 20

After sources download and compilation, cs2cs gives good results:

echo 10 20 | cs2cs +init=epsg:4230 +to +init=epsg:4326
9d59'57.2"E	19d59'55.54"N 102.459

But postgis is still returning the same point.

I have checked /usr/share/proj and epsg file includes both 4230 and 4326 ok

# ED50
<4230> +proj=longlat +ellps=intl +no_defs  <>
# WGS 84
<4326> +proj=longlat +ellps=WGS84 +datum=WGS84 +no_defs  <>

The same problem occurs in three pcs with Ubuntu 10.10 and
ubuntugis-unstable PPA.

I would apreciate very much any hint to solve this problem.

Javier.
-- 
View this message in context: http://old.nabble.com/ST_Transform-problem-%28bad-results%29-tp31303595p31303595.html
Sent from the PostGIS - User mailing list archive at Nabble.com.