Latest Contributions by DavidPike

‎02-26-2020

Hi Blake thanks for your reply, unfortunately the dataset is at work. I'm trying to count the number of exactly overlapping polyline geometries (road segments that have been crossed multiple times in a network analyst route) for some reason I get the correct join count appended but with missing rows. Keep all is enabled and an outer join. I've got around this with a search cursor and geometry token but it seemed strange behaviour.

‎02-25-2020

I'm attempting to do a spatial join of polylines in desktop 10.4.1. The join is against itself to create a join count field for features within the same feature class. I.e. the target is itself. I just want each polygon segment to have a join count I dont really care if theres overlapping segments. The problem is that I seem to be losing records even though I have 'keep all' specified and this is also the default. Thanks for your time.

‎02-25-2020

Hi, You can assign a default value in the fields tab of the feature class, right-click properties on the FC in catalogues or catalog window. If you need this done in python can you screenshot your inputs to the assign default to field method?

‎02-25-2020

The output is the position of the raster with the lowest number for that cell. Ie. Its saying raster 3 has the lowest, faster 2 has the lowest etc. On a cell by cell basis. Dependent on the order you add them to the tool.

‎02-24-2020

Hello, I'm attempting to set up a geoprocessing task to be set up on portal 10.7.1 and desktop 10.5.1. My desktop script tool takes user input by them filling out rows in a table as an input parameter. This is done using the recordset data type (much in the same way as a feature set allows interactive drawing) and a template gdb table located in the datastore. Regardless of any other functions of the tool, when published the input options for the geoprocessing widget are 'table in map' or 'json'. Table in map is useless as it just takes the data from a table in the layers, with no interactive editing allowed, and the json is just a string which my users would have no idea what to do with. Thanks for your time

‎02-24-2020

Create a taster for A and a taster for B. Ensure the same cell size and snap raster are specified in the environments. Turn NoData values to 0 in both tasters (reclassify). Add together both rasters in raster calculator. Reclassify the output raster values. E.g. if A =1, B = 2 there will be a raster cell with a value of 3 at the junction. Turn this value back to 1.

‎02-23-2020

Can you not use the 'group' option and do it for every 'name' or 'objectid' etc. ?

‎02-23-2020

Ah ok I'll throw in a try except block. I guess its erorring on those?

‎02-23-2020

OK, but you did mention by retail type also. do you mean you want a density raster of all of the locations, then separate density rasters for each of the retail types? e.g. product 1 - density of all retail prod 2 - density of supermarkets prod 3 - density of pet stores … ? also you want this as a point density and not a kernel density?

‎02-23-2020

I'm still unsure. Please elaborate

‎02-23-2020

You have a double quote at the end of the selection on the first one and a = f" for the second one

‎02-23-2020

yes, thought the list() was retuning objects.. I've edited the code in my original reply which should work.

‎02-22-2020

import os import arcpy #enter your fgb input and output paths here file_gdb = r'your gdb to loop thru' output_gdb = r'output gdb' arcpy.env.workspace = file_gdb datasets = arcpy.ListDatasets("", "Feature") fd_path_list = [ ] fd_name_dict = { } for feature_dataset in datasets: fd_path = os.path.join(file_gdb, feature_dataset) fd_path_list.append(fd_path) fd_name_dict[fd_path] = str(feature_dataset) print(fd_name_dict) for fd in fd_path_list: fc_path_list_WD = [ ] fc_path_list_WED = [ ] arcpy.env.workspace = fd featureclasses = arcpy.ListFeatureClasses() for fc in featureclasses: if "_WD" in fc: fc_path = os.path.join(fd, fc) fc_path_list_WD.append(fc_path) if "_WED" in fc: fc_path = os.path.join(fd, fc) fc_path_list_WED.append(fc_path) wd_path = os.path.join(output_gdb, fd_name_dict[fd]) print(fd_name_dict[fd]) wd_path = wd_path + "_WD" print(wd_path) if len(fc_path_list_WD) != 0: try: arcpy.Merge_management(fc_path_list_WD, wd_path) except: continue wed_path = os.path.join(output_gdb, fd_name_dict[fd]) wed_path = wed_path + "_WED" print(wed_path) if len(fc_path_list_WED) != 0: try: arcpy.Merge_management(fc_path_list_WED, wed_path) except: continue del fc_path_list_WD del fc_path_list_WED print("done") ‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍

‎02-22-2020

Yes but I need to know what the target shapefiles look like

‎02-22-2020

The only thing is the spatial join. What are the shapefiles and how are they structured? I'll need to define this in the code

Online Status	Offline
Date Last Visited	‎01-14-2026 12:10 PM

My Ideas

Latest Contributions by DavidPike

Re: Spatial Join

Spatial Join

Re: Python help for default value in attribute field

Re: Polyline to Raster - any way to set value preference

Record Set interactive table input

Re: Polyline to Raster - any way to set value preference

Re: How do I use the "Summarize Nearby" tool to assess multiple locations at once?

Re: Loop through feature datasets in a gdb

Re: Points density

Re: Points density

Re: Model tool to Script not running

Re: Loop through feature datasets in a gdb

Re: Loop through feature datasets in a gdb

Re: Loop through feature datasets in a gdb

Re: Loop through feature datasets in a gdb

Re: Define Projection crack-down

Re: Generate polygons that represent historic air ...

Re: Enterprise geodatabase connection slow - machi...

Re: Discrepancy between measure tool and expected ...

Re: Why does the geographic transformation involvi...

GIS Jobs