如何检查R语言SpatRaster中是否存在NA并返回布尔值

在terra包中，anyNA()针对SpatRaster的实现是逐像元返回逻辑值，因此会生成一个与原栅格维度一致的SpatRaster对象，而非单一的TRUE/FALSE结果。要判断整个栅格是否存在NA，可以用以下两种高效方法：

方法1：全局统计（适合大栅格，无需加载全部数据到内存） #

使用global()函数结合自定义统计逻辑，直接计算栅格中NA的总数，再判断是否大于0：

library(terra)
r <- rast("~/Downloads/CHELSA_bio10_1981-2010_V.2.1.tif")

# 统计栅格中NA的数量
na_total <- global(r, function(x) sum(is.na(x)), na.rm = FALSE)

# 判断是否存在NA
has_missing_values <- na_total > 0
has_missing_values

方法2：加载像元值检查（适合小栅格） #

如果栅格数据量较小，可以直接提取所有像元值，再用any(is.na())判断：

# 提取所有像元值并检查是否存在NA
has_missing_values <- any(is.na(values(r)))
has_missing_values

补充说明 #

你之前得到的anyNA(r)返回全FALSE的SpatRaster，说明当前预览的栅格区域没有NA，但如果要确认整个栅格的情况，必须用全局统计或全像元检查的方法——因为anyNA()的逐像元结果不会自动汇总为单一逻辑值。

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