ArcGIS pythonは、長い時系列の複数のラスターデータの平均値を計算します

通常、複数のラスターを平均化する必要があります。たとえば、1年の各月のNDVI値を合計して12で割ると、月平均NDVIまで待機します。このプロセスはラスター計算機で実装できますが、時系列が長いと面倒ですが、現時点ではPythonコードが最適です。

次の図は、ラスターデータを追加する、つまり、対応するラスターを追加して新しいラスターデータを生成する原理を示しています。平均値をグリッド数で割る必要があります。

スタンドアロンスクリプトで次のコードを記述します。

贈り物としてのPythonの平均値コード：

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