r - 如何在R-raster中获得网格周围的等高线？

Question

在 R 中有一个栅格，我如何在网格周围绘制轮廓线(不连接中心或其他任何东西，真正遵循网格的边界)具有某些值(或由某些掩码标识)？

以下示例显示了如何获取值为 0.6 的区域周围的等高线:如何执行相同的操作，但使用遵循网格边界的线？
该函数应返回一个对象以添加到绘图中(作为 SpatialLinesDataFrame 的 rasterToContour )，并且相邻的网格应包含在一条等高线中(即，只应绘制多边形的外边界)。我找不到 rasterToPolygons 的解决方案(有关视觉方面的信息，请参阅 here，但在这里对我没有帮助)。

set.seed(2)
r <- raster(nrow=10, ncol=10)
r[] <- runif(ncell(r))
r[r>0.6] <- 0.6
rc <- rasterToContour(r, levels=c(0.6))
plot(r)
plot(rc, add=TRUE)

Answer 1

我会使用 clump() 的组合和 rasterToPolygons() :

library(raster)
library(rgeos)  ## For dissolve = TRUE in rasterToPolygons()

## Recreate your data
set.seed(2)
r <- raster(nrow = 10, ncol = 10)
r[] <- runif(ncell(r))   
plot(r)

## Compute and then plot polygons surrounding cells with values greater than 0.6
SP <- rasterToPolygons(clump(r > 0.6), dissolve = TRUE)
plot(SP, add = TRUE)