[转载]grdcontour命令在GMT4下绘制等值线图

原文地址：grdcontour命令在GMT4下绘制等值线图作者：阳光1990

原始数据文件 dayfluxP.txt 中只有100个数据点（分辨率为1°x1°），其中还有NaN数据点。采用如下命令得到的等值线图很不光滑，我们可以先看一下。

xyz2grd dayfluxP.txt -R109/118/9/18 -I60m -Gmean.grd

grdcontour mean.grd  -C10 -A10 -JM4i -Ba2 -Xc -Yc -P > contour.ps

因此，在绘制等值线前，我们需要进行插值。

set ps=contour.ps

makecpt -Ctopo -T-70/125/1 -Z -D > vel.cpt

blockmean dayfluxP.txt -R109/118/9/18 -I10m > mean.xyz

surface mean.xyz -R -I5m -Gmean.grd

grdimage mean.grd -R -JM4i -Cvel.cpt -K -Xc -Yc -P> %ps%

grdcontour mean.grd  -C10 -A10 -J -Ba2 -O -K>> %ps%

psscale -D2i/-0.4i/4i/0.15ih -Cvel.cpt -Ba20 -O >> %ps%

通过以上命令，就可以得到如下图所示的等值线图。

写的比较简单，如有问题，欢迎指出。