Python利用GDAL模块实现读取栅格数据并对指定数据加以筛选掩膜

目录

1 代码分段讲解

1.1 模块与路径准备1.2 栅格图像文件名读取与配对1.3 输出文件名称准备1.4 栅格文件数据与信息读取1.5 差值计算与QA波段筛选1.6 结果栅格文件写入与保存

2 完整代码

本文介绍基于Python语言中gdal模块，对遥感影像数据进行栅格读取与计算，同时基于QA波段对像元加以筛选、掩膜的操作。

本文所要实现的需求具体为：现有自行计算的全球叶面积指数（LAI）.tif格式栅格产品（下称“自有产品”），为了验证其精确度，需要与已有学者提出的成熟产品——GLASS全球LAI.hdf格式栅格产品（下称“GLASS产品”）进行做差对比；其中，自有产品除了LAI波段外，还有一个<strong>质量评估波段