通过arcpy模块,利用.shp裁剪影像栅格的python代码
时间: 2023-11-28 07:49:20 浏览: 196
ArcGIS根据SHP裁剪栅格的Python代码
可以使用arcpy中的ExtractByMask函数来实现影像栅格与.shp面数据的裁剪操作。下面是一个简单的示例代码:
```python
import arcpy
# 设置工作空间
arcpy.env.workspace = r"C:\path\to\workspace"
# 定义输入的影像栅格和.shp面数据
raster_path = "input_raster.tif"
shapefile_path = "input_shapefile.shp"
# 读取.shp面数据
shapefile = arcpy.FeatureClassToFeatureClass_conversion(shapefile_path, arcpy.env.workspace, "temp.shp")
# 使用ExtractByMask进行裁剪操作
out_raster = arcpy.sa.ExtractByMask(raster_path, shapefile)
# 保存输出影像
out_raster.save("output_raster.tif")
```
需要注意的是,上述代码中的路径需要根据实际情况进行修改。另外,使用ExtractByMask函数时,输入的影像栅格和.shp面数据需要具有相同的坐标系。
阅读全文