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; } } }

Zonal stats: fail to open raster dataset for a multi-band tiff.

RichardThurau

Hi,
Trying to extract some values from some tiny AOI polygons from multiple bands of a .tif.

AOI = "X:\\aoi.shp"
Img1_1 = "X:\\Glovertown_in_1.tif\\Band_1"
tab_1_1 = "X:\\Tab1_1"

Tab1_1 = ZonalStatisticsAsTable(AOI, "OBJECTID", Img1_1, "tab_1_1", "", "MINIMUM")

error:
ExecuteError: ERROR 999999: Error executing function.
Failed to open raster dataset
Failed to execute (ZonalStatisticsAsTable).

Can anyone tell me what's going wrong here?

Thanks

Rich

タグ付きの投稿をもっと見る

Accepted answers

All comments

ChrisSnyder

Try

Tab1_1 = ZonalStatisticsAsTable(AOI, "OBJECTID", Img1_1, tab_1_1, "", "MINIMUM")

instead of

Tab1_1 = ZonalStatisticsAsTable(AOI, "OBJECTID", Img1_1, "tab_1_1", "", "MINIMUM")

RichardThurau

Chris,
I tried that as well, right after I posted. No Dice.

The tools says it supports floating rasters for the minimum statistic. Hoping not to have to create a bunch of extra output just to convert to int, which is my current hunch.

Thanks for the reply.

Rich

ChrisSnyder

What if you convert band1 of your .tif to a grid... Does it work then? Maybe ZonalStats doesn't like using multiband .tif files as input?

RichardThurau

Chris,
That is possible but this is a small piece of code for a larger project.
1. If I have to convert everything to grid, it's probably a deal breaker. I'll just work it out in SML in ERDAS. I'm trying to use python because it would allow me to automate a bit further down the line.

2. I've seen other people (code examples in these forums) work successfully with .tif

Thanks

RT

ChrisSnyder

I ran a little test, and I was able to use Band_1 from a multiband .tif in the Zonal Statistics as Table tool via ArcToolbox.

arcpy.ZonalStatisticsAsTable(roadshedFC, "ROADSHED_ID", r"C:\temp\test.tif\Band_1", outZnStatsTbl, "", "ALL")

Maybe it's just the map algebra is wonkey...

Does the toolbox tool work for you?

RichardThurau

Thanks for your persistence Chris,
Yes, of course (sarcasm) the tool runs fine. Unfortunately, this isn't going to be very useful in the grand scheme.
Here's the tool, python snipet:

arcpy.ZonalStatisticsAsTable_sa("X:/aoi.shp","OBJECTID","X:/image_in_1.tif/Band_1","x:/zonalst_glovert1","DATA","MINIMUM")

But we have made progress. What's wrong with this code?

ChrisSnyder

Try using a different field than OBJECTID from your AOI polygon. Used to be that OBJECTID wasn't an allowed field for any SA tools. Now it seems to be okay, but sometimes the ESRI tools revert to past behaviors...

What if you use the "ALL" default for the statistic type? Better yet, just use "" to force the default value.

RichardThurau

Chris, no dice.
Trying default values as you suggested. Created an IMG, didn't work. Brought the IMG into a GDB, didn't work. All the same error "failure to open raster dataset".

Where is Python disconnected from the tool? I had a similar problem with Mosaic a while back, where python wouldn't run, but the tool does. I ended up having to do something completely different just to work around.

Anyone have any foundational wisdom here?

Chris probably needs some relief! I know I do.

RT

Luke_Pinner

I've had trouble with ZonalStats using objectid fields previously. Try using another field of SmallInteger or Integer type, no floating points (i.e. Double) or text.

You also need to make sure that the ArcGIS ScratchWorkspace environment setting (or the Windows %TEMP% environment variable) points to a path without spaces as ArcGIS writes temporary grids there and will fail with a very unhelpful error message if the path has spaces.

RichardThurau

Luke,
Thanks. I did create an ID field that did not work. I also passed the image layers to Int() prior to opening. I've tried on multiple images in the area, using .tif and .img file types and different bands.

I'll have to check my scratch workspace tomorrow. Thanks. Not sure why that would work in a model and not python though. Maybe I'll post a new thread with that difference in the title.

RT

KatSuda

Richard, did you ever find what the problem was?

I am having similar trouble trying to composite bands. I've tried using them as tiffs and as grids, but I get error 000271: Cannot open the input datasets.

Kat

RichardThurau

Kat,
Honestly, I don't remember if this was ever resolved. I'm thinking it was not, as I use tab areas to calculate my own stats usually (although zonal stats in 10.1 is superior because to calculates partial pixel values). I'd suggest you start a new thread and post your version info and code.

This community is pretty great at either giving an answer quickly, or, if they can't, you will know quickly that you have to work around.

Best

Rich

KatSuda

Thanks, Rich. I'll post the question.

SalmanAhmed3

Just don't change the ArcGIS default output location when running the tool. You can change the name of the output file though. But let it save in the default directory given by ArcGIS (somewhere in User-Documents-ArcGIS-...). Works for me with this trick.

EricIronside_1-1622050916668.png

dave583.jpg

会話に参加する

サインインしたメンバーは投稿、更新のフォローなどができます。初めてですか？無料アカウントを登録してください。

クイックリンク

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