topic Re: Export raster statistics? in Data Management Questions

Export raster statistics?

StefanHaglund1 — Fri, 08 Oct 2010 12:38:30 GMT

Hi!

Does anyone know if there is a tool out there that lets me export the statistics of a raster datatset?

Basically, I would like to write out the statistics of several rasters to a table or textfile. Think of how it looks under the Statistics heading in the Rater Dataset Properties window.

There is the Get raster properties tool, which I could call several times to get most of what I want, but I was hoping for something smoother than that.

I could try and come up with a script that does this but I've looked through the geoprocessor and can't seem to find anything useful there either.

Any suggestions?
Nooski

Re: Export raster statistics?

RobertBerger — Fri, 08 Oct 2010 15:12:37 GMT

Hi Nooski,

I don't know of any other tool/script that easily lets you export the statistics. In ArcMap you can export the current statistics to an xml file (when using statistics from current display extent) but I don't think there is a way to script this.
Get Raster Properties might be your best bet.

Robert

Re: Export raster statistics?

Luke_Pinner — Sat, 11 Dec 2021 19:48:46 GMT

Something like:

# GetRasterProperties.py 
# Description: Write raster properties to a file. 
# Note: Written for ArcGIS 9.3. 

# Create the Geoprocessor object 
import arcgisscripting
gp = arcgisscripting.create()

# Get input parameters
InRaster = gp.GetParameterAsText(0)
OutFile = gp.GetParameterAsText(1)

#Property Types
PropertyTypes=['MINIMUM','MAXIMUM','MEAN','STD','UNIQUEVALUECOUNT','TOP','LEFT','RIGHT','BOTTOM','CELLSIZEX','CELLSIZEY','VALUETYPE','COLUMNCOUNT','ROWCOUNT','BANDCOUNT']

#Output file
OutFile=open(OutFile,'w')

# Process: GetRasterProperties 
for PropertyType in PropertyTypes:
    try:
        Property = gp.GetRasterProperties(InRaster, PropertyType)
        gp.AddMessage('%s=%s'%(PropertyType,Property))
        OutFile.write('%s=%s\n'%(PropertyType,Property))
    except:
        # Print error message if an error occurs 
        gp.AddMessage(gp.GetMessages(2))

OutFile.close()

Re: Export raster statistics?

BryceStath — Tue, 21 Dec 2010 18:59:03 GMT

This could be really useful if it writes the statistic I need (mean) for many rasters to one file.

How would I run this? I've tried running it stand alone in PyWin, no dice, and by creating a script within Toolbox within an ArcMap session that has the rasters I want to generate stats for.

thanks

Re: Export raster statistics?

StefanHaglund1 — Sat, 11 Dec 2021 19:48:48 GMT

All,

Here's a follow up.

With help from Luke I managed to get a crude script for 10. It looks like this:

import arcpy
from arcpy import env
env.workspace = "C:\Path"
outfile = open("C:\RastStats.txt", 'w')
PropTypes = ["MINIMUM", "MAXIMUM", "MEAN"]
rastlist = arcpy.ListRasters("*", "")
for raster in rastlist:
    rastname = str(raster)
    outfile.write(rastname + "\n")
    for PropType in PropTypes:
        print PropType
        Prop = arcpy.GetRasterProperties_management(raster, PropType)
        textinfo = (PropType, Prop)
        textstring = str(textinfo)
        outfile.write(textstring + "\n")
outfile.close()

There is still some formatting to do, the result looks like this:
--------
raster1
('MINIMUM', <Result '0'>)
('MAXIMUM', <Result '2203.59692382813'>)
('MEAN', <Result '400.089714329223'>)
--------

The lines ('MINIMUM', <Result '0'>) and so on is what comes out of the GetRasterProperties tool in the script, it would be good to have it formatted nicer, like:
--------
raster1:
MINIMUM = 0
MAXIMUM = 2203.59692382813
MEAN = 400.089714329223
--------

What I have now is enough for me at the moment, but I wonder why ESRI (or arcpy) formats the results like it does? If you run it as a tool dialog you get "MEAN = 400.089714329223".

Also, the script is pretty slow, maybe the nested for-loops slows things down?

Re: Export raster statistics?

BryceStath — Sat, 11 Dec 2021 19:48:51 GMT

All,

Here's a follow up.

