Latest Contributions by DarrenWiens2

Yes, you do need closed polylines. If ET GeoWizards Polylines to Polygons tool doesn't work, go to the Build Polygons tool in the Polygon toolset and follow those instructions. ... View more

What is it doing? ... View more

"Centroids" can mean different things and are handled differently by different tools. Many tools treat centroids according to the center of the feature's mass (e.g. a donut's centroid would be in the center of the donuthole). Some tools have the option to force the centroid inside (e.g. Feature to Point's Inside option). Likewise, compare the difference between the centroid and trueCentroid properties of the Polygon object. ... View more

You need to convert the line to a polygon (perhaps using Feature To Polygon, I use ET GeoWizards) and then Clip or Intersect. edit: is your ROW a polyline or polygon feature class? If polygon, just Clip or Intersect. ... View more

If you want to join the total area back to the individual polygons, you can add a new dummy field to each table and join on that. Or, rather than Summary Statistics, use Dissolve to combine your polygons to get the total area, then Spatial Join back to individual polygons. If that's not what you're after, please elaborate. ... View more

The most straightforward way would be something like the following, although you may run into trouble with /" (untested): import shutil index = arcpy.GetParameterAsText(0) # get the first input parameter entered into the tool dialog shutil.copy("c:/chris/Test1/CitizenComments.mxd" , "c:/chris/test2/" + index + "CitizenComments.mxd") To better work with paths, use os.path.join: import shutil, os index = arcpy.GetParameterAsText(0) # get the first input parameter entered into the tool dialog shutil.copy("c:/chris/Test1/CitizenComments.mxd" , os.path.join("c:/chris/test2", index + "CitizenComments.mxd") ... View more

Did you specifically follow this instruction? I had the same problem and this part from your link solved it for me: Customize > ArcMap Options > Raster > Raster Dataset > and Un-check "Use world file to define the coordinates of the raster" ... View more

I think you're looking for the Multiple Attributes option. Value Field = ownership category, Variation by Symbol Size = acreage. ... View more

If you're dead set against creating an intermediate feature class, you can do so using Python: >>> points1 = 'Layer1' # input layer 1 ... points2 = 'Layer2' # input layer 2 ... sr = arcpy.Describe(points1).spatialReference # CRS of layer 1 ... lines = [] # temp line list ... dict1 = {row[1]:row[0].centroid for row in arcpy.da.SearchCursor(points1,['SHAPE@','NAME'],spatial_reference=sr)} # load layer 1 to dictionary ... dict2 = {row[1]:row[0].centroid for row in arcpy.da.SearchCursor(points2,['SHAPE@','NAME'],spatial_reference=sr)} # load layer 2 to dictionary ... for k1,v1 in dict1.iteritems(): # loop dictionary 1 ...    if k1 in dict2: # see if key exists in dictionary 2 ...        new_line = arcpy.Polyline(arcpy.Array([v1,dict2[k1]]),sr) # create line ...        lines.append(new_line) # add to line list ... arcpy.CopyFeatures_management(lines,r'in_memory\lines') # write lines to feature class ... View more

dict.setdefault(row[0]+'-'+row[1],[row[0],row[1],row[2]] + ([0]*len(time_dict))) The above line translates to: In the dictionary "dict", create a key using a value like 'Mossel Bay-Asla Park A' (the first two columns in the Excel sheet). Set the value to a list like the following ['Mossel Bay','Asla Park A','Alma Clinic',0,0,0,0,0,0] (enough zeroes to accommodate all of the possible time bins). Your error is probably telling you that in row[0]+'-'+row[1] you are trying to concatenate a number with a string. The '*' in the SearchCursor means all fields are read in order, so row[0] is first field, row[1] is second field, etc. This all means your fields are in a different order than mine (TOWN first, SETTLEMENT second, NAME third, TIME fourth, BUILDINGS fifth). ... View more

I've removed the long comments - I think this should be correct: >>> times = set([]) ... with arcpy.da.SearchCursor("Alma_Clinic_SAA_Stats$","TIME") as cursor: ...    for row in cursor: ...        times.add(row[0]) ... time_dict = {} ... for i,item in enumerate(sorted(times)): ...    time_dict[item] = i + 3 ... dict = {} ... with arcpy.da.SearchCursor("Alma_Clinic_SAA_Stats$","*") as cursor: ...    for row in cursor: ...        dict.setdefault(row[0]+'-'+row[1],[row[0],row[1],row[2]] + ([0]*len(time_dict))) ...        dict[row[0]+'-'+row[1]][time_dict[row[3]]] = row[4] ... final_rows = [] ... for k,v in dict.iteritems(): ...    final_rows.append(v) ... print final_rows ... View more

I believe you're running into the error because your expression does not return anything for certain cases, where the value is <= 299. Add an 'else' statement to always return a value. if ([Summary_Statistics_2.SUM_DrainArea_SF]) > 299 then n = "MAJOR" else n = "" end if ... View more

Why do you bother asking for a parameter (O_4), if you're immediately going to overwrite it with "O"? ... View more

I see this is in the web dev place, but in Desktop, you would run the Intersect tool, with your line feature class as the only input, and POINT as the output type. ... View more