Latest Contributions by DarrenWiens2

Others are reporting the same issue, currently with no solution. ... View more

I'm confused by this response - the link does not lead to anything specifically to do with arcpy, rather a tool, which like any tool can be called via arcpy, but as far as I can tell it does not have anything to do with setting locator properties. ... View more

Hi Dan, Thanks for this example. Just to be the devil's advocate, I'll add that I'm still not exactly sold on the practicality of using it. As far as I can tell, you save yourself two lines in the two_arg loop (over making 3 individual calls), and nothing at all in the one-arg loop. result = getattr(arcpy, one_arg[0])(*args) can be rewritten as: result = arcpy.AddXY_management(outputs) and... result_bak = getattr(arcpy, one_arg[1])(result,backups) can be rewritten: result_bak = arcpy.CopyFeatures_management(outputs,backups) Unless you're looping through tools, there is no benefit I can see and it reduces readability (for me, at least). I do admit that this is a slick way to loop through tools, however, you need a special script that requires multiple tools with the exact same combination of parameters in order to take advantage of it. I'd be interested to see a spreadsheet listing tools vs parameter type to see which tools would be compatible with each other. ... View more

Hmmm not sure why it doesn't work for you. Is your line feature class projected? If not, that may affect the calculated distances between points, which would affect the angle calculations. ... View more

Are your expected areas either 1,000,000 (meters) or 10,764,000 (feet) times too large, by chance? ... View more

In ArcMap, copy/paste the code into the Python window (under Geoprocessing menu). Change the name of the layer in the SearchCursor (keep within quotes). It will output to an in_memory layer called lines and add it to your map. You can export that layer to one on disk, or change the script to save to a new location (e.g. substitute r'C:\someFolder\someShapefile.shp' for r'in_memory\lines'). "SHAPE@" is a geometry token - don't change it. You can think of it as "the geometry". ... View more

Can you provide an example where this could be useful? I came up with the following (maybe you don't know if you want to 'AddXY_management' or 'RepairGeometry_management' or any other one-parameter tool), but am having a hard time thinking how you could write much code beyond a simple tool call that applies to multiple tools: import arcpy tool = arcpy.GetParameterAsText(0) input = arcpy.GetParameterAsText(1) arguments = [input] result = getattr(arcpy, tool)(*arguments) ... View more

Please post the script you've got so far. ... View more

You may also be interested in tracing the existing border to get started on your new polygon. If not, perhaps you could provide a picture of your exact problem. ... View more

Although this doesn't answer your current question, you should really look into updating your workflow to take advantage of topological editing‌ (editing shared borders). You can also use topological tools to align a mismatched edge. ... View more

Not exactly sure I've captured what you want, but this script will cycle through all the lines in a polyline feature class, and output lines up until it reaches an acute angle. >>> import math ... lines = [] ... with arcpy.da.SearchCursor("YOUR_POLYLINE_FC_HERE","SHAPE@") as cursor: ...    for row in cursor: ...        for part in row[0]: ...            linePnts = [] ...            writePnt = 1 ...            pt_count = 1 ...            for pnt in part: ...                pnt = arcpy.PointGeometry(pnt) ...                if pt_count > 1: ...                    if pt_count >2: ...                        distAB = oneBack.distanceTo(twoBack) ...                        distBC = pnt.distanceTo(oneBack) ...                        distAC = pnt.distanceTo(twoBack) ...                        angB = math.degrees(math.acos((((distAB*distAB)+(distBC*distBC))-(distAC*distAC))/(2*distAB*distBC))) ...                        if angB < 90: ...                            writePnt = 0 ...                    twoBack = oneBack ...                oneBack = pnt ...                pt_count += 1 ...                if writePnt == 1: ...                    linePnts.append(arcpy.Point(pnt.centroid.X,pnt.centroid.Y)) ...            line = arcpy.Polyline(arcpy.Array(linePnts)) ...            lines.append(line) ... arcpy.CopyFeatures_management(lines,r'in_memory\lines') ... View more

Yes, this is a job for Page Definition Query. Your index layer (species) will cycle through species, and setting the Page Definition Query to match, will only show the matching species. ... View more

This wouldn't take too long to do per folder simply through the GUI and batch processing. One batch for Add Field (area), one batch for Calculate Field, and finally Summary Statistics to add up the areas. Of course, if you have many folders, you'd want to automate it with a model or script. ... View more