<div dir="ltr"><div><div>Thank you all for the suggestions!<br></div>Cheers,<br></div>madi<br><div class="gmail_extra"><br><div class="gmail_quote">On Mon, Aug 8, 2016 at 9:01 PM, Blumentrath, Stefan <span dir="ltr"><<a href="mailto:Stefan.Blumentrath@nina.no" target="_blank">Stefan.Blumentrath@nina.no</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">





<div link="blue" vlink="purple" lang="NO-BOK">
<div>
<p><span lang="EN-US">Hei Madi,<u></u><u></u></span></p>
<p><span lang="EN-US"><u></u> <u></u></span></p>
<p><span lang="EN-US">Given that you are quite familiar with GRASS GIS here a GRASS / GDAL solution:<u></u><u></u></span></p>
<p><span lang="EN-US"><u></u> <u></u></span></p>
<p style="margin-left:36.0pt">
<u></u><span lang="EN-US"><span>1.<span style="font:7.0pt "Times New Roman"">     
</span></span></span><u></u><span lang="EN-US">Patch all tiles into one VRT with gdalbuildvrt
<u></u><u></u></span></p>
<p style="margin-left:36.0pt">
<u></u><span lang="EN-US"><span>2.<span style="font:7.0pt "Times New Roman"">     
</span></span></span><u></u><span lang="EN-US">Link the VRT and the individual tiles to GRASS (r.external)<u></u><u></u></span></p>
<p style="margin-left:36.0pt">
<u></u><span lang="EN-US"><span>3.<span style="font:7.0pt "Times New Roman"">     
</span></span></span><u></u><span lang="EN-US">Loop over the tiles and adjust you region accordingly and apply what GRASS module you like.<u></u><u></u></span></p>
<p><span lang="EN-US"><u></u> <u></u></span></p>
<p><span lang="EN-US">So something like this:<u></u><u></u></span></p>
<p><span lang="EN-US"><u></u> <u></u></span></p>
<p><span lang="EN-US">gdalbuildvrt ...<u></u><u></u></span></p>
<p><span lang="EN-US"><u></u> <u></u></span></p>
<p><span lang="EN-US">r.external input=vrt ...<u></u><u></u></span></p>
<p><span lang="EN-US">r.external input=tile ...<u></u><u></u></span></p>
<p><span lang="EN-US"><u></u> <u></u></span></p>
<p><span lang="EN-US">for tile in $(g.list type=raster ...)<u></u><u></u></span></p>
<p><span lang="EN-US">do<u></u><u></u></span></p>
<p><span lang="EN-US">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)<u></u><u></u></span></p>
<p><span lang="EN-US">r.*<u></u><u></u></span></p>
<p><span lang="EN-US">done<u></u><u></u></span></p>
<p><span lang="EN-US"><u></u> <u></u></span></p>
<p><span lang="EN-US">Cheers<u></u><u></u></span></p>
<p><span lang="EN-US">Stefan<u></u><u></u></span></p>
<p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Calibri",sans-serif;color:#1f497d" lang="EN-US"><u></u> <u></u></span></p>
<p class="MsoNormal"><b><span style="font-size:11.0pt;font-family:"Calibri",sans-serif" lang="EN-US">From:</span></b><span style="font-size:11.0pt;font-family:"Calibri",sans-serif" lang="EN-US"> gdal-dev [mailto:<a href="mailto:gdal-dev-bounces@lists.osgeo.org" target="_blank">gdal-dev-bounces@<wbr>lists.osgeo.org</a>]
<b>On Behalf Of </b>Margherita Di Leo<br>
<b>Sent:</b> 8. august 2016 15:43<br>
<b>To:</b> Ari Jolma <<a href="mailto:ari.jolma@gmail.com" target="_blank">ari.jolma@gmail.com</a>><br>
<b>Cc:</b> <a href="mailto:gdal-dev@lists.osgeo.org" target="_blank">gdal-dev@lists.osgeo.org</a><br>
<b>Subject:</b> Re: [gdal-dev] Tiling with overlap<u></u><u></u></span></p><div><div class="h5">
<p class="MsoNormal"><span lang="EN-US"><u></u> <u></u></span></p>
<div>
<p class="MsoNormal">Ari,<u></u><u></u></p>
<div>
<p class="MsoNormal"><u></u> <u></u></p>
<div>
<p class="MsoNormal">On Mon, Aug 8, 2016 at 3:01 PM, Ari Jolma <<a href="mailto:ari.jolma@gmail.com" target="_blank">ari.jolma@gmail.com</a>> wrote:<u></u><u></u></p>
<blockquote style="border:none;border-left:solid #cccccc 1.0pt;padding:0cm 0cm 0cm 6.0pt;margin-left:4.8pt;margin-right:0cm">
<div>
<p class="MsoNormal">08.08.2016, 15:19, Margherita Di Leo kirjoitti:<br>
<br>
<u></u><u></u></p>
<blockquote style="margin-top:5.0pt;margin-bottom:5.0pt">
<div>
<div>
<div>
<p class="MsoNormal" style="margin-bottom:12.0pt">Dear all,<u></u><u></u></p>
</div>
<p class="MsoNormal">is there a simple way to create tiles with a certain overlap?<u></u><u></u></p>
</div>
</div>
</blockquote>
<p class="MsoNormal"><br>
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.<u></u><u></u></p>
</div>
</blockquote>
<div>
<p class="MsoNormal"><u></u> <u></u></p>
</div>
<div>
<p class="MsoNormal" style="margin-bottom:12.0pt">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.<br>
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"?<u></u><u></u></p>
</div>
<div>
<p class="MsoNormal">Thank you<u></u><u></u></p>
</div>
<div>
<p class="MsoNormal"><br>
<br>
[1] <a href="http://topotools.cr.usgs.gov/gmted_viewer/" target="_blank">http://topotools.cr.usgs.gov/<wbr>gmted_viewer/</a>
<u></u><u></u></p>
</div>
</div>
<p class="MsoNormal"><br clear="all">
<br>
-- <u></u><u></u></p>
<div>
<div>
<div>
<div>
<div>
<div>
<div>
<div>
<div>
<p class="MsoNormal"><span style="color:#666666">Margherita Di Leo</span><u></u><u></u></p>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
</div></div></div>
</div>

</blockquote></div><br><br clear="all"><br>-- <br><div class="gmail_signature" data-smartmail="gmail_signature"><div dir="ltr"><div><div dir="ltr"><div><div dir="ltr"><div><div dir="ltr"><div><font color="#666666">Margherita Di Leo</font></div></div></div></div></div></div></div></div></div>
</div></div>