Skip to content
This repository has been archived by the owner on Apr 19, 2022. It is now read-only.

Commit

Permalink
updates
Browse files Browse the repository at this point in the history
  • Loading branch information
seisman committed Jun 20, 2017
1 parent 07b2020 commit 35cc589
Show file tree
Hide file tree
Showing 6 changed files with 89 additions and 38 deletions.
1 change: 0 additions & 1 deletion source/explain_-I.rst_
Original file line number Diff line number Diff line change
Expand Up @@ -10,7 +10,6 @@
- ``+n`` 表明 ``<xinc>`` 和 ``<yinc>`` 不是网格间隔,而是X和Y方向的节点数。
此时会根据节点数、网格区域范围以及网格配准方式重新计算网格间隔。


注意:

1. 若 ``<yinc>`` 设置为0,则表示其与 ``<xinc>`` 相同
Expand Down
30 changes: 21 additions & 9 deletions source/gmtinfo.rst
Original file line number Diff line number Diff line change
Expand Up @@ -74,18 +74,30 @@ gmtinfo
$ gmt info input.dat -C -I2
0 2 0 4 0 3.2

上面的例子中, ``-I`` 选项后只接了一个增量,此时第二列数据也会使用同样的增量(出于向后兼容性考虑)。若只想要为第一列指定增量,第二列不指定,可以使用 ``-Ip`` 选项::
上面的例子中, ``-I`` 选项后只接了一个增量,此时第二列数据也会使用同样的
增量(出于向后兼容性考虑)。若只想要为第一列指定增量,第二列不指定,可以使用
``-Ip`` 选项::

$ gmt info input.dat -C -Ip2
0 2 1 3 0 3.2

对于某些算法,比如FFT或者surface命令, ``-R`` 选项的不同取值会影响到运算的速度, ``-If`` 会输出优化后的 ``-R`` 以得到最快FFT结果, ``-Is`` 会输出优化后的 ``-R`` 以得到最快的surface结果。
对于某些算法,比如FFT或者surface命令, ``-R`` 选项的不同取值会影响到运算的
速度, ``-If`` 会输出优化后的 ``-R`` 以得到最快FFT结果, ``-Is`` 会输出优化
后的 ``-R`` 以得到最快的surface结果。

``-L``
确定多个表数据( ``-Af`` )或多段数据( ``-As`` )的共同范围。若与 ``-I`` 选项连用,则最终的结果会向内近似以使得极值位于真实数据范围内。

``-S``
见官方文档。

``-T<dz>[/<col>]``
以 ``-Tzmin/zmax/dz`` 的形式报告第 ``<col>`` 列的最值,默认报告第1列。该选项必须 ``-I`` 选项连用。
确定多个表数据( ``-Af`` )或多段数据( ``-As`` )的共同范围。若与 ``-I``
选项连用,则最终的结果会向内近似以使得极值位于真实数据范围内。

``-S[xy]``
Add extra space for error bars. Useful together with **-I** option
and when later plotting with :doc:`psxy` **-E**. **-Sx** leaves space
for horizontal error bars using the values in third
(2) column. **-Sy** leaves space for vertical error
bars using the values in third (2) column. **-S**
or **-Sxy** leaves space for both error bars using the values in
third and fourth (2 and 3) columns.

``-T<dz>[+c<col>]``
以 ``-Tzmin/zmax/dz`` 的形式报告第 ``<col>`` 列的最值, ``<col>`` 的默认值为0。
该选项必须 ``-I`` 选项连用。
20 changes: 16 additions & 4 deletions source/gmtwhich.rst
Original file line number Diff line number Diff line change
Expand Up @@ -6,18 +6,21 @@ gmtwhich
:官方文档: :ref:`gmt:gmtwhich`
:简介: 返回指定文件的完整路径

GMT会依次在三个目录中去寻找命令行中指定的文件::
GMT会依次在如下目录中去寻找命令行中指定的文件::

当前目录 > $GMT_USERDIR > $GMT_DATADIR
当前目录 > $GMT_USERDIR > $GMT_DATADIR > $GMT_CACHEDIR

该命令会报告文件的完整路径,使得用户可以确认自己在使用的究竟是哪个数据文件。

