topic Re: Python script to calculate NDVI anomalies. in ArcGIS Pro Questions

Python script to calculate NDVI anomalies.

LuizVianna — Thu, 07 Jun 2018 19:11:38 GMT

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

This is a python script to calculate NDVI anomalies.

Thanks.

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

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

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

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

Re: Python script to calculate NDVI anomalies.

DanPatterson_Retired — Sat, 11 Dec 2021 06:16:42 GMT

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

And the error message?

The last line?

#print ("Decade "+d+" of "+y+" ok")

# ---- maybe try this instead, it handles type conversion to string

print("Decade {} of {} ok".format(d, y))‍‍‍‍‍

Re: Python script to calculate NDVI anomalies.

LuizVianna — Fri, 08 Jun 2018 12:33:48 GMT

Thank you Dan!

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

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

The errors are:

For the arcpy.gp.Minus_3d approach:

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

For the "NDVIAnom = NDVINormY - NDVIMean

NDVIAnom.save(NDVIAnomOut)" approach

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

Re: Python script to calculate NDVI anomalies.

JoeBorgione — Sat, 11 Dec 2021 06:16:45 GMT

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

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

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

#Define workspaces
arcpy.env.workspace="D:\\Modis250\\NDVI_10dias_norm" #Workspace with decendial (ten days period) NDVI

#Filenames standard is "norm_ndvi_yyyy_dec_d.tif" where "yyyy" is the year and "d" is the ten days period: d=1 = january 1-10; d=2 = january 11-20; ...

for n in range (1,37):
    d=str(1)
    y = str(2004)
    NDVINormAll=arcpy.ListRasters("norm_ndvi_*_"+d+".tif", "tif")
    #List of decendial NDVI for all years to calculate the dcendial mean

    NDVINormY=arcpy.ListRasters("norm_ndvi_"+y+"_dec_"+d+".tif", "tif")
    #List of decendial NDVI rasters for the last year (y)

    NDVIMeanOut = "D:\\Modis250\\NDVI_10dias_mean\\"+"mean_dec_"+d+".tif" 
    #Define the decendial NDVI mean output (could be temporary)

    NDVIMean = arcpy.gp.CellStatistics_sa(NDVINormAll, NDVIMeanOut,"MEAN", "NODATA")
    #Calculate the decendial NDVI mean (could be temporary)

    NDVIAnomOut = "D:\\Modis250\\NDVI_10dias_anom\\"+"anom_"+y+"_dec_"+d+".tif" 
    Define the last year decendial NDVI anomaly output

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

    print ("Decade "+d+" of "+y+" ok")‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍

Re: Python script to calculate NDVI anomalies.

LuizVianna — Fri, 08 Jun 2018 14:40:43 GMT

Thank you Joe!

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

I changed lines 24 and line 30 to:

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

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

and when I run the code, after line 33

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

the error is:

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

If I try:

NDVIAnom = NDVINormY - NDVIMean
NDVIAnom.save(NDVIAnomOut)

the error is

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

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

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

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

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

Re: Python script to calculate NDVI anomalies.

XanderBakker — Fri, 08 Jun 2018 14:55:07 GMT

Hi Joe Borgione ,

What does not make sense is:

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

Re: Python script to calculate NDVI anomalies.

JoeBorgione — Sat, 11 Dec 2021 06:16:48 GMT

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

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

Xander Bakker‌

Luiz Vianna‌

for n in range(1,37):
  print '{}'.format(n)‍‍

Re: Python script to calculate NDVI anomalies.

XanderBakker — Sat, 11 Dec 2021 06:16:50 GMT

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

Look at the following structure:

year_range = range(2004, 2011) # to get the values from 2004 until 2010
for year in year_range:
    print year

    # for each year:
    # - create list of all NDVI for current year: lst_ndvi_ras
    # - calculate the mean of all NDVI for current year: mean_ndvi_year

    for ndvi_ras in lst_ndvi_ras:
        # calculate anomaly and save result with unique name‍‍‍‍‍‍‍‍‍‍

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

Re: Python script to calculate NDVI anomalies.

LuizVianna — Fri, 08 Jun 2018 16:34:58 GMT

Hi Xander!

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

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

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

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

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

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

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

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

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

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

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

Thank you!

Re: Python script to calculate NDVI anomalies.

LuizVianna — Fri, 08 Jun 2018 16:40:40 GMT

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

Re: Python script to calculate NDVI anomalies.

XanderBakker — Sat, 11 Dec 2021 06:16:53 GMT

