從 R 中的 shapefile 創建柵格

Question

我有一個研究區域的 shapefile，我想創建一個柵格，其中屬於該 shapefile 的每個單元格的 value=1，其他單元格的 value=NA。

我想創建一個與我將在分析中使用的環境柵格具有相同原點和像元大小的新柵格。

這是我正在使用的代碼，但我總是收到此錯誤：“rep(value, length.out = 2) 中的錯誤：嘗試復制‘S4’類型的 object”

#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)

（我還沒有弄清楚如何將值分配給柵格）

有人對我為什么會收到此錯誤有任何建議嗎？ 關於如何分配我想要的值？

Answer 1

你可以做

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

要創建此處不需要的新柵格，您可以執行

y <- raster(x)