使用R的raster包加载WorldClim栅格为何出现大量NA值？

问题根源 #

你遇到的大量NA值主要来自两个方面：

• WorldClim全球栅格数据本身包含大量无数据区域（比如海洋、非目标大陆），直接调用getValues()会返回所有栅格单元，自然会看到大量NA
• 当前先掩膜再裁剪的处理顺序效率低下，且无法彻底移除研究区外的无效NA值

解决方案步骤 #

1. 统一投影坐标系 #

首先确保研究区矢量数据与栅格数据的投影一致，否则mask/crop操作会失效：

library(raster)
library(sf)
library(rnaturalearth)

# 加载栅格并查看投影
temp <- raster("wc2.1_10m_bio_1.tif")
print(crs(temp))

# 加载研究区矢量并转换投影
Asia <- ne_countries(scale="medium", type="map_units", returnclass="sf", country=c('Myanmar','Thailand','Laos','Cambodia','Vietnam'))
study_area <- select(Asia, geometry, name_long)
study_area <- st_transform(study_area, crs = crs(temp))

2. 先裁剪再掩膜 #

先裁剪到研究区的大致范围，减少后续计算的数据量，再用矢量边界掩膜得到精确的研究区栅格：

# 定义研究区范围
ex <- extent(90, 112, 4, 30)

# 先裁剪栅格到研究区范围
temp_crop <- crop(temp, ex)
preci_crop <- crop(raster("wc2.1_10m_bio_12.tif"), ex)
elev_crop <- crop(raster("wc2.1_10m_elev.tif"), ex)

# 用研究区矢量进行掩膜
temp_masked <- mask(temp_crop, study_area)
preci_masked <- mask(preci_crop, study_area)
elev_masked <- mask(elev_crop, study_area)

3. 构建干净栅格栈并验证 #

# 构建栅格栈
envi_stack <- stack(elev_masked, temp_masked, preci_masked)
names(envi_stack) <- c("elev", "temp", "preci")

# 查看各图层非NA值数量（验证处理效果）
cellStats(envi_stack, function(x) sum(!is.na(x)))

# 提取无NA的环境变量值（用于SDM模型输入）
env_values <- na.omit(getValues(envi_stack))

4. 可选：用terra包提升处理效率 #

10m分辨率的WorldClim数据体积庞大，raster包处理速度慢且易内存溢出，推荐用terra包替代：

library(terra)

# 加载栅格
temp <- rast("wc2.1_10m_bio_1.tif")
preci <- rast("wc2.1_10m_bio_12.tif")
elev <- rast("wc2.1_10m_elev.tif")

# 转换矢量格式并统一投影
study_area <- vect(study_area)
study_area <- project(study_area, crs(temp))

# 批量裁剪+掩膜
envi_stack <- c(elev, temp, preci) |>
  crop(ex) |>
  mask(study_area)

names(envi_stack) <- c("elev", "temp", "preci")

# 提取无NA值
env_values <- na.omit(values(envi_stack))

关键注意事项 #

• WorldClim 10m分辨率单文件体积超1GB，处理前确保有16GB以上可用内存
• SDM模型通常只需要有效区域的环境变量数值，用na.omit()提取即可，无需保留全栅格的NA值
• 处理后用plot(envi_stack)可视化，确认有效区域是否覆盖目标研究区

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