export dbf to csv file

dbf2csv.py

#!/usr/bin/python
import sys, pyproj, shapefile

input_file = sys.argv[1]
output_file = sys.argv[2]
source_spatialref = sys.argv[3]
target_spatialref = sys.argv[4]

sf = shapefile.Reader(input_file)
shapes = sf.shapes()
f = open(output_file, 'w')

def quote(token):
  if ',' in token: return '"' + token + '"' 
  else: return token

f.write(','.join([quote(x[0]) for x in sf.fields] + ['lat', 'lng']) + '\n')

# Examples: 'ESRI:102667', 'WGS84' from spatialreference.org
source_proj = pyproj.Proj(init=source_spatialref, preserve_units=True)
target_proj = pyproj.Proj(proj='latlong', datum=target_spatialref, ellps=target_spatialref)

for idx, record in enumerate(sf.records()):
  lng,lat = pyproj.transform(source_proj, target_proj, shapes[idx].points[0][0], shapes[idx].points[0][1])
  f.write(','.join(map(lambda x: quote(str(x).strip()), record) + [str(lat), str(lng)]) + '\n')

f.close()

pip install pyproj
pip install shapefile