segy2grd
- 官方文档:
- 简介:
segy2grd 用于读取 IEEE SEGY 文件并创建二进制网格文件。它可以执行简单的映射(相当于 xyz2grd -Z)，或者在单个网格单元包含 SEGY 文件中多个样本值的情况下进行更复杂的平均值计算。
如果某些节点没有填充数据，segy2grd 会进行报告。这些未被约束的节点会被设置为用户指定的值[默认值为 NaN]。对于包含多个值的节点，将被设为这些值的平均值。
语法
gmt segy2grd segyfile -Ggrdfile -Iincrement -Rregion [ -A[n|z] ] [ -L[nsamp] ] [ -M[ntraces] ] [ -Qx|yvalue ] [ -S[header] ]
必须选项
- segyfile
segyfile 是一个 IEEE 浮点数格式的 SEGY 文件。所有的道集都假定从 0 时间/深度域开始。
- -Goutgrid
输出的网格文件名
- -Ix_inc/[y_inc]
x_inc [以及可选的 y_inc] 是网格间距。可以在数值后添加 m 表示分钟，或添加 s 表示秒。
可选选项
- -A[n|z]
将属于同一节点(相同坐标)的多个值相加(等同于 -Az)。 添加 n 则仅统计分配到每个节点的数据点数量。 [默认情况下(不使用 -A 选项)会计算平均值]。该选项不用于简单映射。
- -L[nsamp]
使用 nsamp 覆盖每道的采样点数
- -M[ntraces]
设定要读取的道数。默认尝试读取10000道。 -M0 将读取二进制头中指定的道数， -Mntraces 将尝试仅读取 n 道。
- -Qx|yvalue
- 根据不同指令可用于更改两种不同的设置:
-Qxx-scale 对道头中的坐标应用缩放因子 x-scale ，使其匹配 -R 中指定的坐标范围
-Qys_int 当 SEGY 文件中的采样间隔不正确时，指定采样间隔为 s_int 。可重复使用。
- -S[header]
设置可变间距。 header 可以是 c 表示 CDP(共同深度点)，o 表示偏移量， 或 bnumber 表示从第 number 字节开始的4字节浮点数。 如果未设置 -S，则假定样本按照 -I 提供的 x_inc , y_inc 进行均匀间隔。
示例
从一个等间距的 SEGY 文件 test.segy 创建网格文件，读取每道上 18-25 秒时间域(或千米深度域)的数据，第一道会被假定位于 X=198 的位置:
gmt segy2grd test.segy -I0.1/0.1 -Gtest.nc -R198/208/18/25
从 SEGY 文件 test.segy 创建网格文件，根据 CDP 编号定位道集，其中每千米有 10 个 CDP，采样间隔为 0.1。由于网格间隔大于 SEGY 文件的采样间隔，因此单个样本将在网格单元内进行平均:
gmt segy2grd test.segy -Gtest.nc -R0/100/0/10 -I0.5/0.2 -Qx0.1 -Qy0.1