import arcpy
from arcpy import env
from arcpy.sa import *
env.workspace = "H:/yiruzhou/py2/BEIFEN"
rasterlist = arcpy.ListRasters("*", "tif")
##声明一个数组
#a = ["Anshun","Bijie","Duyun","Dushan","Guiyang","Libo","Liping","Luodian","Panxian","Puan","Qianxi",“Renhuai”,“Rongjiang”,“Sinan”,“Songtao”,“Tianzhu”,“Tongren”,“Tongzi”,“Wangmo”,“Xingren”,“Zhengan”,“Zhijin”,“Zunyi”,“sansui”]
a=["Panxian","Puan","Qianxi","Renhuai","Rongjiang","Sinan","Songtao","Tianzhu","Tongren","Tongzi","Wangmo","Xingren","Zhengan","Zhijin","Zunyi","sansui"]
##遍历之
for sitename in a:
for raster in rasterlist:
print str(sitename+raster)
outR = "H:/yiruzhou/py2/BEIFEN/"+raster
outExtractByMask = ExtractByMask(outR, sitename+'.shp')
outExtractByMask.save("H:/yiruzhou/py2/"+sitename+"/LTD/"+raster)
print("OK!")
网友评论