GIS软件工程实习报告(最短路径分析)

AE开发之基于几何网络的最短路径分析1、实习目的本次实习目的在于熟练掌握ArcGIS Engine开发工具并能够通过C#语言在VS2010开发环境中完成查询几何网络的最短路径分析的功能。

2、实习时间2015年5月23日星期六3、实习内容3.1实验环境操作系统：Win dows 2007二次开发平台：VS 2010开发环境、ArcGIS Desktop10.0、AE开发组件3.2实验任务完成基于几何网络分析的最短路径查询功能，即实现通过在几何网络地图中指定起始点，能够查询经过起始点的最短路线，并能够通过缩放功能在地图窗口居中显示。

3.3实验步骤331新建项目选择文件新建项目，如图选择项目类型中Visual C#，再选择Win dows Application ，记为“ FindShortPath ”，点击确定。

Vi-5.u-ai BaskVtsual C?丄VTiual C#Wnd OWE Wet WPF.ES8®VtsuaJ C#思于宙ata育wirdow琵体泗户畀酝应用程序的：E目v^ ual aArcGJS CloudReporiEjng SharePoint Silwirliglhr WCFVkiuJ C#VrtualC# Fi ndShortP firth]ffiHCL）：d:\doicurnervt或yi翼Idl studig 201D\Pr-pjcct5■翩WE3S3T绘兴方宝宕粗⑷：FindShortPisthWindow敷迪用邕了WPF窗览H应用呈甲WindoMrs 觀男WPF用户邮注WorkflowWPF自走:込陛ft黑Vkud C#*曲^方妄缠目云口"⑹•••332添加控件od rjdShDrtPjlh.-^n - MEkvnc^r Viuul-ME.i IWQ QIQ^l -adJBl Hltl ISUMl BfllAl IKlOl£WFl 勺RilN 』■口|W) ■场；F3.3.3控件绑定因为添加的控件只是单独存在，但是程序需要各控件间协同工作，因此要进 行控件绑定.3 RndUiprtPfltijaiA df 』IVc^pdTiifll!J n n^lirtMCJWCB W HndrMa^flaiqrMr.-EiP ijl Errt 呼*rtyt b *枫 J 卫曲㈣ dl 5uc« Ji 吟 Mg 丄 口R EFE L E .W F^nnLa^gnr 住f 魁 FormLj^iu切疥E ・fE ・EXE白 Mu PHmlndflAstxi^frdl A-ln^crdNar din atu 直 ncJ 利X £4 札fnEjuDwfaui^BH ： JxTDC^irlJ-d IAmCilS LD J O FngriH Tunk/rCorfa 口IIS iCuO Enitii magag ■加jiSJMaipCArril'Oll(ilEzlid J334创建几何网络1.在ArcCataLog 中新建个人地理数据库一“甘地的个人地理数据库”，然后 新建要素数据集“ road ”，接着，鼠标右键选择要素数据集“ road ”，选择“导入 要素类”，导入需要创建几何网络的要素类 “主要公路”，输出要素类的名字改为 “road ”，如下图所示：2.鼠标右键选择要素数据集“ road ”，选择新建“几何网络”，贝U 出现“新建 几何网络”对话框，作如下图所示的一系列设置，最终得到几何网络“road Net ”。

