从 R 中的 shapefile 创建栅格答案

【问题标题】：Creating raster from shapefile in R从 R 中的 shapefile 创建栅格
【发布时间】：2021-02-10 16:30:45
【问题描述】：

我有一个研究区域的 shapefile，我想创建一个栅格，其中属于该 shapefile 的每个单元格的 value=1，其他单元格的 value=NA。

我想创建这个新栅格，其原点和像元大小与我将在分析中使用的环境栅格相同。

这是我正在使用的代码，但我总是收到此错误：“rep(value, length.out = 2) 中的错误：尝试复制 'S4' 类型的对象”

#Create raster from shapefile

library(raster)
library(rgdal)

x <- readOGR(dsn="D:/Mestrado/Tese/AreaEstudo/final", layer="areaestudofinaluniforme")
pad <- raster("D:/Mestrado/Tese/Dados modelação/Mais recente/Variáveis maxent/bio_11.asc")

ext <- extent(pad)
reso <- res(pad)
ori <- origin(pad)

newraster <- raster(x, ext, reso, ori)

（我还没有弄清楚如何将值分配给栅格）

有人对我为什么会收到此错误有任何建议吗？以及如何分配我想要的值？

【问题讨论】：

见gis.stackexchange.com/questions/332619/r-raster-to-shapefile

标签： r r-raster

【解决方案1】：

你可以的

r <- rasterize(x, r, field=1)

要创建一个新的栅格，这里不需要，你可以这样做

y <- raster(x)

【讨论】：