我有一個26層的光柵堆棧,代表了一年中不同的星期。在這個區域內,有854個不規則形狀的多邊形,我對每個層的(mean和sd)的彙總值感興趣。彙總跨多邊形的值,對於柵格堆棧中的每個圖層 - 如何將值與shapefile合併?
我相信我已經成功地提取了這些值。但是,當我嘗試將它們與多邊形文件合併時,我不成功。我相信extract()
函數已經剝奪了它們的標識條件(Name_2)的值。我嘗試將extract()
函數的「大型列表」輸出轉換爲數據框,但這也不起作用。
有什麼建議嗎?
謝謝。
Brazil0 <- getData('GADM', country='BRA', level=0)
Brazil1 <- getData('GADM', country='BRA', level=1)
Brazil2 <- getData('GADM', country='BRA', level=2)
MinasGerais<-subset(Brazil2, NAME_1 =='Minas Gerais')
e <- extent(MinasGerais)
MG <- as(Brazil2[Brazil2$NAME_1 == 'Minas Gerais',], 'SpatialPolygons')
row.names(MG) <- as.character(1:length(MG))
setwd("H:/Brazilian Arbovirus Project/ratsers_EpiWeek")
Precip_2017 <- stack("Precip_2017.tif")
Precip_crop <- crop(Precip_2017, e, snap="out")
crop<- setValues(Precip_crop, NA)
Precip_raster <- rasterize(MinasGerais, crop)
Precip_MG<- mask(x=Precip_crop, mask=Precip_raster)
MG_Municip_mean_2 <- extract(Precip_MG,
MG,
method='simple',
match.ID=FALSE,
FUN=mean,
#sp=TRUE,
small=TRUE)
output=data.frame(MG_Municip_mean_2)
MG_1 <- merge(MinasGerais, MG_Municip_mean_2, by='NAME_2', all=TRUE)
writeOGR(MG_1, getwd(), "MG_Municipalities", driver="ESRI Shapefile",
check_exists=TRUE,
overwrite_layer=TRUE)
當你說你「相信你已經提取這些值」,有你在調試代碼證實了這一點,還是你只是假設它的工作? – FluffyKitten
您可以使用dput(Precip_2017)創建類似的數據來添加一些示例數據嗎?你使用哪些軟件包? – minem