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

Python script to calculate NDVI anomalies.

LuizVianna

Why the last line of this code is not working? It is the line of "Minus_3d" calculation?

This is a python script to calculate NDVI anomalies.

Thanks.

# Import arcpy module
import arcpy
from arcpy import env
from arcpy.sa import *
arcpy.env.overwriteOutput = True

# Check out any necessary licenses
arcpy.CheckOutExtension("spatial")

#Define workspaces
arcpy.env.workspace="D:\\Modis250\\NDVI_10dias_norm" #Workspace with decendial (ten days period) NDVI
#Filenames standard is "norm_ndvi_yyyy_dec_d.tif" where "yyyy" is the year and "d" is the ten days period: d=1 = january 1-10; d=2 = january 11-20; ...

for n in range (1,37):
d=str(1)
y = str(2004)
NDVINormAll=arcpy.ListRasters("norm_ndvi_*_"+d+".tif", "tif")#List of decendial NDVI for all years to calculate the dcendial mean
NDVINormY=arcpy.ListRasters("norm_ndvi_"+y+"_dec_"+d+".tif", "tif")#List of decendial NDVI rasters for the last year (y)
NDVIMeanOut = "D:\\Modis250\\NDVI_10dias_mean\\"+"mean_dec_"+d+".tif" #Define the decendial NDVI mean output (could be temporary)
NDVIMean = arcpy.gp.CellStatistics_sa(NDVINormAll, NDVIMeanOut,"MEAN", "NODATA")#Calculate the decendial NDVI mean (could be temporary)
NDVIAnomOut = "D:\\Modis250\\NDVI_10dias_anom\\"+"anom_"+y+"_dec_"+d+".tif" #Define the last year decendial NDVI anomaly output
NDVIAnom=arcpy.gp.Minus_3d(NDVINormY, NDVIMean, NDVIAnomOut) #Calculate the last year decendial NDVI anomaly and record the output tifs
print ("Decade "+d+" of "+y+" ok")

Batch_calculate_decadal_anomaly.py.zip

Find more posts tagged with

Accepted answers

XanderBakker

I tried it too and had some different results in Pro (including ArcGIS Pro closing very efficiently). I also obtained an error with the Cell statistics in some cases and decided to remove the minus and include the paths in the list of rasters that are used in the Cell Statistics. This seems to be working (you will need to change the paths back again).

<SPAN class="keyword token">def</SPAN> <SPAN class="token function">main</SPAN><SPAN class="punctuation token">(</SPAN><SPAN class="punctuation token">)</SPAN><SPAN class="punctuation token">:</SPAN>
&nbsp;&nbsp;&nbsp; <SPAN class="keyword token">import</SPAN> arcpy
&nbsp;&nbsp;&nbsp; <SPAN class="keyword token">import</SPAN> os
&nbsp;&nbsp;&nbsp; arcpy<SPAN class="punctuation token">.</SPAN>env<SPAN class="punctuation token">.</SPAN>overwriteOutput <SPAN class="operator token">=</SPAN> <SPAN class="token boolean">True</SPAN>

&nbsp;&nbsp;&nbsp; <SPAN class="comment token"># Check out SA license</SPAN>
&nbsp;&nbsp;&nbsp; arcpy<SPAN class="punctuation token">.</SPAN>CheckOutExtension<SPAN class="punctuation token">(</SPAN><SPAN class="string token">"spatial"</SPAN><SPAN class="punctuation token">)</SPAN>

&nbsp;&nbsp;&nbsp; <SPAN class="comment token"># settings</SPAN>
<SPAN class="comment token">##&nbsp;&nbsp;&nbsp; ws_norm = r'D:\Modis250\NDVI_10dias_norm'</SPAN>
<SPAN class="comment token">##&nbsp;&nbsp;&nbsp; ws_mean = r'D:\Modis250\NDVI_10dias_mean'</SPAN>
<SPAN class="comment token">##&nbsp;&nbsp;&nbsp; ws_anom = r'D:\Modis250\NDVI_10dias_anom'</SPAN>
&nbsp;&nbsp;&nbsp; ws_norm <SPAN class="operator token">=</SPAN> r<SPAN class="string token">'C:\GeoNet\Modis250\NDVI_10dias_norm'</SPAN>
&nbsp;&nbsp;&nbsp; ws_mean <SPAN class="operator token">=</SPAN> r<SPAN class="string token">'C:\GeoNet\Modis250\NDVI_10dias_mean'</SPAN>
&nbsp;&nbsp;&nbsp; ws_anom <SPAN class="operator token">=</SPAN> r<SPAN class="string token">'C:\GeoNet\Modis250\NDVI_10dias_anom'</SPAN>


&nbsp;&nbsp;&nbsp; <SPAN class="comment token"># create a list of all raster in norm folder</SPAN>
&nbsp;&nbsp;&nbsp; arcpy<SPAN class="punctuation token">.</SPAN>env<SPAN class="punctuation token">.</SPAN>workspace <SPAN class="operator token">=</SPAN> ws_norm
&nbsp;&nbsp;&nbsp; lst_ndvi_ras <SPAN class="operator token">=</SPAN> arcpy<SPAN class="punctuation token">.</SPAN>ListRasters<SPAN class="punctuation token">(</SPAN><SPAN class="punctuation token">)</SPAN>
&nbsp;&nbsp;&nbsp; <SPAN class="keyword token">print</SPAN> <SPAN class="punctuation token">(</SPAN><SPAN class="string token">"lst_ndvi_ras: {}"</SPAN><SPAN class="punctuation token">.</SPAN>format<SPAN class="punctuation token">(</SPAN>lst_ndvi_ras<SPAN class="punctuation token">)</SPAN><SPAN class="punctuation token">)</SPAN>

