Hi, I am trying to calculate the average using multiple rasters using years from 1981 to 2016. For example, average for week 1 for the month of 9 from years 1981 to 2016 is as follows: Normal_Week_1_1981_2016_9.tif =average [Week_1_Sum1981_9.tif . . Week_1_Sum2016_9.tif] My following code outputs all the tif files at the "out" path but at the end it gives me the error below. What I am doing wrong here.? Just noticed the output raster values are also wrong. Any help is appreciated. Thanks in advance. arcpy . CheckOutExtension ( "Spatial" ) arcpy . OverwriteOutput = True years = [ 1981 , 1982 , 1983 , 1984 , 1985 , 1986 , 1987 , 1988 , 1989 , 1990 , 1991 , 1992 , 1993 , 1994 , 1995 , 1996 , 1997 , 1998 , 1999 , 2000 , 2001 , 2002 , 2003 , 2004 , 2005 , 2006 , 2007 , 2008 , 2009 , 2010 , 2011 , 2012 , 2013 , 2014 , 2015 , 2016 ] out = r 'H:\PRISM_800m_weekly_sum\Normal_out' path = r 'H:\PRISM_800m_weekly_sum_Normal' lstFiles = [ ] for root , dirs , files in os . walk ( path ) : arcpy . env . workspace = root files = arcpy . ListRasters ( "*" , "TIF" ) for file in files : ( filename , extension ) = os . path . splitext ( file ) if ( extension == ".tif" ) : lstFiles . append ( os . path . join ( root , file ) ) print ( lstFiles ) new_list1 = [ ] new_list2 = [ ] new_list3 = [ ] new_list4 = [ ] new_list5 = [ ] months = [ 9 , 10 , 11 ] for month in months : for file in lstFiles : if fnmatch . fnmatch ( file , "H:\\PRISM_800m_weekly_sum_Normal" + "*" + "\\" + "Week_1" + "_*_" + str ( month ) + ".tif" ) : new_list1 . append ( file ) new_list1 . sort ( ) print ( new_list1 ) env . workspace = out finras1 = out + "\\" + "Normal_Week_1" + "_1981_2016" + "_" + str ( month ) + ".tif" calc1 = arcpy . sa . CellStatistics ( new_list1 , statistics_type = "MEAN" ) arcpy . CopyRaster_management ( calc1 , finras1 , "" , "" , "" , "" , "" , "32_BIT_FLOAT" ) for file in lstFiles : if fnmatch . fnmatch ( file , "H:\\PRISM_800m_weekly_sum_Normal" + "*" + "\\" + "Week_2" + "_*_" + str ( month ) + ".tif" ) : new_list2 . append ( file ) new_list2 . sort ( ) print ( new_list2 ) env . workspace = out finras2 = out + "\\" + "Normal_Week_2" + "_1981_2016" + "_" + str ( month ) + ".tif" calc2 = arcpy . sa . CellStatistics ( new_list2 , statistics_type = "MEAN" ) arcpy . CopyRaster_management ( calc2 , finras2 , "" , "" , "" , "" , "" , "32_BIT_FLOAT" ) for file in lstFiles : if fnmatch . fnmatch ( file , "H:\\PRISM_800m_weekly_sum_Normal" + "*" + "\\" + "Week_3" + "_*_" + str ( month ) + ".tif" ) : new_list3 . append ( file ) new_list3 . sort ( ) print ( new_list3 ) env . workspace = out finras3 = out + "\\" + "Normal_Week_3" + "_1981_2016" + "_" + str ( month ) + ".tif" calc3 = arcpy . sa . CellStatistics ( new_list3 , statistics_type = "MEAN" ) arcpy . CopyRaster_management ( calc3 , finras3 , "" , "" , "" , "" , "" , "32_BIT_FLOAT" ) for file in lstFiles : if fnmatch . fnmatch ( file , "H:\\PRISM_800m_weekly_sum_Normal" + "*" + "\\" + "Week_4" + "_*_" + str ( month ) + ".tif" ) : new_list4 . append ( file ) new_list4 . sort ( ) print ( new_list4 ) env . workspace = out finras4 = out + "\\" + "Normal_Week_4" + "_1981_2016" + "_" + str ( month ) + ".tif" calc4 = arcpy . sa . CellStatistics ( new_list4 , statistics_type = "MEAN" ) arcpy . CopyRaster_management ( calc4 , finras4 , "" , "" , "" , "" , "" , "32_BIT_FLOAT" ) for file in lstFiles : if fnmatch . fnmatch ( file , "H:\\PRISM_800m_weekly_sum_Normal" + "*" + "\\" + "Week_5" + "_*_" + str ( month ) + ".tif" ) : new_list5 . append ( file ) new_list5 . sort ( ) print ( new_list5 ) env . workspace = out finras5 = out + "\\" + "Normal_Week_5" + "_1981_2016" + "_" + str ( month ) + ".tif" calc5 = arcpy . sa . CellStatistics ( new_list5 , statistics_type = "MEAN" ) arcpy . CopyRaster_management ( calc5 , finras5 , "" , "" , "" , "" , "" , "32_BIT_FLOAT" )

Glad that the output is showing the grouping that you want with your modification. I have added back in some of your earlier code, but I may have missed something. import arcpy # make sure this directory exists out = r 'H:\PRISM_800m_weekly_sum\Normal_out' for week in range ( 1 , 5 + 1 ) : # for weeks 1 to 5 for month in range ( 9 , 11 + 1 ) : # for months 9 to 11 rasterlist = [ ] # initialize a list for the rasters for year in range ( 1981 , 2016 + 1 ) : # for years 1981 to 2016 rasterfile = "H:\\PRISM_800m_weekly_sum_Normal\\{}\\Week_{}_Sum{}_{}.tif" . format ( year , week , year , month ) # probably should add some error checking to see that file exists before adding it to rasterlist if arcpy . Exists ( rasterfile ) : rasterlist . append ( rasterfile ) # if found, append file to the raster list else : print "Not found: {}" . format ( rasterfile ) # same indentation as "for year in..." # rasterlist.sort() # you probably don't need to sort the list finras1 = out + "\\Normal_Week_{}_1981_2016_{}.tif" . format ( week , month ) print "Processing file: {}" . format ( finras1 ) calc1 = arcpy . sa . CellStatistics ( rasterlist , statistics_type = "MEAN" , ignore_nodata = "DATA" ) arcpy . CopyRaster_management ( calc1 , finras1 , "" , "" , "" , "" , "" , "32_BIT_FLOAT" ) print "Done." Also, keep in mind what Dan Patterson said in his comments in another of your questions regarding working with rasters. Remember that indentation is important, and it is often good to have some print statements in your code to help you know where things are going right or wrong. Hope this helps.

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

