[gdal-dev] Problem calling gdalwarp from within python script

Even Rouault even.rouault at mines-paris.org
Fri Jan 24 13:17:43 PST 2014


Le vendredi 24 janvier 2014 21:46:50, Simon Richard Proud a écrit :
> Does anyone have any other ideas as to what would be causing gdalwarp to
> output entire blank images when called from a python script? I can get it
> to run using os.system and subprocess.call (using both shell=True and
> shell=False) but always just get a blank image (all pixels=0). This does
> not happen if I run exactly the same command from a terminal window.

I can't see any rational reason for that... Just in case : are you sure you 
are running the same GDAL lib and binary in the terminal and from Python ?

> 
> Thanks,
> Simon
> 
> ________________________________________
> Fra: gdal-dev-bounces at lists.osgeo.org [gdal-dev-bounces at lists.osgeo.org]
> på vegne af Sean Gillies [sean at mapbox.com] Sendt: 20. januar 2014
> 11:56
> Til: gdal-dev at lists.osgeo.org
> Emne: Re: [gdal-dev] Problem calling gdalwarp from within python script
> 
> On 1/20/14, 1:57 AM, Even Rouault wrote:
> > Selon Simon Richard Proud <srp at geo.ku.dk>:
> >> Hi,
> >> Unfortunately that produces a file not found error with filename
> >> "65535". If I remove both arguments "-srcnodata" and "65535" then I go
> >> back to getting the usage guidelines.
> > 
> > You need to separate also the numerical values after -te and -tr
> > 
> > ['gdalwarp','-q','-multi','-r', 'cubicspline','-t_srs', '+proj=geos
> > +lon_0=0.0 +h=35785831', '-te', '-5568748.2758', '-5568748.4774',
> > '5568748.2758', '5568748.2758','-tr', '3000.40316582', '3000.40316582',
> > '-srcnodata', '65535','TEST.tiff', 'TEST2.tiff']
> 
> The shlex.split() function is handy for turning strings that you would
> have passed to os.system() into a list of args for subprocess functions.
> 
>      >>> import shlex
>      >>> shlex.split('gdalwarp -q -multi -r cubicspline -te
> 
> -5568748.2758 -5568748.4774 5568748.2758 5568748.2758')
>      ['gdalwarp', '-q', '-multi', '-r', 'cubicspline', '-te',
> '-5568748.2758', '-5568748.4774', '5568748.2758', '5568748.2758']
> 
> ...
> 
> --
> Sean Gillies
> sean at mapbox.com
> _______________________________________________
> gdal-dev mailing list
> gdal-dev at lists.osgeo.org
> http://lists.osgeo.org/mailman/listinfo/gdal-dev
> 
> 
> _______________________________________________
> gdal-dev mailing list
> gdal-dev at lists.osgeo.org
> http://lists.osgeo.org/mailman/listinfo/gdal-dev

-- 
Geospatial professional services
http://even.rouault.free.fr/services.html


More information about the gdal-dev mailing list