Solve is not producing routing when more than 100 incidents are added
Ah gotcha, that makes sense.
So I have reworked the code to more reflect the the nax documentation provided. I'm getting to the solve and now it is failing:
[0, 'Executing network analysis service Job Id: "jef9f9f2c394d4fb5881f9fafabf52911".'],
[0, 'Submitted.'],
[0, 'Executing...'],
[0, 'No "Facilities" found for "Location 15" in "Incidents".'],
[0, 'No "Facilities" found for "Location 14" in "Incidents".'],
... # Starts at 15 goes down to 1 then from 16 to 100
[0, 'No "Facilities" found for "Location 100" in "Incidents".'],
[0, 'ERROR 030212: Solve did not find a solution.'],
[0, 'No solution found.'],
[0, 'Failed to execute (FindClosestFacilities).'],
[0, 'Failed.']
Updated Code:
# Get Current Working Directory
path = os.getcwd()
print('Current Working Directory: ', path)
# Build more paths
na_gdb = os.path.join(path, 'Network Adequacy.gdb')
arcpy.env.workspace = na_gdb
closest_facility_path = os.path.join(na_gdb, 'Closest Facility')
routes = os.path.join(closest_facility_path, 'Routes')
incidents = os.path.join(closest_facility_path, 'Incidents')
facilities = os.path.join(closest_facility_path, 'Facilities')
# Build intermediary paths
mce_incidents_xy = os.path.join(na_gdb, 'Incidents_XYTableToPoint')
mce_facilities_xy = os.path.join(na_gdb, 'Facilities_XYTableToPoint')
# Build output paths
staging = os.path.join(path, 'Staging Data')
incidents_routing = os.path.join(staging, 'Incidents_Routing_ExportTable.csv')
facilities_routing = os.path.join(staging, 'Facilities_Routing_ExportTable.csv')
all_routing = os.path.join(staging, 'Routes_ExportTable.csv')
# Delete previous versions of files
print('Deleting old directories..')
delete_list = [na_gdb, incidents_routing, facilities_routing, all_routing]
for delete in delete_list:
if arcpy.Exists(delete):
arcpy.Delete_management(delete)
# Create geodatabase
print('Create fresh gdb')
arcpy.management.CreateFileGDB(
path,
'Network Adequacy.gdb',
)
#%%
for spec in spec_list:
for mce in mce_list:
# Build output paths
mce_incidents_csv = os.path.join(path, f'01 Incidents Output/{mce}_Incidents.csv')
mce_facilities_csv = os.path.join(path, f'02 Facilities Output/{mce}_{spec}_Facilities.csv')
mce_routing_csv = os.path.join(path, f'03 Routing Output/{mce}_{spec}_Routes.csv')
# Read incident and facilities csv
mce_incidents_df = pd.read_csv(mce_incidents_csv, dtype=str)
mce_facilities_df = pd.read_csv(mce_facilities_csv, dtype=str)
# Calculate routes
print(f'routing {len(mce_incidents_df)} incidents to {len(mce_facilities_df)} facilities...')
if len(mce_facilities_df) == 0:
# If no facilities:
print(f'No facilities for {mce} - {spec}')
incident_column = mce_incidents_df['Incident_Match_ID']
# Build report for routing
mce_routing_df = pd.DataFrame({
'MCE_Name': mce,
'Specialty': spec,
'Incident_Match_ID': incident_column,
'Facility_Match_ID': np.nan,
f'Driving_Minutes': np.nan,
f'Driving_Miles': np.nan
})
print('Export No Facilities Report csv')
mce_routing_df.to_csv(mce_routing_csv, index=False)
else:
arcpy.management.XYTableToPoint(
path + f'/01 Incidents Output/{mce}_Incidents.csv',
mce_incidents_xy,
'X',
'Y',
None,
coordinate_system
)
# Add facility data to geodatabase
arcpy.management.XYTableToPoint(
path + f'/02 Facilities Output/{mce}_{spec}_Facilities.csv',
mce_facilities_xy,
'X',
'Y',
None,
coordinate_system
)
# Instantiate a ClosestFacility solver object
closest_facility = arcpy.nax.ClosestFacility(portal)
# Set properties
closest_facility.travelMode = travel_mode
closest_facility.timeUnits = arcpy.nax.TimeUnits.Minutes
closest_facility.defaultImpedanceCutoff = 15
closest_facility.defaultTargetFacilityCount = 1
closest_facility.routeShapeType = arcpy.nax.RouteShapeType.TrueShapeWithMeasures
# Load inputs
closest_facility.load(arcpy.nax.ClosestFacilityInputDataType.Facilities, mce_facilities_xy)
closest_facility.load(arcpy.nax.ClosestFacilityInputDataType.Incidents, mce_incidents_xy)
# Solve the analysis
result = closest_facility.solve()
# Export the results to a feature class
if result.solveSucceeded:
result.export(arcpy.nax.ClosestFacilityOutputDataType.Routes, all_routing)
else:
print("Solve failed")
print(result.solverMessages(arcpy.nax.MessageSeverity.All))
# Export Routes table
print('Exporting Routes')
arcpy.conversion.ExportTable(
routes,
all_routing,
)
# Export Incidents table
print('Exporting Incidents Table')
arcpy.conversion.ExportTable(
incidents,
incidents_routing,
)
# Export Facilities table
print('Exporting Facilities Table')
arcpy.conversion.ExportTable(
facilities,
facilities_routing,
)
# Read CSV into dataframes
print('Building Routing Table')
routes_df = pd.read_csv(all_routing)
# Set variables for columns
incident_column = routes_df['IncidentID']
facility_column = routes_df['FacilityID']
route_column = routes_df['Name']
minutes_column = routes_df['Total_TravelTime']
miles_column = routes_df['Total_Miles']
# Create final dataframe
df = pd.DataFrame({
'MCE_Name': mce,
'Specialty': spec,
'Incident_Match_ID': incident_column,
'Facility_Match_ID': facility_column,
'Route': route_column,
'Driving_Minutes': minutes_column,
'Driving_Miles': miles_column
})
# Save as csv
print(f'Writing {mce} - {spec} to csv..')
df.to_csv(mce_routing_csv, index=False)
So from the error message it seems like the incidents (100) and the facilities (1) are being added, but the routes are failing to generate.
I went and tried to replicate this in the Pro GUI and got the same results, I then tried in the GUI with Output Geometry to Straight Line and it output that fine so that is leading me to believe the Incidents are not being placed on the network. Is there a tolerance setting or should I be doing something else to generate the distance a point is from joining the road network? I remember messing with this about 10 years ago and specifically remember making a line feature class of the closest road network vertex with each facility and incident so that the tolerance could be customized.
Just closing the loop here, have worked extensively trying to get this to work, there were existing licensing issues that needed to be resolved, a memory upgrade and republishing all the routing geoprocessing services and slowly stripping back other things that are interacting with it so even just using the NA REST endpoints with my previous data converted to JSON still does produce an error so there is something fundamentally wrong with the install or the server and is probably not replicable by anyone.
However, if anyone comes across this in the future and is having trouble with a StreetMap Premium producing routes:
1. Turn on debug logging
2. Run Analysis in Pro (like an end user)
3. If still failing and nothing in the logs; try in Portal with the same data but uploaded.
4. If still failing and nothing in the logs; access the developer tools in your browser and you can extract the json of the data from the network tool under submitJob and you can just place that in the html REST endpoint located here: https://machine.domain.com/webAdapter/rest/services/Routing/NetworkAnalysis/GPServer/FindClosestFaci...
Then at least you know nothing else is getting in the way and the entire thing is borked.
