create lines from points by position of the points and orientation

Question

I want to create straight lines based on the points. I used function "Points to Line" which accomplish the task only when the points are sequentially ordered. If they are not ordered, then it creates a zig zag lines (I have included the figures for illustration - Line 2). I have hundreds of lines, and manually editing is not an ideal solution.

How can one create a straight line without manually editing the lines?

I believe one has to take into account of the orientation and the position of the points to address this problem.

[![enter image description here][3]][3]

I also tried running the EMST (Euclidean minimum spanning tree) function in the spanning tree tool. I could connect several lines in the same direction by using "Unsplit Line", but don't know what to do with the lines that connect parallel lines. These are unnecessary lines.

Answer 1

For this task, you'll want to create a dictionary of line numbers as the dictionary's key and point XYs as their value. Then create your new line feature class with the same spatial reference as your point feature class. Iterate through your dictionary, sort your XYs, and then create a polyline geometry object from your sorted points.

This method assumes that all points that should be connected in a single line have a shared ID value.

Given:

Try this script:

import os
import arcpy

#point feature class
pntFc = r"C:\Path\To\PntFc"
#id field
idFld = "Line"
#line feature class
lineFc = r"Path\To\Output\LineFc"

#create point dictionary
di = {}
with arcpy.da.SearchCursor (pntFc, [idFld, "SHAPE@XY"]) as curs:
     for line, xy in curs:
         try: di [line] += [xy]
         except KeyError: di [line] = [xy]

#create new line feature class
sr = arcpy.Describe (pntFc).spatialReference
linePath, lineName = os.path.split (lineFc)
arcpy.CreateFeatureclass_management (linePath, lineName, "POLYLINE", spatial_reference = sr)

#add line ID field to line feature class
fldType = [f.type for f in arcpy.ListFields (pntFc) if f.name == idFld][0]
arcpy.AddField_management (lineFc, idFld, fldType)

#insert cursor
with arcpy.da.InsertCursor (lineFc, [idFld, "SHAPE@"]) as iCurs:
    #iterate dictionary
    for lineId in di:
        #get point xys
        points = di [lineId]
        #create array
        array = arcpy.Array ()
        for x, y in sorted (points):
            array.add (arcpy.Point (x, y))
        #create geometry
        geom = arcpy.Geometry ("polyline", array, sr)
        #insert row
        row = (lineId, geom)
        iCurs.insertRow (row)

Result: