Find closest facility to multiple incidents (service)

View on GitHubSample viewer app

Find routes from several locations to the respective closest facility.

screenshot

Use case

Quickly and accurately determining the most efficient route between a location and a facility is a frequently encountered task. For example, a city's fire department may need to know which fire stations in the vicinity offer the quickest routes to multiple fires. Solving for the closest fire station to the fire's location using an impedance of "travel time" would provide this information.

How to use the sample

Click 'Solve Routes' to solve and display the route from each incident (fire) to the nearest facility (fire station). Click 'Reset' to clear the results.

How it works

  1. Create a ClosestFacilityTask using a URL from an online service.
  2. Get the default set of ClosestFacilityParameters from the task: closestFacilityTask.createDefaultParametersAsync().
  3. Create ServiceFeatureTables for both facilities incidents.
  4. Add all facilities to the task parameters: closestFacilityParameters.setFacilitiesWithFeatureTable(facilitiesFeatureTable, parameters).
  5. Add all incidents to the task parameters: closestFacilityParameters.setIncidentsWithFeatureTable(incidentsFeatureTable, parameters).
  6. Get ClosestFacilityResult by solving the task with the provided parameters: closestFacilityTask.solveClosestFacilityAsync(closestFacilityParameters).
  7. Find the closest facility for each incident by iterating over the list of Incidents.
  8. Display the route as a Graphic using the closestFacilityRoute.routeGeometry().

Relevant API

  • ClosestFacilityParameters
  • ClosestFacilityResult
  • ClosestFacilityRoute
  • ClosestFacilityTask
  • Facility
  • Graphic
  • GraphicsOverlay
  • Incident

Tags

incident, network analysis, route, search

Sample Code

FindClosestFacilityToMultipleIncidentsService.cppFindClosestFacilityToMultipleIncidentsService.cppFindClosestFacilityToMultipleIncidentsService.hFindClosestFacilityToMultipleIncidentsService.qml
Use dark colors for code blocksCopy
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
154
155
156
157
158
159
160
161
162
163
164
165
166
167
168
169
170
171
172
173
174
175
176
177
178
179
180
181
182
183
184
185
186
187
188
189
190
191
192
193
194
195
196
197
198
199
200
201
202
203
204
205
206
207
208
209
210
211
212
213
214
215
216
217
// [WriteFile Name=FindClosestFacilityToMultipleIncidentsService, Category=Routing]
// [Legal]
// Copyright 2019 Esri.
//
// Licensed under the Apache License, Version 2.0 (the "License");
// you may not use this file except in compliance with the License.
// You may obtain a copy of the License at
// http://www.apache.org/licenses/LICENSE-2.0
//
// Unless required by applicable law or agreed to in writing, software
// distributed under the License is distributed on an "AS IS" BASIS,
// WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied.
// See the License for the specific language governing permissions and
// limitations under the License.
// [Legal]

#ifdef PCH_BUILD
#include "pch.hpp"
#endif // PCH_BUILD

// sample headers
#include "FindClosestFacilityToMultipleIncidentsService.h"

// ArcGIS Maps SDK headers
#include "ClosestFacilityParameters.h"
#include "ClosestFacilityResult.h"
#include "ClosestFacilityRoute.h"
#include "ClosestFacilityTask.h"
#include "Envelope.h"
#include "Error.h"
#include "FeatureLayer.h"
#include "GeometryEngine.h"
#include "Graphic.h"
#include "GraphicListModel.h"
#include "GraphicsOverlay.h"
#include "GraphicsOverlayListModel.h"
#include "LayerListModel.h"
#include "Map.h"
#include "MapQuickView.h"
#include "MapTypes.h"
#include "PictureMarkerSymbol.h"
#include "Polyline.h"
#include "QueryParameters.h"
#include "ServiceFeatureTable.h"
#include "SimpleLineSymbol.h"
#include "SimpleRenderer.h"
#include "SymbolTypes.h"

// Qt headers
#include <QFuture>
#include <QUuid>

using namespace Esri::ArcGISRuntime;

FindClosestFacilityToMultipleIncidentsService::FindClosestFacilityToMultipleIncidentsService(QObject* parent /* = nullptr */):
  QObject(parent),
  m_map(new Map(BasemapStyle::ArcGISStreetsRelief, this)),
  m_task(new ClosestFacilityTask(QUrl("https://sampleserver6.arcgisonline.com/arcgis/rest/services/NetworkAnalysis/SanDiego/NAServer/ClosestFacility"), this)),
  m_resultsOverlay(new GraphicsOverlay(this))
{
  // enable busy indicator while loading
  m_busy = true;

  createSymbols();

  createFeatureLayers();

  connect(m_task, &ClosestFacilityTask::doneLoading, this, [this](const Error& e)
  {
    if (!e.isEmpty())
    {
      qDebug() << e.message();
      return;
    }
    setupRouting();
  });
  m_task->load();
}

FindClosestFacilityToMultipleIncidentsService::~FindClosestFacilityToMultipleIncidentsService() = default;

void FindClosestFacilityToMultipleIncidentsService::init()
{
  // Register the map view for QML
  qmlRegisterType<MapQuickView>("Esri.Samples", 1, 0, "MapView");
  qmlRegisterType<FindClosestFacilityToMultipleIncidentsService>("Esri.Samples", 1, 0, "FindClosestFacilityToMultipleIncidentsServiceSample");
}

MapQuickView* FindClosestFacilityToMultipleIncidentsService::mapView() const
{
  return m_mapView;
}

// Set the view (created in QML)
void FindClosestFacilityToMultipleIncidentsService::setMapView(MapQuickView* mapView)
{
  if (!mapView || mapView == m_mapView)
    return;

  m_mapView = mapView;
  m_mapView->setMap(m_map);

  emit mapViewChanged();
}

