That seemed to be the problem. I figured the versions wouldn't make the difference, but I guess if 2163 got shifted that would make sense. Thanks for the input, hopefully this helps anyone else that might have stored old 2163 data and is trying to use the new version.<br>
<br><div class="gmail_quote">On Wed, Oct 20, 2010 at 3:07 PM, Paragon Corporation <span dir="ltr"><<a href="mailto:lr@pcorp.us">lr@pcorp.us</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin: 0pt 0pt 0pt 0.8ex; border-left: 1px solid rgb(204, 204, 204); padding-left: 1ex;">
<div>
<div dir="ltr" align="left"><span><font color="#0000ff" face="Arial" size="2">It is most likely the version of proj4. I know when I
went from proj 4.5 to 4.6, (or was it 4.6 to 4.7). 2163 got shifted.
Frank explained is because it was missing information and in prior proj I think
it guessed and in newer it was decided not to make up stuff.</font></span></div>
<div dir="ltr" align="left"><span><font color="#0000ff" face="Arial" size="2"></font></span> </div>
<div dir="ltr" align="left"><span><font color="#0000ff" face="Arial" size="2">That was a couple years ago. but it is the same national
atlas projection I am fond of.</font></span></div>
<div dir="ltr" align="left"><span><font color="#0000ff" face="Arial" size="2"></font></span> </div>
<div dir="ltr" align="left"><span><font color="#0000ff" face="Arial" size="2">Do a </font></span></div>
<div dir="ltr" align="left"><span><font color="#0000ff" face="Arial" size="2"></font></span> </div>
<div dir="ltr" align="left"><span><font color="#0000ff" face="Arial" size="2">SELECT postgis_full_version();</font></span></div>
<div dir="ltr" align="left"><span><font color="#0000ff" face="Arial" size="2"></font></span> </div>
<div dir="ltr" align="left"><span><font color="#0000ff" face="Arial" size="2">on both -- I bet they are running different versions of
proj.</font></span></div>
<div dir="ltr" align="left"><span><font color="#0000ff" face="Arial" size="2"></font></span> </div>
<div dir="ltr" align="left"><span><font color="#0000ff" face="Arial" size="2">Details here.</font></span></div>
<div dir="ltr" align="left"><span><font color="#0000ff" face="Arial" size="2"></font></span> </div>
<div dir="ltr" align="left"><span><font color="#0000ff" face="Arial" size="2"><a href="http://postgis.refractions.net/pipermail/postgis-users/2008-January/018400.html" target="_blank">http://postgis.refractions.net/pipermail/postgis-users/2008-January/018400.html</a></font></span></div>
<div dir="ltr" align="left"><span><font color="#0000ff" face="Arial" size="2"></font></span> </div>
<div dir="ltr" align="left"><span><font color="#0000ff" face="Arial" size="2">It is a wonder my limited memory actually remembers these
things.</font></span></div>
<div dir="ltr" align="left"><span><font color="#0000ff" face="Arial" size="2"></font></span> </div>
<div dir="ltr" align="left"><span><font color="#0000ff" face="Arial" size="2">Thanks,</font></span></div>
<div dir="ltr" align="left"><span><font color="#0000ff" face="Arial" size="2">Regina</font></span></div>
<div dir="ltr" align="left"><span><font color="#0000ff" face="Arial" size="2"></font></span> </div><font color="#0000ff" face="Arial" size="2"></font><br>
<div dir="ltr" align="left" lang="en-us">
<hr>
<font face="Tahoma" size="2"><b>From:</b>
<a href="mailto:postgis-users-bounces@postgis.refractions.net" target="_blank">postgis-users-bounces@postgis.refractions.net</a>
[mailto:<a href="mailto:postgis-users-bounces@postgis.refractions.net" target="_blank">postgis-users-bounces@postgis.refractions.net</a>] <b>On Behalf Of </b>Josh
Abts<br><b>Sent:</b> Wednesday, October 20, 2010 3:28 PM<br><b>To:</b> PostGIS
Users Discussion<br><b>Subject:</b> Re: [postgis-users] ST_Transform Incorrect
Results?<br></font><br></div><div><div></div><div class="h5">
<div></div>Sure thing!<br><br>SRID 2163<br>Server 1: "+proj=laea +lat_0=45
+lon_0=-100 +x_0=0 +y_0=0 +a=6370997 +b=6370997 +units=m +no_defs "<br>Server 2:
"+proj=laea +lat_0=45 +lon_0=-100 +x_0=0 +y_0=0 +a=6370997 +b=6370997 +units=m
+no_defs "<br><br>SRID 4326<br>Server 1: "+proj=longlat +ellps=WGS84
+datum=WGS84 +no_defs "<br>Server 2: "+proj=longlat +ellps=WGS84 +datum=WGS84
+no_defs "<br><br>Unless I am missing something, they both seem to be exactly
the same...<br><br>I found that earlier as well (forgot to include that as
well). But that really makes me think it is something with the version of
proj4 (or some data file it might be using?), but maybe someone has other
ideas?<br><br>
<div class="gmail_quote">On Wed, Oct 20, 2010 at 2:24 PM, Paul Ramsey <span dir="ltr"><<a href="mailto:pramsey@opengeo.org" target="_blank">pramsey@opengeo.org</a>></span> wrote:<br>
<blockquote style="border-left: 1px solid rgb(204, 204, 204); margin: 0pt 0pt 0pt 0.8ex; padding-left: 1ex;" class="gmail_quote">Could you show the proj4text strings for your coordinate
systems? The<br>srtext strings are actually just there for
decoration.<br><br>P<br>
<div>
<div></div>
<div><br>On Wed, Oct 20, 2010 at 12:16 PM, Josh Abts <<a href="mailto:josh.abts@vistracks.com" target="_blank">josh.abts@vistracks.com</a>>
wrote:<br>> As a followup, the following query reproduces the problem
exactly:<br>><br>> select
astext(st_transform(st_setsrid(st_makepoint(977119.259431679,<br>>
-252548.899746532), 2163) ,4326))<br>><br>> On Server 1 the result is:
POINT(-88.11068 42.2810949683921)<br>><br>> On Server 2 the result is:
POINT(-88.11068 42.0896736466198)<br>><br>> It was also mentioned on
#postgis to include the differing srtext entries I<br>> had found on the
two servers....<br>><br>> For SRID 4326:<br>><br>> Server 1 -
"GEOGCS["WGS 84",DATUM["WGS_1984",SPHEROID["WGS<br>>
84",6378137,298.257223563,AUTHORITY["EPSG","7030"]],TOWGS84[0,0,0,0,0,0,0],AUTHORITY["EPSG","6326"]],PRIMEM["Greenwich",0,AUTHORITY["EPSG","8901"]],UNIT["degree",0.01745329251994328,AUTHORITY["EPSG","9122"]],AUTHORITY["EPSG","4326"]]"<br>
><br>>
Server 2 - "GEOGCS["WGS 84",DATUM["WGS_1984",SPHEROID["WGS<br>>
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"]]"<br>
><br>>
For SRID 2163:<br>><br>> Server 1 - "PROJCS["US National Atlas Equal
Area",GEOGCS["Unspecified datum<br>> based upon the Clarke 1866
Authalic<br>>
Sphere",DATUM["Not_specified_based_on_Clarke_1866_Authalic_Sphere",SPHEROID["Clarke<br>>
1866 Authalic<br>>
Sphere",6370997,0,AUTHORITY["EPSG","7052"]],AUTHORITY["EPSG","6052"]],PRIMEM["Greenwich",0,AUTHORITY["EPSG","8901"]],UNIT["degree",0.01745329251994328,AUTHORITY["EPSG","9122"]],AUTHORITY["EPSG","4052"]],PROJECTION["Lambert_Azimuthal_Equal_Area"],PARAMETER["latitude_of_center",45],PARAMETER["longitude_of_center",-100],PARAMETER["false_easting",0],PARAMETER["false_northing",0],UNIT["metre",1,AUTHORITY["EPSG","9001"]],AUTHORITY["EPSG","2163"]]"<br>
><br>>
Server 2 - "PROJCS["unnamed",GEOGCS["unnamed<br>>
ellipse",DATUM["unknown",SPHEROID["unnamed",6370997,0]],PRIMEM["Greenwich",0],UNIT["degree",0.0174532925199433]],PROJECTION["Lambert_Azimuthal_Equal_Area"],PARAMETER["latitude_of_center",45],PARAMETER["longitude_of_center",-100],PARAMETER["false_easting",0],PARAMETER["false_northing",0],UNIT["Meter",1],AUTHORITY["EPSG","2163"]]"<br>
><br>><br>>
On Wed, Oct 20, 2010 at 11:14 AM, Josh Abts <<a href="mailto:josh.abts@vistracks.com" target="_blank">josh.abts@vistracks.com</a>>
wrote:<br>>><br>>> Hello,<br>>><br>>> I have been
having a problem where I am getting mixed results from two<br>>> servers
using the same dataset. The command I am running is as
follows:<br>>><br>>> select
astext(st_transform(start_point,<br>>>
4326)),astext(st_transform(end_point, 4326)) from mytable where
gid=123<br>>><br>>> The two columns are originally using SRID 2163
so I am transforming them<br>>> into 4326. But the results on one
server (correct) are different than the<br>>> results on a second server
(incorrect). They are off by .2 latitude<br>>>
roughly.<br>>><br>>> For example one result set is as
follows:<br>>><br>>> Server 1 (correct): "POINT(-88.11068
42.2810949683921)";"POINT(-88.10731<br>>> 42.2844149305216)"<br>>>
Server 2 (incorrect): "POINT(-88.11068
42.0896736466198)";"POINT(-88.10731<br>>>
42.0929914612221)"<br>>><br>>> As you can see the longitudes are
correctly transformed, but the latitude<br>>> is
not.<br>>><br>>> I have also included the postgis_full_version
output in case a version<br>>> issue might be the
culprit.<br>>><br>>> Server 1: "POSTGIS="1.3.2"
GEOS="3.0.0rc4-CAPI-1.3.3" PROJ="Rel. 4.5.0, 22<br>>> Oct 2006"
USE_STATS"<br>>> Server 2: "POSTGIS="1.5.2" GEOS="3.2.2-CAPI-1.6.2"
PROJ="Rel. 4.7.1, 23<br>>> September 2009" LIBXML="2.7.5" USE_STATS
(procs from 1.5 r5385 need<br>>> upgrade)"<br>>><br>>> So
you can see I am using a newer version of PostGIS as well as Proj4
but<br>>> I can't seem to figure out why this would result in an
incorrect<br>>> St_Transform result.<br>>><br>>> I also
checked further into the spatial_ref_sys table on both servers and<br>>>
Server 1 has slightly more detailed srtext but they have the same
proj4text<br>>> lines. I tried copying the srtext from Server 1 to
Server 2 and there was<br>>> no change.<br>>><br>>> The only
other difference would be the system platforms themselves which<br>>>
are below.<br>>><br>>> Server 1: "PostgreSQL 8.3.4, compiled by
Visual C++ build 1400" (this is a<br>>> windows 2003 server)<br>>>
Server 2: "PostgreSQL 8.4.5 on x86_64-pc-linux-gnu, compiled by
GCC<br>>> gcc-4.4.real (Ubuntu 4.4.1-4ubuntu9) 4.4.1,
64-bit"<br>>><br>>> Anyone have any ideas of what I could look
into to see why I might be<br>>> getting this slight variation in
results? I can't find anything that would<br>>> seem to be the
problem or cause of a varying transformation...<br>>><br>>>
Thanks!<br>>> Josh<br>><br>><br>><br>> --<br>> Joshua
Abts<br>> IT Specialist<br>> VisTracks, Inc.<br>> <a href="mailto:josh.abts@vistracks.com" target="_blank">josh.abts@vistracks.com</a><br>>
(630) 596-5426<br>><br></div></div>>
_______________________________________________<br>> postgis-users mailing
list<br>> <a href="mailto:postgis-users@postgis.refractions.net" target="_blank">postgis-users@postgis.refractions.net</a><br>>
<a href="http://postgis.refractions.net/mailman/listinfo/postgis-users" target="_blank">http://postgis.refractions.net/mailman/listinfo/postgis-users</a><br>><br>><br>_______________________________________________<br>
postgis-users
mailing list<br><a href="mailto:postgis-users@postgis.refractions.net" target="_blank">postgis-users@postgis.refractions.net</a><br><a href="http://postgis.refractions.net/mailman/listinfo/postgis-users" target="_blank">http://postgis.refractions.net/mailman/listinfo/postgis-users</a><br>
</blockquote></div><br><br clear="all"><br>-- <br>Joshua Abts<br>IT Specialist<br>VisTracks, Inc.<br><a href="mailto:josh.abts@vistracks.com" target="_blank">josh.abts@vistracks.com</a><br>(630)
596-5426<br></div></div></div>
<br>_______________________________________________<br>
postgis-users mailing list<br>
<a href="mailto:postgis-users@postgis.refractions.net">postgis-users@postgis.refractions.net</a><br>
<a href="http://postgis.refractions.net/mailman/listinfo/postgis-users" target="_blank">http://postgis.refractions.net/mailman/listinfo/postgis-users</a><br>
<br></blockquote></div><br><br clear="all"><br>-- <br>Joshua Abts<br>IT Specialist<br>VisTracks, Inc.<br><a href="mailto:josh.abts@vistracks.com">josh.abts@vistracks.com</a><br>(630) 596-5426<br>