使用 R 从 .asc 文件（ESRI ascii 网格）中删除标题（6 行）并导出

Question

我有超过 800 个 .asc 文件（ESRI ascii 网格），每个文件都有一个由 6 行组成的标题，然后是由空格分隔的栅格数据。这里以一个小文件为例。我使用 read.asciigrid（sp 包）阅读它。

new("SpatialGridDataFrame" , 数据 = 结构(list(mydata.asc = c(4, 4, 4, 4, 3, 4, 4, 4, 1, 1, 1, 1, 1, 4, 4, 4, 4, 3, 4 , 4, 4, 1, 1, 1, 1, 1, 4, 4, 4, 4, 3, 4, 4, 4, 1, 1, 1, 1, 1, 4, 4, 4, 4, 3, 4, 4, 4, 6, 1, 1, 1, 1, 4, 4, 4, 4, 3, 4, 4, 4, 6, 1, 1, 1, 1, 4, 4, 4, 4, 3, 4, 4, 4, 6, 1, 1, 1, 1, 4, 4, 4, 4, 4, 4, 4, 4, 6, 6, 1, 1, 1, 4, 4, 4, 4, 4, 4, 4, 4, 6, 6, 1, 1, 1, 4, 4, 4, 4, 4, 4, 4, 4, 4, 1, 1, 1, 1, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 6, 6, 6, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 6)), .Names = "mydata.asc", row.names = c(NA, -143L), class= "data.frame") , grid = new("GridTopology" , cellcenter.offset = c(394984.42630274, 2671265.4912109) , 像元大小 = c(25, 25) , cells.dim = c(13L, 11L) ) , bbox = 结构 (c(394971.92630274, 2671252.9912109, 395296.92630274, 2671527.9912109), .Dim = c(2L, 2L), .Dimnames = list(NULL, c("min", "max"))) , proj4string = new("CRS" , projargs = NA_character_ ) )

这是使用文本编辑器查看文件时的样子。

这是我想做的步骤

1) 读入文件 2）删除前6行（标题） 3) 将文件另存为具有相同文件名但位于不同位置的 .asc 文件

当然，我想对 800 个文件执行此操作，但如果我能弄清楚如何对一个文件执行此操作，我应该能够编写一个函数来遍历所有文件。

感谢您的帮助。

-al

更新： 感谢@Luca Braglia，这是对我有用的最终代码。

设置工作目录

setwd("c:/temp/hdr/ascii")
newdir <- "c:/temp/hdr/ascii_no_hdr/"

files <- dir(pattern="*.asc")

for (my.file in files){
  i <- read.table(my.file,skip=6,sep="")
  write.table(i,file=paste(newdir,my.file,sep=""),sep="",row.names=FALSE,col.names=FALSE)
}

我不想要 col 和 row 的名称。一段非常简单有效的代码。

Answer 1

您可以列出所有文件，在 for 循环中读取所有文件（使用 read.table 的 skip 选项）

## you are in the directory with your asc files
files <- dir(pattern="*.asc")

# loop
for (my.file in  files) {
     i <- read.table(my.file, skip = 6, sep = " ")
     # change names here if you don't want V1, V2 ...
     write.table(i, file = paste("new_dir", my.file, sep = "/"), 
                 sep = " ", row.names = FALSE)

}