On 8 August 2016 at 21:01, Blumentrath, Stefan <[email protected]> wrote: > Hei Madi, >
Hi, > > Given that you are quite familiar with GRASS GIS here a GRASS / GDAL > solution: > > > > 1. Patch all tiles into one VRT with gdalbuildvrt > > 2. Link the VRT and the individual tiles to GRASS (r.external) > > 3. Loop over the tiles and adjust you region accordingly and apply what > GRASS module you like. > > > > So something like this: > > > > gdalbuildvrt ... > > > > r.external input=vrt ... > > r.external input=tile ... > > > > for tile in $(g.list type=raster ...) > > do > > g.region -p raster=$tile n=n+x s=s-x w=w-x e=e+x (where x is your desired > overlap in map units) > > r.* > > done > I didn't answer before because I was thinking to go OT but in GRASS you can use t.tile with overlaps https://grass.osgeo.org/grass70/manuals/r.tile.html > > Cheers > > Stefan > -- ciao Luca www.lucadelu.org _______________________________________________ gdal-dev mailing list [email protected] http://lists.osgeo.org/mailman/listinfo/gdal-dev
