Xander Bakker Python script to calculate NDVI anomalies https://geonet.esri.com/people/xander_bakker I´d like to automate all my workflow, from downloaded original raster datasets to anomaly calculator: 3 years data and folder structure can be downloaded here. Question 1: What is the best option: to run separated scripts or create a unique code? Question 2: How to translate the model builder to python script (Python script 2)? Question 3: If the best option is to run separated, how to make it in a logical sequence (script 1 -> script 2 -> script 3) Summary workflow: * Select and download NDVI filtered raster datasets;* Store downloaded original raster datasets at X:\Modis250\Originais\ * Organizing files and folders: (Python script 1) * Rename original raster datasets (MODIS.ndvi.YYYYjjj.yL6000.BOKU.tif) to filename standard “ndvi_YYYY_dec_dd.tif”, where “YYYY” is the year, “jjj” is the Julian day and “dd” is the decade (ten days period). * Transfer renamed files from X:\Modis250\Originais to X:\Modis250\NDVI10dias * Pre-processing raster datasets: (Python script 2)* Project renamed raster datasets ("ndvi_YYYY_dec_dd.tif") from Lat/Lon WGS 1984 to Lat/Lon SIRGAS2000 using "SAD_1969_To_WGS_1984_14 + SAD_1969_To_SIRGAS_2000_1" geographic transformation. * Overwrite spurious data values: NDVI = Con("%NDVI%" 9998,9999,"%NDVI%")) * Standardize raster datasets: NDVI = (NDVI – Mean(NDVI))/Std(NDVI) * Scale raster datasets to 0 – 100: NDVI = ((NDVI-Min(NDVI))/(Max(NDVI)-Min(NDVI)) * 100 * Mask standardized raster datasets by the area of interest (X:\Modis250\SRTM\mde_sc_90.tif) and output to X:\Modis250\NDVI_10dias_norm * Processing data (Python script 3)* Calculate NDVI mean for each decadal: mean_dec_dd.tif = MEAN (ndvi_YYY1_dec_dd.tif, ndvi_YYY2_dec_dd.tif, ndvi_YYY3_dec_dd.tif,…) * Calculate NDVI anomaly for each year by decade: anom_yyyy_dec_dd.tif = ndvi_YYYY_dec_dd.tif - mean_dec_dd.tif Python scripts 1 and 3 are working. The scrip 2 (Pre-processing) is a model in a toolbox. Python script 1: written by myself, maybe need to be improved! #Rename original files import arcpy , datetime , glob , os arcpy . env . overwriteOutput = True arcpy . env . workspace = r 'O:\Modis250\Originais' infolder = 'O:\Modis250\Originals' outfolder = 'O:\Modis250\NDVI_10dias' list1 = arcpy . ListRasters ( "*.tif" ) list1 . sort ( ) for raster in list1 : source_path = os . path . join ( infolder , raster ) oldFilename = raster dd = int ( oldFilename [ 15 : 18 ] ) if dd 11 : dd = 1 elif dd 21 : dd = 2 elif dd 31 : dd = 3 elif dd 41 : dd = 4 elif dd 51 : dd = 5 elif dd 61 : dd = 6 elif dd 71 : dd = 7 elif dd 81 : dd = 8 elif dd 91 : dd = 9 elif dd 101 : dd = 10 elif dd 111 : dd = 11 elif dd 121 : dd = 12 elif dd 131 : dd = 13 elif dd 141 : dd = 14 elif dd 151 : dd = 15 elif dd 161 : dd = 16 elif dd 171 : dd = 17 elif dd 181 : dd = 18 elif dd 191 : dd = 19 elif dd 201 : dd = 20 elif dd 211 : dd = 21 elif dd 221 : dd = 22 elif dd 231 : dd = 23 elif dd 241 : dd = 24 elif dd 251 : dd = 25 elif dd 261 : dd = 26 elif dd 271 : dd = 27 elif dd 281 : dd = 28 elif dd 291 : dd = 29 elif dd 301 : dd = 30 elif dd 311 : dd = 31 elif dd 321 : dd = 32 elif dd 331 : dd = 33 elif dd 341 : dd = 34 elif dd 351 : dd = 35 else : dd = 36 dd = str ( dd ) newFilename = "ndvi_" + oldFilename [ 11 : 15 ] + "_dec_" + dd + ".tif" destination_path = os . path . join ( outfolder , newFilename ) os . rename ( source_path , destination_path ) Model Builder -> Python script 2 Python script 3: Python script to calculate NDVI anomalies def main ( ) : import arcpy import os arcpy . env . overwriteOutput = True # Check out SA license arcpy . CheckOutExtension ( "spatial" ) # settings ws_norm = r 'D:\Modis250\NDVI_10dias_norm' ws_mean = r 'D:\Modis250\NDVI_10dias_mean' ws_anom = r 'D:\Modis250\NDVI_10dias_anom' # create a list of all raster in norm folder arcpy . env . workspace = ws_norm lst_ndvi_ras = arcpy . ListRasters ( ) lst_ndvi_ras . sort ( ) print ( "lst_ndvi_ras: {}" . format ( lst_ndvi_ras ) ) # get list of years decades = GetListOfDecades ( lst_ndvi_ras ) decades . sort ( ) print ( "decades: {}" . format ( decades ) ) # loop through each year for decade in decades : # Get rasters for given year lst_ndvi_decade = GetListOfRasterForGivenDecade ( lst_ndvi_ras , decade ) print ( "lst_ndvi_decade: {}" . format ( lst_ndvi_decade ) ) # calculate mean raster for decade lst_ndvi_decade_paths = [ os . path . join ( ws_norm , r ) for r in lst_ndvi_decade ] print ( "lst_ndvi_decade_paths: {}" . format ( lst_ndvi_decade_paths ) ) mean_ras = arcpy . sa . CellStatistics ( lst_ndvi_decade_paths , "MEAN" , "DATA" ) # store mean ras for decade # mean_dec_d_.tif out_name_mean = "mean_dec_{}.tif" . format ( decade ) ### out_name_path_mean = os . path . join ( ws_mean , out_name_mean ) print ( "out_name_path_mean: {}" . format ( out_name_path_mean ) ) mean_ras . save ( out_name_path_mean ) # loop through each raster in decade for ndvi_ras in lst_ndvi_decade : # calculate anom for each raster in decade # anom_2002_dec_1.tif = norm_ndvi_2002_dec_1.tif - mean_dec_1_.tif; print ( "ndvi_ras: {}" . format ( ndvi_ras ) ) # minus_ras = arcpy.sa.Minus(ndvi_ras, out_name_path_mean) out_name_path_norm = os . path . join ( ws_norm , ndvi_ras ) minus_ras = arcpy . Raster ( out_name_path_mean ) - arcpy . Raster ( out_name_path_norm ) # store anom raster out_name_anom = ndvi_ras . replace ( "norm_ndvi" , "anom" ) out_name_path_anom = os . path . join ( ws_anom , out_name_anom ) minus_ras . save ( out_name_path_anom ) print ( "out_name_path_anom: {}" . format ( out_name_path_anom ) ) # return SA license arcpy . CheckInExtension ( "spatial" ) def GetListOfDecades ( lst ) : # norm_ndvi_yyyy_dec_d.tif lst1 = list ( set ( [ r . split ( "_" ) [ 4 ] for r in lst ] ) ) return [ r . split ( '.' ) [ 0 ] for r in lst1 ] def GetListOfRasterForGivenDecade ( lst , decade ) : lst1 = [ ] search = "_{}.tif" . format ( decade ) for r in lst : if search in r . lower ( ) : lst1 . append ( r ) return lst1 if __name__ == '__main__' : main ( )

