Question

Source Code

地球掩膜数据

地球掩膜数据指明了地表的水陆地理特征。该网格数据的Z值只有5个值。这5个值从0到4依次指代海洋（ocean）、陆地（land）、湖泊（lake）、湖中岛（island-in-lake）和岛中湖（pond-in-island-in-lake）。

该数据的主要作用是对网格数据进行掩膜操作，以孤立出仅位于陆地或海洋的数据。该掩膜数据可以用 grdlandmask 模块生成，但对于高分辨率、全精度 GSHHG数据而言，这一过程通常非常耗时。因而GMT官方提供了计算好的掩膜数据供用户使用。该数据也分不同的分辨率，从15秒到1度不等。

下表列出了GMT提供的地球掩膜数据 earth_mask 的代码、配准方式、分辨率以及文件大小。

代码	配准方式	分辨率	大小
01d	g,p	1 度	20 KB
30m	g,p	30 分	24 KB
20m	g,p	20 分	30 KB
15m	g,p	15 分	35 KB
10m	g,p	10 分	57 KB
06m	g,p	6 分	106 KB
05m	g,p	5 分	117 KB
04m	g,p	4 分	149 KB
03m	g,p	3 分	209 KB
02m	g,p	2 分	393 KB
01m	g,p	1 分	998 KB
30s	g,p	30 秒	2.6 MB
15s	g,p	15 秒	7.5 MB

• 代码中 d、m和 s 表示度、分和秒

• 配准方式中 g 和 p 分别表示网格线配准和像素配准

使用方法

当需要使用地球掩膜数据时，可以通过 @earth_mask_rru[_reg] 的形式调用。其中 rru 是网格文件的分辨率代码（上表第一列），reg 是要使用的配准方式代码（上表第二列）。配准方式代码是可选的，若不指定，则默认返回像素配准的网格数据。如果命令中使用了 -R 选项，则只会读取该区域内的掩膜数据。

查看30分的地球掩膜数据的信息:

gmt grdinfo @earth_mask_30m

如果用户希望仅仅获得海洋或者水陆的掩膜，需要使用 grdmath 模块对掩膜数据进行调整。下面给出一些例子。

利用 grdmath 将掩膜数据 earth_mask_01m_g 中大于（GT）0 的为 1, 其余为 0。这就意味着海洋部分改为0，其余均为1:

gmt grdmath @earth_mask_01m_g 0 GT = ocean_land.grd

生成一个像素配准的掩膜数据，其陆地上为NaN，海洋中为1，可以用:

gmt grdmath @earth_mask_01m_p 0 LE 0 NAN = land_NaN.grd

生成一个像素配准的掩膜数据，其湖泊中值为1，其余地方均为NaN:

gmt grdmath @earth_mask_01m 2 EQ 0 NAN = lakes_and_NaN.grd

技术细节

• 数据会自动下载并保存到 ~/.gmt/server/earth/earth_mask 目录下