GMT 绘制台站分布图

set ps=test.ps
set J=M4i
set R=73/135.5/10/54
rem gmt gmtset  FONT_ANNOT_PRIMARY 8p FONT_TITLE 8p 
gmt gmtset FORMAT_GEO_MAP ddd:mm
gmt gmtset  MAP_FRAME_TYPE fancy
gmt gmtset  MAP_FRAME_WIDTH 0.1c
gmt gmtset GRID_CROSS_SIZE_PRIMARY = 0c 
gmt gmtset TICK_LENGTH -0.0c
gmt gmtset ANNOT_FONT_PRIMARY = Helvetica 
rem 写入文件头
gmt psxy -J%J% -R%R% -T -K > %ps%
rem 真正的绘图命令
rem gmt grdcut topo30.grd -Gchina.grd  -R$R
rem gmt grdgradient china.grd  -A50 -Gtemp.grd  -Ne0.6 -M
gmt grdimage china.grd -Itemp.grd -J -R  -Ba10f5WSne   -Cwiki-scotland.cpt -K >%ps%
gmt psxy shengjie.xy -R -J   -K -O  -W0.3/255/0/0 >>%ps% 
gmt psxy station.txt -R -J -St0.2 -O -K -P -G0/0/255 >>%ps%
gmt psxy earthquake.txt -R -J -Sc0.2 -O -K -P -G255/0/0 >>%ps%
 gawk "{ print $1,$2,8,0,36,3,$3 }" station.txt| pstext -J -R -K -O  -D0.0c/0.1c >>%ps%
  gawk "{ print $1,$2,8,0,36,3,$3 }" zz.txt| pstext -J -R -K -O  -D0.0c/0.1c >>%ps%
 gawk "{ print $1,$2,8,0,36,3,$3 }" earthquake.txt| pstext -J -R -K -O  -D0.0c/0.2c >>%ps%
 gmt pslegend mylegend.txt -J -R -Dx8.3/0.2/1.6/1.2/BL -F -O -K -C2p/2p  >> %ps%
rem gmt pslegend mylegend.txt -J -R -Dx8.3/0.2/1.6/1.2/BL -F+gwhite -O -K -V >> %ps%
rem 写入文件尾
gmt psxy -J%J% -R%R% -T -O >> %ps%
gmt psconvert %ps% -Tj -E300  -P -C-sFONTPATH=c:\windows\Fonts
del gmt.conf gmt.history