I have five subfolders with Rasters name with the same parameter. I am trying to figure out a way to Loop average of all raster and save output raster in outfolder . For instance Subfolder1 Subfolder2 Subfolder3 Subfolder4 Subfolder5 Outfolder r001_NPP.tif r001_npp.tif r001_NPP.tif r_001_NPP.tif r001_NPP.tif r001_MEAN_NPP.tif I have tried to write code for this, Below i post my code what i have so far import arcpy , os , glob arcpy . CheckOutExtension ( "Spatial" ) ws1 = glob . glob ( r 'D:\subfolder1\*.tif' ) ws2 = glob . glob ( r 'D:\subfolder2\*.tif' ) ws3 = glob . glob ( r 'D:\subfolder3\*.tif' ) ws4 = glob . glob ( r 'D:\subfolder4\*.tif' ) ws5 = glob . glob ( r 'D:\subfolder5\*.tif' ) outws = r 'D:\outfolder' for r in ws1 , ws2 , ws3 , ws4 : basename = os . path . basename ( r ) . split ( "_" ) [ 0 ] r1 = arcpy . sa . Raster ( os . path . join ( ws1 , basename + "r_NPP.tif" ) ) r2 = arcpy . sa . Raster ( os . path . join ( ws2 , basename + "r_NPP.tif" ) ) r3 = arcpy . sa . Raster ( os . path . join ( ws3 , basename + "r_NPP.tif" ) ) r4 = arcpy . sa . Raster ( os . path . join ( ws4 , basename + "r_NPP.tif" ) ) r5 = arcpy . sa . Raster ( os . path . join ( ws5 , basename + "r_NPP.tif" ) ) result = ( r1 + r2 + r3 + r4 + r5 ) / 5 outname = basename + "r_Mean_NPP.tif" result . save ( os . path . join ( outws , outname ) )

I can't reproduce the error you mentioned. I did encounter some other errors that should be corrected like: * include import os at beginning of script and * change "out_workplace" into "out_workspace" import arcpy import os arcpy . CheckOutExtension ( "spatial" ) # generate list of rasters ##workspace = r"D:\TEST" ##out_workplace = r"D:\TEST\MEAN_NPP" workspace = r "C:\GeoNet\MeanRasters\Case20161009" out_workspace = r "C:\GeoNet\MeanRasters\Case20161009\MEAN_NPP" for ras_num in [ str ( i ) . zfill ( 3 ) for i in range ( 1 , 13 ) ] : print ras_num rasters = [ ] walk = arcpy . da . Walk ( workspace , datatype = "RasterDataset" , type = "TIF" ) for dirpath , dirnames , filenames in walk : for filename in filenames : if ras_num in filename : rasters . append ( os . path . join ( dirpath , filename ) ) # calculate mean if len ( rasters ) > 0 : print "Calculating mean of rasters: " , ', ' . join ( rasters ) ras_mean = arcpy . sa . CellStatistics ( rasters , "MEAN" ) ras_mean . save ( os . path . join ( out_workspace , 'r{0}_MEAN_NPP.TIF' . format ( ras_num ) ) ) else : print "no rasters for ras_num" , ras_num This created the mean raster for ras_num 001 to 005, see attached ZIP with resulting rasters

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

Perform Raster Calculation from Multiple Sub-folder Using ArcPy

ShouvikJha

I have five subfolders with Rasters name with the same parameter. I am trying to figure out a way to Loop average of all raster and save output raster in outfolder . For instance

Subfolder1&nbsp;&nbsp; Subfolder2&nbsp;&nbsp; Subfolder3&nbsp;&nbsp; Subfolder4&nbsp;&nbsp;&nbsp; Subfolder5&nbsp;&nbsp; Outfolder
r001_NPP.tif r001_npp.tif r001_NPP.tif r_001_NPP.tif r001_NPP.tif r001_MEAN_NPP.tif‍‍‍‍‍‍

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

I have tried to write code for this, Below i post my code what i have so far

<SPAN class="keyword token">import</SPAN> arcpy<SPAN class="punctuation token">,</SPAN> os<SPAN class="punctuation token">,</SPAN> glob
arcpy<SPAN class="punctuation token">.</SPAN>CheckOutExtension<SPAN class="punctuation token">(</SPAN><SPAN class="string token">"Spatial"</SPAN><SPAN class="punctuation token">)</SPAN>

ws1 <SPAN class="operator token">=</SPAN> glob<SPAN class="punctuation token">.</SPAN>glob<SPAN class="punctuation token">(</SPAN>r<SPAN class="string token">'D:\subfolder1\*.tif'</SPAN><SPAN class="punctuation token">)</SPAN>
ws2 <SPAN class="operator token">=</SPAN> glob<SPAN class="punctuation token">.</SPAN>glob<SPAN class="punctuation token">(</SPAN>r<SPAN class="string token">'D:\subfolder2\*.tif'</SPAN><SPAN class="punctuation token">)</SPAN>
ws3 <SPAN class="operator token">=</SPAN> glob<SPAN class="punctuation token">.</SPAN>glob<SPAN class="punctuation token">(</SPAN>r<SPAN class="string token">'D:\subfolder3\*.tif'</SPAN><SPAN class="punctuation token">)</SPAN>
ws4 <SPAN class="operator token">=</SPAN> glob<SPAN class="punctuation token">.</SPAN>glob<SPAN class="punctuation token">(</SPAN>r<SPAN class="string token">'D:\subfolder4\*.tif'</SPAN><SPAN class="punctuation token">)</SPAN>
ws5 <SPAN class="operator token">=</SPAN> glob<SPAN class="punctuation token">.</SPAN>glob<SPAN class="punctuation token">(</SPAN>r<SPAN class="string token">'D:\subfolder5\*.tif'</SPAN><SPAN class="punctuation token">)</SPAN>
outws <SPAN class="operator token">=</SPAN> r<SPAN class="string token">'D:\outfolder'</SPAN>

