nav[aria-label="Primary Navigation"] { padding: 0; & ul { list-style: none; width: 100%; display: flex; flex-direction: row; justify-content: start; align-items: start; gap: 30px; padding: 0; & li { margin: 0; } & ul li { list-style: none; } } }

Getting the Mean value out of the Max values of a raster data set

anTonialcaraz

Hi,

As the title says, I'm trying to calculate in python the mean value out of all maximum values in a list of rasters (696) in a GDB.

First thing would be to calculate all max values but I'm not sure how to "hold" all those values to then calculate the mean.

Any help will be greatly appreciate it.

Running ArcGis 10.5.1

Finde mehr Beiträge mit dem Tag

Accepted answers

JakeSkinner

Change int to float:

<SPAN class="comment token"># Iterate through each raster and append max val to list</SPAN>
<SPAN class="keyword token">for</SPAN> raster <SPAN class="keyword token">in</SPAN> arcpy<SPAN class="punctuation token">.</SPAN>ListRasters<SPAN class="punctuation token">(</SPAN><SPAN class="string token">"*"</SPAN><SPAN class="punctuation token">)</SPAN><SPAN class="punctuation token">:</SPAN>
&nbsp;&nbsp;&nbsp; val <SPAN class="operator token">=</SPAN> arcpy<SPAN class="punctuation token">.</SPAN>GetRasterProperties_management<SPAN class="punctuation token">(</SPAN>raster<SPAN class="punctuation token">,</SPAN> <SPAN class="string token">"MAXIMUM"</SPAN><SPAN class="punctuation token">)</SPAN>
&nbsp;&nbsp;&nbsp; maxVal <SPAN class="operator token">=</SPAN> float<SPAN class="punctuation token">(</SPAN>val<SPAN class="punctuation token">.</SPAN>getOutput<SPAN class="punctuation token">(</SPAN><SPAN class="number token">0</SPAN><SPAN class="punctuation token">)</SPAN><SPAN class="punctuation token">)</SPAN>
&nbsp;&nbsp;&nbsp; maxVals<SPAN class="punctuation token">.</SPAN>append<SPAN class="punctuation token">(</SPAN>maxVal<SPAN class="punctuation token">)</SPAN><SPAN class="line-numbers-rows"><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN></SPAN>

All comments

XanderBakker

Hi anToni alcaraz

I'm not entirely sure what you want to achieve. When you mention you want to calculate the maximum, do you refer to a single value for the raster (like the metadata property of the raster) or the maximum for each pixel over the series of rasters or the maximum in a certain neighborhood (moving window) for each raster. And the mean would be using the values of each pixel for the series of maximum rasters?

You can create temp rasters using the in_memory workspace if they are not to big in case you do not want to create a physical result for the intermediate max rasters.

MartinSirkovsky1

Basically you would iterate over the rasters, find the maximum, add the result to array, and at the end you would calculate mean of the values in array.

You could use Get Raster Properties function to get the maximum of the raster.

JakeSkinner

You could try the following:

<SPAN class="keyword token">import</SPAN> arcpy
arcpy<SPAN class="punctuation token">.</SPAN>env<SPAN class="punctuation token">.</SPAN>workspace <SPAN class="operator token">=</SPAN> r<SPAN class="string token">"C:\TEMP\PYTHON\Data.gdb"</SPAN>

<SPAN class="comment token"># Create list</SPAN>
maxVals <SPAN class="operator token">=</SPAN> <SPAN class="punctuation token">[</SPAN><SPAN class="punctuation token">]</SPAN>