云南.广西1斷龍如地酸枣列恥(j3ca_czcinrjpts.shp 应中 3B#.5hp 遛甲国嘶界月hp.画申西刻注・hp凹主甕耐hp 曰主豪旦蛰曲p且示沁主妥埜豁.出口已到也诰过歸。

k 庇NvAC#+.1♦)IftittaSC. \1T X «TS \JL^L ui ttziHy品J 轶甲血放攜县畀* L»b*r E 投窮) '申皿■:彌芒日*椿出事;fr 晳湖南贵州:U 凹福建 广东夕筋 香港.竽;e72fi学丘□■iSi 认工也目录“於咲星旳〒也至此，得到几何网络“ road_Net”，如下图所示：0 riwiMTupjnud 1- ArcMap - JUrlrha宜血聊H E U 凹柱im 补如酸icq <<a^®：ui •土:亦口 ■口㈣"MKi335代码实现最短路径分析① 设置 ToolStriplT 甘坦一实习:Xvf c n②添加代码private lActiveView m 」pActiveView;private IMap m_ipMap;〃 地图控件中地图 private IGraphicsCo ntai ner pGC;〃 图形对象private bool clicked = false; int clickedco unt = 0;private double m_dblPathCost = 0;private IGeometricNetwork m_ipGeometricNetwork; private IPoi ntCollection m 」pPoi nts;〃 输入点集合private IPoi ntToEID m_ipPo in tToEID;private lEnum NetEID m_ipE num NetEID_J un ctio ns;星迪吐出 ^-^44.C4K卜厂2- 4口己 J 吞 才“ 1 M.ata.uu" £ IH 3FS SB S« wd ：g 种-"P ：4ESS亀.耳 flO ■£» 驚 2 鬥佝|叩祁 ~ -|8U- ■■- irtf-上 *-a BK鬲 rcd»d,FlM_lunciiMc*w g gd■ 0 flttliKprivate lEnum NetEID m_ipE num NetEID_Edges;private IPolyli ne m_ipPolyli ne;private IMapC on trol3 map ctrlMai nMap = null;private void Form1_Load(object sen der, Eve ntArgs e) {// 对象初始化mapctrlMai nMap(IM apCon trol3)this.axM apCon troll.Object;m」pActiveView = axMapCo ntroll.ActiveView;m」pMap = m」pActiveView.FocusMap;clicked = false;pGC = m」pMap as IGraphicsC on tai ner;}private void Findpath_Click(object sender, EventArgs e) {mapctrlMa in Map.Curre ntTool = n ull;// 设置鼠标样式axMapCo ntroll.MousePoi nteresriC on trolsMousePoi nter.esriPoi nterCrosshair;if (axM apCo yerCou nt == 0)return; }m 」pActiveView = axMapCo ntroll.ActiveView; m 」pMap = m 」pActiveView.FocusMap; clicked = false;pGC = m 」pMap as IGraphicsC on tai ner; ILayer ipLayer = m_ipMap.get_Layer(O); IFeatureLayer ipFeatureLayer = ipLayer as IFeatureLayer;IFeatureDatasetipFDSipFeatureLayer.FeatureClass.FeatureDataset;Ope nF eatureDatasetNetwork(ipFDS); clicked = true; clickedco unt = 0;pGC.DeleteAIIEIeme nts();private void SolvePath_Click(object sender, EventArgs e) {axMapCo ntrol1.MousePoi nteresriCo ntrolsMousePoi nter.esriPoi nterDefault;MessageBox.Show (”请先加载几何网络数据！"）；if (SolvePathGa n("Weight"))〃{先解析路径IPolyline ipPolyResult =PathPolyLine();〃最短路径clicked = false;if (i pPolyResult.lsEmpty){MessageBox.Show(”没有路径可到!!");}else{IRgbColor color = new RgbColorClass(); color.Red = 255;color.Blue = 64;color.Green = 128;Lin eEleme ntClass eleme ntLin eEleme ntClass();ILin eSymbol lin esymbolSimpleLi neSymbolClass();li nesymbol.Color = color as IColor;lin esymbol.Width = 3;eleme nt.Geometry = m_ipPolyli ne;最后返回newneweleme nt.Symbol = lin esymbol;pGC.AddEleme nt(eleme nt, 0);m_ipActiveView.PartialRefresh(esriViewDrawPhase.esriViewGraphics,null, null);axMapC on trol1.FlashShape(eleme nt.Geometry);MessageBox.Show(”路径经过”+ m」pEnumNetEID_Edges.Count + "条线"+ "\r\n" + "经过节点数为"+m」pEnumNetEID_Junctions.Count + "\r\n" + "路线长度为" + m」pPolyline.Length.ToString("#######.##") + "\r\n"," 几何路径信息");}// 路径缩放功能private void SuoFa ng_Click(object sen der, Eve ntArgs e)if (m_ipPolyli ne == nu II){MessageBox.Show(”当前没有执行路径查询！！！请确认! );else {this.axMapC on trolI.Exte nt = m」pPolyli ne.E nv elope;}}#region 自定义路径查询函数public void OpenF eatureDatasetNetwork(IFeatureDataset FeatureDataset){CloseWorkspace();if (!I nitializeNetworkA ndMap (FeatureDataset))Console.WriteLine(”打开network 出错");}// 关闭工作空间private void CloseWorkspace(){m_ipGeometricNetwork = nu II;m」pPoi nts = n ull;m_ipPo in tToEID = null;m_ipE num NetEID_J un ctio ns = nu II;m_ipE num NetEID_Edges = nu II;m_ipPolyli ne = n ull;// 初始化几何网络和地图bool In itializeNetworkA ndMap(IFeatureDatasetFeatureDataset)IFeatureClassC ontainer ipFeatureClassC ontainer; IFeatureClass ipFeatureClass; IGeoDataset ipGeoDataset; ILayer ipLayer;IFeatureLayer ipFeatureLayer;IEn velope ipE nvelope, ipMaxE nv elope; double dblSearchTol; INetworkCollect ionipNetworkCollecti on = FeatureDatasetas INetworkCollectio n;int count = ipNetworkCollectio n.GeometricNetworkCo unt;m_ipGeometricNetworkipNetworkCollecti on. get_GeometricNetwork(O);INetwork ipNetwork = m_work; if (m 」pMap != n ull)ipFeatureClassC ontainer = m_ipGeometricNetwork asIFeatureClassC ontainer;private // 获取第一个几何网络工作空间count = ipFeatureClassC ontain er.ClassCo unt;for (int i = 0; i < count; i++){ipFeatureClassipFeatureClassC on tai ner.get_Class(i);ipFeatureLayer = new FeatureLayerClass();ipFeatureLayer.FeatureClass = ipFeatureClass;for (int j = 0; j < m_yerCo un t; j++){if (m_ipMap.get_Layer(j).Name.ToUpper()== .ToUpper()){con ti nue;}}}}count = m_yerC ount;ipMaxE nv elope = new Env elopeClass();for (int i = 0; i < count; i++) ipLayer = m_ipMap.get_Layer(i);ipFeatureLayer = ipLayer as IFeatureLayer; ipGeoDataset =ipFeatureLayer as IGeoDataset; ipE nvelope = ipGeoDataset.Exte nt;ipMaxE nv elope.U nio n(ipE nv elope);}m_ipPo in tToEID = new Poi ntToEIDCIass(); m_ipPointToEID.SourceMap = m」pMap;m_ipPo in tToEID.GeometricNetwork m_ipGeometricNetwork;double dbIWidth = ipMaxE nv elope.Width;double dblHeight = ipMaxE nv elope.Height;if (dblWidth > dblHeight)dblSearchTol = dblWidth / 100;elsedblSearchTol = dblHeight / 100;m_ipPo in tToEID.S napTolera nee = dblSearchTol; return true; }// 返回路径的几何体public IPolyline PathPolyLine()IEIDI nfo ipEIDI nfo;IGeometry ipGeometry;if (m_ipPolyli ne != n ull) return m_ipPolyli ne;m_ipPolyli ne = new Polyli neClass();as IGeometryCollectio n ipNewGeometryColl = m_ipPolyli ne IGeometryCollectio n;〃引用传递ISpatialRefere nee ipSpatialRefere neem_ipMap.SpatialRefere nee;IEIDHelper ipEIDHelper = new EIDHelper();ipEIDHelper.GeometricNetwork = m_ipGeometricNetwork;ipEIDHelper.OutputSpatialRefere neeipSpatialRefere nee;ipEIDHelper.Retur nGeometries = true;lEnu mEIDI nfo ipE numEIDI nfoipEIDHelper.CreateE nu mEIDI nfo(m」pEnum NetEID_Edges);int count = ipE nu mEID In fo.Co unt;ipE nu mEIDI nfo.Reset();for (int i = 0; i < count; i++){stri ng[] info;info = new stri ng[co unt];ipEIDI nfo = ipE numEIDInfo.Next();ipGeometry = ipEIDIipNewGeometryColl.AddGeometryCollecti onasnfo.Geometry;(ipGeometry IGeometryCollectio n);in fo[i] = m_i pPolyli ne丄en gth.ToString();}return m_ipPolyli ne;}public bool SolvePathGa n(stri ngWeightName) {try{int intJun cti on UserClassID;int intJun cti on UserID;int intJun cti on UserSubID;int intJun cti on ID;IPo int ipF oundJunctionPoint;ITraceFlowSolverGENipTraceFlowSolver new TraceFlowSolver() as ITraceFlowSolverGEN;if (m_ipGeometricNetwork == nu II) { return false;INetwork ipNetwork = m_work; ipNetSolver.SourceNetwork = ipNetwork;Jun cti onF lag[]pJun cti onF lagList = JunctionF lag[i ntCo un t];for (i nt i = 0; i < in tCo unt; i++) { INetFlagipNetFlag = new JunctionFlag() INetFlag;INetSolver ipNetSolver ipTraceFlowSolver asINetSolver;INetEleme ntsipNetEleme nts = ipNetwork asINetEleme nts; if (m 」pPoi nts == nu II) { MessageBox.Show (” 请选择点!!"); return false; } int points 有值吗？ Illiin tCo unt = m 」pPoi nts.Poi ntCou nt;〃 这定义一个边线旗数组最近点*************** IIIIIIIIIInewasm 」pPoi nts.get_Poi nt(i);//m 」pPoi ntToEID.GetNearestJ un ctio n(ipJ un ctio nPoi nt, out intJun cti on ID, out ipF oundJunctionPoin t);ipNetEleme nts.QueryIDs(i ntJ un ctio nID,esriEleme ntType.esriETJ un cti on,out intJunction UserClassID,intJun cti on UserID, out intJunction UserSubID);erClassIDintJun cti on UserClassID;erID = in tJu nctio nUserlD; erSublD = in tJu nctio nUserSublD; IJun cti onF lagpTemp= (IJ un ctio nF lag)(ipNetFlagas IJunctionF lag);pJun cti onF lagList[i] = pTemp; }ipTraceFlowSolver.PutJ un ctio nOrigi ns(refpJunctionF lagList);INetSchema ipNetSchema = ipNetwork as INetSchema; INetWeightipNetWeightipNetSchema.get_WeightByName(WeightName);IPoi ntipJunctionPoint查找输入点的最近的网络点outINetSolverWeights ipNetSolverWeights ipTraceFlowSolver as INetSolverWeights;ipNetSolverWeights.FromToEdgeWeightipNetWeight;// 开始边线的权重ipNetSolverWeights.ToFromEdgeWeightipNetWeight;// 终止边线的权重object]] vaRes = new object[intCount - 1];// 通过findpath得到边线和交汇点的集合ipTraceFlowSolver.Fi ndPath(esriFlowMethod.esriFMC onn ected, esriShortestPathObjF n.esriSPObjF nMin Sum,out m」pEnum NetEID_J un ctio ns, out m_ipE num NetEID_Edges, i ntCou nt - 1, ref vaRes);m_dblPathCost = 0;for (int i = 0; i < vaRes.Le ngth; i++){double m_Va = (double)vaRes[i];m_dblPathCost = m_dblPathCost + m_Va;}m_ipPolyli ne = n ull;return true;catch (Excepti on ex) {Console.WriteLine(ex.Message); return false;}#en dregi onprivate void axMapCo ntrol1_O nM ouseDow n(o bjectsen der, IMap Co ntrolEve nts2_ On MouseDow nEve nt e)if (clicked != true)return;IPoi nt ipNew;if (m」pPoi nts == n ull){m」pPoi nts = new Multipoi nt();ipNewm」pActiveView.Scree nDisplay.DisplayTra nsformatio n.ToMapPoi nt(e.x,e.y)；object o = Type.Miss ing;m_ipPoi nts.AddPoi nt(ipNew, ref o, ref o);IEleme nt eleme nt;ITextEleme nt texteleme nt = new TextEleme ntClass();eleme nt = texteleme nt as IEleme nt;clickedco un t++;texteleme nt.Text = clickedcou nt.ToStri ng();eleme nt.Geometrym」pActiveView.Scree nDisplay.DisplayTra nsformatio n.ToMapPoi nt(e.x, e.y);pGC.AddEleme nt(eleme nt, 0);m_ipActiveView.PartialRefresh(esriViewDrawPhase.esriViewGraphics, null, null);3.4实验结果在得到的窗口中加载“甘地的个人地理数据库”（在压缩包中有）中的几何网络“ road_Net ”，将几何网络地图加载到地图显示窗口中，如下图所示:即1 —3融圉■Pha宜按豳：|申gdShapefiledDrUE.Nilsfi±服芬器颓攥集•显示粪3? |口1 gyppnrt显打开三］取消亘亘I亘LvnEvent e)点击图标添加起始点，效果如下图最后点击图标 得到最短路径的缩放，效果如下图:这样就完成了最短路径的查询分析和显示的功能4、实习总结在整个过程中所出现的问题较多，体现出个人基础的不足，程序的出错有时 候出现在标点符号以及英文字母的大小写上面等等问题， 这些都说明了不熟练的 问题。