<SPAN class="keyword token">for</SPAN> r <SPAN class="keyword token">in</SPAN> ws1<SPAN class="punctuation token">,</SPAN> ws2<SPAN class="punctuation token">,</SPAN> ws3<SPAN class="punctuation token">,</SPAN> ws4<SPAN class="punctuation token">:</SPAN>
&nbsp;&nbsp;&nbsp; basename <SPAN class="operator token">=</SPAN> os<SPAN class="punctuation token">.</SPAN>path<SPAN class="punctuation token">.</SPAN>basename<SPAN class="punctuation token">(</SPAN>r<SPAN class="punctuation token">)</SPAN><SPAN class="punctuation token">.</SPAN>split<SPAN class="punctuation token">(</SPAN><SPAN class="string token">"_"</SPAN><SPAN class="punctuation token">)</SPAN><SPAN class="punctuation token">[</SPAN><SPAN class="number token">0</SPAN><SPAN class="punctuation token">]</SPAN>
&nbsp;&nbsp;&nbsp; r1 <SPAN class="operator token">=</SPAN> arcpy<SPAN class="punctuation token">.</SPAN>sa<SPAN class="punctuation token">.</SPAN>Raster<SPAN class="punctuation token">(</SPAN>os<SPAN class="punctuation token">.</SPAN>path<SPAN class="punctuation token">.</SPAN>join<SPAN class="punctuation token">(</SPAN>ws1<SPAN class="punctuation token">,</SPAN> basename <SPAN class="operator token">+</SPAN> <SPAN class="string token">"r_NPP.tif"</SPAN><SPAN class="punctuation token">)</SPAN><SPAN class="punctuation token">)</SPAN>
&nbsp;&nbsp;&nbsp; r2 <SPAN class="operator token">=</SPAN> arcpy<SPAN class="punctuation token">.</SPAN>sa<SPAN class="punctuation token">.</SPAN>Raster<SPAN class="punctuation token">(</SPAN>os<SPAN class="punctuation token">.</SPAN>path<SPAN class="punctuation token">.</SPAN>join<SPAN class="punctuation token">(</SPAN>ws2<SPAN class="punctuation token">,</SPAN> basename <SPAN class="operator token">+</SPAN> <SPAN class="string token">"r_NPP.tif"</SPAN><SPAN class="punctuation token">)</SPAN><SPAN class="punctuation token">)</SPAN>
&nbsp;&nbsp;&nbsp; r3 <SPAN class="operator token">=</SPAN> arcpy<SPAN class="punctuation token">.</SPAN>sa<SPAN class="punctuation token">.</SPAN>Raster<SPAN class="punctuation token">(</SPAN>os<SPAN class="punctuation token">.</SPAN>path<SPAN class="punctuation token">.</SPAN>join<SPAN class="punctuation token">(</SPAN>ws3<SPAN class="punctuation token">,</SPAN> basename <SPAN class="operator token">+</SPAN> <SPAN class="string token">"r_NPP.tif"</SPAN><SPAN class="punctuation token">)</SPAN><SPAN class="punctuation token">)</SPAN>
&nbsp;&nbsp;&nbsp; r4 <SPAN class="operator token">=</SPAN> arcpy<SPAN class="punctuation token">.</SPAN>sa<SPAN class="punctuation token">.</SPAN>Raster<SPAN class="punctuation token">(</SPAN>os<SPAN class="punctuation token">.</SPAN>path<SPAN class="punctuation token">.</SPAN>join<SPAN class="punctuation token">(</SPAN>ws4<SPAN class="punctuation token">,</SPAN> basename <SPAN class="operator token">+</SPAN> <SPAN class="string token">"r_NPP.tif"</SPAN><SPAN class="punctuation token">)</SPAN><SPAN class="punctuation token">)</SPAN>
&nbsp;&nbsp;&nbsp; r5 <SPAN class="operator token">=</SPAN> arcpy<SPAN class="punctuation token">.</SPAN>sa<SPAN class="punctuation token">.</SPAN>Raster<SPAN class="punctuation token">(</SPAN>os<SPAN class="punctuation token">.</SPAN>path<SPAN class="punctuation token">.</SPAN>join<SPAN class="punctuation token">(</SPAN>ws5<SPAN class="punctuation token">,</SPAN> basename <SPAN class="operator token">+</SPAN> <SPAN class="string token">"r_NPP.tif"</SPAN><SPAN class="punctuation token">)</SPAN><SPAN class="punctuation token">)</SPAN>
&nbsp;&nbsp;&nbsp; 

&nbsp;&nbsp;&nbsp; result <SPAN class="operator token">=</SPAN> <SPAN class="punctuation token">(</SPAN>r1 <SPAN class="operator token">+</SPAN> r2 <SPAN class="operator token">+</SPAN> r3 <SPAN class="operator token">+</SPAN> r4 <SPAN class="operator token">+</SPAN> r5<SPAN class="punctuation token">)</SPAN> <SPAN class="operator token">/</SPAN> <SPAN class="number token">5</SPAN>
&nbsp;&nbsp;&nbsp; outname <SPAN class="operator token">=</SPAN> basename <SPAN class="operator token">+</SPAN> <SPAN class="string token">"r_Mean_NPP.tif"</SPAN>

