20241212
要求：
（1）输入原始影像是反射率数据，每步都是可调用的单独的程序
（2）无效像元对应的值必须为0，data ignore value = 0
（3）输入影像必须为bsq
（4）实测数据的格式为3列：经度，纬度，含量


1 提取准确水体范围（1）自动实现ndwi阈值（会出现提取错误，例如孤岛现象）（2）矢量shp文件转栅格
输出的水域掩膜和影像有相同行列数和分辨率

2 otsu找过曝/耀斑严重区域（是否通过dn值自动识别）

3 去耀斑的方法：添加水域掩膜

4 修正模型参数：（1）通过实测数据（避过过曝/耀斑严重的位置）修正模型参数，【任总决定不反演】，现在还是通过envi做的（2）模型是否合适（3）csv规定格式  json
        （4）输出修正模型信息（怎么输出？）（5）随机选择几个点做精度评价，需要将实测数据分为反演集合和测试集合
输出的模型信息为json格式：模型类型、模型参数、模型精度（带坐标）

5 根据输入（原始影像）间隔（先像素后距离）和取点策略（仅实现一种） 输出一系列是水体且【避过过曝/耀斑严重的位置】的像素 的csv坐标（修改）

6 根据上步坐标，获取光谱和其真实坐标位置
7 对第4、6步数据进行计算，输出第5步关键点的反演值
8 插值：和arcgis的结果有差异