我在同一个文件夹中有50个geotiff文件。 所有这些都代表了世界部分地区的高程数据。 我想合并某些geotiff文件,我发现R中的马赛克可以帮助我们。 我已将这些geotiff移动到同一个文件夹中,我在下面写了一个R脚本:
setwd()
a<-lista.files(pattern="*.tiff",file.name=TRUE)
combind<-merge(a,fun=mean)
但是,此脚本返回错误: as.data.frame(y)中的错误
请问如何改进我的剧本?
答案 0 :(得分:2)
您可以使用功能强大的GDAL
功能。根据我的经验,这些比纯R代码快得多。
我的方法是使用library(gdalUtils)
:
首先,构建虚拟光栅文件(vrt):
library(gdalUtils)
setwd(...)
gdalbuildvrt(gdalfile = "*.tif", # uses all tiffs in the current folder
output.vrt = "dem.vrt")
然后,将虚拟栅格复制到实际的物理文件:
gdal_translate(src_dataset = "dem.vrt",
dst_dataset = "dem.tif",
output_Raster = TRUE # returns the raster as Raster*Object
# if TRUE, you should consider to assign
# the whole function to an object like dem <- gddal_tr..
options = c("BIGTIFF=YES", "COMPRESSION=LZW"))
另一个纯粹的(可能较慢的)raster
包解决方案将是:
f <- list.files(path = "your/path", pattern = ".tif$", full.names = TRUE)
rl <- lapply(f, raster)
do.call(merge, c(rl, tolerance = 1))
您必须调整tolerance
,因为光栅文件可能不具有相同的来源。