&nbsp;&nbsp;&nbsp; result<SPAN class="punctuation token">.</SPAN>save<SPAN class="punctuation token">(</SPAN>os<SPAN class="punctuation token">.</SPAN>path<SPAN class="punctuation token">.</SPAN>join<SPAN class="punctuation token">(</SPAN>outws<SPAN class="punctuation token">,</SPAN> outname<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>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN></SPAN>

Encuentra más publicaciones etiquetadas con

Accepted answers

XanderBakker

I can't reproduce the error you mentioned. I did encounter some other errors that should be corrected like:

include import os at beginning of script and
change "out_workplace" into "out_workspace"

<SPAN class="keyword token">import</SPAN> arcpy
<SPAN class="keyword token">import</SPAN> os
arcpy<SPAN class="punctuation token">.</SPAN>CheckOutExtension<SPAN class="punctuation token">(</SPAN><SPAN class="string token">"spatial"</SPAN><SPAN class="punctuation token">)</SPAN>

<SPAN class="comment token"># generate list of rasters</SPAN>
<SPAN class="comment token">##workspace = r"D:\TEST"</SPAN>
<SPAN class="comment token">##out_workplace = r"D:\TEST\MEAN_NPP"</SPAN>
workspace <SPAN class="operator token">=</SPAN> r<SPAN class="string token">"C:\GeoNet\MeanRasters\Case20161009"</SPAN>
out_workspace <SPAN class="operator token">=</SPAN> r<SPAN class="string token">"C:\GeoNet\MeanRasters\Case20161009\MEAN_NPP"</SPAN>

<SPAN class="keyword token">for</SPAN> ras_num <SPAN class="keyword token">in</SPAN> <SPAN class="punctuation token">[</SPAN>str<SPAN class="punctuation token">(</SPAN>i<SPAN class="punctuation token">)</SPAN><SPAN class="punctuation token">.</SPAN>zfill<SPAN class="punctuation token">(</SPAN><SPAN class="number token">3</SPAN><SPAN class="punctuation token">)</SPAN> <SPAN class="keyword token">for</SPAN> i <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">13</SPAN><SPAN class="punctuation token">)</SPAN><SPAN class="punctuation token">]</SPAN><SPAN class="punctuation token">:</SPAN>
&nbsp;&nbsp;&nbsp; <SPAN class="keyword token">print</SPAN> ras_num
&nbsp;&nbsp;&nbsp; rasters <SPAN class="operator token">=</SPAN> <SPAN class="punctuation token">[</SPAN><SPAN class="punctuation token">]</SPAN>
&nbsp;&nbsp;&nbsp; walk <SPAN class="operator token">=</SPAN> arcpy<SPAN class="punctuation token">.</SPAN>da<SPAN class="punctuation token">.</SPAN>Walk<SPAN class="punctuation token">(</SPAN>workspace<SPAN class="punctuation token">,</SPAN> datatype<SPAN class="operator token">=</SPAN><SPAN class="string token">"RasterDataset"</SPAN><SPAN class="punctuation token">,</SPAN> type<SPAN class="operator token">=</SPAN><SPAN class="string token">"TIF"</SPAN><SPAN class="punctuation token">)</SPAN>
&nbsp;&nbsp;&nbsp; <SPAN class="keyword token">for</SPAN> dirpath<SPAN class="punctuation token">,</SPAN> dirnames<SPAN class="punctuation token">,</SPAN> filenames <SPAN class="keyword token">in</SPAN> walk<SPAN class="punctuation token">:</SPAN>
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; <SPAN class="keyword token">for</SPAN> filename <SPAN class="keyword token">in</SPAN> filenames<SPAN class="punctuation token">:</SPAN>
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; <SPAN class="keyword token">if</SPAN> ras_num <SPAN class="keyword token">in</SPAN> filename<SPAN class="punctuation token">:</SPAN>
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; rasters<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>dirpath<SPAN class="punctuation token">,</SPAN> filename<SPAN class="punctuation token">)</SPAN><SPAN class="punctuation token">)</SPAN>

&nbsp;&nbsp;&nbsp; <SPAN class="comment token"># calculate mean</SPAN>
&nbsp;&nbsp;&nbsp; <SPAN class="keyword token">if</SPAN> len<SPAN class="punctuation token">(</SPAN>rasters<SPAN class="punctuation token">)</SPAN> <SPAN class="operator token">></SPAN> <SPAN class="number token">0</SPAN><SPAN class="punctuation token">:</SPAN>
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; <SPAN class="keyword token">print</SPAN> <SPAN class="string token">"Calculating mean of rasters: "</SPAN><SPAN class="punctuation token">,</SPAN> <SPAN class="string token">', '</SPAN><SPAN class="punctuation token">.</SPAN>join<SPAN class="punctuation token">(</SPAN>rasters<SPAN class="punctuation token">)</SPAN>
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; ras_mean <SPAN class="operator token">=</SPAN> arcpy<SPAN class="punctuation token">.</SPAN>sa<SPAN class="punctuation token">.</SPAN>CellStatistics<SPAN class="punctuation token">(</SPAN>rasters<SPAN class="punctuation token">,</SPAN> <SPAN class="string token">"MEAN"</SPAN><SPAN class="punctuation token">)</SPAN>
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; ras_mean<SPAN class="punctuation token">.</SPAN>save<SPAN class="punctuation token">(</SPAN>os<SPAN class="punctuation token">.</SPAN>path<SPAN class="punctuation token">.</SPAN>join<SPAN class="punctuation token">(</SPAN>out_workspace<SPAN class="punctuation token">,</SPAN> <SPAN class="string token">'r{0}_MEAN_NPP.TIF'</SPAN><SPAN class="punctuation token">.</SPAN>format<SPAN class="punctuation token">(</SPAN>ras_num<SPAN class="punctuation token">)</SPAN><SPAN class="punctuation token">)</SPAN><SPAN class="punctuation token">)</SPAN>
&nbsp;&nbsp;&nbsp; <SPAN class="keyword token">else</SPAN><SPAN class="punctuation token">:</SPAN>
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; <SPAN class="keyword token">print</SPAN> <SPAN class="string token">"no rasters for ras_num"</SPAN><SPAN class="punctuation token">,</SPAN> ras_num‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍<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>

This created the mean raster for ras_num 001 to 005, see attached ZIP with resulting rasters

MEAN_NPP.zip

All comments

DanPatterson_Retired

No problem Xander, I will look for it

ShouvikJha

xander_bakker‌, Thank you for the suggestion. I am extremely sorry. Next time surely i will follow your suggestion. Sorry for the inconvenience.

XanderBakker

Hi shouvik jha ,

When you ask a new question, could you please create a new question? In case it has a relation with a previous question, include the link to that question. I branched the question to a new thread: Scatter plot of two rasters . However, you replied to a different level, so the last comments have not been passed to the new thread. You should not mix questions in the same thread. In Dan would resolved this question, and you would mark it correct, people looking at this thread due to the title will be completely confused.

Please move your posts to the new thread.

Sorry Dan Patterson for the inconvenience.

DanPatterson_Retired

If the rasters aren't huge convert to points, then Ordinary Least Squares (OLS)—Help | ArcGIS for Desktop I think the report allows for a graph

ShouvikJha

Thank you. I am not getting your points. i am trying plot two variable of raster like this ...

XanderBakker

I'm glad it's working now...

ShouvikJha

xander_bakker‌. Thank you very much. Above Code running perfectly.

Regarding error issue which you got different error massage, because i corrected that line before running script.

XanderBakker

I can't reproduce the error you mentioned. I did encounter some other errors that should be corrected like:

include import os at beginning of script and
change "out_workplace" into "out_workspace"

<SPAN class="keyword token">import</SPAN> arcpy
<SPAN class="keyword token">import</SPAN> os
arcpy<SPAN class="punctuation token">.</SPAN>CheckOutExtension<SPAN class="punctuation token">(</SPAN><SPAN class="string token">"spatial"</SPAN><SPAN class="punctuation token">)</SPAN>

<SPAN class="comment token"># generate list of rasters</SPAN>
<SPAN class="comment token">##workspace = r"D:\TEST"</SPAN>
<SPAN class="comment token">##out_workplace = r"D:\TEST\MEAN_NPP"</SPAN>
workspace <SPAN class="operator token">=</SPAN> r<SPAN class="string token">"C:\GeoNet\MeanRasters\Case20161009"</SPAN>
out_workspace <SPAN class="operator token">=</SPAN> r<SPAN class="string token">"C:\GeoNet\MeanRasters\Case20161009\MEAN_NPP"</SPAN>

<SPAN class="keyword token">for</SPAN> ras_num <SPAN class="keyword token">in</SPAN> <SPAN class="punctuation token">[</SPAN>str<SPAN class="punctuation token">(</SPAN>i<SPAN class="punctuation token">)</SPAN><SPAN class="punctuation token">.</SPAN>zfill<SPAN class="punctuation token">(</SPAN><SPAN class="number token">3</SPAN><SPAN class="punctuation token">)</SPAN> <SPAN class="keyword token">for</SPAN> i <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">13</SPAN><SPAN class="punctuation token">)</SPAN><SPAN class="punctuation token">]</SPAN><SPAN class="punctuation token">:</SPAN>
&nbsp;&nbsp;&nbsp; <SPAN class="keyword token">print</SPAN> ras_num
&nbsp;&nbsp;&nbsp; rasters <SPAN class="operator token">=</SPAN> <SPAN class="punctuation token">[</SPAN><SPAN class="punctuation token">]</SPAN>
&nbsp;&nbsp;&nbsp; walk <SPAN class="operator token">=</SPAN> arcpy<SPAN class="punctuation token">.</SPAN>da<SPAN class="punctuation token">.</SPAN>Walk<SPAN class="punctuation token">(</SPAN>workspace<SPAN class="punctuation token">,</SPAN> datatype<SPAN class="operator token">=</SPAN><SPAN class="string token">"RasterDataset"</SPAN><SPAN class="punctuation token">,</SPAN> type<SPAN class="operator token">=</SPAN><SPAN class="string token">"TIF"</SPAN><SPAN class="punctuation token">)</SPAN>
&nbsp;&nbsp;&nbsp; <SPAN class="keyword token">for</SPAN> dirpath<SPAN class="punctuation token">,</SPAN> dirnames<SPAN class="punctuation token">,</SPAN> filenames <SPAN class="keyword token">in</SPAN> walk<SPAN class="punctuation token">:</SPAN>
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; <SPAN class="keyword token">for</SPAN> filename <SPAN class="keyword token">in</SPAN> filenames<SPAN class="punctuation token">:</SPAN>
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; <SPAN class="keyword token">if</SPAN> ras_num <SPAN class="keyword token">in</SPAN> filename<SPAN class="punctuation token">:</SPAN>
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; rasters<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>dirpath<SPAN class="punctuation token">,</SPAN> filename<SPAN class="punctuation token">)</SPAN><SPAN class="punctuation token">)</SPAN>

