ArcGIS python calcula o valor médio de múltiplos dados raster em uma longa série de tempo

Normalmente, precisamos calcular a média de vários rasters. Por exemplo, se somarmos o valor NDVI de cada mês do ano e dividirmos por 12, esperaremos até o NDVI médio mensal. Embora esse processo possa ser implementado em uma calculadora raster, Quando a série temporal é longa, é mais problemática.Neste momento, o código python é a melhor escolha.

A figura abaixo mostra o princípio de adicionar dados raster, ou seja, adicionar raster correspondentes para gerar novos dados raster. O valor médio deve ser dividido pelo número de grades.

Escreva o seguinte código em um script autônomo:

Código de valor médio Python como um presente:

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