topic Re: Calculating distance between many users of a handful of facilities, a small hangup in ArcGIS Network Analyst Questions

Calculating distance between many users of a handful of facilities, a small hangup

chriscoutts — Sat, 07 Feb 2015 18:23:16 GMT

I know the ROUTENAME issue has been addressed many times in these forums over the past couple years. I have scoured them and am hung up on one point.

In the shortest distance thread it is stated that "You might have to load multiple copies of the same park because multiple people use the same park." In my case, these are facilities.

In the Re: Finding Distances Between Assigned Locations thread it is stated that "...your 368 locations are repeating. You could do this with some relates...", but there is not direction on how to use relates to generate repeated "locations."

I need to figure out how to make hundreds of duplicate copies of the same facility. Once duplicated, I will then give each one of the hundreds of copies of the same facility a unique ID. Then I can match the facility to the thousands of users of it using ROUTENAME.

It is just duplicating the facility that is hanging me up. Please help. Thanks.

Re: Calculating distance between many users of a handful of facilities, a small hangup

DarrenWiens2 — Sat, 11 Dec 2021 19:00:25 GMT

If I'm reading this right, your main question is how to duplicate geometries. Here's how, using Python:

>>> duplicates = 100
>>> points = []
>>> with arcpy.da.SearchCursor("YOUR_POINT_LAYER_NAME",["SHAPE@"]) as cursor:
...    for row in cursor:
...        for i in range(duplicates):
...            points.append(row[0])
...           
>>> arcpy.CopyFeatures_management(points,'in_memory\out_points')
<Result 'in_memory\\out_points'>
>>> arcpy.SpatialJoin_analysis("out_points","YOUR_POINT_LAYER_NAME",r'in_memory\sp_join_output',"JOIN_ONE_TO_ONE","KEEP_ALL",'#',"ARE_IDENTICAL_TO")
<Result 'in_memory\\sp_join_output'>

Re: Calculating distance between many users of a handful of facilities, a small hangup

chriscoutts — Sat, 07 Feb 2015 19:45:09 GMT

I really appreciate the prompt response. Thanks very much.

In one feature class (facilities) I have 19 individual points. I am trying to make hundreds of copies of each point (facility). Will the script you provided tackle this? If so, I'm a novice GIS user, could you tell me where to plug in this script? Thanks again, Darren.

Re: Calculating distance between many users of a handful of facilities, a small hangup

DarrenWiens2 — Sat, 07 Feb 2015 19:56:33 GMT

This script will make 100 copies of each feature in the feature class "Your point layer name" (of course, change that to your later name). Change the number of duplicates if you want.

This script can go into the Python window, accessed under the geoprocessing menu in ArcMap.

Re: Calculating distance between many users of a handful of facilities, a small hangup

chriscoutts — Sat, 07 Feb 2015 20:09:09 GMT

Beautiful. I am going to give it a run, but it may be a couple days before I can get back to it. Thanks so much.

Re: Calculating distance between many users of a handful of facilities, a small hangup

chriscoutts — Sat, 07 Feb 2015 20:28:15 GMT

Runtime error Traceback (most recent call last): File "<string>", line 1, in <module> File "c:\program files\arcgis\desktop10.1\arcpy\arcpy\analysis.py", line 481, in SpatialJoin raise e ExecuteError: ERROR 000732: Target Features: Dataset out_points does not exist or is not supported ERROR 000732: Join Features: Dataset ex_centroid does not exist or is not supported

This is what is spit back at me. Running now, but will go back through the code to make sure I put in exactly what you provided.

Re: Calculating distance between many users of a handful of facilities, a small hangup

DarrenWiens2 — Sat, 07 Feb 2015 20:46:39 GMT

That was my fault - I ran my test on a feature class called "ex_centroid". I've updated the script to say "YOUR_POINT_LAYER_NAME" in the Spatial Join call (so, change that to your own).

Re: Calculating distance between many users of a handful of facilities, a small hangup

XanderBakker — Sun, 08 Feb 2015 20:13:36 GMT

You may be interested in this thread: Copy features a number of times based on a numeric field value

Re: Calculating distance between many users of a handful of facilities, a small hangup

chriscoutts — Tue, 10 Feb 2015 20:09:10 GMT

This is EXACTLY what I need to do. I have never used Python script before, but I've invested no less than 8 hours trying to use the "Copy features..." code you provided. Using Geoprocessing=>Python does not work. Lots of >>> and ... symbols that appear to be present in some codes (such as above), but not in yours.

I also created a .py file with the code, created a toolbox, and used the script wizard in an attempt to run it as well. No luck. Of course, I adapted your code to my own file location paths. That is where it goes wrong early on. A lot of red text results make the claim that the files do not exist even though I have confirmed many times that they do.

I would love to run this script if you think you can guide me through getting it running without burdening yourself too much.

Thank you.

Re: Calculating distance between many users of a handful of facilities, a small hangup

chriscoutts — Tue, 10 Feb 2015 20:16:29 GMT

I was really excited about getting this working, but I am failing right in line 1 somehow. Here's a screenshot:

I am really put some time into trying to understand where I was failing before bugging you, but I am stuck.

Re: Calculating distance between many users of a handful of facilities, a small hangup

DarrenWiens2 — Tue, 10 Feb 2015 20:48:47 GMT

You need matching quotes around "SHAPE@"

Re: Calculating distance between many users of a handful of facilities, a small hangup

chriscoutts — Wed, 25 Mar 2015 21:06:26 GMT

This is the code that eventually worked. The one in the Copy features a number of times based on a numeric field value thread did not

def main():

import arcpy

import os

fc_in = r"C:\Users\ccoutts\Desktop\Props\CRC_MultiGRant\Map1\Facilities1.shp" # this one exists

fld_count = "obs"

fc_out = r"C:\Users\ccoutts\Desktop\Props\CRC_MultiGRant\Map1\Facilities2.shp" # this one will be created

sr = arcpy.Describe(fc_in).spatialReference

# create the empty output featureclass

path, name = os.path.split(fc_out)

arcpy.CreateFeatureclass_management(path, name, "POINT", fc_in, "SAME_AS_TEMPLATE", "SAME_AS_TEMPLATE", sr)

# insert the features into the output fc

with arcpy.da.SearchCursor(fc_in, '*') as curs_in:

flds_in = curs_in.fields

idx_cnt = flds_in.index(fld_count)

with arcpy.da.InsertCursor(fc_out, '*') as curs_out:

for row in curs_in:

cnt = row[idx_cnt]

for i in range(0, cnt):

curs_out.insertRow(row)

main()