Hi, I wish to call gdalwarp from within a python script but I cannot get it to work.
The command I wish to run is: 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 I'm running on Ubuntu 12.04. The above command works perfectly when called from a terminal window and produces a good output image. But when I run from python I get problems. First, I tried to run with subprocess.call(["gdalwarp","-q -multi ... TEST2.tiff"]) then all I see is the default usage instructions. No file is created. Same if I use popen. If I use: os.system("gdalwarp -q -multi ... TEST2.tiff") then it does run and I see the expected text on the screen. However, in this case the output file is blank - it is the right size and contains the right georeferencing information but every pixel in the image has a value of zero. Does anyone have an idea what could be the problem here? Thanks! -- View this message in context: http://osgeo-org.1560.x6.nabble.com/Problem-calling-gdalwarp-from-within-python-script-tp5098902.html Sent from the GDAL - Dev mailing list archive at Nabble.com. _______________________________________________ gdal-dev mailing list gdal-dev@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/gdal-dev