grdfiltergmt - 云端在线

程序：

您的姓名

grdfilter - 在空间（或时间）域中过滤网格

概要

过滤器 英格丽 距离标志 x宽度[/宽2][修饰符] 外格 [ 增量 [ i|p|r ]
[ 地区 ] [ ] [水平] ] [ -f]

请注意： 选项标志和相关参数之间不允许有空格。

商品描述

过滤器 将使用选定的卷积之一在时域中过滤网格文件
或非卷积各向同性或矩形滤波器并使用笛卡尔计算距离
或球形几何形状。输出网格文件可以选择生成为子区域
输入的（通过 -R）和/或新的增量（通过 -I）或注册（通过 -T）。 在这个
方式，输入数据中可能有“额外空间”，这样边缘就不会被使用，并且
输出可以在输入边缘的二分之一宽度内。如果滤波器是低通的，
那么输出的采样频率可能低于输入。

所需 争论

英格丽 要过滤的点的网格文件。 （参见下面的网格文件格式）。

-D距离标志
距离 旗 告诉网格 (x,y) 如何与过滤器相关 宽度 如下：

旗 = p: 网格 (px,py) 与 宽度 奇数个像素；笛卡尔距离。

旗 = 0: 网格 (x,y) 与单位相同 宽度, 笛卡尔距离。

旗 = 1：以度为单位的网格（x，y）， 宽度 以公里为单位，笛卡尔距离。

旗 = 2：以度为单位的网格（x，y）， 宽度 以 km 为单位，dx 按 cos(middle y) 缩放，笛卡尔
距离。

上述选项最快，因为它们只允许计算权重矩阵
一次。 接下来的三个选项较慢，因为它们会重新计算每个选项的权重
纬度。

旗 = 3：以度为单位的网格（x，y）， 宽度 以公里为单位，dx 由余弦（y）缩放，笛卡尔
距离计算。

旗 = 4：以度为单位的网格（x，y）， 宽度 以公里为单位，球面距离计算。

旗 = 5：墨卡托网格 (x,y) -Jm1 个图像单位， 宽度 以公里为单位，球面距离
计算。

-FX宽度[/宽2][修饰符]
设置过滤器类型。在卷积和非卷积滤波器之间进行选择。使用任何
过滤器代码 x （如下所列）后跟全直径 宽度。这给出了一个
各向同性滤光片；追加/宽2 对于矩形过滤器（需要 -DP or -D0).
某些过滤器允许使用可选参数和修饰符。

卷积滤波器（及其代码）是：

(b) 棚车：所有的重量都是相等的。

(c) 余弦拱：权重遵循余弦拱曲线。

(g) 高斯：权重由高斯函数给出，其中 宽度 是6倍
传统的高斯西格玛。

(f) 自定义：权重由过滤器权重网格中预先计算的值给出
文件 重量，其尺寸必须为奇数；还要求 -D0 和输出间距
必须与输入间距匹配或者是整数倍。

(o) 运算符：权重由过滤器权重网格中预先计算的值给出
文件 重量，其尺寸必须为奇数；还要求 -D0 和输出间距
必须与输入间距匹配或者是整数倍。假设权重总和为
零，因此不会进行权重总和和标准化的累积。

非卷积滤波器（及其代码）是：

(m) 中值：返回中值。选择另一个分位数追加 +q分位数 in
0-1范围[默认为0.5，即中位数]。

(p) 最大似然概率（模态估计器）：返回模态值。 如果更多
发现不止一种模式，我们返回它们的平均值。附加 +l or +u 如果你宁愿
想要返回最低或最高的模态值。

(h）直方图模式（另一种模式估计器）：返回以模态值为中心
直方图中的主峰。追加/箱宽 指定分箱
间隔。使用修饰符 +c 将垃圾箱集中在多个 箱宽 [默认有
bin 边缘是 箱宽]。如果发现不止一种模式，我们
返回它们的平均值。附加 +l or +u 如果您想退回
模态值的最低或最高。

(l) Lower：返回所有值的最小值。

(L) Lower：仅返回所有正值中的最小值。

(u) Upper：返回所有值中的最大值。

(U) 上部：仅返回最大值或所有负值。

在案件 L|U 可能没有数据通过初始符号测试； 在
在这种情况下，过滤器将返回 NaN。

-G外格
外格 是过滤器的输出网格文件。 （参见下面的网格文件格式）。

不是必须的 争论

-I新奇[单元][=|+][/inc[单元][=|+]]
x_公司 [和可选 y_公司] 是输出增量。 附加 m 来表示弧线
分钟，或 s 来指示弧秒。如果新的 x_公司, y_公司 不是整数
旧的倍数（在输入数据中），过滤将相当大
慢点。 [默认值：与输入相同。]

-Ni|p|r
确定输入网格中的 NaN 值如何影响过滤后的输出：追加 i 至
在计算过滤值时忽略所有 NaN [默认]， r 与...相同 i
除非输入节点为 NaN，则输出节点将设置为 NaN（仅
如果两个网格共同注册则适用），并且 p 这将强制过滤值
如果在过滤器圆内发现任何具有 NaN 值的网格节点，则为 NaN。

-R 西, 东, 南及 北 定义输出点的区域。 [默认：
与输入相同。]

-T 切换输出网格的节点注册，以便与
输入网格 [默认提供与输入网格相同的注册]。

-V[水平] （更多的 ...）
选择详细级别 [c]。

-f[我|o]信息 （更多的 ...）
指定输入和/或输出列的数据类型。

-^ or 只是 -
打印一条关于命令语法的短消息，然后退出（注意：在 Windows 上
只用 -).

