需求
為每個(gè)輸入像元位置計(jì)算其周圍指定鄰域內(nèi)的值的統(tǒng)計(jì)數(shù)據(jù)；
批量得到tif柵格影像每個(gè)像元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!")
