POST
|
Thanks, Rachel. I had been running this in the foreground, but I do believe I have the 64bit background processing installed. But it ran without a hitch once I installed the updated version of PuLP. Thanks for your help!
... View more
11-10-2017
04:23 PM
|
0
|
0
|
1479
|
POST
|
Deelesh Mandloi, jsandhu-esristaff I have used that linked resource (Allocation Tools for ArcGIS) previously (with ArcMap 10.1). I'm trying to use it again(with 10.5), but keep getting an error: [Errno 13] Permission denied: '7476-pulp.mps'. The error is thrown at the point where the LP model is solved (line 479, prob.solve()). There was a similar issue posted in the comments on the page you link a year ago with no response. Do you have any suggestions on how to address this issue? Thanks.
... View more
11-07-2017
11:02 AM
|
0
|
15
|
1479
|
POST
|
Thanks, Mark. I've already tried both of these. I set both the processing extent and raster cell size to match the input raster in the Fill tool environment, and was already working locally.
... View more
11-05-2013
05:03 PM
|
0
|
0
|
8186
|
POST
|
I'm trying to fill a DEM. I keep getting multiple Error 010429 messages. I can't find great documentation on this problem. I've tried running the tool on multiple input DEMs, saved in various locations, with outputs written to multiple formats. The messages associated with the error varies, but includes things like "Error in GRID ID: CalcPackedBytes: invalid depth: -1," and "Error in GRID IO: S_VXfer: Read error". It is often accompanied by Error 010244: Could not set the analysis window. I can't find any documentation on these. Has anyone else encountered these problems, and found a solution? Thanks,
... View more
11-04-2013
09:08 AM
|
0
|
11
|
19738
|
POST
|
In case anyone is interested, I figured this out. You can access the sublayers in a Network Analyst layer in a python script launched from command line using the GetNAClassNames function. I also got the multiprocessing thing running successfully, but it keeps bailing out part way through for no apparent reason that I can figure out. It sort of looks like one of the jobs is failing to add a field to the output table, and this ends the whole thing. Frustrating.
... View more
07-13-2013
07:50 PM
|
0
|
0
|
176
|
POST
|
It sounds like you are saying you have two input tables, each with x/y coordinates in each row, and you are wanting to see if the coordinates for rows in one table match the coordinates for rows in the other table. Is that correct? If so, I don't know of an existing function to do it. But you could likely proceed using the following steps: 1. Join the two tables 2. Create a Boolean variable to be used to indicate if the coordinates match 3. Use a search cursor to run through the rows, and compare the x's and y's for each row. Set the Boolean variable to false if the don't match. Alternatively, if you don't care about slight precisional differences, you could calculate the Euclidean distance between the coordinates, and set the Boolean to false only if it exceeds some threshold. 4. Report out the value of that variable at the end. If the value of the variable is false at the end, you know that the coordinates don't match at some point.
... View more
07-13-2013
08:25 AM
|
0
|
0
|
349
|
POST
|
Hello, all. I've got a python code that I'm trying to get ready for use with multiprocessing (using StacyR's helpful code here). I've got a list of roads in an area threatened by wildfire, and I'm trying to see which are the most critical to protect by seeing how many residents will loose access to critical facilities if a given road is set as non-traversable. So, I'm running through a list of the pertinent roads (around 20,000, which is why I'm trying to get multiprocessing to work), and in each run, I set the selected road as non-traversable, solve the closest facility problem, and generate a summary statistics table to determine the number of pop who no longer have access. The code below is the function that is called for each road in the list of 20,000. It solves the analysis, does some calculations, generates a summary statistics table, and writes it out to a workspace. It all works perfect if I run it as an in-process script tool in ArcMap, but when I try to run it through the command line, or not in process (out process?) through ArcMap, it fails when I try to use MakeTableView (in the 4th line) to get the information from the closest facilities layer routes. Does anyone have any idea what could be causing this? EDIT I forgot to note: when it hits that line, it says path\analysis_layer.lyr\Routes does not exist (where path\analysis_layer.lyr is the path and name of the closest facility layer saved to disk), even though it just used that layer file to solve the closest facility problem in the previous line. END EDIT
def findRoadImportance(current_road,analysis_layer,table_view,OWS):
arcpy.UpdateAnalysisLayerAttributeParameter_na(analysis_layer, "Drop_list", "Dropped Road", current_road)
arcpy.Solve_na(analysis_layer,"SKIP", "TERMINATE", "")
arcpy.MakeTableView_management(analysis_layer + "\\Routes", table_view, "", "", "Name Name VISIBLE NONE;Total_Length Total_Length VISIBLE NONE")
arcpy.DeleteRows_management(analysis_layer + "\\Routes")
arcpy.AddField_management(table_view, "Pop", "LONG")
arcpy.AddField_management(table_view, "PW","DOUBLE")
arcpy.CalculateField_management(table_view, "Pop", "getPop( !Name!)", "PYTHON_9.3", "def getPop(field):\\n return field[:field.find(\" \")]")
arcpy.CalculateField_management(table_view,"PW", "calcMean( !Pop!, !Total_Length!)", "PYTHON_9.3",
"def calcMean(field1,field2):\\n return field1 * field2")
out_table = OWS + "\\" + str(int(current_road)) + "_result.dbf"
arcpy.Statistics_analysis(table_view, out_table, [["Pop","SUM"],["PW","SUM"]])
arcpy.AddField_management(out_table, "Road_ID","LONG")
arcpy.CalculateField_management(out_table, "Road_ID", current_road)
arcpy.AddMessage(" Output table written with drop road ID %s..."%(current_road))
return
Also, I'm wondering if multiprocessing will actually work with this approach. Since I'm using UpdateAnalysisLayerAttributeParameter to iteratively set roads as non-traversable in the same underlying Closest Facilities layer, can this truly be split apart between cores? Or will the results be messed up because changes made in the underlying layer by one job will influence the results in other jobs? Thanks for any help you all can give, matt
... View more
07-12-2013
02:25 PM
|
0
|
2
|
2133
|
POST
|
Awesome, thanks, Melinda. That is just what I was looking for. MB is getting really good in the last release or two. matt
... View more
07-03-2013
08:30 AM
|
0
|
0
|
227
|
POST
|
Hey, all, I'm looking to rank how critical a given subset of road segments are within a transportation network by comparing how many people would loose access to a set of facilities if one road were closed versus another. I'm trying to do this using a closest facility analysis with a road network in which I've created a parameterized attribute that allows me (in the ArcMap interface) to supply the id for a road segment that is then dropped from the analysis. Once I got this working, I wanted to apply the concept in Model Builder, and iterate through a list of roads to remove, and then re-solve, the closest facility problem for each iteration. The problem is I can't see how to specify the value to use for my parameterized attribute through the Model Builder interface. Does anyone know if this is possible, and if it is, how to make it work? Thanks for the help, matt
... View more
07-02-2013
04:22 PM
|
0
|
2
|
856
|
POST
|
Thanks, Curtis. That is what I have done. I'm just leary of outputs I can't explain. I always feel like if I can't explain what happens in one place, I can't be sure I know what is going on elsewhere.
... View more
08-22-2012
12:29 PM
|
0
|
0
|
204
|
POST
|
I calculated the local standard deviation of 1,000 input rasters using the cell statistics tool, and many of the cells returned values of -1.#IND. I haven't been able to find anywhere what it means, or how to address it. When I export it as an ASCII grid to use in another program (R), this value imports as a -1, which is problematic. It looks like these values may be returned for cells where the input values are identical across the local cell values. (I haven't checked everyone of them, though). Does anyone know what this is? I've tried searching the ESRI support, as well as Python, and haven't been able to find anything yet.
... View more
08-21-2012
05:21 PM
|
0
|
3
|
464
|
POST
|
In case this is helpful to anyone: I have been using Path Distance to model pedestrian evacuations, using reclassified landcover as a cost surface, and elevation as the surface layer and vertical factor, all with 1 meter resolution data. When I was running it with ArcGIS Desktop 10, it would run for about 24 hours, and result in dramatic underestimates of travel times. I still had my copies of 9.3.1, so I removed 10 reinstalled that, and the whole process took about an hour, and is returning more reasonable results. So I am suspecting the issue here is something to do with changes made at the 10 release. I like 10 overall, but Path Distance doesn't seem to work well in it.
... View more
11-04-2011
09:20 AM
|
0
|
0
|
1440
|
POST
|
Any news on this? I'm trying to figure out some issues with an analysis I'm running for a moderately sized study area, but with 1 meter pixels. It is an area of pretty low relief, but we are getting what seem to be dramatically low estimates for pedestrian travel times. It worked, however, before I upgraded to 10. Were there changes in the Path Distance tool at the 10 release that are causing this issue?
... View more
11-02-2011
06:48 AM
|
0
|
0
|
1440
|
POST
|
Sander, If you are doing cost distance analysis (rather than path distance), I would suggest this: input your end point in the model as a starting point. This will calculate the least cost distance from every point in your study area to the hill top. Then once you have your points around the hill defined, you can simply do a zonal function to determine the cost at each point. This will only work, however, if you are using isotropic approaches to model cost distances. This means that the direction of travel doesn't matter, and is implemented using the Cost Distance tool. You can still do something like this using Path Distance, but you have to do a bit more work.
... View more
11-02-2011
06:32 AM
|
0
|
0
|
256
|
POST
|
Jay, Thanks for the help. This was the first time I've ever used a geometric network. It worked out great. Did just what I needed. Thanks.
... View more
10-05-2011
06:49 AM
|
0
|
0
|
341
|
Online Status |
Offline
|
Date Last Visited |
11-11-2020
02:23 AM
|