&nbsp;&nbsp;&nbsp; <SPAN class="comment token"># calculate mean</SPAN>
&nbsp;&nbsp;&nbsp; <SPAN class="keyword token">if</SPAN> len<SPAN class="punctuation token">(</SPAN>rasters<SPAN class="punctuation token">)</SPAN> <SPAN class="operator token">></SPAN> <SPAN class="number token">0</SPAN><SPAN class="punctuation token">:</SPAN>
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; <SPAN class="keyword token">print</SPAN> <SPAN class="string token">"Calculating mean of rasters: "</SPAN><SPAN class="punctuation token">,</SPAN> <SPAN class="string token">', '</SPAN><SPAN class="punctuation token">.</SPAN>join<SPAN class="punctuation token">(</SPAN>rasters<SPAN class="punctuation token">)</SPAN>
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; ras_mean <SPAN class="operator token">=</SPAN> arcpy<SPAN class="punctuation token">.</SPAN>sa<SPAN class="punctuation token">.</SPAN>CellStatistics<SPAN class="punctuation token">(</SPAN>rasters<SPAN class="punctuation token">,</SPAN> <SPAN class="string token">"MEAN"</SPAN><SPAN class="punctuation token">)</SPAN>
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; ras_mean<SPAN class="punctuation token">.</SPAN>save<SPAN class="punctuation token">(</SPAN>os<SPAN class="punctuation token">.</SPAN>path<SPAN class="punctuation token">.</SPAN>join<SPAN class="punctuation token">(</SPAN>out_workspace<SPAN class="punctuation token">,</SPAN> <SPAN class="string token">'r{0}_MEAN_NPP.TIF'</SPAN><SPAN class="punctuation token">.</SPAN>format<SPAN class="punctuation token">(</SPAN>ras_num<SPAN class="punctuation token">)</SPAN><SPAN class="punctuation token">)</SPAN><SPAN class="punctuation token">)</SPAN>
&nbsp;&nbsp;&nbsp; <SPAN class="keyword token">else</SPAN><SPAN class="punctuation token">:</SPAN>
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; <SPAN class="keyword token">print</SPAN> <SPAN class="string token">"no rasters for ras_num"</SPAN><SPAN class="punctuation token">,</SPAN> ras_num‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍<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>

This created the mean raster for ras_num 001 to 005, see attached ZIP with resulting rasters

MEAN_NPP.zip

curtvprice

Hey no fair. 😉

ShouvikJha

Thank you very much all of you for your kind cooperation.

Xander Bakker I ran the above code , but i am getting error. I think not having problem with code. Might be we need more explanation what i am trying to do, please follow the below

