r - 是否有一个 R 例程来绘制测深栅格图，仅选择从 0 到 -3000 的测深值，不包括高于零的值？

Question

我有一个那不勒斯湾的栅格，其中包含海洋扩张的水深（即负值）值和陆地的高度数据（即正值）。我想从地图中排除零（高度）以上的测深正值，我该怎么做？

标签： r

Answer 1

这是一个最小的、独立的、可重现的示例

library(anglr)
library(raster)
geb <- anglr::gebco
plot(geb)

解决方案

x <- reclassify(geb, cbind(0, Inf, NA))
plot(x)