<SPAN class="comment token"># Iterate through each raster and append max val to list</SPAN>
<SPAN class="keyword token">for</SPAN> raster <SPAN class="keyword token">in</SPAN> arcpy<SPAN class="punctuation token">.</SPAN>ListRasters<SPAN class="punctuation token">(</SPAN><SPAN class="string token">"*"</SPAN><SPAN class="punctuation token">)</SPAN><SPAN class="punctuation token">:</SPAN>
&nbsp;&nbsp;&nbsp; val <SPAN class="operator token">=</SPAN> arcpy<SPAN class="punctuation token">.</SPAN>GetRasterProperties_management<SPAN class="punctuation token">(</SPAN>raster<SPAN class="punctuation token">,</SPAN> <SPAN class="string token">"MAXIMUM"</SPAN><SPAN class="punctuation token">)</SPAN>
&nbsp;&nbsp;&nbsp; maxVal <SPAN class="operator token">=</SPAN> int<SPAN class="punctuation token">(</SPAN>val<SPAN class="punctuation token">.</SPAN>getOutput<SPAN class="punctuation token">(</SPAN><SPAN class="number token">0</SPAN><SPAN class="punctuation token">)</SPAN><SPAN class="punctuation token">)</SPAN>
&nbsp;&nbsp;&nbsp; maxVals<SPAN class="punctuation token">.</SPAN>append<SPAN class="punctuation token">(</SPAN>maxVal<SPAN class="punctuation token">)</SPAN>

<SPAN class="comment token"># Sum list</SPAN>
sum <SPAN class="operator token">=</SPAN> <SPAN class="number token">0</SPAN>
<SPAN class="keyword token">for</SPAN> val <SPAN class="keyword token">in</SPAN> maxVals<SPAN class="punctuation token">:</SPAN>
&nbsp;&nbsp;&nbsp; sum <SPAN class="operator token">+=</SPAN> val

<SPAN class="comment token"># Find Mean</SPAN>
meanVal <SPAN class="operator token">=</SPAN> sum <SPAN class="operator token">/</SPAN> len<SPAN class="punctuation token">(</SPAN>maxVals<SPAN class="punctuation token">)</SPAN>
<SPAN class="keyword token">print</SPAN><SPAN class="punctuation token">(</SPAN>meanVal<SPAN class="punctuation token">)</SPAN><SPAN class="line-numbers-rows"><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN></SPAN>

DanPatterson_Retired

Can't let a numpy opportunity go by...

<SPAN class="keyword token">import</SPAN> arcpy
arcpy<SPAN class="punctuation token">.</SPAN>env<SPAN class="punctuation token">.</SPAN>workspace <SPAN class="operator token">=</SPAN> r<SPAN class="string token">"C:\Data\rasters"</SPAN>
<SPAN class="keyword token">import</SPAN> numpy <SPAN class="keyword token">as</SPAN> np
r_max <SPAN class="operator token">=</SPAN> np<SPAN class="punctuation token">.</SPAN>finfo<SPAN class="punctuation token">(</SPAN>np<SPAN class="punctuation token">.</SPAN>float<SPAN class="punctuation token">)</SPAN><SPAN class="punctuation token">.</SPAN>min
<SPAN class="keyword token">for</SPAN> raster <SPAN class="keyword token">in</SPAN> arcpy<SPAN class="punctuation token">.</SPAN>ListRasters<SPAN class="punctuation token">(</SPAN><SPAN class="string token">"*"</SPAN><SPAN class="punctuation token">)</SPAN><SPAN class="punctuation token">:</SPAN>
&nbsp;&nbsp;&nbsp; ras <SPAN class="operator token">=</SPAN> arcpy<SPAN class="punctuation token">.</SPAN>RasterToNumPyArray<SPAN class="punctuation token">(</SPAN>raster<SPAN class="punctuation token">)</SPAN>
&nbsp;&nbsp;&nbsp; r_max <SPAN class="operator token">=</SPAN> max<SPAN class="punctuation token">(</SPAN>r_max<SPAN class="punctuation token">,</SPAN> np<SPAN class="punctuation token">.</SPAN>nanmax<SPAN class="punctuation token">(</SPAN>ras<SPAN class="punctuation token">)</SPAN><SPAN class="punctuation token">)</SPAN>
<SPAN class="keyword token">print</SPAN><SPAN class="punctuation token">(</SPAN>r_max<SPAN class="punctuation token">)</SPAN>


<SPAN class="number token">255</SPAN><SPAN class="line-numbers-rows"><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN></SPAN>

anTonialcaraz

Thanks Xander and Martin for your suggestions.

Sorry Xander I may have not explained the issue correctly.

