如何从SpatialPixelsDataFrame中去除含NA值的像素？

下面提供两种无需依赖过时函数的方法，最终都能返回SpatialPixelsDataFrame类型对象：

方法1：直接在sp对象上处理（推荐） #

利用complete.cases()函数筛选出@data槽中无NA值的行，直接对原SpatialPixelsDataFrame进行子集化操作，无需转换为sf类：

library(sp)
library(raster)

# 生成示例数据（和你的代码一致）
m <- raster::raster(ncol=100, nrow=100)
m[] <- runif(ncell(m))
m[m < 0.2] <- NA
spdf <- as(m, "SpatialPixelsDataFrame")

# 筛选无NA的像素
spdf_clean <- spdf[complete.cases(spdf@data), ]

# 验证结果：检查@data中是否还有NA
any(is.na(spdf_clean@data)) # 应返回FALSE
# 验证类型
class(spdf_clean) # 应返回"SpatialPixelsDataFrame"

方法2：通过sf类中转（使用stats原生na.omit） #

如果需要通过sf类处理，可直接使用stats::na.omit（而非过时的sp.na.omit）处理sf对象，再转换回SpatialPixelsDataFrame：

library(sp)
library(sf)
library(raster)

# 生成示例数据
m <- raster::raster(ncol=100, nrow=100)
m[] <- runif(ncell(m))
m[m < 0.2] <- NA
spdf <- as(m, "SpatialPixelsDataFrame")

# 转sf -> 移除NA -> 转回SpatialPixelsDataFrame
sf_obj <- st_as_sf(spdf)
sf_clean <- na.omit(sf_obj)
spdf_clean <- as(sf_clean, "SpatialPixelsDataFrame")

# 验证结果
any(is.na(spdf_clean@data)) # FALSE
class(spdf_clean) # "SpatialPixelsDataFrame"

关键说明 #

• complete.cases()会逐行检查@data中的所有列，仅保留无任何NA值的行；
• 方法1无需额外依赖sf包，操作更直接高效；
• 方法2遵循警告提示，使用原生na.omit处理sf对象，避免过时函数的警告。

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