gma 教程 | 栅格处理 | 矢量裁剪栅格

目标

【裁剪】。按矢量边界裁剪栅格。

环境

系统： Window 10+ (X64)
Python 版本： 3.8.8+
gma 版本： 1.0.7+

gma 安装和详细功能帮助见：地理与气象分析库

函数

gma.rasp.Clip(InFile, OutFile, CutLineFile, LayerID = 0, FeatureID = None, InNoData = None, OutNoData = None, OutFormat = ‘GTiff’)

参数：

 InFile：str。输入栅格路径。

 OutFile：str。输出文件路径。

 CutLineFile：str。裁剪矢量文件路径。

可选参数：

 LayerID = int 。输入裁剪矢量的图层 ID。默认第一个图层（0）。

 FeatureID = list 。掩膜图层的要素 ID。默认掩膜全部要素（None）。

 InNoData = float。输入栅格的无效值。默认不指定（None）无效值。

 OutNoData = float。输出栅格的无效值。默认不指定（None）无效值。

 OutFormat = str。输出栅格文件格式。默认为 GTiff，其他格式详见 ToOtherFormat 函数。

案例

数据

栅格数据： ESA 2020年陆表覆盖数据
矢量数据： 河南省18个地级行政单位矢量边界

裁剪

import gma# 定义原始栅格和裁剪矢量
InFile = r'D:\CSDN\gma\栅格处理\0.1原始影像\ESA_Henan_2020.tif'
CutLineFile = r'D:\CSDN\gma\栅格处理\0.2原始矢量\Henan_City.shp'# 定义裁剪结果路径并执行裁剪
OutFile = r'D:\CSDN\gma\栅格处理\Clip\ESA_Henan_2020_Clip.tif'
gma.rasp.Clip(InFile, OutFile, CutLineFile, OutNoData=0)

裁剪结果如下：

如果需要按矢量文件中的某个要素裁剪，可增加 FeatureID 参数，示例如下：

如果有需求，可先查看矢量属性表，获取要素信息，方法如下：

# 打开矢量数据集
VEDataSet = gma.Open(CutLineFile)
# 获取矢量图层
LY = VEDataSet.GetLayer(0)
# 获取属性表
ATable = LY.GetAttributeTable()
print(ATable)

我们裁剪河南省洛阳市，其 FeatureID 为 2，配置参数如下：

### 按 图层ID 裁剪
OutFileFID = r'D:\CSDN\gma\栅格处理\Clip\ESA_Henan_2020_Clip_FID2.tif'
gma.rasp.Clip(InFile, OutFileFID, CutLineFile, FeatureID = [2], OutNoData=0)

执行结果如下所示：

疑问与反馈

邮箱：858998723@qq.com