>>> raster_as_numpy_array = arcpy.RasterToNumPyArray('C:\Users\xxxxx\Downloads\q47121g52be.tif\q47121g52be.tif')
>>> raster_geometric_mean = scipy.stats.mstats.gmean(raster_as_numpy_array)
>>> raster_geometric_mean
array([ 0., 0., 0., ..., 0., 0., 0.], dtype=float16)
So I tried the above on integer and float rasters and am not sure what I did wrong there...how do you get to the gm as a single reported value? If you can't tell I am fairly new at this point to Python -thanks you guys for all your help!