基于arcpy, 利用sdjzu.shp对GF1_PMS2_fuse2m_clip.tif进行裁切,使用得裁切后影像的形状与sdjzu.shp完全一致。
时间: 2023-07-15 15:14:11 浏览: 165
arcmap中用clip裁切数据
3星 · 编辑精心推荐
可以使用以下代码实现:
```python
import arcpy
# 设置工作空间
arcpy.env.workspace = r"path/to/workspace"
# 定义输入数据路径
input_raster = r"path/to/GF1_PMS2_fuse2m_clip.tif"
input_shapefile = r"path/to/sdjzu.shp"
# 定义输出数据路径
output_raster = r"path/to/output.tif"
# 获取裁切范围
extent = arcpy.Describe(input_shapefile).extent
# 裁切栅格数据
arcpy.Clip_management(input_raster, str(extent), output_raster, input_shapefile, "#", "ClippingGeometry")
print("裁切完成!")
```
其中,`arcpy.Describe(input_shapefile).extent`可以获取到shapefile的范围,`arcpy.Clip_management()`可以进行裁剪操作,`ClippingGeometry`参数指定了裁剪的几何形状。
阅读全文