I have a raster dataset that is made up of 280 tiles, 51Gb total. I need to smooth/blur the data, retaining the same resolution. I would like to blur across tile boundaries, where possible, so that points on the edge of each raster contain information from the neighbouring raster(s). I would like the output to be in the same format or smaller tiles.
I there an easy way to do this? I see that it's possible to blur using a KernelFilteredSource in VRT, but it's not clear how that can be made to work across tile boundaries.
I guess the other option is to write a script to load each tile, plus the margins from the surrounding 8 tiles, blur the result, trim, and re-save. I'm just wondering if there's a more efficient way to do it before I try that.