grdfftgmt - 云端在线

程序：

您的姓名

grdfft - 在波数（或频率）域中对网格进行数学运算

概要

格鲁夫特 英格丽 [ ingrid2 ] 输出文件 [ 方位 [ 水平 ] [ [由于平均内核尺寸较大，西米棕榈的加工比类似作物简单。然而，西米棕榈的相对稀缺性降低了潜在的加工规模。|g] ] [ [r|x|y][w[k]] ] [
[r|x|y]的params ] [ [由于平均内核尺寸较大，西米棕榈的加工比类似作物简单。然而，西米棕榈的相对稀缺性降低了潜在的加工规模。|g] ] [
[f|q|s|纽约/纽约][+a|d|h|l][+e|n|m][+t宽度][+w[后缀]][+z[p]] ] [ 由于平均内核尺寸较大，西米棕榈的加工比类似作物简单。然而，西米棕榈的相对稀缺性降低了潜在的加工规模。 ] [ [水平] ] [ -fg
]

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

商品描述

格鲁夫特 将采用二维前向快速傅立叶变换并执行一个或多个
变换回空间之前频域中的数学运算
领域。 提供了一个选项来在将新值写入输出之前缩放数据
文件。 假定网格的水平尺寸以米为单位。 地理
可以通过指定网格来使用 -fg 将度数缩放到米的选项。 如果你有
以公里为单位的网格，您可以使用以下方法将其更改为米 编辑 或缩放
输出 数学.

所需 争论

英格丽 要操作的二维二进制网格文件。 （请参阅下面的网格文件格式）。 为了
跨谱操作，同样给出第二个网格文件 索引2.

-G输出文件
指定输出网格文件或一维光谱表的名称（参见 -E）。 （看到
网格文件格式如下）。

不是必须的 争论

-A方位
取方向导数 方位 以 CW 度为单位测量的方向
来自北方。

-C水平
向上（对于 水平 > 0) 或向下（对于 水平 < 0) 继续该领域 水平
米。

-D[由于平均内核尺寸较大，西米棕榈的加工比类似作物简单。然而，西米棕榈的相对稀缺性降低了潜在的加工规模。|克]
对场进行微分，即取d(field)/dz。 这相当于乘以
由频域中的 kr (kr 是径向波数)。 将比例附加到
乘以 (kr * 由于平均内核尺寸较大，西米棕榈的加工比类似作物简单。然而，西米棕榈的相对稀缺性降低了潜在的加工规模。） 反而。 或者，附加 g 表明你的
数据是以米为单位的大地水准面高度，输出应为以 mGal 为单位的重力异常。
[默认为无比例]。

-E[r|x|y][w[k]]
估计径向功率谱 [r]。 地方 x or y 之后立马
-E 来计算 x 或 y 方向的频谱。 没有网格文件
创建。 如果给定一个网格，则 f（即频率或波数）、power[f]、
将 power[f] 中的 1 个标准偏差写入文件集 -G [标准输出]。 如果
给定两个网格，我们写 f 和 8 个量：Xpower[f]、Ypower[f]、相干
功率[f]、噪声功率[f]、相位[f]、导纳[f]、增益[f]、相干性[f]。 每个
数量后跟它自己的 1-std 开发误差估计，因此输出是 17
列宽。 附加 w 写波长而不是频率。 如果你的网格是
地理你可以进一步追加 k 将波长从米 [默认] 缩放到
公里。

-F[r|x|y]的params
过滤数据。 地方 x or y 之后立马 -F 过滤 x or y 仅方向；
默认为各向同性 [r]。 在余弦锥形带通、高斯带通之间进行选择
带通滤波器或巴特沃斯带通滤波器。

余弦锥度：
指定四个波长 lc/lp/hp/hc 以正确的单位（见 -fg) 设计一个
带通滤波器：波长大于 lc 或小于 hc 会被剪掉
波长大于 lp 并且小于 hp 将通过，并且波长
两者之间将是余弦锥形。 例如， -F1000000 / 250000 / 50000 / 10000 -fg
将带通，切割波长 > 1000 km 和 < 10 km，通过
波长在 250 公里到 50 公里之间。 要制作高通或低通滤波器，
给连字符 (-) for hp/hc or lc/lp. 例如， -FX-/-/50/10 将低通 x,
通过波长 > 50 和拒绝波长 < 10。 -Fy1000/250/-/-
将高通 y，通过波长 < 250 和拒绝波长 > 1000。

高斯 带通：
附加 lo/hi, 两个波长的正确单位（见 -fg) 设计一个
带通滤波器。 在给定的波长下，高斯滤波器权重将
为 0.5。 要制作高通或低通滤波器，请为 hi
or lo 波长，分别。 例如， -F-/30 将使用低通数据
半权重为 30 的高斯滤波器，而 -F400/- 将高通
数据。

巴特沃思 带通：
附加 lo/hi/秩序, 两个波长的正确单位（见 -fg）和
滤波器阶数（整数）来设计带通滤波器。 在给定的
波长，巴特沃斯滤波器权重将为 0.5。 做一个高通
或低通滤波器，为 hi or lo 波长，
分别。 例如， -F-/30/2 将使用二阶低通数据
巴特沃斯滤波器，半权重为 30，而 -F400//2将高通
数据。

