postgis / proj 900913到4326投影问题与Y坐标

Question

这是一个老人，但我似乎无法找到解决方案。

当我想在900913坐标上对4326系统进行st_transform时，y坐标会移动。

示例：

选择 AsText（ 转变（ 转变（ GeomFromText（'POINT（449760.25168159 6790560.4594059）'，900913）， 4326 ） 900913 ） ）

这里原始的900913说明点是st_stransformed到4326并返回到900913.结果不是原点，y不同。 （我将在稍后插入结果，我在这里没有。）

我尝试改变4316的proj4text，添加+ nadgrids = @ null就像我在某处读到

srid 4326的proj4text目前是：

“从spatial_ref_sys中选择proj4text，其中srid = 4326”  + proj = longlat + ellps = WGS84 + datum = WGS84 + no_defs

srid 900913的proj4text目前是：

 “从spatial_ref_sys中选择proj4text，其中srid = 900913”  + proj = merc + a = 6378137 + b = 6378137 + lat_ts = 0.0 + lon_0 = 0.0 + x_0 = 0.0 + y_0 = 0 + k = 1.0 + units = m  + nadgrids = @ null + no_defs

我也尝试过从900913到另一个投影到4326的投影，但我得到与从900913到4326的直接转换完全相同的点。

任何想法？

EJ

Answer 1

EPSG：900913是不明确的预测。您应该使用EPSG：3857，它应该是完全相同的球形墨卡托，而是标准的。

任何重投都是一种有损操作。来回转换坐标会产生最后二进制数字的噪声，这通常是亚毫米误差。

Answer 2

您使用的是什么版本的PostGIS？我已尝试使用此配置查询：POSTGIS="1.5.2" GEOS="3.2.2-CAPI-1.6.2" PROJ="Rel. 4.7.1, 23 September 2009"并且它可以正常工作。

尽管如此，我还有一些涉及900913坐标的问题。