grdblend¶

官方文件:	grdblend
簡介:	將多個部分重疊的網格檔合併成一個網格檔

grdblend 模塊會讀取多個網格檔以及一個合併參數文件，並將多個網格檔合併成一個網格檔。

合併過程中會遇到如下幾種情況：

某個節點只有一個值：使用該值填充該節點
某個節點有多個值：對多個值計算其權重均值並賦值給該節點
某個節點沒有值：該節點的值由 -N 選項控制，默認爲NaN

輸入網格數據中，若網格節點配準方式或網格間隔與輸出數據不同，則會自動調用 grdsample 對這些網格數據進行重採樣。

語法¶

gmt grdblend [ blendfile | grid1 grid2 … ] -Goutgrid -Iincrement -Rregion [ -Cf|l|o|u[+n|p] ] [ -Nnodata ] [ -Q ] [ -Zscale ] [ -V[level] ] [ -W[z] ] [ -fflags ] [ -nflags ] [ -rreg ] [ --PAR=value ]

必選選項¶

grid1 grid2

要合併的網格檔

blendfile

合併參數文件。

文件中包含了要合併的網格檔的文件名列表，其格式爲:

網格檔名（必須）  -R選項（可選）  相對權重（可選）

在考慮權重的情況下，-R 範圍外的部分會給零權重，-R 範圍內的部分則會使用指定的相對權重，在邊界處會加上2D餘弦taper權重。
若相對權重爲負值，則整個反過來，即範圍內的會給零權重，範圍外的會給權重的絕對值。
-R 選項可以取值 -，此時會直接使用該網格檔的數據範圍。
若未指定權重，則默認權重爲1
也可以不指定 blendfile 而只將所有網格檔名列在命令行中，此時會使用網格數據的真實範圍，且所有網格檔的權重都是1

-Goutgrid: 輸出的網格檔名

-Ixinc[unit][+e|n][/yinc[unit][+e|n]]

指定X和Y方向的网格间隔

xinc X方向的网格间隔
yinc Y方向的网格间隔
unit 网格间隔的单位。对于地理坐标系而言默认值为度，也可设置为其他距离单位
+e 微调X和Y方向范围的最大值，使得其是网格间隔的整数倍（默认会微调网格间隔以适应给定的数据范围）
+n 表明 xinc 和 yinc 不是网格间隔，而是X和Y方向的节点数。此时会根据节点数、网格区域范围以及网格配准方式重新计算网格间隔。

注意：

若 yinc 设置为0，则表示其与 xinc 相同
若使用 -Rgrdfile 选项，则网格间隔和配准方式已经根据网格文件自动初始化，此时依然可以使用 -I 和 -r 覆盖相应的值

-Rxmin/xmax/ymin/ymax[+r][+uunit] (more …): 指定数据范围

可選選項¶

-Cf|l|o|u[+n|p]

合併網格檔時，若多個網格出現重疊，則忽略權重並按照下面的規則設置重疊區域的值：

-Cf 將第一個訪問該節點的網格數據作爲該節點的值
-Co 將最後一個訪問該節點的網格數據作爲該節點的值
-Cl 將所有網格檔在該節點的值的最小值作爲該節點的值
-Cu 將所有網格檔在該節點的值的最大值作爲該節點的值

對於 -Cf 和 -Co 而言，網格檔的順序決定了這些節點的值。在這種模式下，權重和餘弦tapering會被忽略。

使用 +n 或 +p 則首先將網格節點值初始化爲第一個網格檔的值，對於接下來的其它網格檔，僅當其值小於等於0或大於等於0時才考慮是否更新該節點的值。

-Nnodata: 將無數據的節點的值設置爲 nodata，默認值爲NaN

-Q: 創建一個無頭段的二進制網格檔以供其它程序使用

-V[level] (more …): 设置 verbose 等级 [c]

-W[z]

不合並數據，僅輸出每個節點所使用的權重。

-Wz 表示輸出 Z值乘以權重的和。

-Zscale: 在輸出前先將數據乘以比例因子 scale，默認值爲 1

-f[i|o]colinfo (more …): 指定输入或输出列的数据类型

-n[b|c|l|n][+a][+bBC][+c][+tthreshold] (more …): 设置网格文件的插值方式

-r[g|p] (more …): 设置网格配置方式 [默认为网格线配准]

-^ 或 -: 显示简短的帮助信息，包括模块简介和基本语法信息（Windows下只能使用 -）
-+ 或 +: 显示帮助信息，包括模块简介、基本语法以及模块特有选项的说明
-? 或无参数: 显示完整的帮助信息，包括模块简介、基本语法以及所有选项的说明
--PAR=value: 临时修改GMT参数的值，可重复多次使用。参数列表见 gmt.conf

示例¶

假設要合併幾個網格檔，則可以設置合併參數文件 blend.job 的內容爲:

piece_1.nc -R0/30/-90/90 1
piece_2.nc -R25/50/-90/90 1.5
piece_3.nc -R45/80/-90/90 0.9
piece_4.nc -R80/160/-90/90 1

執行如下命令即可實現數據合併:

gmt grdblend blend.job -Gblend.nc -R0/160/-90/90 -I1m/1m -V

將所有網格檔 MB_*.nc 以相同權重合並:

gmt grdblend MB_*.nc -Gblend.nc -R0/360/-90/90 -I1m/1m -V

grdblend¶

語法¶

必選選項¶

可選選項¶

示例¶

相關模塊¶