Error in multiple raster processing

DuminduJayasekera

Hi,

I am trying to calculate the average using multiple rasters using years from 1981 to 2016. For example, average for week 1 for the month of 9 from years 1981 to 2016 is as follows:

Normal_Week_1_1981_2016_9.tif =average [Week_1_Sum1981_9.tif

Week_1_Sum2016_9.tif]

My following code outputs all the tif files at the "out" path but at the end it gives me the error below. What I am doing wrong here.? Just noticed the output raster values are also wrong.

Any help is appreciated.

Thanks in advance.

arcpy<SPAN class="punctuation token">.</SPAN>CheckOutExtension<SPAN class="punctuation token">(</SPAN><SPAN class="string token">"Spatial"</SPAN><SPAN class="punctuation token">)</SPAN>
arcpy<SPAN class="punctuation token">.</SPAN>OverwriteOutput <SPAN class="operator token">=</SPAN> <SPAN class="token boolean">True</SPAN>

years <SPAN class="operator token">=</SPAN> <SPAN class="punctuation token">[</SPAN><SPAN class="number token">1981</SPAN><SPAN class="punctuation token">,</SPAN><SPAN class="number token">1982</SPAN><SPAN class="punctuation token">,</SPAN><SPAN class="number token">1983</SPAN><SPAN class="punctuation token">,</SPAN><SPAN class="number token">1984</SPAN><SPAN class="punctuation token">,</SPAN><SPAN class="number token">1985</SPAN><SPAN class="punctuation token">,</SPAN><SPAN class="number token">1986</SPAN><SPAN class="punctuation token">,</SPAN><SPAN class="number token">1987</SPAN><SPAN class="punctuation token">,</SPAN><SPAN class="number token">1988</SPAN><SPAN class="punctuation token">,</SPAN><SPAN class="number token">1989</SPAN><SPAN class="punctuation token">,</SPAN><SPAN class="number token">1990</SPAN><SPAN class="punctuation token">,</SPAN><SPAN class="number token">1991</SPAN><SPAN class="punctuation token">,</SPAN><SPAN class="number token">1992</SPAN><SPAN class="punctuation token">,</SPAN><SPAN class="number token">1993</SPAN><SPAN class="punctuation token">,</SPAN><SPAN class="number token">1994</SPAN><SPAN class="punctuation token">,</SPAN><SPAN class="number token">1995</SPAN><SPAN class="punctuation token">,</SPAN><SPAN class="number token">1996</SPAN><SPAN class="punctuation token">,</SPAN><SPAN class="number token">1997</SPAN><SPAN class="punctuation token">,</SPAN><SPAN class="number token">1998</SPAN><SPAN class="punctuation token">,</SPAN><SPAN class="number token">1999</SPAN><SPAN class="punctuation token">,</SPAN><SPAN class="number token">2000</SPAN><SPAN class="punctuation token">,</SPAN><SPAN class="number token">2001</SPAN><SPAN class="punctuation token">,</SPAN><SPAN class="number token">2002</SPAN><SPAN class="punctuation token">,</SPAN><SPAN class="number token">2003</SPAN><SPAN class="punctuation token">,</SPAN><SPAN class="number token">2004</SPAN><SPAN class="punctuation token">,</SPAN><SPAN class="number token">2005</SPAN><SPAN class="punctuation token">,</SPAN><SPAN class="number token">2006</SPAN><SPAN class="punctuation token">,</SPAN><SPAN class="number token">2007</SPAN><SPAN class="punctuation token">,</SPAN><SPAN class="number token">2008</SPAN><SPAN class="punctuation token">,</SPAN><SPAN class="number token">2009</SPAN><SPAN class="punctuation token">,</SPAN><SPAN class="number token">2010</SPAN><SPAN class="punctuation token">,</SPAN><SPAN class="number token">2011</SPAN><SPAN class="punctuation token">,</SPAN><SPAN class="number token">2012</SPAN><SPAN class="punctuation token">,</SPAN><SPAN class="number token">2013</SPAN><SPAN class="punctuation token">,</SPAN><SPAN class="number token">2014</SPAN><SPAN class="punctuation token">,</SPAN><SPAN class="number token">2015</SPAN><SPAN class="punctuation token">,</SPAN><SPAN class="number token">2016</SPAN><SPAN class="punctuation token">]</SPAN>

out <SPAN class="operator token">=</SPAN> r<SPAN class="string token">'H:\PRISM_800m_weekly_sum\Normal_out'</SPAN>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 

path <SPAN class="operator token">=</SPAN> r<SPAN class="string token">'H:\PRISM_800m_weekly_sum_Normal'</SPAN>
lstFiles <SPAN class="operator token">=</SPAN> <SPAN class="punctuation token">[</SPAN><SPAN class="punctuation token">]</SPAN>&nbsp; 
<SPAN class="keyword token">for</SPAN> root<SPAN class="punctuation token">,</SPAN> dirs<SPAN class="punctuation token">,</SPAN> files <SPAN class="keyword token">in</SPAN> os<SPAN class="punctuation token">.</SPAN>walk<SPAN class="punctuation token">(</SPAN>path<SPAN class="punctuation token">)</SPAN><SPAN class="punctuation token">:</SPAN>&nbsp; 
&nbsp;&nbsp;&nbsp; arcpy<SPAN class="punctuation token">.</SPAN>env<SPAN class="punctuation token">.</SPAN>workspace <SPAN class="operator token">=</SPAN> root&nbsp; 
&nbsp;&nbsp;&nbsp; files <SPAN class="operator token">=</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="string token">"TIF"</SPAN><SPAN class="punctuation token">)</SPAN>&nbsp; 
&nbsp;&nbsp;&nbsp; <SPAN class="keyword token">for</SPAN> file <SPAN class="keyword token">in</SPAN> files<SPAN class="punctuation token">:</SPAN>&nbsp; 
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; <SPAN class="punctuation token">(</SPAN>filename<SPAN class="punctuation token">,</SPAN> extension<SPAN class="punctuation token">)</SPAN> <SPAN class="operator token">=</SPAN> os<SPAN class="punctuation token">.</SPAN>path<SPAN class="punctuation token">.</SPAN>splitext<SPAN class="punctuation token">(</SPAN>file<SPAN class="punctuation token">)</SPAN>&nbsp; 
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; <SPAN class="keyword token">if</SPAN><SPAN class="punctuation token">(</SPAN>extension <SPAN class="operator token">==</SPAN> <SPAN class="string token">".tif"</SPAN><SPAN class="punctuation token">)</SPAN><SPAN class="punctuation token">:</SPAN>&nbsp; 
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; lstFiles<SPAN class="punctuation token">.</SPAN>append<SPAN class="punctuation token">(</SPAN>os<SPAN class="punctuation token">.</SPAN>path<SPAN class="punctuation token">.</SPAN>join<SPAN class="punctuation token">(</SPAN>root<SPAN class="punctuation token">,</SPAN> file<SPAN class="punctuation token">)</SPAN><SPAN class="punctuation token">)</SPAN>&nbsp; 
<SPAN class="keyword token">print</SPAN> <SPAN class="punctuation token">(</SPAN>lstFiles<SPAN class="punctuation token">)</SPAN>&nbsp; 

