use*_*029 3 foreach r doparallel terra
我以前从未遇到过此问题,但是当我尝试在 R 中使用 foreach 循环时遇到此错误:“{ 中的错误:任务 1 失败 - “NULL 值作为符号地址传递”。
我几乎不可能生成一个小的、可重复的示例(我已经尝试过!),因为我试图从巨大的栅格中提取数据并从该数据创建 csv 文件。但是,这是我的代码。
bi_2021 <- rast('G:\\GridMet_Yearly\\bi_2021.nc')
cl <- makeCluster(2)
registerDoParallel(cl)
r = 1
foreach (r=1:10, .packages = c('tidyverse','lubridate')) %dopar% {
rc <- row_char[r]
cc <- col_char[r]
ce <- cell_char[r]
rn <- row_num[r]
cn <- col_num[r]
fname <- paste0('G:/GridMet_Cells_RawData/row',rc,'_col',cc,'_cell',ce,'.csv')
data_df <- data.frame(read_csv(fname, show_col_types = FALSE)) # read previous data in
data_df <- data[which(year(data$Date) < 2021),]
# add rows for 2021 daily data
data_df[15342:15673,] <- NA
data_df$Date[15342:15673] <- seq(as.Date('2021-01-01'),as.Date('2021-11-28'),'days')
data_df$bi[15342:15673] <- as.numeric(bi_2021[rn,cn,][1:332]) # THIS IS THE LINE IT DOESN'T
LIKE!
write_csv(final_df,paste0('G://GridMet_Cells_RawData2//row',rc,'_col',cc,'_cell',ce,'.csv'))
rm(data_df,cc,ce,cn,fname,rc,rn)}
Run Code Online (Sandbox Code Playgroud)
这是数据示例。
Date bi erc etr_alfalfa fm100 fm1000 etr_grass pp rhmax rhmin shum srad ud tmin tmax vpd us
<date> <dbl> <dbl> <dbl> <dbl> <dbl> <dbl> <dbl> <dbl> <dbl> <dbl> <dbl> <dbl> <dbl> <dbl> <dbl> <dbl>
1 1979-01-01 22 16 0.1 18.2 24.5 0.1 0 74.4 38.7 0.14 78.1 232 -40.8 -26.8 20 2.5
Run Code Online (Sandbox Code Playgroud)
由于某种原因,它不允许我从栅格中提取数据,并将其放入每个变量的特定行中。关于为什么要这样做以及如何解决它有什么想法吗?当它没有放入 foreach 循环中时它工作得很好。谢谢。
您不能直接在并行化中使用 SpatRaster。请参阅此github 问题。某些方法内置了支持(predict、app和lapp)tapp,并且在其他情况下有不同的方法来解决此问题(也许请参阅此问题)。
例如,如果您使用此行,它可能会起作用
bi_2021 <- rast('G:\\GridMet_Yearly\\bi_2021.nc')
Run Code Online (Sandbox Code Playgroud)
循环内foreach。