void FindClosestFacilityToMultipleIncidentsService::createSymbols()
{
  m_facilitySymbol = new PictureMarkerSymbol(QUrl("https://static.arcgis.com/images/Symbols/SafetyHealth/FireStation.png"), this);
  m_facilitySymbol->setWidth(30.0f);
  m_facilitySymbol->setHeight(30.0f);

  m_incidentSymbol = new PictureMarkerSymbol(QUrl("https://static.arcgis.com/images/Symbols/SafetyHealth/esriCrimeMarker_56_Gradient.png"), this);
  m_incidentSymbol->setWidth(30.0f);
  m_incidentSymbol->setHeight(30.0f);

  m_routeSymbol = new SimpleLineSymbol(SimpleLineSymbolStyle::Solid, Qt::blue, 2.0f, this);
}

void FindClosestFacilityToMultipleIncidentsService::createFeatureLayers()
{
  m_facilitiesFeatureTable = new ServiceFeatureTable(QUrl("https://services2.arcgis.com/ZQgQTuoyBrtmoGdP/ArcGIS/rest/services/San_Diego_Facilities/FeatureServer/0"), this);
  m_facilitiesFeatureLayer = new FeatureLayer(m_facilitiesFeatureTable, this);
  m_facilitiesFeatureLayer->setRenderer(new SimpleRenderer(m_facilitySymbol, this));

  m_incidentsFeatureTable = new ServiceFeatureTable(QUrl("https://services2.arcgis.com/ZQgQTuoyBrtmoGdP/ArcGIS/rest/services/San_Diego_Incidents/FeatureServer/0"), this);
  m_incidentsFeatureLayer = new FeatureLayer(m_incidentsFeatureTable, this);
  m_incidentsFeatureLayer->setRenderer(new SimpleRenderer(m_incidentSymbol, this));

  // connect to the doneLoading signal which calls the slot to set the view point geometry
  connect(m_facilitiesFeatureTable, &ServiceFeatureTable::doneLoading, this, &FindClosestFacilityToMultipleIncidentsService::setViewpointGeometry);
  connect(m_incidentsFeatureTable, &ServiceFeatureTable::doneLoading, this, &FindClosestFacilityToMultipleIncidentsService::setViewpointGeometry);

  m_facilitiesFeatureTable->load();
  m_incidentsFeatureTable->load();
}

void FindClosestFacilityToMultipleIncidentsService::setupRouting()
{
  m_task->createDefaultParametersAsync().then(this, [this](const ClosestFacilityParameters& defaultParameters)
  {
    QueryParameters params;
    params.setWhereClause("1=1");
    m_facilityParams = defaultParameters;
    m_facilityParams.setFacilitiesWithFeatureTable(m_facilitiesFeatureTable, params);
    m_facilityParams.setIncidentsWithFeatureTable(m_incidentsFeatureTable, params);

    m_busy = false;
    m_solveButtonEnabled = true;
    emit busyChanged();
    emit solveButtonChanged();
  });
}

void FindClosestFacilityToMultipleIncidentsService::setViewpointGeometry(const Error& e)
{
  if (!e.isEmpty())
  {
    qDebug() << e.message();
    return;
  }

  // proceed only if both layers have been loaded
  if ((m_facilitiesFeatureTable->loadStatus() != LoadStatus::Loaded) || (m_incidentsFeatureTable->loadStatus() != LoadStatus::Loaded))
    return;

  // return if set viewpoint future has already been created
  if (m_setViewpointFuture.isValid())
    return;

  m_mapView->map()->operationalLayers()->append(m_facilitiesFeatureLayer);
  m_mapView->map()->operationalLayers()->append(m_incidentsFeatureLayer);

  m_setViewpointFuture = m_mapView->setViewpointGeometryAsync(GeometryEngine::unionOf(m_facilitiesFeatureLayer->fullExtent(), m_incidentsFeatureLayer->fullExtent()), 20);
}

void FindClosestFacilityToMultipleIncidentsService::solveRoute()
{
  m_busy = true;
  m_solveButtonEnabled = false;
  emit busyChanged();
  emit solveButtonChanged();

  m_task->solveClosestFacilityAsync(m_facilityParams).then(this, [this]( const ClosestFacilityResult& closestFacilityResult)
  {
    if (closestFacilityResult.isEmpty())
    {
      qDebug() << "Empty result";
      return;
    }

    // finding the closest facility for each incident to create a route graphic between each pair
    for (int incidentIndex = 0; incidentIndex < m_incidentsFeatureTable->numberOfFeatures(); incidentIndex++)
    {
      const auto indexes = closestFacilityResult.rankedFacilityIndexes(incidentIndex);
      const int closestFacilityIndex = indexes.first();
      const ClosestFacilityRoute route = closestFacilityResult.route(closestFacilityIndex, incidentIndex);
      Graphic* m_routeGraphic = new Graphic(route.routeGeometry(), m_routeSymbol, this);

      m_resultsOverlay->graphics()->append(m_routeGraphic);
    }

    m_mapView->graphicsOverlays()->append(m_resultsOverlay);
    m_resetButtonEnabled = true;
    emit resetButtonChanged();
    m_busy = false;
    emit busyChanged();
  });
}

void FindClosestFacilityToMultipleIncidentsService::resetGO()
{
  m_mapView->graphicsOverlays()->clear();
  m_resetButtonEnabled = false;
  m_solveButtonEnabled = true;
  emit solveButtonChanged();
  emit resetButtonChanged();
}

Your browser is no longer supported. Please upgrade your browser for the best experience. See our browser deprecation post for more details.