Pyshp का उपयोग करके .csv फ़ाइल को .shp में कनवर्ट करें?


10

मैं यह समझने की कोशिश कर रहा हूं कि मैं अजगर में सीएसवी मॉड्यूल का उपयोग कैसे कर सकता हूं कि अजगर स्क्रिप्ट के रूप में एक ही फ़ोल्डर में एक सीएसवी फ़ाइल खोलें, और फिर शेपफाइल मॉड्यूल pyshp का उपयोग करके एक आकृति बनाते हैं।

सीएसवी फ़ाइल इस तरह दिखती है, लेकिन रिकॉर्ड की हजार पंक्तियों की एक जोड़ी हो सकती है:

id_nr;date;target;start_lat;start_lon
1;2012-05-21;navpoint 25x;55.123654;13.456954
1;2012-05-23;navpoint 11f;55.143654;12.456954

जवाबों:


14

Pyshp मॉड्यूल हैंग करने के लिए थोड़ा मुश्किल है, लेकिन जब आप इसे ले जाते हैं तो वास्तव में उपयोगी होता है। मैंने एक स्क्रिप्ट लिखी है जो उदाहरण डेटा के सीएसवी में पढ़ता है और सही डेटाटाइप्स की विशेषताओं के रूप में संग्रहीत डेटा के साथ एक आकृति को लिखता है। Pyshp / xbase डेटाटाइपिंग हमेशा मेरे लिए मुश्किल रहा है जब तक कि मुझे xbase प्रारूप के लिए यह उपयोगकर्ता गाइड नहीं मिला और इस प्रश्न के परिणामस्वरूप मैंने अपने ब्लॉग पर संबंधित pyshp डेटासेट के बारे में एक छोटा सा नोट लिखा है, जिसका एक हिस्सा मैंने नीचे चिपकाया है। :

  • C ASCII वर्ण है
  • N एक डबल सटीक पूर्णांक है जो लंबाई में लगभग 18 वर्णों तक सीमित है
  • D, YYYYMMDD प्रारूप में तारीखों के लिए है, जिसमें कोई रिक्त स्थान या हाइफ़न वर्गों के बीच नहीं है।
  • F, N के समान लंबाई सीमा वाले फ्लोटिंग पॉइंट नंबरों के लिए है
  • L तार्किक डेटा के लिए है जो कि आकृति पूर्णांक की विशेषता तालिका में 1 (सच) या 0 (गलत) के रूप में एक छोटे पूर्णांक के रूप में संग्रहीत है। वे मान जो प्राप्त कर सकते हैं वे हैं 1, 0, y, n, Y, N, T, F या पायथन बिल्डिंस True andse

पूरी लिस्टिंग इस प्रकार है:

import shapefile as shp
import csv

out_file = 'GPS_Pts.shp'

#Set up blank lists for data
x,y,id_no,date,target=[],[],[],[],[]

#read data from csv file and store in lists
with open('input.csv', 'rb') as csvfile:
    r = csv.reader(csvfile, delimiter=';')
    for i,row in enumerate(r):
        if i > 0: #skip header
            x.append(float(row[3]))
            y.append(float(row[4]))
            id_no.append(row[0])
            date.append(''.join(row[1].split('-')))#formats the date correctly
            target.append(row[2])

#Set up shapefile writer and create empty fields
w = shp.Writer(shp.POINT)
w.autoBalance = 1 #ensures gemoetry and attributes match
w.field('X','F',10,8)
w.field('Y','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)

आशा है कि ये आपकी मदद करेगा।


बहुत अच्छी स्क्रिप्ट है। मुझे एक त्रुटि मिली क्योंकि इसे पाठ के रूप में नहीं पढ़ा था इसलिए मैंने इस पंक्ति को बदल दिया: csvfile के रूप में खुले ('input.csv', 'rt') के साथ:
विरुद्ध

1
मुझे लगता है कि आप लूप से पहले हेड (आर) का उपयोग करके प्रदर्शन को बेहतर बना सकते हैं, अगर किसी स्टेटमेंट का उपयोग करके चेक करने के बजाय हेडर को छोड़ दें।
रोजीको

@sgrieve - यह स्क्रिप्ट विशिष्ट पूर्व-निर्धारित क्षेत्रों के साथ एक सीएसवी को परिवर्तित करती है। मैं किसी भी csv को एक फीचर क्लास में बदलने के लिए एक सामान्य स्क्रिप्ट चाहता हूँ। शायद इसे प्राप्त करने के लिए उपयोगी आर्कपी कार्य हैं?
वाटरमैन

2

एक विकल्प के रूप में आपको सूचियों में डेटा रखने की आवश्यकता नहीं है।

# import libraries
import shapefile, csv

# create a point shapefile
output_shp = shapefile.Writer(shapefile.POINT)
# for every record there must be a corresponding geometry.
output_shp.autoBalance = 1
# create the field names and data type for each.
# you can insert or omit lat-long here
output_shp('Date','D')
output_shp('Target','C',50)
output_shp('ID','N')
# count the features
counter = 1
# access the CSV file
with open('input.csv', 'rb') as csvfile:
    reader = csv.reader(csvfile, delimiter=',')
    # skip the header
    next(reader, None)
    #loop through each of the rows and assign the attributes to variables
    for row in reader:
        id= row[0]
        target= row[1]
        date = row[2]
        # create the point geometry
        output_shp.point(float(longitude),float(latitude))
        # add attribute data
        output_shp.record(id, target, date)
        print "Feature " + str(counter) + " added to Shapefile."
        counter = counter + 1
# save the Shapefile
output_shp.save("output.shp")

आप इस कार्यान्वयन का एक कार्यशील उदाहरण यहां पा सकते हैं ।

हमारी साइट का प्रयोग करके, आप स्वीकार करते हैं कि आपने हमारी Cookie Policy और निजता नीति को पढ़ और समझा लिया है।
Licensed under cc by-sa 3.0 with attribution required.