new_list1<SPAN class="operator token">=</SPAN><SPAN class="punctuation token">[</SPAN><SPAN class="punctuation token">]</SPAN>
new_list2<SPAN class="operator token">=</SPAN><SPAN class="punctuation token">[</SPAN><SPAN class="punctuation token">]</SPAN>
new_list3<SPAN class="operator token">=</SPAN><SPAN class="punctuation token">[</SPAN><SPAN class="punctuation token">]</SPAN>
new_list4<SPAN class="operator token">=</SPAN><SPAN class="punctuation token">[</SPAN><SPAN class="punctuation token">]</SPAN>
new_list5<SPAN class="operator token">=</SPAN><SPAN class="punctuation token">[</SPAN><SPAN class="punctuation token">]</SPAN>

months<SPAN class="operator token">=</SPAN><SPAN class="punctuation token">[</SPAN><SPAN class="number token">9</SPAN><SPAN class="punctuation token">,</SPAN><SPAN class="number token">10</SPAN><SPAN class="punctuation token">,</SPAN><SPAN class="number token">11</SPAN><SPAN class="punctuation token">]</SPAN>
<SPAN class="keyword token">for</SPAN> month <SPAN class="keyword token">in</SPAN> months<SPAN class="punctuation token">:</SPAN>
&nbsp;&nbsp;&nbsp;&nbsp; <SPAN class="keyword token">for</SPAN> file <SPAN class="keyword token">in</SPAN> lstFiles<SPAN class="punctuation token">:</SPAN>
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; <SPAN class="keyword token">if</SPAN> fnmatch<SPAN class="punctuation token">.</SPAN>fnmatch<SPAN class="punctuation token">(</SPAN>file<SPAN class="punctuation token">,</SPAN> <SPAN class="string token">"H:\\PRISM_800m_weekly_sum_Normal"</SPAN> <SPAN class="operator token">+</SPAN> <SPAN class="string token">"*"</SPAN> <SPAN class="operator token">+</SPAN> <SPAN class="string token">"\\"</SPAN> <SPAN class="operator token">+</SPAN> <SPAN class="string token">"Week_1"</SPAN> <SPAN class="operator token">+</SPAN> <SPAN class="string token">"_*_"</SPAN> <SPAN class="operator token">+</SPAN> str<SPAN class="punctuation token">(</SPAN>month<SPAN class="punctuation token">)</SPAN> <SPAN class="operator token">+</SPAN> <SPAN class="string token">".tif"</SPAN><SPAN class="punctuation token">)</SPAN><SPAN class="punctuation token">:</SPAN>
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; new_list1<SPAN class="punctuation token">.</SPAN>append<SPAN class="punctuation token">(</SPAN>file<SPAN class="punctuation token">)</SPAN>
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; new_list1<SPAN class="punctuation token">.</SPAN>sort<SPAN class="punctuation token">(</SPAN><SPAN class="punctuation token">)</SPAN>
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; <SPAN class="keyword token">print</SPAN> <SPAN class="punctuation token">(</SPAN>new_list1<SPAN class="punctuation token">)</SPAN>
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 
&nbsp;&nbsp;&nbsp;&nbsp; env<SPAN class="punctuation token">.</SPAN>workspace <SPAN class="operator token">=</SPAN> out
&nbsp;&nbsp;&nbsp;&nbsp; finras1 <SPAN class="operator token">=</SPAN> out <SPAN class="operator token">+</SPAN> <SPAN class="string token">"\\"</SPAN> <SPAN class="operator token">+</SPAN> <SPAN class="string token">"Normal_Week_1"</SPAN> <SPAN class="operator token">+</SPAN> <SPAN class="string token">"_1981_2016"</SPAN> <SPAN class="operator token">+</SPAN> <SPAN class="string token">"_"</SPAN> <SPAN class="operator token">+</SPAN> str<SPAN class="punctuation token">(</SPAN>month<SPAN class="punctuation token">)</SPAN> <SPAN class="operator token">+</SPAN> <SPAN class="string token">".tif"</SPAN>
&nbsp;&nbsp;&nbsp;&nbsp; calc1 <SPAN class="operator token">=</SPAN> arcpy<SPAN class="punctuation token">.</SPAN>sa<SPAN class="punctuation token">.</SPAN>CellStatistics<SPAN class="punctuation token">(</SPAN>new_list1<SPAN class="punctuation token">,</SPAN> statistics_type <SPAN class="operator token">=</SPAN> <SPAN class="string token">"MEAN"</SPAN><SPAN class="punctuation token">)</SPAN>
&nbsp;&nbsp;&nbsp;&nbsp; arcpy<SPAN class="punctuation token">.</SPAN>CopyRaster_management<SPAN class="punctuation token">(</SPAN>calc1<SPAN class="punctuation token">,</SPAN>finras1<SPAN class="punctuation token">,</SPAN><SPAN class="string token">""</SPAN><SPAN class="punctuation token">,</SPAN><SPAN class="string token">""</SPAN><SPAN class="punctuation token">,</SPAN><SPAN class="string token">""</SPAN><SPAN class="punctuation token">,</SPAN><SPAN class="string token">""</SPAN><SPAN class="punctuation token">,</SPAN><SPAN class="string token">""</SPAN><SPAN class="punctuation token">,</SPAN><SPAN class="string token">"32_BIT_FLOAT"</SPAN><SPAN class="punctuation token">)</SPAN>

