ArcGIS网络分析最短路径分析源代码（VB6.0）精品资料.doc

《ArcGIS网络分析最短路径分析源代码（VB6.0）精品资料.doc》由会员分享，可在线阅读，更多相关《ArcGIS网络分析最短路径分析源代码（VB6.0）精品资料.doc（30页珍藏版）》请在淘文阁 - 分享文档赚钱的网站上搜索。

1、ArcGIS网络分析最短路径分析源代码(VB6.0)12Copyright1995-2005ESRI34AllrightsreservedunderthecopyrightlawsoftheUnitedStates.56Youmayfreelyredistributeandusethissamplecode,withorwithoutmodification.78Disclaimer:THESAMPLECODEISPROVIDEDASISANDANYEXPRESSORIMPLIED9WARRANTIES,INCLUDINGTHEIMPLIEDWARRANTIESOFMERCHANTABILIT

2、YANDFITNESS10FORAPARTICULARPURPOSEAREDISCLAIMED.INNOEVENTSHALLESRIOR11CONTRIBUTORSBELIABLEFORANYDIRECT,INDIRECT,INCIDENTAL,SPECIAL,EXEMPLARY,12ORCONSEQUENTIALDAMAGES(INCLUDING,BUTNOTLIMITEDTO,PROCUREMENTOF13SUBSTITUTEGOODSORSERVICES;LOSSOFUSE,DATA,ORPROFITS;ORBUSINESS14INTERRUPTION)SUSTAINEDBYYOUORA

3、THIRDPARTY,HOWEVERCAUSEDANDONANY15THEORYOFLIABILITY,WHETHERINCONTRACT,STRICTLIABILITY,ORTORTARISINGINANY16WAYOUTOFTHEUSEOFTHISSAMPLECODE,EVENIFADVISEDOFTHEPOSSIBILITYOF17SUCHDAMAGE.1819Foradditionalinformationcontact:EnvironmentalSystemsResearchInstitute,Inc.2021Attn:ContractsDept.2223380NewYorkStre

4、et2425Redlands,California,U.S.A.923732627Email:contracts2829OptionExplicit3031vbversionofthePathFinderobject3233本地变量34Privatem_ipGeometricNetworkAsesriGeoDatabase.IGeometricNetwork35Privatem_ipMapAsesriCarto.IMap36Privatem_ipPointsAsesriGeometry.IPointCollection37Privatem_ipPointToEIDAsesriNetworkAn

5、alysis.IPointToEID38返回结果变量39Privatem_dblPathCostAsDouble40Privatem_ipEnumNetEID_JunctionsAsesriGeoDatabase.IEnumNetEID41Privatem_ipEnumNetEID_EdgesAsesriGeoDatabase.IEnumNetEID42Privatem_ipPolylineAsesriGeometry.IPolyline434445OptionallysettheMap(e.g.thecurrentmapinArcMap),46otherwiseadefaultmapwill

6、bemade(forIPointToEID).4748PublicPropertySetMap(MapAsesriCarto.IMap)49Setm_ipMap=Map50EndProperty5152PublicPropertyGetMap()AsesriCarto.IMap53SetMap=m_ipMap54EndProperty5556EitherOpenAccessNetworkorOpenFeatureDatasetNetwork57needstobecalled.5859PublicSubOpenAccessNetwork(AccessFileNameAsString,Featur

7、eDatasetNameAsString)6061DimipWorkspaceFactoryAsesriGeoDatabase.IWorkspaceFactory62DimipWorkspaceAsesriGeoDatabase.IWorkspace63DimipFeatureWorkspaceAsesriGeoDatabase.IFeatureWorkspace64DimipFeatureDatasetAsesriGeoDatabase.IFeatureDataset6566AfterthisSubexits,wellhaveanINetworkinterface67andanIMapint

8、erfaceinitializedforthenetworkwellbeusing.6869closedownthelastoneifopened70CloseWorkspace7172openthemdb73SetipWorkspaceFactory=NewesriDataSourcesGDB.AccessWorkspaceFactory74SetipWorkspace=ipWorkspaceFactory.OpenFromFile(AccessFileName,0)7576gettheFeatureWorkspace77SetipFeatureWorkspace=ipWorkspace78

9、79opentheFeatureDataset80SetipFeatureDataset=ipFeatureWorkspace.OpenFeatureDataset(FeatureDatasetName)8182initializeNetworkandMap(m_ipNetwork,m_ipMap)83IfNotInitializeNetworkAndMap(ipFeatureDataset)ThenErr.Raise0,OpenAccessNetwork,ErrorinitializingNetworkandMap8485EndSub8687PublicSubOpenFeatureDatas

