使用Arcgis做路径规划方式(使用python脚本调用)

#调用argis接口 根据点查找路径,并将路径保存到shp文件中
import arcpy
from arcpy import env

try:
    #Check out the Network Analyst extension license
    arcpy.CheckOutExtension("Network")

    #Set environment settings  配置argis的当前基础数据的数据文件地址
    env.workspace = "C:/Documents/ArcGIS/routeana0.gdb"
        env.overwriteOutput = True
    
    #Set local variables  配置好的网络数据集(NetworkDataset)位置  即基础数据
    inNetworkDataset = "train/train_ND"
    #输出的经过点的图层名称
    outNALayerName = "StationRoute"
    impedanceAttribute = "Length"
    inAddressLocator = "train_station_copy"
    #所有经过点的图层位置
    allinFeatures = "s2sdis"
    #输入经过点数据表位置
    inAddressTable = "G:/ArcGIS/python/StopAddresses.csv"
    #输入经过点数据表 字段
    inAddressFields = "Name Name VISIBLE NONE"
    #输出经过点图层名称
    outStops = "GeocodedStops"
    #输出经过点图层保存地址
    outLayerFile = "G:/ArcGIS/python" + "/" + outNALayerName + ".lyr"
    searchTolerance = "3000 Meters"
    
    #Create a new Route layer. For this scenario, the default value for all the
    #remaining parameters statisfies the analysis requirements
    outNALayer = arcpy.na.MakeRouteLayer(inNetworkDataset, outNALayerName,
                                         impedanceAttribute)
    
    #Get the layer object from the result object. The route layer can now be
    #referenced using the layer object.  得到输出图层集
    outNALayer = outNALayer.getOutput(0)
    
    #Get the names of all the sublayers within the route layer.
    #得到输出图层名称
    subLayerNames = arcpy.na.GetNAClassNames(outNALayer)
    #Stores the layer names that we will use later  根据名称获取对应图层
    stopsLayerName = subLayerNames["Stops"]
    routesLayerName = subLayerNames["Routes"]
    
    #需要分析多次 筛选多批次的经过点  数据量大时可分批执行
    for i in range(int(start),int(end)) :
        筛选当前批次的经过点到临时图层station2station0
        inFeatures = arcpy.Select_analysis(allinFeatures, "station2station0", '"sort" = '+str(i))
        
        #Load the geocoded address locations as stops mapping the address field from
        #geocoded stop features as Name property using field mappings.
        #网络数据集中清除原来的点，并加入输出经过点图层中的点
        arcpy.na.AddLocations(outNALayer, stopsLayerName, inFeatures, "Name Name #",
                                  searchTolerance, append = "CLEAR",exclude_restricted_elements = "EXCLUDE")
    
        #Solve the route layer, ignore any invalid locations such as those that
        #can not be geocoded
        #执行查找路径方法  可能存在找不到路径的情况 捕捉异常,不打断循环
        try:
            result = arcpy.na.Solve(outNALayer,"SKIP")
        except Exception as e:
            # If an error occurred, print line number and error message
            import traceback, sys
            tb = sys.exc_info()[2]
            print "An error occured on line %i" % tb.tb_lineno
            # 删除临时图层，避免下次创建图层时出错
            arcpy.Delete_management("station2station0")
            print str(e)
            print str(i)
            continue
    
        # 删除临时图层，避免下次创建图层时出错
        arcpy.Delete_management("station2station0")
        #Get the polygons sublayer from the service area layer
        #获取结果中路径的图层
        routes_sublayer = arcpy.mapping.ListLayers(outNALayer,
                                            routesLayerName)[0]
        #将获得的路径汇总到时s2s的shp中
        arcpy.Append_management([routes_sublayer], "s2s", "TEST")

    print("Script completed successfully")

except Exception as e:
    # If an error occurred, print line number and error message
    import traceback, sys
    tb = sys.exc_info()[2]
    print ("An error occured on line %i" % tb.tb_lineno)
    print (str(e))