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

博文

GMT绘图--grdfilter滤波

已有 5231 次阅读 2016-1-11 11:03 |个人分类:GMT|系统分类:科研笔记| grdfilter

GMT绘图----grdfilter滤波
Grdfilter 对 grid文件进行空间(时间)域滤波
grdfilter ingrid -Ddistance_flag -Fxwidth[/width2][modifiers]  -Goutgrid  [ -Iincrement ]  [ -Ni|p|r ]  [ -Rregion ] [ -T ] [ -V[level] ] [ -f<flags> ]
使用选取的褶积或非褶积各向同性或矩形滤波器对grid文件进行滤波并计算笛卡尔或球面几何。两种几何下的距离计算公式不同。
-D distance_flag。指明网格(x,y)和滤波宽度的关系
    flag = p: grid (px,py) with width an odd number of pixels; Cartesian distances.
    flag = 0: grid (x,y) same units as width, Cartesian distances.
    flag = 1: grid (x,y) in degrees, width in kilometers, Cartesian distances.
    flag = 2: grid (x,y) in degrees, width in km, dx scaled by cos(middle y), Cartesian distances.
    以上几个选项计算过程是快的,下面几个由于要在每个纬度上重算权重,所以计算起来较慢。
    flag = 3: grid (x,y) in degrees, width in km, dx scaled by cosine(y), Cartesian distance calculation.
    flag=4: 网格以度为单位,滤波器宽度为km,球面距离。
    flag = 5: grid (x,y) in Mercator -Jm1 img units, width in km, Spherical distance calculation.

-Fxwidth[/width2][modifiers]。指定滤波器类型及滤波宽度。一般为各向同性的滤波,只有在附上width2时才为矩形滤波。
   目前GMT支持的滤波有褶积滤波和非褶积滤波:
   褶积滤波有:
    -b Boxcar等权重滤波,即直径范围内的节点数据的权重是相同的。
    -c 余弦曲线滤波,权重遵循余弦曲线。
    -g 高斯滤波,权重由高斯函数决定。
 此外,还有自定义滤波函数、非卷积类型的滤波
-Goutgrid 输出文件
-Ixinc[unit][=|+][/yinc[unit][=|+]] 输出文件中的网格间隔。若单位为分或秒,需要附上m或s。若新的网格   间隔x_inc, y_inc不是输入文件中旧间隔的整数倍,滤波过程会非常慢。默认是与输入相同。
-Ni|p|r 指明输入grid文件中NaN值对滤波的影响方式。i:忽略所有NaN值(默认);r:忽略所有NaN值,不 过在为NaN值的节点上,输出的值中该节点仍为NaN (only applies if both grids are co-registered);

  p:在滤波半径内若存在NaN值则该节点的滤波后的值将强制设置为NaN



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

上一篇:GMT绘图 输出特殊字符
收藏 IP: 124.114.234.*| 热度|

0

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

数据加载中...

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

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

Powered by ScienceNet.cn

Copyright © 2007- 中国科学报社

返回顶部