Maybe something like the untested code below could work:

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

    # get list of years
    decades = GetListOfDecades(lst_ndvi_ras)

    # loop through each year
    for decade in decades:
        # Get rasters for given year
        lst_ndvi_decade = GetListOfRasterForGivenDecade(lst_ndvi_ras, decade)

        # calculate mean raster for decade
        mean_ras = arcpy.sa.CellStatistics(lst_ndvi_year, "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)
        mean_ras.save(out_name_path_mean)

        # loop through each raster in decade
        for ndvi_ras in lst_ndvi_ras:
            # calculate anom for each raster in decade
            # anom_2002_dec_1.tif = norm_ndvi_2002_dec_1.tif - mean_dec_1_.tif;
            minus_ras = arcpy.sa.Minus(ndvi_ras, 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)


    # 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()‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍

Re: Python script to calculate NDVI anomalies.

LuizVianna — Fri, 08 Jun 2018 18:00:15 GMT

Thank you, Xander!

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

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

I expected for 'D:\Modis250\NDVI_10dias_mean'

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

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

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

Re: Python script to calculate NDVI anomalies.

XanderBakker — Fri, 08 Jun 2018 18:07:04 GMT

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

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

Re: Python script to calculate NDVI anomalies.

XanderBakker — Sat, 11 Dec 2021 06:16:56 GMT

Something like this:

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

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

# get list of years
decades = GetListOfDecades(lst_ndvi_ras)

# loop through each year
for decade in decades:
    # Get rasters for given year
    lst_ndvi_decade = GetListOfRasterForGivenDecade(lst_ndvi_ras, decade)

    # calculate mean raster for decade
    mean_ras = arcpy.sa.CellStatistics(lst_ndvi_year, "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)
    mean_ras.save(out_name_path_mean)

    # loop through each raster in decade
    for ndvi_ras in lst_ndvi_ras:
        # calculate anom for each raster in decade
        # anom_2002_dec_1.tif = norm_ndvi_2002_dec_1.tif - mean_dec_1_.tif;
        minus_ras = arcpy.sa.Minus(ndvi_ras, 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)


# return SA license
arcpy.CheckInExtension("spatial")‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍

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

Re: Python script to calculate NDVI anomalies.

LuizVianna — Fri, 08 Jun 2018 18:45:23 GMT

Xander,

Is "lst_ndvi_year" in line 27 correct? Is´nt it "lst_ndvi_decade"?

Re: Python script to calculate NDVI anomalies.

XanderBakker — Fri, 08 Jun 2018 19:02:59 GMT

Good catch , you're absolutely right (that's why I don't like creating code without data to test on).

Additionally on line 36 it should be lst_ndvi_decade instead of lst_ndvi_ras

Re: Python script to calculate NDVI anomalies.

LuizVianna — Fri, 08 Jun 2018 19:19:49 GMT

I corrected the code and executed. The error now is:

Traceback (most recent call last):
File "<string>", line 65, in <module>
File "<string>", line 39, in main
File "c:\program files\arcgis\pro\Resources\arcpy\arcpy\sa\Functions.py", line 4271, in Minus
in_raster_or_constant2)
File "c:\program files\arcgis\pro\Resources\arcpy\arcpy\sa\Utils.py", line 53, in swapper
result = wrapper(*args, **kwargs)
File "c:\program files\arcgis\pro\Resources\arcpy\arcpy\sa\Functions.py", line 4268, in Wrapper
["Minus", in_raster_or_constant1, in_raster_or_constant2])
RuntimeError: ERROR 000732: Input Raster: Dataset norm_ndvi_2002_dec_31.tif does not exist or is not supported

Re: Python script to calculate NDVI anomalies.

XanderBakker — Fri, 08 Jun 2018 20:31:16 GMT

Did it create anything (for other decades) or is this the first decade it processed? Was the mean raster created correctly? Does the mentioned raster exist and is it valid (does it load into Desktop)?

Re: Python script to calculate NDVI anomalies.

LuizVianna — Tue, 12 Jun 2018 17:44:18 GMT

It created one tif of mean and crashed. The error is:

I've made the same folder structure and sample files available for a period of 3 years in the GoogleDrive (https://drive.google.com/file/d/1G4baOl1NP7G1rEfBpu9SZr4cxEOd1-Jo/view?usp=sharing).

Best.

Re: Python script to calculate NDVI anomalies.

XanderBakker — Sat, 11 Dec 2021 06:16:59 GMT

Thanks for sharing the data. I made some changes to the code and that seems to be working:

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()
    print "lst_ndvi_ras:", lst_ndvi_ras

    # get list of years
    decades = GetListOfDecades(lst_ndvi_ras)
    print "decades:", 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:", lst_ndvi_decade

        # calculate mean raster for decade
        mean_ras = arcpy.sa.CellStatistics(lst_ndvi_decade, "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:", 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:", ndvi_ras
            minus_ras = arcpy.sa.Minus(ndvi_ras, 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 "out_name_path_anom:", 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()‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍

Here an example "anom" raster (anom_2002_dec_1.tif):