合并多个栅格图块时栅格损坏

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

有人已经建议了错误来源： https://stackoverflow.com/a/67607770/9101903< /a> 可能是坏的或损坏的光栅图块。但似乎没有解决方案来识别或忽略这个特定的栅格图块。

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

Answer 1

如果有损坏的文件，您应该能够像这样找到它

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

当发生错误时，tile i is损坏

您也可以尝试

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