ArcPy 使用Shp批量裁剪栅格 Clip_maintain clipping extent

 小记一下 Clip_management 的使用方法

inRaster = "d.tif" #输入栅格路径
outRaster = r"C:/Users/xhz36/Documents/ArcGIS/Default.gdb/c_Clip" #输出栅格路径
clipShp = "zlshp" #裁剪shp
nodata = 256 #Nodata值
 
def clipRasterbyShp(inRaster,outRaster,clipShp):
    ext = arcpy.Describe(clipShp).Extent
    arcpy.Clip_management(inRaster,"%.12f %.12f %.12f %.12f"%(ext.XMin,ext.YMin,ext.XMax,ext.YMax),outRaster,clipShp,"NoData","ClippingGeometry","NO_MAINTAIN_EXTENT")
 
clipRasterbyShp(inRaster,outRaster,clipShp)

Clip_management (in_raster, rectangle, out_raster, {in_template_dataset}, {nodata_value}, {clipping_geometry}, {maintain_clipping_extent})

输入的extent字符串，精度一般在9~12位小数，若精度较低，在贴边时可能会产生，Clip Shp超出Raster 范围的错误。

参数说明：
in_raster 输入栅格数据集：Mosaic Dataset; Mosaic Layer; Raster Dataset; Raster Layer

rectangle 用于定义裁剪所需的最小外接矩形的四个坐标按以下顺序定义：X 最小值、Y 最小值、X 最大值和 Y 最大值。

如果指定的裁剪范围没有与输入栅格数据集对齐，则“裁剪”工具可确保使用适当的对齐。这可能使输出的实际范围与此工具中指定的范围略有不同。

 Envelope 
out_raster 输出栅格数据集。确保输出格式可以支持适当的像素深度。

以文件格式存储栅格数据集时，需要指定文件扩展名，具体如下：

.bil - Esri BIL 
.bip - Esri BIP 
.bmp - BMP 
.bsq - Esri BSQ 
.dat - ENVI DAT 
.gif - GIF 
.img - ERDAS IMAGINE 
.jpg - JPEG 
.jp2 - JPEG 2000 
.png - PNG 
.tif - TIFF 
无扩展名 - Esri Grid
以地理数据库形式存储栅格数据集时，不应向栅格数据集的名称添加文件扩展名。

将栅格数据集存储到 JPEG 文件、JPEG 2000 文件、TIFF 文件或地理数据库时，可以指定压缩类型和压缩质量。

in_template_dataset 输入的矢量或栅格模板
(可选) 可以将现有的栅格或矢量图层作为裁剪范围使用。该裁剪输出包括了与最小外接矩形相交的所有像素。

将某要素类作为输出范围并且要基于面要素进行栅格裁剪时，请使用 clipping_geometry 选项。选中此项可以提高输出的像素深度。因此，需要确保输出格式可以支持适当的像素深度。

nodata_value  (可选) 具有指定值的所有像素将在输出栅格数据集中被设置为 NoData。

clipping_geometry  (可选) 如果使用要素类作为输出范围，则可以通过该要素类的范围或面周长来裁剪栅格。

NONE —基于要素类的最小外接矩形对栅格数据集进行裁剪。 
ClippingGeometry —基于面形状的周长对栅格数据集进行裁剪。 
使用裁剪几何时，可能提高输出的像素深度。因此，需要确保输出格式可以支持适当的像素深度。

maintain_clipping_extent  (可选) 选择剪裁后的输出栅格数据集的范围。

NO_MAINTAIN_EXTENT —输出栅格将与输入栅格保持相同的像元对齐方式，输出范围将相应进行调整。这是默认设置。
MAINTAIN_EXTENT —输出范围将与指定的剪裁范围完全匹配，行数和列数将相应进行调整。由于行数和列数会根据输出进行调整，因此，将对输入值重采样。
 Boolean