GIS软件工程实习报告最短路径分析.docx

1、GIS软件工程实习报告最短路径分析AE开发之基于几何网络的最短路径分析1、实习目的本次实习目的在于熟练掌握ArcGIS Engine开发工具并能够通过C#语言在VS 2010开发环境中完成查询几何网络的最短路径分析的功能。2、实习时间2015年5月23日星期六3、实习内容3.1实验环境操作系统：Windows 2007二次开发平台：VS 2010开发环境、ArcGIS Desktop10.0、AE开发组件3.2实验任务完成基于几何网络分析的最短路径查询功能，即实现通过在几何网络地图中指定起始点，能够查询经过起始点的最短路线，并能够通过缩放功能在地图窗口居中显示。3.3实验步骤3.3.1新建项目

2、选择文件新建项目，如图选择项目类型中Visual C#，再选择Windows Application，记为“FindShortPath”，点击确定。3.3.2添加控件3.3.3控件绑定因为添加的控件只是单独存在，但是程序需要各控件间协同工作，因此要进行控件绑定。3.3.4创建几何网络1.在ArcCataLog中新建个人地理数据库“甘地的个人地理数据库”，然后新建要素数据集“road”，接着，鼠标右键选择要素数据集“road”，选择“导入要素类”，导入需要创建几何网络的要素类“主要公路”，输出要素类的名字改为“road”，如下图所示:2.鼠标右键选择要素数据集“road”，选择新建“几何网络”，

3、则出现“新建几何网络”对话框，作如下图所示的一系列设置，最终得到几何网络“road_Net”。 至此，得到几何网络“road_Net”，如下图所示:3.3.5代码实现最短路径分析1设置ToolStrip12添加代码private IActiveView m_ipActiveView; private IMap m_ipMap;/地图控件中地图 private IGraphicsContainer pGC;/图形对象 private bool clicked = false; int clickedcount = 0; private double m_dblPathCost = 0; priv

4、ate IGeometricNetwork m_ipGeometricNetwork; private IPointCollection m_ipPoints;/输入点集合 private IPointToEID m_ipPointToEID; private IEnumNetEID m_ipEnumNetEID_Junctions; private IEnumNetEID m_ipEnumNetEID_Edges; private IPolyline m_ipPolyline; private IMapControl3 mapctrlMainMap = null;private void F

5、orm1_Load(object sender, EventArgs e) /对象初始化 mapctrlMainMap = (IMapControl3)this.axMapControl1.Object; m_ipActiveView = axMapControl1.ActiveView; m_ipMap = m_ipActiveView.FocusMap; clicked = false; pGC = m_ipMap as IGraphicsContainer; private void Findpath_Click(object sender, EventArgs e) mapctrlMa

6、inMap.CurrentTool = null; /设置鼠标样式 axMapControl1.MousePointer = esriControlsMousePointer.esriPointerCrosshair; if (axMapControl1.LayerCount = 0) MessageBox.Show(请先加载几何网络数据！); return; m_ipActiveView = axMapControl1.ActiveView; m_ipMap = m_ipActiveView.FocusMap; clicked = false; pGC = m_ipMap as IGraph

