thanks, That does the trick. Oz. On Mon, Mar 30, 2009 at 7:02 PM, Even Rouault <even.roua...@mines-paris.org>wrote:
> Oz, > > See the hereafter a correction to insert in your code that will hopefully > fix > your issue. > > Le Monday 30 March 2009 17:25:13 Oz Nahum, vous avez écrit : > > Hi, > > I have encountered a weired problem with gdal_rasterize. > > > > when I run the script from here > > > http://trac.osgeo.org/gdal/wiki/FAQRaster#HowcanIcreateablankrasterbasedona > >vectorfilesextentsforusewithgdal_rasterize > > > > When I run: > > #!/usr/bin/env python > > > > from osgeo import gdal > > from osgeo import osr > > from osgeo import ogr > > import numpy > > #from numpy import ones > > import stat, sys, os, string, commands > > > > shp = 'AOI' > > tiff = 'AOI.tif' > > px = 1 > > tiff_width = 80 > > tiff_height = 80 > > > > # Import vector shapefile > > vector = ogr.GetDriverByName('ESRI Shapefile') > > src_ds = vector.Open(shp +'.shp') > > src_lyr = src_ds.GetLayerByIndex(index=0) > > src_extent = src_lyr.GetExtent() > > > > # Create new raster layer with 4 bands > > raster = gdal.GetDriverByName('GTiff') > > dst_ds = raster.Create( tiff, tiff_width, tiff_height, 1, gdal.GDT_Byte) > > > > # Create raster GeoTransform based on upper left corner and pixel > > resolution raster_transform = [src_extent[0], px, 0.0, src_extent[3], > 0.0, > > -px] dst_ds.SetGeoTransform( raster_transform ) > > > > # Get projection of shapefile and assigned to raster > > srs = osr.SpatialReference() > > srs.ImportFromWkt(src_lyr.GetSpatialRef().__str__()) > > dst_ds.SetProjection( srs.ExportToWkt() ) > > > > # Create blank raster with fully opaque alpha band > > zeros = numpy.zeros( (tiff_height, tiff_width), numpy.uint8 ) > > dst_ds.GetRasterBand(1).WriteArray( zeros ) > > # Properly close the dataset before calling gdal_rasterize, so that the TIF > # file is properly written. > dst_ds = None > > > > > commandString = 'gdal_rasterize -burn 255 -l AOI AOI.shp AOI.tif' > > commandOutput = commands.getoutput(commandString) > > print commandOutput > > #end of code > > > > I get the error this error: > > > > augmented to the gdal_rasterize command I see this error. When I run it > via > > a shell script, I don't see the error: > > ERROR 4: `AOI.tif' not recognised as a supported file format. > > > > > > #!/bin/sh > > python create_empty_raster.py > > gdal_rasterize -burn 255 -l AOI AOI.shp AOI.tif > > > > Runs ok. > > --> This runs ok as at the end of the create_empty_raster.py script, the > GDAL > dataset is automatically closed. > > > > > So I am confused, why I can't call gdal_rasterize from python ? > > Thanks, for any help > > > > Oz Nahum > > Universtät Tuebingen > > Applied Environmental Sciences > > > -- ---- Imagine there's no countries It isn't hard to do Nothing to kill or die for And no religion too Imagine all the people Living life in peace ---- You all must read 'The God Delusion' http://en.wikipedia.org/wiki/The_God_Delusion --- when one person suffers from a delusion it is called insanity. When many people suffer from a delusion it is called religion." Robert Pirsig, Zen and the Art of Motorcycle Maintenance
_______________________________________________ gdal-dev mailing list gdal-dev@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/gdal-dev