批量去除栅格图像异常值

#coding:utf-8
import arcpy
from arcpy import env
from arcpy.sa import *
import sys
reload(sys)
sys.setdefaultencoding('utf8')

arcpy.CheckOutExtension("ImageAnalyst")  # 检查许可
arcpy.CheckOutExtension("spatial")  # 检查许可
arcpy.env.overwriteOutput=True
env.workspace = r"I:\硕士论文数据\ET数据\ET裁剪_MOD16A3GF"  # 数据存放目录

outfile = r'I:\大论文数据处理\ET_年和季_去掉Nodata\year'  # 输出数据存放路径

for file1 in arcpy.ListFiles("*ET_500m*.tif"):  # 循环列表下的tif文件
    s = file1   
    C_outfile = outfile + '\\' + s  # 输出的文件名，完整路径 
    inRaster = Raster(s)

    outSetNull = SetNull(inRaster>65500,inRaster)*0.1
    outSetNull = SetNull(outSetNull==181,outSetNull)
    #以上两条件都为去除条件
    outSetNull.save(C_outfile)

     # 保存数据

print ('over')