ArcGIS python calcule la valeur moyenne de plusieurs données raster dans une longue série chronologique

En général, nous devons calculer la moyenne de plusieurs rasters. Par exemple, si nous additionnons la valeur NDVI de chaque mois de l'année et la divisons par 12, nous attendrons la moyenne mensuelle du NDVI. Bien que ce processus puisse être implémenté dans un calculateur raster, Lorsque la série chronologique est longue, c'est plus gênant. Pour le moment, le code Python est le meilleur choix.

La figure suivante illustre le principe de l'ajout de données raster, c'est-à-dire l'ajout de rasters correspondants pour générer de nouvelles données raster. La valeur moyenne doit être divisée par le nombre de grilles.

Écrivez le code suivant dans un script autonome:

Code de valeur moyenne Python en cadeau:

import arcpy
from arcpy.sa import *
arcpy.CheckOutExtension("spatial")
arcpy.gp.overwriteOutput=1
 
#custom
arcpy.env.workspace="G:\\Phenology of 30 Years\\GIMMS 3g\\15Length\\1Length\\"
#custom
outpath="G:\\Phenology of 30 Years\\GIMMS 3g\\15Length\\2mean_len\\"
#custom
outfilename="mean";n=30;
 
Sum=0
 
files=arcpy.ListRasters()
 
for file in files:
    Sum=Sum+Raster(file)
(Sum/n).save(outpath+outfilename)
 
print("Done,please close")