Python routine to calculate MODIS NDVI anomaly: Xander Bakker modisndvi calculation ndvi Summary workflow: * Select and download NDVI filtered raster datasets;* Store downloaded original raster datasets at X:\Modis250\Originais\ * Organizing files and folders: (Python script 1) * Rename original raster datasets (MODIS.ndvi.YYYYjjj.yL6000.BOKU.tif) to filename standard “ndvi_YYYY_dec_dd.tif”, where “YYYY” is the year, “jjj” is the Julian day and “dd” is the decade (ten days period). * Transfer renamed files from X:\Modis250\Originais to X:\Modis250\NDVI10dias def main ( ) : import arcpy import os import shutil arcpy . env . overwriteOutput = True ws_in = r 'D:\Modis250\Originais' ws_out = r 'D:\Modis250\NDVI_10dias' arcpy . env . workspace = ws_in lst_ras = arcpy . ListRasters ( "*.tif" ) for ras_name in sorted ( lst_ras ) : # MODIS.ndvi.YYYYjjj.yL6000.BOKU.tif in_name_path = os . path . join ( ws_in , ras_name ) y = ras_name [ 11 : 15 ] j = int ( ras_name [ 15 : 18 ] ) d = GetDecadeFromJulian ( int ( y ) , j ) out_name = "ndvi_{}_dec_{}.tif" . format ( y , d ) out_name_path = os . path . join ( ws_out , out_name ) #os.rename(in_name_path, out_name_path) shutil . move ( in_name_path , out_name_path ) print ( out_name + " renamed" ) def GetDecadeFromJulian ( year , j ) : import datetime date = datetime . datetime ( year -1 , 12 , 31 ) + datetime . timedelta ( j ) m = ( date . month - 1 ) * 3 if date . day 10 : d = 1 elif date . day 20 : d = 2 else : d = 3 return m + d if __name__ == '__main__' : main ( ) * Pre-processing raster datasets: (Python script 2)* Project renamed raster datasets ("ndvi_YYYY_dec_dd.tif") from Lat/Lon WGS 1984 to Lat/Lon SIRGAS2000 using "SAD_1969_To_WGS_1984_14 + SAD_1969_To_SIRGAS_2000_1" geographic transformation. * Overwrite spurious data values: NDVI = Con("%NDVI%" 9998,9999,"%NDVI%")) * Standardize raster datasets: NDVI = (NDVI – Mean(NDVI))/Std(NDVI) * Scale raster datasets to 0 – 100: NDVI = ((NDVI-Min(NDVI))/(Max(NDVI)-Min(NDVI)) * 100 * Mask standardized raster datasets by the area of interest (X:\Modis250\SRTM\mde_sc_90.tif) and output to X:\Modis250\NDVI_10dias_norm def main ( ) : import arcpy import os arcpy . CheckOutExtension ( "spatial" ) arcpy . env . overwriteOutput = True # Set Geoprocessing environments ws_scr = r 'D:\Modis250\Temp' #Scratch workspace ws_in = r 'D:\Modis250\NDVI_10dias' #Input workspace ws_out = r 'D:\Modis250\NDVI_10dias_norm' #Output workspace mask = r 'D:\Modis250\SRTM\mde_sc_90.tif' #Raster used as mask arcpy . env . scratchWorkspace = ws_scr arcpy . env . workspace = ws_in lst_ras = arcpy . ListRasters ( "*.tif" ) lst_ras . sort ( ) # Pre-processing original rasters for raster1 in lst_ras : print ( "Processing: " + raster1 ) #outRaster1= Original raster projected from CGS-WGS1964 to CGS-SIRGAS-2000 outRaster1 = os . path . join ( ws_scr , "rep_{}" . format ( raster1 ) ) # Temporary file arcpy . ProjectRaster_management ( raster1 , outRaster1 , "GEOGCS['GCS_SIRGAS_2000',DATUM['D_SIRGAS_2000',SPHEROID['GRS_1980',6378137.0,298.257222101]],PRIMEM['Greenwich',0.0],UNIT['Degree',0.0174532925199433]]" , "NEAREST" , "1.8000795216263E-03 1.80029836029412E-03" , "'SAD_1969_To_WGS_1984_14 + SAD_1969_To_SIRGAS_2000_1'" , "" , "GEOGCS['GCS_WGS_1984',DATUM['D_WGS_1984',SPHEROID['WGS_1984',6378137.0,298.257223563]],PRIMEM['Greenwich',0.0],UNIT['Degree',0.0174532925199433]]" ) print ( "rep_{}" . format ( raster1 ) + " projected" ) #outRaster2 = outRaster1 masked by Area of Interest (AOI) outRaster2 = arcpy . sa . ExtractByMask ( outRaster1 , mask ) # Temporary file outRaster2 . save ( os . path . join ( ws_scr , "mask_{}" . format ( raster1 ) ) ) print ( "mask_{}" . format ( raster1 ) + " masked" ) #outRaster3 = Overwrited spurious data values from outRaster2 outRaster3 = arcpy . sa . Con ( outRaster2 1 , 0 , arcpy . sa . Con ( outRaster2 > 9998 , 9999 , outRaster2 ) ) # Temporary file outRaster3 . save ( os . path . join ( ws_scr , "aju_{}" . format ( raster1 ) ) ) print ( "aju_{}" . format ( raster1 ) + " ajusted" ) #outRaster4 = Standardized outRaster3 #ras_mean = arcpy.GetRasterProperties_management(outRaster3, "MEAN") #ras_std = arcpy.GetRasterProperties_management(outRaster3, "STD") ras_mean = float ( arcpy . GetRasterProperties_management ( outRaster3 , "MEAN" ) . getOutput ( 0 ) ) ras_std = float ( arcpy . GetRasterProperties_management ( outRaster3 , "STD" ) . getOutput ( 0 ) ) outRaster4 = ( ( outRaster3 - ras_mean ) / ras_std ) # Temporary file outRaster4 . save ( os . path . join ( ws_scr , "ras_std_{}" . format ( raster1 ) ) ) print ( "ras_std_{}" . format ( raster1 ) + " standardized" ) #outRaster = Final raster scaled to 0 -100 #ras_max = arcpy.GetRasterProperties_management(outRaster4, "MAXIMUM") #ras_min = arcpy.GetRasterProperties_management(outRaster4, "MINIMUM") ras_max = float ( arcpy . GetRasterProperties_management ( outRaster4 , "MAXIMUM" ) . getOutput ( 0 ) ) ras_min = ras_min = float ( arcpy . GetRasterProperties_management ( outRaster4 , "MINIMUM" ) . getOutput ( 0 ) ) outRaster = ( ( outRaster4 - ras_min ) / ( ras_max - ras_min ) * 100 ) # Final result outRaster . save ( os . path . join ( ws_out , "norm_{}" . format ( raster1 ) ) ) print ( "norm_{}" . format ( raster1 ) + " processed" ) arcpy . CheckInExtension ( "spatial" ) if __name__ == '__main__' : main ( ) * Processing data (Python script 3)* Calculate NDVI mean for each decadal: mean_dec_dd.tif = MEAN (ndvi_YYY1_dec_dd.tif, ndvi_YYY2_dec_dd.tif, ndvi_YYY3_dec_dd.tif,…) * Calculate NDVI anomaly for each year by decade: anom_yyyy_dec_dd.tif = ndvi_YYYY_dec_dd.tif - mean_dec_dd.tif def main ( ) : import arcpy import os arcpy . env . overwriteOutput = True # Check out SA license arcpy . CheckOutExtension ( "spatial" ) # settings ws_norm = r 'D:\Modis250\NDVI_10dias_norm' ws_mean = r 'D:\Modis250\NDVI_10dias_mean' ws_anom = r 'D:\Modis250\NDVI_10dias_anom' ## ws_norm = r'C:\GeoNet\Modis250\NDVI_10dias_norm' ## ws_mean = r'C:\GeoNet\Modis250\NDVI_10dias_mean' ## ws_anom = r'C:\GeoNet\Modis250\NDVI_10dias_anom' # create a list of all raster in norm folder arcpy . env . workspace = ws_norm lst_ndvi_ras = arcpy . ListRasters ( ) #print ("Raster list: {}".format(lst_ndvi_ras)) # get list of years decades = GetListOfDecades ( lst_ndvi_ras ) decades . sort ( ) #print ("Decades list: {}".format(decades)) # loop through each year for decade in decades : # Get rasters for given year lst_ndvi_decade = GetListOfRasterForGivenDecade ( lst_ndvi_ras , decade ) print ( "Decade by year list: {}" . format ( lst_ndvi_decade ) ) # calculate mean raster for decade lst_ndvi_decade_paths = [ os . path . join ( ws_norm , r ) for r in lst_ndvi_decade ] #print("lst_ndvi_decade_paths: {}".format(lst_ndvi_decade_paths)) mean_ras = arcpy . sa . CellStatistics ( lst_ndvi_decade_paths , "MEAN" , "DATA" ) # store mean ras for decade # mean_dec_d_.tif out_name_mean = "mean_dec_{}.tif" . format ( decade ) ### out_name_path_mean = os . path . join ( ws_mean , out_name_mean ) print ( "Mean raster: {}" . format ( out_name_mean ) ) mean_ras . save ( out_name_path_mean ) # loop through each raster in decade for ndvi_ras in lst_ndvi_decade : # calculate anom for each raster in decade # anom_2002_dec_1.tif = norm_ndvi_2002_dec_1.tif - mean_dec_1_.tif; print ( "Processing: {}" . format ( ndvi_ras ) ) # minus_ras = arcpy.sa.Minus(ndvi_ras, out_name_path_mean) out_name_path_norm = os . path . join ( ws_norm , ndvi_ras ) minus_ras = arcpy . Raster ( out_name_path_norm ) - arcpy . Raster ( out_name_path_mean ) # store anom raster out_name_anom = ndvi_ras . replace ( "norm_ndvi" , "anom" ) out_name_path_anom = os . path . join ( ws_anom , out_name_anom ) minus_ras . save ( out_name_path_anom ) print ( "Anom: {}" . format ( out_name_path_anom ) ) # return SA license arcpy . CheckInExtension ( "spatial" ) def GetListOfDecades ( lst ) : # norm_ndvi_yyyy_dec_d.tif lst1 = list ( set ( [ r . split ( "_" ) [ 4 ] for r in lst ] ) ) return [ r . split ( '.' ) [ 0 ] for r in lst1 ] def GetListOfRasterForGivenDecade ( lst , decade ) : lst1 = [ ] search = "_{}.tif" . format ( decade ) for r in lst : if search in r . lower ( ) : lst1 . append ( r ) return lst1 if __name__ == '__main__' : main ( ) Batch file to run on prompt: @ECHO OFF REM Runs NDVI anomaly scripts ECHO Backup original files xcopy D:\Modis250\Originais\*.tif D:\Modis250\Originais_bkp /a /c /i /y /z PAUSE ECHO Selecting decade and renaming original files to "ndvi_YYYY_dec_dd.tif" standard. C:\Python27\ArcGIS10.5\Python.exe D:\Modis250\Python\Script1.py ECHO NDVI files renamed, hit enter to pre-processing data. PAUSE ECHO Projecting, removing spurious data, standardizing and scaling ndvi decade files to 0-100. C:\Python27\ArcGIS10.5\Python.exe D:\Modis250\Python\Script2.py ECHO NDVI pre-processed, hit enter to calculate anomalies. PAUSE ECHO Calculating decade anomalies. C:\Python27\ArcGIS10.5\Python.exe D:\Modis250\Python\Script3.py Del D:\Temp /s Del D:\NDVI_10dias /s ECHO Done. PAUSE

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 routine to calculate MODIS NDVI anomaly

LuizVianna

Xander Bakker‌

Python script to calculate NDVI anomalies

https://geonet.esri.com/people/xander_bakker‌

I´d like to automate all my workflow, from downloaded original raster datasets to anomaly calculator:

3 years data and folder structure can be downloaded here.

Question 1: What is the best option: to run separated scripts or create a unique code?

Question 2: How to translate the model builder to python script (Python script 2)?

Question 3: If the best option is to run separated, how to make it in a logical sequence (script 1 -> script 2 -> script 3)

Summary workflow:

