====== 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"