裁剪 (数据管理)
摘要
创建栅格数据集的空间子集。
插图
 |
用法
-
使用此工具可以基于模板范围提取部分栅格数据集。该裁剪输出包含与模板范围相交的所有像素。
-
通过使用以 x 和 y 坐标的最小值和最大值确定的包络矩形或使用输出范围文件来指定剪切区域。如果指定的裁剪范围与输入栅格数据集没有对齐，“裁剪”工具可确保使用适当的对齐。这可能使输出的实际范围与此工具中指定的范围略有不同。
-
可以将已有的栅格或矢量图层作为剪切范围使用。如果使用要素类作为输出范围，则可以通过要素类的最小外接矩形或要素的面几何来裁剪栅格。使用裁剪几何时，可能提高输出的像素深度。因此，需要确保输出格式可以支持适当的像素深度。
-
在使用 ArcMap 时，也可以将选中的要素作为剪切范围使用。如果选定要素类中的某一要素并选中将输入要素用于裁剪几何，则输出将裁减掉选定区域。如果选中要素类中的某一要素，但未选中将输入要素用于裁剪几何，则输出将裁剪掉该要素的最小外接矩形。
-
可将输出保存为 BIL、BIP、BMP、BSQ、DAT、GIF、GRID、IMG、JPEG、JPEG 2000、PNG、TIFF 格式或任意地理数据库栅格数据集。
-
如果指定的裁剪范围与输入栅格数据集没有对齐，“裁剪”工具可确保使用适当的对齐。这可能使输出的实际范围与此工具中指定的范围略有不同。
-
如果使用要素类作为输出范围，则可以通过该要素类的范围或面周长来裁剪栅格。
-
该范围值必须与栅格数据集在同一空间坐标系内，并使用相同的单位。
语法
| 参数 | 说明 | 数据类型 |
in_raster |
输入栅格数据集。 | Raster layer |
rectangle |
裁剪所需的定义最小外接矩形的四个坐标。定义顺序如下：X 最小值、Y 最小值、X 最大值和 Y 最大值。 如果指定的裁剪范围与输入栅格数据集没有对齐，则“裁剪”工具可确保使用适当的对齐。这可能使输出的实际范围与此工具中指定的范围略有不同。 | Envelope |
out_raster |
输出栅格数据集。确保输出格式可以支持适当的像素深度。 以文件格式存储栅格数据集时，需要指定文件扩展名，具体如下：
以地理数据库形式存储栅格数据集时，不应向栅格数据集的名称添加文件扩展名。 将栅格数据集存储到 JPEG 文件、JPEG 2000 文件、TIFF 文件或地理数据库时，可以指定压缩类型和压缩质量。 | Raster dataset |
in_template_dataset （可选） |
可作为剪切范围使用的已有栅格或矢量图层。该裁剪输出包括与最小外接矩形相交的所有像素。 将某要素类作为输出范围且想要基于面要素进行栅格裁剪时，请选择“裁剪几何”选项。使用“裁剪几何”时，可能提高输出的像素深度。因此，需要确保输出格式可以支持适当的像素深度。 | Raster Layer; Feature Layer |
nodata_value （可选） |
具有指定值的所有像素将在输出栅格数据集中被设置为 NoData。 | Double |
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()