topic Re: retrieve raster statistics and export to an .xls in Python Questions

retrieve raster statistics and export to an .xls

HeatherHook — Tue, 22 May 2018 15:49:37 GMT

I need to retrieve minimum, maximum, mean, and standard deviation statistical values from a directory of raster files and store the values in an excel spreadsheet.

My approach as been cursors and GetRasterProperties_management.

Any insight would be greatly appreciated.

Re: retrieve raster statistics and export to an .xls

DanPatterson_Retired — Tue, 22 May 2018 16:06:00 GMT

Get Raster Properties..

returns a result object, but what is your code that provides the shell around dealing with the result and putting it into a format that you can use for Table to Excel ?

Re: retrieve raster statistics and export to an .xls

HeatherHook — Tue, 22 May 2018 16:19:09 GMT

That's a very good question. I was using a text document file with each object separated by a coma space.

Re: retrieve raster statistics and export to an .xls

DanPatterson_Retired — Tue, 22 May 2018 16:29:33 GMT

actually I was referring to the code that cycles through the rasters and uses the tools's arcpy code access to get the stats. Or do you not have that bit yet? since you refer to cursors

Re: retrieve raster statistics and export to an .xls

HeatherHook — Tue, 22 May 2018 16:46:38 GMT

for loop. but I really haven't gotten that bit to work properly yet.

Re: retrieve raster statistics and export to an .xls

HeatherHook — Tue, 22 May 2018 16:49:27 GMT

i really feel like there is a more efficient approach I just have no idea what it could be. I also thought there would already be an existing tool for this objective but my searches have come up empty handed.

Re: retrieve raster statistics and export to an .xls

Anonymous User — Wed, 23 May 2018 07:15:17 GMT

I did something similar the other day. Instead of keeping it all inside GIS, I did a standalone script and exported from the script to a .txt file, which I then imported in Excel.

The iterating and printing part of the script looks like this:

#!/usr/bin/env python
# -*- coding: utf-8 -*-

import arcpy, os, sys, string, glob

arcpy.env.overwriteOutput = True

shpWS = "E:\MidlGIS\I7010"

arcpy.env.workspace = shpWS

shplist = arcpy.ListFeatureClasses()

print (len(shplist))

for fc in shplist:
count = arcpy.GetCount_management(fc)
ord = fc
ord = ord[6:]
ord = ord[:-4]
print fc[1:5]+"____"+ ord +'____'+str(count)