利用 arcpy 将栅格按照要素裁剪

可以通过 ArcToolbox 中的 Extract by Mask 工具来实现。这个工具可从栅格数据集中提取指定掩模下的像元值。

以下是使用 arcpy 进行栅格裁剪操作的步骤：

1. 打开 ArcMap，导入需要处理的栅格和矢量数据，并确保它们都以正确的投影系统提供。

2. 打开 ArcToolbox 界面，找到 Spatial Analyst Tools 工具集中的 Extraction 工具箱。然后选择 Extract by Mask 函数。

3. 在 Extract by Mask 函数中输入如下参数：

  - Input raster：需要裁剪的栅格数据路径，例如："C:/data/raster.tif"。
  - Output workspace：输出栅格数据集存储目录路径，例如："C:/output"。
  - Extract mask：作为裁剪边界的矢量图层（shapefile或feature class）路径，例如： "C:/data/polygon.shp"。
  - Invert mask：是否应该反转提取掩模，即提取裁剪范围之外的区域。默认为“False”。
  
  import arcpy

  # Set input and output file paths
  in_raster = "D:/path/to/input/raster.tif"
  out_workspace = "D:/path/to/output"
  mask_shapefile = "D:/path/to/mask/boundary.shp"

  # Extract the Raster Dataset using the mask boundary shapefile
  arcpy.gp.ExtractByMask_sa(in_raster, mask_shapefile, out_workspace)
  ```

4. 确定所有参数后，单击“运行”按钮执行 Extract by Mask 函数。

5. 提取成功后，在输出目录下找到新的栅格数据集并在 ArcMap 中打开它，即可看到已经按掩模切割后的栅格结果。

注意：要素掩蔽裁剪操作仅适用于参照空间范围内的输入栅格数据。如果掩蔽范围与输入栅格数据范围不匹配，则需要先对输入栅格进行裁剪或者进行重采样以保持一致。此外，执行该函数时需要注意性能和计算资源限制，因为该函数需要对整个栅格数据集进行处理且可能涉及大量像元处理。