Select and download NDVI filtered raster datasets;
1. Store downloaded original raster datasets at X:\Modis250\Originais\
2. Organizing files and folders: (Python script 1)
3. Rename original raster datasets (MODIS.ndvi.YYYYjjj.yL6000.BOKU.tif) to filename standard “ndvi_YYYY_dec_dd.tif”, where “YYYY” is the year, “jjj” is the Julian day and “dd” is the decade (ten days period).
4. Transfer renamed files from X:\Modis250\Originais to X:\Modis250\NDVI10dias
Pre-processing raster datasets: (Python script 2)
1. Project renamed raster datasets ("ndvi_YYYY_dec_dd.tif") from Lat/Lon WGS 1984 to Lat/Lon SIRGAS2000 using "SAD_1969_To_WGS_1984_14 + SAD_1969_To_SIRGAS_2000_1" geographic transformation.
2. Overwrite spurious data values: NDVI = Con("%NDVI%"<1,0,Con("%NDVI%">9998,9999,"%NDVI%"))
3. Standardize raster datasets: NDVI = (NDVI – Mean(NDVI))/Std(NDVI)
4. Scale raster datasets to 0 – 100: NDVI = ((NDVI-Min(NDVI))/(Max(NDVI)-Min(NDVI)) * 100
5. Mask standardized raster datasets by the area of interest (X:\Modis250\SRTM\mde_sc_90.tif) and output to X:\Modis250\NDVI_10dias_norm
Processing data (Python script 3)
1. Calculate NDVI mean for each decadal: mean_dec_dd.tif = MEAN (ndvi_YYY1_dec_dd.tif, ndvi_YYY2_dec_dd.tif, ndvi_YYY3_dec_dd.tif,…)
2. Calculate NDVI anomaly for each year by decade: anom_yyyy_dec_dd.tif = ndvi_YYYY_dec_dd.tif - mean_dec_dd.tif

Python scripts 1 and 3 are working. The scrip 2 (Pre-processing) is a model in a toolbox.

Python script 1: written by myself, maybe need to be improved!

<SPAN class="comment token">#Rename original files</SPAN>
<SPAN class="keyword token">import</SPAN> arcpy<SPAN class="punctuation token">,</SPAN> datetime<SPAN class="punctuation token">,</SPAN> glob<SPAN class="punctuation token">,</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>
arcpy<SPAN class="punctuation token">.</SPAN>env<SPAN class="punctuation token">.</SPAN>workspace<SPAN class="operator token">=</SPAN>r<SPAN class="string token">'O:\Modis250\Originais'</SPAN>
infolder<SPAN class="operator token">=</SPAN><SPAN class="string token">'O:\Modis250\Originals'</SPAN>
outfolder<SPAN class="operator token">=</SPAN><SPAN class="string token">'O:\Modis250\NDVI_10dias'</SPAN>
list1<SPAN class="operator token">=</SPAN>arcpy<SPAN class="punctuation token">.</SPAN>ListRasters<SPAN class="punctuation token">(</SPAN><SPAN class="string token">"*.tif"</SPAN><SPAN class="punctuation token">)</SPAN>
list1<SPAN class="punctuation token">.</SPAN>sort<SPAN class="punctuation token">(</SPAN><SPAN class="punctuation token">)</SPAN>
&nbsp;&nbsp;&nbsp; 
<SPAN class="keyword token">for</SPAN> raster <SPAN class="keyword token">in</SPAN> list1<SPAN class="punctuation token">:</SPAN>
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;source_path <SPAN class="operator token">=</SPAN> os<SPAN class="punctuation token">.</SPAN>path<SPAN class="punctuation token">.</SPAN>join<SPAN class="punctuation token">(</SPAN>infolder<SPAN class="punctuation token">,</SPAN> raster<SPAN class="punctuation token">)</SPAN>
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;oldFilename<SPAN class="operator token">=</SPAN>raster
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;dd <SPAN class="operator token">=</SPAN> int<SPAN class="punctuation token">(</SPAN>oldFilename<SPAN class="punctuation token">[</SPAN><SPAN class="number token">15</SPAN><SPAN class="punctuation token">:</SPAN><SPAN class="number token">18</SPAN><SPAN class="punctuation token">]</SPAN><SPAN class="punctuation token">)</SPAN>
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;<SPAN class="keyword token">if</SPAN> dd <SPAN class="operator token"><</SPAN> <SPAN class="number token">11</SPAN><SPAN class="punctuation token">:</SPAN> dd <SPAN class="operator token">=</SPAN> <SPAN class="number token">1</SPAN>
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;<SPAN class="keyword token">elif</SPAN> dd <SPAN class="operator token"><</SPAN> <SPAN class="number token">21</SPAN><SPAN class="punctuation token">:</SPAN> dd <SPAN class="operator token">=</SPAN> <SPAN class="number token">2</SPAN>
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;<SPAN class="keyword token">elif</SPAN> dd <SPAN class="operator token"><</SPAN> <SPAN class="number token">31</SPAN><SPAN class="punctuation token">:</SPAN> dd <SPAN class="operator token">=</SPAN> <SPAN class="number token">3</SPAN>
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;<SPAN class="keyword token">elif</SPAN> dd <SPAN class="operator token"><</SPAN> <SPAN class="number token">41</SPAN><SPAN class="punctuation token">:</SPAN> dd <SPAN class="operator token">=</SPAN> <SPAN class="number token">4</SPAN>
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;<SPAN class="keyword token">elif</SPAN> dd <SPAN class="operator token"><</SPAN> <SPAN class="number token">51</SPAN><SPAN class="punctuation token">:</SPAN> dd <SPAN class="operator token">=</SPAN> <SPAN class="number token">5</SPAN>
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;<SPAN class="keyword token">elif</SPAN> dd <SPAN class="operator token"><</SPAN> <SPAN class="number token">61</SPAN><SPAN class="punctuation token">:</SPAN> dd <SPAN class="operator token">=</SPAN> <SPAN class="number token">6</SPAN>
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;<SPAN class="keyword token">elif</SPAN> dd <SPAN class="operator token"><</SPAN> <SPAN class="number token">71</SPAN><SPAN class="punctuation token">:</SPAN> dd <SPAN class="operator token">=</SPAN> <SPAN class="number token">7</SPAN>
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;<SPAN class="keyword token">elif</SPAN> dd <SPAN class="operator token"><</SPAN> <SPAN class="number token">81</SPAN><SPAN class="punctuation token">:</SPAN> dd <SPAN class="operator token">=</SPAN> <SPAN class="number token">8</SPAN>
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;<SPAN class="keyword token">elif</SPAN> dd <SPAN class="operator token"><</SPAN> <SPAN class="number token">91</SPAN><SPAN class="punctuation token">:</SPAN> dd <SPAN class="operator token">=</SPAN> <SPAN class="number token">9</SPAN>
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;<SPAN class="keyword token">elif</SPAN> dd <SPAN class="operator token"><</SPAN> <SPAN class="number token">101</SPAN><SPAN class="punctuation token">:</SPAN> dd <SPAN class="operator token">=</SPAN> <SPAN class="number token">10</SPAN>
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;<SPAN class="keyword token">elif</SPAN> dd <SPAN class="operator token"><</SPAN> <SPAN class="number token">111</SPAN><SPAN class="punctuation token">:</SPAN> dd <SPAN class="operator token">=</SPAN> <SPAN class="number token">11</SPAN>
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;<SPAN class="keyword token">elif</SPAN> dd <SPAN class="operator token"><</SPAN> <SPAN class="number token">121</SPAN><SPAN class="punctuation token">:</SPAN> dd <SPAN class="operator token">=</SPAN> <SPAN class="number token">12</SPAN>
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;<SPAN class="keyword token">elif</SPAN> dd <SPAN class="operator token"><</SPAN> <SPAN class="number token">131</SPAN><SPAN class="punctuation token">:</SPAN> dd <SPAN class="operator token">=</SPAN> <SPAN class="number token">13</SPAN>
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;<SPAN class="keyword token">elif</SPAN> dd <SPAN class="operator token"><</SPAN> <SPAN class="number token">141</SPAN><SPAN class="punctuation token">:</SPAN> dd <SPAN class="operator token">=</SPAN> <SPAN class="number token">14</SPAN>
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;<SPAN class="keyword token">elif</SPAN> dd <SPAN class="operator token"><</SPAN> <SPAN class="number token">151</SPAN><SPAN class="punctuation token">:</SPAN> dd <SPAN class="operator token">=</SPAN> <SPAN class="number token">15</SPAN>
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;<SPAN class="keyword token">elif</SPAN> dd <SPAN class="operator token"><</SPAN> <SPAN class="number token">161</SPAN><SPAN class="punctuation token">:</SPAN> dd <SPAN class="operator token">=</SPAN> <SPAN class="number token">16</SPAN>
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;<SPAN class="keyword token">elif</SPAN> dd <SPAN class="operator token"><</SPAN> <SPAN class="number token">171</SPAN><SPAN class="punctuation token">:</SPAN> dd <SPAN class="operator token">=</SPAN> <SPAN class="number token">17</SPAN>
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;<SPAN class="keyword token">elif</SPAN> dd <SPAN class="operator token"><</SPAN> <SPAN class="number token">181</SPAN><SPAN class="punctuation token">:</SPAN> dd <SPAN class="operator token">=</SPAN> <SPAN class="number token">18</SPAN>
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;<SPAN class="keyword token">elif</SPAN> dd <SPAN class="operator token"><</SPAN> <SPAN class="number token">191</SPAN><SPAN class="punctuation token">:</SPAN> dd <SPAN class="operator token">=</SPAN> <SPAN class="number token">19</SPAN>
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;<SPAN class="keyword token">elif</SPAN> dd <SPAN class="operator token"><</SPAN> <SPAN class="number token">201</SPAN><SPAN class="punctuation token">:</SPAN> dd <SPAN class="operator token">=</SPAN> <SPAN class="number token">20</SPAN>
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;<SPAN class="keyword token">elif</SPAN> dd <SPAN class="operator token"><</SPAN> <SPAN class="number token">211</SPAN><SPAN class="punctuation token">:</SPAN> dd <SPAN class="operator token">=</SPAN> <SPAN class="number token">21</SPAN>
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;<SPAN class="keyword token">elif</SPAN> dd <SPAN class="operator token"><</SPAN> <SPAN class="number token">221</SPAN><SPAN class="punctuation token">:</SPAN> dd <SPAN class="operator token">=</SPAN> <SPAN class="number token">22</SPAN>
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;<SPAN class="keyword token">elif</SPAN> dd <SPAN class="operator token"><</SPAN> <SPAN class="number token">231</SPAN><SPAN class="punctuation token">:</SPAN> dd <SPAN class="operator token">=</SPAN> <SPAN class="number token">23</SPAN>
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;<SPAN class="keyword token">elif</SPAN> dd <SPAN class="operator token"><</SPAN> <SPAN class="number token">241</SPAN><SPAN class="punctuation token">:</SPAN> dd <SPAN class="operator token">=</SPAN> <SPAN class="number token">24</SPAN>
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;<SPAN class="keyword token">elif</SPAN> dd <SPAN class="operator token"><</SPAN> <SPAN class="number token">251</SPAN><SPAN class="punctuation token">:</SPAN> dd <SPAN class="operator token">=</SPAN> <SPAN class="number token">25</SPAN>
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;<SPAN class="keyword token">elif</SPAN> dd <SPAN class="operator token"><</SPAN> <SPAN class="number token">261</SPAN><SPAN class="punctuation token">:</SPAN> dd <SPAN class="operator token">=</SPAN> <SPAN class="number token">26</SPAN>
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;<SPAN class="keyword token">elif</SPAN> dd <SPAN class="operator token"><</SPAN> <SPAN class="number token">271</SPAN><SPAN class="punctuation token">:</SPAN> dd <SPAN class="operator token">=</SPAN> <SPAN class="number token">27</SPAN>
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;<SPAN class="keyword token">elif</SPAN> dd <SPAN class="operator token"><</SPAN> <SPAN class="number token">281</SPAN><SPAN class="punctuation token">:</SPAN> dd <SPAN class="operator token">=</SPAN> <SPAN class="number token">28</SPAN>
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;<SPAN class="keyword token">elif</SPAN> dd <SPAN class="operator token"><</SPAN> <SPAN class="number token">291</SPAN><SPAN class="punctuation token">:</SPAN> dd <SPAN class="operator token">=</SPAN> <SPAN class="number token">29</SPAN>
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;<SPAN class="keyword token">elif</SPAN> dd <SPAN class="operator token"><</SPAN> <SPAN class="number token">301</SPAN><SPAN class="punctuation token">:</SPAN> dd <SPAN class="operator token">=</SPAN> <SPAN class="number token">30</SPAN>
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;<SPAN class="keyword token">elif</SPAN> dd <SPAN class="operator token"><</SPAN> <SPAN class="number token">311</SPAN><SPAN class="punctuation token">:</SPAN> dd <SPAN class="operator token">=</SPAN> <SPAN class="number token">31</SPAN>
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;<SPAN class="keyword token">elif</SPAN> dd <SPAN class="operator token"><</SPAN> <SPAN class="number token">321</SPAN><SPAN class="punctuation token">:</SPAN> dd <SPAN class="operator token">=</SPAN> <SPAN class="number token">32</SPAN>
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;<SPAN class="keyword token">elif</SPAN> dd <SPAN class="operator token"><</SPAN> <SPAN class="number token">331</SPAN><SPAN class="punctuation token">:</SPAN> dd <SPAN class="operator token">=</SPAN> <SPAN class="number token">33</SPAN>
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;<SPAN class="keyword token">elif</SPAN> dd <SPAN class="operator token"><</SPAN> <SPAN class="number token">341</SPAN><SPAN class="punctuation token">:</SPAN> dd <SPAN class="operator token">=</SPAN> <SPAN class="number token">34</SPAN>
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;<SPAN class="keyword token">elif</SPAN> dd <SPAN class="operator token"><</SPAN> <SPAN class="number token">351</SPAN><SPAN class="punctuation token">:</SPAN> dd <SPAN class="operator token">=</SPAN> <SPAN class="number token">35</SPAN>
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;<SPAN class="keyword token">else</SPAN><SPAN class="punctuation token">:</SPAN> dd <SPAN class="operator token">=</SPAN> <SPAN class="number token">36</SPAN>
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;dd <SPAN class="operator token">=</SPAN> str<SPAN class="punctuation token">(</SPAN>dd<SPAN class="punctuation token">)</SPAN>
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;newFilename<SPAN class="operator token">=</SPAN><SPAN class="string token">"ndvi_"</SPAN> <SPAN class="operator token">+</SPAN> oldFilename<SPAN class="punctuation token">[</SPAN><SPAN class="number token">11</SPAN><SPAN class="punctuation token">:</SPAN><SPAN class="number token">15</SPAN><SPAN class="punctuation token">]</SPAN> <SPAN class="operator token">+</SPAN> <SPAN class="string token">"_dec_"</SPAN> <SPAN class="operator token">+</SPAN> dd <SPAN class="operator token">+</SPAN> <SPAN class="string token">".tif"</SPAN>
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;destination_path<SPAN class="operator token">=</SPAN>os<SPAN class="punctuation token">.</SPAN>path<SPAN class="punctuation token">.</SPAN>join<SPAN class="punctuation token">(</SPAN>outfolder<SPAN class="punctuation token">,</SPAN> newFilename<SPAN class="punctuation token">)</SPAN>
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;os<SPAN class="punctuation token">.</SPAN>rename<SPAN class="punctuation token">(</SPAN>source_path<SPAN class="punctuation token">,</SPAN> destination_path<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>

Model Builder -> Python script 2

Python script 3: Python script to calculate NDVI anomalies

<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>
&nbsp;&nbsp;&nbsp; ws_norm <SPAN class="operator token">=</SPAN> r<SPAN class="string token">'D:\Modis250\NDVI_10dias_norm'</SPAN>
&nbsp;&nbsp;&nbsp; ws_mean <SPAN class="operator token">=</SPAN> r<SPAN class="string token">'D:\Modis250\NDVI_10dias_mean'</SPAN>
&nbsp;&nbsp;&nbsp; ws_anom <SPAN class="operator token">=</SPAN> r<SPAN class="string token">'D:\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;&nbsp;&nbsp;lst_ndvi_ras<SPAN class="punctuation token">.</SPAN>sort<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;&nbsp;&nbsp;decades<SPAN class="punctuation token">.</SPAN>sort<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">"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>

Trouver plus de publications taguées avec

Accepted answers

LuizVianna

Python routine to calculate MODIS NDVI anomaly:

Xander Bakker‌

modis‌ndvi calculation ndvi

Summary workflow:

Select and download NDVI filtered raster datasets;
1. Store downloaded original raster datasets at X:\Modis250\Originais\
2. Organizing files and folders: (Python script 1)
3. Rename original raster datasets (MODIS.ndvi.YYYYjjj.yL6000.BOKU.tif) to filename standard “ndvi_YYYY_dec_dd.tif”, where “YYYY” is the year, “jjj” is the Julian day and “dd” is the decade (ten days period).
4. Transfer renamed files from X:\Modis250\Originais to X:\Modis250\NDVI10dias

<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;&nbsp; <SPAN class="keyword token">import</SPAN> arcpy
&nbsp;&nbsp;&nbsp;&nbsp; <SPAN class="keyword token">import</SPAN> os
&nbsp;&nbsp;&nbsp;&nbsp; <SPAN class="keyword token">import</SPAN> shutil

&nbsp;&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;&nbsp; ws_in <SPAN class="operator token">=</SPAN> r<SPAN class="string token">'D:\Modis250\Originais'</SPAN>
&nbsp;&nbsp;&nbsp;&nbsp; ws_out <SPAN class="operator token">=</SPAN> r<SPAN class="string token">'D:\Modis250\NDVI_10dias'</SPAN>

&nbsp;&nbsp;&nbsp;&nbsp; arcpy<SPAN class="punctuation token">.</SPAN>env<SPAN class="punctuation token">.</SPAN>workspace <SPAN class="operator token">=</SPAN> ws_in
&nbsp;&nbsp;&nbsp;&nbsp; lst_ras <SPAN class="operator token">=</SPAN> arcpy<SPAN class="punctuation token">.</SPAN>ListRasters<SPAN class="punctuation token">(</SPAN><SPAN class="string token">"*.tif"</SPAN><SPAN class="punctuation token">)</SPAN>

&nbsp;&nbsp;&nbsp;&nbsp; <SPAN class="keyword token">for</SPAN> ras_name <SPAN class="keyword token">in</SPAN> sorted<SPAN class="punctuation token">(</SPAN>lst_ras<SPAN class="punctuation token">)</SPAN><SPAN class="punctuation token">:</SPAN>
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; <SPAN class="comment token"># MODIS.ndvi.YYYYjjj.yL6000.BOKU.tif</SPAN>
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; in_name_path <SPAN class="operator token">=</SPAN> os<SPAN class="punctuation token">.</SPAN>path<SPAN class="punctuation token">.</SPAN>join<SPAN class="punctuation token">(</SPAN>ws_in<SPAN class="punctuation token">,</SPAN> ras_name<SPAN class="punctuation token">)</SPAN>
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; y <SPAN class="operator token">=</SPAN> ras_name<SPAN class="punctuation token">[</SPAN><SPAN class="number token">11</SPAN><SPAN class="punctuation token">:</SPAN><SPAN class="number token">15</SPAN><SPAN class="punctuation token">]</SPAN>
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; j <SPAN class="operator token">=</SPAN> int<SPAN class="punctuation token">(</SPAN>ras_name<SPAN class="punctuation token">[</SPAN><SPAN class="number token">15</SPAN><SPAN class="punctuation token">:</SPAN><SPAN class="number token">18</SPAN><SPAN class="punctuation token">]</SPAN><SPAN class="punctuation token">)</SPAN>
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; d <SPAN class="operator token">=</SPAN> GetDecadeFromJulian<SPAN class="punctuation token">(</SPAN>int<SPAN class="punctuation token">(</SPAN>y<SPAN class="punctuation token">)</SPAN><SPAN class="punctuation token">,</SPAN> j<SPAN class="punctuation token">)</SPAN>

&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; out_name <SPAN class="operator token">=</SPAN> <SPAN class="string token">"ndvi_{}_dec_{}.tif"</SPAN><SPAN class="punctuation token">.</SPAN>format<SPAN class="punctuation token">(</SPAN>y<SPAN class="punctuation token">,</SPAN> d<SPAN class="punctuation token">)</SPAN>
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; out_name_path <SPAN class="operator token">=</SPAN> os<SPAN class="punctuation token">.</SPAN>path<SPAN class="punctuation token">.</SPAN>join<SPAN class="punctuation token">(</SPAN>ws_out<SPAN class="punctuation token">,</SPAN> out_name<SPAN class="punctuation token">)</SPAN>
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; <SPAN class="comment token">#os.rename(in_name_path, out_name_path)</SPAN>
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; shutil<SPAN class="punctuation token">.</SPAN>move<SPAN class="punctuation token">(</SPAN>in_name_path<SPAN class="punctuation token">,</SPAN> out_name_path<SPAN class="punctuation token">)</SPAN>
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; <SPAN class="keyword token">print</SPAN> <SPAN class="punctuation token">(</SPAN>out_name <SPAN class="operator token">+</SPAN> <SPAN class="string token">" renamed"</SPAN><SPAN class="punctuation token">)</SPAN>

<SPAN class="keyword token">def</SPAN> <SPAN class="token function">GetDecadeFromJulian</SPAN><SPAN class="punctuation token">(</SPAN>year<SPAN class="punctuation token">,</SPAN> j<SPAN class="punctuation token">)</SPAN><SPAN class="punctuation token">:</SPAN>
&nbsp;&nbsp;&nbsp;&nbsp; <SPAN class="keyword token">import</SPAN> datetime
&nbsp;&nbsp;&nbsp;&nbsp; date <SPAN class="operator token">=</SPAN> datetime<SPAN class="punctuation token">.</SPAN>datetime<SPAN class="punctuation token">(</SPAN>year<SPAN class="number token">-1</SPAN><SPAN class="punctuation token">,</SPAN> <SPAN class="number token">12</SPAN><SPAN class="punctuation token">,</SPAN> <SPAN class="number token">31</SPAN><SPAN class="punctuation token">)</SPAN> <SPAN class="operator token">+</SPAN> datetime<SPAN class="punctuation token">.</SPAN>timedelta<SPAN class="punctuation token">(</SPAN>j<SPAN class="punctuation token">)</SPAN>
&nbsp;&nbsp;&nbsp;&nbsp; m <SPAN class="operator token">=</SPAN> <SPAN class="punctuation token">(</SPAN>date<SPAN class="punctuation token">.</SPAN>month <SPAN class="operator token">-</SPAN> <SPAN class="number token">1</SPAN><SPAN class="punctuation token">)</SPAN> <SPAN class="operator token">*</SPAN> <SPAN class="number token">3</SPAN>
&nbsp;&nbsp;&nbsp;&nbsp; <SPAN class="keyword token">if</SPAN> date<SPAN class="punctuation token">.</SPAN>day <SPAN class="operator token"><=</SPAN> <SPAN class="number token">10</SPAN><SPAN class="punctuation token">:</SPAN>
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; d <SPAN class="operator token">=</SPAN> <SPAN class="number token">1</SPAN>
&nbsp;&nbsp;&nbsp;&nbsp; <SPAN class="keyword token">elif</SPAN> date<SPAN class="punctuation token">.</SPAN>day <SPAN class="operator token"><=</SPAN> <SPAN class="number token">20</SPAN><SPAN class="punctuation token">:</SPAN>
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; d <SPAN class="operator token">=</SPAN> <SPAN class="number token">2</SPAN>
&nbsp;&nbsp;&nbsp;&nbsp; <SPAN class="keyword token">else</SPAN><SPAN class="punctuation token">:</SPAN>
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; d <SPAN class="operator token">=</SPAN> <SPAN class="number token">3</SPAN>
&nbsp;&nbsp;&nbsp;&nbsp; <SPAN class="keyword token">return</SPAN> m <SPAN class="operator token">+</SPAN> d

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

Pre-processing raster datasets: (Python script 2)
1. Project renamed raster datasets ("ndvi_YYYY_dec_dd.tif") from Lat/Lon WGS 1984 to Lat/Lon SIRGAS2000 using "SAD_1969_To_WGS_1984_14 + SAD_1969_To_SIRGAS_2000_1" geographic transformation.
2. Overwrite spurious data values: NDVI = Con("%NDVI%"<1,0,Con("%NDVI%">9998,9999,"%NDVI%"))
3. Standardize raster datasets: NDVI = (NDVI – Mean(NDVI))/Std(NDVI)
4. Scale raster datasets to 0 – 100: NDVI = ((NDVI-Min(NDVI))/(Max(NDVI)-Min(NDVI)) * 100
5. Mask standardized raster datasets by the area of interest (X:\Modis250\SRTM\mde_sc_90.tif) and output to X:\Modis250\NDVI_10dias_norm

<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;&nbsp; <SPAN class="keyword token">import</SPAN> arcpy
&nbsp;&nbsp;&nbsp;&nbsp; <SPAN class="keyword token">import</SPAN> os
&nbsp;&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;&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;&nbsp; <SPAN class="comment token"># Set Geoprocessing environments</SPAN>
&nbsp;&nbsp;&nbsp;&nbsp; ws_scr <SPAN class="operator token">=</SPAN> r<SPAN class="string token">'D:\Modis250\Temp'</SPAN> <SPAN class="comment token">#Scratch workspace</SPAN>
&nbsp;&nbsp;&nbsp;&nbsp; ws_in <SPAN class="operator token">=</SPAN> r<SPAN class="string token">'D:\Modis250\NDVI_10dias'</SPAN> <SPAN class="comment token">#Input workspace</SPAN>
&nbsp;&nbsp;&nbsp;&nbsp; ws_out <SPAN class="operator token">=</SPAN> r<SPAN class="string token">'D:\Modis250\NDVI_10dias_norm'</SPAN> <SPAN class="comment token">#Output workspace</SPAN>
&nbsp;&nbsp;&nbsp;&nbsp; mask <SPAN class="operator token">=</SPAN> r<SPAN class="string token">'D:\Modis250\SRTM\mde_sc_90.tif'</SPAN> <SPAN class="comment token">#Raster used as mask</SPAN>
&nbsp;&nbsp;&nbsp;&nbsp; arcpy<SPAN class="punctuation token">.</SPAN>env<SPAN class="punctuation token">.</SPAN>scratchWorkspace <SPAN class="operator token">=</SPAN> ws_scr
&nbsp;&nbsp;&nbsp;&nbsp; arcpy<SPAN class="punctuation token">.</SPAN>env<SPAN class="punctuation token">.</SPAN>workspace <SPAN class="operator token">=</SPAN> ws_in
&nbsp;&nbsp;&nbsp;&nbsp; lst_ras <SPAN class="operator token">=</SPAN> arcpy<SPAN class="punctuation token">.</SPAN>ListRasters<SPAN class="punctuation token">(</SPAN><SPAN class="string token">"*.tif"</SPAN><SPAN class="punctuation token">)</SPAN>
&nbsp;&nbsp;&nbsp;&nbsp; lst_ras<SPAN class="punctuation token">.</SPAN>sort<SPAN class="punctuation token">(</SPAN><SPAN class="punctuation token">)</SPAN>

&nbsp;&nbsp;&nbsp; <SPAN class="comment token"># Pre-processing original rasters</SPAN>
&nbsp;&nbsp;&nbsp;&nbsp; <SPAN class="keyword token">for</SPAN> raster1 <SPAN class="keyword token">in</SPAN> lst_ras<SPAN class="punctuation token">:</SPAN>
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; <SPAN class="keyword token">print</SPAN> <SPAN class="punctuation token">(</SPAN><SPAN class="string token">"Processing: "</SPAN> <SPAN class="operator token">+</SPAN> raster1<SPAN class="punctuation token">)</SPAN>
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; <SPAN class="comment token">#outRaster1= Original raster projected from CGS-WGS1964 to CGS-SIRGAS-2000&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; </SPAN>
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; outRaster1 <SPAN class="operator token">=</SPAN> os<SPAN class="punctuation token">.</SPAN>path<SPAN class="punctuation token">.</SPAN>join<SPAN class="punctuation token">(</SPAN>ws_scr<SPAN class="punctuation token">,</SPAN> <SPAN class="string token">"rep_{}"</SPAN><SPAN class="punctuation token">.</SPAN>format<SPAN class="punctuation token">(</SPAN>raster1<SPAN class="punctuation token">)</SPAN><SPAN class="punctuation token">)</SPAN>&nbsp; <SPAN class="comment token"># Temporary file</SPAN>
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; arcpy<SPAN class="punctuation token">.</SPAN>ProjectRaster_management<SPAN class="punctuation token">(</SPAN>raster1<SPAN class="punctuation token">,</SPAN> outRaster1<SPAN class="punctuation token">,</SPAN><SPAN class="string token">"GEOGCS['GCS_SIRGAS_2000',DATUM['D_SIRGAS_2000',SPHEROID['GRS_1980',6378137.0,298.257222101]],PRIMEM['Greenwich',0.0],UNIT['Degree',0.0174532925199433]]"</SPAN><SPAN class="punctuation token">,</SPAN> <SPAN class="string token">"NEAREST"</SPAN><SPAN class="punctuation token">,</SPAN> <SPAN class="string token">"1.8000795216263E-03 1.80029836029412E-03"</SPAN><SPAN class="punctuation token">,</SPAN> <SPAN class="string token">"'SAD_1969_To_WGS_1984_14 + SAD_1969_To_SIRGAS_2000_1'"</SPAN><SPAN class="punctuation token">,</SPAN> <SPAN class="string token">""</SPAN><SPAN class="punctuation token">,</SPAN> <SPAN class="string token">"GEOGCS['GCS_WGS_1984',DATUM['D_WGS_1984',SPHEROID['WGS_1984',6378137.0,298.257223563]],PRIMEM['Greenwich',0.0],UNIT['Degree',0.0174532925199433]]"</SPAN><SPAN class="punctuation token">)</SPAN>
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; <SPAN class="keyword token">print</SPAN> <SPAN class="punctuation token">(</SPAN><SPAN class="string token">"rep_{}"</SPAN><SPAN class="punctuation token">.</SPAN>format<SPAN class="punctuation token">(</SPAN>raster1<SPAN class="punctuation token">)</SPAN> <SPAN class="operator token">+</SPAN> <SPAN class="string token">" projected"</SPAN><SPAN class="punctuation token">)</SPAN>

&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; <SPAN class="comment token">#outRaster2 = outRaster1 masked by Area of Interest (AOI)</SPAN>
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; outRaster2 <SPAN class="operator token">=</SPAN> arcpy<SPAN class="punctuation token">.</SPAN>sa<SPAN class="punctuation token">.</SPAN>ExtractByMask<SPAN class="punctuation token">(</SPAN>outRaster1<SPAN class="punctuation token">,</SPAN> mask<SPAN class="punctuation token">)</SPAN>&nbsp; <SPAN class="comment token"># Temporary file</SPAN>
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; outRaster2<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>ws_scr<SPAN class="punctuation token">,</SPAN> <SPAN class="string token">"mask_{}"</SPAN><SPAN class="punctuation token">.</SPAN>format<SPAN class="punctuation token">(</SPAN>raster1<SPAN class="punctuation token">)</SPAN><SPAN class="punctuation token">)</SPAN><SPAN class="punctuation token">)</SPAN>
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; <SPAN class="keyword token">print</SPAN> <SPAN class="punctuation token">(</SPAN><SPAN class="string token">"mask_{}"</SPAN><SPAN class="punctuation token">.</SPAN>format<SPAN class="punctuation token">(</SPAN>raster1<SPAN class="punctuation token">)</SPAN> <SPAN class="operator token">+</SPAN> <SPAN class="string token">" masked"</SPAN><SPAN class="punctuation token">)</SPAN>

&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; <SPAN class="comment token">#outRaster3 = Overwrited spurious data values from outRaster2</SPAN>
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; outRaster3 <SPAN class="operator token">=</SPAN> arcpy<SPAN class="punctuation token">.</SPAN>sa<SPAN class="punctuation token">.</SPAN>Con<SPAN class="punctuation token">(</SPAN>outRaster2 <SPAN class="operator token"><</SPAN> <SPAN class="number token">1</SPAN><SPAN class="punctuation token">,</SPAN> <SPAN class="number token">0</SPAN><SPAN class="punctuation token">,</SPAN> arcpy<SPAN class="punctuation token">.</SPAN>sa<SPAN class="punctuation token">.</SPAN>Con<SPAN class="punctuation token">(</SPAN>outRaster2 <SPAN class="operator token">></SPAN> <SPAN class="number token">9998</SPAN><SPAN class="punctuation token">,</SPAN> <SPAN class="number token">9999</SPAN><SPAN class="punctuation token">,</SPAN> outRaster2<SPAN class="punctuation token">)</SPAN><SPAN class="punctuation token">)</SPAN>&nbsp; <SPAN class="comment token"># Temporary file</SPAN>
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; outRaster3<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>ws_scr<SPAN class="punctuation token">,</SPAN> <SPAN class="string token">"aju_{}"</SPAN><SPAN class="punctuation token">.</SPAN>format<SPAN class="punctuation token">(</SPAN>raster1<SPAN class="punctuation token">)</SPAN><SPAN class="punctuation token">)</SPAN><SPAN class="punctuation token">)</SPAN>
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; <SPAN class="keyword token">print</SPAN> <SPAN class="punctuation token">(</SPAN><SPAN class="string token">"aju_{}"</SPAN><SPAN class="punctuation token">.</SPAN>format<SPAN class="punctuation token">(</SPAN>raster1<SPAN class="punctuation token">)</SPAN> <SPAN class="operator token">+</SPAN> <SPAN class="string token">" ajusted"</SPAN><SPAN class="punctuation token">)</SPAN>

