【必要なデータ】 K-NET,KiK-NET観測地点の緯度、経度をそれぞれのホームページから頂いた。k-net.dat, kik-net.datにはそれぞれ緯度、経度のデータが２列になって入っている。県境のデータはhttp://www.geocities.jp/ne_o_t/index.htm から頂いた。
【プログラム】 ▲がK-NET、△がKIK-NETです。
#!/bin/bash
# 日本地図を書く
pscoast -R126/147/25/46 -Jm1.5 -B1 -P -W -Dh -X1 -Y1 -K > k_kik_location.ps
# 県境を書く
psxy ken.txt -W1ta -R -Jm -B -P -M -O -K >> k_kik_location.ps
# k-net観測点のプロット
awk '{print $1,$2}' k-net.dat | psxy -R -Jm -B -P -M -ST.1 -O -K -G0/0/0 >> k_kik_location.ps
# kik-net観測点のプロット
awk '{print $1,$2}' kik-net.dat | psxy -R -Jm -B -P -M -ST.1 -O -K >> k_kik_location.ps
