[postgis-users] Coordinate changes using transform from 1.3.1 to 1.3.6

Tue May 26 05:23:04 PDT 2009

Hello,

the SRID definition in spatial_ref_sys is exactly the same. Remember that I
am restoring the same database dump in two different DBMS.

The information:

First DBMS:

SELECT version(), postgis_full_version();

"PostgreSQL 8.2.13 on i686-pc-mingw32, compiled by GCC gcc.exe (GCC) 3.4.2
(mingw-special)";"POSTGIS="1.3.1" GEOS="3.0.0rc4-CAPI-1.3.3" PROJ="Rel.
4.5.0, 22 Oct 2006" USE_STATS DBPROC="0.3.0" RELPROC="1.3.1" (needs proc
upgrade)"

spatial_ref_sys

4326;"EPSG";4326;"GEOGCS["WGS 84",DATUM["WGS_1984",SPHEROID["WGS
84",6378137,298.257223563,AUTHORITY["EPSG","7030"]],AUTHORITY["EPSG","6326"]],PRIMEM["Greenwich",0,AUTHORITY["EPSG","8901"]],UNIT["degree",0.01745329251994328,AUTHORITY["EPSG","9122"]],AUTHORITY["EPSG","4326"]]";"+proj=longlat
+ellps=WGS84 +datum=WGS84 +no_defs "

23029;"EPSG";23029;"PROJCS["ED50 / UTM zone
29N",GEOGCS["ED50",DATUM["European_Datum_1950",SPHEROID["International
1924",6378388,297,AUTHORITY["EPSG","7022"]],AUTHORITY["EPSG","6230"]],PRIMEM["Greenwich",0,AUTHORITY["EPSG","8901"]],UNIT["degree",0.01745329251994328,AUTHORITY["EPSG","9122"]],AUTHORITY["EPSG","4230"]],PROJECTION["Transverse_Mercator"],PARAMETER["latitude_of_origin",0],PARAMETER["central_meridian",-9],PARAMETER["scale_factor",0.9996],PARAMETER["false_easting",500000],PARAMETER["false_northing",0],UNIT["metre",1,AUTHORITY["EPSG","9001"]],AUTHORITY["EPSG","23029"]]";"+proj=utm
+zone=29 +ellps=intl +units=m +no_defs "

Second DBMS:

SELECT version(), postgis_full_version();

"PostgreSQL 8.2.13 on i686-pc-mingw32, compiled by GCC gcc.exe (GCC) 3.4.2
(mingw-special)";"POSTGIS="1.3.6" GEOS="3.1.0-CAPI-1.5.0" PROJ="Rel. 4.6.1,
21 August 2008" USE_STATS DBPROC="0.3.0" RELPROC="1.3.6" (needs proc
upgrade)"

spatial_ref_sys

4326;"EPSG";4326;"GEOGCS["WGS 84",DATUM["WGS_1984",SPHEROID["WGS
84",6378137,298.257223563,AUTHORITY["EPSG","7030"]],AUTHORITY["EPSG","6326"]],PRIMEM["Greenwich",0,AUTHORITY["EPSG","8901"]],UNIT["degree",0.01745329251994328,AUTHORITY["EPSG","9122"]],AUTHORITY["EPSG","4326"]]";"+proj=longlat
+ellps=WGS84 +datum=WGS84 +no_defs "

23029;"EPSG";23029;"PROJCS["ED50 / UTM zone
29N",GEOGCS["ED50",DATUM["European_Datum_1950",SPHEROID["International
1924",6378388,297,AUTHORITY["EPSG","7022"]],AUTHORITY["EPSG","6230"]],PRIMEM["Greenwich",0,AUTHORITY["EPSG","8901"]],UNIT["degree",0.01745329251994328,AUTHORITY["EPSG","9122"]],AUTHORITY["EPSG","4230"]],PROJECTION["Transverse_Mercator"],PARAMETER["latitude_of_origin",0],PARAMETER["central_meridian",-9],PARAMETER["scale_factor",0.9996],PARAMETER["false_easting",500000],PARAMETER["false_northing",0],UNIT["metre",1,AUTHORITY["EPSG","9001"]],AUTHORITY["EPSG","23029"]]";"+proj=utm
+zone=29 +ellps=intl +units=m +no_defs "

Thank you very much in advance.

Best regards,
  Miguel

On Tue, May 26, 2009 at 1:39 PM, Mark Cave-Ayland <
mark.cave-ayland at siriusit.co.uk> wrote:

> Miguel R. Luaces wrote:
>
>  Dear all,
>>
>> we have a PostgreSQL 8.1 database with postgis 1.2 installed. We have some
>> data in the database stored using EPSG:23029.
>>
>> We make a backup with pg_dump and we restore the database in PostgreSQL
>> 8.2 with Postgis 1.3.1 without any problem.
>>
>> But, if we restore the backup in PostgreSQL 8.2 with Postgis 1.3.6, the
>> function transform() returns a different result. More precisely:
>>
>> The same query:
>>
>> select astext(the_geom), astext(transform(the_geom)) from monumentos where
>> id = ,
>>
>> In Postgis 1.3.1:
>>
>> "Castelo de San Antón";"POINT(549556.009066526
>> 4801796.3531214)";"POINT(-8.38839299999999 43.3657910268222)"
>>
>> In Postgis 1.3.6:
>>
>> "Castelo de San Antón";"POINT(549556.009066526
>> 4801796.3531214)";"POINT(-8.38839299999999 43.3666118995462)"
>>
>> The second coordinate of the transformation to 4326 changes by 0.004
>> degrees (approx).
>>
>> Any idea?
>>
>> Thank you very much in advance.
>>
>> Best regards,
>>  Miguel
>>
>
>
> Hi Miguel,
>
> What versions of PostgreSQL/PostGIS are you using on both servers, and with
> which versions of the PROJ.4 library?
>
> SELECT version(), postgis_full_version();
>
> Also can you confirm that the spatial_ref_sys entries for 23029 and 4326
> are identical between the two servers?
>
>
> ATB,
>
> Mark.
>
> --
> Mark Cave-Ayland - Senior Technical Architect
> PostgreSQL - PostGIS
> Sirius Corporation plc - control through freedom
> http://www.siriusit.co.uk
> t: +44 870 608 0063
>
> _______________________________________________
> postgis-users mailing list
> postgis-users at postgis.refractions.net
> http://postgis.refractions.net/mailman/listinfo/postgis-users
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.osgeo.org/pipermail/postgis-users/attachments/20090526/d050972e/attachment.html>