&nbsp;&nbsp;&nbsp;&nbsp; <SPAN class="keyword token">for</SPAN> file <SPAN class="keyword token">in</SPAN> lstFiles<SPAN class="punctuation token">:</SPAN>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; <SPAN class="keyword token">if</SPAN> fnmatch<SPAN class="punctuation token">.</SPAN>fnmatch<SPAN class="punctuation token">(</SPAN>file<SPAN class="punctuation token">,</SPAN> <SPAN class="string token">"H:\\PRISM_800m_weekly_sum_Normal"</SPAN> <SPAN class="operator token">+</SPAN> <SPAN class="string token">"*"</SPAN> <SPAN class="operator token">+</SPAN> <SPAN class="string token">"\\"</SPAN> <SPAN class="operator token">+</SPAN> <SPAN class="string token">"Week_2"</SPAN> <SPAN class="operator token">+</SPAN> <SPAN class="string token">"_*_"</SPAN> <SPAN class="operator token">+</SPAN> str<SPAN class="punctuation token">(</SPAN>month<SPAN class="punctuation token">)</SPAN> <SPAN class="operator token">+</SPAN> <SPAN class="string token">".tif"</SPAN><SPAN class="punctuation token">)</SPAN><SPAN class="punctuation token">:</SPAN>
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; new_list2<SPAN class="punctuation token">.</SPAN>append<SPAN class="punctuation token">(</SPAN>file<SPAN class="punctuation token">)</SPAN>
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; new_list2<SPAN class="punctuation token">.</SPAN>sort<SPAN class="punctuation token">(</SPAN><SPAN class="punctuation token">)</SPAN>
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; <SPAN class="keyword token">print</SPAN> <SPAN class="punctuation token">(</SPAN>new_list2<SPAN class="punctuation token">)</SPAN>
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 
&nbsp;&nbsp;&nbsp;&nbsp; env<SPAN class="punctuation token">.</SPAN>workspace <SPAN class="operator token">=</SPAN> out
&nbsp;&nbsp;&nbsp;&nbsp; finras2 <SPAN class="operator token">=</SPAN> out <SPAN class="operator token">+</SPAN> <SPAN class="string token">"\\"</SPAN> <SPAN class="operator token">+</SPAN> <SPAN class="string token">"Normal_Week_2"</SPAN> <SPAN class="operator token">+</SPAN> <SPAN class="string token">"_1981_2016"</SPAN> <SPAN class="operator token">+</SPAN> <SPAN class="string token">"_"</SPAN> <SPAN class="operator token">+</SPAN> str<SPAN class="punctuation token">(</SPAN>month<SPAN class="punctuation token">)</SPAN> <SPAN class="operator token">+</SPAN> <SPAN class="string token">".tif"</SPAN>
&nbsp;&nbsp;&nbsp;&nbsp; calc2 <SPAN class="operator token">=</SPAN> arcpy<SPAN class="punctuation token">.</SPAN>sa<SPAN class="punctuation token">.</SPAN>CellStatistics<SPAN class="punctuation token">(</SPAN>new_list2<SPAN class="punctuation token">,</SPAN> statistics_type <SPAN class="operator token">=</SPAN> <SPAN class="string token">"MEAN"</SPAN><SPAN class="punctuation token">)</SPAN>
&nbsp;&nbsp;&nbsp;&nbsp; arcpy<SPAN class="punctuation token">.</SPAN>CopyRaster_management<SPAN class="punctuation token">(</SPAN>calc2<SPAN class="punctuation token">,</SPAN>finras2<SPAN class="punctuation token">,</SPAN><SPAN class="string token">""</SPAN><SPAN class="punctuation token">,</SPAN><SPAN class="string token">""</SPAN><SPAN class="punctuation token">,</SPAN><SPAN class="string token">""</SPAN><SPAN class="punctuation token">,</SPAN><SPAN class="string token">""</SPAN><SPAN class="punctuation token">,</SPAN><SPAN class="string token">""</SPAN><SPAN class="punctuation token">,</SPAN><SPAN class="string token">"32_BIT_FLOAT"</SPAN><SPAN class="punctuation token">)</SPAN>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 
&nbsp;&nbsp;&nbsp;&nbsp; <SPAN class="keyword token">for</SPAN> file <SPAN class="keyword token">in</SPAN> lstFiles<SPAN class="punctuation token">:</SPAN>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; <SPAN class="keyword token">if</SPAN> fnmatch<SPAN class="punctuation token">.</SPAN>fnmatch<SPAN class="punctuation token">(</SPAN>file<SPAN class="punctuation token">,</SPAN> <SPAN class="string token">"H:\\PRISM_800m_weekly_sum_Normal"</SPAN> <SPAN class="operator token">+</SPAN> <SPAN class="string token">"*"</SPAN> <SPAN class="operator token">+</SPAN> <SPAN class="string token">"\\"</SPAN> <SPAN class="operator token">+</SPAN> <SPAN class="string token">"Week_3"</SPAN> <SPAN class="operator token">+</SPAN> <SPAN class="string token">"_*_"</SPAN> <SPAN class="operator token">+</SPAN> str<SPAN class="punctuation token">(</SPAN>month<SPAN class="punctuation token">)</SPAN> <SPAN class="operator token">+</SPAN> <SPAN class="string token">".tif"</SPAN><SPAN class="punctuation token">)</SPAN><SPAN class="punctuation token">:</SPAN>
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; new_list3<SPAN class="punctuation token">.</SPAN>append<SPAN class="punctuation token">(</SPAN>file<SPAN class="punctuation token">)</SPAN>
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; new_list3<SPAN class="punctuation token">.</SPAN>sort<SPAN class="punctuation token">(</SPAN><SPAN class="punctuation token">)</SPAN>
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; <SPAN class="keyword token">print</SPAN> <SPAN class="punctuation token">(</SPAN>new_list3<SPAN class="punctuation token">)</SPAN>
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 
&nbsp;&nbsp;&nbsp;&nbsp; env<SPAN class="punctuation token">.</SPAN>workspace <SPAN class="operator token">=</SPAN> out
&nbsp;&nbsp;&nbsp;&nbsp; finras3 <SPAN class="operator token">=</SPAN> out <SPAN class="operator token">+</SPAN> <SPAN class="string token">"\\"</SPAN> <SPAN class="operator token">+</SPAN> <SPAN class="string token">"Normal_Week_3"</SPAN> <SPAN class="operator token">+</SPAN> <SPAN class="string token">"_1981_2016"</SPAN> <SPAN class="operator token">+</SPAN> <SPAN class="string token">"_"</SPAN> <SPAN class="operator token">+</SPAN> str<SPAN class="punctuation token">(</SPAN>month<SPAN class="punctuation token">)</SPAN> <SPAN class="operator token">+</SPAN> <SPAN class="string token">".tif"</SPAN>
&nbsp;&nbsp;&nbsp;&nbsp; calc3 <SPAN class="operator token">=</SPAN> arcpy<SPAN class="punctuation token">.</SPAN>sa<SPAN class="punctuation token">.</SPAN>CellStatistics<SPAN class="punctuation token">(</SPAN>new_list3<SPAN class="punctuation token">,</SPAN> statistics_type <SPAN class="operator token">=</SPAN> <SPAN class="string token">"MEAN"</SPAN><SPAN class="punctuation token">)</SPAN>
&nbsp;&nbsp;&nbsp;&nbsp; arcpy<SPAN class="punctuation token">.</SPAN>CopyRaster_management<SPAN class="punctuation token">(</SPAN>calc3<SPAN class="punctuation token">,</SPAN>finras3<SPAN class="punctuation token">,</SPAN><SPAN class="string token">""</SPAN><SPAN class="punctuation token">,</SPAN><SPAN class="string token">""</SPAN><SPAN class="punctuation token">,</SPAN><SPAN class="string token">""</SPAN><SPAN class="punctuation token">,</SPAN><SPAN class="string token">""</SPAN><SPAN class="punctuation token">,</SPAN><SPAN class="string token">""</SPAN><SPAN class="punctuation token">,</SPAN><SPAN class="string token">"32_BIT_FLOAT"</SPAN><SPAN class="punctuation token">)</SPAN>
&nbsp;&nbsp;&nbsp;&nbsp; 
&nbsp;&nbsp;&nbsp;&nbsp; <SPAN class="keyword token">for</SPAN> file <SPAN class="keyword token">in</SPAN> lstFiles<SPAN class="punctuation token">:</SPAN>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; <SPAN class="keyword token">if</SPAN> fnmatch<SPAN class="punctuation token">.</SPAN>fnmatch<SPAN class="punctuation token">(</SPAN>file<SPAN class="punctuation token">,</SPAN> <SPAN class="string token">"H:\\PRISM_800m_weekly_sum_Normal"</SPAN> <SPAN class="operator token">+</SPAN> <SPAN class="string token">"*"</SPAN> <SPAN class="operator token">+</SPAN> <SPAN class="string token">"\\"</SPAN> <SPAN class="operator token">+</SPAN> <SPAN class="string token">"Week_4"</SPAN> <SPAN class="operator token">+</SPAN> <SPAN class="string token">"_*_"</SPAN> <SPAN class="operator token">+</SPAN> str<SPAN class="punctuation token">(</SPAN>month<SPAN class="punctuation token">)</SPAN> <SPAN class="operator token">+</SPAN> <SPAN class="string token">".tif"</SPAN><SPAN class="punctuation token">)</SPAN><SPAN class="punctuation token">:</SPAN>
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; new_list4<SPAN class="punctuation token">.</SPAN>append<SPAN class="punctuation token">(</SPAN>file<SPAN class="punctuation token">)</SPAN>
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; new_list4<SPAN class="punctuation token">.</SPAN>sort<SPAN class="punctuation token">(</SPAN><SPAN class="punctuation token">)</SPAN>
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; <SPAN class="keyword token">print</SPAN> <SPAN class="punctuation token">(</SPAN>new_list4<SPAN class="punctuation token">)</SPAN>
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 
&nbsp;&nbsp;&nbsp;&nbsp; env<SPAN class="punctuation token">.</SPAN>workspace <SPAN class="operator token">=</SPAN> out
&nbsp;&nbsp;&nbsp;&nbsp; finras4 <SPAN class="operator token">=</SPAN> out <SPAN class="operator token">+</SPAN> <SPAN class="string token">"\\"</SPAN> <SPAN class="operator token">+</SPAN> <SPAN class="string token">"Normal_Week_4"</SPAN> <SPAN class="operator token">+</SPAN> <SPAN class="string token">"_1981_2016"</SPAN> <SPAN class="operator token">+</SPAN> <SPAN class="string token">"_"</SPAN> <SPAN class="operator token">+</SPAN> str<SPAN class="punctuation token">(</SPAN>month<SPAN class="punctuation token">)</SPAN> <SPAN class="operator token">+</SPAN> <SPAN class="string token">".tif"</SPAN>
&nbsp;&nbsp;&nbsp;&nbsp; calc4 <SPAN class="operator token">=</SPAN> arcpy<SPAN class="punctuation token">.</SPAN>sa<SPAN class="punctuation token">.</SPAN>CellStatistics<SPAN class="punctuation token">(</SPAN>new_list4<SPAN class="punctuation token">,</SPAN> statistics_type <SPAN class="operator token">=</SPAN> <SPAN class="string token">"MEAN"</SPAN><SPAN class="punctuation token">)</SPAN>
&nbsp;&nbsp;&nbsp;&nbsp; arcpy<SPAN class="punctuation token">.</SPAN>CopyRaster_management<SPAN class="punctuation token">(</SPAN>calc4<SPAN class="punctuation token">,</SPAN>finras4<SPAN class="punctuation token">,</SPAN><SPAN class="string token">""</SPAN><SPAN class="punctuation token">,</SPAN><SPAN class="string token">""</SPAN><SPAN class="punctuation token">,</SPAN><SPAN class="string token">""</SPAN><SPAN class="punctuation token">,</SPAN><SPAN class="string token">""</SPAN><SPAN class="punctuation token">,</SPAN><SPAN class="string token">""</SPAN><SPAN class="punctuation token">,</SPAN><SPAN class="string token">"32_BIT_FLOAT"</SPAN><SPAN class="punctuation token">)</SPAN>
&nbsp;&nbsp;&nbsp;&nbsp; 
&nbsp;&nbsp;&nbsp;&nbsp; 
&nbsp;&nbsp;&nbsp;&nbsp; <SPAN class="keyword token">for</SPAN> file <SPAN class="keyword token">in</SPAN> lstFiles<SPAN class="punctuation token">:</SPAN>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; <SPAN class="keyword token">if</SPAN> fnmatch<SPAN class="punctuation token">.</SPAN>fnmatch<SPAN class="punctuation token">(</SPAN>file<SPAN class="punctuation token">,</SPAN> <SPAN class="string token">"H:\\PRISM_800m_weekly_sum_Normal"</SPAN> <SPAN class="operator token">+</SPAN> <SPAN class="string token">"*"</SPAN> <SPAN class="operator token">+</SPAN> <SPAN class="string token">"\\"</SPAN> <SPAN class="operator token">+</SPAN> <SPAN class="string token">"Week_5"</SPAN> <SPAN class="operator token">+</SPAN> <SPAN class="string token">"_*_"</SPAN> <SPAN class="operator token">+</SPAN> str<SPAN class="punctuation token">(</SPAN>month<SPAN class="punctuation token">)</SPAN> <SPAN class="operator token">+</SPAN> <SPAN class="string token">".tif"</SPAN><SPAN class="punctuation token">)</SPAN><SPAN class="punctuation token">:</SPAN>
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; new_list5<SPAN class="punctuation token">.</SPAN>append<SPAN class="punctuation token">(</SPAN>file<SPAN class="punctuation token">)</SPAN>
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; new_list5<SPAN class="punctuation token">.</SPAN>sort<SPAN class="punctuation token">(</SPAN><SPAN class="punctuation token">)</SPAN>
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; <SPAN class="keyword token">print</SPAN> <SPAN class="punctuation token">(</SPAN>new_list5<SPAN class="punctuation token">)</SPAN>&nbsp;&nbsp;&nbsp;&nbsp; 

