Latest Contributions by DarrenWiens2

Have you read this help page: Joining attributes of one table to another? Also, is this what you mean by join (adding matching columns onto records based on a matching attribute value)? Sometimes people say "join" and mean something else entirely. ... View more

Can you compare statewide point density variance? It looks like the main difference is that point density is more variable (more clumps and voids) in the optimized surface than the more even original surface. ... View more

Here's how you can do this using arcpy and geometry objects: >>> radius = 200000 # buffer radius >>> polys = [] >>> sr = arcpy.Describe("YOUR_POINT_LAYER_NAME_HERE").spatialReference >>> with arcpy.da.SearchCursor("YOUR_POINT_LAYER_NAME_HERE","SHAPE@",'#',sr) as cursor: ...     for row in cursor: ...         pt1 = row[0].centroid ...         pt2 = arcpy.Point(row[0].centroid.X - radius*2, row[0].centroid.Y) ...         pt3 = arcpy.Point(row[0].centroid.X - radius*2, row[0].centroid.Y - radius*2) ...         triangle = arcpy.Polygon(arcpy.Array([pt1,pt2,pt3]),sr) ...         buffer = row[0].buffer(radius) ...         wedge = triangle.intersect(buffer,4) ...         polys.append(wedge) ...         >>> arcpy.CopyFeatures_management(polys,'in_memory\wedges') <Result 'in_memory\\wedges'> It helps that all wedges are the same shape, so the triangle shape used to clip (intersect) the buffer is hard-coded. You could fairly easily alter the script to respond to attribute data. Sorry for swooping in on your fun, Xander ... View more

It matters whether or not you run scripts within an ArcMap session, and within that session whether the script is run as a tool or from the Python window. Your error is actually due to not importing the arcpy library. It is auto-loaded in ArcMap, but not standalone PyScripter. Make ' import arcpy ' (no quotes) the first line in a standalone arcpy script. You are correct, though, that the ' CURRENT ' mxd will not be available standalone; reference it by file path. ... View more

That was my fault - I ran my test on a feature class called "ex_centroid". I've updated the script to say "YOUR_POINT_LAYER_NAME" in the Spatial Join call (so, change that to your own). ... View more

This script will make 100 copies of each feature in the feature class "Your point layer name" (of course, change that to your later name). Change the number of duplicates if you want. This script can go into the Python window, accessed under the geoprocessing menu in ArcMap. ... View more

If I'm reading this right, your main question is how to duplicate geometries. Here's how, using Python: >>> duplicates = 100 >>> points = [] >>> with arcpy.da.SearchCursor("YOUR_POINT_LAYER_NAME",["SHAPE@"]) as cursor: ...    for row in cursor: ...        for i in range(duplicates): ...            points.append(row[0]) ...           >>> arcpy.CopyFeatures_management(points,'in_memory\out_points') <Result 'in_memory\\out_points'> >>> arcpy.SpatialJoin_analysis("out_points","YOUR_POINT_LAYER_NAME",r'in_memory\sp_join_output',"JOIN_ONE_TO_ONE","KEEP_ALL",'#',"ARE_IDENTICAL_TO") <Result 'in_memory\\sp_join_output'> ... View more

The line containing writerow is not within a loop, so it only writes once, after the code exits the for loop. Indent further to include in the loop and write multiple rows. ... View more

Okay, what does "it doesn't work" mean? Is there an error? If so, what is it? ... View more

Does it work if you add a '\' in the path? arcpy.mapping.ExportToPDF(mxd, r"D:\Temp\" + str(pageName) + ".pdf") ... View more

You shouldn't need elevations for your structures - simply subtract a reasonable value from your DEM at each structure. 1.) Polygon to Raster, as Jake suggests. 2.) Subtract a value (e.g. 10m) from your DEM: Con (IsNull("structures"), "DEM", "DEM"-10) ... View more

Here's my five minute solution: 1.) Add and calculate a field in your fiber feature class to equal 1 for all rows. 2.) Convert to raster (requires Spatial Analyst, I believe), using new column as value. Now have fiber raster layer, all with value = 1. 3.) Run Cost Distance tool (requires Spatial Analyst), using new raster as cost surface and hub as source. 4.) Extract resulting cost distance values to node points. This will only transfer distances to nodes if they happen to fall on the fiber raster. If not, need to edit nodes so they fall on raster and rerun. This was done somewhat naively, not knowing if there are geometrical consequences calculating continuous distances between square pixels, but the values seem reasonable. ... View more

If you want to update for each row, place cursor.updateRow(row) at the same indentation level as the if/else. ... View more