使用terra库转换不规则XYZ文件为栅格时Z值丢失问题排查

1. Z值丢失的原因及修复 #

你使用terra::rasterize时仅得到值1，是因为该函数默认会将矢量要素的存在标记为1，并未指定要提取的属性字段。你的SpatVector中保留了z列属性，只需在rasterize中添加field="z"参数，即可将原始Z值写入栅格：

final_raster <- terra::rasterize(v_etd_sf, r1, field = "z")

2. 合理设置ncols和nrows的方法 #

直接指定ncols=700和nrows=700属于主观猜测，更科学的方式是基于目标分辨率和数据空间范围自动计算行列数：

步骤： #

• 获取矢量数据的空间范围（extent）
• 确定适合数据精度的目标分辨率（例如WGS84坐标系下的0.0001度，需根据你的数据调整）
• 基于范围和分辨率自动生成模板栅格

代码示例：

# 获取矢量数据的空间范围
ext <- terra::ext(v_etd_sf)
# 设置目标分辨率（单位：度）
target_res <- 0.0001
# 创建匹配范围与分辨率的模板栅格
r1 <- terra::rast(ext, res = target_res, crs = terra::crs(v_etd_sf))

完整修正后的代码 #

# Purpose: convert from an irregular cell sizes xyz file to raster
library(terra)
library(data.table)

# 加载XYZ文件
etd <- fread("./exemplo.xyz")
data.table::setnames(etd, 1:3, c("x", "y", "z")) 

# 直接转换为SpatVector（无需中转sf）
v_etd <- terra::vect(etd, geom = c("x", "y"), crs = 4326)

# 创建模板栅格（基于目标分辨率）
ext <- terra::ext(v_etd)
target_res <- 0.0001  # 根据数据精度调整此值
r_template <- terra::rast(ext, res = target_res, crs = terra::crs(v_etd))

# 栅格化并保留Z值
final_raster <- terra::rasterize(v_etd, r_template, field = "z")

# 验证结果
plot(final_raster)
range(final_raster$z, na.rm = TRUE)  # 显示Z值的实际范围

额外说明 #

• terra::rast(type="xyz")仅支持规则网格的XYZ数据（x和y的间隔完全均匀），你的数据是不规则分布的点，因此该方法报错属于正常情况。
• 如果点数据存在空白区域需要填充，可考虑使用插值方法（如terra::interpolate结合IDW或克里金算法），具体需根据数据特性选择合适的插值方式。

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