With help from Luke I managed to get a crude script for 10. It looks like this:

import arcpy
from arcpy import env
env.workspace = "C:\Path"
outfile = open("C:\RastStats.txt", 'w')
PropTypes = ["MINIMUM", "MAXIMUM", "MEAN"]
rastlist = arcpy.ListRasters("*", "")
for raster in rastlist:
    rastname = str(raster)
    outfile.write(rastname + "\n")
    for PropType in PropTypes:
        print PropType
        Prop = arcpy.GetRasterProperties_management(raster, PropType)
        textinfo = (PropType, Prop)
        textstring = str(textinfo)
        outfile.write(textstring + "\n")
outfile.close()

This looks promising, but I'm still on 9.3
I'll be trying to get this to work, but could use help if anyone feels like translating this.

Re: Export raster statistics?

Luke_Pinner — Thu, 03 Feb 2011 20:36:56 GMT

This looks promising, but I'm still on 9.3
I'll be trying to get this to work, but could use help if anyone feels like translating this.

See the 3rd post - (http://forums.arcgis.com/showthread.php?p=45823#poststop) it's for 9.3 and will do what you want with minimal changes. Basically instead of looping through statistic types, loop through rasters and get the mean - ie mean = gp.GetRasterProperties(someraster, 'MEAN')
.

Re: Export raster statistics?

BryceStath — Sat, 11 Dec 2021 19:48:54 GMT

I've been working on this (see code below).... and it runs although it doesn't want to produce an output file for each grid within gp.Workspace = "C:\\working\\Temp\\testing\\grids

Lets say I have two sample grids in there
f1p060_10c2
f1p060_10c4

I end up only getting in my output folder
C:\working\Temp\testing\outputs3\f1p060_10c2_stats.dbf

So it appears that my while loop isn't looping.
The end goal here is to have all of my raster stats in one file. If I can get this process to work, to write the stats from each input raster to one file, that's good enough. I have an excel script that will push all files in one directory into one master file.
Thanks for any help.

# BatchZonalStatisticsAsTable.py
# Description: Summarizes values of a raster within the zones of another dataset and reports the results to a table
# Requirements: None
# Author: David Coley
# Date: 05/10/06

#import win32com.client, sys, os, string
import arcgisscripting
gp = arcgisscripting.create()
gp.CheckOutExtension("Spatial")
gp.CellSize = 1 #can set as argument
gp.OverWriteOutput = 1

try:
    gp.Workspace = "C:\\working\\Temp\\testing\\grids" 
    InZonalLayer = "C:\\working\\Temp\\testing\\TreatmentFranksTract2010.mdb\\Controls"
    InZonalValueField = "PIN" #can set as argumment
    #set output workspace which holds tables in a gdb
    OutWorkSpace = "C:\\working\\Temp\\testing\\outputs3"
    #Get the raster datasets in the input workspace and loop through them from the start

    InputRasters = gp.ListRasters()
    InputRasters.reset()
    InputRaster = InputRasters.next()
    OutTable = "C:\\working\\Temp\\testing\\outputs3\\" + InputRaster + "_stats.dbf" 

    while InputRaster:
     
        # Process: Zonal Statistics as Table 
            gp.ZonalStatisticsAsTable_sa(InZonalLayer, InZonalValueField, InputRaster, OutTable, "DATA")
        # print InputRaster
            InputRaster = InputRasters.next()
        #print OutTable
except:
    gp.AddMessage(gp.GetMessages(2))
    print gp.GetMessages(2)

Re: Export raster statistics?

ArvindBhuta — Tue, 23 Aug 2011 14:31:25 GMT

What parameters must I input in the parameters tab to get this to work? I am new to scripts and getting them to work.

Re: Export raster statistics?

curtvprice — Tue, 17 Jan 2012 16:47:32 GMT

What parameters must I input in the parameters tab to get this to work? I am new to scripts and getting them to work.

This script is hard coded - it does not take parameters.

Re: Export raster statistics?

curtvprice — Sat, 11 Dec 2021 19:48:57 GMT

So it appears that my while loop isn't looping.

You skipped the first record, so you're only getting the second.

    InputRasters = gp.ListRasters()
    ## InputRasters.reset()  NOT NECESSARY, the enumeration is new
    ## InputRaster = InputRasters.next()  # DELETE - this skips to the second record