[gdal-dev] How to specify extent of a VRTWarpedDataset in Python?

[Antonio Valentino]

Hi Even, hi Anton,

Il giorno Wed, 01 Jun 2011 11:32:55 +0200
Even Rouault <even.rouault at mines-paris.org> ha scritto:

> Selon Anton Korosov <anton.korosov at nersc.no>:
> 
> > 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