7、icsContainer; ILayer ipLayer = m_ipMap.get_Layer(0); IFeatureLayer ipFeatureLayer = ipLayer as IFeatureLayer; IFeatureDataset ipFDS = ipFeatureLayer.FeatureClass.FeatureDataset; OpenFeatureDatasetNetwork(ipFDS); clicked = true; clickedcount = 0; pGC.DeleteAllElements(); private void SolvePath_Click(

8、object sender, EventArgs e) axMapControl1.MousePointer = esriControlsMousePointer.esriPointerDefault; if (SolvePathGan(Weight)/先解析路径 IPolyline ipPolyResult = PathPolyLine();/最后返回最短路径 clicked = false; if (ipPolyResult.IsEmpty) MessageBox.Show(没有路径可到!); else IRgbColor color = new RgbColorClass(); colo

9、r.Red = 255; color.Blue = 64; color.Green = 128; LineElementClass element = new LineElementClass(); ILineSymbol linesymbol = new SimpleLineSymbolClass(); linesymbol.Color = color as IColor; linesymbol.Width = 3; element.Geometry = m_ipPolyline; element.Symbol = linesymbol; pGC.AddElement(element, 0)

10、; m_ipActiveView.PartialRefresh(esriViewDrawPhase.esriViewGraphics, null, null); axMapControl1.FlashShape(element.Geometry); MessageBox.Show(路径经过 + m_ipEnumNetEID_Edges.Count + 条线 + rn + 经过节点数为 + m_ipEnumNetEID_Junctions.Count + rn + 路线长度为 + m_ipPolyline.Length.ToString(#.#) + rn, 几何路径信息); /路径缩放功能 p

11、rivate void SuoFang_Click(object sender, EventArgs e) if (m_ipPolyline = null) MessageBox.Show(当前没有执行路径查询！请确认！); else this.axMapControl1.Extent = m_ipPolyline.Envelope; #region 自定义路径查询函数 public void OpenFeatureDatasetNetwork(IFeatureDataset FeatureDataset) CloseWorkspace(); if (!InitializeNetworkAnd

12、Map(FeatureDataset) Console.WriteLine(打开network出错); /关闭工作空间 private void CloseWorkspace() m_ipGeometricNetwork = null; m_ipPoints = null; m_ipPointToEID = null; m_ipEnumNetEID_Junctions = null; m_ipEnumNetEID_Edges = null; m_ipPolyline = null; /初始化几何网络和地图 private bool InitializeNetworkAndMap(IFeatur

13、eDataset FeatureDataset) IFeatureClassContainer ipFeatureClassContainer; IFeatureClass ipFeatureClass; IGeoDataset ipGeoDataset; ILayer ipLayer; IFeatureLayer ipFeatureLayer; IEnvelope ipEnvelope, ipMaxEnvelope; double dblSearchTol; INetworkCollection ipNetworkCollection = FeatureDataset as INetwork

14、Collection; int count = ipNetworkCollection.GeometricNetworkCount; /获取第一个几何网络工作空间 m_ipGeometricNetwork = ipNetworkCollection.get_GeometricNetwork(0); INetwork ipNetwork = m_ipGeometricNetwork.Network; if (m_ipMap != null) ipFeatureClassContainer = m_ipGeometricNetwork as IFeatureClassContainer; coun

15、t = ipFeatureClassContainer.ClassCount; for (int i = 0; i count; i+) ipFeatureClass = ipFeatureClassContainer.get_Class(i); ipFeatureLayer = new FeatureLayerClass(); ipFeatureLayer.FeatureClass = ipFeatureClass; for (int j = 0; j m_ipMap.LayerCount; j+) if (m_ipMap.get_Layer(j).Name.ToUpper() = ipFe

16、atureLayer.Name.ToUpper() continue; count = m_ipMap.LayerCount; ipMaxEnvelope = new EnvelopeClass(); for (int i = 0; i dblHeight) dblSearchTol = dblWidth / 100; else dblSearchTol = dblHeight / 100; m_ipPointToEID.SnapTolerance = dblSearchTol; return true; /返回路径的几何体 public IPolyline PathPolyLine() IE

17、IDInfo ipEIDInfo; IGeometry ipGeometry; if (m_ipPolyline != null) return m_ipPolyline; m_ipPolyline = new PolylineClass(); IGeometryCollection ipNewGeometryColl = m_ipPolyline as IGeometryCollection;/引用传递 ISpatialReference ipSpatialReference = m_ipMap.SpatialReference; IEIDHelper ipEIDHelper = new E

18、IDHelper(); ipEIDHelper.GeometricNetwork = m_ipGeometricNetwork; ipEIDHelper.OutputSpatialReference = ipSpatialReference; ipEIDHelper.ReturnGeometries = true; IEnumEIDInfo ipEnumEIDInfo = ipEIDHelper.CreateEnumEIDInfo(m_ipEnumNetEID_Edges); int count = ipEnumEIDInfo.Count; ipEnumEIDInfo.Reset(); for

19、 (int i = 0; i count; i+) string info; info = new stringcount; ipEIDInfo = ipEnumEIDInfo.Next(); ipGeometry = ipEIDInfo.Geometry; ipNewGeometryColl.AddGeometryCollection(ipGeometry as IGeometryCollection); infoi = m_ipPolyline.Length.ToString(); return m_ipPolyline; public bool SolvePathGan(string W

20、eightName) try int intJunctionUserClassID; int intJunctionUserID; int intJunctionUserSubID; int intJunctionID; IPoint ipFoundJunctionPoint; ITraceFlowSolverGEN ipTraceFlowSolver = new TraceFlowSolver() as ITraceFlowSolverGEN; INetSolver ipNetSolver = ipTraceFlowSolver as INetSolver; if (m_ipGeometri

21、cNetwork = null) return false; INetwork ipNetwork = m_ipGeometricNetwork.Network; ipNetSolver.SourceNetwork = ipNetwork; INetElements ipNetElements = ipNetwork as INetElements; if (m_ipPoints = null) MessageBox.Show(请选择点!); return false; int intCount = m_ipPoints.PointCount;/这里的points有值吗？ /定义一个边线旗数组

22、 /*最近点*/ IJunctionFlag pJunctionFlagList = new JunctionFlagintCount; for (int i = 0; i intCount; i+) INetFlag ipNetFlag = new JunctionFlag() as INetFlag; IPoint ipJunctionPoint = m_ipPoints.get_Point(i); /查找输入点的最近的网络点 m_ipPointToEID.GetNearestJunction(ipJunctionPoint, out intJunctionID, out ipFoundJ

23、unctionPoint); ipNetElements.QueryIDs(intJunctionID, esriElementType.esriETJunction, out intJunctionUserClassID, out intJunctionUserID, out intJunctionUserSubID); ipNetFlag.UserClassID = intJunctionUserClassID; ipNetFlag.UserID = intJunctionUserID; ipNetFlag.UserSubID = intJunctionUserSubID; IJuncti

24、onFlag pTemp = (IJunctionFlag)(ipNetFlag as IJunctionFlag); pJunctionFlagListi = pTemp; ipTraceFlowSolver.PutJunctionOrigins(ref pJunctionFlagList); INetSchema ipNetSchema = ipNetwork as INetSchema; INetWeight ipNetWeight = ipNetSchema.get_WeightByName(WeightName); INetSolverWeights ipNetSolverWeigh

25、ts = ipTraceFlowSolver as INetSolverWeights; ipNetSolverWeights.FromToEdgeWeight = ipNetWeight;/开始边线的权重 ipNetSolverWeights.ToFromEdgeWeight = ipNetWeight;/终止边线的权重 object vaRes = new objectintCount - 1; /通过findpath得到边线和交汇点的集合 ipTraceFlowSolver.FindPath(esriFlowMethod.esriFMConnected, esriShortestPath

26、ObjFn.esriSPObjFnMinSum,out m_ipEnumNetEID_Junctions, out m_ipEnumNetEID_Edges, intCount - 1, ref vaRes); m_dblPathCost = 0; for (int i = 0; i vaRes.Length; i+) double m_Va = (double)vaResi; m_dblPathCost = m_dblPathCost + m_Va; m_ipPolyline = null; return true; catch (Exception ex) Console.WriteLin

27、e(ex.Message); return false; #endregion private void axMapControl1_OnMouseDown(object sender, IMapControlEvents2_OnMouseDownEvent e) if (clicked != true) return; IPoint ipNew; if (m_ipPoints = null) m_ipPoints = new Multipoint(); ipNew = m_ipActiveView.ScreenDisplay.DisplayTransformation.ToMapPoint(

28、e.x, e.y); object o = Type.Missing; m_ipPoints.AddPoint(ipNew, ref o, ref o); IElement element; ITextElement textelement = new TextElementClass(); element = textelement as IElement; clickedcount+; textelement.Text = clickedcount.ToString(); element.Geometry = m_ipActiveView.ScreenDisplay.DisplayTransformation.ToMapPoint(e.x, e.y); pGC.AddElement(element, 0); m_ipActiveView.PartialRefresh(esriViewDrawPhase.esriViewGraphics, null, null); 3.4实验结果在得到的窗口中加载“甘地的个人地理数据库”（在压缩包中有）中的几何网络“road_Net”，将几何网络地图加载到地图显示窗口中，如下图所示：点击图标添加起始点，效果如下图：然后点击图标，计算最短路径，效果如下图：最后点击图标得到最短路径的缩