&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; <SPAN class="comment token">#outRaster4 = Standardized outRaster3</SPAN>
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; <SPAN class="comment token">#ras_mean = arcpy.GetRasterProperties_management(outRaster3, "MEAN")</SPAN>
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; <SPAN class="comment token">#ras_std = arcpy.GetRasterProperties_management(outRaster3, "STD")</SPAN>
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; ras_mean <SPAN class="operator token">=</SPAN> float<SPAN class="punctuation token">(</SPAN>arcpy<SPAN class="punctuation token">.</SPAN>GetRasterProperties_management<SPAN class="punctuation token">(</SPAN>outRaster3<SPAN class="punctuation token">,</SPAN> <SPAN class="string token">"MEAN"</SPAN><SPAN class="punctuation token">)</SPAN><SPAN class="punctuation token">.</SPAN>getOutput<SPAN class="punctuation token">(</SPAN><SPAN class="number token">0</SPAN><SPAN class="punctuation token">)</SPAN><SPAN class="punctuation token">)</SPAN>
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; ras_std <SPAN class="operator token">=</SPAN> float<SPAN class="punctuation token">(</SPAN>arcpy<SPAN class="punctuation token">.</SPAN>GetRasterProperties_management<SPAN class="punctuation token">(</SPAN>outRaster3<SPAN class="punctuation token">,</SPAN> <SPAN class="string token">"STD"</SPAN><SPAN class="punctuation token">)</SPAN><SPAN class="punctuation token">.</SPAN>getOutput<SPAN class="punctuation token">(</SPAN><SPAN class="number token">0</SPAN><SPAN class="punctuation token">)</SPAN><SPAN class="punctuation token">)</SPAN>
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; outRaster4 <SPAN class="operator token">=</SPAN> <SPAN class="punctuation token">(</SPAN><SPAN class="punctuation token">(</SPAN>outRaster3<SPAN class="operator token">-</SPAN>ras_mean<SPAN class="punctuation token">)</SPAN><SPAN class="operator token">/</SPAN>ras_std<SPAN class="punctuation token">)</SPAN>&nbsp; <SPAN class="comment token"># Temporary file</SPAN>
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; outRaster4<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>ws_scr<SPAN class="punctuation token">,</SPAN> <SPAN class="string token">"ras_std_{}"</SPAN><SPAN class="punctuation token">.</SPAN>format<SPAN class="punctuation token">(</SPAN>raster1<SPAN class="punctuation token">)</SPAN><SPAN class="punctuation token">)</SPAN><SPAN class="punctuation token">)</SPAN>
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; <SPAN class="keyword token">print</SPAN> <SPAN class="punctuation token">(</SPAN><SPAN class="string token">"ras_std_{}"</SPAN><SPAN class="punctuation token">.</SPAN>format<SPAN class="punctuation token">(</SPAN>raster1<SPAN class="punctuation token">)</SPAN> <SPAN class="operator token">+</SPAN> <SPAN class="string token">" standardized"</SPAN><SPAN class="punctuation token">)</SPAN>

