Esportazione di shapefile da python

Il seguente script esporta in formato “pseudo csv” i dati dello shapefile e la “bounding box” di ogni singola figura censita

"""
Esportazione dei dati di uno shapefile
"""
import shapelib, dbflib, string
 
SHAPEFILE="nomeshp"
 
shp = shapelib.ShapeFile(SHAPEFILE)
dbf = dbflib.DBFFile(SHAPEFILE)
 
sCount = shp.info()[0]
dCount = dbf.record_count()
 
if sCount <> dCount:
 raise "Errore: numero oggeti in shape %s (%d) diverso da record  dbf (%d) " %(SHAPEFILE,sCount,dCount)
 
 
#stampo l'header!
header = dbf.read_record(0)
#print header.keys()
intestazione = string.join(header.keys(),"|")
intestazione += "|XMIN|YMIN|XMAX|YMAX"
 
print intestazione
 
#scandisco riga per riga il file dbf
for i in range (dbf.record_count() ):
 rc = dbf.read_record(i)
 sc = shp.read_object(i)
 row = ""
 for value in rc.values():
  row += "%s|" %(value)
 #aggiungo la extent/envelope
 (min,max) = sc.extents()
 xmin = min[0]
 ymin = min[1]
 xmax = max[0]
 ymax = max[1]
 
 row += "%f|%f|%f|%f" %(xmin,ymin,xmax,ymax)
 print row
 #print "\n"