Hi, I'm trying to run Mean annual cell statistics for groups of monthly rasters. The output I'm getting is the same as the last raster in the list ("camp_t_srf_sep" in this example) so i assume the script is not actually iterating through all rasters. Any help will be greatly appreciate it. # Example of list of rasters: # camp_t_srf_apr # camp_t_srf_aug # camp_t_srf_dec # camp_t_srf_feb # camp_t_srf_jan # camp_t_srf_jul # camp_t_srf_jun # camp_t_srf_mar # camp_t_srf_may # camp_t_srf_nov # camp_t_srf_oct # camp_t_srf_sep import arcpy from arcpy import env from arcpy . sa import * arcpy . CheckOutExtension ( "Spatial" ) env . overwriteOutput = True env . extent = arcpy . Extent ( - 180.0 , - 90.0 , 180.0 , 90.0 ) env . workspace = arcpy . GetParameterAsText ( 0 ) StageAge = arcpy . GetParameterAsText ( 1 ) rasterlist = arcpy . ListRasters ( StageAge + "_*" ) for raster in rasterlist : output_mean = CellStatistics ( raster , "MEAN" , "DATA" ) output_mean . save ( env . workspace + "\\" + raster [ : 11 ] + "ann" )

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; } } }

CELL STATISTICS with LIST RASTERS

anTonialcaraz

Hi,

I'm trying to run Mean annual cell statistics for groups of monthly rasters.

The output I'm getting is the same as the last raster in the list ("camp_t_srf_sep" in this example) so i assume the script is not actually iterating through all rasters.

Any help will be greatly appreciate it.

<SPAN class="comment token"># Example of list of rasters:</SPAN>
 <SPAN class="comment token"># camp_t_srf_apr</SPAN>
 <SPAN class="comment token"># camp_t_srf_aug</SPAN>
 <SPAN class="comment token"># camp_t_srf_dec</SPAN>
 <SPAN class="comment token"># camp_t_srf_feb</SPAN>
 <SPAN class="comment token"># camp_t_srf_jan</SPAN>
 <SPAN class="comment token"># camp_t_srf_jul</SPAN>
 <SPAN class="comment token"># camp_t_srf_jun</SPAN>
 <SPAN class="comment token"># camp_t_srf_mar</SPAN>
 <SPAN class="comment token"># camp_t_srf_may</SPAN>
 <SPAN class="comment token"># camp_t_srf_nov</SPAN>
 <SPAN class="comment token"># camp_t_srf_oct</SPAN>
 <SPAN class="comment token"># camp_t_srf_sep</SPAN>

<SPAN class="keyword token">import</SPAN> arcpy
<SPAN class="keyword token">from</SPAN> arcpy <SPAN class="keyword token">import</SPAN> env
<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>overwriteOutput <SPAN class="operator token">=</SPAN> <SPAN class="token boolean">True</SPAN>
env<SPAN class="punctuation token">.</SPAN>extent <SPAN class="operator token">=</SPAN> arcpy<SPAN class="punctuation token">.</SPAN>Extent<SPAN class="punctuation token">(</SPAN><SPAN class="operator token">-</SPAN><SPAN class="number token">180.0</SPAN><SPAN class="punctuation token">,</SPAN> <SPAN class="operator token">-</SPAN><SPAN class="number token">90.0</SPAN><SPAN class="punctuation token">,</SPAN> <SPAN class="number token">180.0</SPAN><SPAN class="punctuation token">,</SPAN> <SPAN class="number token">90.0</SPAN><SPAN class="punctuation token">)</SPAN>

env<SPAN class="punctuation token">.</SPAN>workspace <SPAN class="operator token">=</SPAN> arcpy<SPAN class="punctuation token">.</SPAN>GetParameterAsText<SPAN class="punctuation token">(</SPAN><SPAN class="number token">0</SPAN><SPAN class="punctuation token">)</SPAN>
StageAge <SPAN class="operator token">=</SPAN> arcpy<SPAN class="punctuation token">.</SPAN>GetParameterAsText<SPAN class="punctuation token">(</SPAN><SPAN class="number token">1</SPAN><SPAN class="punctuation token">)</SPAN>


