R Studio中提取TIFF影像矢量框值后NA无法去除，均值计算报错求助

问题重现 #

尝试从TIFF地图文件的矢量框提取数值并计算均值，但提取结果始终带NA，且计算均值时报错：

# 加载栅格
DEM <- raster("del15N_130_raster.tif")
DEM <- setMinMax(DEM)
plot(DEM)

# 加载矢量并提取数值
AntC <- read_sf("shp/Antartica_typeC.shp")
plot(AntC, add = T)
vals <- terra::extract(DEM, AntC, na.rm=TRUE )
vals_mean <- mean(vals$layer)

报错信息：

Warning message:
In mean.default(vals$layer) :
argument is not numeric or logical: returning NA

尝试过na.omit(DEM)、mean(vals$layer, na.rm=TRUE)等方法，问题依旧。

解决方案 #

1. 统一使用terra包处理栅格 #

避免raster与terra包混合使用导致的兼容问题，直接用terra加载栅格：

# 加载为SpatRaster对象（terra的标准栅格格式）
DEM <- terra::rast("del15N_130_raster.tif")
plot(DEM)

2. 检查并匹配坐标系 #

坐标系不匹配是提取全NA的核心原因之一，先验证两者坐标系是否一致：

# 查看矢量坐标系
st_crs(AntC)
# 查看栅格坐标系
terra::crs(DEM)

若不一致，将矢量投影转换为栅格的坐标系：

AntC <- sf::st_transform(AntC, terra::crs(DEM))

3. 直接用extract内置函数计算均值 #

无需手动提取所有值再计算，直接在terra::extract中指定fun=mean并开启na.rm=TRUE：

# 一次性提取每个矢量框的均值，自动忽略NA
mean_vals <- terra::extract(DEM, AntC, fun = mean, na.rm = TRUE)
# 查看结果，每行对应一个矢量要素的均值
print(mean_vals)

4. 排查剩余NA的可能原因 #

若仍得到NA，检查以下两点：

• 矢量框是否完全在栅格范围外：用terra::intersect(DEM, AntC)验证重叠区域
• 栅格对应区域是否本身全为NA：用terra::extract(DEM, AntC)查看原始提取值，确认是否存在有效数值

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