PostGISのspatial_ref_sys用世界測地系データ

GIS

　PostGISに入っているspatial_ref_sysテーブルに測地成果2000の測地系が入っていないので、Proj4から必要な部分だけ抜き出してInsert文を作ってみました。

一応

test=# select transform('SRID=2449;POINT(100 100)',4612);
transform

                                                                                                      • -

SRID=4612;POINT(137.167775884152 36.0009013187198)
(1 row)

と
http://vldb.gsi.go.jp/sokuchi/surveycalc/xy2blf.html
緯度、経度への換算
で確かめました

9/16日修正・追加

この日記の時点で未完成だったので、もうちょっとまともなSQL文を作りました。


測地成果2000のSQL文

このSQL文では、世界測地系（JGD2000)の平面直角座標 I系 〜 XIX系および世界測地系（JGD2000)の緯度経度の空間参照IDが定義しています。
※PostGIS 0.9の時点では、世界測地系はspatial_ref_sysテーブルに含まれていません.


使い方例
    

        #psql -d geometry_databese -f proj_jpn.txt

10/23追加
SRIDについては
id:hogemanさんのエントリが参考になります
http://d.hatena.ne.jp/hogeman/20041022