jgdtrans#
Coordinate Transformer by Gridded Correction Parameter which Geospatial Information Authority of Japan (GIAJ, formerly GSIJ) distributing [1].
国土地理院が公開しているパラメータファイル(par ファイル)による座標変換(順逆変換)を提供します。
Features:
Supports offline transformation (no web API)
Supports both original forward and backward transformation
Supports verified backward transformation
Supports all TKY2JGD, PatchJGD, PatchJGD(H), HyokoRev, SemiDynaEXE and POS2JGD (geonetF3 and ITRF2014)
Clean implementation
Example:
import jgdtrans
with open('SemiDyna2023.par') as fp:
tf = jgdtrans.load(fp, format='SemiDynaEXE')
# Geospatial Information Authority of Japan
origin = (36.10377479, 140.087855041, 2.34)
# forward transformation
result = tf.forward(*origin)
# prints Point(latitude=36.103773017086695, longitude=140.08785924333452, altitude=2.4363138578103)
print(result)
# backward transformation
point = tf.backward(*result)
# prints Point(latitude=36.10377479, longitude=140.087855041, altitude=2.34)
print(point)
# backward transformation compatible to GIAJ web app/APIs
p = tf.backward_compat(*result)
# prints Point(latitude=36.10377479000002, longitude=140.087855041, altitude=2.339999999578243)
print(p)
This package does not contain parameter files, the GIAJ distributes such files [2].
このパッケージはパラメータファイルを提供しません。公式サイトよりダウンロードしてください [2]。
We use TKY2JGD for Windows Ver.1.3.79 [3] as a reference.