10、etNetwork(FeatureDatasetAsesriGeoDatabase.IFeatureDataset)88closedownthelastoneifopened89CloseWorkspace9091weassumethatthecallerhaspassedavalidFeatureDataset9293initializeNetworkandMap(m_ipNetwork,m_ipMap)94IfNotInitializeNetworkAndMap(FeatureDataset)ThenErr.Raise0,OpenFeatureDatasetNetwork,Errorini

11、tializingNetworkandMap9596EndSub9798Thecollectionofpointstotravelthroughmustbeset.99100PublicPropertySetStopPoints(PointsAsesriGeometry.IPointCollection)101Setm_ipPoints=Points102EndProperty103104PublicPropertyGetStopPoints()AsesriGeometry.IPointCollection105SetStopPoints=m_ipPoints106EndProperty107

12、108Calculatethepath109110PublicSubSolvePath(WeightNameAsString)111112DimipNetworkAsesriGeoDatabase.INetwork113DimipTraceFlowSolverAsesriNetworkAnalysis.ITraceFlowSolver114DimipNetSolverAsesriNetworkAnalysis.INetSolver115DimipNetFlagAsesriNetworkAnalysis.INetFlag116DimipaNetFlag()AsesriNetworkAnalysi

13、s.IEdgeFlag117DimipEdgePointAsesriGeometry.IPoint118DimipNetElementsAsesriGeoDatabase.INetElements119DimintEdgeUserClassIDAsLong120DimintEdgeUserIDAsLong121DimintEdgeUserSubIDAsLong122DimintEdgeIDAsLong123DimipFoundEdgePointAsesriGeometry.IPoint124DimdblEdgePercentAsDouble125DimipNetWeightAsesriGeoD

14、atabase.INetWeight126DimipNetSolverWeightsAsesriNetworkAnalysis.INetSolverWeights127DimipNetSchemaAsesriGeoDatabase.INetSchema128DimintCountAsLong129DimiAsLong130DimvaRes()AsVariant131132makesureweareready133Debug.AssertNotm_ipPointsIsNothing134Debug.AssertNotm_ipGeometricNetworkIsNothing135136insta

15、ntiateatraceflowsolver137SetipTraceFlowSolver=NewesriNetworkAnalysis.TraceFlowSolver138139gettheINetSolverinterface140SetipNetSolver=ipTraceFlowSolver141142setthesourcenetworktosolveon143SetipNetwork=m_ipGeometricNetwork.Network144SetipNetSolver.SourceNetwork=ipNetwork145146makeedgeflagsfromthepoint

16、s147148theINetElementsinterfaceisneededtogetUserID,UserClassID,149andUserSubIDfromanelementid150SetipNetElements=ipNetwork151152getthecount153intCount=m_ipPoints.PointCount154Debug.AssertintCount1155156dimensionourIEdgeFlagarray157ReDimipaNetFlag(intCount)158159Fori=0TointCount-1160makeanewEdgeFlag1

17、61SetipNetFlag=NewesriNetworkAnalysis.EdgeFlag162SetipEdgePoint=m_ipPoints.Point(i)163lookuptheEIDforthecurrentpoint(thiswillpopulateintEdgeIDanddblEdgePercent)164m_ipPointToEID.GetNearestEdgeipEdgePoint,intEdgeID,ipFoundEdgePoint,dblEdgePercent165Debug.AssertintEdgeID0elsePoint(eid)notfound166ipNet

18、Elements.QueryIDsintEdgeID,esriETEdge,intEdgeUserClassID,intEdgeUserID,intEdgeUserSubID167Debug.Assert(intEdgeUserClassID0)And(intEdgeUserID0)elsePointnotfound168ipNetFlag.UserClassID=intEdgeUserClassID169ipNetFlag.UserID=intEdgeUserID170ipNetFlag.UserSubID=intEdgeUserSubID171SetipaNetFlag(i)=ipNetF

19、lag172Next173174addtheseedgeflags175ipTraceFlowSolver.PutEdgeOriginsintCount,ipaNetFlag(0)176177settheweight(costfield)tosolveon178179gettheINetSchemainterface180SetipNetSchema=ipNetwork181SetipNetWeight=ipNetSchema.WeightByName(WeightName)182Debug.AssertNotipNetWeightIsNothing183184settheweight(use

20、thesameforbothdirections)185SetipNetSolverWeights=ipTraceFlowSolver186SetipNetSolverWeights.FromToEdgeWeight=ipNetWeight187SetipNetSolverWeights.ToFromEdgeWeight=ipNetWeight188189initializearrayforresultstonumberofsegmentsinresult190ReDimvaRes(intCount-1)191192solveit193ipTraceFlowSolver.FindPathesr

