<div dir="ltr">Hi, <br>I have encountered a weired problem with gdal_rasterize.<br><br>when I run the script from here <a href="http://trac.osgeo.org/gdal/wiki/FAQRaster#HowcanIcreateablankrasterbasedonavectorfilesextentsforusewithgdal_rasterize">http://trac.osgeo.org/gdal/wiki/FAQRaster#HowcanIcreateablankrasterbasedonavectorfilesextentsforusewithgdal_rasterize</a><br>
<br>When I run:<br>#!/usr/bin/env python<br><br>from osgeo import gdal<br>from osgeo import osr<br>from osgeo import ogr<br>import numpy<br>#from numpy import ones<br>import stat, sys, os, string, commands<br><br>shp = 'AOI'<br>
tiff = 'AOI.tif'<br>px = 1<br>tiff_width = 80<br>tiff_height = 80<br><br># Import vector shapefile<br>vector = ogr.GetDriverByName('ESRI Shapefile')<br>src_ds = vector.Open(shp +'.shp')<br>src_lyr = src_ds.GetLayerByIndex(index=0)<br>
src_extent = src_lyr.GetExtent()<br><br># Create new raster layer with 4 bands<br>raster = gdal.GetDriverByName('GTiff')<br>dst_ds = raster.Create( tiff, tiff_width, tiff_height, 1, gdal.GDT_Byte)<br><br># Create raster GeoTransform based on upper left corner and pixel resolution<br>
raster_transform = [src_extent[0], px, 0.0, src_extent[3], 0.0, -px]<br>dst_ds.SetGeoTransform( raster_transform )<br><br># Get projection of shapefile and assigned to raster<br>srs = osr.SpatialReference()<br>srs.ImportFromWkt(src_lyr.GetSpatialRef().__str__())<br>
dst_ds.SetProjection( srs.ExportToWkt() )<br><br># Create blank raster with fully opaque alpha band<br>zeros = numpy.zeros( (tiff_height, tiff_width), numpy.uint8 )<br>dst_ds.GetRasterBand(1).WriteArray( zeros )<br><br>commandString = 'gdal_rasterize -burn 255 -l AOI AOI.shp AOI.tif'<br>
commandOutput = commands.getoutput(commandString)<br>print commandOutput<br>#end of code<br><br>I get the error this error:<br><br>augmented to the gdal_rasterize command I see this error. When I run it via a shell script, I don't see the error:<br>
ERROR 4: `AOI.tif' not recognised as a supported file format.<br><br><br>#!/bin/sh<br>python create_empty_raster.py<br>gdal_rasterize -burn 255 -l AOI AOI.shp AOI.tif<br><br clear="all">Runs ok. <br><br>So I am confused, why I can't call gdal_rasterize from python ?<br>
Thanks, for any help<br><br>Oz Nahum<br>Universtät Tuebingen<br>Applied Environmental Sciences<br>-- <br><br>----<br> Imagine there's no countries<br> It isn't hard to do<br> Nothing to kill or die for<br>
And no religion too<br> Imagine all the people<br> Living life in peace <br><br>
</div>