<div dir="ltr"><div><div><div><div><div><div><div><div>Hello,<br><br></div>I'm working on an application with requirements described below. Please have a look at what I came up with so far and tell me if I'm on an awfully wrong path.<br>
<br></div><span style="font-family:courier new,monospace">Requirements<br>============</span><br><br></div>Application takes as input a number of image files +/- geo-referenced<br></div><div>Input may have different resolutions, SRS, bands<br>
</div>The user can alter projection and geo-reference parameters<br></div>Output consists of tiles resulting from all those images; <br>The arrangement of the tiles can be customized (the way tiles are arranged and named)<br>
<br>gdalbuildvrt, gdalwarp and gdal2tiles.py indicate that these requirements can be satisfied by gdal.<br><br><br></div><span style="font-family:courier new,monospace">Implementation<br>==============</span><br></div><ul>
<li>for each input dataset create a virtual dataset using <span style="color:rgb(0,0,0)">GDALAutoCreateWarpedVRT()</span><span style="color:rgb(0,0,0)">; this takes into consideration the changes made by the user and target SRS</span></li>
<li><span style="color:rgb(0,0,0)">follow the steps in </span>gdalbuildvrt (CreateVRTNonSeparate) to generate a single virtual file; this involves GDALProxyPoolDatasetCreate, VRTAddSimpleSource (is there a simpler way?)</li>
<li>create base layer tiles from virtual dataset; use GDALDataset::RasterIO()<br></li><li>create tiles higher in the pyramid (gdal2tiles.py) from existing ones</li></ul><p>Any advice / idea / pointers to higher level functions that do some of the work is highly appreciated.</p>
<p><br></p><p>Thanks,</p><p>Nick<br></p></div></div>