使用R的terra包加载asc.gz文件创建栅格栈时坐标读取异常

你的代码出现坐标读取问题，核心原因是手动读取ASCII文件并转成矩阵时丢失了栅格的头文件信息。ASC格式的栅格文件开头包含ncols、nrows、xllcorner、yllcorner、cellsize等关键空间参数，用read.table读取时会把这些参数当成数据行处理，导致生成的栅格没有正确的地理坐标、分辨率和范围，自然绘图显示异常。

正确处理方式 #

terra包的rast()函数支持直接读取压缩的.asc.gz文件，会自动解析头文件里的空间参数，不需要手动读表转矩阵。

方法1：直接批量读取生成栅格栈（推荐） #

这是最简洁高效的方式，一行代码就能完成：

# 替换为你的文件所在目录
filenames <- list.files(path = "./your_data_dir", pattern = "\\.asc\\.gz$", full.names = TRUE)
# 直接读取所有文件，自动生成栅格栈
my_raster <- rast(filenames)

方法2：循环逐个读取添加（适合需要分步处理的场景） #

如果必须用循环，也要直接用rast()读取文件，而不是手动解析：

filenames <- list.files(path = "./your_data_dir", pattern = "\\.asc\\.gz$", full.names = TRUE)
# 初始化空栅格栈
my_raster <- rast()

for (file_path in filenames) {
  # 直接读取压缩ASC文件，保留空间参数
  current_rast <- rast(file_path)
  # 添加到栅格栈
  my_raster <- c(my_raster, current_rast)
}

验证方法 #

读取完成后，可以用以下代码确认空间参数是否正确：

# 查看栅格范围
ext(my_raster)
# 查看分辨率
res(my_raster)
# 查看投影信息
crs(my_raster)

绘制单个图层时，坐标显示就会正常：

# 绘制第一个图层
plot(my_raster[[1]])

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