arcpy批量裁剪遥感影像数据

可以使用arcpy中的ExtractByMask函数来批量裁剪遥感影像数据。下面是一个代码,用于将一个目录(input_folder)内的所有遥感影像数据根据指定的掩码文件(mask_file)进行裁剪处理:


import arcpy, os

# 设置工作空间和输入数据路径
arcpy.env.workspace = r"C:\data\Imagery.gdb"
input_folder = r"C:\input"
mask_file = r"C:\data\Masks.gdb\LandCover"

# 循环遍历目录内所有影像,并且执行ExtractByMask操作
for file in os.listdir(input_folder):
    if file.endswith(".img") or file.endswith(".tif"):
        input_raster = os.path.join(input_folder, file)
        output_raster = os.path.join(arcpy.env.workspace, "Clip_" + file)
        arcpy.sa.ExtractByMask(input_raster, mask_file).save(output_raster)

代码中,首先设置工作空间和输入数据路径,然后循环遍历指定目录下的所有影像文件。如果该文件是以".img"或者".tif"结尾的遥感影像,则将其加载到ExtractByMask函数中,同时指定掩码文件。函数执行完成后,使用Save函数将结果保存到输出路径中。

需要注意的是,代码仅适用于掩码文件(mask_file)为单一目标范围的情况下。如果需要对多个目标范围进行裁剪操作,则需要修改代码,例如使用循环或者使用ExtractMultipleLayersByMask函数,来实现批量裁剪所有目标范围。

  • 1
    点赞
  • 2
    收藏
    觉得还不错? 一键收藏
  • 打赏
    打赏
  • 0
    评论

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包

打赏作者

认真学GIS

你的鼓励将是我创作的最大动力

¥1 ¥2 ¥4 ¥6 ¥10 ¥20
扫码支付:¥1
获取中
扫码支付

您的余额不足,请更换扫码支付或充值

打赏作者

实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

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

余额充值