Guest User

Untitled

a guest
Nov 18th, 2018
85
0
Never
Not a member of Pastebin yet? Sign Up, it unlocks many cool features!
text 1.55 KB | None | 0 0
  1. import arcpy
  2. pnt = arcpy.Point()
  3. pntX = 0;
  4. pntY = 0;
  5. array = arcpy.Array()
  6. mxd = arcpy.mapping.MapDocument("CURRENT")
  7. spatial_reference = mxd.activeDataFrame.spatialReference
  8. pntFc = r"C:UsersAdministratorDesktopCESIMStest.shp"
  9. df = arcpy.mapping.ListDataFrames(mxd)[0]
  10. for lyr in arcpy.mapping.ListLayers(mxd, "", df):
  11. if lyr.name == "testLine":
  12. for row in arcpy.da.SearchCursor(lyr, ["SHAPE@"]):
  13. partnum = 0
  14. for part in row[0]:
  15. for pnt in part:
  16. if pnt:
  17. pointGeom = arcpy.Point()
  18. pointGeom.x = pnt.X
  19. pointGeom.y = pnt.Y
  20. array.add(pointGeom)
  21. print("{}, {}".format(pnt.X, pnt.Y))
  22. partnum += 1
  23. print "one"
  24. if lyr.name == "test":
  25. for row in arcpy.da.SearchCursor(lyr, ["SHAPE@XY"]):
  26. x, y = row[0]
  27. pntX =x
  28. pntY =y
  29. pnt.x = x
  30. pnt.y = y
  31. print "two"
  32. plyLne = arcpy.Polyline(array)
  33. pnt_geometry = arcpy.PointGeometry(arcpy.Point(pntX,pntY), spatial_reference)
  34. distance = plyLne.distanceTo(pnt_geometry)
  35. split_lines =
  36. arcpy.management.SplitLineAtPoint(plyLne,pnt,arcpy.Geometry(),distance)
  37. ptGeoms = [arcpy.Geometry ("POINT", split_lines, sr)]
  38. with arcpy.da.InsertCursor (pntFc, "SHAPE@") as iCurs:
  39. iCurs.insertRow (ptGeoms)
  40. arcpy.RefreshActiveView ()
  41. print "done"
  42.  
  43. distance = plyLne.distanceTo(pnt_geometry)
Add Comment
Please, Sign In to add comment