数据：青藏高原250米逐日水体制图（2000-2022）

网址：https://data.tpdc.ac.cn/zh-hans/data/763cfc76-4f18-4c11-b45b-ec23d5fc769d

数据格式：

Tif文件内容：每个tif文件里有365或者366个band，每天是一个band

第一步（波段分离）：

这里我只选择了两个图幅进行波段分离，因为只需要这两个图幅范围内的信息

输出文件名为：

下边是python程序

第二步（空间图幅合并）：

这里是用的gdal进行图幅合并，基于linux系统写的bash脚本，我的系统gdal使用时需要加载对应的module

因为只需要两个图幅，这里是把这两个图幅的每天的数据对应进行合并，合并后文件为：

Linux的sh脚本为：

第三步（提取tif文件到文本）：

这里我是用一个流域的不规则多边形（shp格式）对合并后的图幅进行裁剪，裁剪出流域范围，然后统计流域范围内各个值的网格数（值是通过arcgis预览获取的）进行后续的分析

输出结果为：

这里total pixels和zero pixels有问题，它是不规则多边形所在的矩形框内所有的网格（如下图黑色区域）。但是有值的没问题，结果是不规则多边形内的网格数。

其它：（获取不规则多边形内部的网格数）：

首先用arcgis的extract by mask工具，用流域的polygon.shp裁剪merge.tif，然后用下边的python程序获取不规则多边形内部的各网格数