<div dir="ltr"><div dir="ltr"><div class="gmail_default" style="font-family:tahoma,sans-serif">Thank you, Jeff. This is helpful. </div><div class="gmail_default" style="font-family:tahoma,sans-serif"><br></div><div class="gmail_default" style="font-family:tahoma,sans-serif">Shira</div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Thu, Jul 23, 2020 at 5:04 AM Jeff McKenna <<a href="mailto:jmckenna@gatewaygeomatics.com">jmckenna@gatewaygeomatics.com</a>> wrote:<br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex">For the history of this: I first reported this shift back in 2018 <br>
(<a href="http://lists.maptools.org/pipermail/proj/2018-February/008082.html" rel="noreferrer" target="_blank">http://lists.maptools.org/pipermail/proj/2018-February/008082.html</a>). <br>
It was fixed for 3857 in the PROJ 5.2.0 release.<br>
<br>
My advice to you is to change to using EPSG:3857 as most software now no <br>
longer supports the non-standard "900913".<br>
<br>
-jeff<br>
<br>
<br>
<br>
-- <br>
Jeff McKenna<br>
MapServer Consulting and Training Services<br>
co-founder of FOSS4G<br>
<a href="http://gatewaygeo.com/" rel="noreferrer" target="_blank">http://gatewaygeo.com/</a><br>
<br>
<br>
<br>
<br>
<br>
On 2020-07-20 3:25 p.m., Shira Bezalel wrote:<br>
> We recently recompiled PostGIS in one of our Postgres 12.3 environments <br>
> with PROJ 6 (6.3.2) from PROJ 4 (4.9.3).<br>
> <br>
> postgis_full_version() output:<br>
> POSTGIS="3.0.1 ec2a9aa" [EXTENSION] PGSQL="120" GEOS="3.7.1-CAPI-1.11.1 <br>
> 27a5e771" PROJ="6.3.2" LIBXML="2.9.4" LIBJSON="0.12.1"<br>
> <br>
> We're now noticing that using ST_Transform(geometry, srid)**to transform <br>
> some points in California from SRID 4326 to 900913 is causing a shift <br>
> along the Y axis of about 20-30 meters. The difference is evidenced by a <br>
> comparison of the old and new ST_Y() values.<br>
> <br>
> If I force the PROJ 4 definition to be used, the transformation works <br>
> fine with no shift. (To force the PROJ 4 definition, I used the method <br>
> described in <br>
> <a href="http://blog.cleverelephant.ca/2019/02/proj4-postgis.html" rel="noreferrer" target="_blank">http://blog.cleverelephant.ca/2019/02/proj4-postgis.html</a> and set <br>
> authname, authsrid and srtext to null so that the proj4text column is used.)<br>
> <br>
> Is this shift expected with PROJ 6? I realize 900913 is deprecated and <br>
> we plan to migrate to 3857 as a longer term solution, but looking for <br>
> a shorter term fix if possible since much data and code is impacted by <br>
> this. Note we are seeing no shift issues with SRID 3857.<br>
> <br>
> Let me know if I can provide more information.<br>
> <br>
> Thanks much,<br>
> Shira<br>
> <br>
<br>
_______________________________________________<br>
postgis-users mailing list<br>
<a href="mailto:postgis-users@lists.osgeo.org" target="_blank">postgis-users@lists.osgeo.org</a><br>
<a href="https://lists.osgeo.org/mailman/listinfo/postgis-users" rel="noreferrer" target="_blank">https://lists.osgeo.org/mailman/listinfo/postgis-users</a></blockquote></div></div>