AnsweredAssumed Answered

deflection angle script

Question asked by jfritz on Nov 29, 2018
Latest reply on Nov 30, 2018 by jfritz

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 

Outcomes