arcpy遥感影像批量判断满足条件的栅格

实现内容：批量判断某个文件夹下的所有影像，依次遍历每个栅格值，大于25并且小于100的栅格则输出为1，否之输出为0。最后每个影像对应一张0-1的栅格图。（已测试结果正确性）

# coding=utf-8
import arcpy
from arcpy import env
from arcpy.sa import *


# Set environment settings
# 输入工作空间文件夹（即存放需批处理tif影像的文件夹）
env.workspace = "D:/毕业论文/3.arcmap数据处理/暴雨数量/tif"
# 注意此处‘/’的方向


# Set local variables
# 遍历工作空间中的tif格式数据
rasters = arcpy.ListRasters("*", "tif")

# Check out the ArcGIS Spatial Analyst extension license
arcpy.CheckOutExtension("Spatial")
# 循环判断每个影像是否满足条件
for raster in rasters:
      outCellStatistics = Con(raster, 1, 0, "value <= 250 and value > 100 ")
      outCellStatistics.save("D:/d100_250/"+raster)