&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; <SPAN class="comment token">#outRaster = Final raster scaled to 0 -100</SPAN>
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; <SPAN class="comment token">#ras_max = arcpy.GetRasterProperties_management(outRaster4, "MAXIMUM")</SPAN>
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; <SPAN class="comment token">#ras_min = arcpy.GetRasterProperties_management(outRaster4, "MINIMUM")</SPAN>
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; ras_max <SPAN class="operator token">=</SPAN> float<SPAN class="punctuation token">(</SPAN>arcpy<SPAN class="punctuation token">.</SPAN>GetRasterProperties_management<SPAN class="punctuation token">(</SPAN>outRaster4<SPAN class="punctuation token">,</SPAN> <SPAN class="string token">"MAXIMUM"</SPAN><SPAN class="punctuation token">)</SPAN><SPAN class="punctuation token">.</SPAN>getOutput<SPAN class="punctuation token">(</SPAN><SPAN class="number token">0</SPAN><SPAN class="punctuation token">)</SPAN><SPAN class="punctuation token">)</SPAN>
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; ras_min <SPAN class="operator token">=</SPAN> ras_min <SPAN class="operator token">=</SPAN> float<SPAN class="punctuation token">(</SPAN>arcpy<SPAN class="punctuation token">.</SPAN>GetRasterProperties_management<SPAN class="punctuation token">(</SPAN>outRaster4<SPAN class="punctuation token">,</SPAN> <SPAN class="string token">"MINIMUM"</SPAN><SPAN class="punctuation token">)</SPAN><SPAN class="punctuation token">.</SPAN>getOutput<SPAN class="punctuation token">(</SPAN><SPAN class="number token">0</SPAN><SPAN class="punctuation token">)</SPAN><SPAN class="punctuation token">)</SPAN>
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; outRaster <SPAN class="operator token">=</SPAN> <SPAN class="punctuation token">(</SPAN><SPAN class="punctuation token">(</SPAN>outRaster4<SPAN class="operator token">-</SPAN>ras_min<SPAN class="punctuation token">)</SPAN><SPAN class="operator token">/</SPAN><SPAN class="punctuation token">(</SPAN>ras_max<SPAN class="operator token">-</SPAN>ras_min<SPAN class="punctuation token">)</SPAN> <SPAN class="operator token">*</SPAN> <SPAN class="number token">100</SPAN><SPAN class="punctuation token">)</SPAN>&nbsp; <SPAN class="comment token"># Final result</SPAN>
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; outRaster<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>ws_out<SPAN class="punctuation token">,</SPAN> <SPAN class="string token">"norm_{}"</SPAN><SPAN class="punctuation token">.</SPAN>format<SPAN class="punctuation token">(</SPAN>raster1<SPAN class="punctuation token">)</SPAN><SPAN class="punctuation token">)</SPAN><SPAN class="punctuation token">)</SPAN>
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; <SPAN class="keyword token">print</SPAN> <SPAN class="punctuation token">(</SPAN><SPAN class="string token">"norm_{}"</SPAN><SPAN class="punctuation token">.</SPAN>format<SPAN class="punctuation token">(</SPAN>raster1<SPAN class="punctuation token">)</SPAN> <SPAN class="operator token">+</SPAN> <SPAN class="string token">" processed"</SPAN><SPAN class="punctuation token">)</SPAN>

