GMT对陆地和海洋数据的识别和赋值操作

小二

很多情况下，我们的数据同时分布在海洋和陆地上的，有时候需要剔除海洋数据或者陆地数据，或者对海洋或陆地的数据进行不同的操作，GMT可以通过内置的海岸线和水文数据库实现上面的操作。

例如为了实现海洋区域数据的准确统计，我们要设置陆地数据为0，GMT可以通过两个程序的组合来实现：

grdlandmask -R -Dh -I20m -V -N1/0 -Glandmask.nc

grdmath  landmask.nc globe_hz18_ku_seaice_bscat_2.nc MUL = globe_hz18_ku_seaice_bscat_22.nc

首先是grdlandmath依照海岸线识别陆海，并设置陆地和水域数值，-N1/0 表示将陆地的数据设置为1，水域的数据设置为0，结果输出一个同网格密度的二进制文件（netCDF格式），该网格文件中的陆地数据全部为0，海洋网格点数据全部为1；然后对grdlandmask输出的网格文件和待处理的实际数据进行乘法操作，新输出的网格文件中陆地的数据就全部设置为0了，而海洋数据保持不变。接下来即可对其进行下一步的数学统计操作。

figure： distribution of sigma after landmask setting land values to 0.

−N可以设置为wet/dry二个类型，也可以设置ocean/land/lake/island/pond五个参数，下图书设置-N1/0/0/0/0，将内陆湖泊、岛屿和水塘都屏蔽之后的结果。

figure： distribution of sigma after landmask setting all continental  values to 0.

如果只提取出海洋范围内的数据，也可以通过下面更简单的方式实现：

gmtselect 098.nc -R0/360/-82.4/-50 -bi -bo -Nk/s/s/s/s > 098_lm.nc

其中-N的含义和grdlandmask的设置相同，gmtselect的作用顾名思义是数据的选择，它提供多种数据选择策略，包括基于内置海岸线，依点位置进行数据的选择提取方法。得到的结果文件中只含有通过区域的数据，那些不符合筛选条件的数据被过滤，因此结果文件已经不是格网数据，不能对其进行格网操作，但是数学统计是可以继续执行的。