[GIS] Python: how to create a point shapefile from a text file

pythonshapefile

I'm writing a Python code to read the points in a polygon shapefile and save them to a point shapefile.

I'm not sure about my approach so if you have a better idea (even if it is totally different from mine) please let me know.

So first I made a text file and stored the points' (x,y) in that .txt file. then I tried to make a point shapefile from the text file but it gave an error.

Here is the code (just the last part):

creat point shape-file from text file 
import fileinput
import string
import os
env.overwriteOutput=True
outpath="C:/roadpl"
newfc="newpoint.shp" 
arcpy.CreateFeatureclass_management(outpath, newfc, "Point")
infile="C:/roadpl/roadL5.txt"
cursor=arcpy.da.InsertCursor(newfc, ["SHAPE@"])
array=arcpy.Array()
for line in fileinput.input(infile):
    point.X, point.Y = line.split() 
    line_array.add(point)
cursor.insertRow([arcpy.Point(array)])
fileinput.close()
del cursor

Here is the error:

Traceback (most recent call last):
  File "C:\roadpl\P_Code_L5", line 49, in <module>
    point.X, point.Y  = line.split()
  File "C:\Program Files (x86)\ArcGIS\Desktop10.1\arcpy\arcpy\arcobjects\_base.py", line 87, in _set
    return setattr(self._arc_object, attr_name, cval(val))
RuntimeError: Point: Input value is not numeric

I made the input file (.txt file) from a polygon shapefile using:

cur = arcpy.SearchCursor("C:/roadpl/L5P.shp", ["SHAPE@"])
f=open("C:/roadpl/roadL5.txt", "w")

for row in cur:
    geom = row.shape
    row = geom.getPart()
    for part in row:
        for point in part:
          print >>f, (point.X, point.Y)
f.close()

Best Answer

Your X and Y values are text and need to be converted to numeric values. You can tweak your code as follows:

for line in fileinput.input(infile):
    X,Y=line.split()
    cursor.insertRow([arcpy.Point(float(X),float(Y))])

EDIT
I had been focusing on the fact that your input was text but on closer inspection of your code I think I have spotted another error. As I understand it, you want to import each pair of coordinates as a separate point. At the moment you are adding all the points into a single array and then trying to define that whole array as a point. A point can only have two coordinates. So, see the revised code above.

Related Question