Message&nbsp;&nbsp;&nbsp;&nbsp; File Name&nbsp;&nbsp;&nbsp;&nbsp; Line&nbsp;&nbsp;&nbsp;&nbsp; Position&nbsp;&nbsp;&nbsp;&nbsp; 
Traceback&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 
&nbsp;&nbsp;&nbsp; <SPAN class="operator token"><</SPAN>module<SPAN class="operator token">></SPAN>&nbsp;&nbsp;&nbsp;&nbsp; <SPAN class="operator token"><</SPAN>module1<SPAN class="operator token">></SPAN>&nbsp;&nbsp;&nbsp;&nbsp; <SPAN class="number token">20</SPAN>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 
&nbsp;&nbsp;&nbsp; CellStatistics&nbsp;&nbsp;&nbsp;&nbsp; C<SPAN class="punctuation token">:</SPAN>\Program Files\ArcGIS\Desktop10<SPAN class="number token">.3</SPAN>\ArcPy\arcpy\sa\Functions<SPAN class="punctuation token">.</SPAN>py&nbsp;&nbsp;&nbsp;&nbsp; <SPAN class="number token">2951</SPAN>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 
&nbsp;&nbsp;&nbsp; swapper&nbsp;&nbsp;&nbsp;&nbsp; C<SPAN class="punctuation token">:</SPAN>\Program Files\ArcGIS\Desktop10<SPAN class="number token">.3</SPAN>\ArcPy\arcpy\sa\Utils<SPAN class="punctuation token">.</SPAN>py&nbsp;&nbsp;&nbsp;&nbsp; <SPAN class="number token">53</SPAN>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 
&nbsp;&nbsp;&nbsp; Wrapper&nbsp;&nbsp;&nbsp;&nbsp; C<SPAN class="punctuation token">:</SPAN>\Program Files\ArcGIS\Desktop10<SPAN class="number token">.3</SPAN>\ArcPy\arcpy\sa\Functions<SPAN class="punctuation token">.</SPAN>py&nbsp;&nbsp;&nbsp;&nbsp; <SPAN class="number token">2947</SPAN>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 
RuntimeError<SPAN class="punctuation token">:</SPAN> ERROR <SPAN class="number token">999998</SPAN><SPAN class="punctuation token">:</SPAN> Unexpected Error<SPAN class="punctuation token">.</SPAN>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 
‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍<SPAN class="line-numbers-rows"><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN></SPAN>

I am giving here what I am trying to do, I have three sub-folder (Name : NPP_2008, NPP_2009, NPP_2010),

r001_MEAN_NPP.tif (mean of all r001_NPP)<CODE>
<CODE>r002_NPP.tif r002_npp.tif&nbsp;&nbsp;&nbsp; r002_NPP.tif&nbsp;&nbsp; <CODE>r002_MEAN_NPP.tif <SPAN style="color: #800000;">(mean of all r002_NPP)</SPAN>
So on ..............

<CODE><CODE><SPAN style="color: #0000ff;">NPP_2008&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; NPP_2009&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; NPP_2010&nbsp;</SPAN>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; <SPAN style="color: #008000;">MEAN_NPP (output Folder)</SPAN>
<SPAN style="color: #ff00ff;">r001_NPP.tif&nbsp; r001_npp.tif&nbsp;&nbsp; r001_NPP.tif&nbsp;&nbsp; </SPAN>

Every sub-folder hold monthly raster data for different year with the name of e.g r001_NPP, r002_NPP……..r012_NPP. So I want calculate mean of those raster which is r001_NPP series in every sub-folder, similarly program will loop through till r012_NPP series of all raster to calculate mean and save it another subfolder (Name : Mean_NPP) with output name 001_MEAN_NPP(which is mean of all r001_NPP), 002_MEAN NPP…..012_MEAN_NPP.

Here with i have attached reference raster which i am using for my raster calculation

Reference_RASTER.rar

XanderBakker

Absolutely Dan, you are right!

DanPatterson_Retired

Xander, thanks.... line 7, perhaps...

now before people get sloppy and forget to use 'raw' or unix based folder formatting... lets get those little 'r' s in there, especially if people get sloppy and/or use lower case, or upper case gets translated to lower case.... sooo. a 'case' in point...

<SPAN class="operator token">>></SPAN><SPAN class="operator token">></SPAN> w <SPAN class="operator token">=</SPAN> <SPAN class="string token">'D:\Test\Folder'</SPAN>
<SPAN class="operator token">>></SPAN><SPAN class="operator token">></SPAN> <SPAN class="keyword token">print</SPAN><SPAN class="punctuation token">(</SPAN>w<SPAN class="punctuation token">)</SPAN>
D<SPAN class="punctuation token">:</SPAN>\Test\Folder
<SPAN class="operator token">>></SPAN><SPAN class="operator token">></SPAN>
<SPAN class="operator token">>></SPAN><SPAN class="operator token">></SPAN> <SPAN class="comment token"># now let's get sloppy</SPAN>
<SPAN class="operator token">>></SPAN><SPAN class="operator token">></SPAN>
<SPAN class="operator token">>></SPAN><SPAN class="operator token">></SPAN> w <SPAN class="operator token">=</SPAN> <SPAN class="string token">'D:\test\folder'</SPAN>
<SPAN class="operator token">>></SPAN><SPAN class="operator token">></SPAN> <SPAN class="keyword token">print</SPAN><SPAN class="punctuation token">(</SPAN>w<SPAN class="punctuation token">)</SPAN>
D<SPAN class="punctuation token">:</SPAN>&nbsp;&nbsp;&nbsp;&nbsp; est 
older‍‍‍‍‍‍‍‍‍‍<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>

XanderBakker

I'm looking at the code and trying to figure what you want to do exactly, since the indentation is not correct.

First of all I assume you want to perform this calculation for each month (001, .. , 012). To do so, you will need increase the range to (1, 13) (see line 9 below).

I also assume that you do not want to have the same output results for each month and that the list of rasters should be filtered by the month value. Therefore I included line 14.

Please be aware that in you original code you create the raster r001_MEAN_NPP.TIF and since that raster is stored in a subfolder of the input workspace, it will be included in the result for the second month. I don't think that is what you want. Also remember your closing brackets on line 19 (there were two missing in your code)

<SPAN class="keyword token">import</SPAN> arcpy
arcpy<SPAN class="punctuation token">.</SPAN>CheckOutExtension<SPAN class="punctuation token">(</SPAN><SPAN class="string token">"spatial"</SPAN><SPAN class="punctuation token">)</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>

