R栅格分析（坡度计算与for循环）缺失值报错问题求解

你之前新增的判断仅覆盖了treefile无任何数据的场景，即便瓦片有有效数据，仍可能存在部分点位的坡度值为NA：比如点位落在DEM空值区域、DEM边缘导致坡度计算返回空值，此时循环中执行slope <= 20判断时，R遇到NA值就会触发缺失值报错。

先执行以下代码确认坡度列是否存在空值：

sum(is.na(treefile@data$slope))

如果返回结果大于0，即可确认是该问题导致的报错。

方案1：修改原有循环逻辑 #

在坡度判断前增加空值校验，你可以根据业务需求补充坡度为空时的处理规则：

if( nrow(treefile@data)==0 ){
  out1 = r.crop
  out1[out1] = NA
  out1 = stack(out1,out1,out1)
  return(out1)
} else {
  for( i in 1:nrow(treefile@data)){
    # 增加坡度非空判断
    if( treefile@data[i,"species"] == "NS" && !is.na(treefile@data[i,"slope"])){
      if( treefile@data[i,"slope"] <= 20 ){
        treefile@data[i,"species"] = "U_NS_f"
      } else {
        treefile@data[i,"species"] = "U_NS_s"
      }
    }
    # 可选：补充坡度为空时的NS树种处理逻辑，例如下方默认赋值为U_NS_f
    # else if(treefile@data[i,"species"] == "NS" && is.na(treefile@data[i,"slope"])){
    #   treefile@data[i,"species"] = "U_NS_f"
    # }
  }
}

方案2：使用向量化操作替代循环（更高效，适合大数据量） #

R中针对数据框的向量化操作远快于for循环，建议用该写法替换原有循环：

if( nrow(treefile@data)==0 ){
  out1 = r.crop
  out1[out1] = NA
  out1 = stack(out1,out1,out1)
  return(out1)
} else {
  # 筛选出树种为NS且坡度不为空的点位
  ns_valid_mask <- treefile@data$species == "NS" & !is.na(treefile@data$slope)
  # 按坡度赋值
  treefile@data[ns_valid_mask & treefile@data$slope <=20, "species"] <- "U_NS_f"
  treefile@data[ns_valid_mask & treefile@data$slope >20, "species"] <- "U_NS_s"
  
  # 可选：处理坡度为空的NS树种
  # treefile@data[treefile@data$species == "NS" & is.na(treefile@data$slope), "species"] <- "U_NS_f"
}

如果坡度为空的点位属于无效数据，你可以在提取坡度后直接过滤掉这些点位，从根源避免后续空值问题：

treefile@data$slope = round(raster::extract(r.temp, treefile))
# 过滤坡度为空的点位
treefile <- treefile[!is.na(treefile@data$slope), ]

内容的提问来源于stack exchange，提问作者Max