Thank you all for the suggestions! Cheers, madi On Mon, Aug 8, 2016 at 9:01 PM, Blumentrath, Stefan < [email protected]> wrote:
> Hei Madi, > > > > 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 > > > > Cheers > > Stefan > > > > *From:* gdal-dev [mailto:[email protected]] *On Behalf Of > *Margherita > Di Leo > *Sent:* 8. august 2016 15:43 > *To:* Ari Jolma <[email protected]> > *Cc:* [email protected] > *Subject:* Re: [gdal-dev] Tiling with overlap > > > > Ari, > > > > On Mon, Aug 8, 2016 at 3:01 PM, Ari Jolma <[email protected]> wrote: > > 08.08.2016, 15:19, Margherita Di Leo kirjoitti: > > Dear all, > > is there a simple way to create tiles with a certain overlap? > > > Do you mean by "overlap" the geographic area the tiles should cover? > gdal2tiles.py creates the tiles for the area the source dataset covers. So > you can control the area through the source dataset. > > > > Thanks for your answer, sorry I'm not sure I understand the methodology > you are proposing. Let me explain what I'm trying to do so that the problem > is clearer. I have a dataset that is a global coverage, e.g. GMTED [1]. > That is a tiled dataset. However for my processing chain I need tiles that > overlap on each other for a certain number of pixels. This means that the > area at the border of each tile should be covered by the tile itself and by > the adjacent one to certain extent. > My method would be the following. I would create a virtual mosaic and then > create the tiles with this rule. From the manual of gdal2tiles.py it is > not clear to me if i can obtain this result. What do you mean by "you can > control the area through the source dataset"? > > Thank you > > > > [1] http://topotools.cr.usgs.gov/gmted_viewer/ > > > > -- > > Margherita Di Leo > -- Margherita Di Leo
_______________________________________________ gdal-dev mailing list [email protected] http://lists.osgeo.org/mailman/listinfo/gdal-dev
