CSV TO SHAPEFILE - Использование pyshp - PullRequest
Новая
       28

CSV TO SHAPEFILE - Использование pyshp

0 голосов
/ 05 сентября 2018

Мне нужна помощь.

Я получил несколько CSV-файлов с географической информацией, и мне нужно преобразовать их в SHP ... Я пытаюсь использовать следующий код, но он не работает ..

Я впервые работаю с python. Буду очень признателен, если вы поможете мне. 

import shapefile as shp
import csv

out_file = 'test.shp'

#Set up blank lists for data
lat,lon=[],[]

#read data from csv file and store in lists
with open('input.csv', 'rt') as csvfile:
    r = csv.reader(csvfile, delimiter=';')
    for i,row in enumerate(r):
        if i > 0: #skip header
            lat.append((row[0]))
            lon.append((row[1]))

#Set up shapefile writer and create empty fields
w = shp.Writer(shp.POINT)
w.autoBalance = 1 #ensures gemoetry and attributes match
w.field('lat','F',10,8)
w.field('lon','F',10,8)
w.field('Date','D')
w.field('Target','C',50)
w.field('ID','N')

#loop through the data and write the shapefile
for j,k in enumerate(x):
    w.point(k,y[j]) #write the geometry
    w.record(k,y[j],date[j], target[j], id_no[j]) #write the attributes

#Save shapefile
w.save(out_file)

Вот что мне возвращает шпион: 

Python 3.6.5 |Anaconda, Inc.| (default, Mar 29 2018, 13:32:41) [MSC v.1900 64 bit (AMD64)]
Type "copyright", "credits" or "license" for more information.

IPython 6.4.0 -- An enhanced Interactive Python.
Traceback (most recent call last):

  File "<ipython-input-1-dece8279ec6d>", line 1, in <module>
    runfile('C:/Users/gcnxq/Downloads/py/csv.py', wdir='C:/Users/gcnxq/Downloads/py')

  File "C:\Users\gcnxq\AppData\Local\Continuum\anaconda3\lib\site-packages\spyder\utils\site\sitecustomize.py", line 705, in runfile
    execfile(filename, namespace)

  File "C:\Users\gcnxq\AppData\Local\Continuum\anaconda3\lib\site-packages\spyder\utils\site\sitecustomize.py", line 102, in execfile
    exec(compile(f.read(), filename, 'exec'), namespace)

  File "C:/Users/gcnxq/Downloads/py/csv.py", line 18, in <module>
    w = shp.Writer(shp.POINT)

  File "C:\Users\gcnxq\Downloads\py\shapefile.py", line 1018, in __init__
    self.shp = self.__getFileObj(os.path.splitext(target)[0] + '.shp')

  File "C:\Users\gcnxq\AppData\Local\Continuum\anaconda3\lib\ntpath.py", line 224, in splitext
    p = os.fspath(p)

TypeError: expected str, bytes or os.PathLike object, not int

1 Ответ

0 голосов
/ 20 октября 2018

Мне удалось сделать это, используя arcpy модули MakeXYEventLayer_management, а затем CopyFeatures_management. Пример ниже 

import arcpy 
arcpy.env.workspace = r"C:/home"
csvFile = r"C:/home/some.csv"
in_x_field = "long"
in_y_field = "lat"
out_layer = "csvLayer"
spatial_reference = 4326

csvLayer=arcpy.MakeXYEventLayer_management(csvFile,in_x_field,in_y_field,out_layer,spatial_reference)
arcpy.CopyFeatures_management(csvLayer,"shapename.shp")
...