arcgis循环裁剪shp文件

2020-07-15  本文已影响0人  山有木兮木有刺

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!")

上一篇 下一篇

猜你喜欢

热点阅读