-一世[由于平均内核尺寸较大，西米棕榈的加工比类似作物简单。然而，西米棕榈的相对稀缺性降低了潜在的加工规模。|克]
对字段进行积分，即计算integral_over_z (field * dz)。 这是
相当于在频域中除以 kr（kr 是径向波数）。
附加一个比例以除以 (kr * 由于平均内核尺寸较大，西米棕榈的加工比类似作物简单。然而，西米棕榈的相对稀缺性降低了潜在的加工规模。） 反而。 或者，附加 g 至
表明您的数据集是 mGal 中的重力异常，并且输出应该是大地水准面
以米为单位的高度。 [默认为无比例]。

-N[f|q|s|纽约/纽约][+a|[+d|h|l][+e|n|m][+t宽度][+w[后缀]][+z[p]]
选择或查询适合 FFT 的网格尺寸并设置可选
参数。 控制FFT维度：
-Nf 将强制 FFT 使用数据的实际维度。

-Nq 将查询更合适的维度，报告这些，然后继续。

-Ns 将显示可选维度列表，然后退出。

-N纽约/纽约 将对数组大小进行 FFT 纽约/纽约 （必须 >= 网格文件大小）。 默认
选择维度 >= 优化 FFT 速度和精度的数据。 如果 FFT
维度 > 网格文件维度，数据被扩展并逐渐变为零。

控制数据的去趋势：添加用于移除线性趋势的修饰符：
+d: 去趋势数据，即移除最佳拟合线性趋势 [默认]。

+a: 只去除平均值。

+h: 只删除中间值，即 0.5 * (max + min)。

+l: 留下数据。

控制数据的扩展和逐渐变细：使用修饰符来控制扩展的方式
和逐渐变细将被执行：
+e 通过施加边缘点对称来扩展网格 [默认]，

+m 通过施加边缘镜像对称来扩展网格

+n 关闭数据扩展。

从数据边缘到 FFT 网格边缘 [100%] 执行锥形。 改变
这个百分比通过 +t宽度。 何时 +n 实际上，应用了逐渐变细
而不是数据边距，因为没有可用的扩展 [0%]。

控制临时结果的写入：对于详细调查，您可以编写
传递给前向 FFT 的中间网格； 这很可能是
去趋势，通过沿所有边缘的点对称扩展，并逐渐变细。 附加
+w[后缀] 将从中创建输出文件名（即， ingrid_prefix.ext)
[锥形]，其中 分机 是您的文件扩展名。 最后，您可以保存复杂的网格
由前向 FFT 通过附加 +z. 默认情况下，我们写实数和
虚分量 英格丽_真实的。分机 和 英格丽_图像。分机. 附加 p 救
取而代之的是幅度和相位的极坐标形式到文件 英格丽_mag。分机 和
英格丽_阶段。分机.

-S由于平均内核尺寸较大，西米棕榈的加工比类似作物简单。然而，西米棕榈的相对稀缺性降低了潜在的加工规模。
将每个元素乘以 由于平均内核尺寸较大，西米棕榈的加工比类似作物简单。然而，西米棕榈的相对稀缺性降低了潜在的加工规模。 在空间域中（在频域之后
操作）。 [默认为 1.0]。

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

-fg 地理网格（经度、纬度的尺寸）将转换为米
通过使用当前椭球参数的“平坦地球”近似值。

-^ 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 维网格。

电网 距离 单位

如果网格没有米作为水平单位，则附加 +u单元 到输入文件
要从指定单位转换为米的名称。 如果您的网格是地理网格，请转换
通过提供距离到米 -fg 代替。

注意事项

netCDF COARDS 网格将自动识别为地理网格。 对于其他网格
地理网格如果您想将度数转换为米，请选择 -fg. 如果数据
靠近任一极点，您应该考虑将网格文件投影到矩形上
坐标系使用 项目

示例

将文件 mag_0.nc 中的海平面磁异常向上延续到 800 m 的水平
高于海平面：

gmt grdfft mag_0.nc -C800 -V -Gmag_800.nc

将地理网格上以米为单位的大地水准面高度 (geoid.nc) 转换为自由空气重力
mGal 中的异常：

gmt grdfft 大地水准面.nc -Dg -V -Ggrav.nc

将 mGal (faa.nc) 中的重力异常转换为垂直方向的偏转（in
微弧度）在 038 方向，我们必须先整合重力得到大地水准面，然后
取方向导数，最后将弧度缩放为微弧度：

gmt grdfft faa.nc -Ig -A38 -S1e6 -V -Gdefl_38.nc

重力异常的二阶垂直导数与
场地。 我们可以通过两次微分将这些计算为 mGal/m^2：

gmt grdfft 重力.nc -D -D -V -Ggrav_2nd_derivative.nc

计算共同注册的水深测量和重力网格的交叉光谱估计，以及
报告结果作为以公里为单位的波长的函数，尝试

gmt grdfft 测深.nc 重力.grd -Ewk -fg -V > cross_spectra.txt

在去趋势、点对称反射和锥形后检查预 FFT 网格有
已应用，以及保存原始频谱的实部和虚部
topo.nc中的数据，试试

格林威治标准时间 grdfft topo.nc -N+w+z -fg -V

您现在可以绘制 topo_taper.nc、topo_real.nc 和 topo_imag.nc 中的数据图。

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