import ogr path=r"C:\Users\Rich\Documents\seasonal_water_yield_workspace\aggregated_results_cz.shp" v=ogr.Open(path) l=v.GetLayer() with open('agg_results_cz.csv', 'wb') as out_csv: for feature in l: line = str(feature.GetFID()) line += ','+str(feature.GetField('vri_sum')) line += ','+str(feature.GetField('qb'))+'\n' out_csv.write(line)