如何用NA替换对应于多边形的值？

Question

我有一个栅格和一个 shapefile：

library(cartography)
library(sf)
library(raster)

r <- raster(matrix(rnorm(10*12), nrow=10), xmn = -180, xmx= 180, ymn = -90, ymx= 90)
mtq <- st_read(system.file("gpkg/mtq.gpkg", package="cartography"), quiet = TRUE)

我想将栅格 r 与 shapefile mtq 相交，并将所有多边形的相应像素设为 NA（用 NA 替换栅格中的像素值）和return 光栅。

Answer 1

您可能正在寻找 mask；它同时存在于旧 {raster} 和闪亮的新 {terra}.

中

请注意，我不得不稍微重写您的 r 对象，因为它与 {cartography} 中的马提尼克矢量对象不太兼容。

编辑：如果正如评论中所指出的那样，您正在寻找用 NA 替换多边形内部（而不是外部）的值，我的答案仍然是 raster::mask()，只是稍作调整掩蔽对象的（您需要光栅范围内多边形的倒数）。

library(cartography)
library(sf)
library(raster)

mtq <- st_read(system.file("gpkg/mtq.gpkg", package="cartography"), quiet = TRUE) %>% 
  dplyr::summarise() # dissolve internal boundaries

r <- raster(matrix(rnorm(10*12), nrow=10), 
            xmn = st_bbox(mtq)["xmin"], 
            xmx= st_bbox(mtq)["xmax"], 
            ymn = st_bbox(mtq)["ymin"], 
            ymx= st_bbox(mtq)["ymax"], 
            crs = st_crs(mtq))

plot(r) # original raster - full extent + range

# the masking object:
mask <- st_bbox(r) %>% # take extent of your raster... 
  st_as_sfc() %>% # make it a sf object
  st_set_crs(st_crs(mtq)) %>% # in CRS of your polygon 
  st_difference(mtq) %>% # intersect with the polygon object
  st_as_sf() # interpret as sf (and not sfc) object

result <- r %>% 
  mask(mask)

plot(result)
plot(st_geometry(mtq), add = T)

如何用NA替换对应于多边形的值？

How to replace values that correspond to polygons by NA?

r

raster

sf

r-raster