ArcGIS中最短路径的实现.docxVIP

ArcGIS中最短路径的实现

5、然后再右键新建一个GeometricNetwork，选择从已存在的图元中建立几何网络。

6、打开ArcMap，把刚才建立的“PersonalGeoDatabaseFeatureClass”添加到地图中，这样几何网络就建立好了。

这样我们就建立好一个几何网络了。我们现在要通过编程来实现最短路径，用到的接口主要有INetworkCollection，IGeometricNetwork，IPointToEID，ITraceFlowSolverGEN（它实现了ITraceFlowSolver的接口），INetSchema，IEIDHelper等。主要步骤如下：

1、获取几何网络工作空间

2、定义一个边线旗数组，把离点串最近的网络元素添加进数组

3、设置开始和结束边线的权重

4、进行路径分析

5、得到路径分析的结果

上次介绍了用几何网络实现的“最短路径”，这次用网络数据集实现真正的最短路径功能，跟上次一样，先处理下数据。

1、先打开ArcCatalog，连接到目标文件夹，假定该文件下有一个名为road的道路图层。

2、在road图层上右键新建一个网络数据集，并按照其默认设置直至完成。

3、打开该地图的工作空间，把刚才新建的网络数据集添加工作空间中。

4、在网络分析菜单中选择新建最近设施点。

这时在工作空间里，可以看到多了一个名为“ClosestFacility”的图层。它下面还有4个子图层，名字分别为“Facilities”，“Incidents”，“Barriers”，“Routes”。“Facilities”就是设施点图层，也就是目的点，“Incidents”的意思就是出发点，“Barriers”是障碍点，意思就是地图某条道路附近有一个障碍点，如果障碍点与道路距离在容限范围内，则表示此道路不通，“Routes”就是最终的结果。这样我们编程实现最短路径的思路就出现了：

1、添加出发点。

2、添加目的点。

3、生成最优路径，获取结果。

图片看不清楚？请点击这里查看原图（大图）。这里的添加出发点或者目的点，是往“Facilities”或“Incidents”图层上添加元素。获取结果也是从“Routes”中获取Polyline。往“Facilities”或“Incidents”图层上添加元素用到的主要方法是INALocator的QueryLocationByPoint函数，生成路径主要接口是INASolver和它的Solve方法。获取结果是按属性查找，因为“Routes”类其实就是一个图层类，只不过只是存在于内存。

1CMapControlDefaultm_map;

2IPointCollectionPtrm_ipPointCollection;

3

4ILayerPtripLayer=m_map.GetLayer(0);//网络数据集

5INALayerPtripNaLayer=ipLayer;

6if(NULL==ipNaLayer)

7{

8return;

9}

10

11INAContextPtripNaContext;

12HRESULThr=ipNaLayer-get_Context(ipNaContext);

13INAClassLoaderPtripNAClassLoader(CLSID_NAClassLoader);

14INALocatorPtripNALocator=NULL;

15hr=ipNaContext-get_Locator(ipNALocator);

16ipNALocator-put_SnapToleranceUnits(esriMeters);

17ipNALocator-put_SnapTolerance(200);

18ipNaContext;

19hr=ipNAClassLoader-putref_Locator(ipNALocator);

20

21INamedSetPtripNamedSet=NULL;

22ipNaContext-get_NAClasses(ipNamedSet);

23

24CStringszName=Facilities;

25BSTRbstrName=szName.Allo