very nice explanation Mike<br><br><div class="gmail_quote">On Thu, Aug 12, 2010 at 11:39 AM, Mike Toews <span dir="ltr"><<a href="mailto:mwtoews@gmail.com">mwtoews@gmail.com</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex;">
Hi,<br>
<br>
Your coordinates may be flipped. Was it 59N 18E? If so, use x,y<br>
notation: 'POINT(18 59)', which results in 'POINT(18.0000000000006<br>
58.9999999999905)', which is close enough.<br>
<br>
Also keep in mind that you are outside the projection bounds:<br>
<a href="http://spatialreference.org/ref/epsg/3021/" target="_blank">http://spatialreference.org/ref/epsg/3021/</a> (just a bit too far east).<br>
Whenever you are outside the projection bounds, the likelihood of<br>
storage precision errors increase. To understand why this is, you can<br>
think of taking the tangent of a two angles that are nearly a<br>
right-angle (89.9991 and 89.9992) which have very different results<br>
due to nature of the geometry.<br>
<font color="#888888"><br>
-Mike<br>
</font><div><div></div><div class="h5"><br>
On 11 August 2010 13:10, Erik Rehn <<a href="mailto:erik@slagkryssaren.com">erik@slagkryssaren.com</a>> wrote:<br>
> Hello Postgis Users!<br>
><br>
> This is my first post on this list so I will start by asking<br>
> a simple (and probably stupid) question. :)<br>
><br>
> While using ST_AsKml() to produce an overlay for Google Earth I<br>
> noticed that all my geometries where shifted slightly south-east.<br>
> I figured this had something to do with the transformation between<br>
> the projection that my geometries are stored in (SRID 3021) and WGS84 (4326)<br>
> that is outputted by ST_AsKml()<br>
><br>
> Just to test I ran this:<br>
><br>
> SELECT ST_AsText(<br>
>    ST_Transform(<br>
>        ST_Transform(<br>
>            ST_GeomFromText('POINT(59 18)',4326),<br>
>        3021),<br>
>    4326));<br>
><br>
> I input a point in WGS84 (59,18), transforms it to 3021 and then back to<br>
> WGS84. The result I get is:<br>
> POINT(58.8672757036296 18.0394763349359)<br>
><br>
> Can anyone explain this? Am I missing something regarding ST_Transform()?<br>
><br>
> Im running Postgis 1.5 on Windows.<br>
><br>
> Thank you for any help!<br>
> /Erik<br>
><br>
> --<br>
> Erik Rehn<br>
> Slagkryssaren<br>
> <a href="mailto:erik@slagkryssaren.com">erik@slagkryssaren.com</a><br>
> <a href="http://www.slagkryssaren.com" target="_blank">www.slagkryssaren.com</a><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>
_______________________________________________<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>
</div></div></blockquote></div><br>