So, there are a couple of things to mention, not so much the question but issues that I see with the script itself.
- Avoid nesting loops. In terms of performance, nested loops can be processing killers if they are used within themselves.
- All cursor functions have SQL clauses that can be used, so having separate select by attribute query layer is unnecessary.
- In my personal opinion, never set functions within the main portion of a script. I typically set them above the main processes and below the imports.
I would also recommend checking out arcpy geometry to get a better understanding of how to utilize feature geometries.
import arcpy
from arcpy import ListFields
from arcpy.da import SearchCursor as Searching , UpdateCursor as Updating
# Retrieves records of an input layer/feature and returns a dictionary of values
def GetRecords( Layer , ["OID@"] + Fields + ['SHAPE@'] , SQLClause ):
SearchLayer = Searching( Layer , Fields )
if SQLClause is not None: SearchLayer = Searching( Layer , Fields , SQLClause )
return { row[0] : row[ 1:] for row in SearchLayer }
# Returns a list of field names
def GetFieldNames( Layer , ExcludeFields ):
fieldnames = [ field.name for field in ListFields( Layer ) ]
if type( ExcludeFields ) is list: fieldnames = list( set( fieldnames ).difference( set( ExcludeFields ) ) )
return fieldnames
# Checks the distance between two geometries
def CheckProximity( InputFeatureRecords , DistanceFeatureRecords , SetDistance ):
WithinSetLimits = { }
for IF_id , IF_Attributes in InputFeatureRecords.items():
for DF_id , DF_Attributes in DistanceFeatureRecords.items():
Point = IF_Attributes[ 0 ]
Line = DF_Attributes[ 0 ]
if Point.distanceTo( Line ) <= SetDistance: WithinSetLimits[ IF_id ] = DF_Attributes[ 0 ]
else: WithinSetLimits[ IF_id ] = "Not Within 100 Feet of Stream or River"
return WithinSetLimits
Poles = '<PoleLayer>'
Streams = '<StreamLayer>'
PoleFields = GetFieldNames( infeature , '<[ list of field names to exclude ]>' )
StreamFields = [ "Stream_ID" , "StreamPermanence" ]
PoleRecords = GetRecords( infeature , ["OID@"] + Fields + ['SHAPE@'] )
StreamRecords = GetRecords( infeature , ["OID@"] + Fields + ['SHAPE@'] )
CheckProximity( PoleRecords , StreamRecords , 100 )
The sample above isn't a solution but I just thought to provide something useful to give you a rough idea.
