[postgis-users] ST_Transform Incorrect Results?
Josh Abts
josh.abts at vistracks.com
Wed Oct 20 09:14:12 PDT 2010
Hello,
I have been having a problem where I am getting mixed results from two
servers using the same dataset. The command I am running is as follows:
select astext(st_transform(start_point,
4326)),astext(st_transform(end_point, 4326)) from mytable where gid=123
The two columns are originally using SRID 2163 so I am transforming them
into 4326. But the results on one server (correct) are different than the
results on a second server (incorrect). They are off by .2 latitude
roughly.
For example one result set is as follows:
Server 1 (correct): "POINT(-88.11068 42.2810949683921)";"POINT(-88.10731
42.2844149305216)"
Server 2 (incorrect): "POINT(-88.11068 42.0896736466198)";"POINT(-88.10731
42.0929914612221)"
As you can see the longitudes are correctly transformed, but the latitude is
not.
I have also included the postgis_full_version output in case a version issue
might be the culprit.
Server 1: "POSTGIS="1.3.2" GEOS="3.0.0rc4-CAPI-1.3.3" PROJ="Rel. 4.5.0, 22
Oct 2006" USE_STATS"
Server 2: "POSTGIS="1.5.2" GEOS="3.2.2-CAPI-1.6.2" PROJ="Rel. 4.7.1, 23
September 2009" LIBXML="2.7.5" USE_STATS (procs from 1.5 r5385 need
upgrade)"
So you can see I am using a newer version of PostGIS as well as Proj4 but I
can't seem to figure out why this would result in an incorrect St_Transform
result.
I also checked further into the spatial_ref_sys table on both servers and
Server 1 has slightly more detailed srtext but they have the same proj4text
lines. I tried copying the srtext from Server 1 to Server 2 and there was
no change.
The only other difference would be the system platforms themselves which are
below.
Server 1: "PostgreSQL 8.3.4, compiled by Visual C++ build 1400" (this is a
windows 2003 server)
Server 2: "PostgreSQL 8.4.5 on x86_64-pc-linux-gnu, compiled by GCC
gcc-4.4.real (Ubuntu 4.4.1-4ubuntu9) 4.4.1, 64-bit"
Anyone have any ideas of what I could look into to see why I might be
getting this slight variation in results? I can't find anything that would
seem to be the problem or cause of a varying transformation...
Thanks!
Josh
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.osgeo.org/pipermail/postgis-users/attachments/20101020/3f27f746/attachment.html>
More information about the postgis-users
mailing list