栅格重分类
from arcpy.sa import Con
import arcpy
from arcpy.sa import *
arcpy.CheckOutExtension('spatial')
arcpy.env.workspace = r"D:\Factor\MOD11A2_ST\xiangytjdata_tif"
output_path = "D:\Factor\MOD11A2_ST\Reclassify_xytj_MOD11A2"
rasterlist = arcpy.ListRasters("*","tif")
for raster in rasterlist:
inRaster = arcpy.Raster(raster)
out = output_path + str(raster)
outCon =Con(inRaster < 289,1,
Con((inRaster > 289) & (inRaster <= 291), 2,
Con((inRaster > 291) & (inRaster <= 293), 3,
Con((inRaster > 293) & (inRaster <= 295), 4,
Con((inRaster > 295) & (inRaster <= 297), 5,
Con((inRaster > 297) & (inRaster <= 299), 6,
Con((inRaster > 299) & (inRaster <= 301), 7,
Con((inRaster > 301) & (inRaster <= 303), 8,
Con((inRaster > 303) & (inRaster <= 305), 9,
Con((inRaster > 305) & (inRaster <= 307), 10,
Con((inRaster > 307), 11,)))))))))))
outCon.save(out)
print ('SUCCESS! OK!')