My ultimate goal here is to standardise a number of raster layers (696 of them). My first objective was to get the maximum value (a unique single value) for all those rasters. I've done that; but after examining the maximum value I realise that doing the mean of all the maximum values of all rasters would be a more accurate value (to later standardise them).

What Martin says is what I need to do but I'm not familiar with numpy and arrays and I'm not quite sure how to proceed. I also wonder if there is a function in scipy that would do that directly? Below is just the start of the script but as I said I don't know how to proceed.

<SPAN class="keyword token">import</SPAN> arcpy
<SPAN class="keyword token">import</SPAN> math
<SPAN class="keyword token">import</SPAN> os
<SPAN class="keyword token">import</SPAN> numpy
<SPAN class="keyword token">import</SPAN> scipy
<SPAN class="keyword token">from</SPAN> arcpy <SPAN class="keyword token">import</SPAN> env
arcpy<SPAN class="punctuation token">.</SPAN>env<SPAN class="punctuation token">.</SPAN>overwriteOutput<SPAN class="operator token">=</SPAN><SPAN class="token boolean">True</SPAN>
<SPAN class="keyword token">from</SPAN> arcpy<SPAN class="punctuation token">.</SPAN>sa <SPAN class="keyword token">import</SPAN><SPAN class="operator token">*</SPAN>
arcpy<SPAN class="punctuation token">.</SPAN>CheckOutExtension<SPAN class="punctuation token">(</SPAN><SPAN class="string token">"Spatial"</SPAN><SPAN class="punctuation token">)</SPAN>



env<SPAN class="punctuation token">.</SPAN>workspace <SPAN class="operator token">=</SPAN> r"D<SPAN class="punctuation token">:</SPAN>\000_TEST_MAGNITUDE_summary\MAX_VAL_calculation\Vel_rasters<SPAN class="punctuation token">.</SPAN>gdb



<SPAN class="comment token"># List rasters and get MAX values</SPAN>

rasterlist1 <SPAN class="operator token">=</SPAN> arcpy<SPAN class="punctuation token">.</SPAN>ListRasters<SPAN class="punctuation token">(</SPAN><SPAN class="punctuation token">)</SPAN>

Vel_max <SPAN class="operator token">=</SPAN> arcpy<SPAN class="punctuation token">.</SPAN>GetRasterProperties_management<SPAN class="punctuation token">(</SPAN>rasterlist1<SPAN class="punctuation token">,</SPAN> <SPAN class="string token">"MAXIMUM"</SPAN><SPAN class="punctuation token">)</SPAN>


<SPAN class="comment token"># Create array and add result from above</SPAN>

Vel_array <SPAN class="operator token">=</SPAN> arcpy<SPAN class="punctuation token">.</SPAN>RasterToNumPyArray<SPAN class="punctuation token">(</SPAN>rasterlist1<SPAN class="punctuation token">)</SPAN><SPAN class="line-numbers-rows"><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN></SPAN>

anTonialcaraz

Thank you Jake. That looks promising. Just trying now but I got an error in line 16. Is it '+' or '='?

JakeSkinner

It's actually +=. It's essentially a quicker way of writing:

sum <SPAN class="operator token">=</SPAN> sum <SPAN class="operator token">+</SPAN> val<SPAN class="line-numbers-rows"><SPAN>‍</SPAN></SPAN>

XanderBakker

