ST_Transformどうなの?


これからの流れ
https://jpusergis.wordpress.com/2017/11/15/donotusentv2injp/

当然TKY2JGDとは違います。
コード書けない人はどーしてんのか。

POSTGIS=2.3.2
GEOS=3.6.1-CAPI-1.10.1
PROJ=Rel. 4.9.1
GDAL=GDAL 2.1.3
LIBXML=2.7.8
LIBJSON=0.12
RASTER

緯度 35.003215692
経度 134.997229989
X座標 -110380.5661 m
Y座標 60597.9468 m

TOKYO 4301

SELECT
 ST_AsText(
  ST_Transform(
   ST_GeomFromText('POINT(135 35)', 4301),
   4612
  )
 );

“POINT(134.997204249013 35.0031971806904)”
X座標 -110382.6352 m
Y座標 60595.6111 m

並べると
X座標 -110380.5661 m
X座標 -110382.6352 m
Y座標 60597.9468 m
Y座標 60595.6111 m

緯度 24.467
経度 122.98
緯度 24.472484139
経度 122.980455228
X座標 -168814.9237 m
Y座標 -103351.4901 m

SELECT
 ST_AsText(
  ST_Transform(
   ST_GeomFromText('POINT(122.98 24.467)', 4301),
   4612
  )
 );

“POINT(122.978487384933 24.4711542616074)”

X座標 -168960.7563 m
Y座標 -103552.0726 m

並べると
X座標 -168814.9237 m
X座標 -168960.7563 m
Y座標 -103351.4901 m
Y座標 -103552.0726 m

まあ補正パラメータはなぁ。。。。

カテゴリー: 開発, 設計 タグ: パーマリンク