Latest Contributions by curtvprice

‎05-11-2011

The two nearest cells would be 5 x 5 neighborhood: o o o o o o o o o o o o x o o o o o o o o o o o o a 2 x 2 neighborhood probably looks more like this: o o x o

‎05-11-2011

Thanks, that was most of what I needed. First, though I had to compute which pixels where >30% slope. I did this with the raster calculator which I eventually found hidden under Map Algebra. I was then able to use the Tabulate Area under Zonal tools to finish the job. Here's another approach: 1. Tag high-slope pixels with the polygon ID, others NoData (Raster Calculator or Single Output Map Algebra) (Set the extent, snap, and raster processing cell size to elevgrid before executing) SetNull(slope(elevgrid,"DEGREES") <= 30,polygrid(inpoly,polyid)) 2. calculate area of each zone with the Zonal Geometry As Table tool

‎05-11-2011

Does it matter if rasters have a different pixel type or pixel depth when using the raster calculator? My experience is that no, it doesn't, the intermediate data sets and output will use a data-appropriate pixel depth. If you want to reliably force the output pixel depth of a raster, the only method I know of that works, is to copy the raster with the Copy Raster tool.

‎05-02-2011

Anyone... anyone...?

‎04-29-2011

Which version of ArcGIS?

‎04-29-2011

Does anyone know if the following docell code can be done with map algebra in ArcGIS 10? docell if (flr_ezg eq 1) to_flg_lg = flg_lg_ezg(1,0) if (flr_ezg eq 2) to_flg_lg = flg_lg_ezg(1,1) if (flr_ezg eq 4) to_flg_lg = flg_lg_ezg(0,1) if (flr_ezg eq 😎 to_flg_lg = flg_lg_ezg(-1,1) if (flr_ezg eq 16) to_flg_lg = flg_lg_ezg(-1,0) if (flr_ezg eq 32) to_flg_lg = flg_lg_ezg(-1,-1) if (flr_ezg eq 64) to_flg_lg = flg_lg_ezg(0,-1) if (flr_ezg eq 128) to_flg_lg = flg_lg_ezg(1,-1) end I would like to code in arcobjects c# with map algebra. I have read about numpy array as one way to do it in python. Thanks for any help Bruce Vandenberg Bruce, Neighborhood notation is not supported at Arc 10. We've been asking. Please do contact tech support on this. They may send you to using numpy rasters to do this kind of thing. Another workaround may possibly to use the MultiOutputMapAlgebra tool, but it isn't supported at 10.0 although it is still there in the gp object. You may be able to work out something using the Shift tool to offset a bunch of rasters and then use Pick to select the right one from the list. You can convert the powers of two to 1...8 using this map algebra expression. This could be used with Pick. (Log2(flr_ezg) + 1) Hope this helps you out!

‎04-18-2011

What was the problem with the union processing -- did the tool fail, or did the result not have the results you needed? Have you considered doing the analysis in the raster domain using the Combine tool? This tool would give you a table of every unique combination of inputs at each cell which can then be used with the Pivot Table tool (or exported to a database or spreadsheet) for summarizing the combinatorial analysis results.

‎04-18-2011

I'm wondering if that was a pathname issue. BTW the best way to strip off the extension is to use the python os module. This will work even if the input file name doesn't include an extension: import os ascOut = os.path.splitext(ascIn)[0] + ".img"

‎04-15-2011

It's also kind of dangerous to name datasets or fields anywhere in arcgis that start with a number.

‎04-15-2011

Does this tool require an ArcInfo level license? Yes unfortunately most of the tools do require an arcinfo license. I have a version that does not require an ArcInfo license for as many of the tools that I am testing. I have successfully got it running with versions 9.3 SP 2, 9.3.1 SP 2 and 10.0 SP 1. Let me know if you'd like to test it out -- my email is cprice ------ at ------ usgs.gov

‎04-14-2011

Are you thinking a series of embedded CON statements? That's exactly what I was thinking. You could write python code to generate the map algebra needed. Now that I think about it, I think the approach I'm suggesting is way too complex. If these are integer codes, why not run a combine() to make a temp grid with all unique combinations in the VAT, add a field to the VAT, and then use Calculate Value to populate your output using a Python or VBScript function, then use the Lookup tool (equivalent to the old grid.item syntax) to create your output grid. (Python would be neater because you can use dictionaries and/or lists to elegantly do your lookup.)

‎04-12-2011

Jamal, the best you can do is to create a point feature at the point you are interest in, populate it with the raster value using the Extract Values to Points tool, and then label the point. Hope that helps.

‎04-12-2011

It is best to use NoData to represent NoData. Any raster operations you do should handle NoData correctly and your statistics will not include the NoData areas. So, yes, use SetNull to set your -9999 cells to NoData and use that grid for the mask, extent and snap raster for future operations. If you want to dump the raster out for use outside ArcGIS, you can either convert the NoData using the IsNull function to a value, or the conversion too (Grid To Ascii) often allows you to specify a value for NoData cells.

‎04-12-2011

I'm looking for a geoprocessing tool that would this: given this value in raster x and this value in raster y, raster z (output) gets this value. I'm thinking something like "Reclass using Ascii file" but using a matrix of two paramaters maybe? Any ideas? This sounds like something that could be pretty handily accomplished using map algebra. You could use the Calculate Value tool to build a map algebra expression and pass it to the single-output map algebra tool (9.3.1) or the new map algebra tool (10.0). You could get pretty fancy building the expression by writing a function in python in the code window of the Calculate Value tool.

‎04-12-2011

The intermittent problem I was having with unicode was occurring only in modelbuilder. The problem path strings there are created by the system and the error would occur while the Calculate Value tool was reading the string -- so there was no place to work around it in code. It was a very intermittent problem, fixed by resetting my HKCU registry hive, so I'm very happy that the ESRI support folks were able to reproduce it.

Online Status	Offline
Date Last Visited	‎06-19-2024 12:10 AM

My Ideas

Latest Contributions by curtvprice

Re: Focal Statistics sum

Re: How to compute % area above 30% slope for 100+ polygons in layer

Re: Pixel type and depth in raster calculator

Re: Frequency missing a "case" field

Re: Normalizing Values within Rasters using Raster Calculator

Implementing neighborhood notation in Arc 10?

Re: A grid applied in a sum of polygons included(how to analize polygon in polygon)!

Re: ERROR 000865: Input ASCII raster file: 299_327_Intensity.asc does not exist.

Re: ERROR 000865: Input ASCII raster file: 299_327_Intensity.asc does not exist.

Re: Tabulate Area of Overlapping Polygons

Re: Any geoprocessing "Reclassify" type tool that uses a matrix?

Re: is there a way to label the cells of the raster with theirs values?

Re: Setting "-9999" as a 'NoData' Value (Changes Minimum Value)

Re: Any geoprocessing "Reclassify" type tool that uses a matrix?

Re: Calculate Value tool fails, unicode error from paths with "\U" (rawunicodeescape)

Re: ModelBuilder SQL Expression Using Inline Varia...

PyCharm Setup for ArcGIS Desktop

Re: Error 010240 - Could not save raster dataset.....

ArcInfo Workstation Commands

Re: Is it possible to configure SA in a Python scr...

ArcGIS Desktop Installation Support

R-ArcGIS

Coordinate Reference Systems

Python Snippets