I am trying to mosaic very large rasters in R. My code, taken from this helpful answer, worked perfectly for 3 test rasters (50 MB, 50 MB, 120 MB), but when I tried a larger test on 7 rasters (ranging from 1 MB to 400 MB) I returned an error.
setwd("~/raster_tiles") gdalbuildvrt(gdalfile = "*.tif", # uses all tiffs in the current folder output.vrt = "test.vrt") gdal_translate(src_dataset = "test.vrt", dst_dataset = "test.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")) Error in .local(.Object, ...) : In addition: Warning message: In system(cmd, intern = TRUE) : running command '"C:\Program Files\QGIS 2.18\bin\gdal_translate.exe" -of "GTiff" "geom1to7test.vrt" "geom1to7test.tif"' had status 1 Error in .rasterObjectFromFile(x, objecttype = "RasterBrick", ...) : Cannot create a RasterLayer object from this file. (file does not exist)
I can see a .vrt file in my folder of ~19 MB, but can't seem to return a GeoTIFF. Does it only have to do with the file size when I use gdal_translate on more and bigger rasters? Any suggestions?