QGIS clipping virtual raster into tiles using grid layer (polygons)

Question

Using QGIS

I want to clip a 100gb VRT raster (UAS) into tiles using a predefined polygon grid (polygon gpkg). I want each tile file to be named by the value contained in the "grid_square" field inside the grid gpkg.

What raster tool should I use for this?

Answer 1

Are you willing to try a solution in R ?

library(sf)
library(terra)

grid <- st_read(dsn = "polygon_grid.gpkg")
nr = nrow(terrace_grid)
r <- rast("uas.vrt")

lapply(1:nr, function(t) {
    tile <- grid[t,]
    nm <- tile$grid_square
    cropped <- crop(r, tile)
    cropped_name = paste0(nm, ".tif")
    writeRast(cropped, cropped_name)
}