Latest Contributions by PaulHuffman

‎04-02-2015

I have that problem when my pour points are a little off the local minimums. But you said you used the Snap Pour Point Tool. In an edit session, try displaying your flow accumulation raster behind your pour points. Play with the symbology classification of your accumulation raster so you can see the stream lines distinct from the real low accumulation values. Move the pour points so they will sit right on top of a pixel in the accumulation with the highest value near where you thought the pour point should be. Sometimes, however, with real flat terrain, the flow modeling just doesn't work as expected. Especially with coarser DEM data like 30 m cells. You don't get more accurate data by resampling to 10m, you get each 30 m cell chopped into 9 cells with the same value (basically).

‎03-12-2015

I think this worked for me. I launched this script from IDLE, waited a couple minutes, the IDLE shell window showed a single set of >>> after the restart and wouldn't advance to a blank line, just beeped at me. CPU usage clicked along at 6 to 14%. Then I created the input raster by copying SimHAWS2000 to SimHAWS2000_1, and in a couple minutes both the output rasters came out. # Import arcpy module import arcpy import time from arcpy import env from arcpy.sa import * # Check out any necessary licenses arcpy.CheckOutExtension("spatial") arcpy.env.overwriteOutput = True # Local variables: SimHAWS2000_1 = "G:\\elevation\\lidar_2005\\ThreeWay\\ThreeWayTerrians.gdb\\SimHAWS2000_1" Input_true_raster_or_constant_value = "8" Input_true_raster_or_constant_value__2_ = "-8" Con_SimHAWS20001 = "G:\\elevation\\lidar_2005\\ThreeWay\\ThreeWayTerrians.gdb\\Con_SimHAWS20001" SimHAWS2kCon = "G:\\elevation\\lidar_2005\\ThreeWay\\ThreeWayTerrians.gdb\\SimHAWS2kCon" result = None while result is None: try: time.sleep(60) # Process: Con arcpy.gp.Con_sa(SimHAWS2000_1, Input_true_raster_or_constant_value, Con_SimHAWS20001, SimHAWS2000_1, "VALUE >= 8") # Process: Con (2) arcpy.gp.Con_sa(Con_SimHAWS20001, Input_true_raster_or_constant_value__2_, SimHAWS2kCon, Con_SimHAWS20001, "VALUE <= -8") result = "Both Cons executed" except: pass

‎03-12-2015

I hope to try this later today. I'm a bit concerned that a solution like this could gobble up CPU while it runs and fails. As one of the commenters on stackoverflow offers: "But any proposed answer should be safe, or at least note the pitfalls. This does not offer protection against 100% CPU consumption and endangers future readers. – Brad Koch May 30 '14 at 19:03" Putting a sleep inside the while or try might help reduce the CPU use, but I don't know for sure. I'm guessing sleep uses less than attempting an operation on a big raster and failing. Also needed will be env.overwriteOutput, since I have seen above that a tool can create an output file even though the input incomplete and the tool going to therefore fail.

‎03-11-2015

Thanks. I didn't expect to find this under ArcDesktop Help.

‎03-11-2015

Where can I find the matrix of databases compatible with ArcServer 10.3? Trying to decide if I should stay with the SQL server 2012 I have or get SQL server 2014. Looking at http://server.arcgis.com/en/server/latest/install/windows/system-requirements.htm#ESRI_SECTION1_FCC5755E241C42F5BA664230E496D0B7http:// and http://www.esri.com/~/media/Files/Pdfs/library/brochures/pdfs/arcgis-server-functionality-matrix.pdf but can't find it.

‎03-10-2015

