Latest Contributions by PeterWilson

Background: I've used the Near Analysis tool which has added the following fields: NEAR_X (X Coordinate) NEAR_Y (Y Coordinate) I'd like to use the following fields from the output of using the Near Analysis tool to generate a Structured NumPy Array: HydroID (Unique ID: Long Integer) NEAR_X (X Coordinate) NEAR_Y (Y Coordinate) I've read through the ESRI online help documentation on: Working with NumPy in ArcGIS TableToNumPyArray , and I'm unable to figure out how to generate a Structured NumPy Array that I can use as input to NumPyArrayToFeatureClass. The Structured NumPy Array format that is required as input to NumPyArrayToFeatureClass is: array = numpy.array([(1, (471316.3835861763, 5000448.782036674)),                      (2, (470402.49348005146, 5000049.216449278))], numpy.dtype([('idfield',numpy.int32),('XY', '<f8', 2)])) From what I gathered the coordinates are stored as tuples within the Structured NumPy Array. I'm still learning to code so any assistance in how to covert my data into a Structured NumPy Array that I may use as input into NumPyArrayToFeatureClass will be appreciated. Regards PeterW ... View more

@I'm trying to assign my Watersheds polygon into regions of non overlapping polygons. I'm currently using the ESRI Spatial Analyst Supplemental Tools New Spatial Analyst Supplemental tools, v1.3 | ArcGIS Blog The Zonal Statistics Tool 2 implements the grouping of the overlapping polygons to generate non overlapping polygons to run the Zonal Statistics, but I'm unable to following the logic that that's being used as I'm still new to Python. I've also found the following link on the Stack Overflow website and download the sample code: arcgis desktop - Exploding overlapping to new non-overlapping polygons? - Geographic Information Systems Stack Exchange I've made amendments to the source code from the Stack Overflow website: import os import arcpy from arcpy import GetParameterAsText fc = GetParameterAsText(0) idName = GetParameterAsText(1) dirname = os.path.dirname(arcpy.Describe(fc).catalogPath) desc = arcpy.Describe(dirname) if hasattr(desc, "datasetType") and desc.datasetType=='FeatureClass':     dirname = os.path.dirname(dirname) arcpy.env.workspace = dirname def countOverlaps(fc,idName):     intersect = arcpy.Intersect_analysis(fc,'intersect')     findID = arcpy.FindIdentical_management(intersect,"explFindID","Shape")     arcpy.MakeFeatureLayer_management(intersect,"intlyr")     arcpy.AddJoin_management("intlyr",arcpy.Describe("intlyr").OIDfieldName,findID,"IN_FID","KEEP_ALL")     segIDs = {}     featseqName = "explFindID.FEAT_SEQ"     idNewName = "intersect."+idName     for row in arcpy.SearchCursor("intlyr"):         idVal = row.getValue(idNewName)         featseqVal = row.getValue(featseqName)         segIDs[featseqVal] = []     for row in arcpy.SearchCursor("intlyr"):         idVal = row.getValue(idNewName)         featseqVal = row.getValue(featseqName)         segIDs[featseqVal].append(idVal)     segIDs2 = {}     for row in arcpy.SearchCursor("intlyr"):         idVal = row.getValue(idNewName)         segIDs2[idVal] = []     for x,y in segIDs.iteritems():         for segID in y:             segIDs2[segID].extend([k for k in y if k != segID])     for x,y in segIDs2.iteritems():         segIDs2 = list(set(y))     arcpy.RemoveJoin_management("intlyr",arcpy.Describe(findID).name)     if 'overlaps' not in [k.name for k in arcpy.ListFields(fc)]:         arcpy.AddField_management(fc,'overlaps',"TEXT")     if 'ovlpCount' not in [k.name for k in arcpy.ListFields(fc)]:         arcpy.AddField_management(fc,'ovlpCount',"SHORT")     urows = arcpy.UpdateCursor(fc)     for urow in urows:         idVal = urow.getValue(idName)         if segIDs2.get(idVal):             urow.overlaps = str(segIDs2[idVal]).strip('[]')             urow.ovlpCount = len(segIDs2[idVal])         urows.updateRow(urow) countOverlaps(fc,idName) def explodeOverlaps(fc,idName):     countOverlaps(fc,idName)     arcpy.AddField_management(fc,'expl',"SHORT")     urows = arcpy.UpdateCursor(fc,'"overlaps" IS NULL')     for urow in urows:         urow.expl = 1         urows.updateRow(urow)     i=1     lyr = arcpy.MakeFeatureLayer_management(fc)     while int(arcpy.GetCount_management(arcpy.SelectLayerByAttribute_management(lyr,"NEW_SELECTION",'"expl" IS NULL')).getOutput(0)) > 0:         ovList=[]         urows = arcpy.UpdateCursor(fc,'"expl" IS NULL','','','ovlpCount D')         for urow in urows:             ovVal = urow.overlaps             idVal = urow.getValue(idName)             intList = ovVal.replace(' ','').split(',')             for x in intList:                 intList[intList.index(x)] = int(x)             if idVal not in ovList:                 urow.expl = i             urows.updateRow(urow)             ovList.extend(intList)         i+=1 explodeOverlaps(fc,idName) but I'm receiving the following error message and I'm unable to resolve it. I'd really appreciate assistance in resolving the following as I'm still new  to python. Regards Peter Wilson ... View more

I've recently upgraded my Arc Hydro installation from 10.2.0.40 to 10.2.0.85 and my model stopped working where "Batch Subwatershed Delineation" was inserted. The output parameter for SubwatershedPoint  seems to be affected with the following error message: "Error 000735 Configuration for out_subwatershedpoint_features cannot be found in the configuration. : Value is required. The parameter is usually automatically populated with the string variable "SubwatershedPoint" and the output is saved to the output File Geodatabase Feature Dataset. I've tried installing previous versions as far back as 10.2.0.55 and the error message still persists. As soon as I re-install 10.2.0.40 the error message dissappears. Unnfortunately I need to be able to use the newer version due to changes that have been made. If anyone can assist in resoloving the following, it will truly be appreicated. I've attached a jpeg of the Batch Subwatershed Delineation with the error message for BatchSubwatershedPoint. Regards Peter Wilson ... View more