Not a member of Pastebin yet?
Sign Up,
it unlocks many cool features!
- import arcpy
- pnt = arcpy.Point()
- pntX = 0;
- pntY = 0;
- array = arcpy.Array()
- mxd = arcpy.mapping.MapDocument("CURRENT")
- spatial_reference = mxd.activeDataFrame.spatialReference
- pntFc = r"C:UsersAdministratorDesktopCESIMStest.shp"
- df = arcpy.mapping.ListDataFrames(mxd)[0]
- for lyr in arcpy.mapping.ListLayers(mxd, "", df):
- if lyr.name == "testLine":
- for row in arcpy.da.SearchCursor(lyr, ["SHAPE@"]):
- partnum = 0
- for part in row[0]:
- for pnt in part:
- if pnt:
- pointGeom = arcpy.Point()
- pointGeom.x = pnt.X
- pointGeom.y = pnt.Y
- array.add(pointGeom)
- print("{}, {}".format(pnt.X, pnt.Y))
- partnum += 1
- print "one"
- if lyr.name == "test":
- for row in arcpy.da.SearchCursor(lyr, ["SHAPE@XY"]):
- x, y = row[0]
- pntX =x
- pntY =y
- pnt.x = x
- pnt.y = y
- print "two"
- plyLne = arcpy.Polyline(array)
- pnt_geometry = arcpy.PointGeometry(arcpy.Point(pntX,pntY), spatial_reference)
- distance = plyLne.distanceTo(pnt_geometry)
- split_lines =
- arcpy.management.SplitLineAtPoint(plyLne,pnt,arcpy.Geometry(),distance)
- ptGeoms = [arcpy.Geometry ("POINT", split_lines, sr)]
- with arcpy.da.InsertCursor (pntFc, "SHAPE@") as iCurs:
- iCurs.insertRow (ptGeoms)
- arcpy.RefreshActiveView ()
- print "done"
- distance = plyLne.distanceTo(pnt_geometry)
Add Comment
Please, Sign In to add comment