I was able to test this, and it looks like Exists evaluates as true before the raster is completed. I ran a little script when the input raster, strm2M, was about 38% complete - import arcpy import time from arcpy import env from arcpy.sa import * #arcpy.CheckOutExtension("3D") arcpy.CheckOutExtension("Spatial") env.workspace = r"G:\elevation\lidar_2005\ThreeWay\ThreeWayTerrians.gdb" while not arcpy.Exists("strm2M"): time.sleep(10) #print "File not found" print "File found" arcpy.RasterToPolyline_conversion("strm2M", "strm2Mline", "ZERO", 50, "SIMPLIFY") This script bombed at line 14 with a "cannot acquire a lock" and "Error 010213: Error in reading raster Band_1". When I ran the script after the con tool that was making strm2M completed, it hit an error that "strm2Mline exists. Can't overwrite". So the previous run created a new feature class but wasn't able to put any features in it. The attribute table had no records. When I deleted the interrupted feature class strm2Mline, the script ran successfully. Looks like arcpy.Exists doesn't work for this. I need something like arcpy.ReallyCompletedReleasedLockThisTime..

‎03-10-2015

An associate pointed out to me that if you have installed the 64 Bit Background Geoprocessing Tools, the Fill tool is running asynchronously in the background. You can start another geoprocessing tool or script from the ArcMap or ArcCatalog session you used to do the Fill, and it will wait to run until your Fill process has finished. But in my case, I haven't turned on background processing, and, besides, I had started the first script from IDLE, not the python window or tool within ArcMap. Rather than asking the OS about processes, I think all I need to do is use arcpy.Exists by adding something like the following to the top of my second script: import arcpy import time from arcpy import env from arcpy.sa import * arcpy.CheckOutExtension("Spatial") env.workspace = r"G:\elevation\lidar_2005\ThreeWay\ThreeWayTerrians.gdb" while not arcpy.Exists("Fill_TWtoPom12"): time.sleep(10) #rest of the script ...... ...... But I wonder if Exists might jump the gun and evaluate as true before the raster is fully complete. I think I saw an example yesterday of a raster that showed up in Catalog, was able to be added to a map, but had map draw errors because it wasn't completely done. Maybe I can test this later today.

‎03-09-2015

Sorry if I didn't make this clear. Step one was finished, I started step two in the late afternoon and was surprised that it was taking so long. I wanted to put steps 3, 4, 5 into a script, get it in the que somehow so it would wait until the Fill operation completed before starting in on steps 3, 4 , and 5, sometime during the night, so it be all finished in the morning. I had this other idea today that rather than reach into the OS to see if the process was still running, maybe the test could be an arcpy Exists. However, is it possible that a grid can pass the Exist test when it's not really finished? I have a Fill raster now that was able to be added to ArcMap, but will not draw, and it seems like the script that created it is still running.

‎03-09-2015

I was processing some big rasters from Lidar data last week. There were five basic steps I wanted to complete: 1. Terrain to raster 2. Fill. 3. Flow Direction, 4. Flow Accumulation, 5. A Con to select the higher values of the Flow Accumulation to get to stream lines. While step 2 was running, I started thinking "Hey, this is going to take a while. I could put steps 3,4, and 5 in a script or model, and get them ready to run unattended once this Fill completes." But the Fill was still running when I went home. This got me wishing that I could launch the script running the remaining steps, but make the script wait until the Fill completed. Like I used to do in Unix environments. I installed the wmi module for python, http://stackoverflow.com/questions/1632234/list-running-processes-on-64-bit-windows but when I ran it, my process list was empty. They say it runs the same for Windows 7 64 like I am running, but I didn't figure this out. I took a look at wmic http://stackoverflow.com/questions/16326529/python-get-process-names-cpu-mem-usage-and-peak-mem-usage-in-windows . I could get a list of my processes at cmd, so it will work in python, but I didn't figure out a way recognize my Fill process or a way to test for the absence of my Fill process. I suppose the way to do it is to have the script run in circles until the process is not found, then let the script continue, like I did once with Unix shell scripts. Does anyone have an example of how to do this in Windows.

‎03-05-2015

