基于arcpy, 利用sdjzu.shp对GF1_PMS2_fuse2m_clip.tif进行裁切,使用得裁切后影像的形状与sdjzu.shp完全一致。
时间: 2023-07-16 13:15:03 浏览: 111
GRB2网格预报数据发布影像服务.zip_.grb数据_Arcpy_GRB2格式数据_arcgis 气象_气象数据
您可以使用以下代码进行裁剪:
```python
import arcpy
# 设置工作空间和输出路径
arcpy.env.workspace = r"path/to/workspace"
output_path = r"path/to/output"
# 定义输入和裁切要素
input_raster = "GF1_PMS2_fuse2m_clip.tif"
clip_feature = "sdjzu.shp"
# 获取裁剪要素的范围
desc = arcpy.Describe(clip_feature)
ext = desc.extent
# 定义输出影像的名称和路径
output_raster = output_path + "\output.tif"
# 调用裁剪函数
arcpy.Clip_management(input_raster, str(ext), output_raster, clip_feature, "#", "ClippingGeometry")
print("裁剪完成!")
```
在代码中,首先设置了工作空间和输出路径,然后定义了输入的影像和裁剪要素。接着,使用 `arcpy.Describe()` 函数获取了裁剪要素的范围,并将其转换为字符串格式传递给 `arcpy.Clip_management()` 函数。最后,将输出影像的名称和路径定义为变量 `output_raster`,并调用裁剪函数。
运行完毕后,输出影像的形状应该与裁剪要素完全一致。
阅读全文