I am using ArcMap 10.4 (with Spatial Analyst.)
I have a feature layer with ~11.000 partly overlapping but similar features. For each feature I want to extract some values from ~25 underlying rasters.
I set up a model that iterates the raster folder and applies the needed values to the features. It works as intended, except for Zonal Statistics as a table doesn't handle overlapping features.
I have tried the following in order of solving this problem:
- Set up the iterators in submodels, but they only run once. In order of fixing that I tried making them preconditions for the forthcoming tools. That didn't solve produce intended result.
- I tried simply deleting every second feature, by
MOD(layerFID,2)=1, which was close to correct, but isn't really a solution.
- I set up an iterator that handles the overlapping features, by splitting the into separate shapefiles, performs Zonal Statistics as a table on them, and merges them again. It works, despite of being very slow. And then I can't iterate the raster folder.
- I tried running Zonal Statistics 2 (and 3), with and without "arcpy.env.overwriteOutput = False". I yields odd results. As far as I know COUNT should be the same for all the zones when the area and the raster is the same.
I don't see any decent solutions that doesn't involve two iterators.