&nbsp;&nbsp;&nbsp; <SPAN class="comment token"># get list of years</SPAN>
&nbsp;&nbsp;&nbsp; decades <SPAN class="operator token">=</SPAN> GetListOfDecades<SPAN class="punctuation token">(</SPAN>lst_ndvi_ras<SPAN class="punctuation token">)</SPAN>
&nbsp;&nbsp;&nbsp; <SPAN class="keyword token">print</SPAN> <SPAN class="punctuation token">(</SPAN><SPAN class="string token">"decades: {}"</SPAN><SPAN class="punctuation token">.</SPAN>format<SPAN class="punctuation token">(</SPAN>decades<SPAN class="punctuation token">)</SPAN><SPAN class="punctuation token">)</SPAN>

&nbsp;&nbsp;&nbsp; <SPAN class="comment token"># loop through each year</SPAN>
&nbsp;&nbsp;&nbsp; <SPAN class="keyword token">for</SPAN> decade <SPAN class="keyword token">in</SPAN> decades<SPAN class="punctuation token">:</SPAN>
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; <SPAN class="comment token"># Get rasters for given year</SPAN>
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; lst_ndvi_decade <SPAN class="operator token">=</SPAN> GetListOfRasterForGivenDecade<SPAN class="punctuation token">(</SPAN>lst_ndvi_ras<SPAN class="punctuation token">,</SPAN> decade<SPAN class="punctuation token">)</SPAN>
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; <SPAN class="keyword token">print</SPAN><SPAN class="punctuation token">(</SPAN><SPAN class="string token">"lst_ndvi_decade: {}"</SPAN><SPAN class="punctuation token">.</SPAN>format<SPAN class="punctuation token">(</SPAN>lst_ndvi_decade<SPAN class="punctuation token">)</SPAN><SPAN class="punctuation token">)</SPAN>

&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; <SPAN class="comment token"># calculate mean raster for decade</SPAN>
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; lst_ndvi_decade_paths <SPAN class="operator token">=</SPAN> <SPAN class="punctuation token">[</SPAN>os<SPAN class="punctuation token">.</SPAN>path<SPAN class="punctuation token">.</SPAN>join<SPAN class="punctuation token">(</SPAN>ws_norm<SPAN class="punctuation token">,</SPAN> r<SPAN class="punctuation token">)</SPAN> <SPAN class="keyword token">for</SPAN> r <SPAN class="keyword token">in</SPAN> lst_ndvi_decade<SPAN class="punctuation token">]</SPAN>
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; <SPAN class="keyword token">print</SPAN><SPAN class="punctuation token">(</SPAN><SPAN class="string token">"lst_ndvi_decade_paths: {}"</SPAN><SPAN class="punctuation token">.</SPAN>format<SPAN class="punctuation token">(</SPAN>lst_ndvi_decade_paths<SPAN class="punctuation token">)</SPAN><SPAN class="punctuation token">)</SPAN>
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; mean_ras <SPAN class="operator token">=</SPAN> arcpy<SPAN class="punctuation token">.</SPAN>sa<SPAN class="punctuation token">.</SPAN>CellStatistics<SPAN class="punctuation token">(</SPAN>lst_ndvi_decade_paths<SPAN class="punctuation token">,</SPAN> <SPAN class="string token">"MEAN"</SPAN><SPAN class="punctuation token">,</SPAN> <SPAN class="string token">"DATA"</SPAN><SPAN class="punctuation token">)</SPAN>

&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; <SPAN class="comment token"># store mean ras for decade</SPAN>
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; <SPAN class="comment token"># mean_dec_d_.tif</SPAN>
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; out_name_mean <SPAN class="operator token">=</SPAN> <SPAN class="string token">"mean_dec_{}.tif"</SPAN><SPAN class="punctuation token">.</SPAN>format<SPAN class="punctuation token">(</SPAN>decade<SPAN class="punctuation token">)</SPAN>&nbsp; <SPAN class="comment token">###</SPAN>
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; out_name_path_mean <SPAN class="operator token">=</SPAN>&nbsp; os<SPAN class="punctuation token">.</SPAN>path<SPAN class="punctuation token">.</SPAN>join<SPAN class="punctuation token">(</SPAN>ws_mean<SPAN class="punctuation token">,</SPAN> out_name_mean<SPAN class="punctuation token">)</SPAN>
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; <SPAN class="keyword token">print</SPAN> <SPAN class="punctuation token">(</SPAN><SPAN class="string token">"out_name_path_mean: {}"</SPAN><SPAN class="punctuation token">.</SPAN>format<SPAN class="punctuation token">(</SPAN>out_name_path_mean<SPAN class="punctuation token">)</SPAN><SPAN class="punctuation token">)</SPAN>
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; mean_ras<SPAN class="punctuation token">.</SPAN>save<SPAN class="punctuation token">(</SPAN>out_name_path_mean<SPAN class="punctuation token">)</SPAN>

