ジオイド

日本のジオイド2011

ここから日本のジオイドのデータ(gsigeo2011_ver2_1_asc.zip)をダウンロードできる。

その中の説明書（asc取扱説明書.pdf）によると,

「ジオイド・モデルは、北緯20度から50度、東経120度から150度の範囲について、緯度間隔1分、経度間隔1.5分の格子位置におけるジオイド高として与えられています。 テキスト（ASCII）形式のジオイド・モデル‘gsigeome_ver5.asc’は、１行目がヘッダ行であり、２行目以降に格子点でのジオイド高が格納されています。」

# ジオイド高（楕円体を基準にしたジオイドの高さ）の, GRASSへのインポート
$ grass    # start GRASS in latlon projection
$ unzip gsigeoid_ver5_asc.zip 
$ cd gsigeoid_ver5_asc/program
$ cat <<EOF > dummy
north:     50.0083333
south:     19.9916667
east:     150.0125
west:     119.9875
rows:    1801
cols:    1201
null:     999.0000
EOF
$ cat gsigeome_ver5.asc | \
 awk 'NR>1{for (i=1; i<=NF; i++)print $i}' | awk '$1>-999{print $1}' | \
 awk '{printf"%s ",$1}NR%1201==0{printf"\n"}' | \
 awk '{a[NR]=$0}END{for (i=0;i<NR;i++)print a[NR-i]}'  >> dummy
$ g.region w=119 e=151 s=19 n=51 nsres=0.0166666667 ewres=0.025
$ r.in.ascii input=dummy output=t --o

# このファイルは, 南から北に向けてデータが入っている。
# しかも1行は28データ(1ラインの一部)なので, 
# 単純に行順をひっくり返せば良いというものではない。
# いったん1ライン1データにして, ゴミ（改行だけの行）
# を消し, あらためて1ライン1行にまとめなおしてから行をひっくり返す。