gis - postgis / proj 900913 から 4326 への Y 座標の投影の問題

Question

これは古き良きものですが、解決策が見つからないようです。

900913 座標で 4326 システムに st_transform を実行したい場合、y 座標がシフトします。

例：

SELECT AsText( Transform( Transform( GeomFromText( 'POINT( 449760.25168159 6790560.4594059 )', 900913), 4326 ), 900913 ) )

ここで、元の 900913 開始点は st_stransformed で 4326 になり、900913 に戻ります。結果は元の点ではなく、y が異なります。(後で結果を挿入します。ここにはありません)。

どこかで読んだように +nadgrids=@null を追加して、4326 の proj4text を変更してみました

現在、srid 4326 の proj4text は次のとおりです。

"srid=4326 の spatial_ref_sys から proj4text を選択" +proj=longlat +ellps=WGS84 +datum=WGS84 +no_defs

現在、srid 900913 の proj4text は次のとおりです。

"srid=900913 の spatial_ref_sys から proj4text を選択" +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