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