Latest Contributions by curtvprice

I'm having the same issue with beta 1 using script tools written in vbscript that access the geoprocessor through gpdispatch. You can add me to the list! Still having this problem. Error popup message follows: Could not create object named "esriGeoprocessing.GPDispatch.1". ArcMap 10.0 PR 1 (Build 2348) ... View more

A mole from the geoprocessing team sent me this neat example. This method could be used to generate $$XMAP and $$YMAP using the environment raster extent and cell size. Sure would be nice to include a script tool in the Arc 10 toolbox to both serve as a simple example of how to use NumpyArrayToRaster, and keep the functionality in the standard product. Otherwise the handfull of geeks that use this functionality will have to do it ourselves. Here's the (now trivial, with the example below) tool I will probably write if you don't release one: CreateLocationRaster <Out_raster> {ROWMAP | COLMAP | XMAP | YMAP} {Template_raster} ## NumPyArrayToRaster (in_array, {lower_left_corner}, ##       {x_cell_size}, {y_cell_size}, {value_to_nodata}) import os, sys, arcpy import numpy as NUM # Set path sPath = sys.path[0] outPath = os.path.join(os.path.dirname(sPath),"output") # Set environment arcpy.env.workspace = outPath arcpy.env.overwriteOutput = 1 #------------------------------------------ # Create and populate an array myArray = NUM.zeros((10,16),dtype='int16') # Fill row with row number for i in range(0,10):   myArray = i + 1 print myArray # Create a raster rowRas = arcpy.NumPyArrayToRaster(myArray) rowRas.save("rowRas") print "...saved rowras" #------------------------------------------ # Fill col with col number for j in range(0,16):   myArray[:,j] = j+1 print myArray # Create a raster colRas = arcpy.NumPyArrayToRaster(myArray) colRas.save("colRas") print "...saved colras" And here's another (certifiably cool) way to do it, sent by a USGS colleague... using the NumPy "indices" method:     import numpy as np     import arcpy     # Setup some rasters of rows and columns     # (like old $$ROWMAP and $$COLMAP)     nprows = np.indices((10,10))[0]     npcols = np.indices((10,10))[1]     # Convert the numpy arrays to ESRI rasters (ie Raster objects)     # called 'rows' and 'cols'     rows = arcpy.NumPyArrayToRaster(nprows)     cols = arcpy.NumPyArrayToRaster(npcols) ... View more

*** Do not write below this line *** ESRI Incident #796389 (Barbara S) - Contact - Python diagnostic message is incomplete/wrong when tool name is ambiguous Hello, Curtis, My name is Barbara and I will be assisting you with this incident report. I have duplicated this on 9.3.x, 10B2, and on a weekly build for version 10. I have filed a defect report (#NIM054814) for this behavior: [#NIM054814  Error diagnostics for the Usage geoprocessing command offers a Choices list with no choices.] Thank you for bringing this to our attention. Best regards, Barbara S. ... View more

Any map algebra expression with a dot in my experience is problematic. It's possible you may have better luck if you use Gp.Workspace so you don't need to put the ".gdb" in the path to the raster. ... View more

When I do this at the ArcMap command line (9.4b2): >>> atanb6  = Ln(upareag + 1) / Tan((Raster("nedslp") + 1) * (3.14 / 180.)) I get a layer in my map. If I hit the back button and do it again with a slight change (like I am used to in Grid), nothing happens (apparently). >>> atanb6  = Ln(upareag + 1) / Tan((Raster("nedslp") + 1) * (3.14 / 181.)) I need to change the name of the output: >>> atanb7  = Ln(upareag + 1) / Tan((Raster("nedslp") + 1) * (3.14 / 181.)) This works - a new layer is created in the map with the output. I would expect the layer to be re-written with a new result, not nothing to happen. What's going on, and does this deserve an error message of some type? ... View more

The jury is out - raster calculator has some huge shortcomings in terms of validation and its lack of integration with the rest of the geoprocessing environment. I do hope the raster team can post a video or an early patch of the graphical interface to create command line map algebra, or a video showing how raster calculator type stuff can be done in the python command line window. With the intellisense etc I would think it may be *easier* to use than the raster calculator. We will see though.... I need to play with it a bit before I decide whether this is a net negative. The biggest issue is the loss of a map algebra tool within modelbuilder, but there are workarounds. May I suggest you post a use-case of a raster calculator expression you would normally do and maybe someone from the raster team can walk you through how to do it in 9.4? ... View more

That was it! Those links work.  Thanks! ... View more

I don't seem to have access to helpdev.esri.com. Is there another site that is accessible to beta testers outside ESRI? ... View more

I had a nice discussion with a product team member (who shall remain anonymous) about this a couple of weeks ago. He suggested getting conversant with the numpy raster object -- arcpy will have hooks to easily go back and forth, which may be a really handy way to create function-based rasters that don't lend themselves easily to map algebra, or are slower to manipulate within the map algebra environment. Is this "hook" to numpy rasters going to be available in beta 2? If not what is a workaround to try in the meantime? ... View more

I'm having the same issue with beta 1 using script tools written in vbscript that access the geoprocessor through gpdispatch. You can add me to the list! ... View more