【ArcGIS自定义脚本工具】批量裁剪栅格（多对多）

Salierib

已于 2022-06-19 16:55:27 修改

阅读量1.2k

点赞数 2

分类专栏： ArcGIS自定义脚本编程文章标签： python ArcGIS

于 2021-02-13 11:23:09 首次发布

本文链接：https://blog.csdn.net/qq_37948866/article/details/113727873

版权

ArcGIS自定义脚本编程专栏收录该内容

24 篇文章 88 订阅

订阅专栏

该博客介绍了如何利用ArcGIS的Python脚本进行批量裁剪栅格影像操作。通过输入待裁剪的栅格文件路径、裁剪掩模和输出目录，脚本将自动生成对应裁剪结果，并保存到指定文件夹中。同时，脚本会检查输出文件夹是否存在，避免重复裁剪并提供进度反馈。

摘要由CSDN通过智能技术生成

系列文章目录： ArcGIS自定义脚本编程

一、功能介绍

在这里插入图片描述

二、脚本代码

#!/usr/bin/python
# -*- coding: UTF-8 -*-


import arcpy
from arcpy import env
import os
import time
import sys
reload(sys)
sys.setdefaultencoding('utf8')

tifs = arcpy.GetParameterAsText(0)
masks = arcpy.GetParameterAsText(1)
out_dir = arcpy.GetParameterAsText(2)
tifs = tifs.split(";")
masks = masks.split(";")
names = [os.path.splitext(os.path.basename(mask))[0] for mask in masks]

size = len(tifs) * len(masks)
num = 1
for i, mask in enumerate(masks):
    # create a new folder named by mask's name
    new_folder = out_dir + os.sep + names[i]
    if not os.path.exists(new_folder):
        os.mkdir(new_folder)
    else:
        arcpy.AddMessage("Folder {0} already exists. Please check it.".format(new_folder))
    for tif in tifs:
        s = time.time()
        cliped_tif = os.path.join(new_folder, "{0}_{1}".format(names[i],os.path.split(tif)[1]))
        if not os.path.exists(cliped_tif):
            arcpy.Clip_management(tif, "#", cliped_tif, mask, "#", "ClippingGeometry")
            e = time.time()
            arcpy.AddMessage("{0}/{1} | {2} Completed, time used {3}s".format(num, size, cliped_tif, e - s))
        else:
            e = time.time()
            arcpy.AddMessage("{0}/{1} | {2} already exists.".format(num, size, cliped_tif))
        num += 1