裁剪 (Data Management)
摘要
创建包括栅格数据集、镶嵌数据集或影像服务图层在内的栅格空间子集。
插图
 |
用法
-
使用此工具可以基于模板范围提取部分栅格数据集。该裁剪输出包含与模板范围相交的所有像素。如果要提取要素数据集的一部分，可以使用分析工具箱中的裁剪工具。
-
通过使用以 x 和 y 坐标的最小值和最大值确定的包络矩形或使用输出范围文件来指定剪切区域。如果指定的裁剪范围与输入栅格数据集没有对齐，“裁剪”工具可确保使用适当的对齐。这可能使输出的实际范围与此工具中指定的范围略有不同。
-
可以将已有的栅格或矢量图层作为剪切范围使用。如果使用要素类作为输出范围，则可以通过要素类的最小外接矩形或要素的面几何来裁剪栅格。使用裁剪几何时，可能提高输出的像素深度。因此，需要确保输出格式可以支持适当的像素深度。
-
在使用 ArcMap 时，也可以将选中的要素作为剪切范围使用。如果选择要素类中的某一要素并选中选择范围（将 clipping_geometry 设置为 ClippingGeometry），则输出将裁减掉选定区域。如果选择要素类中的某一要素，但未选中选择范围，则输出会将该要素的最小外接矩形裁剪掉。
-
可将输出结果保存为 BIL、BIP、BMP、BSQ、DAT、GIF、Esri Grid、IMG、JPEG、JPEG 2000、PNG、TIFF 格式或任意地理数据库栅格数据集。
-
该范围值必须与栅格数据集在同一空间坐标系内，并使用相同的单位。
语法
| 参数 | 说明 | 数据类型 |
in_raster |
输入栅格数据集。 | Mosaic Dataset; Mosaic Layer; Raster Dataset; Raster Layer |
rectangle |
用于定义裁剪所需的最小外接矩形的四个坐标按以下顺序定义：X 最小值、Y 最小值、X 最大值和 Y 最大值。 如果指定的裁剪范围没有与输入栅格数据集对齐，则“裁剪”工具可确保使用适当的对齐。这可能使输出的实际范围与此工具中指定的范围略有不同。 | Envelope |
out_raster |
输出栅格数据集。确保输出格式可以支持适当的像素深度。 以文件格式存储栅格数据集时，需要指定文件扩展名，具体如下：
以地理数据库形式存储栅格数据集时，不应向栅格数据集的名称添加文件扩展名。 将栅格数据集存储到 JPEG 文件、JPEG 2000 文件、TIFF 文件或地理数据库时，可以指定压缩类型和压缩质量。 | Raster Dataset |
in_template_dataset (可选) |
可以将现有的栅格或矢量图层作为裁剪范围使用。该裁剪输出包括了与最小外接矩形相交的所有像素。 将某要素类作为输出范围并且要基于面要素进行栅格裁剪时，请使用 clipping_geometry 选项。选中此项可以提高输出的像素深度。因此，需要确保输出格式可以支持适当的像素深度。 | Raster Layer; Feature Layer |
nodata_value (可选) |
具有指定值的所有像素将在输出栅格数据集中被设置为 NoData。 | String |
clipping_geometry (可选) |
如果使用要素类作为输出范围，则可以通过该要素类的范围或面周长来裁剪栅格。
使用裁剪几何时，可能提高输出的像素深度。因此，需要确保输出格式可以支持适当的像素深度。 | Boolean |
代码实例
这是“裁剪”工具的 Python 示例。
import arcpy
arcpy.Clip_management("c:/data/image.tif","1952602 294196 1953546 296176",
"c:/data/clip.gdb/clip01", "#", "#", "NONE")
这是“裁剪”工具的 Python 脚本的示例。
##====================================
##Clip
##Usage: Clip_management in_raster rectangle out_raster {in_template_dataset} {nodata_value} {NONE | ClippingGeometry}
try:
import arcpy
arcpy.env.workspace = r"C:/Workspace"
##Clip Raster Dataset by known extent - Left Bottom Right Top
arcpy.Clip_management("image.tif","1952602.23 294196.279 1953546.23 296176.279","clip.gdb/clip", "#", "#", "NONE")
##Clip Raster Dataset with feature geometry
arcpy.Clip_management("image.tif", "#", "clip.tif","feature.shp", "0", "ClippingGeometry")
except:
print "Clip example failed."
print arcpy.GetMessages()