Hello, 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.
Requirements ============ Application takes as input a number of image files +/- geo-referenced Input may have different resolutions, SRS, bands The user can alter projection and geo-reference parameters Output consists of tiles resulting from all those images; The arrangement of the tiles can be customized (the way tiles are arranged and named) gdalbuildvrt, gdalwarp and gdal2tiles.py indicate that these requirements can be satisfied by gdal. Implementation ============== - for each input dataset create a virtual dataset using GDALAutoCreateWarpedVRT(); this takes into consideration the changes made by the user and target SRS - follow the steps in gdalbuildvrt (CreateVRTNonSeparate) to generate a single virtual file; this involves GDALProxyPoolDatasetCreate, VRTAddSimpleSource (is there a simpler way?) - create base layer tiles from virtual dataset; use GDALDataset::RasterIO() - create tiles higher in the pyramid (gdal2tiles.py) from existing ones Any advice / idea / pointers to higher level functions that do some of the work is highly appreciated. Thanks, Nick
_______________________________________________ gdal-dev mailing list [email protected] http://lists.osgeo.org/mailman/listinfo/gdal-dev
