[发明专利]一种基于OLCI的内陆富营养化湖泊水体蓝藻生物量比例遥感估算方法

权利要求书

1.一种基于OLCI的内陆富营养化湖泊水体蓝藻生物量比例遥感估算方法，其特征在于包括以下步骤：

S1、从OLCI遥感影像上提取所述内陆富营养化湖泊水体研究区域，经过预处理，得到所述区域对应水体的遥感反射率；所述遥感影像由像元组成，所述预处理包括辐射定标、几何校正、大气校正，所述遥感反射率第一波段至第十六波段的遥感反射率值依次记为OA1至OA16；其中，第六波段中心波长560nm、第七波段中心波长620nm、第八波段中心波长665nm、第十一波段中心波长709nm、第十六波段中心波长778nm；

S2、求像元处OA6和OA10的比值；如果所述比值大于事先设定的阈值，则认为所述像元的蓝藻生物量比例大于90％，记为极高梯度，输出结果；小于该阈值时，认为所述像元的蓝藻生物量比例小于90％，进入S3；其中，该阈值为1.77；

S3、对目标内陆富营养化湖泊水体研究区域对应水体实地测量，得到采样区域的遥感反射率，通过光学仪器测量得到浮游藻类的吸收系数a_ph在全波段的实测值，作为所述采样区域的实测值；

S4、由OA16计算得到所述像元778nm处的后向散射值b_b(778)，由OA11/OA8、b_b(778)，得到浮游藻类在665nm处的吸收系数a_ph(665)，由a_ph(665)的反演值与实测值的比值，获取γ值，由OA11/OA8、b_b(778)，以及γ值得到校正后的浮游藻类在665nm处的吸收系数a_ph(665)；由OA11/OA7、b_b(778)，得到浮游藻类在620nm处的吸收系数a_ph(620)，由a_ph(620)的反演值与实测值的比值，获取δ值，由OA11/OA7、b_b(778)，以及δ值得到校正后的浮游藻类在620nm处的吸收系数a_ph(620)；由OA11/OA8获取值，由OA7/OA11获取值，结合a_ph(620)、a_ph(665)、和分离得到藻蓝蛋白色素在620nm处的吸收系数a_pc(620)，和叶绿素a色素在665nm的吸收系数a_chl(665)；

S5、构建归一化指数：[a_pc(620)-a_chl(665)]/[a_pc(620)+a_chl(665)]，用于指征所述像元蓝藻生物量比例；

S6、将S2至S5所述的步骤逐像元作用于所述内陆富营养化湖泊水体研究区域的OLCI遥感影像上，得到S1所述区域对应水体的蓝藻生物量比例的空间分布。

2.如权利要求1所述的一种基于OLCI的内陆富营养化湖泊水体蓝藻生物量比例遥感估算方法，其特征在于，步骤S1中大气校正采用的是SeaDAS 7.3的MUMM算法，在该算法中，_Rrs(779)与_Rrs(865)的比值设定为1.9。

3.如权利要求1所述的一种基于OLCI的内陆富营养化湖泊水体蓝藻生物量比例遥感估算方法，其特征在于，步骤S1中还包括大气校正精度分析的步骤：首先将实测点的遥感反射率的实测值通过光谱响应函数响应到遥感影像的各个波段，然后将实测点的遥感反射率和OLCI对应位置的遥感反射率进行比对，得到指数对对比结果进行评价，所述指数包括均方根误差和平均相对误差；平均相对误差小于30％时，认为大气校正效果较好，将结果输出至步骤S2。

4.如权利要求1所述的一种基于OLCI的内陆富营养化湖泊水体蓝藻生物量比例遥感估算方法，其特征在于，步骤S2中所述阈值取1.77。

5.如权利要求1所述的一种基于OLCI的内陆富营养化湖泊水体蓝藻生物量比例遥感估算方法，其特征在于，步骤S4中γ的取值为0.85，δ的取值为0.96。