损坏的栅格，同时合并多个栅格图块

Question

我有大约 1000 个栅格图块。我正在尝试将它们合并为一个栅格图块。

library(raster)
raster_tile_path <- list.files("file_path", full.names = TRUE,
                                        pattern = ".tif" ))

merge_tile  <- lapply(raster_tile_path, raster)
merge_tile  <- do.call(merge, c(merge_tile, tolerance = 1))

运行后，我收到此错误：

Error in rgdal::getRasterData(con, offset = offs, region.dim = reg, band = object@data@band) : Failure during raster IO

有人已经提出了错误来源：可能是一个坏的或损坏的栅格图块。但似乎没有解决方案来识别或忽略这个特定的栅格图块。

有人知道如何解决这个问题吗？

Answer 1

如果有损坏的文件，你应该可以这样找到它

for (i in 1:length(merge_tile)) {
    x = merge_tile[[i]] * 1
}

发生错误时，图块 i 已损坏

你也可以试试

library(terra)
v <- vrt(raster_tile_path)
writeRaster(v, "file.tif")

损坏的栅格，同时合并多个栅格图块

Corrupt raster, while merging multiple raster tiles

r

raster

r-raster

r-mosaic