&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; <SPAN class="comment token"># loop through each raster in decade</SPAN>
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; <SPAN class="keyword token">for</SPAN> ndvi_ras <SPAN class="keyword token">in</SPAN> lst_ndvi_decade<SPAN class="punctuation token">:</SPAN>
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; <SPAN class="comment token"># calculate anom for each raster in decade</SPAN>
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; <SPAN class="comment token"># anom_2002_dec_1.tif = norm_ndvi_2002_dec_1.tif - mean_dec_1_.tif;</SPAN>
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; <SPAN class="keyword token">print</SPAN> <SPAN class="punctuation token">(</SPAN><SPAN class="string token">"ndvi_ras: {}"</SPAN><SPAN class="punctuation token">.</SPAN>format<SPAN class="punctuation token">(</SPAN>ndvi_ras<SPAN class="punctuation token">)</SPAN><SPAN class="punctuation token">)</SPAN>
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; <SPAN class="comment token"># minus_ras = arcpy.sa.Minus(ndvi_ras, out_name_path_mean)</SPAN>
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; out_name_path_norm <SPAN class="operator token">=</SPAN> os<SPAN class="punctuation token">.</SPAN>path<SPAN class="punctuation token">.</SPAN>join<SPAN class="punctuation token">(</SPAN>ws_norm<SPAN class="punctuation token">,</SPAN> ndvi_ras<SPAN class="punctuation token">)</SPAN>
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; minus_ras <SPAN class="operator token">=</SPAN> arcpy<SPAN class="punctuation token">.</SPAN>Raster<SPAN class="punctuation token">(</SPAN>out_name_path_mean<SPAN class="punctuation token">)</SPAN> <SPAN class="operator token">-</SPAN> arcpy<SPAN class="punctuation token">.</SPAN>Raster<SPAN class="punctuation token">(</SPAN>out_name_path_norm<SPAN class="punctuation token">)</SPAN>

&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; <SPAN class="comment token"># store anom raster</SPAN>
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; out_name_anom <SPAN class="operator token">=</SPAN> ndvi_ras<SPAN class="punctuation token">.</SPAN>replace<SPAN class="punctuation token">(</SPAN><SPAN class="string token">"norm_ndvi"</SPAN><SPAN class="punctuation token">,</SPAN> <SPAN class="string token">"anom"</SPAN><SPAN class="punctuation token">)</SPAN>
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; out_name_path_anom <SPAN class="operator token">=</SPAN> os<SPAN class="punctuation token">.</SPAN>path<SPAN class="punctuation token">.</SPAN>join<SPAN class="punctuation token">(</SPAN>ws_anom<SPAN class="punctuation token">,</SPAN> out_name_anom<SPAN class="punctuation token">)</SPAN>
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; minus_ras<SPAN class="punctuation token">.</SPAN>save<SPAN class="punctuation token">(</SPAN>out_name_path_anom<SPAN class="punctuation token">)</SPAN>
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; <SPAN class="keyword token">print</SPAN> <SPAN class="punctuation token">(</SPAN><SPAN class="string token">"out_name_path_anom: {}"</SPAN><SPAN class="punctuation token">.</SPAN>format<SPAN class="punctuation token">(</SPAN>out_name_path_anom<SPAN class="punctuation token">)</SPAN><SPAN class="punctuation token">)</SPAN>


&nbsp;&nbsp;&nbsp; <SPAN class="comment token"># return SA license</SPAN>
&nbsp;&nbsp;&nbsp; arcpy<SPAN class="punctuation token">.</SPAN>CheckInExtension<SPAN class="punctuation token">(</SPAN><SPAN class="string token">"spatial"</SPAN><SPAN class="punctuation token">)</SPAN>

