diff --git a/source/module/index.rst b/source/module/index.rst index 643b5126e9..0b5a816d71 100644 --- a/source/module/index.rst +++ b/source/module/index.rst @@ -96,6 +96,7 @@ - :doc:`sac` - :doc:`sample1d` - :doc:`segy` + - :doc:`segy2grd` - :doc:`solar` - :doc:`spectrum1d` - :doc:`sph2grd` @@ -159,7 +160,6 @@ - :doc:`gmt:greenspline` - :doc:`gmt:movie` - :doc:`gmt:plot3d` - - :doc:`gmt:supplements/segy/segy2grd` - :doc:`gmt:supplements/segy/segyz` - :doc:`gmt:supplements/spotter/backtracker` - :doc:`gmt:supplements/spotter/gmtpmodeler` @@ -379,7 +379,7 @@ :doc:`gmt:supplements/gsfml/fzprofiler` Bash script to plot fracture zone cross-profiles :doc:`gmt:supplements/gsfml/mlconverter` Convert chrons to ages using selected magnetic timescale **其他模块** .. _module_other: -:doc:`gmt:supplements/segy/segy2grd` Converting SEGY data to a GMT grid +:doc:`segy2grd` 将 SEGY 文件转换为网格文件 :doc:`grdgdal` 在 GMT 中运行 GDAL 命令 :doc:`gshhg` 从 GSHHG 或 WDBII 数据文件中提取数据 :doc:`img2google` 由测深墨卡托 img 网格创建谷歌地球 KML 文件 @@ -487,6 +487,7 @@ sac sample1d segy + segy2grd solar spectrum1d sph2grd diff --git a/source/module/segy2grd.rst b/source/module/segy2grd.rst new file mode 100644 index 0000000000..c06685e6ac --- /dev/null +++ b/source/module/segy2grd.rst @@ -0,0 +1,97 @@ +.. index:: ! segy2grd +.. include:: common_SYN_OPTs.rst_ + +segy2grd +============= + +:官方文档: :doc:`gmt:supplements/segy/segy2grd` +:简介: + +**segy2grd** 用于读取 IEEE SEGY 文件并创建二进制网格文件。它可以执行简单的映射(相当于 xyz2grd -Z),或者在单个网格单元包含 SEGY 文件中多个样本值的情况下进行更复杂的平均值计算。 + +如果某些节点没有填充数据,**segy2grd** 会进行报告。这些未被约束的节点会被设置为用户指定的值[默认值为 NaN]。对于包含多个值的节点,将被设为这些值的平均值。 + +语法 +-------- + +**gmt segy2grd** +*segyfile* +|-G|\ *grdfile* +|SYN_OPT-I| +|SYN_OPT-R| +[ |-A|\ [**n**\|\ **z**] ] +[ |-L|\ [*nsamp*] ] +[ |-M|\ [*ntraces*] ] +[ |-Q|\ **x**\|\ **y**\ *value* ] +[ |-S|\ [*header*] ] + +必须选项 +------------------ + +*segyfile* + *segyfile* 是一个 IEEE 浮点数格式的 SEGY 文件。所有的道集都假定从 0 时间/深度域开始。 + +.. _-G: + +**-G**\ outgrid + 输出的网格文件名 + +.. _-I: + +**-I**\ *x_inc*\ /[*y_inc*] + *x_inc* [以及可选的 *y_inc*] 是网格间距。可以在数值后添加 **m** 表示分钟,或添加 **s** 表示秒。 + +可选选项 +------------------ + +.. _-A: + +**-A**\ [**n**\|\ **z**] + 将属于同一节点(相同坐标)的多个值相加(等同于 **-Az**)。 + 添加 **n** 则仅统计分配到每个节点的数据点数量。 + [默认情况下(不使用 **-A** 选项)会计算平均值]。该选项不用于简单映射。 + +.. _-L: + +**-L**\ [*nsamp*] + 使用 *nsamp* 覆盖每道的采样点数 + +.. _-M: + +**-M**\ [*ntraces*] + 设定要读取的道数。默认尝试读取10000道。 + **-M**\ 0 将读取二进制头中指定的道数, + **-M**\ *ntraces* 将尝试仅读取 *n* 道。 + +.. _-Q: + +**-Q**\ **x**\|\ **y**\ *value* + 根据不同指令可用于更改两种不同的设置: + **-Qx**\ *x-scale* 对道头中的坐标应用缩放因子 *x-scale* ,使其匹配 -R 中指定的坐标范围 + + **-Qy**\ *s_int* 当 SEGY 文件中的采样间隔不正确时,指定采样间隔为 *s_int* 。可重复使用。 + +.. _-S: + +**-S**\ [*header*] + 设置可变间距。 *header* 可以是 **c** 表示 CDP(共同深度点),**o** 表示偏移量, + 或 **b**\ *number* 表示从第 *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 + +相关模块 +-------- + +:doc:`grd2xyz`, +:doc:`segy`, +:doc:`xyz2grd`