基于Python实现批量化焦点统计

需求

为每个输入像元位置计算其周围指定邻域内的值的统计数据；
批量得到tif栅格影像每个像元3*3矩阵窗口的平均值。

代码

import arcpy
from arcpy import env
from arcpy.sa import *
env.workspace = "D:/GPP/tem_opt/result/test"
rasterlist = arcpy.ListRasters("*", "tif")   
for raster in rasterlist:
   print str(raster)
   out = "D:/GPP/tem_opt/result/test2/"+raster
   arcpy.gp.FocalStatistics_sa(raster,out, "Rectangle 3 3 CELL", 
                               "MEAN", "DATA");
print("OK!")

官方说明文档

焦点统计