地理データ

都道府県庁データ (緯度　経度　県庁名)

pref_capital.txt

Bashスクリプト　(県庁所在地を地図上にプロット。psxyは最初の2列しかデータを読まないのでこれで問題なし)

#!/bin/bash
PROJ='M15'
RANGE='127/146/26/46'
MARK='c0.2'
OUT='japanmap.eps'
OUTpng='japanmap.png'
gmtdefaults -D > .gmtdefaults4
gmtset PLOT_DEGREE_FORMAT +D D_FORMAT %.1f
gmtset PAPER_MEDIA A4+
pscoast -J$PROJ -R$RANGE -Ba4f2 -Di -W2 -P -K > $OUT
psxy japan.txt -J$PROJ -R$RANGE -S$MARK -W3 -Gblack -O >> $OUT
convert -quality 100 $OUT $OUTpng