The monthly maximum composition (MVC) data for 12 months of the year is placed in "F:\\vegetation Change\\data\\gimms data\\1mvc\\" with the data name format mvc_198801,mvc_198802 .... mvc_198812. The processing year is 1981-2006 and the code is:
Import arcpyarcpy. Checkoutextension ("spatial") for I in Range (1981,2007): a1 = "F:\\vegetation Change\\data\\gimms data\\1mvc\\" + "MVC_" +s TR (i) + "a2" = "F:\\vegetation change\\data\\gimms data\\1mvc\\" + "mvc_" +str (i) + "a3 =" F:\\vegetation Change\\da Ta\\gimms data\\1mvc\\ "+" mvc_ "+str (i) +" a4 = "f:\\vegetation change\\data\\gimms data\\1mvc\\" + "mvc_" +str (i) + "04" A5 = "F:\\vegetation change\\data\\gimms data\\1mvc\\" + "mvc_" +str (i) + "A6 =" F:\\vegetation Change\\data\\gimms Da ta\\1mvc\\ "+" mvc_ "+str (i) +" a7 = "F:\\vegetation change\\data\\gimms data\\1mvc\\" + "mvc_" +str (i) + "A8 =" F:\\v Egetation Change\\data\\gimms data\\1mvc\\ "+" mvc_ "+str (i) +" A9 = "F:\\vegetation Change\\data\\gimms Data\\1MVC\\" + "Mvc_" +str (i) + "a10" = "f:\\vegetation change\\data\\gimms data\\1mvc\\" + "mvc_" +str (i) + "ten" A11 = "f:\\vegetation Change\\data\\gimms data\\1mvc\\ "+" mvc_ "+str (i) +" one "A12 =" f:\\vegetation change\\data\\gimms data\\1mvc\\ "+" mvc_ "+ STR (i) + "12" out = "F:\\vegetation Change\\data\\gimms data\\2max\\max_" +str (i) Arcpy.gp.CellStatistics_sa ((a1,a2,a3,a4,a5,a6,a7 , A8,a9,a10,a11,a12), out, "MAXIMUM", "DATA") print (str (i) + "have done")
ArcGIS Python implements the annual maximum for MODIS ndvi batch