point cloud to polyline feature

TanyaBurman1

Hey Everyone,

 

I have a txt file with x,y,z and 0/1 columns. There are about 35000 point coordinates. When displayed, a pattern is clearly visible. How should i convert it to a line feature so that it can be used as a barrier in IDW interpolation. Please note: I already tried point to line feature tool but that doesn't follow the point pattern.

 

Any help would be appreciated.

 

Thanks

XanderBakker

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)