gis - 带有Y坐标的postgis/proj 900913至4326投影问题

Question

这是个老歌，但我似乎找不到解决方法。

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

例：



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



这里原始的900913声明点被st_s转换为4326然后返回900913。结果不是原始点，y有所不同。 （我稍后会插入结果，我在这里没有）。

我尝试更改4326的proj4text，添加+ nadgrids = @ null，就像我在某处阅读的一样

srid 4326的proj4text当前为：


“从spatial_ref_sys中选择proj4text，其中srid = 4326”
+ proj = longlat + ellps = WGS84 +基准= 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的直接转换完全相同的点。

任何人有任何想法吗？

电子杂志

Answer 1

EPSG：900913是不确定的投影。您应该改用EPSG：3857，它应该是完全相同的Spherical Mercator，但是标准化的。

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