Python：arcpy批处理

按照一个矢量文件批量裁剪多个raster数据：

import os
import arcpy
import glob

filepath=r'F:RasterData'
Clppath=r'F:Clp_RasterData'
mask=r'F:SHPrec.shp'

#筛选tif文件
files=glob.glob(os.path.join(filepath,'*.tif'))
for f in files:
	outname=os.path.join(Clppath, os.path.basename(f).split(".")[0] + ".tif")
	out_extract=arcpy.sa.ExtractByMask(f, mask)
	out_extract.save(outname)
print('OVER')	

按照多个矢量文件裁剪一个raster数据：
主要区别在于文件的搜索

import os
import arcpy
import glob

filepath=r'F:ratserdataDEM.tif'
shppath=r'F:SHP'
respath=r'F:SplitDEM'

shps=glob.glob(os.path.join(shppath, "*.shp"))
for s in shps:
	outname = os.path.join(respath, os.path.basename(s).split(".")[0] + ".tif")
	out_extract = arcpy.sa.ExtractByMask(filepath, s)
	out_extract.save(outname)
print('OVER')

将所有图像都拼接到第一个图像中去：

import os
import arcpy
import glob

filepath=r'F:ratserdata'
files=glob.glob(os.path.join(filepath, "*.tif"))
r0=files[0]

for f in files[1:]:
	# f向r0拼接，留在r0上
	#拼接方法还有：LAST，FIRST，MAXIMUM，MINIMUM
	arcpy.Mosaic_management(f,r0,'BLEND')
print('OVER')