zoukankan      html  css  js  c++  java
  • [转载]grdcontour命令在GMT4下绘制等值线图

    原始数据文件 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

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


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

    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%

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

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

    写的比较简单,如有问题,欢迎指出。
  • 相关阅读:
    java
    java
    java
    java
    java
    java
    java
    java
    sed命令的用法
    linux系统产生随机数的6种方法
  • 原文地址:https://www.cnblogs.com/gisalameda/p/12840533.html
Copyright © 2011-2022 走看看