Arcpy栅格掩膜

遥感影像太多，要获取影像一部分研究区就要进行栅格裁剪，arcpy代码如下，可以成功运行，做个记录

小提示：一次不要运行太大量的数据，要不然容易崩。

import arcpy
from arcpy import env
import os
path1= "F:\TM" #遥感影像文件夹所在路径，很多景影像下载后在不同的文件夹里，这个路径是文件夹的路径              
fileList1 = os.listdir(path1)   #获取path1路径下文件夹的名字     
path2="F:\LayerClip\TMClip"    #裁剪后存放的路径
fileList2 = os.listdir(path2)  #获取path2路径下文件夹的名字 
i=0
for f in fileList1:
     inpath=path1+'\\'+fileList1[i]    #遥感影像的路径
     env.workspace =inpath
     rasterList = arcpy.ListRasters("*","TIF")
     mask = "F:\study_area\LT5_studyarea.shp"   
     output_path=path2+'\\'+fileList2[i]+'\\d'            
     for raster in rasterList:
          out = output_path + raster[0:-4] + "_NL" + ".tif" #对生成文件进行命名
          arcpy.Clip_management(raster,"#",out,mask,"0","ClippingGeometry")
          print out+" clip has done!"
     i=i+1