GMT6.1如何设置海平面为负值?

How to set the elevations of sea to a common negative value in GMT6.1?

我想绘制某个区域的地势(涉及海洋和陆地),但我不想显示海洋测深的变化,所以我必须将它们设置为共同的负值以确保所有测深图颜色相同,如浅蓝色。不知道怎么办?

使用 grdclip 很容易。使用命令行 GMT 时,您可以尝试:

gmt grdclip grid -Ggrid_clipped -Sb0/-1000 -V

将输入网格的所有值设置为 grid< 0 到 -1000(根据需要调整值)并将剪切的网格写入输出文件 grid_clipped.

如果您使用的是 PyGMT,您可以尝试:

import pygmt 

grid = pygmt.grdclip(grid, below = [0, -1000])

你可以直接将裁剪的网格交给另一个 PyGMT 函数进行绘图 fig.grdimage(grid = grid).