ArcGIS 10.7中裁剪栅格的核心方法有两种：使用图形用户界面（GUI）工具和利用ArcPy或模型构建器进行批处理。下面将分别详细介绍这两种方法的操作步骤。 ### 一、使用图形用户界面工具进行单次裁剪 ArcGIS 10.7提供了多个用于裁剪栅格的工具，其中“按掩膜提取”（Extract by Mask）是最常用、最灵活的工具。以下是其详细操作步骤： 1. **打开ArcMap 10.7**，将待裁剪的栅格数据（例如.tif/.img文件）和作为裁剪范围的掩膜数据（通常是.shp格式的多边形矢量数据，如行政区划）加载到地图中[ref_2][ref_5]。 2. **搜索并打开工具**。 * 在菜单栏中点击**Geoprocessing（地理处理）** → **ArcToolbox**。 * 或者直接在右侧的搜索框中输入“按掩膜提取”，点击搜索。 3. **定位并打开裁剪工具**。 * 在ArcToolbox中导航至：`Spatial Analyst Tools（空间分析工具）` → `Extraction（提取）` → `Extract by Mask（按掩膜提取）`[ref_1][ref_5]。 * 双击打开“按掩膜提取”对话框。 4. **设置工具参数**[ref_1][ref_2]。 * **Input raster（输入栅格）**：选择你想要裁剪的栅格图层。 * **Input raster or feature mask data（输入栅格或要素掩膜数据）**：选择用于定义裁剪范围的.shp矢量图层。 * **Output raster（输出栅格）**：指定裁剪后文件的保存路径和名称。建议使用英文命名，避免路径中出现中文。 * （可选）可以根据需要设置输出像元大小等高级选项。 5. **运行并检查结果**。 * 点击“OK”按钮运行工具。处理完成后，裁剪结果会自动添加到当前地图中。 | 方法 | 工具名称 | 工具箱位置 | 适用场景 | | :--- | :--- | :--- | :--- | | **推荐** | **按掩膜提取 (Extract by Mask)** | Spatial Analyst Tools -> Extraction | 最常用，使用矢量或多边形作为精确边界进行裁剪[ref_1][ref_2]。 | | 其他方法 | 裁剪 (Clip) | Data Management Tools -> Raster -> Raster Processing | 常用于使用矩形范围进行简单裁剪，但“按掩膜提取”功能更通用[ref_5]。 | | 其他方法 | 分割栅格 (Split Raster) | Data Management Tools -> Raster -> Raster Processing | 用于将一个大栅格按指定行数列数或按多边形要素拆分成多个小文件[ref_2]。 | ### 二、使用模型构建器进行批量裁剪 当需要对多个栅格文件进行相同区域的裁剪时，手动操作效率低下。可以使用**模型构建器**实现自动化批量处理[ref_3]。 1. **打开模型构建器**。 * 在菜单栏中点击 **Geoprocessing（地理处理）** → **ModelBuilder（模型构建器）**，打开一个新的模型窗口[ref_3]。 2. **添加迭代器和裁剪工具**。 * 在模型构建器工具箱中，找到 **Iterators（迭代器）**，将其中的“**Iterate Rasters（迭代栅格）**”拖拽到模型画布中。 * 配置“迭代栅格”参数：将`Workspace or Raster Catalog`设置为存放所有待裁剪栅格文件的文件夹路径[ref_3]。 3. **连接并配置裁剪工具**。 * 从ArcToolbox中，将`Spatial Analyst Tools -> Extraction -> Extract by Mask`工具拖到模型画布中。 * 使用连接线（Connect Tool）将“迭代栅格”输出的 `Raster` 连接到“按掩膜提取”工具的 `Input raster` 上。 * 双击“按掩膜提取”工具，在其对话框中设置好固定的`Input raster or feature mask data`（掩膜.shp文件）[ref_3]。 4. **设置输出路径**。 * 这是关键步骤，需用变量构造动态输出路径。右键点击模型的空白处，选择 `Create Variable` -> `From Parameter` -> 选择 `Output raster`。 * 双击这个新创建的`Output raster`变量，在值（Value）中使用 **%Name%** 等内联变量模型表达式来构造输出名，例如：`E:\output_folder\%Name%_clipped.tif`，这能确保每个输出文件以输入栅格名称加后缀的形式命名[ref_3]。 5. **保存并运行模型**。 * 保存模型。点击模型构建器工具栏上的“运行”（Run）按钮。模型会自动遍历输入文件夹中的所有栅格，并使用同一掩膜分别裁剪，输出到指定文件夹。 ### 三、使用ArcPy脚本进行高级批处理 对于更复杂的批量任务（如处理嵌套子文件夹、定制命名规则），编写Python脚本是最强大的方法。核心函数是 **`arcpy.sa.ExtractByMask`**。 以下是一个基础示例脚本，演示如何批量处理一个文件夹中的所有栅格： ```python # -*- coding: utf-8 -*- import arcpy from arcpy.sa import ExtractByMask # 导入按掩膜提取函数 import os # 设置工作环境（可选） arcpy.env.workspace = r"C:\YourInputFolder" # 设置包含待裁剪栅格的文件夹路径 arcpy.env.overwriteOutput = True # 允许覆盖同名输出文件 # 定义输入参数 mask_shp = r"C:\YourMaskFolder\boundary.shp" # 裁剪范围的矢量文件路径[ref_1] output_folder = r"C:\YourOutputFolder" # 输出结果文件夹路径 # 检查Spatial Analyst扩展模块是否可用 if arcpy.CheckExtension("Spatial") == "Available": arcpy.CheckOutExtension("Spatial") else: raise Exception("Spatial Analyst license is not available.") # 遍历工作空间中的所有栅格文件 # arcpy.ListRasters() 可以指定文件类型，如 "*.tif"，这里列出所有栅格 raster_list = arcpy.ListRasters() for raster in raster_list: print(f"正在处理: {raster}") # 构建输出栅格路径和名称 # 这里使用原文件名加 "_Clipped" 后缀 out_raster_name = os.path.splitext(raster)[0] + "_Clipped.tif" out_raster_path = os.path.join(output_folder, out_raster_name) # 执行按掩膜提取操作[ref_1] try: outExtractByMask = ExtractByMask(raster, mask_shp) outExtractByMask.save(out_raster_path) print(f"成功保存至: {out_raster_path}") except Exception as e: print(f"处理 {raster} 时出错: {e}") # 归还扩展模块许可 arcpy.CheckInExtension("Spatial") print("批量裁剪处理完成。") ``` * **应用场景**：此脚本适用于处理单个文件夹内的大量栅格，例如一个年份内的所有月度NDVI数据。 * **扩展应用**：如参考资料所述，在处理多级目录（如`年/月/日.tif`）时，可以使用`os.walk()`函数递归遍历所有子文件夹，实现更复杂的组织结构处理[ref_1]。