0

我有从一个大列创建 141 x 131 矩阵并将其导出为文本文件的代码。到目前为止，它运作良好。我将使用 ASCII 到栅格工具将这些文件导入 ArcGIS。为了让 ArcGIS 导入它们，需要一个五行标题，标题后面的数据如下所示：

NCOLS 131
NROWS 141
XLLCENTER -150
YLLCENTER -70
CELLSIZE 1
NODATA_VALUE -32768
265 41.7 1.8 12 16 18.3 16.5 18.3 19.7 24.5 26.1 28.4...

一旦 R 保存文本文件，手动粘贴此标题很简单，但我需要处理 900 多个文件，因此我希望 R 在导出文本文件时添加标题。这是可能的还是我可以将标题行添加为矩阵的前 5 行然后导出它？

编辑：下面的这一部分是关于我在评论中提出的另一个问题。

这是我的整个代码：

batch_tpose <- function(filename) {
dat <- matrix(scan(file = filename),ncol=3,byrow=T)
dat_trans <- matrix(dat[,3], 141, byrow=FALSE)

header <- "NCOLS 131
NROWS 141
XLLCENTER -150
YLLCENTER -70
CELLSIZE 1
NODATA_VALUE -32768"

writeLines(header, filename)
write.table(dat_trans, file = filename, row.names=FALSE, col.names=FALSE, append=TRUE)

## This works, I just need the output to be something like out/filename...


}

batch_tpose("200RelVort-nwatl-20201008-15.txt")

我希望能够batch_tpose在底部的命令中使用list.files. 现在我只是用一个文件进行测试。这有效，但它将文件写入工作目录。它还会覆盖源文件，因为它位于同一目录中。如何将输出文件写入“out”目录？我已经尝试过out/filename，file = out/filename但这会创建一个名为“文件名”的新文件。我尝试了其他组合，但它们导致了错误。

4

1 回答 1

1

您可以先用（或类似函数）写入标题，readLines然后再写入数据。write.table诀窍就是选项append = TRUE。这里有一个例子：

header <- "NCOLS 131
NROWS 141
XLLCENTER -150
YLLCENTER -70
CELLSIZE 1
NODATA_VALUE -32768"

data <- matrix(rnorm(100), ncol=10)

writeLines(header, "file.txt")
write.table(data, file="file.txt", col.names=FALSE, row.names=FALSE, append=TRUE)

您还可以paste将标题文本与通常的 R 函数一起使用：

ncols=131
nrows=141

header <- paste0("NCOLS ", ncols, 
"\nNROWS ", nrows,
"\nXLLCENTER -150
YLLCENTER -70
CELLSIZE 1
NODATA_VALUE -32768")

换行符可以指定为文本字符串中的换行符，也可以使用换行符“\n”代码。

于 2021-06-02T19:09:26.470 回答