&nbsp;&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">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>

Processing data (Python script 3)
1. Calculate NDVI mean for each decadal: mean_dec_dd.tif = MEAN (ndvi_YYY1_dec_dd.tif, ndvi_YYY2_dec_dd.tif, ndvi_YYY3_dec_dd.tif,…)
2. Calculate NDVI anomaly for each year by decade: anom_yyyy_dec_dd.tif = ndvi_YYYY_dec_dd.tif - mean_dec_dd.tif

<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;&nbsp; <SPAN class="keyword token">import</SPAN> arcpy
&nbsp;&nbsp;&nbsp;&nbsp; <SPAN class="keyword token">import</SPAN> os
&nbsp;&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;&nbsp; <SPAN class="comment token"># Check out SA license</SPAN>
&nbsp;&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;&nbsp; <SPAN class="comment token"># settings</SPAN>
&nbsp;&nbsp;&nbsp;&nbsp; ws_norm <SPAN class="operator token">=</SPAN> r<SPAN class="string token">'D:\Modis250\NDVI_10dias_norm'</SPAN>
&nbsp;&nbsp;&nbsp;&nbsp; ws_mean <SPAN class="operator token">=</SPAN> r<SPAN class="string token">'D:\Modis250\NDVI_10dias_mean'</SPAN>
&nbsp;&nbsp;&nbsp;&nbsp; ws_anom <SPAN class="operator token">=</SPAN> r<SPAN class="string token">'D:\Modis250\NDVI_10dias_anom'</SPAN>
<SPAN class="comment token">##&nbsp;&nbsp;&nbsp; ws_norm = r'C:\GeoNet\Modis250\NDVI_10dias_norm'</SPAN>
<SPAN class="comment token">##&nbsp;&nbsp;&nbsp; ws_mean = r'C:\GeoNet\Modis250\NDVI_10dias_mean'</SPAN>
<SPAN class="comment token">##&nbsp;&nbsp;&nbsp; ws_anom = r'C:\GeoNet\Modis250\NDVI_10dias_anom'</SPAN>


&nbsp;&nbsp;&nbsp;&nbsp; <SPAN class="comment token"># create a list of all raster in norm folder</SPAN>
&nbsp;&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;&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;&nbsp; <SPAN class="comment token">#print ("Raster list: {}".format(lst_ndvi_ras))</SPAN>

&nbsp;&nbsp;&nbsp;&nbsp; <SPAN class="comment token"># get list of years</SPAN>
&nbsp;&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;&nbsp; decades<SPAN class="punctuation token">.</SPAN>sort<SPAN class="punctuation token">(</SPAN><SPAN class="punctuation token">)</SPAN>
&nbsp;&nbsp;&nbsp;&nbsp; <SPAN class="comment token">#print ("Decades list: {}".format(decades))</SPAN>

