Arcgis下python实现矢量shp批量裁剪栅格影像(亲测可用)

以下代码,需要在arcgis系统的python编程窗口实现,需要替换栅格影像文件夹路径,矢量shp文件路径,及裁剪后的输出路径。

另外,需要注意,栅格格式如果不是tif,也需要修改为对应的栅格类型。

# -*- coding: utf-8 -*-
import arcpy  
arcpy.CheckOutExtension("spatial")  
arcpy.gp.overwriteOutput=0
arcpy.env.workspace = r'D:\1paper\1ETandP\dayquchunull'  #被裁剪栅格影像所在文件夹
rasters = arcpy.ListRasters("*", 'tif')  #栅格数据格式设置
mask= r'D:\1paper\1ETandP\liuyu\touy.shp'   #裁剪范围shp文件路径
outpath=r'D:\1paper\1ETandP\dayquchunullcj'  #影像裁剪结果输出文件夹

for raster in rasters:  
    print(raster)  
    out=outpath+'\\'+ raster.split(".")[1]+ ".tif"  #重命名,将原文件以.为分割,提取哪一段为新名,0为.前的第1段字符
    arcpy.gp.ExtractByMask_SA(raster, mask, out)  
    print(raster+"  has done")  
print("All done") 
 

  • 9
    点赞
  • 2
    收藏
    觉得还不错? 一键收藏
  • 1
    评论
评论 1
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值