#調(diào)用argis接口 根據(jù)點(diǎn)查找路徑,并將路徑保存到shp文件中
import arcpy
from arcpy import env

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

    #Set environment settings  配置argis的當(dāng)前基礎(chǔ)數(shù)據(jù)的數(shù)據(jù)文件地址
    env.workspace = "C:/Documents/ArcGIS/routeana0.gdb"
        env.overwriteOutput = True
    
    #Set local variables  配置好的網(wǎng)絡(luò)數(shù)據(jù)集(NetworkDataset)位置  即基礎(chǔ)數(shù)據(jù)
    inNetworkDataset = "train/train_ND"
    #輸出的經(jīng)過點(diǎn)的圖層名稱
    outNALayerName = "StationRoute"
    impedanceAttribute = "Length"
    inAddressLocator = "train_station_copy"
    #所有經(jīng)過點(diǎn)的圖層位置
    allinFeatures = "s2sdis"
    #輸入經(jīng)過點(diǎn)數(shù)據(jù)表位置
    inAddressTable = "G:/ArcGIS/python/StopAddresses.csv"
    #輸入經(jīng)過點(diǎn)數(shù)據(jù)表 字段
    inAddressFields = "Name Name VISIBLE NONE"
    #輸出經(jīng)過點(diǎn)圖層名稱
    outStops = "GeocodedStops"
    #輸出經(jīng)過點(diǎn)圖層保存地址
    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  根據(jù)名稱獲取對應(yīng)圖層
    stopsLayerName = subLayerNames["Stops"]
    routesLayerName = subLayerNames["Routes"]
    
    #需要分析多次 篩選多批次的經(jīng)過點(diǎn)  數(shù)據(jù)量大時(shí)可分批執(zhí)行
    for i in range(int(start),int(end)) :
        篩選當(dāng)前批次的經(jīng)過點(diǎn)到臨時(shí)圖層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.
        #網(wǎng)絡(luò)數(shù)據(jù)集中清除原來的點(diǎn)，并加入輸出經(jīng)過點(diǎn)圖層中的點(diǎn)
        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
        #執(zhí)行查找路徑方法  可能存在找不到路徑的情況 捕捉異常,不打斷循環(huán)
        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
            # 刪除臨時(shí)圖層，避免下次創(chuàng)建圖層時(shí)出錯(cuò)
            arcpy.Delete_management("station2station0")
            print str(e)
            print str(i)
            continue
    
        # 刪除臨時(shí)圖層，避免下次創(chuàng)建圖層時(shí)出錯(cuò)
        arcpy.Delete_management("station2station0")
        #Get the polygons sublayer from the service area layer
        #獲取結(jié)果中路徑的圖層
        routes_sublayer = arcpy.mapping.ListLayers(outNALayer,
                                            routesLayerName)[0]
        #將獲得的路徑匯總到時(shí)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))