rasterlist <SPAN class="operator token">=</SPAN> arcpy<SPAN class="punctuation token">.</SPAN>ListRasters<SPAN class="punctuation token">(</SPAN>StageAge <SPAN class="operator token">+</SPAN> <SPAN class="string token">"_*"</SPAN><SPAN class="punctuation token">)</SPAN>
<SPAN class="keyword token">for</SPAN> raster <SPAN class="keyword token">in</SPAN> rasterlist<SPAN class="punctuation token">:</SPAN>
&nbsp;&nbsp;&nbsp; output_mean <SPAN class="operator token">=</SPAN> CellStatistics<SPAN class="punctuation token">(</SPAN>raster<SPAN class="punctuation token">,</SPAN> <SPAN class="string token">"MEAN"</SPAN><SPAN class="punctuation token">,</SPAN> <SPAN class="string token">"DATA"</SPAN><SPAN class="punctuation token">)</SPAN>
&nbsp;&nbsp;&nbsp; output_mean<SPAN class="punctuation token">.</SPAN>save<SPAN class="punctuation token">(</SPAN>env<SPAN class="punctuation token">.</SPAN>workspace <SPAN class="operator token">+</SPAN> <SPAN class="string token">"\\"</SPAN> <SPAN class="operator token">+</SPAN> raster<SPAN class="punctuation token">[</SPAN><SPAN class="punctuation token">:</SPAN><SPAN class="number token">11</SPAN><SPAN class="punctuation token">]</SPAN> <SPAN class="operator token">+</SPAN> <SPAN class="string token">"ann"</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>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN></SPAN>

Find more posts tagged with

Accepted answers

DanPatterson_Retired

You can use your for loop to build the list... just don't have cellstatistics in the loop

All comments

anTonialcaraz

I know it works like as below. Just thought there would be a more elegant and economic way to put it using ListRasters

rasterlist5 <SPAN class="operator token">=</SPAN> <SPAN class="punctuation token">[</SPAN>StageAge <SPAN class="operator token">+</SPAN> <SPAN class="string token">"_t_srf_apr"</SPAN><SPAN class="punctuation token">,</SPAN> StageAge <SPAN class="operator token">+</SPAN> <SPAN class="string token">"_t_srf_aug"</SPAN><SPAN class="punctuation token">,</SPAN> StageAge <SPAN class="operator token">+</SPAN> <SPAN class="string token">"_t_srf_dec"</SPAN><SPAN class="punctuation token">,</SPAN> StageAge <SPAN class="operator token">+</SPAN> <SPAN class="string token">"_t_srf_feb"</SPAN><SPAN class="punctuation token">,</SPAN> StageAge <SPAN class="operator token">+</SPAN> <SPAN class="string token">"_t_srf_jan"</SPAN><SPAN class="punctuation token">,</SPAN> StageAge <SPAN class="operator token">+</SPAN> <SPAN class="string token">"_t_srf_jul"</SPAN><SPAN class="punctuation token">,</SPAN> StageAge <SPAN class="operator token">+</SPAN> <SPAN class="string token">"_t_srf_jun"</SPAN><SPAN class="punctuation token">,</SPAN> \
 StageAge <SPAN class="operator token">+</SPAN> <SPAN class="string token">"_t_srf_mar"</SPAN><SPAN class="punctuation token">,</SPAN> StageAge <SPAN class="operator token">+</SPAN> <SPAN class="string token">"_t_srf_may"</SPAN><SPAN class="punctuation token">,</SPAN> StageAge <SPAN class="operator token">+</SPAN> <SPAN class="string token">"_t_srf_nov"</SPAN><SPAN class="punctuation token">,</SPAN> StageAge <SPAN class="operator token">+</SPAN> <SPAN class="string token">"_t_srf_oct"</SPAN><SPAN class="punctuation token">,</SPAN> StageAge <SPAN class="operator token">+</SPAN> <SPAN class="string token">"_t_srf_sep"</SPAN><SPAN class="punctuation token">]</SPAN>

output_mean <SPAN class="operator token">=</SPAN> CellStatistics<SPAN class="punctuation token">(</SPAN>rasterlist5<SPAN class="punctuation token">,</SPAN> <SPAN class="string token">"MEAN"</SPAN><SPAN class="punctuation token">,</SPAN> <SPAN class="string token">"DATA"</SPAN><SPAN class="punctuation token">)</SPAN>

output_mean<SPAN class="punctuation token">.</SPAN>save<SPAN class="punctuation token">(</SPAN>OUTWorkspace <SPAN class="operator token">+</SPAN> <SPAN class="string token">"\\"</SPAN> <SPAN class="operator token">+</SPAN> raster<SPAN class="punctuation token">[</SPAN><SPAN class="punctuation token">:</SPAN><SPAN class="number token">11</SPAN><SPAN class="punctuation token">]</SPAN> <SPAN class="operator token">+</SPAN> <SPAN class="string token">"ann"</SPAN><SPAN class="punctuation token">)</SPAN>‍‍‍‍<SPAN class="line-numbers-rows"><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN></SPAN>

DanPatterson_Retired

because you are cycling through each raster one at a time.

for raster in rasterlist:

give it the whole list... examine the syntax in the code sample in the link I provided

anTonialcaraz

Thank you Dan.

yes, the print statement gives me the list of rasters I want to calculate the mean from.

I thought I was giving the script the list of rasters with the variable rasterlist. I don't quite understand why it only processes the last raster in the list.

DanPatterson_Retired

Throw a print statement in... just to make sure you are getting the correct names, it will silently do nothing otherwise

And also.. Cell statistics expects a list/stack of rasters, not one at a time

http://desktop.arcgis.com/en/arcmap/latest/tools/spatial-analyst-toolbox/cell-statistics.htm

anTonialcaraz

just to add that I'm using 10.2.1 at the moment.

many thanks

Join the Conversation

Signed in members can post, follow updates, and more. New here? Register a free account.

Quick Links

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