批量polygon切一个栅格数据

有一个文件夹里有N个POLYGON，其实是将接图表里的图幅split出来，需要用这些polygon去切一个栅格数据。当然有种最笨的方法，就是依次用extract  by mask去裁减那个栅格数据。或者在arctoolboxt的extract  by mask工具选右键batch，还是很慢。

 借助了python，参考了arcgis的帮助，自己写了段小代码，居然实现了，哈哈。

import arcgisscripting
import os
gp = arcgisscripting.create()
gp.CheckOutExtension ("Spatial")

 

try:
  InRaster="F:\\temp\\aa"
  filenames = os.listdir("F:/temp/polygon")
  for name in filenames:
     filenames[filenames.index(name)]=name[:-4]
     filename=set(filenames) 
  for name1 in filename:
       if len(name1)>0:
           InFuture="F:\\temp\\polygon\\"+name1+".shp"
           utRaster="F:\\temp\\"+name1
           print InFuture
           gp.ExtractByMask_sa(InRaster, InFuture, OutRaster)
except:
 print gp.GetMessages()

来自：http://www.gisall.com/?uid-13588-action-viewspace-itemid-2260