https://community.esri.com/t5/python-questions/extend-line-not-within-up-to-polygon-feature/m-p/1186583#M64833 <P>Thank you, I will try it.</P> Mon, 27 Jun 2022 04:25:18 GMT JohannesBierer 2022-06-27T04:25:18Z https://community.esri.com/t5/python-questions/extend-line-not-within-up-to-polygon-feature/m-p/1185381#M64807 <P>Hello,</P><P>I've found a script of JoshuaBixby from here:</P><P><A href="https://community.esri.com/t5/python-questions/extend-line-up-to-polygon-feature/td-p/705806/page/2" target="_blank">https://community.esri.com/t5/python-questions/extend-line-up-to-polygon-feature/td-p/705806/page/2</A></P><P>it work well to extend lines within a polygon up to the border of the polygon feature. But it couldn't work with lines which are somehow between the polygons (see attached image - in black original line).</P><P>Could you think of a solution for this?</P><P>I tested "Have their center in" in line 21 instead of "within" but it doesn't work ...</P><P>I use notebook in ArcGIS Pro 2.9</P><LI-CODE lang="python"># -*- coding: utf-8 -*- import arcpy import math arcpy.env.workspace = FGDB arcpy.env.overwriteOutput = True dist = 100 fc_line = "line feature class" fc_line_1 = "Sperren_C_1" arcpy.CopyFeatures_management(fc_line, fc_line_1) fc_line_FL = arcpy.MakeFeatureLayer_management(fc_line_1, "fc_line") fc_poly = "Polygon Feature Class" fc_poly_FL = arcpy.MakeFeatureLayer_management(fc_poly, "fc_poly") with arcpy.da.SearchCursor(fc_poly_FL, "SHAPE@") as scur: for poly, in scur: boundary = poly.boundary() arcpy.SelectLayerByLocation_management(fc_line_FL, "WITHIN", poly) with arcpy.da.UpdateCursor(fc_line_FL, "SHAPE@") as ucur: for line, in ucur: arr, = line.getPart() SR = line.spatialReference p1, p2 = arr[0], arr[1] angle = math.atan2(p2.Y - p1.Y, p2.X - p1.X) p = arcpy.Point(p1.X - dist * math.cos(angle), p1.Y - dist * math.sin(angle)) arr.insert(0, p) pn1, pn = arr[len(arr)-2], arr[len(arr)-1] angle = math.atan2(pn.Y - pn1.Y, pn.X - pn1.X) p = arcpy.Point(pn.X + dist * math.cos(angle), pn.Y + dist * math.sin(angle)) arr.append(p) line = arcpy.Polyline(arr, SR, True, False) line = line.cut(boundary)[1] arr, = line.getPart() p1, pn = arr[0], arr[len(arr)-1] p1.Z, pn.Z = 0, 0 arr.insert(0, p1) arr.append(pn) line = arcpy.Polyline(arr, SR, True, False) ucur.updateRow([line])</LI-CODE><P>&nbsp;</P> Thu, 23 Jun 2022 08:12:44 GMT https://community.esri.com/t5/python-questions/extend-line-not-within-up-to-polygon-feature/m-p/1185381#M64807 JohannesBierer 2022-06-23T08:12:44Z https://community.esri.com/t5/python-questions/extend-line-not-within-up-to-polygon-feature/m-p/1185979#M64821 <P>Maybe you can store the length of the line (either original line length or after the extension) and then check it after it is cut by the boundary to see if there is a difference and go from there?</P> Fri, 24 Jun 2022 13:06:30 GMT https://community.esri.com/t5/python-questions/extend-line-not-within-up-to-polygon-feature/m-p/1185979#M64821 Anonymous User 2022-06-24T13:06:30Z https://community.esri.com/t5/python-questions/extend-line-not-within-up-to-polygon-feature/m-p/1186583#M64833 <P>Thank you, I will try it.</P> Mon, 27 Jun 2022 04:25:18 GMT https://community.esri.com/t5/python-questions/extend-line-not-within-up-to-polygon-feature/m-p/1186583#M64833 JohannesBierer 2022-06-27T04:25:18Z