multiplying rasters range problem

arcgisarcgipro_6

Hello 

I multiply 5 different rasters in raster calculator. However, I can not get the minimum and maximum range as it expected mathematically from the output raster layer. I would appreciate if you post the solution here. 

Thanks,

Gigi

danretired

That is because no location/cell (ie in your stack of 5 rasters) has a combination of 5, 5, 5, 5, 5

In order to get 3125 all the stacked rasters have to have a value of 5 at one location

a = np.array([[0, 1, 2], [3, 4, 5]])  # ---- one raster
a5 = np.vstack([a, a, a, a, a]).reshape(5, 2, 3)  # 5 identical stacked rasters
a
array([[0, 1, 2],
       [3, 4, 5]])

a5
array([[[0, 1, 2],
        [3, 4, 5]],
       [[0, 1, 2],
        [3, 4, 5]],
       [[0, 1, 2],
        [3, 4, 5]],
       [[0, 1, 2],
        [3, 4, 5]],
       [[0, 1, 2],
        [3, 4, 5]]])

np.cumprod(a5, axis=0)  # ---- their product looking down through a "cell"

array([[[   0,    1,    2],
        [   3,    4,    5]],

       [[   0,    1,    4],
        [   9,   16,   25]],

       [[   0,    1,    8],
        [  27,   64,  125]],

       [[   0,    1,   16],
        [  81,  256,  625]],

       [[   0,    1,   32],
        [ 243, 1024, 3125]]], dtype=int32)