使用raster包rasterFromXYZ转换XYZ数据报错x单元格大小不规则如何解决

Error in rasterFromXYZ(DATA) : x cell sizes are not regular

该报错的本质是rasterFromXYZ()要求输入的XYZ数据对应规则格网，即x坐标的相邻间隔、y坐标的相邻间隔都必须完全一致。90%以上的报错场景是浮点精度误差导致的坐标微小偏差，而非数据本身的格网不规则。

• 方案1：调整digits参数放宽精度阈值
  绝大多数场景下通过调整坐标匹配的精度参数即可解决问题，digits用来控制坐标保留的小数位数，过滤浮点运算带来的微小偏差，可根据你的坐标精度调整数值：

  library(raster)
  # 示例：坐标精确到小数点后4位就设置digits为4，精度更高可调整为5或6
  rst <- rasterFromXYZ(DATA, digits = 4)
  

• 方案2：手动构造规则栅格模板强制转换
  如果调整digits无效，可以手动计算理论分辨率，通过rasterize函数完成转换，不需要依赖函数自动识别坐标间隔：

  library(raster)
  # 提取坐标和值列，默认第一列为x、第二列为y、第三列为属性值
  x <- DATA[,1]
  y <- DATA[,2]
  z <- DATA[,3]
  # 取出现频次最高的间隔作为规则格网的分辨率
  res_x <- median(diff(sort(unique(x))))
  res_y <- median(diff(sort(unique(y))))
  # 生成空规则栅格模板
  rst_template <- raster(
    xmn = min(x) - res_x/2,
    xmx = max(x) + res_x/2,
    ymn = min(y) - res_y/2,
    ymx = max(y) + res_y/2,
    res = c(res_x, res_y)
  )
  # 将XYZ值映射到栅格模板上
  rst <- rasterize(DATA[,1:2], rst_template, field = DATA[,3])
  

• 方案3：排查原始数据的坐标异常
  如果以上两种方案都无效，说明原始数据确实存在非规则坐标点，需要先做预处理：
  1. 提取所有唯一x坐标，排序后计算相邻差值，定位间隔异常的坐标点
  2. 修正异常坐标、剔除重复/错位的无效点后，再重新调用转换函数
     排查代码示例：

  x_unique <- sort(unique(DATA[,1]))
  # 输出所有x坐标的相邻间隔，定位异常值
  print(diff(x_unique))
  

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