<SPAN class="comment token"># generate list of rasters</SPAN>
workspace <SPAN class="operator token">=</SPAN> r<SPAN class="string token">"D:\TEST"</SPAN>
out_workplace <SPAN class="operator token">=</SPAN> <SPAN class="string token">"D:\TEST\MEAN_NPP"</SPAN>

<SPAN class="keyword token">for</SPAN> ras_num <SPAN class="keyword token">in</SPAN> <SPAN class="punctuation token">[</SPAN>str<SPAN class="punctuation token">(</SPAN>i<SPAN class="punctuation token">)</SPAN><SPAN class="punctuation token">.</SPAN>zfill<SPAN class="punctuation token">(</SPAN><SPAN class="number token">3</SPAN><SPAN class="punctuation token">)</SPAN> <SPAN class="keyword token">for</SPAN> i <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">13</SPAN><SPAN class="punctuation token">)</SPAN><SPAN class="punctuation token">]</SPAN><SPAN class="punctuation token">:</SPAN>
&nbsp;&nbsp;&nbsp; rasters <SPAN class="operator token">=</SPAN> <SPAN class="punctuation token">[</SPAN><SPAN class="punctuation token">]</SPAN>
&nbsp;&nbsp;&nbsp; walk <SPAN class="operator token">=</SPAN> arcpy<SPAN class="punctuation token">.</SPAN>da<SPAN class="punctuation token">.</SPAN>Walk<SPAN class="punctuation token">(</SPAN>workspace<SPAN class="punctuation token">,</SPAN> datatype<SPAN class="operator token">=</SPAN><SPAN class="string token">"RasterDataset"</SPAN><SPAN class="punctuation token">,</SPAN> type<SPAN class="operator token">=</SPAN><SPAN class="string token">"TIF"</SPAN><SPAN class="punctuation token">)</SPAN>
&nbsp;&nbsp;&nbsp; <SPAN class="keyword token">for</SPAN> dirpath<SPAN class="punctuation token">,</SPAN> dirnames<SPAN class="punctuation token">,</SPAN> filenames <SPAN class="keyword token">in</SPAN> walk<SPAN class="punctuation token">:</SPAN>
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; <SPAN class="keyword token">for</SPAN> filename <SPAN class="keyword token">in</SPAN> filenames<SPAN class="punctuation token">:</SPAN>
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; <SPAN class="keyword token">if</SPAN> ras_num <SPAN class="keyword token">in</SPAN> filename<SPAN class="punctuation token">:</SPAN>
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; rasters<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>dirpath<SPAN class="punctuation token">,</SPAN> filename<SPAN class="punctuation token">)</SPAN><SPAN class="punctuation token">)</SPAN>

