pygmt - 如何在 GMT6.1 中将海平面设置为一个常见的负值？

Question

我想画一些区域的浮雕（涉及海洋和陆地），但我不想显示海洋测深的变化，所以我必须将它们设置为一个共同的负值，以确保所有测深映射到相同颜色，如浅蓝色。我不知道该怎么做？

Answer 1

这很容易使用grdclip。使用命令行 GMT 时，您可以尝试：

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

它将输入网格的所有值grid< 0 设置为 -1000 （根据需要调整值）并将剪裁的网格写入 outfile grid_clipped。

如果您使用的是 PyGMT，您可以尝试：

import pygmt 

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

您可以在其中直接将剪裁的网格交给另一个 PyGMT 函数进行绘图，例如 fig.grdimage(grid = grid).