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