nav[aria-label="Primary Navigation"] { padding: 0; & ul { list-style: none; width: 100%; display: flex; flex-direction: row; justify-content: start; align-items: start; gap: 30px; padding: 0; & li { margin: 0; } & ul li { list-style: none; } } }

calculate zonal statistics using overlapping polygons

HannahPomella1

Hello! I'm trying to calculate zonal statistics with a grid as input and overlapping circles in a feature class as areas. As I understud the software is converting the feature class with the polygons in a raster file in the background and clips the polygons. Can I avoid this somehow or is there an other possibility to extract this spatial information from the grid? What I would need is the number of cells with specific values within a 2000m radius around points. Thanks!

Encuentra más publicaciones etiquetadas con

Accepted answers

DarrenWiens2

What type of feature class are you working with? Where clauses are data-specific (see the help for Make Feature Layer, where clause section). My example used a shapefile - fields are enclosed by double quotes. File GDB fields have no quotes. Personal GDB fields are enclosed by []. A nice way to figure it out is to run Make Feature Layer once through the normal tool dialog to figure out the syntax. Once you get it to work once, you can export the result from the results window to a Python snippet and inspect the syntax more fully in a text editor or IDE

All comments

ChrisDonohue__GISP

So if I understand your dilemma correctly, the issue is the overlapping polygons (circles) not being processed correctly by the Zonal Statistics tool. If that is the case, the help for the tool offers this explanation and resolution for this issue:

If the zone feature input has overlapping polygons, the zonal analysis will not be performed for each individual polygon. Since the feature input is converted to a raster, each location can only have one value.
An alternative method is to process the zonal statistics iteratively for each of the polygon zones and collate the results.

Source: ArcGIS Help 10.1

So the solution then would involve iterating the Zonal Statistics tool for each individual polygon. Common ways to accomplish this include Python and Modelbuilder (using an Iterator).

Chris Donohue, GISP

(EDIT: fixed typos]

DarrenWiens2

In addition to Chris' comment about processing zones individually, it sounds like you may be more interested in the Tabulate Area tool (which will list the cells in each class), rather than Zonal Statistics (which performs some numerical summarization of your classes, like mean, range, etc.).

curtvprice

Esri has provided an enhanced Zonal Statistics and Tabulate Area tools that handle overlapping polygons in the Spatial Analyst Supplemental tools.

I have posted a toolbox on ArcGIS Online that does both zonal statistics of tabulations of categorical data from points. (It does the buffers for you). It also includes a tool that will weed points if you don't want to run them one at a time, which can be slow for thousands of points:

NAWQA Area-Characterization Toolbox

For a lot more discussion on this, see:

Re: Zonal Statistics as Table: Overlapping Polygons: Solution Required

Hope this is helpful!

HannahPomella1

hi, thank you very much for your aswers and sorry for my late replay.

@Darren Wiens You are right, i need the Tabluate Area tool to get the results I want. But there is still the problem with the overlapping polygons. For sure I can calculate single circle areas and summarise then the results in one table, but I have to calculate >1000 points, so a tool would be nice to have.

@Curtis Price thank you for your tipps. I tried now the spatial analyst supplemental tool on ArcGIS 10.3.

If I use the Tabulate Area 2 tool on three not overlapping circles it's working perfectly. Then I tried it on three overlapping circles and it does the same as the standart Tabulate Area tool. Its not calculating the entire area. Maybe somebody has a hint what could be the problem?

I tried also the NAWQA Area-Characterization Toolbox, but its alwas producing a error "the extent my input featurs is not exact"... No idea what could be wrong with them.

Thanks!

Hannah

DarrenWiens2

Here is how you can do it with Python. Requires Spatial Analyst. This example uses an overlapping zone feature class called "buffer_1000" (must have "FID" field), class raster called "raster" (must have "Value" field), and outputs one table for each feature in the zone feature class called "outTable1", "outTable2", etc. You can then run Merge to consolidate these output tables. Of course, change the feature class names to suit your data.

>>> zoneFC = "buffer_1000"
... raster = "raster"
... with arcpy.da.SearchCursor(zoneFC,"OID@") as cursor:
...    for row in cursor:
...        arcpy.MakeFeatureLayer_management(zoneFC,"tempFL","\"FID\"=" + str(row[0]))
...        rast = arcpy.sa.TabulateArea("tempFL","FID",raster,"Value",r"in_memory\outTable" + str(row[0]))

HannahPomella1

Thank you very much for the code.

I changed the file names but he reclaims a syntax error in row 6

>>> zoneFC="POIzone_Random"

... raster="Allminerals_Extension.bmp"

... with arcpy.da.SearchCursor(zoneFC,"POI_Nr") as cursor:

... for row in cursor:

... arcpy.MakeFeatureLayer_management(zoneFC,"tempFL","\"POI_Nr\"=" + str(row[0]))

... rast=arcpy.sa.TabulateArea("tempFL","POI_Nr",raster,"Value",r"E:\Projekte\Coop_Meyer\outTable" + str(row[0]) >>>

...

Parsing error SyntaxError: invalid syntax (line 6)

DarrenWiens2

There was a copy/paste error in there - I've updated above. The final line should end with this, not >>>: str(row[0]))

rast = arcpy.sa.TabulateArea("tempFL","FID",raster,"Value",r"in_memory\outTable" + str(row[0]))

HannahPomella1

Hi!

I was using a personal GDB. Now I export the data to shapefile and run the code again.

Its working at least a bit 🙂

It creates a table with one entry and then stops with a error.

Únete a la Conversación

Los miembros registrados pueden publicar, seguir actualizaciones y más. ¿Nuevo aquí? Regístrate gratis.

Enlaces Rápidos

Community Help
Find useful guides, FAQs, and documents to help you navigate and make the most of Esri Community.