21、iFMConnected,esriSPObjFnMinSum,m_ipEnumNetEID_Junctions,m_ipEnumNetEID_Edges,intCount-1,vaRes(0)194195computetotalcost196m_dblPathCost=0197Fori=LBound(vaRes)ToUBound(vaRes)198m_dblPathCost=m_dblPathCost+vaRes(i)199Next200201clearthelastpolylineresult202Setm_ipPolyline=Nothing203204EndSub205206Proper

22、tytogetthecost207208PublicPropertyGetPathCost()AsDouble209PathCost=m_dblPathCost210EndProperty211212Propertytogettheshape213214PublicPropertyGetPathPolyLine()AsesriGeometry.IPolyline215216DimipEIDHelperAsesriNetworkAnalysis.IEIDHelper217DimcountAsLong,iAsLong218DimipEIDInfoAsesriNetworkAnalysis.IEID

23、Info219DimipEnumEIDInfoAsesriNetworkAnalysis.IEnumEIDInfo220DimipGeometryAsesriGeometry.IGeometry221DimipNewGeometryCollAsesriGeometry.IGeometryCollection222DimipSpatialReferenceAsesriGeometry.ISpatialReference223224ifthelineisalreadycomputedsincethelastpath,justreturnit225IfNotm_ipPolylineIsNothing

24、Then226SetPathPolyLine=m_ipPolyline227ExitProperty228EndIf229230Setm_ipPolyline=NewesriGeometry.Polyline231SetipNewGeometryColl=m_ipPolyline232233apathshouldbesolvedfirst234Debug.AssertNotm_ipEnumNetEID_EdgesIsNothing235236makeanEIDHelperobjecttotranslateedgestogeometricfeatures237SetipEIDHelper=New

25、esriNetworkAnalysis.EIDHelper238SetipEIDHelper.GeometricNetwork=m_ipGeometricNetwork239SetipSpatialReference=m_ipMap.SpatialReference240SetipEIDHelper.OutputSpatialReference=ipSpatialReference241ipEIDHelper.ReturnGeometries=True242243getthedetailsusingtheIEIDHelperclasses244SetipEnumEIDInfo=ipEIDHel

26、per.CreateEnumEIDInfo(m_ipEnumNetEID_Edges)245count=ipEnumEIDInfo.count246247settheiteratortobeginning248ipEnumEIDInfo.Reset249250Fori=1Tocount251252getthenextEIDandacopyofitsgeometry(itmakesaClone)253SetipEIDInfo=ipEnumEIDInfo.Next254SetipGeometry=ipEIDInfo.Geometry255256ipNewGeometryColl.AddGeomet

27、ryCollectionipGeometry257258NextEID259260returnthemergedgeometryasaPolyline261SetPathPolyLine=m_ipPolyline262263EndProperty264265Private266267PrivateSubCloseWorkspace()268makesureweletgoofeverythingandstartwithnewresults269Setm_ipGeometricNetwork=Nothing270Setm_ipPoints=Nothing271Setm_ipPointToEID=N

28、othing272Setm_ipEnumNetEID_Junctions=Nothing273Setm_ipEnumNetEID_Edges=Nothing274Setm_ipPolyline=Nothing275EndSub276277PrivateFunctionInitializeNetworkAndMap(FeatureDatasetAsesriGeoDatabase.IFeatureDataset)AsBoolean278279DimipNetworkCollectionAsesriGeoDatabase.INetworkCollection280DimipNetworkAsesri

29、GeoDatabase.INetwork281DimcountAsLong,iAsLong282DimipFeatureClassContainerAsesriGeoDatabase.IFeatureClassContainer283DimipFeatureClassAsesriGeoDatabase.IFeatureClass284DimipGeoDatasetAsesriGeoDatabase.IGeoDataset285DimipLayerAsesriCarto.ILayer286DimipFeatureLayerAsesriCarto.IFeatureLayer287DimipEnve

30、lopeAsesriGeometry.IEnvelope,ipMaxEnvelopeAsesriGeometry.IEnvelope288DimdblSearchTolAsDouble289DimdblWidthAsDouble,dblHeightAsDouble290291OnErrorGoToTrouble292293getthenetworks294SetipNetworkCollection=FeatureDataset295296eventhoughaFeatureDatasetcanhavemanynetworks,welljust297assumethefirstone(othe

31、rwiseyouwouldpassthenetworknamein,etc.)298299getthecountofnetworks300count=ipNetworkCollection.GeometricNetworkCount301302Debug.Assertcount0thenException.Create(Nonetworksfound);303304getthefirstGeometricNewtork(0-based)305Setm_ipGeometricNetwork=ipNetworkCollection.GeometricNetwork(0)306307gettheNe

