# coding=utf-8
import sys
reload(sys)
sys.setdefaultencoding('utf-8')
import arcpy
import glob
import os
arcpy.CheckOutExtension('Spatial')
#指定先前拼接后的遥感影像所在目录
inws = r"D:\data\MOD06\cloudfraction\20185"
# 指定裁剪后的影响存放目录
outws = r"D:\data\output\COF"
# 指定shp范围边界文件,即目标区域的边界
mask = r"D:\data\国界\国界.shp"
# 利用glob包,将inws下的所有tif文件读存放到rasters中
rasters = glob.glob(os.path.join(inws, "*.tif"))
# 循环rasters中的所有影像,进行按掩模提取操作
for ras in rasters:
outname = os.path.join(outws, os.path.basename(ras).split(".")[0] + "_clp.tif") # 指定输出文件的命名方式(以被裁剪文件名+_clip.tif命名)
out_extract = arcpy.sa.ExtractByMask(ras, mask) # 执行按掩模提取操作
out_extract.save(outname) # 保存数据在这里插入代码片
Arcpy 批量掩膜裁剪栅格
最新推荐文章于 2024-07-01 22:04:01 发布