Could you try this and share what it prints (it's basically the same code as Jake shared earlier)?

<SPAN class="keyword token">import</SPAN> arcpy
<SPAN class="keyword token">import</SPAN> os

<SPAN class="comment token"># settings</SPAN>
ws <SPAN class="operator token">=</SPAN> r<SPAN class="string token">"D:\000_TEST_MAGNITUDE_summary\MAX_VAL_calculation\Vel_rasters.gdb"</SPAN>
arcpy<SPAN class="punctuation token">.</SPAN>env<SPAN class="punctuation token">.</SPAN>workspace <SPAN class="operator token">=</SPAN> ws

<SPAN class="comment token"># List rasters</SPAN>
raster_names <SPAN class="operator token">=</SPAN> arcpy<SPAN class="punctuation token">.</SPAN>ListRasters<SPAN class="punctuation token">(</SPAN><SPAN class="punctuation token">)</SPAN>

<SPAN class="comment token"># loop through rasters</SPAN>
max_values <SPAN class="operator token">=</SPAN> <SPAN class="punctuation token">[</SPAN><SPAN class="punctuation token">]</SPAN>
<SPAN class="keyword token">for</SPAN> raster_name <SPAN class="keyword token">in</SPAN> raster_names<SPAN class="punctuation token">:</SPAN>
&nbsp;&nbsp;&nbsp; raster_path <SPAN class="operator token">=</SPAN> os<SPAN class="punctuation token">.</SPAN>path<SPAN class="punctuation token">.</SPAN>join<SPAN class="punctuation token">(</SPAN>ws<SPAN class="punctuation token">,</SPAN> raster_name<SPAN class="punctuation token">)</SPAN>
&nbsp;&nbsp;&nbsp; vel_max <SPAN class="operator token">=</SPAN> arcpy<SPAN class="punctuation token">.</SPAN>GetRasterProperties_management<SPAN class="punctuation token">(</SPAN>raster_path<SPAN class="punctuation token">,</SPAN> <SPAN class="string token">"MAXIMUM"</SPAN><SPAN class="punctuation token">)</SPAN>
&nbsp;&nbsp;&nbsp; <SPAN class="keyword token">print</SPAN><SPAN class="punctuation token">(</SPAN><SPAN class="string token">" - {}: {}"</SPAN><SPAN class="punctuation token">.</SPAN>format<SPAN class="punctuation token">(</SPAN>raster_name<SPAN class="punctuation token">,</SPAN> vel_max<SPAN class="punctuation token">)</SPAN><SPAN class="punctuation token">)</SPAN>
&nbsp;&nbsp;&nbsp; max_values<SPAN class="punctuation token">.</SPAN>append<SPAN class="punctuation token">(</SPAN>float<SPAN class="punctuation token">(</SPAN>vel_max<SPAN class="punctuation token">.</SPAN>getOutput<SPAN class="punctuation token">(</SPAN><SPAN class="number token">0</SPAN><SPAN class="punctuation token">)</SPAN><SPAN class="punctuation token">)</SPAN><SPAN class="punctuation token">)</SPAN>

mean_val <SPAN class="operator token">=</SPAN> float<SPAN class="punctuation token">(</SPAN>sum<SPAN class="punctuation token">(</SPAN><SPAN class="punctuation token">[</SPAN>max_values<SPAN class="punctuation token">]</SPAN><SPAN class="punctuation token">)</SPAN><SPAN class="punctuation token">)</SPAN> <SPAN class="operator token">/</SPAN> len<SPAN class="punctuation token">(</SPAN>max_values<SPAN class="punctuation token">)</SPAN>
<SPAN class="keyword token">print</SPAN><SPAN class="punctuation token">(</SPAN><SPAN class="string token">"Mean value: {}"</SPAN><SPAN class="punctuation token">.</SPAN>format<SPAN class="punctuation token">(</SPAN>mean_val<SPAN class="punctuation token">)</SPAN><SPAN class="punctuation token">)</SPAN><SPAN class="line-numbers-rows"><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN></SPAN>

DanPatterson_Retired

Mean of max's? no problem

<SPAN class="keyword token">import</SPAN> arcpy
arcpy<SPAN class="punctuation token">.</SPAN>env<SPAN class="punctuation token">.</SPAN>workspace <SPAN class="operator token">=</SPAN> r<SPAN class="string token">"C:\Data\rasters"</SPAN>
<SPAN class="keyword token">import</SPAN> numpy <SPAN class="keyword token">as</SPAN> np
r_max_mean <SPAN class="operator token">=</SPAN> <SPAN class="punctuation token">[</SPAN><SPAN class="punctuation token">]</SPAN>
<SPAN class="keyword token">for</SPAN> raster <SPAN class="keyword token">in</SPAN> arcpy<SPAN class="punctuation token">.</SPAN>ListRasters<SPAN class="punctuation token">(</SPAN><SPAN class="string token">"*"</SPAN><SPAN class="punctuation token">)</SPAN><SPAN class="punctuation token">:</SPAN>
&nbsp;&nbsp;&nbsp; ras <SPAN class="operator token">=</SPAN> arcpy<SPAN class="punctuation token">.</SPAN>RasterToNumPyArray<SPAN class="punctuation token">(</SPAN>raster<SPAN class="punctuation token">)</SPAN>
&nbsp;&nbsp;&nbsp; r_max_mean<SPAN class="punctuation token">.</SPAN>append<SPAN class="punctuation token">(</SPAN> np<SPAN class="punctuation token">.</SPAN>nanmax<SPAN class="punctuation token">(</SPAN>ras<SPAN class="punctuation token">)</SPAN> <SPAN class="punctuation token">)</SPAN>
<SPAN class="keyword token">print</SPAN><SPAN class="punctuation token">(</SPAN><SPAN class="string token">"max_list\n{}\nmean_max {}"</SPAN><SPAN class="punctuation token">.</SPAN>format<SPAN class="punctuation token">(</SPAN>r_max_mean<SPAN class="punctuation token">,</SPAN> np<SPAN class="punctuation token">.</SPAN>nanmean<SPAN class="punctuation token">(</SPAN>r_max_mean<SPAN class="punctuation token">)</SPAN><SPAN class="punctuation token">)</SPAN><SPAN class="punctuation token">)</SPAN>


max_list
<SPAN class="punctuation token">[</SPAN><SPAN class="number token">255</SPAN><SPAN class="punctuation token">,</SPAN> <SPAN class="number token">99</SPAN><SPAN class="punctuation token">,</SPAN> <SPAN class="number token">99</SPAN><SPAN class="punctuation token">,</SPAN> <SPAN class="number token">99</SPAN><SPAN class="punctuation token">,</SPAN> <SPAN class="number token">99</SPAN><SPAN class="punctuation token">,</SPAN> <SPAN class="number token">99</SPAN><SPAN class="punctuation token">,</SPAN> <SPAN class="number token">99</SPAN><SPAN class="punctuation token">,</SPAN> <SPAN class="number token">99</SPAN><SPAN class="punctuation token">,</SPAN> <SPAN class="number token">99</SPAN><SPAN class="punctuation token">,</SPAN> <SPAN class="number token">99</SPAN><SPAN class="punctuation token">,</SPAN> <SPAN class="number token">99</SPAN><SPAN class="punctuation token">,</SPAN> <SPAN class="number token">99</SPAN><SPAN class="punctuation token">,</SPAN> <SPAN class="number token">99</SPAN><SPAN class="punctuation token">,</SPAN> <SPAN class="number token">99</SPAN><SPAN class="punctuation token">,</SPAN> <SPAN class="number token">99</SPAN><SPAN class="punctuation token">,</SPAN> <SPAN class="number token">99</SPAN><SPAN class="punctuation token">,</SPAN> <SPAN class="number token">99</SPAN><SPAN class="punctuation token">,</SPAN> <SPAN class="number token">99</SPAN><SPAN class="punctuation token">,</SPAN>
 <SPAN class="number token">99</SPAN><SPAN class="punctuation token">,</SPAN> <SPAN class="number token">99</SPAN><SPAN class="punctuation token">,</SPAN> <SPAN class="number token">99</SPAN><SPAN class="punctuation token">,</SPAN> <SPAN class="number token">99</SPAN><SPAN class="punctuation token">,</SPAN> <SPAN class="number token">99</SPAN><SPAN class="punctuation token">,</SPAN> <SPAN class="number token">99</SPAN><SPAN class="punctuation token">,</SPAN> <SPAN class="number token">99</SPAN><SPAN class="punctuation token">,</SPAN> <SPAN class="number token">99</SPAN><SPAN class="punctuation token">,</SPAN> <SPAN class="number token">99</SPAN><SPAN class="punctuation token">,</SPAN> <SPAN class="number token">99</SPAN><SPAN class="punctuation token">,</SPAN> <SPAN class="number token">99</SPAN><SPAN class="punctuation token">,</SPAN> <SPAN class="number token">99</SPAN><SPAN class="punctuation token">,</SPAN> <SPAN class="number token">99</SPAN><SPAN class="punctuation token">,</SPAN> <SPAN class="number token">99</SPAN><SPAN class="punctuation token">]</SPAN>

mean_max <SPAN class="number token">103.875</SPAN><SPAN class="line-numbers-rows"><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN></SPAN>

anTonialcaraz

I see. At first run I got an error in line 13 below. As the rasters are floating point I removed 'int' from the line (as below).

env<SPAN class="punctuation token">.</SPAN>workspace <SPAN class="operator token">=</SPAN> r<SPAN class="string token">"D:\000_TEST_MAGNITUDE_summary\MAX_VAL_calculation\Vel_rasters.gdb"</SPAN>



<SPAN class="comment token"># Create list</SPAN>

maxVals <SPAN class="operator token">=</SPAN> <SPAN class="punctuation token">[</SPAN><SPAN class="punctuation token">]</SPAN>

<SPAN class="comment token"># Iterate through each raster and append max val to list</SPAN>

<SPAN class="keyword token">for</SPAN> raster <SPAN class="keyword token">in</SPAN> arcpy<SPAN class="punctuation token">.</SPAN>ListRasters<SPAN class="punctuation token">(</SPAN><SPAN class="string token">"*"</SPAN><SPAN class="punctuation token">)</SPAN><SPAN class="punctuation token">:</SPAN>
&nbsp;&nbsp;&nbsp; val <SPAN class="operator token">=</SPAN> arcpy<SPAN class="punctuation token">.</SPAN>GetRasterProperties_management<SPAN class="punctuation token">(</SPAN>raster<SPAN class="punctuation token">,</SPAN> <SPAN class="string token">"MAXIMUM"</SPAN><SPAN class="punctuation token">)</SPAN>
&nbsp;&nbsp;&nbsp; maxVal <SPAN class="operator token">=</SPAN> val<SPAN class="punctuation token">.</SPAN>getOutput<SPAN class="punctuation token">(</SPAN><SPAN class="number token">0</SPAN><SPAN class="punctuation token">)</SPAN>
&nbsp;&nbsp;&nbsp; maxVals<SPAN class="punctuation token">.</SPAN>append<SPAN class="punctuation token">(</SPAN>maxVal<SPAN class="punctuation token">)</SPAN>

<SPAN class="comment token"># Sum list</SPAN>

sum <SPAN class="operator token">=</SPAN> <SPAN class="number token">0</SPAN>
<SPAN class="keyword token">for</SPAN> val <SPAN class="keyword token">in</SPAN> maxVals<SPAN class="punctuation token">:</SPAN>
&nbsp;&nbsp;&nbsp; sum <SPAN class="operator token">=</SPAN> sum <SPAN class="operator token">+</SPAN> val

<SPAN class="comment token"># Find Mean</SPAN>

meanVal <SPAN class="operator token">=</SPAN> sum <SPAN class="operator token">/</SPAN> len<SPAN class="punctuation token">(</SPAN>maxVals<SPAN class="punctuation token">)</SPAN>


<SPAN class="keyword token">print</SPAN><SPAN class="punctuation token">(</SPAN>meanVal<SPAN class="punctuation token">)</SPAN>‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍<SPAN class="line-numbers-rows"><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN></SPAN>

But now I'm getting this error which I guess it may be related.

DanPatterson_Retired

scipy relies on numpy, see my second example, it returns a list of the maximums, then simply takes the mean of that

anTonialcaraz

Thank you Xander. It seems there is a minor issue in line 19.

I also tried with float(len(max_values) but getting the same error. Other than that, I'm sure is perfect. Very nice the way it prints the max values for each raster. Both Jake's and Dan's are working and giving the same result. Many thanks.

anTonialcaraz

That works perfect Dan. All the slightly different approaches above give same result. Really appreciate your help once again.

Nimm an der Diskussion teil

Angemeldete Mitglieder können Beiträge verfassen, Updates folgen und mehr. Neu hier? Registriere ein kostenloses Konto.

Schnellzugriffe

Community Help
Find useful guides, FAQs, and documents to help you navigate and make the most of Esri Community.