postgis/proj 900913至4326投影問題與Y座標

Question

這是一個oldie，但我似乎無法找到解決方案。

當我想要一個900913座標系的一個st_transform到一個4326系統時，y座標轉換。

例如：

SELECT AsText（ 變換（ 變換（ GeomFromText（ 'POINT（449760.25168159 6790560.4594059）'，900913），）） ）

這裏的原始900913開始點是st_stransformed到4326並返回到900913.結果不是原始點，y不同。 （我會稍後插入結果，我沒有在這裏）。

我試圖改變proj4text爲4326，加入+ nadgrids = @空像我某處

的proj4text爲SRID 4326讀目前：

「選擇從spatial_ref_sys proj4text其中SRID = 4326」 + PROJ = longlat + ellps = WGS84 +基準= WGS84 + no_defs

爲SRID 900913的proj4text目前：

「塞萊ct proj4text from spatial_ref_sys where 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

你使用的是什麼版本的PostGIS？我試過這個配置的查詢：POSTGIS="1.5.2" GEOS="3.2.2-CAPI-1.6.2" PROJ="Rel. 4.7.1, 23 September 2009"，它工作正常。

儘管如此，我遇到了一些涉及900913座標的問題。

Answer 2

EPSG：900913是不明確的投影。你應該使用EPSG：3857來代替，它應該是完全相同的球形墨卡託，但是是標準的。

任何重投影是有損操作。來回轉換座標會使最後的二進制數字產生噪音，這通常是亞毫米誤差。