windows下GMT绘制等值线

set plot=ave.ps
gmt gmtinfo 2017050600.txt

rem 查看文件信息，方便制作cpt，文件格式为 x y z

gmt makecpt -Cjet -T20/60/5 > std.cpt
pscontour 2017050600.txt -R110/118/20/26 -JM6i -Cstd.cpt -Bafg -P -I  >>%plot%

::pscontour 2017050600.txt -R110/118/20/26 -JM6i -Bafg -Wthin -Cstd.cpt -Lthinnest,- -Gd1i>> %plot%

::pscontour 2017050600.txt -R110/118/20/26 -JM6i -Bafg -Wthin -Cstd.cpt >> %plot%

rem 三种类型

psconvert %plot% -A -Tg -P

rem 保存为png图片

del *.history
del *.conf
del *.ps
pause

结果：