使用terra::rasterize栅格化甲虫观测数据框时全量数据无值问题

问题背景 #

使用terra包栅格化10条甲虫观测数据子集时，结果正常且统计正确；但处理18万条全量数据时，生成的栅格无有效值（全为NA），同时出现警告：

In class(object) <- "environment" :
  Setting class(x) to "environment" sets attribute to NULL; result will no longer be an S4 object

排查与解决方案 #

1. 检查坐标与分类列的有效性

全量数据可能存在异常值或缺失值，导致栅格化失败：

# 检查x/y坐标的缺失与异常值
sum(is.na(dist_data$x))
sum(is.na(dist_data$y))
summary(dist_data[,c("x","y")])

# 检查Taxon列的缺失
sum(is.na(dist_data$Taxon))

# 过滤异常数据
dist_data_clean <- dist_data[!is.na(dist_data$x) & !is.na(dist_data$y) & !is.na(dist_data$Taxon) & is.finite(dist_data$x) & is.finite(dist_data$y),]

2. 跳过sf中转，直接创建terra矢量对象

警告信息提示S4对象转换异常，可直接用terra::vect()从数据框创建矢量，避免sf转换环节：

# 直接从数据框创建矢量对象
v <- vect(dist_data_clean, geom = c("x", "y"), crs = 4326)

# 创建匹配范围的栅格
e <- ext(v)
r <- rast(e, res = 1, crs = 4326)

# 执行栅格化
new_ras <- rasterize(v, r, "Taxon")
summary(new_ras)

3. 确认栅格范围与矢量范围匹配

若栅格范围未覆盖所有观测点，会导致所有点被判定为在栅格外，生成全NA栅格：

# 查看矢量与栅格的范围
ext(v)
ext(r)

# 手动指定栅格范围为矢量的范围
r <- rast(ext(v), res = 1, crs = 4326)

4. 分批次栅格化大数据量

18万条数据量较大，分物种批次处理可避免内存或处理异常：

# 按物种分组
sp_groups <- split(dist_data_clean, dist_data_clean$Taxon)

# 初始化空栅格
r <- rast(ext(dist_data_clean), res = 1, crs = 4326)
final_rast <- init(r, NA)

# 循环处理每个物种
for(sp_name in names(sp_groups)){
  sp_points <- vect(sp_groups[[sp_name]], geom = c("x", "y"), crs = 4326)
  sp_rast <- rasterize(sp_points, r, field = rep(sp_name, nrow(sp_groups[[sp_name]])))
  final_rast <- cover(final_rast, sp_rast)
}

summary(final_rast)

5. 更新依赖包版本

警告信息可能源于terra或sf包的版本兼容问题，更新至最新版：

install.packages(c("terra", "sf"), dependencies = TRUE)

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