&nbsp;&nbsp;&nbsp;&nbsp; <SPAN class="comment token"># loop through each year</SPAN>
&nbsp;&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;&nbsp;&nbsp; <SPAN class="comment token"># Get rasters for given year</SPAN>
&nbsp;&nbsp;&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;&nbsp;&nbsp; <SPAN class="keyword token">print</SPAN><SPAN class="punctuation token">(</SPAN><SPAN class="string token">"Decade by year list: {}"</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;&nbsp;&nbsp; <SPAN class="comment token"># calculate mean raster for decade</SPAN>
&nbsp;&nbsp;&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;&nbsp;&nbsp; <SPAN class="comment token">#print("lst_ndvi_decade_paths: {}".format(lst_ndvi_decade_paths))</SPAN>
&nbsp;&nbsp;&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;&nbsp;&nbsp; 

&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; <SPAN class="comment token"># store mean ras for decade</SPAN>
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; <SPAN class="comment token"># mean_dec_d_.tif</SPAN>
&nbsp;&nbsp;&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;&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;&nbsp;&nbsp; <SPAN class="keyword token">print</SPAN> <SPAN class="punctuation token">(</SPAN><SPAN class="string token">"Mean raster: {}"</SPAN><SPAN class="punctuation token">.</SPAN>format<SPAN class="punctuation token">(</SPAN>out_name_mean<SPAN class="punctuation token">)</SPAN><SPAN class="punctuation token">)</SPAN>
&nbsp;&nbsp;&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;&nbsp;&nbsp; <SPAN class="comment token"># loop through each raster in decade</SPAN>
&nbsp;&nbsp;&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;&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;&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;&nbsp;&nbsp;&nbsp; <SPAN class="keyword token">print</SPAN> <SPAN class="punctuation token">(</SPAN><SPAN class="string token">"Processing: {}"</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;&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;&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;&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_norm<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_mean<SPAN class="punctuation token">)</SPAN>

&nbsp;&nbsp;&nbsp;&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;&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;&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;&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;&nbsp;&nbsp;&nbsp; <SPAN class="keyword token">print</SPAN> <SPAN class="punctuation token">(</SPAN><SPAN class="string token">"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;&nbsp; <SPAN class="comment token"># return SA license</SPAN>
&nbsp;&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;&nbsp; <SPAN class="comment token"># norm_ndvi_yyyy_dec_d.tif</SPAN>
&nbsp;&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;&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;&nbsp; lst1 <SPAN class="operator token">=</SPAN> <SPAN class="punctuation token">[</SPAN><SPAN class="punctuation token">]</SPAN>
&nbsp;&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;&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;&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;&nbsp;&nbsp;&nbsp; lst1<SPAN class="punctuation token">.</SPAN>append<SPAN class="punctuation token">(</SPAN>r<SPAN class="punctuation token">)</SPAN>
&nbsp;&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;&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>‍</SPAN><SPAN>‍</SPAN></SPAN>

Batch file to run on prompt:

@ECHO OFF 
REM Runs NDVI anomaly scripts
ECHO Backup original files
xcopy D:\Modis250\Originais\*.tif D:\Modis250\Originais_bkp /a /c /i /y /z
PAUSE
ECHO Selecting decade and renaming original files to "ndvi_YYYY_dec_dd.tif" standard.
C:\Python27\ArcGIS10.5\Python.exe D:\Modis250\Python\Script1.py
ECHO NDVI files renamed, hit enter to pre-processing data.
PAUSE
ECHO Projecting, removing spurious data, standardizing and scaling ndvi decade files to 0-100.
C:\Python27\ArcGIS10.5\Python.exe D:\Modis250\Python\Script2.py
ECHO NDVI pre-processed, hit enter to calculate anomalies.
PAUSE
ECHO Calculating decade anomalies.
C:\Python27\ArcGIS10.5\Python.exe D:\Modis250\Python\Script3.py
Del D:\Temp /s
Del D:\NDVI_10dias /s
ECHO Done.
PAUSE‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍<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>

All comments

LuizVianna

Error running at ArcGisPro:

Traceback (most recent call last):
File "<string>", line 40, in <module>
File "<string>", line 13, in main
TypeError: 'NoneType' object is not iterable

in code:

<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; ws_in <SPAN class="operator token">=</SPAN> r<SPAN class="string token">'D:\Modis250\Originals'</SPAN>
&nbsp;&nbsp;&nbsp; ws_out <SPAN class="operator token">=</SPAN> r<SPAN class="string token">'D:\Modis250\NDVI_10dias'</SPAN>

&nbsp;&nbsp;&nbsp; arcpy<SPAN class="punctuation token">.</SPAN>env<SPAN class="punctuation token">.</SPAN>workspace <SPAN class="operator token">=</SPAN> ws_in
&nbsp;&nbsp;&nbsp; lst_ras <SPAN class="operator token">=</SPAN> arcpy<SPAN class="punctuation token">.</SPAN>ListRasters<SPAN class="punctuation token">(</SPAN><SPAN class="string token">"*.tif"</SPAN><SPAN class="punctuation token">)</SPAN>

&nbsp;&nbsp;&nbsp; <SPAN class="keyword token">for</SPAN> ras_name <SPAN class="keyword token">in</SPAN> sorted<SPAN class="punctuation token">(</SPAN>lst_ras<SPAN class="punctuation token">)</SPAN><SPAN class="punctuation token">:</SPAN>
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; in_name_path <SPAN class="operator token">=</SPAN> os<SPAN class="punctuation token">.</SPAN>path<SPAN class="punctuation token">.</SPAN>join<SPAN class="punctuation token">(</SPAN>ws_in<SPAN class="punctuation token">,</SPAN> ras_name<SPAN class="punctuation token">)</SPAN>
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; j <SPAN class="operator token">=</SPAN> int<SPAN class="punctuation token">(</SPAN>ras_name<SPAN class="punctuation token">[</SPAN><SPAN class="number token">15</SPAN><SPAN class="punctuation token">:</SPAN><SPAN class="number token">18</SPAN><SPAN class="punctuation token">]</SPAN><SPAN class="punctuation token">)</SPAN>
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; d <SPAN class="operator token">=</SPAN> GetDecadeFromJulian<SPAN class="punctuation token">(</SPAN>j<SPAN class="punctuation token">)</SPAN>

&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; out_name <SPAN class="operator token">=</SPAN> <SPAN class="string token">"ndvi_{}_dec_{}.tif"</SPAN><SPAN class="punctuation token">.</SPAN>format<SPAN class="punctuation token">(</SPAN>j<SPAN class="punctuation token">,</SPAN> d<SPAN class="punctuation token">)</SPAN>
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; out_name_path <SPAN class="operator token">=</SPAN> os<SPAN class="punctuation token">.</SPAN>path<SPAN class="punctuation token">.</SPAN>join<SPAN class="punctuation token">(</SPAN>ws_out<SPAN class="punctuation token">,</SPAN> out_name<SPAN class="punctuation token">)</SPAN>
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; os<SPAN class="punctuation token">.</SPAN>rename<SPAN class="punctuation token">(</SPAN>in_name_path<SPAN class="punctuation token">,</SPAN> out_name_path<SPAN class="punctuation token">)</SPAN>
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 

<SPAN class="comment token">#def GetDecadeFromJulian(j):</SPAN>
&nbsp;&nbsp;&nbsp; <SPAN class="comment token">#d =&nbsp; (j+1 - ((j-1) % 10)) / 10 + 1</SPAN>
&nbsp;&nbsp;&nbsp; <SPAN class="comment token">#return str(d) if d < 37 else '36'</SPAN>
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;
<SPAN class="keyword token">def</SPAN> <SPAN class="token function">GetDecadeFromJulian</SPAN><SPAN class="punctuation token">(</SPAN>year<SPAN class="punctuation token">,</SPAN> j<SPAN class="punctuation token">)</SPAN><SPAN class="punctuation token">:</SPAN>
&nbsp;&nbsp;&nbsp; <SPAN class="keyword token">import</SPAN> datetime
&nbsp;&nbsp;&nbsp; date <SPAN class="operator token">=</SPAN> datetime<SPAN class="punctuation token">.</SPAN>datetime<SPAN class="punctuation token">(</SPAN>year<SPAN class="punctuation token">,</SPAN> <SPAN class="number token">1</SPAN><SPAN class="punctuation token">,</SPAN> <SPAN class="number token">1</SPAN><SPAN class="punctuation token">)</SPAN> <SPAN class="operator token">+</SPAN> datetime<SPAN class="punctuation token">.</SPAN>timedelta<SPAN class="punctuation token">(</SPAN>j<SPAN class="punctuation token">)</SPAN>
&nbsp;&nbsp;&nbsp; m <SPAN class="operator token">=</SPAN> <SPAN class="punctuation token">(</SPAN>date<SPAN class="punctuation token">.</SPAN>month <SPAN class="operator token">-</SPAN> <SPAN class="number token">1</SPAN><SPAN class="punctuation token">)</SPAN> <SPAN class="operator token">*</SPAN> <SPAN class="number token">3</SPAN>
&nbsp;&nbsp;&nbsp; <SPAN class="keyword token">if</SPAN> date<SPAN class="punctuation token">.</SPAN>day <SPAN class="operator token"><=</SPAN> <SPAN class="number token">10</SPAN><SPAN class="punctuation token">:</SPAN>
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; d <SPAN class="operator token">=</SPAN> <SPAN class="number token">1</SPAN>
&nbsp;&nbsp;&nbsp; <SPAN class="keyword token">elif</SPAN> date<SPAN class="punctuation token">.</SPAN>day <SPAN class="operator token"><=</SPAN> <SPAN class="number token">20</SPAN><SPAN class="punctuation token">:</SPAN>
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; d <SPAN class="operator token">=</SPAN> <SPAN class="number token">2</SPAN>
&nbsp;&nbsp;&nbsp; <SPAN class="keyword token">else</SPAN><SPAN class="punctuation token">:</SPAN>
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; d <SPAN class="operator token">=</SPAN> <SPAN class="number token">3</SPAN>
&nbsp;&nbsp;&nbsp; <SPAN class="keyword token">return</SPAN> m <SPAN class="operator token">+</SPAN> d

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

XanderBakker

This would require changing the function to get the decade like this:

<SPAN class="keyword token">def</SPAN> <SPAN class="token function">GetDecadeFromJulian</SPAN><SPAN class="punctuation token">(</SPAN>year<SPAN class="punctuation token">,</SPAN> j<SPAN class="punctuation token">)</SPAN><SPAN class="punctuation token">:</SPAN>
&nbsp;&nbsp;&nbsp; <SPAN class="keyword token">import</SPAN> datetime
&nbsp;&nbsp;&nbsp; date <SPAN class="operator token">=</SPAN> datetime<SPAN class="punctuation token">.</SPAN>datetime<SPAN class="punctuation token">(</SPAN>year<SPAN class="punctuation token">,</SPAN> <SPAN class="number token">1</SPAN><SPAN class="punctuation token">,</SPAN> <SPAN class="number token">1</SPAN><SPAN class="punctuation token">)</SPAN> <SPAN class="operator token">+</SPAN> datetime<SPAN class="punctuation token">.</SPAN>timedelta<SPAN class="punctuation token">(</SPAN>j<SPAN class="punctuation token">)</SPAN>
&nbsp;&nbsp;&nbsp; m <SPAN class="operator token">=</SPAN> <SPAN class="punctuation token">(</SPAN>date<SPAN class="punctuation token">.</SPAN>month <SPAN class="operator token">-</SPAN> <SPAN class="number token">1</SPAN><SPAN class="punctuation token">)</SPAN> <SPAN class="operator token">*</SPAN> <SPAN class="number token">3</SPAN>
&nbsp;&nbsp;&nbsp; <SPAN class="keyword token">if</SPAN> date<SPAN class="punctuation token">.</SPAN>day <SPAN class="operator token"><=</SPAN> <SPAN class="number token">10</SPAN><SPAN class="punctuation token">:</SPAN>
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; d <SPAN class="operator token">=</SPAN> <SPAN class="number token">1</SPAN>
&nbsp;&nbsp;&nbsp; <SPAN class="keyword token">elif</SPAN> date<SPAN class="punctuation token">.</SPAN>day <SPAN class="operator token"><=</SPAN> <SPAN class="number token">20</SPAN><SPAN class="punctuation token">:</SPAN>
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; d <SPAN class="operator token">=</SPAN> <SPAN class="number token">2</SPAN>
&nbsp;&nbsp;&nbsp; <SPAN class="keyword token">else</SPAN><SPAN class="punctuation token">:</SPAN>
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; d <SPAN class="operator token">=</SPAN> <SPAN class="number token">3</SPAN>
&nbsp;&nbsp;&nbsp; <SPAN class="keyword token">return</SPAN> m <SPAN class="operator token">+</SPAN> d<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>

Both year and Julian (j) should be provided as numeric value.

LuizVianna

Question: how often do you plan to run this script? Does the input data change?<SPAN class="line-numbers-rows"><SPAN>‍</SPAN></SPAN>

The first process will be runned one time for all the past years files (from 2000 to today). After that, the script will be runned every ten days period (decade). The new input (original file "MODIS.ndvi.YYYYjjj.yL6000.BOKU.tif") will be pre-processed and the result will be added to the list of standardized and scaled (0-100) raster datasets, to calculate the anomaly for the last decade in current year. It is a continous process.

As the name of the original file have the "YYYYjjj" format, is it possible to rename the files considering the leap years, according to the table below?

Decade	Date	julian day in common year	julian day in leap year
1	01/01 to 01/10	10	10
2	01/11 to 01/20	20	20
3	01/21 to 01/31	31	31
4	02/01 to 02/10	41	41
5	02/11 to 02/20	51	51
6	02/21 to 02/28-29	59	60
7	03/01 to 03/10	69	70
8	03/11 to 03/20	79	80
9	03/21 to 03/31	90	91
10	04/01 to 04/10	100	101
11	04/11 to 04/20	110	111
12	04/21 to 04/30	120	121
13	05/01 to 05/10	130	131
14	05/11 to 05/20	140	141
15	05/21 to 05/31	151	152
16	06/01 to 06/10	161	162
17	06/11 to 06/20	171	172
18	06/21 to 06/30	181	182
19	07/01 to 07/10	191	192
20	07/11 to 07/20	201	202
21	07/21 to 07/31	212	213
22	08/01 to 08/10	222	223
23	08/11 to 08/20	232	233
24	08/21 to 08/31	243	244
25	09/01 to 09/10	253	254
26	09/11 to 09/20	263	264
27	09/21 to 09/30	273	274
28	10/01 to 10/10	283	284
29	10/11 to 10/20	293	294
30	10/21 to 10/31	304	305
31	11/01 to 11/10	314	315
32	11/11 to 11/20	324	325
33	11/21 to 11/30	334	335
34	12/01 to 12/10	344	345
35	12/11 to 12/20	354	355
36	12/21 to 12/31	365	366

LuizVianna

The script 1 worked, however the output name is wrong: The names should be "ndvi_yyyy_dec_d.tif" where "yyyy" is the year (from the input file name (MODIS.ndvi.YYYYjjj.yL6000.BOKU.tif ) and "d" the decade. Why the ".2" after the decade number in the outputs? If there are more then one input file for decade, de code shoud overwrite the first.

XanderBakker

Let's try and do this in parts.

Question 1: What is the best option: to run separated scripts or create a unique code?

You can work with a single script that does all 3 steps, but it might be better to break it up in parts(as you did). This way you can check the intermediate results and run parts separately and also reuse parts for other workflows.

Question 2: How to translate the model builder to python script (Python script 2)?

You can export the Model to a Python script, but that probably results in some really ugly code. It is better to rewrite the model into a new script. This also allows for some additional checks that might not be easy to do inside a model.

Question 3: If the best option is to run separated, how to make it in a logical sequence (script 1 -> script 2 -> script 3)

Depends. As I mentioned above, it may be easier to run parts separately and do some (manual) checks between the steps to be sure things go right. When you have tested thoroughly, you can call the scripts from a "master" script and run the entire process. Question: how often do you plan to run this script? Does the input data change?

I think that the first script could be rewritten like this (untested):

<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; ws_in <SPAN class="operator token">=</SPAN> <SPAN class="string token">'O:\Modis250\Originals'</SPAN>
&nbsp;&nbsp;&nbsp; ws_out <SPAN class="operator token">=</SPAN> <SPAN class="string token">'O:\Modis250\NDVI_10dias'</SPAN>

&nbsp;&nbsp;&nbsp; arcpy<SPAN class="punctuation token">.</SPAN>env<SPAN class="punctuation token">.</SPAN>workspace <SPAN class="operator token">=</SPAN> ws_in
&nbsp;&nbsp;&nbsp; lst_ras <SPAN class="operator token">=</SPAN> arcpy<SPAN class="punctuation token">.</SPAN>ListRasters<SPAN class="punctuation token">(</SPAN><SPAN class="string token">"*.tif"</SPAN><SPAN class="punctuation token">)</SPAN>

&nbsp;&nbsp;&nbsp; <SPAN class="keyword token">for</SPAN> ras_name <SPAN class="keyword token">in</SPAN> sorted<SPAN class="punctuation token">(</SPAN>lst_ras<SPAN class="punctuation token">)</SPAN><SPAN class="punctuation token">:</SPAN>
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; in_name_path <SPAN class="operator token">=</SPAN> os<SPAN class="punctuation token">.</SPAN>path<SPAN class="punctuation token">.</SPAN>join<SPAN class="punctuation token">(</SPAN>ws_in<SPAN class="punctuation token">,</SPAN> ras_name<SPAN class="punctuation token">)</SPAN>
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; j <SPAN class="operator token">=</SPAN> int<SPAN class="punctuation token">(</SPAN>ras_name<SPAN class="punctuation token">[</SPAN><SPAN class="number token">15</SPAN><SPAN class="punctuation token">:</SPAN><SPAN class="number token">18</SPAN><SPAN class="punctuation token">]</SPAN><SPAN class="punctuation token">)</SPAN>
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; d <SPAN class="operator token">=</SPAN> GetDecadeFromJulian<SPAN class="punctuation token">(</SPAN>j<SPAN class="punctuation token">)</SPAN>

&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; out_name <SPAN class="operator token">=</SPAN> <SPAN class="string token">"ndvi_{}_dec_{}.tif"</SPAN><SPAN class="punctuation token">.</SPAN>format<SPAN class="punctuation token">(</SPAN>j<SPAN class="punctuation token">,</SPAN> d<SPAN class="punctuation token">)</SPAN>
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; out_name_path <SPAN class="operator token">=</SPAN> os<SPAN class="punctuation token">.</SPAN>path<SPAN class="punctuation token">.</SPAN>join<SPAN class="punctuation token">(</SPAN>ws_out<SPAN class="punctuation token">,</SPAN> out_name<SPAN class="punctuation token">)</SPAN>
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; os<SPAN class="punctuation token">.</SPAN>rename<SPAN class="punctuation token">(</SPAN>in_name_path<SPAN class="punctuation token">,</SPAN> out_name_path<SPAN class="punctuation token">)</SPAN>

<SPAN class="keyword token">def</SPAN> <SPAN class="token function">GetDecadeFromJulian</SPAN><SPAN class="punctuation token">(</SPAN>j<SPAN class="punctuation token">)</SPAN><SPAN class="punctuation token">:</SPAN>
&nbsp;&nbsp;&nbsp; d <SPAN class="operator token">=</SPAN>&nbsp; <SPAN class="punctuation token">(</SPAN>j<SPAN class="operator token">+</SPAN><SPAN class="number token">1</SPAN> <SPAN class="operator token">-</SPAN> <SPAN class="punctuation token">(</SPAN><SPAN class="punctuation token">(</SPAN>j<SPAN class="number token">-1</SPAN><SPAN class="punctuation token">)</SPAN> <SPAN class="operator token">%</SPAN> <SPAN class="number token">10</SPAN><SPAN class="punctuation token">)</SPAN><SPAN class="punctuation token">)</SPAN> <SPAN class="operator token">/</SPAN> <SPAN class="number token">10</SPAN> <SPAN class="operator token">+</SPAN> <SPAN class="number token">1</SPAN>
&nbsp;&nbsp;&nbsp; <SPAN class="keyword token">return</SPAN> str<SPAN class="punctuation token">(</SPAN>d<SPAN class="punctuation token">)</SPAN> <SPAN class="keyword token">if</SPAN> d <SPAN class="operator token"><</SPAN> <SPAN class="number token">37</SPAN> <SPAN class="keyword token">else</SPAN> <SPAN class="string token">'36'</SPAN>

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

Rejoindre la conversation

Les membres connectés peuvent publier, suivre les mises à jour, et plus encore. Nouveau ici ? Inscrivez-vous gratuitement.

Liens rapides

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