GoGPT Best VPN GoSearch

OnWorks 网站图标

grdfftgmt - 云端在线

通过 Ubuntu Online、Fedora Online、Windows 在线模拟器或 MAC OS 在线模拟器在 OnWorks 免费托管服务提供商中运行 grdfftgmt

这是命令 grdfftgmt 可以使用我们的多个免费在线工作站之一在 OnWorks 免费托管服务提供商中运行,例如 Ubuntu Online、Fedora Online、Windows 在线模拟器或 MAC OS 在线模拟器

程序:

您的姓名


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


免费服务器和工作站

下载 Windows 和 Linux 应用程序

Linux 命令

Ad




×
广告
❤️在这里购物、预订或购买——免费,有助于保持服务免费。