У меня есть 50 файлов geotiff в той же папке. Все они представляли данные о высотах в части мира. Я хотел бы объединить некоторые файлы geotiff, и я обнаружил, что мозаика в R может нам помочь. Я переместил эти geotiff в ту же папку и написал R-скрипт, показанный ниже:
setwd()
a<-lista.files(pattern="*.tiff",file.name=TRUE)
combind<-merge(a,fun=mean)
Однако этот скрипт вернул ошибку: ошибка в as.data.frame(y)
Могу я спросить, как я могу улучшить свой сценарий?
Вы можете использовать мощные функции 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, поскольку растровые файлы, вероятно, не будут иметь одинаковое происхождение.
options = c("BIGTIFF=YES")
- person loki; 08.05.2018
gdalbuildvrt(gdalfile = ".tif$" ...) не удалось. Я отредактировал свой ответ. Возможно, вы могли бы проверить. Просто замените ".tif$" на ".tif". Я также добавил некоторые параметры.
- person loki; 09.05.2018
do.call(merge, a), а дляmergeнет аргументаfun. - person akrun schedule 08.05.2018Reduce(function(...) merge(...), a)- person akrun schedule 08.05.2018?mergeизraster, пример работает сdo.call(merge, x), и я также не нахожу аргументаfun- person akrun schedule 08.05.2018