Hi Even, hi Anton, Il giorno Wed, 01 Jun 2011 11:32:55 +0200 Even Rouault <[email protected]> ha scritto:
> Selon Anton Korosov <[email protected]>: > > > Hi! > > we are developing software for automatic processing of satellite > > data and decided to base it mostly on GDAL-Python bindings. Can I > > ask you about options to create Warped VRT datasets in Python? > > > > 1. The command > > > > gdalwarp -t_srs '+proj=utm +zone=24 +datum=WGS84' -of VRT -tr 1000 > > 1000 -te -1000000 7000000 2000000 10000000 input.tif output.vrt > > > > creates a VRT file of subclass VRTWarpedDataset. The size, > > resolution and spatial reference is specified by the options of > > that command. Is there any way I can create such VRTWarpedDataset > > using GDAL-Python bindings? > > Yes, with Dataset.Create() by specifying a SUBCLASS=VRTWarpedDataset > creation option. ... But I don't think it will help you much since the > GDALInitializeWarpedVRT() that is then used by gdalwarp to fill the > dataset isn't mapped to SWIG. maybe the gdal.ReprojectImage could help, but I don't have the time to test it right now. best regards -- Antonio Valentino _______________________________________________ gdal-dev mailing list [email protected] http://lists.osgeo.org/mailman/listinfo/gdal-dev
