Need script to select each feature, one at a time

johnodonnell

Hello all-

I'm having a tough time getting started on what should be a relatively simple script. I have a feature class consisting of sinks in a geometric network. I need my script to initially select the first record, perform some functions, finish those functions, then move to the next record, and repeat until the end. I've set up a search cursor to cycle through the records, but my problem is making a feature layer from the "current" record in the search cursor. Any input on how I could do this? Thanks, John

johnodonnell

I got it working! For anyone interested, the following script will take a point feature class representing sinks in a geometric network, select each one at a time, trace upstream from each, and update a field in an upstream feature class with a unique ID of the sink it traces to. Cheers!

import arcpy
arcpy.env.overwriteOutput=True
arcpy.env.workspace = "D:\FacID_Python\SO_CSS_Master_9x_GDB.gdb"

sinkFeatures = "D:\\FacID_Python\\SO_CSS_Master_9x_GDB.gdb\\SO_CSS_Master\\SManhole_Master"
dCatchBasin = "D:\\FacID_Python\\SO_CSS_Master_9x_GDB.gdb\\SO_CSS_Master\\DCatchBasin_Master"
geometricNetwork = "D:\\FacID_Python\\SO_CSS_Master_9x_GDB.gdb\\SO_CSS_Master\\SO_CSS_Master_Net"
mxd = arcpy.mapping.MapDocument("D:\FacID_Python\Working.mxd")

sinkList = []

rows = arcpy.SearchCursor(sinkFeatures)
row = rows.next()

while row:
    featureID = row.OBJECTID
    if featureID not in sinkList:
        sinkList.append(featureID)
    row = rows.next()
del rows

print sinkList

sinkCount = len(sinkList)

for each in range(sinkCount):
    print "%s of %s" % (each +1, sinkCount)
    sinkWhere = '"OBJECTID" = %s' % sinkList[each]

    if arcpy.Exists("TempSinkLayer"):
        arcpy.Delete_management("TempSinkLayer")


    arcpy.MakeFeatureLayer_management(sinkFeatures, "TempSinkLayer", sinkWhere)

    cursor = arcpy.SearchCursor("TempSinkLayer")
    sink = cursor.next()
    while sink:
        facID = sink.getValue("FACILITYID")
        print facID
        sink = cursor.next()
    del sink
    del cursor


    try:
        outGroupLayer = str("Trace" + facID)
        arcpy.SetFlowDirection_management(geometricNetwork, "WITH_DIGITIZED_DIRECTION")
        arcpy.TraceGeometricNetwork_management(geometricNetwork, outGroupLayer, "TempSinkLayer", "TRACE_UPSTREAM", "", "", "", "", "", "", "" ,"", "", "", "", "", "", "")
    except:
        print "Can't trace"

    try:
        dCatchLayer = str(facID + "dCatch")
        arcpy.SelectData_management(outGroupLayer, "DCatchBasin_Master")
        arcpy.MakeFeatureLayer_management("DCatchBasin_Master", dCatchLayer)
        result = int(arcpy.GetCount_management(dCatchLayer).getOutput(0))
        print result

        dCatchs = arcpy.UpdateCursor(dCatchLayer)
        dCatch = dCatchs.next()
        fieldID = "DnComboFacID"

        while dCatch:
            dCatch.setValue(fieldID, facID)
            dCatchs.updateRow(dCatch)
            dCatch = dCatchs.next()

        del dCatch
        del dCatchs

        arcpy.Delete_management(dCatchLayer)
        arcpy.Delete_management("DCatchBasin_Master")
        arcpy.SelectLayerByAttribute_management("DCatchBasin_Master", "CLEAR_SELECTION")
        ##for df in arcpy.mapping.ListDataFrames(mxd):
            ##for lyr in arcpy.mapping.ListLayers(mxd, "", df):
                ##if lyr.name == "Trace" + facID:
                    ##arcpy.mapping.RemoveLayer(df, lyr)
                ##else:
                    ##pass

    except:
        print "not working"

PS- the code works within PyScripter, but for some reason gets wonky results when run within ArcMap.