Iterate different stops for different routes
Hi Hazel,
Are you by any chance familiar with python? You could use excel to create the permutations or even have them written down in notepad, but taking these approaches is going to result in you having to manually input this information into your model.
I would highly suggest using python to accomplish this task. Within python creating the permutations is a one line call. For example,
I'm not sure how well you understand python, but I have included code below that I was able to utilize to create the results you're seeking. You can review this if you want to utilize this workflow to accomplish your task. I pulled the example from the page listed below and modified it.
Make Route Layer
http://desktop.arcgis.com/en/desktop/latest/tools/network-analyst-toolbox/make-route-layer.htm
# Name: MakeRouteLayer_Workflow.py
# Description: Find a best route to visit the stop locations and save the
# route to a layer file. The stop locations are geocoded from a
# text file containing the addresses.
# Requirements: Network Analyst Extension
# Import system modules
import arcpy
import itertools
import os
from arcpy import env
try:
dataPath = os.path.join(os.path.dirname(__file__), "data")
# Check out the Network Analyst extension license
arcpy.CheckOutExtension("Network")
# Set environment settings
env.workspace = os.path.join(dataPath, "SanFrancisco.gdb")
env.overwriteOutput = True
# Set local variables
inNetworkDataset = "Transportation/Streets_ND"
outNALayerName = "BestRoute"
impedanceAttribute = "TravelTime"
outStops = "Analysis/Stores"
outLayerFile = os.path.join(dataPath, "output", outNALayerName + ".lyr")
if arcpy.Exists(os.path.dirname(outLayerFile)):
arcpy.management.Delete(os.path.dirname(outLayerFile))
os.makedirs(os.path.dirname(outLayerFile))
# Create a new Route layer. For this scenario, the default value for all the
# remaining parameters statisfies the analysis requirements
outNALayer = arcpy.na.MakeRouteLayer(inNetworkDataset, outNALayerName, impedanceAttribute)
# Get the layer object from the result object. The route layer can now be
# referenced using the layer object.
outNALayer = outNALayer.getOutput(0)
# Get the names of all the sublayers within the route layer.
subLayerNames = arcpy.na.GetNAClassNames(outNALayer)
# Stores the layer names that we will use later
stopsLayerName = subLayerNames["Stops"]
# Load the geocoded address locations as stops mapping the address field from
# geocoded stop features as Name property using field mappings.
fieldMappings = arcpy.na.NAClassFieldMappings(outNALayer, stopsLayerName)
fieldMappings["Name"].mappedFieldName = "NAME"
# Get a list of all of the OIDs for the stops we need to create a route against
oidfieldname = arcpy.Describe(outStops).oidfieldname
with arcpy.da.SearchCursor(outStops, ["OID@"], "{} < {}".format(*[oidfieldname, 6])) as cursor:
oids = [row[0] for row in cursor]
# Get all of the permutations for all of the stops except the start/end
permutations = list(itertools.permutations(oids[1:-1], len(oids)-2))
stopsLayer = arcpy.management.MakeFeatureLayer(outStops, "LoadOrder")
for n,stops in enumerate(permutations):
# Add the start stop
arcpy.management.SelectLayerByAttribute(stopsLayer, "NEW_SELECTION", "{} = {}".format(*[oidfieldname, oids[0]]))
arcpy.na.AddLocations(outNALayer, stopsLayerName, stopsLayer, fieldMappings, "", exclude_restricted_elements="EXCLUDE",
append="CLEAR")
# Add the in-between stops in the permutation order
for i, stop in enumerate(stops):
arcpy.management.SelectLayerByAttribute(stopsLayer, "NEW_SELECTION", "{} = {}".format(*[oidfieldname, stop]))
arcpy.na.AddLocations(outNALayer, stopsLayerName, stopsLayer, fieldMappings, "", exclude_restricted_elements="EXCLUDE",
append="APPEND")
# Add the end stop
arcpy.management.SelectLayerByAttribute(stopsLayer, "NEW_SELECTION", "{} = {}".format(*[oidfieldname, oids[-1]]))
arcpy.na.AddLocations(outNALayer, stopsLayerName, stopsLayer, fieldMappings, "", exclude_restricted_elements="EXCLUDE",
append="APPEND")
# Solve the route layer, ignore any invalid locations such as those that
# can not be geocoded
arcpy.na.Solve(outNALayer,"SKIP")
# Save the solved route layer as a layer file on disk with relative paths
arcpy.management.SaveToLayerFile(outNALayer,outLayerFile.replace(".lyr", "_{0}.lyr".format(str(n+1).zfill(len(
permutations)))),"RELATIVE")
print("Script completed successfully")
except Exception as e:
# If an error occurred, print line number and error message
import traceback, sys
tb = sys.exc_info()[2]
print("An error occured on line %i" % tb.tb_lineno)
print(str(e))
As a side note, if you wanted you could take lines 52-83 from my code, convert it into a script tool and then use this within a Model to implement the advanced logic needed here to execute the task. This iteration logic would not be natively possible to implement within ModelBuilder.
