After running the script below the "point" layer that is created does not have any features and of course no angles calculated. Any ideas on why this might be would be appreciated.
Thanks,
Jim
import math
02.... sr = arcpy.Describe("Z77").spatialReference
03.... arcpy.CreateFeatureclass_management('in_memory','points',"POINT",spatial_reference=sr)
04.... arcpy.AddField_management('points','Angle',"DOUBLE")
05.... insCur = arcpy.da.InsertCursor("points",('SHAPE@','Angle'))
06.... with arcpy.da.SearchCursor("Z77",["SHAPE@"]) as cursor:
07.... for row in cursor:
08.... for part in row[0]:
09.... pt_count = 1
10.... for pnt in part:
11.... pnt = arcpy.PointGeometry(pnt,sr)
12.... if pt_count > 1:
13.... if pt_count >2:
14.... distAB = oneBack.distanceTo(twoBack)
15.... distBC = pnt.distanceTo(oneBack)
16.... distAC = pnt.distanceTo(twoBack)
17.... angB = math.degrees(math.acos((((distAB*distAB)+(distBC*distBC))-(distAC*distAC))/(2*distAB*distBC)))
18.... insCur.insertRow((oneBack,angB))
19.... twoBack = oneBack
20.... oneBack = pnt
21.... pt_count += 1
22.... del insCur