|||
在实际绘图中,对数据网格化后,会产生位于由原始数据点构成的测网外部的点的数据,这部分数据是不可靠的,往往需要屏蔽掉。本部分的目的是,根据测网外边界构成的闭合多边形产生网格化数据文件(其Z值为赋值),由此文件和根据原始数据生成的网格化文件做积,生成新的网格化文件,该文件特点是:在多边形外部可以为0或者NAN(依赖于不同的赋值情况,也可为其它值),在多边形内部为某类值(如原始值,或其倍数),在多边形边上为某值(如原始值,或其倍数)。
grdmask:生成可以形成闭合多边形的网格。可以设置多边形外部、内部和边上点对应的数值,一般与grdmath组合使用。
grdmask pte.txt -Gfactor.nc -R100/119/28/45 -I0.05/0.05 -NNAN/1/1
grdmath factor.nc test.grd MUL = test_m.grd
grdimage test_m.grd -R100/119/28/45 –Ba5 -Ccolor.cpt –Q >gmt.ps
说明:
grdmask语句中,pte.txt为构成多边形的节点文件,-G指示输出文件名,-NNAN/1/1设置多边形外部点处的值为NAN,多边形的边上及其内部为1。
grdmath语句,两个网格化文件做乘积,由于factor.nc文件的特点,该语句的结果即生成了在多边形外部为NAN,在多边形边上和多边形内部为原test.grd文件中对应点处的数值,没有变化。
Archiver|手机版|科学网 ( 京ICP备07017567号-12 )
GMT+8, 2024-11-26 09:58
Powered by ScienceNet.cn
Copyright © 2007- 中国科学报社