如何在R中给栅格中的特定数字一个不同的颜色（或形状）？

Question

我有一个栅格，该栅格的值为9999几个像素。 在绘制整个栅格之前，我想给这些像素某种颜色。 因此，图例未考虑此值

  Library(raster)
  filename <- system.file("external/test.grd", package="raster")
  r <- raster(filename)
  plot(r) ### normal plot which takes into account all pixels

  r[r>1000]=9999 

   plot(r)

在此，图例（图）不应考虑9999 ，而应给此值指定颜色（或形状）并按正常方式绘制r 。 9999的图例也可以分开

Answer 1

您可以创建不带9999个像元的栅格副本，以及仅带9999个像元的栅格副本，并将其覆盖：

library(raster)
filename <- system.file("external/test.grd", package="raster")
r <- raster(filename)
r[r>1000]=9999 

# raster without 9999
r2 = reclassify(r, matrix(c(1000, Inf, NA), ncol=3))
plot(r2, colNA = NA)

# raster with only 9999
r3 = reclassify(r, matrix(c(-Inf, 1000, NA, 1000, Inf, 9999), ncol=3, byrow=T))
plot(r3, add=T, col='black', colNA = NA, legend=F)