[Python] Rasterizzare un layer vettoriale usando la libreria GDAL

[Maurizio Minissale]

Ciao a tutti,

quello che sto cercando di fare è ottenere un layer raster partendo da un
vettoriale.
Più precisamente ho un vettoriale contenente poligoni e vorrei rasterizzare
l'informazione di uno specifico campo, ovvero il campo 'Index' (che ho
aggiunto al layer vettoriale originale) contenente valori interi compresi
tra 0 e 9.

Ho utilizzato l'istruzione:

err = gdal.RasterizeLayer(target_ds, (3,2,1), source_layer,
burn_values=(0,0,0), options=["ATTRIBUTE= Index"])

che effettua la rasterizzazione partendo da un DEM (target_ds) target, ma
mi ritrovo una matrice piena di 0 e non i valori del campo.

Qualcuno ha qualche suggerimento in merito?
Purtroppo la documentazione a riguardo non è molto chiara...

Grazie in anticipo,
Maurizio.
-------------- parte successiva --------------
Un allegato HTML è stato rimosso...
URL: <http://lists.python.it/pipermail/python/attachments/20130114/e0660f24/attachment.html>