&nbsp;&nbsp;&nbsp;&nbsp; env<SPAN class="punctuation token">.</SPAN>workspace <SPAN class="operator token">=</SPAN> out
&nbsp;&nbsp;&nbsp;&nbsp; finras5 <SPAN class="operator token">=</SPAN> out <SPAN class="operator token">+</SPAN> <SPAN class="string token">"\\"</SPAN> <SPAN class="operator token">+</SPAN> <SPAN class="string token">"Normal_Week_5"</SPAN> <SPAN class="operator token">+</SPAN> <SPAN class="string token">"_1981_2016"</SPAN> <SPAN class="operator token">+</SPAN> <SPAN class="string token">"_"</SPAN> <SPAN class="operator token">+</SPAN> str<SPAN class="punctuation token">(</SPAN>month<SPAN class="punctuation token">)</SPAN> <SPAN class="operator token">+</SPAN> <SPAN class="string token">".tif"</SPAN>
&nbsp;&nbsp;&nbsp;&nbsp; calc5 <SPAN class="operator token">=</SPAN> arcpy<SPAN class="punctuation token">.</SPAN>sa<SPAN class="punctuation token">.</SPAN>CellStatistics<SPAN class="punctuation token">(</SPAN>new_list5<SPAN class="punctuation token">,</SPAN> statistics_type <SPAN class="operator token">=</SPAN> <SPAN class="string token">"MEAN"</SPAN><SPAN class="punctuation token">)</SPAN>
&nbsp;&nbsp;&nbsp;&nbsp; arcpy<SPAN class="punctuation token">.</SPAN>CopyRaster_management<SPAN class="punctuation token">(</SPAN>calc5<SPAN class="punctuation token">,</SPAN>finras5<SPAN class="punctuation token">,</SPAN><SPAN class="string token">""</SPAN><SPAN class="punctuation token">,</SPAN><SPAN class="string token">""</SPAN><SPAN class="punctuation token">,</SPAN><SPAN class="string token">""</SPAN><SPAN class="punctuation token">,</SPAN><SPAN class="string token">""</SPAN><SPAN class="punctuation token">,</SPAN><SPAN class="string token">""</SPAN><SPAN class="punctuation token">,</SPAN><SPAN class="string token">"32_BIT_FLOAT"</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>‍</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><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>

タグ付けされた投稿をもっと見る

Accepted answers

RandyBurton

Glad that the output is showing the grouping that you want with your modification. I have added back in some of your earlier code, but I may have missed something.

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

<SPAN class="comment token"># make sure this directory exists</SPAN>
out <SPAN class="operator token">=</SPAN> r<SPAN class="string token">'H:\PRISM_800m_weekly_sum\Normal_out'</SPAN>

<SPAN class="keyword token">for</SPAN> week <SPAN class="keyword token">in</SPAN> range<SPAN class="punctuation token">(</SPAN><SPAN class="number token">1</SPAN><SPAN class="punctuation token">,</SPAN> <SPAN class="number token">5</SPAN><SPAN class="operator token">+</SPAN><SPAN class="number token">1</SPAN><SPAN class="punctuation token">)</SPAN><SPAN class="punctuation token">:</SPAN> <SPAN class="comment token"># for weeks 1 to 5</SPAN>

&nbsp;&nbsp;&nbsp; <SPAN class="keyword token">for</SPAN> month <SPAN class="keyword token">in</SPAN> range<SPAN class="punctuation token">(</SPAN><SPAN class="number token">9</SPAN><SPAN class="punctuation token">,</SPAN> <SPAN class="number token">11</SPAN><SPAN class="operator token">+</SPAN><SPAN class="number token">1</SPAN><SPAN class="punctuation token">)</SPAN><SPAN class="punctuation token">:</SPAN> <SPAN class="comment token"># for months 9 to 11</SPAN>

&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; rasterlist <SPAN class="operator token">=</SPAN> <SPAN class="punctuation token">[</SPAN><SPAN class="punctuation token">]</SPAN> <SPAN class="comment token"># initialize a list for the rasters</SPAN>
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; <SPAN class="keyword token">for</SPAN> year <SPAN class="keyword token">in</SPAN> range<SPAN class="punctuation token">(</SPAN><SPAN class="number token">1981</SPAN><SPAN class="punctuation token">,</SPAN> <SPAN class="number token">2016</SPAN><SPAN class="operator token">+</SPAN><SPAN class="number token">1</SPAN><SPAN class="punctuation token">)</SPAN><SPAN class="punctuation token">:</SPAN> <SPAN class="comment token"># for years 1981 to 2016</SPAN>

&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; rasterfile <SPAN class="operator token">=</SPAN> <SPAN class="string token">"H:\\PRISM_800m_weekly_sum_Normal\\{}\\Week_{}_Sum{}_{}.tif"</SPAN><SPAN class="punctuation token">.</SPAN>format<SPAN class="punctuation token">(</SPAN>year<SPAN class="punctuation token">,</SPAN> week<SPAN class="punctuation token">,</SPAN> year<SPAN class="punctuation token">,</SPAN> month<SPAN class="punctuation token">)</SPAN>
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; <SPAN class="comment token"># probably should add some error checking to see that file exists before adding it to rasterlist</SPAN>
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; <SPAN class="keyword token">if</SPAN> arcpy<SPAN class="punctuation token">.</SPAN>Exists<SPAN class="punctuation token">(</SPAN>rasterfile<SPAN class="punctuation token">)</SPAN><SPAN class="punctuation token">:</SPAN>
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; rasterlist<SPAN class="punctuation token">.</SPAN>append<SPAN class="punctuation token">(</SPAN>rasterfile<SPAN class="punctuation token">)</SPAN> <SPAN class="comment token"># if found, append file to the raster list</SPAN>
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; <SPAN class="keyword token">else</SPAN><SPAN class="punctuation token">:</SPAN>
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; <SPAN class="keyword token">print</SPAN> <SPAN class="string token">"Not found: {}"</SPAN><SPAN class="punctuation token">.</SPAN>format<SPAN class="punctuation token">(</SPAN>rasterfile<SPAN class="punctuation token">)</SPAN>

&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; <SPAN class="comment token"># same indentation as "for year in..."</SPAN>
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; <SPAN class="comment token"># rasterlist.sort() # you probably don't need to sort the list</SPAN>

&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; finras1 <SPAN class="operator token">=</SPAN> out <SPAN class="operator token">+</SPAN> <SPAN class="string token">"\\Normal_Week_{}_1981_2016_{}.tif"</SPAN><SPAN class="punctuation token">.</SPAN>format<SPAN class="punctuation token">(</SPAN>week<SPAN class="punctuation token">,</SPAN> month<SPAN class="punctuation token">)</SPAN>
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; <SPAN class="keyword token">print</SPAN> <SPAN class="string token">"Processing file: {}"</SPAN><SPAN class="punctuation token">.</SPAN>format<SPAN class="punctuation token">(</SPAN>finras1<SPAN class="punctuation token">)</SPAN>

&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; calc1 <SPAN class="operator token">=</SPAN> arcpy<SPAN class="punctuation token">.</SPAN>sa<SPAN class="punctuation token">.</SPAN>CellStatistics<SPAN class="punctuation token">(</SPAN>rasterlist<SPAN class="punctuation token">,</SPAN> statistics_type <SPAN class="operator token">=</SPAN> <SPAN class="string token">"MEAN"</SPAN><SPAN class="punctuation token">,</SPAN> ignore_nodata<SPAN class="operator token">=</SPAN><SPAN class="string token">"DATA"</SPAN><SPAN class="punctuation token">)</SPAN>
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; arcpy<SPAN class="punctuation token">.</SPAN>CopyRaster_management<SPAN class="punctuation token">(</SPAN>calc1<SPAN class="punctuation token">,</SPAN>finras1<SPAN class="punctuation token">,</SPAN><SPAN class="string token">""</SPAN><SPAN class="punctuation token">,</SPAN><SPAN class="string token">""</SPAN><SPAN class="punctuation token">,</SPAN><SPAN class="string token">""</SPAN><SPAN class="punctuation token">,</SPAN><SPAN class="string token">""</SPAN><SPAN class="punctuation token">,</SPAN><SPAN class="string token">""</SPAN><SPAN class="punctuation token">,</SPAN><SPAN class="string token">"32_BIT_FLOAT"</SPAN><SPAN class="punctuation token">)</SPAN>

<SPAN class="keyword token">print</SPAN> <SPAN class="string token">"Done."</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>‍</SPAN><SPAN>‍</SPAN></SPAN>

Also, keep in mind what Dan Patterson‌ said in his comments in another of your questions regarding working with rasters. Remember that indentation is important, and it is often good to have some print statements in your code to help you know where things are going right or wrong.

Hope this helps.

All comments

DuminduJayasekera

Randy, Thanks again. That is the output I want exactly but I would revise the last line like below. But, how can i modify the above code?

<SPAN class="" style="color: #0077aa; border: 0px; font-weight: inherit;">print</SPAN> <SPAN class="" style="color: #669900; border: 0px; font-weight: inherit;">"\tH:\\PRISM_800m_weekly_sum_Normal\\{}\\Week_{}_Sum{}_{}.tif"</SPAN><SPAN class="" style="color: #999999; border: 0px; font-weight: inherit;">.</SPAN>format<SPAN class="" style="color: #999999; border: 0px; font-weight: inherit;">(</SPAN>year<SPAN class="" style="color: #999999; border: 0px; font-weight: inherit;">,</SPAN> week<SPAN class="" style="color: #999999; border: 0px; font-weight: inherit;">,</SPAN> year<SPAN class="" style="color: #999999; border: 0px; font-weight: inherit;">,</SPAN> month<SPAN class="" style="color: #999999; border: 0px; font-weight: inherit;">)</SPAN>

RandyBurton

Does the following code print a listing of your files in the proper grouping for analysis? (I'm still trying to understand your directory structure and the objective.)

groupNum <SPAN class="operator token">=</SPAN> <SPAN class="number token">1</SPAN>
<SPAN class="keyword token">for</SPAN> week <SPAN class="keyword token">in</SPAN> range<SPAN class="punctuation token">(</SPAN><SPAN class="number token">1</SPAN><SPAN class="punctuation token">,</SPAN> <SPAN class="number token">5</SPAN><SPAN class="operator token">+</SPAN><SPAN class="number token">1</SPAN><SPAN class="punctuation token">)</SPAN><SPAN class="punctuation token">:</SPAN> <SPAN class="comment token"># for weeks 1 to 5</SPAN>
&nbsp;&nbsp;&nbsp; <SPAN class="keyword token">print</SPAN> <SPAN class="string token">"Files in group {}:"</SPAN><SPAN class="punctuation token">.</SPAN>format<SPAN class="punctuation token">(</SPAN>groupNum<SPAN class="punctuation token">)</SPAN>
&nbsp;&nbsp;&nbsp; groupNum <SPAN class="operator token">+=</SPAN> <SPAN class="number token">1</SPAN>
&nbsp;&nbsp;&nbsp; <SPAN class="keyword token">for</SPAN> month <SPAN class="keyword token">in</SPAN> range<SPAN class="punctuation token">(</SPAN><SPAN class="number token">9</SPAN><SPAN class="punctuation token">,</SPAN> <SPAN class="number token">11</SPAN><SPAN class="operator token">+</SPAN><SPAN class="number token">1</SPAN><SPAN class="punctuation token">)</SPAN><SPAN class="punctuation token">:</SPAN> <SPAN class="comment token"># for months 9 to 11</SPAN>
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; <SPAN class="keyword token">if</SPAN> month <SPAN class="operator token">></SPAN> <SPAN class="number token">9</SPAN><SPAN class="punctuation token">:</SPAN>
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; <SPAN class="keyword token">print</SPAN> <SPAN class="string token">"Files in group {}:"</SPAN><SPAN class="punctuation token">.</SPAN>format<SPAN class="punctuation token">(</SPAN>groupNum<SPAN class="punctuation token">)</SPAN>
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; groupNum <SPAN class="operator token">+=</SPAN> <SPAN class="number token">1</SPAN>
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; <SPAN class="keyword token">for</SPAN> year <SPAN class="keyword token">in</SPAN> range<SPAN class="punctuation token">(</SPAN><SPAN class="number token">1981</SPAN><SPAN class="punctuation token">,</SPAN> <SPAN class="number token">2016</SPAN><SPAN class="operator token">+</SPAN><SPAN class="number token">1</SPAN><SPAN class="punctuation token">)</SPAN><SPAN class="punctuation token">:</SPAN> <SPAN class="comment token"># for years 1981 to 2016</SPAN>
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; <SPAN class="keyword token">print</SPAN> <SPAN class="string token">"\tH:\\PRISM_800m_weekly_sum_Normal\\{}\\Week_{}\\Week_{}_Sum{}_{}.tif"</SPAN><SPAN class="punctuation token">.</SPAN>format<SPAN class="punctuation token">(</SPAN>year<SPAN class="punctuation token">,</SPAN> week<SPAN class="punctuation token">,</SPAN> week<SPAN class="punctuation token">,</SPAN> year<SPAN class="punctuation token">,</SPAN> month<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>

If the listing is correct, then we can modify the code to do the analysis.

DuminduJayasekera

Thanks Randy for the reply. I have revised my code and posted above. There is no error now and producing outputs but the output values are NOT correct. I am doing something wrong here.

I have folders names "1981", "1982"....."2016" in "H:\PRISM_800m_weekly_sum_Normal" and each folder (year) has Weekly .tif files (for example Week_1_Sum1981_9.tif, Week_2_Sum1981_9.tif, ...., Week_5_Sum1981_9.tif,.....Week_5_Sum1981_9.tif) for the months from 9 to 11.

Can you help me to fix this?

Thanks in advance.

RandyBurton

Hi Dumindu,

I agree with Dan Patterson‌ that you have some indentation issues with your code.

I may not understand exactly what you are trying to accomplish, but it looks like you want to average the raster files for the first week (week 1) of September (month 9) for the years 1981 to 2016. Then move to the next month and eventually loop through all the weeks. This is suggested by the file you are saving to found in line 19 of your code.

finras1 <SPAN class="operator token">=</SPAN> out <SPAN class="operator token">+</SPAN> <SPAN class="string token">"\\"</SPAN> <SPAN class="operator token">+</SPAN> <SPAN class="string token">"Normal_Week_"</SPAN> <SPAN class="operator token">+</SPAN> str<SPAN class="punctuation token">(</SPAN>week<SPAN class="punctuation token">)</SPAN> <SPAN class="operator token">+</SPAN> <SPAN class="string token">"_1981_2016"</SPAN> <SPAN class="operator token">+</SPAN> <SPAN class="string token">"_"</SPAN> <SPAN class="operator token">+</SPAN> str<SPAN class="punctuation token">(</SPAN>month<SPAN class="punctuation token">)</SPAN> <SPAN class="operator token">+</SPAN> <SPAN class="string token">".tif"</SPAN>
<SPAN class="line-numbers-rows"><SPAN>‍</SPAN></SPAN>

If so, you may have your loops in the wrong order. You may want something like this:

<SPAN class="keyword token">for</SPAN> week <SPAN class="keyword token">in</SPAN> range<SPAN class="punctuation token">(</SPAN><SPAN class="number token">1</SPAN><SPAN class="punctuation token">,</SPAN> <SPAN class="number token">5</SPAN><SPAN class="operator token">+</SPAN><SPAN class="number token">1</SPAN><SPAN class="punctuation token">)</SPAN><SPAN class="punctuation token">:</SPAN> <SPAN class="comment token"># for weeks 1 to 5</SPAN>
&nbsp;&nbsp;&nbsp; <SPAN class="keyword token">for</SPAN> month <SPAN class="keyword token">in</SPAN> range<SPAN class="punctuation token">(</SPAN><SPAN class="number token">9</SPAN><SPAN class="punctuation token">,</SPAN> <SPAN class="number token">11</SPAN><SPAN class="operator token">+</SPAN><SPAN class="number token">1</SPAN><SPAN class="punctuation token">)</SPAN><SPAN class="punctuation token">:</SPAN> <SPAN class="comment token"># for months 9 to 11</SPAN>
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; <SPAN class="keyword token">for</SPAN> year <SPAN class="keyword token">in</SPAN> range<SPAN class="punctuation token">(</SPAN><SPAN class="number token">1981</SPAN><SPAN class="punctuation token">,</SPAN> <SPAN class="number token">2016</SPAN><SPAN class="operator token">+</SPAN><SPAN class="number token">1</SPAN><SPAN class="punctuation token">)</SPAN><SPAN class="punctuation token">:</SPAN> <SPAN class="comment token"># for years 1981 to 2016</SPAN>
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; <SPAN class="keyword token">print</SPAN> <SPAN class="string token">"Week_{}_Sum{}_{}.tif"</SPAN><SPAN class="punctuation token">.</SPAN>format<SPAN class="punctuation token">(</SPAN>week<SPAN class="punctuation token">,</SPAN> year<SPAN class="punctuation token">,</SPAN> month<SPAN class="punctuation token">)</SPAN>
<SPAN class="string token">'''
Result:
Week_1_Sum1981_9.tif
Week_1_Sum1982_9.tif
...
Week_1_Sum2015_9.tif
Week_1_Sum2016_9.tif
Week_1_Sum1981_10.tif
Week_1_Sum1982_10.tif
...
Week_1_Sum2015_10.tif
Week_1_Sum2016_10.tif
Week_1_Sum1981_11.tif
Week_1_Sum1982_11.tif
...
Week_1_Sum2015_11.tif
Week_1_Sum2016_11.tif
Week_2_Sum1981_9.tif
Week_2_Sum1982_9.tif
...
...
Week_5_Sum2015_11.tif
Week_5_Sum2016_11.tif
'''</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>

If this is the case, you may want to change the directory structure of your rasters, with the weeks folders containing the months folders which contain the years folders. This would simplify what you are trying to do with arcpy.env.workspace and arcpy.ListRasters.

It would probably be helpful if you could describe the directory structure you are working with.

DanPatterson_Retired

indentation is wrong in line 17

lines 18 - end should be dedented to line up with the for in line 15

throw in some print statements

try using

arcpy.Exists( some file )

to see if your files actually exist and/or can be found

A continuation of this thread which you might want to close

会話に参加する

サインインしたメンバーは投稿、更新のフォローなどができます。初めてですか？無料アカウントを登録してください。

クイックリンク

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