<SPAN class="keyword token">def</SPAN> <SPAN class="token function">GetListOfDecades</SPAN><SPAN class="punctuation token">(</SPAN>lst<SPAN class="punctuation token">)</SPAN><SPAN class="punctuation token">:</SPAN>
&nbsp;&nbsp;&nbsp; <SPAN class="comment token"># norm_ndvi_yyyy_dec_d.tif</SPAN>
&nbsp;&nbsp;&nbsp; lst1 <SPAN class="operator token">=</SPAN> list<SPAN class="punctuation token">(</SPAN>set<SPAN class="punctuation token">(</SPAN><SPAN class="punctuation token">[</SPAN>r<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">4</SPAN><SPAN class="punctuation token">]</SPAN> <SPAN class="keyword token">for</SPAN> r <SPAN class="keyword token">in</SPAN> lst<SPAN class="punctuation token">]</SPAN><SPAN class="punctuation token">)</SPAN><SPAN class="punctuation token">)</SPAN>
&nbsp;&nbsp;&nbsp; <SPAN class="keyword token">return</SPAN> <SPAN class="punctuation token">[</SPAN>r<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> <SPAN class="keyword token">for</SPAN> r <SPAN class="keyword token">in</SPAN> lst1<SPAN class="punctuation token">]</SPAN>

<SPAN class="keyword token">def</SPAN> <SPAN class="token function">GetListOfRasterForGivenDecade</SPAN><SPAN class="punctuation token">(</SPAN>lst<SPAN class="punctuation token">,</SPAN> decade<SPAN class="punctuation token">)</SPAN><SPAN class="punctuation token">:</SPAN>
&nbsp;&nbsp;&nbsp; lst1 <SPAN class="operator token">=</SPAN> <SPAN class="punctuation token">[</SPAN><SPAN class="punctuation token">]</SPAN>
&nbsp;&nbsp;&nbsp; search <SPAN class="operator token">=</SPAN> <SPAN class="string token">"_{}.tif"</SPAN><SPAN class="punctuation token">.</SPAN>format<SPAN class="punctuation token">(</SPAN>decade<SPAN class="punctuation token">)</SPAN>
&nbsp;&nbsp;&nbsp; <SPAN class="keyword token">for</SPAN> r <SPAN class="keyword token">in</SPAN> lst<SPAN class="punctuation token">:</SPAN>
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; <SPAN class="keyword token">if</SPAN> search <SPAN class="keyword token">in</SPAN> r<SPAN class="punctuation token">.</SPAN>lower<SPAN class="punctuation token">(</SPAN><SPAN class="punctuation token">)</SPAN><SPAN class="punctuation token">:</SPAN>
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; lst1<SPAN class="punctuation token">.</SPAN>append<SPAN class="punctuation token">(</SPAN>r<SPAN class="punctuation token">)</SPAN>
&nbsp;&nbsp;&nbsp; <SPAN class="keyword token">return</SPAN> lst1


<SPAN class="keyword token">if</SPAN> __name__ <SPAN class="operator token">==</SPAN> <SPAN class="string token">'__main__'</SPAN><SPAN class="punctuation token">:</SPAN>
&nbsp;&nbsp;&nbsp; main<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>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN></SPAN>

All comments

XanderBakker

Something like this:

<SPAN class="keyword token">def</SPAN> <SPAN class="token function">GetListOfDecades</SPAN><SPAN class="punctuation token">(</SPAN>lst<SPAN class="punctuation token">)</SPAN><SPAN class="punctuation token">:</SPAN>
&nbsp;&nbsp;&nbsp; <SPAN class="comment token"># norm_ndvi_yyyy_dec_d.tif</SPAN>
&nbsp;&nbsp;&nbsp; lst1 <SPAN class="operator token">=</SPAN> list<SPAN class="punctuation token">(</SPAN>set<SPAN class="punctuation token">(</SPAN><SPAN class="punctuation token">[</SPAN>r<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">4</SPAN><SPAN class="punctuation token">]</SPAN> <SPAN class="keyword token">for</SPAN> r <SPAN class="keyword token">in</SPAN> lst<SPAN class="punctuation token">]</SPAN><SPAN class="punctuation token">)</SPAN><SPAN class="punctuation token">)</SPAN>
&nbsp;&nbsp;&nbsp; <SPAN class="keyword token">return</SPAN> <SPAN class="punctuation token">[</SPAN>r<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> <SPAN class="keyword token">for</SPAN> r <SPAN class="keyword token">in</SPAN> lst1<SPAN class="punctuation token">]</SPAN>

<SPAN class="keyword token">def</SPAN> <SPAN class="token function">GetListOfRasterForGivenDecade</SPAN><SPAN class="punctuation token">(</SPAN>lst<SPAN class="punctuation token">,</SPAN> decade<SPAN class="punctuation token">)</SPAN><SPAN class="punctuation token">:</SPAN>
&nbsp;&nbsp;&nbsp; lst1 <SPAN class="operator token">=</SPAN> <SPAN class="punctuation token">[</SPAN><SPAN class="punctuation token">]</SPAN>
&nbsp;&nbsp;&nbsp; search <SPAN class="operator token">=</SPAN> <SPAN class="string token">"_{}.tif"</SPAN><SPAN class="punctuation token">.</SPAN>format<SPAN class="punctuation token">(</SPAN>decade<SPAN class="punctuation token">)</SPAN>
&nbsp;&nbsp;&nbsp; <SPAN class="keyword token">for</SPAN> r <SPAN class="keyword token">in</SPAN> lst<SPAN class="punctuation token">:</SPAN>
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; <SPAN class="keyword token">if</SPAN> search <SPAN class="keyword token">in</SPAN> r<SPAN class="punctuation token">.</SPAN>lower<SPAN class="punctuation token">(</SPAN><SPAN class="punctuation token">)</SPAN><SPAN class="punctuation token">:</SPAN>
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; lst1<SPAN class="punctuation token">.</SPAN>append<SPAN class="punctuation token">(</SPAN>r<SPAN class="punctuation token">)</SPAN>
&nbsp;&nbsp;&nbsp; <SPAN class="keyword token">return</SPAN> lst1

<SPAN class="keyword token">import</SPAN> arcpy
<SPAN class="keyword token">import</SPAN> os
arcpy<SPAN class="punctuation token">.</SPAN>env<SPAN class="punctuation token">.</SPAN>overwriteOutput <SPAN class="operator token">=</SPAN> <SPAN class="token boolean">True</SPAN>

<SPAN class="comment token"># Check out SA license</SPAN>
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"># settings</SPAN>
ws_norm <SPAN class="operator token">=</SPAN> r<SPAN class="string token">'D:\Modis250\NDVI_10dias_norm'</SPAN>
ws_mean <SPAN class="operator token">=</SPAN> r<SPAN class="string token">'D:\Modis250\NDVI_10dias_mean'</SPAN>
ws_anom <SPAN class="operator token">=</SPAN> r<SPAN class="string token">'D:\Modis250\NDVI_10dias_anom'</SPAN>

<SPAN class="comment token"># create a list of all raster in norm folder</SPAN>
arcpy<SPAN class="punctuation token">.</SPAN>env<SPAN class="punctuation token">.</SPAN>workspace <SPAN class="operator token">=</SPAN> ws_norm
lst_ndvi_ras <SPAN class="operator token">=</SPAN> arcpy<SPAN class="punctuation token">.</SPAN>ListRasters<SPAN class="punctuation token">(</SPAN><SPAN class="punctuation token">)</SPAN>

<SPAN class="comment token"># get list of years</SPAN>
decades <SPAN class="operator token">=</SPAN> GetListOfDecades<SPAN class="punctuation token">(</SPAN>lst_ndvi_ras<SPAN class="punctuation token">)</SPAN>

<SPAN class="comment token"># loop through each year</SPAN>
<SPAN class="keyword token">for</SPAN> decade <SPAN class="keyword token">in</SPAN> decades<SPAN class="punctuation token">:</SPAN>
&nbsp;&nbsp;&nbsp; <SPAN class="comment token"># Get rasters for given year</SPAN>
&nbsp;&nbsp;&nbsp; lst_ndvi_decade <SPAN class="operator token">=</SPAN> GetListOfRasterForGivenDecade<SPAN class="punctuation token">(</SPAN>lst_ndvi_ras<SPAN class="punctuation token">,</SPAN> decade<SPAN class="punctuation token">)</SPAN>

&nbsp;&nbsp;&nbsp; <SPAN class="comment token"># calculate mean raster for decade</SPAN>
&nbsp;&nbsp;&nbsp; mean_ras <SPAN class="operator token">=</SPAN> arcpy<SPAN class="punctuation token">.</SPAN>sa<SPAN class="punctuation token">.</SPAN>CellStatistics<SPAN class="punctuation token">(</SPAN>lst_ndvi_year<SPAN class="punctuation token">,</SPAN> <SPAN class="string token">"MEAN"</SPAN><SPAN class="punctuation token">,</SPAN> <SPAN class="string token">"DATA"</SPAN><SPAN class="punctuation token">)</SPAN>

&nbsp;&nbsp;&nbsp; <SPAN class="comment token"># store mean ras for decade</SPAN>
&nbsp;&nbsp;&nbsp; <SPAN class="comment token"># mean_dec_d_.tif</SPAN>
&nbsp;&nbsp;&nbsp; out_name_mean <SPAN class="operator token">=</SPAN> <SPAN class="string token">"mean_dec_{}_.tif"</SPAN><SPAN class="punctuation token">.</SPAN>format<SPAN class="punctuation token">(</SPAN>decade<SPAN class="punctuation token">)</SPAN>
&nbsp;&nbsp;&nbsp; out_name_path_mean <SPAN class="operator token">=</SPAN>&nbsp; os<SPAN class="punctuation token">.</SPAN>path<SPAN class="punctuation token">.</SPAN>join<SPAN class="punctuation token">(</SPAN>ws_mean<SPAN class="punctuation token">,</SPAN> out_name_mean<SPAN class="punctuation token">)</SPAN>
&nbsp;&nbsp;&nbsp; mean_ras<SPAN class="punctuation token">.</SPAN>save<SPAN class="punctuation token">(</SPAN>out_name_path_mean<SPAN class="punctuation token">)</SPAN>

&nbsp;&nbsp;&nbsp; <SPAN class="comment token"># loop through each raster in decade</SPAN>
&nbsp;&nbsp;&nbsp; <SPAN class="keyword token">for</SPAN> ndvi_ras <SPAN class="keyword token">in</SPAN> lst_ndvi_ras<SPAN class="punctuation token">:</SPAN>
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; <SPAN class="comment token"># calculate anom for each raster in decade</SPAN>
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; <SPAN class="comment token"># anom_2002_dec_1.tif = norm_ndvi_2002_dec_1.tif - mean_dec_1_.tif;</SPAN>
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; minus_ras <SPAN class="operator token">=</SPAN> arcpy<SPAN class="punctuation token">.</SPAN>sa<SPAN class="punctuation token">.</SPAN>Minus<SPAN class="punctuation token">(</SPAN>ndvi_ras<SPAN class="punctuation token">,</SPAN> out_name_path_mean<SPAN class="punctuation token">)</SPAN>

&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; <SPAN class="comment token"># store anom raster</SPAN>
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; out_name_anom <SPAN class="operator token">=</SPAN> ndvi_ras<SPAN class="punctuation token">.</SPAN>replace<SPAN class="punctuation token">(</SPAN><SPAN class="string token">"norm_ndvi"</SPAN><SPAN class="punctuation token">,</SPAN> <SPAN class="string token">"anom"</SPAN><SPAN class="punctuation token">)</SPAN>
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; out_name_path_anom <SPAN class="operator token">=</SPAN> os<SPAN class="punctuation token">.</SPAN>path<SPAN class="punctuation token">.</SPAN>join<SPAN class="punctuation token">(</SPAN>ws_anom<SPAN class="punctuation token">,</SPAN> out_name_anom<SPAN class="punctuation token">)</SPAN>
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; minus_ras<SPAN class="punctuation token">.</SPAN>save<SPAN class="punctuation token">(</SPAN>out_name_path_anom<SPAN class="punctuation token">)</SPAN>


<SPAN class="comment token"># return SA license</SPAN>
arcpy<SPAN class="punctuation token">.</SPAN>CheckInExtension<SPAN class="punctuation token">(</SPAN><SPAN class="string token">"spatial"</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>

... and/or include some print statements to see what is going on.

XanderBakker

You shouldn't run this by parts, since the code is within the "def main()" and "if __name__ == '__main__':" the variable will not be set.

You could try to run this by parts, without the outer parts

LuizVianna

Thank you, Xander!

Unfortunatly I did´n understood this code, so I´ve tried to run the script step by step. However, I dont know how. I copied and paste fom line 1 to 16, executed the code and then I typed "ws_mean", the error was:

Traceback (most recent call last):
File "<string>", line 1, in <module>
NameError: name 'ws_mean' is not defined

I expected for 'D:\Modis250\NDVI_10dias_mean'

If I copy all the script and run, I had this error:

Traceback (most recent call last):
File "<string>", line 65, in <module>
File "<string>", line 27, in main
NameError: name 'lst_ndvi_year' is not defined

I´m working on ArcGis Pro 2.1.2, an I think my Python is 2.7 (I also have ArcGis 10.5).

LuizVianna

The 36 "decades" NDVI means could be temporary... i dont need to store them, just use to calculate the anomaly...

LuizVianna

Hi Xander!

Let me try to explain what I need: There are 423 NDVI tifs in folder 'D: \ Modis250 \ NDVI_10dias_norm'. The tifs are named by year (from 2002 to 2018) and "decade" (36 ten days periods) as: norm_ndvi_yyyy_dec_d.tif, where "yyyy" is the year and the "d" is the "decade".

The loop calculates the mean for each decade, considering all years:

Mean (norm_ndvi_2002_dec_1.tif, norm_ndvi_2003_dec_1.tif, ..., norm_ndvi_2018_dec_1.tif ).

The 36 means are storaged in 'D:\Modis250\NDVI_10dias_mean' like 'mean_dec_d_.tif'.

The loop also calculates the anomaly for a choosing year ("y" in the code), subtracting each decade NDVI from those year by the NDVI mean for each decade:

anom_2002_dec_1.tif = norm_ndvi_2002_dec_1.tif - mean_dec_1_.tif;

anom_2002_dec_2.tif = norm_ndvi_2002_dec_2.tif - mean_dec_2_.tif;

anom_yyyy_dec_d.tif = norm_ndvi_yyyy_dec_d.tif - mean_dec_d_.tif;

The anomaly results shoud be storages in 'D:\Modis250\NDVI_10dias_anom' like anom_yyyy_dec_d.tif

The code seems to be ok until the mean calculation, however the 'Minus' operation is not working.

How can I set the input rasters for the "Minus" operation?

Thank you!

XanderBakker

I can imagine that the idea is to calculate for each raster in the list with "decendial NDVI rasters for the last year (y)" the difference needs to be calculated with the mean. The mean should be calculated outside the look and the loop should go through each item in the list of rasters.

Look at the following structure:

year_range <SPAN class="operator token">=</SPAN> range<SPAN class="punctuation token">(</SPAN><SPAN class="number token">2004</SPAN><SPAN class="punctuation token">,</SPAN> <SPAN class="number token">2011</SPAN><SPAN class="punctuation token">)</SPAN> <SPAN class="comment token"># to get the values from 2004 until 2010</SPAN>
<SPAN class="keyword token">for</SPAN> year <SPAN class="keyword token">in</SPAN> year_range<SPAN class="punctuation token">:</SPAN>
&nbsp;&nbsp;&nbsp; <SPAN class="keyword token">print</SPAN> year

&nbsp;&nbsp;&nbsp; <SPAN class="comment token"># for each year:</SPAN>
&nbsp;&nbsp;&nbsp; <SPAN class="comment token"># - create list of all NDVI for current year: lst_ndvi_ras</SPAN>
&nbsp;&nbsp;&nbsp; <SPAN class="comment token"># - calculate the mean of all NDVI for current year: mean_ndvi_year</SPAN>

&nbsp;&nbsp;&nbsp; <SPAN class="keyword token">for</SPAN> ndvi_ras <SPAN class="keyword token">in</SPAN> lst_ndvi_ras<SPAN class="punctuation token">:</SPAN>
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; <SPAN class="comment token"># calculate anomaly and save result with unique name</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>

I assume you will have a range of years. Each year will have a number of NDVI rasters (is that true?). How are the NDVI rasters for each year stored? Each year has its own folder or the year is part of the name?

JoeBorgione

Good catch Xander! I've indented the entire loop.

Luiz- as Xander mentions, your loop has a malfunction in it as you don't mention the the incrementer (n) to move through the loop.

Xander Bakker‌

Luiz Vianna‌

<SPAN class="keyword token">for</SPAN> n <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">37</SPAN><SPAN class="punctuation token">)</SPAN><SPAN class="punctuation token">:</SPAN>
&nbsp; <SPAN class="keyword token">print</SPAN> <SPAN class="string token">'{}'</SPAN><SPAN class="punctuation token">.</SPAN>format<SPAN class="punctuation token">(</SPAN>n<SPAN class="punctuation token">)</SPAN><SPAN class="line-numbers-rows"><SPAN>‍</SPAN><SPAN>‍</SPAN></SPAN>

XanderBakker

Hi Joe Borgione ,

What does not make sense is:

There is for loop on line 15 that would require indentation afterwards.
During the loop a variable "n" is set that is not used, causing a loop of 37 iteration doing the same thing
the Minus_3d has the output set as parameter "Minus_3d (in_raster_or_constant1, in_raster_or_constant2, out_raster)" and this should not be used in combination as "another_output = Minus_3d (in_raster_or_constant1, in_raster_or_constant2, out_raster)"
The first parameter of the Minus_3d is a list of rasters and should be a single raster

LuizVianna

Thank you Joe!

I´m not expert in Python and I´m just copying some codes from internet and trying to make them work.

I changed lines 24 and line 30 to:

24 NDVIMeanOut = ("D:\\Modis250\\NDVI_10dias_anom\\mean_dec_{}.tif".format(d))

30 NDVIAnomOut = ("D:\\Modis250\\NDVI_10dias_anom\\anom_{}_dec_{}.tif".format(y, d))

and when I run the code, after line 33

33 NDVIAnom=arcpy.gp.Minus_3d(NDVINormY, NDVIMean, NDVIAnomOut)

the error is:

Traceback (most recent call last):
File "<string>", line 24, in <module>
File "c:\program files\arcgis\pro\Resources\arcpy\arcpy\geoprocessing\_base.py", line 506, in <lambda>
return lambda *args: val(*gp_fixargs(args, True))
RuntimeError: Object: Error in executing tool

If I try:

NDVIAnom = NDVINormY - NDVIMean
NDVIAnom.save(NDVIAnomOut)

the error is

Traceback (most recent call last):
File "<string>", line 25, in <module>
TypeError: unsupported operand type(s) for -: 'list' and 'geoprocessing server result object'

New code is: (sorry, I dont know how to format it)

# Import arcpy module
import arcpy
from arcpy import env
from arcpy.sa import *
arcpy.env.overwriteOutput = True

# Check out any necessary licenses
arcpy.CheckOutExtension("spatial")

for n in range (1,37):
d=str(n)
y = str(2004)
NDVINormAll=arcpy.ListRasters("norm_ndvi_*_"+d+".tif", "tif")#List of decendial NDVI for all years to calculate the dcendial mean
NDVINormY=arcpy.ListRasters("norm_ndvi_"+y+"_dec_"+d+".tif", "tif")#List of decendial NDVI rasters for the last year (y)
#NDVIMeanOut = "D:\\Modis250\\NDVI_10dias_mean\\"+"mean_dec_"+d+".tif" #Define the decendial NDVI mean output (could be temporary)
NDVIMeanOut = ("D:\\Modis250\\NDVI_10dias_anom\\mean_dec_{}.tif".format(d))
NDVIMean = arcpy.gp.CellStatistics_sa(NDVINormAll, NDVIMeanOut,"MEAN", "NODATA")#Calculate the decendial NDVI mean (could be temporary)
#NDVIAnomOut = "D:\\Modis250\\NDVI_10dias_anom\\"+"anom_"+y+"_dec_"+d+".tif" #Define the last year decendial NDVI anomaly output
NDVIAnomOut = ("D:\\Modis250\\NDVI_10dias_anom\\anom_{}_dec_{}.tif".format(y, d))
#NDVIAnom=arcpy.gp.Minus_3d(NDVINormY, NDVIMean, NDVIAnomOut) #Calculate the last year decendial NDVI anomaly and record the output tifs
NDVIAnom = NDVINormY - NDVIMean
NDVIAnom.save(NDVIAnomOut)
#print ("Decade "+d+" of "+y+" ok")
print("Decade {} of {} ok".format(d, y))

JoeBorgione

Luiz- I've reformatted your script as Dan suggested to make it a little easier to read; in this case you say your problem is on line 33. The error indicates your arguments are in error. Looking upstream, my guess is Line 30, and your definition of NDVIAnomOut. Try picking out the individual components that make of the string and set them to variables rather than stringing the whole thing inside quotes: that's what I do and then using Dan's approach of '{}'.format(variable_name) works quite nicely...

<SPAN class="comment token"># Import arcpy module</SPAN>
<SPAN class="keyword token">import</SPAN> arcpy
<SPAN class="keyword token">from</SPAN> arcpy <SPAN class="keyword token">import</SPAN> env
<SPAN class="keyword token">from</SPAN> arcpy<SPAN class="punctuation token">.</SPAN>sa <SPAN class="keyword token">import</SPAN> <SPAN class="operator token">*</SPAN> 
arcpy<SPAN class="punctuation token">.</SPAN>env<SPAN class="punctuation token">.</SPAN>overwriteOutput <SPAN class="operator token">=</SPAN> <SPAN class="token boolean">True</SPAN>

<SPAN class="comment token"># Check out any necessary licenses</SPAN>
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">#Define workspaces</SPAN>
arcpy<SPAN class="punctuation token">.</SPAN>env<SPAN class="punctuation token">.</SPAN>workspace<SPAN class="operator token">=</SPAN><SPAN class="string token">"D:\\Modis250\\NDVI_10dias_norm"</SPAN> <SPAN class="comment token">#Workspace with decendial (ten days period) NDVI</SPAN>

<SPAN class="comment token">#Filenames standard is "norm_ndvi_yyyy_dec_d.tif" where "yyyy" is the year and "d" is the ten days period: d=1 = january 1-10; d=2 = january 11-20; ...</SPAN>

<SPAN class="keyword token">for</SPAN> n <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">37</SPAN><SPAN class="punctuation token">)</SPAN><SPAN class="punctuation token">:</SPAN>
&nbsp;&nbsp;&nbsp; d<SPAN class="operator token">=</SPAN>str<SPAN class="punctuation token">(</SPAN><SPAN class="number token">1</SPAN><SPAN class="punctuation token">)</SPAN>
&nbsp;&nbsp;&nbsp; y <SPAN class="operator token">=</SPAN> str<SPAN class="punctuation token">(</SPAN><SPAN class="number token">2004</SPAN><SPAN class="punctuation token">)</SPAN>
&nbsp;&nbsp;&nbsp; NDVINormAll<SPAN class="operator token">=</SPAN>arcpy<SPAN class="punctuation token">.</SPAN>ListRasters<SPAN class="punctuation token">(</SPAN><SPAN class="string token">"norm_ndvi_*_"</SPAN><SPAN class="operator token">+</SPAN>d<SPAN class="operator token">+</SPAN><SPAN class="string token">".tif"</SPAN><SPAN class="punctuation token">,</SPAN> <SPAN class="string token">"tif"</SPAN><SPAN class="punctuation token">)</SPAN>
&nbsp;&nbsp;&nbsp; <SPAN class="comment token">#List of decendial NDVI for all years to calculate the dcendial mean</SPAN>

&nbsp;&nbsp;&nbsp; NDVINormY<SPAN class="operator token">=</SPAN>arcpy<SPAN class="punctuation token">.</SPAN>ListRasters<SPAN class="punctuation token">(</SPAN><SPAN class="string token">"norm_ndvi_"</SPAN><SPAN class="operator token">+</SPAN>y<SPAN class="operator token">+</SPAN><SPAN class="string token">"_dec_"</SPAN><SPAN class="operator token">+</SPAN>d<SPAN class="operator token">+</SPAN><SPAN class="string token">".tif"</SPAN><SPAN class="punctuation token">,</SPAN> <SPAN class="string token">"tif"</SPAN><SPAN class="punctuation token">)</SPAN>
&nbsp;&nbsp;&nbsp; <SPAN class="comment token">#List of decendial NDVI rasters for the last year (y)</SPAN>

&nbsp;&nbsp;&nbsp; NDVIMeanOut <SPAN class="operator token">=</SPAN> <SPAN class="string token">"D:\\Modis250\\NDVI_10dias_mean\\"</SPAN><SPAN class="operator token">+</SPAN><SPAN class="string token">"mean_dec_"</SPAN><SPAN class="operator token">+</SPAN>d<SPAN class="operator token">+</SPAN><SPAN class="string token">".tif"</SPAN> 
&nbsp;&nbsp;&nbsp; <SPAN class="comment token">#Define the decendial NDVI mean output (could be temporary)</SPAN>

&nbsp;&nbsp;&nbsp; NDVIMean <SPAN class="operator token">=</SPAN> arcpy<SPAN class="punctuation token">.</SPAN>gp<SPAN class="punctuation token">.</SPAN>CellStatistics_sa<SPAN class="punctuation token">(</SPAN>NDVINormAll<SPAN class="punctuation token">,</SPAN> NDVIMeanOut<SPAN class="punctuation token">,</SPAN><SPAN class="string token">"MEAN"</SPAN><SPAN class="punctuation token">,</SPAN> <SPAN class="string token">"NODATA"</SPAN><SPAN class="punctuation token">)</SPAN>
&nbsp;&nbsp;&nbsp; <SPAN class="comment token">#Calculate the decendial NDVI mean (could be temporary)</SPAN>

&nbsp;&nbsp;&nbsp; NDVIAnomOut <SPAN class="operator token">=</SPAN> <SPAN class="string token">"D:\\Modis250\\NDVI_10dias_anom\\"</SPAN><SPAN class="operator token">+</SPAN><SPAN class="string token">"anom_"</SPAN><SPAN class="operator token">+</SPAN>y<SPAN class="operator token">+</SPAN><SPAN class="string token">"_dec_"</SPAN><SPAN class="operator token">+</SPAN>d<SPAN class="operator token">+</SPAN><SPAN class="string token">".tif"</SPAN> 
&nbsp;&nbsp;&nbsp; Define the last year decendial NDVI anomaly output

&nbsp;&nbsp;&nbsp; NDVIAnom<SPAN class="operator token">=</SPAN>arcpy<SPAN class="punctuation token">.</SPAN>gp<SPAN class="punctuation token">.</SPAN>Minus_3d<SPAN class="punctuation token">(</SPAN>NDVINormY<SPAN class="punctuation token">,</SPAN> NDVIMean<SPAN class="punctuation token">,</SPAN> NDVIAnomOut<SPAN class="punctuation token">)</SPAN> 
&nbsp;&nbsp;&nbsp; <SPAN class="comment token">#Calculate the last year decendial NDVI anomaly and record the output tifs</SPAN>

&nbsp;&nbsp;&nbsp; <SPAN class="keyword token">print</SPAN> <SPAN class="punctuation token">(</SPAN><SPAN class="string token">"Decade "</SPAN><SPAN class="operator token">+</SPAN>d<SPAN class="operator token">+</SPAN><SPAN class="string token">" of "</SPAN><SPAN class="operator token">+</SPAN>y<SPAN class="operator token">+</SPAN><SPAN class="string token">" ok"</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>

LuizVianna

Thank you Dan!

The problem is in line "NDVIAnom=arcpy.gp.Minus_3d(NDVINormY, NDVIMean, NDVIAnomOut) #Calculate the last year decendial NDVI anomaly and record the output tifs"

I´ve also tried "NDVIAnom = NDVINormY - NDVIMean
NDVIAnom.save(NDVIAnomOut)"

The errors are:

For the arcpy.gp.Minus_3d approach:

Traceback (most recent call last):
File "<string>", line 22, in <module>
File "c:\program files\arcgis\pro\Resources\arcpy\arcpy\geoprocessing\_base.py", line 506, in <lambda>
return lambda *args: val(*gp_fixargs(args, True))
RuntimeError: Object: Error in executing tool

For the "NDVIAnom = NDVINormY - NDVIMean

NDVIAnom.save(NDVIAnomOut)" approach

Traceback (most recent call last):
File "<string>", line 23, in <module>
TypeError: unsupported operand type(s) for -: 'list' and 'geoprocessing server result object'

DanPatterson_Retired

Luiz ... could you format your code Code Formatting the basics please

And the error message?

The last line?

<SPAN class="comment token">#print ("Decade "+d+" of "+y+" ok")</SPAN>

<SPAN class="comment token"># ---- maybe try this instead, it handles type conversion to string</SPAN>

<SPAN class="keyword token">print</SPAN><SPAN class="punctuation token">(</SPAN><SPAN class="string token">"Decade {} of {} ok"</SPAN><SPAN class="punctuation token">.</SPAN>format<SPAN class="punctuation token">(</SPAN>d<SPAN class="punctuation token">,</SPAN> y<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>

Join the Conversation

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

Quick Links

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