选项
----
必须选项
--------

``<files>``
任意数据文件名

可选选项
--------

``-A``
仅考虑用户有读权限的文件

Expand All @@ -26,3 +29,12 @@ GMT会依次在三个目录中去寻找命令行中指定的文件::

``-D``
不报告完整路径,仅打印包含该文件的目录名

``-G``
GMT可以自动下载某些文件:

- 以URL形式给出的文件会自动下载
- 以 ``@filename`` 形式给出的文件会从GMT网站上下载
- GMT提供的全球地形起伏数据 ``earth_relief_*.grd``

使用该选项,若GMT没有在本地缓存目录中找到这些文件,则会尝试下载。
23 changes: 18 additions & 5 deletions source/grdsample.rst
Original file line number Diff line number Diff line change
Expand Up @@ -12,20 +12,33 @@ grdsample
#. 不同的网格间隔或网格节点数( ``-I`` )
#. 不同的网格范围( ``-R`` )

选项
----
默认使用 bicubic 插值方式,可以使用 ``-n`` 选项设置其他插值方式。
该命令可以很安全地将粗网格插值为细网格;反之,将细网格插值为粗网格时,则可能
存在混淆效应,因而需要在插值前使用 ``grdfft`` 或 ``grdfilter`` 对网格文件做
滤波。

必须选项
--------

``<in_grdfile>``
要重采样的2D网格文件

``-G<out_grdfile>``
重采样生成的网格文件

``-I<xinc>[<unit>][=|+]/<yinc>[<unit>][=|+]``
指定X和Y方向的网格间隔,见 :doc:`xyz2grd` 中的介绍。
可选选项
--------

.. include:: explain_-I.rst_

``-R<w>/<e>/<s>/<n>``
指定新网格的数据范围。

若只使用 ``-R`` 选项,则等效于使用 ``grdcut`` 或 ``grdedit -S`` 。

``-T``
交换网格文件的配准方式。即若输入是网格线配准,则输出为像素点配准;若输入是像素点配准,则输出为网格线配准。
交换网格文件的配准方式。即若输入是网格线配准,则输出为像素点配准;若输入
是像素点配准,则输出为网格线配准。

``-n[b|c|l|n][+a][+b<BC>][+c][+t<threshold>]``
重采样时使用的插值算法,见 :ref:`doc:option-n` 一节。
Expand Down
30 changes: 23 additions & 7 deletions source/grdvolume.rst
Original file line number Diff line number Diff line change
Expand Up @@ -6,16 +6,22 @@ grdvolume
:官方文档: :ref:`gmt:grdvolume`
:简介: 计算网格数据中某个等值线所包围的表面积和体积

该命令读取一个2D网格文件,通过指定某条等值线确定某个Z值平面,并计算由该等值线约束的区域网格表面积、网格表面到该平面所包围的体积,以及最大平均高度(体积/面积)。也可以指定一系列等值线,此时该命令会分别计算每个等值线范围内的表面积和体积。
该命令读取一个2D网格文件,通过指定某条等值线确定某个Z值平面,并计算由该等值线
约束的区域网格表面积、网格表面到该平面所包围的体积,以及最大平均高度(体积/面积)。
也可以指定一系列等值线,此时该命令会分别计算每个等值线范围内的表面积和体积。

选项
----
必须选项
---------

``<grdfile>``
输入的2D网格文件名

可选选项
--------

``-C``
指定等值线(Z值平面)。
指定等值线(Z值平面),以计算由该等值线所决定的表面积、体积以及平均高度。
若不使用该选项,则返回整个网格文件的表面积、体积和平均高度。

该选项有四种不同的语法:

Expand All @@ -28,7 +34,9 @@ grdvolume
同时也计算从等值线到 ``<base>`` 的体积

``-S[<unit>]``
将弧度转换为Flat Earth距离,默认为笛卡尔坐标
对于地理网格,默认会将角度转换为 “Flat Earth” 下的距离,默认单位为米。
使用 ``-S`` 选项以指定其他的单位,则输出的表面积单位为此 ``<unit>^2`` ,
而输出的体积单位则是 ``<unit>^2 * z_unit`` 。