-+ or 只是 +
打印广泛的使用（帮助）消息，包括对任何
模块特定选项（但不是 GMT 通用选项），然后退出。

-? or 没有 参数
打印完整的使用（帮助）消息，包括选项的解释，然后
退出。

- 版
打印 GMT 版本并退出。

--显示数据目录
打印 GMT 共享目录的完整路径并退出。

电网 文件 FORMATS

默认情况下，GMT 在 COARDS-complaint netCDF 中将网格写出为单精度浮点数
文件格式。 但是，GMT 能够在许多其他常用网格中生成网格文件
文件格式，也有助于所谓的“打包”网格，写出浮点数
数据为 1 或 2 字节整数。 要指定精度、比例和偏移量，用户应该
添加后缀 =id[/规模/抵消[/南]]， 在哪里 id 是网格的两个字母的标识符
类型和精度，以及 规模 和 抵消 是可选的比例因子和偏移量
应用于所有网格值，以及 南 是用于指示缺失数据的值。 如果
两个字符 id 不提供，如 =/规模 比一个 id=nf 假设。 什么时候
阅读格子，格式一般是自动识别的。 如果不是，相同的后缀
可以添加到输入网格文件名中。 看 转换 和节网格文件格式
GMT 技术参考和说明书以获取更多信息。

当读取包含多个网格的 netCDF 文件时，GMT 将默认读取
可以在该文件中找到的第一个二维网格。 哄GMT阅读另一个
网格文件中的多维变量，追加 ?变量名 到文件名，其中
变量名 是变量的名称。 请注意，您可能需要转义特殊含义
of ? 在 shell 程序中通过在它前面放置一个反斜杠，或者通过放置
引号或双引号之间的文件名和后缀。 这 ?变量名 也可以使用后缀
为输出网格指定一个与默认值不同的变量名：“z”。 看
转换 和 GMT Technical 的 Sections modifiers-for-CF 和 grid-file-format
参考和食谱以获取更多信息，特别是关于如何阅读 3-、
4 维或 5 维网格。

地理 AND 时间 坐标

当输出网格类型为 netCDF 时，坐标将标记为“经度”，
“纬度”或“时间”基于输入数据或网格（如果有）的属性或基于
-f or -R 选项。 例如，两者 -f0x -f1t 和 -R90w/90e/0t/3t 将导致
经度/时间网格。 当 x、y 或 z 坐标为时间时，将存储在网格中
作为 TIME_UNIT 和 TIME_EPOCH 中指定的纪元以来的相对时间 配置文件 文件
或在命令行上。 除此之外 单元 时间变量的属性将指示
这个单位和时代。

示例

假设north_pacific_etopo5.nc是5E到140E的260分钟测深文件
0N 到 50N，并且您想要找到 300 公里半径（600 公里）内的值的中位数
全宽）的输出点，您可以选择从 150E 到 250E 和 10N 到 40N，
您想要每 0.5 度的输出值。使用球面距离计算，
你需要：

gmt grdfilter North_pacific_etopo5.nc -Gfiltered_pa​​cific.nc -Fm600 \
-D4 -R150/250/10/40 -I0.5 -V

如果我们想要一个高通结果，那么可以执行相应的低通结果
使用粗网格间隔进行过滤，因为 grdfilter 会将结果重新采样为相同的结果
分辨率作为输入网格，以便我们可以计算残差，例如

gmt grdfilter Northern_pacific_etopo5.nc -Gresidual_pacific.nc -Fm-600 \
-D4 -R150/250/10/40 -I0.5 -V

此处，residual_pacific.nc 网格将具有与原始网格相同的 5 分钟分辨率。

使用自定义各向异性高斯滤波器 exp 过滤ripples.nc中的数据集
(-0.5*r^2)，其距中心的距离 r 由 (2x^2 + y^2 -2xy)/6 给出，主要
轴与水平面成 63 度角，尝试

gmt grdmath -R-10/10/-10/10 -I1 X 2 POW 2 MUL Y 2 POW ADD XY MUL 2 MUL \
SUB 6 DIV NEG 2 DIV EXP DUP SUM DIV = gfilter.nc
gmt grdfilter Ripples.nc -Ffgfilter.nc -D0 -Gsmooth.nc -V

限制

要使用 -D5 选项输入墨卡托网格必须由 img2mercgrd 使用以下命令创建 -C
选项，因此 y 值的原点是赤道（即 x = y = 0 对应于 lon =
纬度 = 0)。

使用 onworks.net 服务在线使用 grdfiltergmt