1

我目前正在尝试从 CHIRPS 数据库中提取墨西哥的降雨数据。目标是获得 15 年期间每月降雨量的综合数据库。这涉及在从包含特定月份天气状况信息的 .tif 文件中提取数据后合并许多列。
CVE_ENT 和 CVE_MUN 是后来帮助我识别各个城市的两个变量。

在 R 中运行我的代码并查看生成的数据框，一切看起来都很好。但是，一旦我尝试将其提取为 .dta 或 .csv 文件，我就会收到以下错误消息： fwrite(vextractall, file="rainfall55.csv") fwrite(vextractall, file = "rainfall55. csv") ：第 4 列的长度 (1) 与第 1 列的长度 (2456) 不同

该错误发生在多种其他情况下，例如，如果使用 write.dta。

任何可能知道我错过了什么的人？提前非常感谢。

#Data get methods#
tif.raster1 <- raster('chirps-v2.0.1999.01.tif')
crs.LL <- "+proj=longlat +datum=WGS84 +no_defs +ellps=WGS84" 
tif.raster1 <- projectRaster(tif.raster1, crs = crs.LL)
mexico2shp <- readOGR(dsn='GIS Mexican Municipalities', layer='Mexican      Municipalities')
tif.raster1 <- crop(tif.raster1, extent(mexico2shp))
vras.tif1 <- velox(tif.raster1)
iters <- nrow(mexico2shp)
#mapping function#
vextractall <- vras.tif1$extract(mexico2shp, fun=mean)
mexicomm <- as.data.frame(mexico2shp)
vextractall <- as.data.frame(vextractall)
iters <- nrow(mexico2shp)
x <- foreach(a=1:iters) %do% {
if(is.na(vextractall[a,1])) {
 ext <- raster::extract(tif.raster1, mexico2shp[a,], fun=mean)
 vextractall[a,1] <- ext[1,1]
 }
}
vextractall<-as.data.frame(vextractall)
vextractall$CVE_ENT <- mexicomm[,c("CVE_ENT")]
vextractall$CVE_MUN <- mexicomm[,c("CVE_MUN")]
vextractall<-vextractall[,c(ncol(vextractall), 1:(ncol(vextractall)-1))]
vextractall<-vextractall[,c(ncol(vextractall), 1:(ncol(vextractall)-1))]
vextractall <- plyr::rename(vextractall, c("V1"="Milimeters011999"))

tif.raster1 <- raster('chirps-v2.0.1999.02.tif')
tif.raster1 <- crop(tif.raster1, extent(mexico2shp))
vras.tif1 <- velox(tif.raster1)
crs.LL <- "+proj=longlat +datum=WGS84 +no_defs +ellps=WGS84" 
tif.raster1 <- projectRaster(tif.raster1, crs = crs.LL)
vextract2 <- vras.tif1$extract(mexico2shp, fun=mean)
vextract2 <- as.data.frame(vextract2)
iters <- nrow(mexico2shp)
foreach(a=1:iters) %do% {
if(is.na(vextract2[a,1])) {
ext <- raster::extract(tif.raster1, mexico2shp[a,], fun=mean)
vextract2[a,1] <- ext[1,1]
}
 }
vextractall<-as.data.frame(vextractall)
vextract2<-as.data.frame(vextract2)
vextractall$Milimeters021999 <- vextract2
4

0 回答 0