R循环处理多国空间多边形报replacement has length zero错误求助

• 空栅格读取异常：部分国家矢量范围与population/wind栅格无重叠，或是经阈值过滤后区域内所有栅格值被设为NA，后续mask、cellStats操作读取空栅格时触发该报错。
• 无容错机制：循环未做异常捕获，单个国家处理失败就会终止全部流程，无法继续处理剩余国家。
• 延迟读取冲突：raster包默认采用延迟读取逻辑，栅格数据仅在需要时读取本地临时文件，循环运行时间过长导致临时文件被清理、文件连接溢出时也会触发读取失败。

# 预处理：如果栅格体积不大，先全量加载到内存避免读取冲突
population <- readAll(population)
wind <- readAll(wind)

bound <- list(ARG, AUS, AUT, BEL, BGR, BRA, CAN, CHE, CHL, CHN, COL, CRI, CZE, DEU, DNK, ESP, EST, EU27, FIN, FRA, GBR, GRC, HRV, HUN, IDN, IND, IRL, ISL, ISR, ITA, JPN, KOR, LTU, LUX, LVA, MEX, NLD, NOR, NZL, PER, POL, PRT, ROU, RUS, SAU, SVK, SVN, SWE, TUR, USA, ZAF)
result <- data.frame() 

for (i in 1:length(bound)) { 
  country <- bound[[i]]
  iso_code <- country@data$iso3
  cat("正在处理：", iso_code, "\n") # 打印当前处理国家，方便定位出错项
  
  # 加异常捕获，出错时跳过当前国家继续运行
  tryCatch({
    clip1 <- crop(population, extent(country))
    clip2 <- crop(wind, extent(country))
    cat1 <- calc(clip2, fun=function(x){ x[x > 1.1] <- NA; return(x)})
    clip4 <- mask(clip1, cat1)
    
    # 先判断裁剪后栅格是否有有效值，避免空计算
    if (cellStats(!is.na(clip4), sum) == 0) {
      cat(iso_code, "区域无有效值，跳过\n")
      # 可根据需求生成空记录占位
      final <- data.frame(iso3 = iso_code, value = NA, Freq = 0)
    } else {
      ext_cat1 <- cellStats(clip4, 'sum')
      tab <- lapply(ext_cat1, table)
      mat <- as.data.frame(tab)
      final <- cbind(iso3 = iso_code, mat)
    }
    result <- rbind(result, final)
  }, error = function(e) {
    # 打印错误信息但不中断循环
    cat(iso_code, "处理失败，错误信息：", e$message, "\n")
    # 可选择生成错误占位记录
    error_row <- data.frame(iso3 = iso_code, value = NA, Freq = NA)
    result <<- rbind(result, error_row)
  })
}

• 优先使用terra包替代raster包：terra是raster开发者推出的新一代栅格处理包，读取性能和稳定性远高于raster，相同逻辑代码运行速度可提升3-10倍。
• 避免在循环内用rbind累积结果：可以先把每个国家的结果存在列表里，循环结束后一次性合并，大幅提升运行效率。

内容的提问来源于stack exchange，提问作者Mikaël