Latest Contributions by J_B__K_

Hi Robert & all ! At first - you are doing such a great work for the community, your widgets, your approach,... I appreciate it! To my problem - when using graphical search (choose the exact tool, let´s say polygon, then "draw" a polygon), the tool (drawing polygon) stay active after completing search. Is there a way how to disable this? How to turn cursor "back to pan" automatically after drawing a polygon? (v3.5.4) Thank you! ... View more

Michael, thank you for your attempts! Is there anybody who can answer - what is wrong with this short, easy script? What should I do to use this method successfully? ... View more

Yes, it is possible to make a user connection to this AcrGIS Server without any problems, without credentials needed (see the image). It is just a basic, "open", ESRI´s sample server, I have put it into my code just to show the error, just because it is easily accessible... [ATTACH=CONFIG]23992[/ATTACH] ... View more

Hi! Mattias is not alone with this issue 🙂 I´m having the same problem - when using keyboard (arrows and enter), the graphic shows, but when using mouse to click the result, the graphic does not show (location and info pop-up are correct)! Any ideas? btw. Robert´s modification (to use just custom layers, not default geocode service) works great, thank you! ... View more

This works for me... just a simple python script... import arcpy, os from arcpy import env env.overwriteOutput = True # Get arguments: #   Input point feature class #   Output polygon feature class #   Buffer distance #   Boolean type: Maintain fields and field values of the input in the output print "initializing" inPoints   = "C://YourDataDirectory//yourPoint.shp" outPath   = "C://YourDataDirectory" outName   = "SqBuffer" bufDist    = 15000 #in meters outPolys   = outPath+"//"+outName+".shp" print "creating new blank feature"                              arcpy.CreateFeatureclass_management(outPath, outName,"POLYGON","","","", "", "","","","") print "created" # Open searchcursor inRows = arcpy.SearchCursor(inPoints) # Open insertcursor outRows = arcpy.InsertCursor(outPolys) # Create point and array objects pntObj = arcpy.Point() arrayObj = arcpy.Array() print "writing geomety" for inRow in inRows: # One output feature for each input point feature     inShape = inRow.shape     pnt = inShape.getPart(0)     # Need 5 vertices for square buffer: upper right, upper left, lower left,     #   lower right, upper right. Add and subtract distance from coordinates of     #   input point as appropriate.     for vertex in [0,1,2,3,4]:         pntObj.ID = vertex         if vertex in [0,3,4]:             pntObj.X = pnt.X + bufDist         else:             pntObj.X = pnt.X - bufDist         if vertex in [0,1,5]:             pntObj.Y = pnt.Y + bufDist         else:             pntObj.Y = pnt.Y - bufDist         arrayObj.add(pntObj)     # Create new row for output feature     feat = outRows.newRow()         # Assign array of points to output feature     feat.shape = arrayObj     # Insert the feature     outRows.insertRow(feat)     # Clear array of points     arrayObj.removeAll() # Delete inputcursor del outRows print "done Is this what you are looking for? ... View more

