[发明专利]一种水体污染的检测方法、装置、设备以及存储介质

权利要求书

1.一种水体污染的检测方法，其特征在于，包括：

对原始遥感影像数据进行大气校正，得到第一遥感影像数据，其中，所述原始遥感影像数据为与矿区相邻水体所在地表的卫星图像数据；

对所述第一遥感影像数据进行水体校正，得到所述第一遥感影像数据中水体的地表反射率，其中，所述水体包含叶绿素a和悬浮物；

基于所述地表反射率，建立所述叶绿素a与所述水体的光谱波段的第一经验反演模型、以及所述悬浮物与所述水体的光谱波段的第二经验反演模型；

将包含所述水体且经过大气校正后的第二遥感影像数据输入所述第一经验反演模型和所述第二经验反演模型，得到所述叶绿素a的浓度和所述悬浮物的浓度；

根据所述叶绿素a的浓度和所述悬浮物的浓度、识别所述水体与所述矿区的污染关系。

2.根据权利要求1所述的方法，其特征在于，所述对原始遥感影像数据进行大气校正，得到第一遥感影像数据，包括：

获取所述原始遥感影像数据的大气参数值；

将所述大气参数值输入到Flaash模块中，并选择大气模型的类型以对原始遥感影像数据进行大气校正，得到所述第一遥感影像数据。

3.根据权利要求1所述的方法，其特征在于，所述对所述第一遥感影像数据进行水体校正，得到所述第一遥感影像数据中水体的地表反射率，包括：

采用水表面以上的测量方法对与矿区相邻的水体进行采样，得到第一光谱数据；

对所述第一光谱数据进行低通滤波，得到所述水体的第二光谱数据；

将所述第二光谱数据中水面以上的上行辐射率减去天空漫散射光的辐射率，得到离水辐亮度；

将所述离水辐亮度除以水面入射辐照度，得到所述水体的遥感反射率；

根据所述水体的遥感反射率，计算所述第一遥感影像数据中水体的地表反射率。

4.根据权利要求3所述的方法，其特征在于，所述采用水表面以上的测量方法对与矿区相邻的水体进行采样，得到第一光谱数据，包括：

采集仪器的观测平面与太阳入射平面的夹角在90度到135度范围内；

所述采集仪器与水体平面的法线方向的夹角在30度到45度范围内；

选择在湖面风平浪静的情况下，从10：00到16：00间对所述水体进行光谱测量，得到所述第一光谱数据。

5.根据权利要求3所述的方法，其特征在于，所述根据所述水体的遥感反射率，计算所述第一遥感影像数据中水体的地表反射率，包括：

根据以下公式计算所述第一遥感影像数据中水体的地表反射率：

其中R是所述第一遥感影像数据中水体的地表反射率，λ为卫星传感器的光谱通道的波长，λ_min为卫星传感器的光谱通道的起始波长，λ_max为卫星传感器的光谱通道的终止波长，r(λ)为对应λ波长的所述水体的遥感反射率，f(λ)为光谱响应函数。

6.根据权利要求5所述的方法，其特征在于，所述基于所述地表反射率，建立所述叶绿素a与所述水体的光谱波段的第一经验反演模型、以及所述悬浮物与所述水体的光谱波段的第二经验反演模型，包括：

采用单波段阈值法对所述地表反射率进行筛选，提取所述第一遥感影像数据中水体的第三光谱数据；

计算所述第三光谱数据中各波段组合值与所述叶绿素a之间的第一相关系数；

计算所述第三光谱数据中各波段组合值与所述悬浮物之间的第二相关系数；

选择所述第一相关系数为最大值的波段组合值与所述叶绿素a建立所述第一经验反演模型为Y1＝-2.1139*(TM3-TM2)+0.0043；

选择所述第二相关系数为最大值的波段组合值与所述悬浮物建立所述第二经验反演模型为Y2＝2740.9*(TM3-TM1)+10.825；

其中，Y1为叶绿素a的浓度，Y2为悬浮物的浓度，TM1为遥感影像的蓝光波段，TM2为遥感影像的绿光波段，TM3为遥感影像的红光波段。