Add segy2grd (#1428)

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

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`