Latest Contributions by KentaOkuyama

Hello, I really appreciate if anyone can help me out with this. I want to calculate network distance from each subject (n=8000) to a facility (n=290) . Each subject is supposed to go to one specific facility which is in their town.  Therefore, my approach was to 1. create point FC of subjects with the field "townname" which identifies where a subject lives. 2. create point FC of facilities with the field "townname" which identifies where a facility locates at.  3. execute MakreRouteAnalysisLayer by using field mapping However after many attempts, only around 100 lines are created as output.  I want to have 8000 lines as output which each line has distance value to the associated facility. If anyone could give any suggestions on this, I would really appreciate. My codes below:  env.overwriteOutput = True # calculate distance from subject to associated hospital inPath1 = r"C:\Users\kenta\OneDrive\GISData_Import\ADF2012_DISSOLVE.gdb" inPath2 = r"G:\ResearchProjects_data\gen_cohre\geocode.gdb" inPath3 = r"G:\ResearchProjects_data\gen_cohre\Walkabiliy_Indices.gdb" outPath = r"G:\ResearchProjects_data\gen_cohre\Walkability_final.gdb" network = os.path.join(inPath1, "ADF", "ADF_ND") stops_home = os.path.join(inPath3, "hospital") stops_work = os.path.join(inPath2, "subject") layer_name = "hosDistance" out_routes_featureclass = os.path.join(outPath,"hosRoutes") travel_mode = "Driving Time" #Create a new Route layer. result_object = arcpy.na.MakeRouteAnalysisLayer(network, layer_name, travel_mode) #Get the layer object from the result object. #The route layer can now be referenced using the layer object. layer_object = result_object.getOutput(0) #Get the names of all the sublayers within the route layer. sublayer_names = arcpy.na.GetNAClassNames(layer_object) #Stores the layer names that we will use later stops_layer_name = sublayer_names["Stops"] routes_layer_name = sublayer_names["Routes"] #Before loading the subject home and hospital locations as route stops, set #up field mapping. Map the "townname" field from the input data to #the RouteName property in the Stops sublayer, which ensures that each #unique townname will be placed in a separate route.  field_mappings = arcpy.na.NAClassFieldMappings(layer_object, stops_layer_name) field_mappings["RouteName"].mappedFieldName = "townname" #Add the home and hospital locations as Stops. The same field mapping #works for both input feature classes because they both have a field called #"townname" arcpy.na.AddLocations(layer_object, stops_layer_name, stops_home, field_mappings, "") arcpy.na.AddLocations(layer_object, stops_layer_name, stops_work, field_mappings, "", append="APPEND") #Solve the route layer. arcpy.na.Solve(layer_object) # Get the output Routes sublayer and save it to a feature class routes_sublayer = layer_object.listLayers(routes_layer_name)[0] arcpy.management.CopyFeatures(routes_sublayer, out_routes_featureclass) ... View more

I'm trying to plot GPS data in ArcGIS pro for the first time. The raw data is in text format and looks like below $GPRMC,151144.00,A,3335.69893,N,13012.84953,E,5.219,253.95,171218,,,A*68 2136198633 0000004023 $GPRMC,151152.00,A,3335.69887,N,13012.83373,E,2.715,255.40,171218,,,A*65 2136723401 0000004023 $GPGGA,151152.00,3335.69887,N,13012.83373,E,1,05,5.07,331.5,M,27.8,M,,*52 $GPRMC,151153.00,A,3335.70102,N,13012.83402,E,2.412,285.11,171218,,,A*64 2146959215 0000004023 . . . I found that GPX to Feature class was easy and fast, I first converted this text file to .gpx file using website service and succeeded in plotting the data in ArcGIS Pro. 1. GPX to Feature class (tool_gpx.PNG) 2. Successfully mapped (map.PNG) After that, I could convert the point data to line by Point to Lines tool. 3. Converted to point to line (line.PNG) It looks pretty good, but I wonder if there is a way to size or highlight parts of the line that had overlapped. In other words, the car tracked by GPS did not drive straight from the start to the end point, but it traveled back and force sometime.  I would really appreciate if anyone is familiar with GPS and could help answer my question.  ... View more

Hi,  I'm using arcpy.mp to export multiple layers at once, but having a trouble that legend is not updated. Each layer is successfully exported as separate JPEG file by the code below. However, the legend attached to each exported output is all the same (legend of the first layer).  I usually run my script as standalone from Jupyter notebook. When I run the code below in python window within ArcGIS pro, legend was successfully updated in each exported output. If anybody has any idea on why it cannot be done by standalone, I would really appreciate it. ArcGIS Pro 2.2 # This is the code I ran aprx = arcpy.mp.ArcGISProject(r"C:\Users\kenta\Documents\ArcGIS\Projects\CoHRE_Collaboration\CoHRE_Collaboration.aprx") maps = aprx.listMaps("Oki")[0] layers = maps.listLayers("mean*") layout = aprx.listLayouts("Oki_layout")[0] for layer in layers:    layer.visible=False for layer in layers:    layer.visible = True    layout.exportToJPEG(r"C:\Users\kenta\Dropbox\Collaboration\CollaborationOki"+ "\\" + layer.name + ".JPEG")    layer.visible = False ... View more