I have some older lidar data that was collected in some irregular shapes. I'm creating some terrains for this data from the multipoint data and would like to avoid triangulation across concave edges by using an outline polygon feature class as a soft clip. What I have been doing is buffering the multipoints by a positive distance, dissolving any internal polygons that were created, then buffering the new polygon by same distance as step 1, but a negative number. The problem is that the first step, buffering the multipoints, can take a long time for big files. Is there a more efficient method? For some of these multipoints files and in cases where I'm going to use several adjacent multipoint files in the same terrain, I think I might find that just digitizing an outline might be more efficient than my buffering method has been.

‎01-22-2015

It's the Identify Tool! In ArcMap, just point at a pixel in the area of interest, and eventually, an Identify box will pop up with a bunch of helpful attributes, including SRC_DATE2 in this case.

‎01-15-2015

I added a photo from arcgisonline to a map that we are using for engineering at an irrigation dam site. It turned out that the image at World_Imagery captures a very illuminating high water event at the site. It is a different image than at ESRI_Imagery_World_2D or Google Earth. I have records at river flows for each day. If I could date the photo, I could say what flow the river was at when the photo was taken.

‎01-07-2015

So you don't want to have the polygon's attribute table contain an item for the presence of each species because that would be 13000 items in the attribute table? What about a link or relate between the data table and the attribute table so a query on the data table will select the polygons that contain the species of interest? But maybe you need to deliver a shapefile, not a map project.

‎12-15-2014

What's the most efficient workflow for fixing small areas in a filled LIdar surface when you find leaks that are sending the stream flows the wrong way? Like a road bed that is blocking flow but I know there is a culvert. In this case, it's the opposite: I want to build a dike to keep flow accumulation or a cost distance line from jumping out of a ditch. Before, did I do something like create a little polygon the shape of my dike, then set the value to the desired top elevation, then use Polygon to Raster to make a dike raster with cell size the same as my fill surface, then run a con something like FixedRas = Con(DikeRas,DikeRas,FillRas) ? Or does it work better to use the dike polygon as a mask and process the fill raster somehow. Or does it work better to run the dike polygon through Polygon to Points, then set the value in the points, then go back to raster< and run the Con? Or is there a better way that I'm not seeing?

‎11-13-2014

Is your problem that you want to copy the point feature class to another directory location, but when you do the photos don't work any more? Is that because the location of the photos specified in the feature class attribute table are relative? If the photos listed in the table have no path to them, the feature class must be expecting them in the same folder as the point feature class. No, I just tested that. Copied a point shape file with associated geotagged photos to a new folder. The photos didn't come along with the copy, the attribute table still shows no path to the geotagged photos, but the photos still show up in ArcMap with Identify tool. How does that work?

Online Status	Offline
Date Last Visited	‎11-11-2020 02:23 AM

My Ideas

Latest Contributions by PaulHuffman

Re: Spatial Analyst Watershed tool

Re: How can I start a script or model on the condition that a previous step has completed normally.

Re: How can I start a script or model on the condition that a previous step has completed normally.

Re: Where is the database requirements matrix for ArcServer 10.3

Where is the database requirements matrix for ArcServer 10.3

Re: How can I start a script or model on the condition that a previous step has completed normally.

Re: How can I start a script or model on the condition that a previous step has completed normally.

Re: How can I start a script or model on the condition that a previous step has completed normally.

How can I start a script or model on the condition that a previous step has completed normally.

Looking for a more efficient way to create a soft clip for Lidar data.

Re: Is there any way to determine the date for a airphoto on World_Imagery on arcgisonline?

Is there any way to determine the date for a airphoto on World_Imagery on arcgisonline?

Re: Python help needed... select polygons from one shapefile, based on values in another table...

Fiddling with a drainage surface

Re: GeoTagged Photos

Can I not mix coordinate systems in Drone2Map?

Re: Can I use MOD in Query Builder?

Can I use MOD in Query Builder?

Re: Why is there no availability or burrow on clie...

Re: How can I start a script or model on the condi...