r - 使用 R 将 geoTIFF 转换为 ascii/text？

Question

我试图从 here 获得净初级生产力的时间序列全局估计值，将它们下载到位于美国的远程机器(我在肯尼亚，手机上网速度很慢)，将它们从当前格式转换为 ascii，使用如下函数粗化它们:

coarsen = function(input,size_out){
    if (dim(input)[1] != dim(input)[2]) {print("Matrix not square.")} else {    
    dimsize = dim(input)[1]
    fac = dimsize/size_out
    x = suppressWarnings(matrix(c(rep(1, fac), rep(0, dimsize)), dimsize, dimsize/fac))
    t(x) %*% input %*% x}
    }

然后下载它们以便我可以使用它们。这是否可以使用 R 而不是专门构建的 GIS 程序来完成？如果是这样，我该怎么做？首选 R 的解决方案——我将使用一个无法安装程序的集群，但我可以安装 R 包。

Answer 1

光栅 包具有出色的工具，可以读取和写入多种不同格式的网格数据，以及对此类数据执行各种操作。

使用它，您案例中的基本工作流程(对于单个文件)可能是这样的:

library(raster)
f <- "path/to/downloaded/file.tif"
r <- raster(f)
ra <- aggregate(r, fact=2)  ## By default aggregates using mean, but see fun=
writeRaster(ra, "path/to/outfile.asc", format="ascii")

该软件包有据可查。我建议从 ?"raster-package" 开始和 vignette(Raster) (加上 help(package="raster") 中列出的所有帮助文件)，然后如果遇到严重障碍，请返回此处。