32、twork308SetipNetwork=m_ipGeometricNetwork.Network309310TheEIDHelperclassthatconvertspointstoEIDsneedsa311IMap,sowellneedonearoundwithallourlayersadded.312ThisPathfinderobjecthasanoptionalMappropertythanmaybeset313beforeopeningtheNetwork.314Ifm_ipMapIsNothingThen315Setm_ipMap=NewesriCarto.Map316317Ad

33、deachoftheFeatureClassesinthisGeometricNetworkasamapLayer318SetipFeatureClassContainer=m_ipGeometricNetwork319count=ipFeatureClassContainer.ClassCount320Debug.Assertcount0thenException.Create(No(network)featureclassesfound);321322Fori=0Tocount-1323getthefeatureclass324SetipFeatureClass=ipFeatureClas

34、sContainer.Class(i)325makealayer326SetipFeatureLayer=NewesriCarto.FeatureLayer327SetipFeatureLayer.FeatureClass=ipFeatureClass328addlayertothemap329m_ipMap.AddLayeripFeatureLayer330Next331EndIfweneededtomakeaMap332333334Calculatepointsnaptoleranceas1/100ofmapwidth.335count=m_ipMap.LayerCount336Setip

35、MaxEnvelope=NewesriGeometry.Envelope337Fori=0Tocount-1338SetipLayer=m_ipMap.Layer(i)339SetipFeatureLayer=ipLayer340getitsdimensions(forsettingsearchtolerance)341SetipGeoDataset=ipFeatureLayer342SetipEnvelope=ipGeoDataset.Extent343mergewithmaxdimensions344ipMaxEnvelope.UnionipEnvelope345Next346347fin

36、ally,wecansetuptheIPointToEID348Setm_ipPointToEID=NewesriNetworkAnalysis.PointToEID349Setm_ipPointToEID.SourceMap=m_ipMap350Setm_ipPointToEID.GeometricNetwork=m_ipGeometricNetwork351352setsnaptolerance353dblWidth=ipMaxEnvelope.Width354dblHeight=ipMaxEnvelope.Height355356IfdblWidthdblHeightThen357dbl

37、SearchTol=dblWidth/100#358Else359dblSearchTol=dblHeight/100#360EndIf361362m_ipPointToEID.SnapTolerance=dblSearchTol363364InitializeNetworkAndMap=Truegoodtogo365ExitFunction366367Trouble:368InitializeNetworkAndMap=Falsewehadanerror369EndFunction附录资料:不需要的可以自行删除bat文件的基本应用bat是dos下的批处理文件 .cmd是nt内核命令行环境的另

38、一种批处理文件 从更广义的角度来看，unix的shell脚本以及其它操作系统甚至应用程序中由外壳进行解释执行的文本，都具有与批处理文件十分相似的作用，而且同样是由专用解释器以行为单位解释执行，这种文本形式更通用的称谓是脚本语言。所以从某个程度分析，batch, unix shell, awk, basic, perl 等脚本语言都是一样的，只不过应用的范围和解释的平台各有不同而已。甚至有些应用程序仍然沿用批处理这一称呼，而其内容和扩展名与dos的批处理却又完全不同。 = 首先批处理文件是一个文本文件，这个文件的每一行都是一条DOS命令（大部分时候就好象我们在DOS提示符下执行的命令行一样），你

39、可以使用DOS下的Edit或者Windows的记事本(notepad)等任何文本文件编辑工具创建和修改批处理文件。 = 注 = 批处理文件中完全可以使用非dos命令，甚至可以使用不具有可执行特性的普通数据性文件，这缘于windows系统这个新型解释平台的涉入，使得批处理的应用越来越边缘化。所以我们讨论的批处理应该限定在dos环境或者命令行环境中，否则很多观念和设定都需要做比较大的变动。 = 其次，批处理文件是一种简单的程序，可以通过条件语句(if)和流程控制语句(goto)来控制命令运行的流程，在批处理中也可以使用循环语句(for)来循环执行一条命令。当然，批处理文件的编程能力与C语言等编程语

40、句比起来是十分有限的，也是十分不规范的。批处理的程序语句就是一条条的DOS命令（包括内部命令和外部命令），而批处理的能力主要取决于你所使用的命令。 = 注 = 批处理文件(batch file)也可以称之为批处理程序(batch program)，这一点与编译型语言有所不同，就c语言来说，扩展名为c或者cpp的文件可以称之为c语言文件或者c语言源代码，但只有编译连接后的exe文件才可以称之为c语言程序。因为批处理文件本身既具有文本的可读性，又具有程序的可执行性，这些称谓的界限是比较模糊的。 = 第三，每个编写好的批处理文件都相当于一个DOS的外部命令，你可以把它所在的目录放到你的DOS搜索路径(path)中来使得它可以在任意位置运行。一个良好的习惯是在硬盘上建立一个bat或者batch目录（例如C:BATCH），