Very nice - This is actually quite usefull to me - I have done some projects in the past doing network anaylysis on logging roads... And I was looking for a way to apply some factors that would describe how "difficult" the road is to manuver (road type as well as each segment's vertical/horizontal complexity). When I get back to my routing project (months away), I'll try to post some code that uses z values as well. Thanks for posting the code... Glad you've found it usefull! This is exactly the purpose - to make network analysis more precise by considering road´s "curvature"... I'm focusing on this theme in my diploma thesis... Thanks for advice with adding and using fields! ... View more

May be helpful - similar topic (calculating angular changes and average curvature of lines)... http://forums.arcgis.com/threads/49581-get-vertex-coordinates?p=175418&viewfull=1#post175418 ... View more

Sorry for late response - so much another work to do:( Ok, I understand, no problem with posting the code. Just to make it clear - it calculates "average curvature" of each line (angular change devided by length of line). I'm pretty sure there is a room for improvement in the script, I just did not have enough time to work on it... # # This script calculates average "curvature" of roads # # It takes a line feature class (with pre-added fields "angleChng" and "AvCurv"!) # then takes the fist line # then takes first three verticies # and calculates the angular change of line between those three verticies # then it loops through all verticies of the line, calculates angular changes and cumulates them # after going through whole line, it writes down the total angular change, calculates and writes down the average curvature # and then this is repeated for all lines in the feature class # import arcpy, math, datetime, numpy from arcpy import env print "starting" start = datetime.datetime.now() # for calculating time of process lineFC = r"roads.shp" # input fc - with pre-added fields "angleChng" and "AvCurv"! dsc = arcpy.Describe(lineFC) shapeFieldName = dsc.ShapeFieldName updateRows = arcpy.UpdateCursor(lineFC) for updateRow in updateRows: # looping through lines    vertexList = []    shapeObj = updateRow.getValue(shapeFieldName)    partObj = shapeObj.getPart(0)    for pointObj in partObj:       vertexList.append((pointObj.X, pointObj.Y)) # getting coordinates of all line vertices    vertexCount = len(vertexList)    i = 0    angularChange=0    averageCurvature=0    while i + 3 <= vertexCount: # that is until you run out of vertex triplicate pairs to count       x1 = vertexList[0] # coordinates of three vertices       y1 = vertexList[1]       x2 = vertexList[i + 1][0]       y2 = vertexList[i + 1][1]       x3 = vertexList[i + 2][0]       y3 = vertexList[i + 2][1]       u1=x1-x2 # vectors (BA and BC !)       u2=y1-y2       v1=x3-x2       v2=y3-y2       cit=(u1*v1)+(u2*v2) # calculating the angle of vectors       jmen1=math.sqrt((u1*u1)+(u2*u2))       jmen2=math.sqrt((v1*v1)+(v2*v2))       jmen=jmen1*jmen2       cosfi=cit/jmen       fiRad=numpy.arccos(cosfi) # angle of vectors [rad]       fiGrad=fiRad*200/math.pi       fiGrad=200-fiGrad # angle change of one segment of line [grad]       angularChange=angularChange+fiGrad # total angle change of line       i = i + 1    averageCurvature=angularChange/shapeObj.length*1000 # average curvature of line = angle change in one kilometer    updateRow.angleChng=angularChange # adding angle change of line    updateRow.AvCurv=averageCurvature # adding average curvature of line    updateRows.updateRow(updateRow) print "Done in ",datetime.datetime.now() - start, " seconds" del updateRows, updateRow, lineFC Btw.: Maybe a newbie question, but i didn't find any solution yet - Is it possible to add field into fc (with fieldName as a variable) and the to use something like updateRow.fieldName? (if using posted script, it is necessary to have "pre-added" fields - it would be much easier to add them in the beginning of the script...) Thanks! ... View more

Sure you can do it in the cursor... I guess an update cursor so you could then populate your angle statistic field value, right? It'd look something like this then: ... Thanks, my funky statistic 😉 is working perfectly that way... (just had a little problem with "break", but that's already solved).. Thank to all of you for your help! ... View more

Thank you Chris for your answer. I know I can store coordinates of verticies in a list and then just simply read what I want, but I was wondering if there is any "smarter" way... if it is quite a big shapefile/feature class (in my case - roads), it takes some time to write all coordinates into list and it takes quite a lot of memory. What I need is to calculate some characteristics of each line (analytical geometry - vectors and so... in my case I want to get the "angle change" of each line) and store it in new field - it should be better to calculate the characteristic directly, without writing coordinates into list, like "take line - get coordinates of points 1-3 - calculate angle change - store the value - move to points 2-3 and repeat..."... Is this possible? (Or, just to be sure - is there any better way to calculate the "angle change" of each line?:)) ... View more