topic Re: Zonal stats: fail to open raster dataset for a multi-band tiff. in Python Questions

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

RichardThurau — Tue, 06 Dec 2011 14:45:38 GMT

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

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

ChrisSnyder — Tue, 06 Dec 2011 15:21:26 GMT

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")

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

RichardThurau — Tue, 06 Dec 2011 15:26:12 GMT

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

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

ChrisSnyder — Tue, 06 Dec 2011 15:28:05 GMT

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?

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

RichardThurau — Tue, 06 Dec 2011 16:20:39 GMT

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

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

ChrisSnyder — Tue, 06 Dec 2011 16:35:54 GMT

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?

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

RichardThurau — Tue, 06 Dec 2011 17:13:28 GMT

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?

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

ChrisSnyder — Tue, 06 Dec 2011 19:22:36 GMT

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.

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

RichardThurau — Tue, 06 Dec 2011 21:11:45 GMT

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

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

Luke_Pinner — Wed, 07 Dec 2011 23:26:31 GMT

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.

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

RichardThurau — Thu, 08 Dec 2011 01:56:07 GMT

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

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

KatSuda — Sat, 08 Dec 2012 14:30:33 GMT

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

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

RichardThurau — Sat, 08 Dec 2012 15:03:00 GMT

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

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

KatSuda — Sat, 08 Dec 2012 15:25:58 GMT

Thanks, Rich. I'll post the question.

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

SalmanAhmed3 — Mon, 23 May 2016 13:41:50 GMT

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.