光栅：如果坐标完全匹配替换值，如果不是 NA答案

【问题标题】：Raster: if the coordinates exactly match replace the values, if not NA光栅：如果坐标完全匹配替换值，如果不是 NA
【发布时间】：2021-06-24 01:24:01
【问题描述】：

我想用 (Class@data$layer) 中的新值替换实际的栅格 (r) 值。但是，我的条件是如果栅格r 中的坐标与Class 坐标（SpatialPointsDataFrame）完全匹配，则替换为新值（layer），并将所有坐标替换为NA。在我的例子中：

library(sp)  
library(dplyr)
library(raster)

# create raster
r <- matrix(rnorm(n=10000,mean=10),100,100)
r <- raster(r)
extent(r) <- c(544937,545916,7321332,7322524)
projection(r) <- "+proj=utm +zone=22 +ellps=WGS84 +units=m +no_defs"

# Selection of some coordinates that exists in the raster 
ras.coords <- as.data.frame(coordinates(r))
coords.sample <- ras.coords %>% sample_n(100)

# Create a new values
coords.sample$layer <- rpois(n=100, lambda=25)

# Convert  to SPDF
Class <- SpatialPointsDataFrame(coords = coords.sample[,1:2],
                           data = as.data.frame(coords.sample[,3]), 
                           proj4string = crs(r))
names(Class) <-("layer")

#Rewrite raster values, but the values are populated if exist the coordinate correspondence if not replace by NA

r[cellFromXY(r, Class@coords)] <- Class@data$class

显然不起作用，但它表达了我的目标。请帮忙解决我的问题？

【问题讨论】：

标签： r raster r-raster

【解决方案1】：

如果我理解正确，你可以这样做

x <- rasterize(Class, r, "layer")

【讨论】：