[转载]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%
写的比较简单,如有问题,欢迎指出。