R语言使用fread读取ASCII栅格文件如何保留原始格式？

你遇到的问题根源是ESRI ASCII Grid格式文件本身的结构特殊性：前6行为元数据行，每行仅2个字段（参数名+参数值），从第7行开始才是与ncols列数一致的栅格数值行。你直接用fread统一按分隔符读取全文件时，会自动按最大列数（即栅格数据的列数）拆分所有行，导致前6行产生大量空值NA，后续导出自然不符合标准ASC格式要求。

方法1：使用专业空间栅格包读取（推荐） #

直接用terra包的rast()函数读取ASC文件，原生支持格式解析，读取后直接是标准栅格对象，后续批量处理、导出都不会出现格式错乱问题。

library(purrr)
library(terra)

# 批量读取所有asc文件
files <- list.files(pattern = "*.asc")
raster_list <- map(files, rast)

# 后续如果要重命名导出，直接用writeRaster即可，自动保留标准ASC格式
# 示例：导出第一个文件
writeRaster(raster_list[[1]], "新文件名.asc", overwrite=TRUE)

方法2：手动用fread分块读取（如需自定义处理逻辑） #

如果必须手动处理文本内容，可拆分读取元数据和栅格数值两部分：

library(data.table)
library(purrr)

read_asc <- function(file_path) {
  # 读取前6行元数据
  meta <- readLines(file_path, n=6)
  # 读取第7行开始的栅格数值
  grid_vals <- fread(file_path, skip=6, header=FALSE, sep=" ")
  # 导出时先写元数据，再写栅格数值即可还原格式
  output_name <- paste0("重命名_", basename(file_path))
  writeLines(meta, output_name)
  fwrite(grid_vals, output_name, append=TRUE, sep=" ", col.names=FALSE)
  return(TRUE)
}

# 批量处理所有文件
files <- list.files(pattern = "*.asc")
walk(files, read_asc)

• 不要在读取时用fill=TRUE参数，该参数会强制补全行列，是导致前6行产生NA的直接原因之一
• 导出ASC文件时不要输出列名行，也不要额外添加索引列，否则会破坏格式

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