风雨的个人博客分享 http://blog.sciencenet.cn/u/xinxiang

博文

GMT绘图:grdmask、grdmath与grdimage

已有 19623 次阅读 2014-1-23 23:08 |个人分类:GMT|系统分类:科研笔记| GMT, 屏蔽, 多边形

在实际绘图中,对数据网格化后,会产生位于由原始数据点构成的测网外部的点的数据,这部分数据是不可靠的,往往需要屏蔽掉。本部分的目的是,根据测网外边界构成的闭合多边形产生网格化数据文件(其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文件中对应点处的数值,没有变化。






https://blog.sciencenet.cn/blog-387-761732.html

上一篇:Kindle PW安装Koreader解决pdf等文档的重排版问题
下一篇:GMT绘图 输出特殊字符
收藏 IP: 123.138.56.*| 热度|

0

该博文允许注册用户评论 请点击登录 评论 (8 个评论)

数据加载中...

Archiver|手机版|科学网 ( 京ICP备07017567号-12 )

GMT+8, 2023-1-30 12:42

Powered by ScienceNet.cn

Copyright © 2007- 中国科学报社

返回顶部