J'ai une pile raster avec 26 couches qui représentent différentes semaines de l'année. Dans cette zone il y a 854 polygones de forme irrégulière que je suis intéressé par les valeurs récapitulatives pour (moyenne et sd), pour chacune des couches. Je crois que j'ai réussi à extraire ces valeurs. Cependant, lorsque j'essaie de les fusionner avec le fichier polygone, je ne réussis pas. Je crois que la fonction extract()
a dépouillé les valeurs de leurs critères d'identification (Name_2). J'ai essayé de convertir la sortie "Large list" de la fonction extract()
en une trame de données, mais cela n'a pas fonctionné non plus.Résumer des valeurs à travers des polygones, pour chaque couche dans Raster Stack - Comment fusionner des valeurs avec un fichier de formes?
Des suggestions?
Merci.
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)
Quand vous dites que vous « croyez que vous avez extrait ces valeurs », avez-vous confirmé par le débogage du code ou êtes-vous en train de supposer que cela fonctionne? – FluffyKitten
pouvez-vous ajouter quelques exemples de données en utilisant dput (de Precip_2017) ou créer des données similaires? quels paquets utilisez-vous? – minem