ArcGis Python模块批量裁剪代码

import arcpy
… from arcpy import env
… from arcpy.sa import *
… env.workspace = “D:\PVproject4\HC_s1\S1”
… for i in range(50):
… raster1 = ‘DEMno’
… raster2 = str(i+1)
… raster3 = ‘_13.tif.tif’
… raster = ‘%s%s%s’%(raster1,raster2,raster3)
… outExtractByMask = ExtractByMask(raster, “D:\PVproject4\HC_s1\S1\T3t.shp”)
… savepath1 = “D:\PVproject4\HC_s1\T3\ERno”
… savepath = “%s%s”%(savepath1,raster2)
… outExtractByMask.save(savepath)