OK, just in case you're interested, but I'm sure Dan will come up with a better solution:
import arcpy
import os
import math
txt = r"D:\Xander\GeoNet\PointCloud\sample.txt"
fc_pnt = r"D:\Xander\GeoNet\PointCloud\gdb\test.gdb\points05"
fc_lines = r"D:\Xander\GeoNet\PointCloud\gdb\test.gdb\lines05"
max_dist = 25
arcpy.env.overwriteOutput = True
# sr = arcpy.SpatialReference()
# create empty fc
pnt_ws, pnt_name = os.path.split(fc_pnt)
arcpy.CreateFeatureclass_management(pnt_ws, pnt_name, "POINT", "#", "DISABLED", "ENABLED")
fld_X = "X"
fld_Y = "Y"
fld_Z = "Z"
fld_A = "Att"
fld_dist = "Dist"
arcpy.AddField_management(fc_pnt, fld_X, "Double")
arcpy.AddField_management(fc_pnt, fld_Y, "Double")
arcpy.AddField_management(fc_pnt, fld_Z, "Double")
arcpy.AddField_management(fc_pnt, fld_A, "Integer")
arcpy.AddField_management(fc_pnt, fld_dist, "Double")
flds = ("SHAPE@", fld_X, fld_Y, fld_Z, fld_A, fld_dist)
lst_lines = []
with arcpy.da.InsertCursor(fc_pnt, flds) as curs:
# open txt file
arr_line = arcpy.Array()
i = 0
with open(txt, 'r') as f:
for r in f.readlines():
i += 1
r = r.replace('\n', '')
lst = r.split(' ')
if i % 1000 == 0:
print "Processing line: {0}".format(i)
if i == 1:
x = float(lst[3])
y = float(lst[5])
xp = x
yp = y
else:
xp = x
yp = y
x = float(lst[3])
y = float(lst[5])
z = float(lst[10])
a = int(lst[12])
dist = math.hypot(xp-x, yp-y)
pnt = arcpy.Point(x, y, z)
pnt_g = arcpy.PointGeometry(pnt)
curs.insertRow((pnt_g, x, y, z, a, dist, ))
## if dist < max_dist:
## # add to line
## arr_line.add(pnt)
## else:
## # add previous line to list
## if arr_line.count > 1:
## polyline = arcpy.Polyline(arr_line)
## lst_lines.append(polyline)
##
## # create a new line
## arr_line = arcpy.Array()
## arr_line.add(pnt)
if a == 1:
# add to line
arr_line.add(pnt)
else:
# add previous line to list
if arr_line.count > 1:
polyline = arcpy.Polyline(arr_line)
lst_lines.append(polyline)
# create a new line
arr_line = arcpy.Array()
arr_line.add(pnt)
arcpy.CopyFeatures_management(lst_lines, fc_lines)
