问题描述
我在使用 RasterStack 循环时遇到了一点问题。我试图提取过去 10 年中北卡罗来纳州每个县每小时的平均温度。我能够运行一些简单的代码并相对容易地为每一层获取它,但是当我试图循环这个过程以获得光栅堆栈/砖块中每个时间戳的平均值时,我没有太多运气!如果您对此有任何帮助或建议,我将不胜感激。
library(tidyverse)
library(rgdal)
library(sp)
library(raster)
library(doParallel)
library(foreach)
library(ncdf4)
librarcy(exactextractr)
# Brick the NC File
r<-stack("ERA_Temps.nc")
output <- for (i in 1:length(r)) {
wgsraster <- r[[i]]
year <- str_sub(names(wgsraster),2,5)
month <- str_sub(names(wgsraster),7,8)
day <- str_sub(names(wgsraster),10,11)
hour <- str_sub(names(wgsraster),13,14)
NC_Counties <- cbind(shp,exact_extract(wgsraster,shp,c('mean')))
output_df <- data.frame(County.Name = NC_Counties$NAME,GEOID = NC_Counties$GEOID,temp_mean = NC_Counties$mean,year,month,day,hour)
return(output_df)
}
Output <- output_df
示例和更多信息:
r
class : RasterStack
dimensions : 81,81,6561,72 (nrow,ncol,ncell,nlayers)
resolution : 0.25,0.25 (x,y)
extent : -90.125,-69.875,24.875,45.125 (xmin,xmax,ymin,ymax)
crs : +proj=longlat +datum=wgs84 +no_defs
names : X2018.09.14.00.00.00,X2018.09.14.01.00.00,X2018.09.14.02.00.00,X2018.09.14.03.00.00,X2018.09.14.04.00.00,X2018.09.14.05.00.00,X2018.09.14.06.00.00,X2018.09.14.07.00.00,X2018.09.14.08.00.00,X2018.09.14.09.00.00,X2018.09.14.10.00.00,X2018.09.14.11.00.00,X2018.09.14.12.00.00,X2018.09.14.13.00.00,X2018.09.14.14.00.00,...
我可以很容易地获得单个图层所需的内容:
wgsraster <- r[[5]]
year <- str_sub(names(wgsraster),5)
month <- str_sub(names(wgsraster),8)
day <- str_sub(names(wgsraster),11)
hour <- str_sub(names(wgsraster),14)
NC_Counties <- cbind(shp,c('mean')))
output_df <- data.frame(County.Name = NC_Counties$NAME,temp_mean = NC_Counties$exact_extract.wgsraster..shp..c..mean...,hour)
output_df
County.Name GEOID temp_mean year month day hour
1 Saluda 45081 298.3141 2018 09 14 11
2 Williamsburg 45089 300.2301 2018 09 14 11
3 Dorchester 45035 299.3716 2018 09 14 11
4 Berkeley 45015 299.7389 2018 09 14 11
5 Greenwood 45047 298.5175 2018 09 14 11
但是当我尝试运行上面提到的循环时,它要么失败并且没有输出,要么给出以下错误:
Error in .local(x,i,j,...) : you must provide an index
解决方法
暂无找到可以解决该程序问题的有效方法,小编努力寻找整理中!
如果你已经找到好的解决方法,欢迎将解决方案带上本链接一起发送给小编。
小编邮箱:dio#foxmail.com (将#修改为@)