将SpatRaster转换为DataFrame时出现行数不匹配错误的原因

问题背景 #

栅格对象信息 #

library(terra)
    
r
    
class       : SpatRaster 
dimensions  : 6000, 6000, 1  (nrow, ncol, nlyr)
resolution  : 0.0008333333, 0.0008333333  (x, y)
extent      : 89.99958, 94.99958, 4.999583, 9.999583  (xmin, xmax, ymin, ymax)
coord. ref. : lon/lat WGS 84 (EPSG:4326) 
source      : memory 
name        : SSP_2050 
min value   : 0.0100 
max value   : 0.2376 

执行的转换代码 #

r_df <- terra::as.data.frame(r, xy = T, na.rm = F)

报错信息 #

Error in data.frame(..., check.names = FALSE) : 
  arguments imply differing number of rows: 36000000, 43596000

错误原因 #

这个错误的核心是生成的坐标列（x、y）与栅格值列的行数不匹配：

• 从栅格维度计算，6000×6000的栅格应对应36000000行数据，但报错显示其中一组数据的行数为43596000，说明SpatRaster对象内部数据结构存在异常。
• 大概率是栅格在创建、修改或内存存储过程中出现了索引损坏，导致terra::as.data.frame在生成坐标序列时计算出错误的行数。
• 另外，虽表面上范围、分辨率与行列数的数值匹配，但内部存储的范围参数可能存在微小精度误差，触发了错误的行数计算；如果栅格是从其他格式导入后修改过范围/分辨率，也可能导致内部数据结构不一致。

验证方法 #

可单独提取坐标和栅格值，查看两者行数是否一致：

# 提取坐标
coords <- terra::xyFromCell(r, 1:terra::ncell(r))
nrow(coords)
# 提取栅格值
values <- terra::values(r)
nrow(values)

若行数不一致，说明栅格内部结构损坏，建议重新创建或修复该SpatRaster对象。

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