``-T[c|h]``

Expand All @@ -45,9 +53,9 @@ grdvolume

gmt grdvolume hawaii_topo.nc -Sk

计算网格表明与等值线Z=250之间的体积::
计算网格表面与等值线Z=250m之间的体积::

gmt grdvolume peaks.nc -Sk -C250
gmt grdvolume peaks.nc -Se -C250

在等值线100到300范围内,以10为间隔,计算所有等值线所约束的表面积和体积::

Expand All @@ -60,3 +68,11 @@ grdvolume
计算湖内从表面到300米深度范围内水的体积::

gmt grdvolume lake.nc -Cr-300/0

引用
----

如果本命令对你非常有用,可以引用如下文章:

Wessel, P., 1998, An empirical method for optimal robust regional-residual
separation of geophysical data, *Math. Geol.*, **30**\ (4), 391-408.
23 changes: 11 additions & 12 deletions source/index.rst
Original file line number Diff line number Diff line change
Expand Up @@ -46,6 +46,15 @@ GMT |version| 模块列表
- :doc:`gmtget` :列出单个或多个GMT参数的当前值
- :doc:`gmtdefaults` :列出所有GMT参数的当前值

- 信息提取

- :doc:`grdinfo` :从网格文件中提取基本信息
- :doc:`gmtinfo` :从表数据中提取信息
- :doc:`minmax <gmtinfo>`: 已改名为 ``gmtinfo``
- :doc:`gmtwhich` :返回指定文件的完整路径

- 格式转换

- 绘制 1D 数据

- :doc:`pssolar` :计算或/和绘制晨昏线以及民用、航海用以及天文用曙暮光区域
Expand All @@ -56,10 +65,6 @@ GMT |version| 模块列表

- :doc:`gmtlogo` :在图上绘制GMT图形logo

- 信息提取

- :doc:`grdinfo` :从网格文件中提取基本信息

- 1D 数据处理

- :doc:`gmtselect` :根据多个准则筛选数据
Expand All @@ -72,14 +77,11 @@ GMT |version| 模块列表
- :doc:`grdinfo` :从网格文件中提取基本信息
- :doc:`grdpaste` :将两个网格沿着其共同边界拼接成一个文件
- :doc:`grdraster` :从二进制数据中提取子区域并保存为GMT网格文件

- 格式转换
- :doc:`grdsample` :对网格文件做重采样
- :doc:`grdvolume` :计算网格数据中某个等值线所包围的表面积和体积

**以下模块需要重新检查**

- :doc:`gmtinfo` :从表数据中提取信息
- :doc:`minmax <gmtinfo>`: 已改名为 ``gmtinfo``
- :doc:`gmtwhich` :返回指定文件的完整路径
- :doc:`fitcircle` :拟合球面上数据点的平均位置及圆弧
- :doc:`pssac` :在地图上绘制 SAC 格式的地震波形数据
- :doc:`pscoast` :在地图上绘制海岸线、河流、国界线
Expand All @@ -105,13 +107,10 @@ GMT |version| 模块列表
- :doc:`grdblend` :将多个部分重叠的网格文件合并成一个网格文件
- :doc:`grdlandmask` :根据海岸线数据创建陆地-海洋的mask网格文件
- :doc:`grdtrend` :拟合网格的趋势面并计算残差
- :doc:`grdsample` :对网格文件做重采样
- :doc:`grdvolume` :计算网格数据中某个等值线所包围的表面积和体积
- :doc:`grdproject` :对网格数据做地图变换和逆变换
- :doc:`grdmask` :根据多边形数据或点数据创建mask网格文件
- :doc:`kml2gmt` :将Google Earth的KML文件转换为GMT表数据
- :doc:`grd2xyz` :将网格文件转换成表数据
- :doc:`grdconvert` :在不同的网格格式之间互相转换
- :doc:`psconvert` :将GMT生成的PS文件转换为其他图片格式
- :doc:`ps2raster <psconvert>` :将GMT生成的PS文件转换为其他图片格式
- :doc:`xyz2grd` :将XYZ数据或Z数据转换成网格文件
Expand Down

0 comments on commit 35cc589

Please sign in to comment.