Question

官方文档

grdlandmask

简介

根据海岸线数据创建陆地-海洋的mask网格文件

grdlandmask 模块读取指定的海岸线数据，用于确定网格内的每个节点是位于陆地还是水域，并给不同类型的节点赋予不同的值，以生成掩膜文件。生成的掩膜文件可进一步用在 grdmath 中以掩盖掉位于陆地或水域中的数据点。

语法

gmt grdlandmask

-Gmask_grd_file

生成的掩膜网格文件的文件名

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

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

• xinc 和 yinc 为 X 和 Y 方向的网格间隔。对于地理坐标，可以指定网格间隔单位 [默认单位为度]

• +e 微调X和Y方向范围的最大值，使得其是网格间隔的整数倍（默认会微调网格间隔以适应给定的数据范围）

• +n 表明 xinc 和 yinc 不是网格间隔，而是X和Y方向的节点数。此时会根据节点数、网格区域范围以及网格配准方式重新计算网格间隔。

注意：

• 若 yinc 设置为0，则表示其与 xinc 相同

• 若使用 -Rgrdfile 选项，则网格间隔和配准方式已经根据网格文件自动初始化，此时依然可以使用 -I 和 -r 覆盖相应的值

-Rxmin/xmax/ymin/ymax[+r][+uunit] (more …)

指定数据范围

可选选项

grdlandmask 生成的掩膜网格文件属于 分类型 数据，即所有数据只能取几个固定的值，比如 -N0/1 会将水域内的网格值设置为0，将陆地内的网格值设置为1。在这种情况下，对这种数据用标准方法（比如样条）进行插值通常会得到无意义的结果，使用时需要小心。

然而，当你直接绘制该掩膜网格文件时，网格数据会被重新投影使得在投影后的坐标下变成一个矩形。这个过程中涉及到了网格插值，默认使用的插值算法是样条插值，因而可能会在图中产生假象。因而建议在使用 grdimage 绘制此类数据时使用 -nn 选项即 nearest neighbor 插值算法以避免这一问题。

earth_mask 远程数据

用 grdlandmask 模块生成高分辨率、全精度的掩膜数据会非常耗时。GMT官方提供了计算好的掩膜数据供用户直接使用。该数据分不同的分辨率，从15秒到1度不等。详情见 earth_mask: 地球掩膜数据。

示例

将所有陆地上的节点设置为NaN，水域上的节点设置为1:

gmt grdlandmask -R-60/-40/-40/-30 -Dh -I5m -N1/NaN -Gland_mask.nc -V

生成全球1x1度的网格，并将不同性质的区域设置成不同的值:

gmt grdlandmask -R0/360/-90/90 -Dl -I1 -N0/1/2/3/4 -Glevels.nc -V

相关模块

grdmath, grdclip, mask, clip, coast