Solve is not producing routing when more than 100 incidents are added
414
6
a month ago
by
by
by
by
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.
by
