Latest Contributions by sotokan80_

I am using arcpy and ArcGIS 10.4.1. Is there a way to set curb approach programmatically based on the road's hierarchy. i.e in residential roads curb approach will be set to either side and in primary roads curb approach will be set to right side. ... View more

I have a network dataset and i want to find out a way to show oneways direction (an arrow of something similar) without having to use network identify tool. What i want to achieve is to share a shapefile with certain people in order to get info back concerning how accurate is current oneways information. There is a field named oneway with values T or F of B but i have to know for every segment the start and end point. Any ideas? ... View more

I am using ArcGis 10.4.1 and network analyst extension. I am facing an issue when loading barrier lines using arcpy. It appears that the network elements that have restrictions in network dataset attributes have no network location information. If i load the line barriers using ArcMap and a New Route layer everything works ok, if i click "Use Geometry" in Location Position frame, whilst if i use the arcpy created layer if i click  "Use Geometry" in Location Position frame and if i add a new line barrier that is crossing an network element that has restriction i get error: Barrier "whatever" in "Line Barriers" has no associated network location information. If the line barrier doesnt cross elements with restrictions i get no error. Obviously the arcpy created layer has a miconfiguration. Part of the code that i use: ## Create a new Route layer outRouteResultObject = arcpy.na.MakeRouteLayer(inNetworkDataset, outNALayerName, impedanceAttribute, "FIND_BEST_ORDER", ""+EndPoint+"", "", ["Minutes", "Length"], "ALLOW_DEAD_ENDS_AND_INTERSECTIONS_ONLY", RoadRules, "USE_HIERARCHY", "", "TRUE_LINES_WITH_MEASURES","") ##Get the layer object from the result object. The route layer can now be referenced using the layer object. outNALayer = outRouteResultObject.getOutput(0) ##Get the names of all the sublayers within the route layer. subLayerNames = arcpy.na.GetNAClassNames(outNALayer) ##Store the layer names that we will use later stopsLayerName = subLayerNames["Stops"] routesLayerName = subLayerNames["Routes"] ##set up field mapping. Map the "Group_ID" field from the input data to the RouteName property in the Stops sublayer fieldMappingsStore = arcpy.na.NAClassFieldMappings(outNALayer, stopsLayerName) fieldMappingsStore["RouteName"].mappedFieldName = "Group_ID" fieldMappingsStops = arcpy.na.NAClassFieldMappings(outNALayer, stopsLayerName) fieldMappingsStops["RouteName"].mappedFieldName = "Group_ID" fieldMappingsStops["Attr_Minutes"].defaultValue = 2 ##Add locations as Stops. arcpy.na.AddLocations(outNALayer, stopsLayerName, inDepots, fieldMappingsStore, "300 meters", "OBJECTID", [["networkDataset_ND", "MIDDLE"], ["networkDataset_ND_Junctions", "NONE"]], "MATCH_TO_CLOSEST", "CLEAR", "NO_SNAP", "", "", "") arcpy.na.AddLocations(outNALayer, stopsLayerName, inOrders, fieldMappingsStops, "300 meters", "OBJECTID", [["networkDataset_ND", "SHAPE"], ["networkDataset_ND_Junctions", "NONE"]], "MATCH_TO_CLOSEST", "APPEND", "SNAP", "1", "EXCLUDE", searchQuery) #Figure out the layer's impedance attribute solver_props = arcpy.na.GetSolverProperties(outNALayer) impedance = solver_props.impedance pointBarriersLayerName = subLayerNames["PolylineBarriers"]# {u'Barriers': u'Point Barriers', u'Routes': u'Routes', u'PolylineBarriers': u'Line Barriers', u'Stops': u'Stops', u'PolygonBarriers': u'Polygon Barriers'} etc fieldMappingsBR = arcpy.na.NAClassFieldMappings(outNALayer, pointBarriersLayerName) fieldMappingsBR["BarrierType"].defaultValue = 0 #if it is restriction otherwise 1 fieldMappingsBR["Attr_" + impedance].defaultValue = 1.4 #Load restrictions features arcpy.na.AddLocations(outNALayer, pointBarriersLayerName, Barriers_shp, fieldMappingsBR) Can you spot the mistake ... View more

I am facing a frustrating issue, when using oneways, (ArcGIS 10.4.1), as restriction the solved VRP  returns more than one to the same road. An example is shown below. As you can see the driver is forced to return to the same location after 30 minutes of driving. Is there any way to fix this issue? ... View more

I am using ArcGIS 10.4.1 and arcpy. I am facing an issue when creating a table in a gdb using GenerateNearTable_analysis tool. The output table does not show its contents. This only occurs if the workspace is a geodatabase. If i restart arcmap it shows its contents but thats not the point since my code crashes at that point. Any ideas? ... View more

I am using Arcpy and ArcGIS 10.4.1. Is there any way to parse a shapefile's Shape field and get the value in the following format; 0xE6100000010CF091A1F8EB67354050B8D91D27644440. Currently i am storing in native sql server using the following query; insert into [TABLE_NAME] (Shape) SELECT geometry::Point(21.405944385000055, 40.782443744000034 , 4326) i.e # Get data from shapefile with arcpy.da.SearchCursor(table, ['Shape']) as cursor:          for row in cursor:                   Shape = row[0]                   print Shape >>>>(21.405944385000055, 40.782443744000034) And i want it to be 0xE6100000010CF091A1F8EB67354050B8D91D27644440 ... View more

I have a polyline shapefile (1.500.000 road segments) and  a point shapefile (the road segment's name). What is the most efficient way to update the polyline shapefile using info from the point shapefile. Do i have to use a specific tool i.e SelectLayerByLocation_management and a cursor? Is there a more efficient way? ... View more

I am using ArcGIS 10.4.1 and Network Analyst extension. As an IDE i use Spyder and Python version 2.7. I want to have fields in my shapefile that will present coordinates in WGS 84 and in a national Grid format respectively. The initial information is in WGS 84 format. Is it possible to calculate geometry of the fields that will present the coordinates in the national Grid format? Currently i am achieving that by updating the shapefile with a cursor and running for every line the following script which is very time consuming. point_utm = arcpy.PointGeometry(arcpy.Point(float(row[12]),float(row[13])),arcpy.SpatialReference("WGS 1984"), trsf) pointA = point_utm.projectAs(out_coordinate_system, trsf) row[10] = pointA.centroid.X row[11] = pointA.centroid.Y I dont want to use   # Generate the extent coordinates using Add Geometry Properties tool#  arcpy.AddGeometryAttributes_management(in_features, properties, length_unit, area_unit, coordinate_system) because it updates POINT_X and POINT_Y fields thus losing WGS84 coordinates format. ... View more

I am using ArcGis 10.4.1 for desktop and network Analyst extension. Is it possible to use multiprocessing when performing vehicle problem routing analysis? Every day i run a script where the VRP solver finds the best routes for a fleet of 300 motorbikes where each one has to visit up to 400 stops. If it is not feasible to use multiprocessing will upgrading my cpu and ram size reduce the processing time?  Currently i m using a pc equipped  with i3 processor and 16gb ram. ... View more