topic How to determine whether Raster statistics have been calculated in Python Questions

How to determine whether Raster statistics have been calculated

MichaelStjernholm — Fri, 15 Apr 2011 09:03:51 GMT

Hi

I have a script which fails when I use the GetRasterProperties_management function on a dataset where statistics have not yet been calculated.

How can I in a simple way check whether the statistics have been calculated ?

There is strangely enough not a Raster set attribute named "hasStatistics" or similar to get from Describe fuunction
neither does the GetRasterProperties simply rerturn a Null object
so I could use something like :

if arcpy.GetRasterProperties_management (outRaster, "MINIMUM") == None:
CalculateStatistics_management (outRaster)

All help is most welcome

Michael

Re: How to determine whether Raster statistics have been calculated

HongmeiZhu — Fri, 15 Apr 2011 14:44:17 GMT

After raster statistics are calculated, a file with aux.xml extension is created in the same directory as the raster resides. The auxiliary file saves information of statistics for future use. Therefore, if you can see an aux.xml file in the same location, that indicates Raster statistics have been calculated. Or you can simplly put a line after CalculateStatistics_management( ) to inform the calculation has been done successfully.

Print "Done successfully!"

Re: How to determine whether Raster statistics have been calculated

MichaelStjernholm — Mon, 18 Apr 2011 06:42:16 GMT

Thanks Hongmei for reminding me of using the test of whether the aux.xml file exists - I use similar approach in other calculations.

regards michael

Re: How to determine whether Raster statistics have been calculated

NiklasNorrthon — Fri, 10 Dec 2021 22:55:36 GMT

Another way to solve the problem is to just assume that statistics is there and deal with the problem when it isn't.

# Pseudocode:
while True:
    try: 
        arcpy.FunctionThatNeedsRasterStatistics(raster)
        break
    except arcpy.ExecuteError:
        arcpy.CalculateStatistics(raster)

This pattern also has the advantage that it isn't subject to race conditions (not that those are likely for raster statistics, but in other contexts that can be a real issue.

Re: How to determine whether Raster statistics have been calculated

JakeSkinner — Fri, 10 Dec 2021 22:55:38 GMT

I just wanted to note that an aux.xml can exist when statistics do not. For example, the aux.xml will also be used as a pointer to the pyramid file. This is explained here:

http://help.arcgis.com/en/arcgisdesktop/10.0/help/index.html#//009t00000027000000.htm

You may want to go along the lines Niklas mentioned. Ex:

try:
    arcpy.GetRasterProperties_management(input, "STD")
    print "Statistics exist"
except arcpy.ExecuteError:
    arcpy.CalculateStatistics_management(input)
    print "Statistics were calculated"