如何最快检查SpatRaster是否存在至少一个非NA值？

你当前的代码确实会遍历整个栅格计算全局非NA值数量，对大型栅格来说效率偏低。要实现找到第一个非NA值就停止检查的需求，推荐以下几种高效方案：

方法1：极简直接法（优先推荐） #

利用terra::cells的内置特性，指定只返回第一个非NA单元格的索引，内部会在找到目标后立即终止遍历，效率最高：

has_non_na <- function(x) {
  # 获取第一个非NA单元格的索引，存在则返回TRUE
  non_na_cell <- terra::cells(x, na.rm = TRUE, n = 1)
  length(non_na_cell) > 0
}

调用方式：has_non_na(your_spatraster)

方法2：分块迭代检查 #

如果需要自定义遍历粒度（比如适配内存限制），可以分块读取栅格行数据，每块检查完就判断是否存在非NA值，存在则立即返回：

has_non_na <- function(x) {
  block_size <- 10000  # 每次处理的行数，可根据内存调整
  total_rows <- terra::nrow(x)
  # 按块遍历栅格行
  for (i in seq(1, total_rows, block_size)) {
    current_end <- min(i + block_size - 1, total_rows)
    block_vals <- terra::values(x, row = i, nrows = current_end - i + 1)
    if (any(!is.na(block_vals))) {
      return(TRUE)
    }
  }
  return(FALSE)
}

方法3：采样+兜底检查 #

针对非NA值分布较均匀的场景，先通过随机采样快速验证，采样未命中时再逐块兜底检查，兼顾速度与准确性：

has_non_na <- function(x) {
  # 先随机采样100个单元格（数量可按需调整）
  sample_vals <- terra::spatSample(x, size = 100, na.rm = FALSE)
  if (any(!is.na(sample_vals))) {
    return(TRUE)
  }
  # 采样未命中时，逐块兜底检查
  block_size <- 10000
  total_rows <- terra::nrow(x)
  for (i in seq(1, total_rows, block_size)) {
    current_end <- min(i + block_size - 1, total_rows)
    block_vals <- terra::values(x, row = i, nrows = current_end - i + 1)
    if (any(!is.na(block_vals))) {
      return(TRUE)
    }
  }
  return(FALSE)
}

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