Need some help with making multiple feature classes based on spatial relationships.

ColeNelson — Fri, 11 Feb 2022 17:19:46 GMT

Hey everyone, the ESRI Community has been really helpful and would not be as far as I am with this code with you. I do need help with this last bit though.

So the goal of this tool/code is to take the "Lines_Selected" line feature class and divide it up into two line feature classes based on their spatial relationship to my point feature classes.

The Lines_Selected feature class is created in the ForDeveloping.gdb and added to the map frame.

After this is done I want to divide it into two feature classes, based on the lines' spatial relationship to the point feature classes located in the NewHydrant.gdb. They should be called Kuna4WS_StreetsExtracted and Middleton4WS_StreetsExtracted.

However the code I currently have just copies the Lines_Selected feature class from the map frame and saves it in the NewHydrants.gdb, which isn't what I'm wanting. Here is the code that I am using, any help or suggestions is greatly appreciated, thank you!

import arcpy

import os

arcpy.env.overwriteOutput = True

arcpy.env.workspace = r"C:\UpdateDevelopmentPython\ForDeveloping\ForDeveloping.gdb"

EnterWaterSysID = arcpy.GetParameterAsText(0).split(";")

Hydrant_FC = "Hydrants"

WaterSystems = "WaterSystems"

Streetmap = r"C:\UpdateDevelopmentPython\Streetmap Premium\FGDB\StreetMap_Data\NorthAmerica.gdb\Routing\Routing_ND"

qry = "{0} IN ({1})".format(arcpy.AddFieldDelimiters(datasource= Hydrant_FC, field = 'WaterSysID'), ', '.join(EnterWaterSysID))

arcpy.SelectLayerByAttribute_management(in_layer_or_view=Hydrant_FC, where_clause=qry)

arcpy.SelectLayerByLocation_management(Hydrant_FC,"Within",WaterSystems,None,"SUBSET_SELECTION","INVERT")

arcpy.env.workspace = r"C:\UpdateDevelopmentPython\ForDeveloping\ForDeveloping.gdb"

Service_Area = arcpy.na.MakeServiceAreaAnalysisLayer(network_data_source=Streetmap,

layer_name="Service Area",

travel_mode="Driving Distance",

travel_direction="FROM_FACILITIES",

cutoffs=[0.189394],

time_of_day="",

time_zone="LOCAL_TIME_AT_LOCATIONS",

output_type="LINES",

polygon_detail="STANDARD",

geometry_at_overlaps="SPLIT",

geometry_at_cutoffs="RINGS",

polygon_trim_distance="100 Meters",

exclude_sources_from_polygon_generation=[],

accumulate_attributes=["Miles"],

ignore_invalid_locations="SKIP")[0]

Updated_Service_Area = arcpy.na.AddLocations(in_network_analysis_layer=Service_Area,

sub_layer="Facilities",

in_table=Hydrant_FC,

field_mappings="",

search_tolerance="5000 Meters",

sort_field="",

search_criteria=[["Routing_Streets", "SHAPE"],

["Routing_Streets_Override", "NONE"],

["Routing_ND_Junctions", "NONE"]],

match_type="MATCH_TO_CLOSEST",

append="APPEND",

snap_to_position_along_network="NO_SNAP",

snap_offset="5 Meters",

exclude_restricted_elements="EXCLUDE",

search_query=[])[0]

Solve_Succeeded = arcpy.na.Solve(in_network_analysis_layer=Updated_Service_Area,

ignore_invalids="SKIP",

terminate_on_solve_error="TERMINATE",

simplification_tolerance="",

overrides="")

Lines = arcpy.Select_analysis("Lines","Lines_Selected")

arcpy.SelectLayerByAttribute_management(Hydrant_FC, "CLEAR_SELECTION")

Lines_Extracted = r"C:\UpdateDevelopmentPython\ForDeveloping\ForDeveloping.gdb\Lines_Selected"

p = arcpy.mp.ArcGISProject("CURRENT")

m = p.listMaps("Map")[0]

m.addDataFromPath(Lines_Extracted)

gdb = r"C:\UpdateDevelopmentPython\ForDeveloping\NewHydrants.gdb"

HydrantLines = arcpy.ListFeatureClasses()

for fc in HydrantLines:

desc = arcpy.da.Describe(fc)

HydrantFeatures = os.path.join(gdb,desc["baseName"])

arcpy.MakeFeatureLayer_management(fc,HydrantFeatures)

arcpy.SelectLayerByLocation_management("Lines_Selected","WITHIN_A_DISTANCE",HydrantFeatures,"1000 Feet")

HydrantPoints = arcpy.ListFeatureClasses()

for fcs in HydrantPoints:

desc1 = arcpy.da.Describe(fcs)

arcpy.SpatialJoin_analysis("Lines_Selected",fcs, HydrantFeatures)

arcpy.Copy_management("Lines_Selected",HydrantFeatures)

Re: Need some help with making multiple feature classes based on spatial relationships.

DanPatterson — Fri, 11 Feb 2022 17:35:54 GMT

You might have missed the response in your original post

Help with using list layers and select by location... - Esri Community

Re: Need some help with making multiple feature classes based on spatial relationships.

ColeNelson — Fri, 11 Feb 2022 17:42:09 GMT

Lol, I did miss that, I'll give it a try.

topic Need some help with making multiple feature classes based on spatial relationships. in Python Questions

Need some help with making multiple feature classes based on spatial relationships.

Re: Need some help with making multiple feature classes based on spatial relationships.

Re: Need some help with making multiple feature classes based on spatial relationships.