&nbsp;&nbsp;&nbsp; <SPAN class="comment token"># calculate mean</SPAN>
&nbsp;&nbsp;&nbsp; ras_mean <SPAN class="operator token">=</SPAN> CellStatistics<SPAN class="punctuation token">(</SPAN>rasters<SPAN class="punctuation token">,</SPAN> <SPAN class="string token">"MEAN"</SPAN><SPAN class="punctuation token">)</SPAN>
&nbsp;&nbsp;&nbsp; ras_mean<SPAN class="punctuation token">.</SPAN>save<SPAN class="punctuation token">(</SPAN>os<SPAN class="punctuation token">.</SPAN>path<SPAN class="punctuation token">.</SPAN>join<SPAN class="punctuation token">(</SPAN>out_workspace<SPAN class="punctuation token">,</SPAN> <SPAN class="string token">'r{0}_MEAN_NPP.TIF'</SPAN><SPAN class="punctuation token">.</SPAN>format<SPAN class="punctuation token">(</SPAN>ras_num<SPAN class="punctuation token">)</SPAN>‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍<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>

XanderBakker

In addition, line 18 is missing 2 closing brackets "))" at the end of the line. Perhaps you didn't copy the entire line.

DanPatterson_Retired

Here is how you test things... when you get an error, narrow down the code and test it, like in the following example to ensure that you output file names are correct..

<SPAN class="operator token">>></SPAN><SPAN class="operator token">></SPAN> <SPAN class="keyword token">import</SPAN> os
<SPAN class="operator token">>></SPAN><SPAN class="operator token">></SPAN> out_workplace <SPAN class="operator token">=</SPAN> <SPAN class="string token">'D:/Test/MEAN_NPP'</SPAN>&nbsp; <SPAN class="comment token"># note I am using os independent path specs</SPAN>
<SPAN class="operator token">>></SPAN><SPAN class="operator token">></SPAN> num <SPAN class="operator token">=</SPAN> <SPAN class="number token">1</SPAN>&nbsp;&nbsp; <SPAN class="comment token"># obviously if you had this in a loop 'num' would be changing</SPAN>
<SPAN class="operator token">>></SPAN><SPAN class="operator token">></SPAN> p <SPAN class="operator token">=</SPAN> <SPAN class="string token">'r{0}_MEAN_NPP.TIF'</SPAN><SPAN class="punctuation token">.</SPAN>format<SPAN class="punctuation token">(</SPAN>num<SPAN class="punctuation token">)</SPAN>&nbsp;&nbsp; <SPAN class="comment token"># this is the path I am trying to make</SPAN>
<SPAN class="operator token">>></SPAN><SPAN class="operator token">></SPAN> x <SPAN class="operator token">=</SPAN> os<SPAN class="punctuation token">.</SPAN>path<SPAN class="punctuation token">.</SPAN>join<SPAN class="punctuation token">(</SPAN>out_workplace<SPAN class="punctuation token">,</SPAN> p<SPAN class="punctuation token">)</SPAN>&nbsp; <SPAN class="comment token"># joining them together</SPAN>
<SPAN class="operator token">>></SPAN><SPAN class="operator token">></SPAN> x
<SPAN class="string token">'D:/Test/MEAN_NPP/r1_MEAN_NPP.TIF'</SPAN>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; <SPAN class="comment token"># this is what we get</SPAN>
<SPAN class="line-numbers-rows"><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN></SPAN>

Now... is that right?

RebeccaStrauch__GISP

Change line 7 to

out_workplace = r"D:\TEST\MEAN_NPP"

You missed the "r" in that path.

ShouvikJha

Dan , thanks . I have solved the indentation error, and i updated my code but still i could not solve the problem, updated code below

error masage

ile <SPAN class="string token">"<module1>"</SPAN><SPAN class="punctuation token">,</SPAN> line <SPAN class="number token">18</SPAN>
&nbsp;&nbsp;&nbsp; sum<SPAN class="punctuation token">.</SPAN>save<SPAN class="punctuation token">(</SPAN>os<SPAN class="punctuation token">.</SPAN>path<SPAN class="punctuation token">.</SPAN>join<SPAN class="punctuation token">(</SPAN>out_workspace<SPAN class="punctuation token">,</SPAN> <SPAN class="string token">'r{0}_MEAN_NPP.TIF'</SPAN><SPAN class="punctuation token">.</SPAN>format<SPAN class="punctuation token">(</SPAN>ras_num<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;&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="operator token">^</SPAN>
SyntaxError<SPAN class="punctuation token">:</SPAN> invalid syntax<SPAN class="line-numbers-rows"><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN></SPAN>

Updated Code

<SPAN class="keyword token">import</SPAN> arcpy
arcpy<SPAN class="punctuation token">.</SPAN>CheckOutExtension<SPAN class="punctuation token">(</SPAN><SPAN class="string token">"spatial"</SPAN><SPAN class="punctuation token">)</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>

<SPAN class="comment token"># generate list of rasters</SPAN>
workspace <SPAN class="operator token">=</SPAN> r<SPAN class="string token">"D:\TEST"</SPAN>
out_workplace <SPAN class="operator token">=</SPAN> <SPAN class="string token">"D:\TEST\MEAN_NPP"</SPAN>

<SPAN class="keyword token">for</SPAN> ras_num <SPAN class="keyword token">in</SPAN> <SPAN class="punctuation token">[</SPAN>str<SPAN class="punctuation token">(</SPAN>i<SPAN class="punctuation token">)</SPAN><SPAN class="punctuation token">.</SPAN>zfill<SPAN class="punctuation token">(</SPAN><SPAN class="number token">3</SPAN><SPAN class="punctuation token">)</SPAN> <SPAN class="keyword token">for</SPAN> i <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">12</SPAN><SPAN class="punctuation token">)</SPAN><SPAN class="punctuation token">]</SPAN><SPAN class="punctuation token">:</SPAN>
 rasters <SPAN class="operator token">=</SPAN> <SPAN class="punctuation token">[</SPAN><SPAN class="punctuation token">]</SPAN>
walk <SPAN class="operator token">=</SPAN> arcpy<SPAN class="punctuation token">.</SPAN>da<SPAN class="punctuation token">.</SPAN>Walk<SPAN class="punctuation token">(</SPAN>workspace<SPAN class="punctuation token">,</SPAN> datatype<SPAN class="operator token">=</SPAN><SPAN class="string token">"RasterDataset"</SPAN><SPAN class="punctuation token">,</SPAN> type<SPAN class="operator token">=</SPAN><SPAN class="string token">"TIF"</SPAN><SPAN class="punctuation token">)</SPAN>
<SPAN class="keyword token">for</SPAN> dirpath<SPAN class="punctuation token">,</SPAN> dirnames<SPAN class="punctuation token">,</SPAN> filenames <SPAN class="keyword token">in</SPAN> walk<SPAN class="punctuation token">:</SPAN>
&nbsp;&nbsp;&nbsp; <SPAN class="keyword token">for</SPAN> filename <SPAN class="keyword token">in</SPAN> filenames<SPAN class="punctuation token">:</SPAN>
&nbsp;&nbsp;&nbsp;&nbsp; rasters<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>dirpath<SPAN class="punctuation token">,</SPAN> filename<SPAN class="punctuation token">)</SPAN><SPAN class="punctuation token">)</SPAN>

<SPAN class="comment token"># calculate mean</SPAN>
sum <SPAN class="operator token">=</SPAN> CellStatistics<SPAN class="punctuation token">(</SPAN>rasters<SPAN class="punctuation token">,</SPAN> <SPAN class="string token">"MEAN"</SPAN><SPAN class="punctuation token">)</SPAN>
sum<SPAN class="punctuation token">.</SPAN>save<SPAN class="punctuation token">(</SPAN>os<SPAN class="punctuation token">.</SPAN>path<SPAN class="punctuation token">.</SPAN>join<SPAN class="punctuation token">(</SPAN>out_workspace<SPAN class="punctuation token">,</SPAN> <SPAN class="string token">'r{0}_MEAN_NPP.TIF'</SPAN><SPAN class="punctuation token">.</SPAN>format<SPAN class="punctuation token">(</SPAN>ras_num<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>

DanPatterson_Retired

Ok.... add four (4) more spaces to the beginning of line 13 so it is indented properly

Here is the python style guide recommendations PEP 8 -- Style Guide for Python Code | Python.org

ShouvikJha

xander_bakker‌, We need your cooperation to solve above problem. Kindly cooperate with us.

ShouvikJha

Dan . thanks , you are right. I didn't post the error massage. Line no 13 i am getting the Indentation Error massage. I tried to fix it but not solved yet. error massage below

File <SPAN class="string token">"<module1>"</SPAN><SPAN class="punctuation token">,</SPAN> line <SPAN class="number token">13</SPAN>
&nbsp;&nbsp;&nbsp; rasters<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>dirpath<SPAN class="punctuation token">,</SPAN> filename<SPAN class="punctuation token">)</SPAN><SPAN class="punctuation token">)</SPAN>
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; <SPAN class="operator token">^</SPAN>
IndentationError<SPAN class="punctuation token">:</SPAN> expected an indented block
<SPAN class="line-numbers-rows"><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN></SPAN>

DanPatterson_Retired

your code is only minimally updated from what Neil provided and it was untested.

What was the error message?

issues to get you started...

Line 13 is not indented properly
Walk—Help | ArcGIS for Desktop you need to check the correct format
datatype and type in Walk are incorrect
you shouldn't use from X import *

You need to provide some input to the solution and to do that providing commentary as to what went wrong and where would go a long way

ShouvikJha

Thank you. I have updated my suggested by you, problem still persist, i am not understanding how to iterate through each folder, take all series of 001_NPP from all folder and calculate average oh them and save it as 001_NPP_MEAN, below is my updated code and along with code i have attached reference raster also which i am using for my raster calculation,

<SPAN class="keyword token">import</SPAN> arcpy
arcpy<SPAN class="punctuation token">.</SPAN>CheckOutExtension<SPAN class="punctuation token">(</SPAN><SPAN class="string token">"spatial"</SPAN><SPAN class="punctuation token">)</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>

<SPAN class="comment token"># generate list of rasters</SPAN>
workspace <SPAN class="operator token">=</SPAN> r<SPAN class="string token">"D:\TEST"</SPAN>
out_workplace <SPAN class="operator token">=</SPAN> <SPAN class="string token">"D:\TEST\MEAN_NPP"</SPAN>

rasters <SPAN class="operator token">=</SPAN> <SPAN class="punctuation token">[</SPAN><SPAN class="punctuation token">]</SPAN>
walk <SPAN class="operator token">=</SPAN> arcpy<SPAN class="punctuation token">.</SPAN>da<SPAN class="punctuation token">.</SPAN>Walk<SPAN class="punctuation token">(</SPAN>workspace<SPAN class="punctuation token">,</SPAN> datatype<SPAN class="operator token">=</SPAN><SPAN class="string token">"RasterDataset"</SPAN><SPAN class="punctuation token">,</SPAN> type<SPAN class="operator token">=</SPAN><SPAN class="string token">"TIF"</SPAN><SPAN class="punctuation token">)</SPAN>
<SPAN class="keyword token">for</SPAN> dirpath<SPAN class="punctuation token">,</SPAN> dirnames<SPAN class="punctuation token">,</SPAN> filenames <SPAN class="keyword token">in</SPAN> walk<SPAN class="punctuation token">:</SPAN>
&nbsp;&nbsp;&nbsp; <SPAN class="keyword token">for</SPAN> filename <SPAN class="keyword token">in</SPAN> filenames<SPAN class="punctuation token">:</SPAN>
&nbsp;&nbsp;&nbsp; rasters<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>dirpath<SPAN class="punctuation token">,</SPAN> filename<SPAN class="punctuation token">)</SPAN><SPAN class="punctuation token">)</SPAN>

<SPAN class="comment token"># calculate mean</SPAN>
sum <SPAN class="operator token">=</SPAN> CellStatistics<SPAN class="punctuation token">(</SPAN>rasters<SPAN class="punctuation token">,</SPAN> <SPAN class="string token">"MEAN"</SPAN><SPAN class="punctuation token">)</SPAN>
sum<SPAN class="punctuation token">.</SPAN>save<SPAN class="punctuation token">(</SPAN>os<SPAN class="punctuation token">.</SPAN>path<SPAN class="punctuation token">.</SPAN>join<SPAN class="punctuation token">(</SPAN>out_workspace<SPAN class="punctuation token">,</SPAN> <SPAN class="string token">"npp_sum.tif"</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>

Reference_RASTER.rar

curtvprice

I would look into arcpy.da.walk() to generate the list of rasters, and use the Spatial Analyst Cell Statistics tool on the raster list.

Walk—Help | ArcGIS for Desktop

Cell Statistics—Help | ArcGIS for Desktop

Cell Statistics is much more efficient, especially with long lists of rasters -- and it also has a very handy DATA/NODATA option so null is treated as zero in the cell by cell summation.

Untested code:

<SPAN class="keyword token">import</SPAN> arcpy
arcpy<SPAN class="punctuation token">.</SPAN>CheckOutExtension<SPAN class="punctuation token">(</SPAN><SPAN class="string token">"spatial"</SPAN><SPAN class="punctuation token">)</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>

<SPAN class="comment token"># generate list of rasters</SPAN>
workspace <SPAN class="operator token">=</SPAN> r<SPAN class="string token">"D:\"</SPAN>
rasters <SPAN class="operator token">=</SPAN> <SPAN class="punctuation token">[</SPAN><SPAN class="punctuation token">]</SPAN>
walk <SPAN class="operator token">=</SPAN> arcpy<SPAN class="punctuation token">.</SPAN>da<SPAN class="punctuation token">.</SPAN>Walk<SPAN class="punctuation token">(</SPAN>workspace<SPAN class="punctuation token">,</SPAN> datatype<SPAN class="operator token">=</SPAN><SPAN class="string token">"RasterDataset"</SPAN><SPAN class="punctuation token">,</SPAN> type<SPAN class="operator token">=</SPAN><SPAN class="string token">"TIF"</SPAN><SPAN class="punctuation token">)</SPAN>
<SPAN class="keyword token">for</SPAN> dirpath<SPAN class="punctuation token">,</SPAN> dirnames<SPAN class="punctuation token">,</SPAN> filenames <SPAN class="keyword token">in</SPAN> walk<SPAN class="punctuation token">:</SPAN>
&nbsp;&nbsp;&nbsp; <SPAN class="keyword token">for</SPAN> filename <SPAN class="keyword token">in</SPAN> filenames<SPAN class="punctuation token">:</SPAN>
&nbsp;&nbsp;&nbsp; rasters<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>dirpath<SPAN class="punctuation token">,</SPAN> filename<SPAN class="punctuation token">)</SPAN><SPAN class="punctuation token">)</SPAN>

<SPAN class="comment token"># calculate mean</SPAN>
sum <SPAN class="operator token">=</SPAN> CellStatistics<SPAN class="punctuation token">(</SPAN>rasters<SPAN class="punctuation token">,</SPAN> <SPAN class="string token">"MEAN"</SPAN><SPAN class="punctuation token">)</SPAN>
sum<SPAN class="punctuation token">.</SPAN>save<SPAN class="punctuation token">(</SPAN>os<SPAN class="punctuation token">.</SPAN>path<SPAN class="punctuation token">.</SPAN>join<SPAN class="punctuation token">(</SPAN>out_workspace<SPAN class="punctuation token">,</SPAN> <SPAN class="string token">"npp_sum.tif"</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>

Únete a la Conversación

Los miembros registrados pueden publicar, seguir actualizaciones y más. ¿Nuevo aquí? Regístrate gratis.

Enlaces Rápidos

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