python批处理系列代码-栅格数据按掩膜处理_python按阈值掩膜栅格图像

import arcpy
from arcpy import env
from arcpy.sa import *
import os
import os.path
import sys

#arcpy.env.workspace = "D:\\TIF\\IMERGD\\Y"
rootdir = 'D:\\TIF\\IMERGD\\Y'

filenames = os.listdir(rootdir)
for filename in filenames:
        if os.path.splitext(filename)[1] == '.tif':
            inRaster =rootdir + os.sep + filename
            outname = filename.replace('.tif','_cjmm')
            inMaskData = "D:\\TIF\\HB_PROVINCE.shp"
            arcpy.CheckOutExtension("Spatial")
            outExtractByMask = ExtractByMask(inRaster, inMaskData)
            print "ok"
            outExtractByMask.save("D:\\test\\"+outname+'.tif')
            print outname