【技术实现步骤摘要】
高光谱卫星的水体叶绿素a浓度反演方法及系统
[0001]本公开涉及高光谱卫星
，尤其涉及一种高光谱卫星的水体叶绿素a浓度反演方法及系统。
技术介绍
[0002]叶绿素a浓度是水环境监测的重要指标。目前叶绿素a浓度的监测方法有现场直接测量和遥感间接测量。现场直接测量是通过实地采集水样后低温保存，之后在实验室内进行分析进而得到采样点水体的叶绿素a浓度。这种方法由于人工测量的采样点和采样次数是有限的，不能全面评价水体叶绿素a的空间分布情况。遥感间接测量是通过水体光谱特征与叶绿素a浓度的关系，寻找两者之间的关系进行建模反演得到后者的结果。
[0003]由于中国卫星高光谱数据与SeaWIFS、MODIS、MERIS等国外水色卫星数据在数据质量(遥感反射率精度)和载荷指标(波段设置、空间分辨率)等方面的差异，现有技术直接应用于国产高光谱卫星无法获得较好的叶绿素a浓度反演效果。
技术实现思路
[0004]为克服相关技术中存在的问题，本公开提供一种高光谱卫星的水体叶绿素a浓度反演方法及系统。
[0005...
【技术保护点】
【技术特征摘要】
1.一种高光谱卫星的水体叶绿素a浓度反演方法，其特征在于，包括：获取湖库水体本底数据集中采样点的遥感反射率数据，并计算所述采样点在高光谱卫星数据的预设波段的等效遥感反射率；根据所述采样点的等效遥感反射率，分别确定不同水体类型在各预设波段的平均遥感反射率光谱，其中，所述水体类型包括：贫营养水体类型、中营养水体类型和富营养水体类型；确定每个水体类型对应的叶绿素a反演模型；分别计算各个水体类型对应的平均遥感反射率光谱与高光谱影像的水体像元的遥感反射率在所述预设波段的光谱角余弦，以确定所述高光谱影像的水体像元所属的目标水体类型；采用与所述目标水体类型对应的叶绿素a反演模型，确定所述高光谱影像的水体像元对应的目标叶绿素a浓度。2.根据权利要求1所述的高光谱卫星的水体叶绿素a浓度反演方法，其特征在于，采用以下第一计算公式计算所述采样点在高光谱卫星数据波段的等效遥感反射率：其中，表示所述水体在高光谱卫星数据的预设波段的等效遥感反射率，表示所述水体的遥感反射率数据，f
i
(λ)表示第i个高光谱卫星数据的预设波段的光谱响应函数，λ表示波长，λ
min
和λ
max
分别表示f
i
(λ)函数的最小波长和最大波长；采用以下第二计算公式计算不同水体类型在各预设波段的平均遥感反射率光谱：其中，表示任一水体类型在第i个高光谱卫星数据的预设波段的平均遥感反射率，N是该水体类型的湖库水体本底数据集的样本个数。3.根据权利要求1所述的高光谱卫星的水体叶绿素a浓度反演方法，其特征在于，确定每个水体类型对应的叶绿素a反演模型，包括：获取湖库水体本底数据集中各个水体类型采样点的水质数据，其中，所述水质数据包括叶绿素a浓度数据；对各水体类型，根据归属该类采样点的等效遥感反射率和对应的叶绿素a浓度数据，确定每个水体类型对应的叶绿素a反演模型。4.根据权利要求1所述的高光谱卫星的水体叶绿素a浓度反演方法，其特征在于，确定所述高光谱影像的水体像元所属的目标水体类型，包括：当任一水体类型对应的平均遥感反射率光谱与高光谱影像的水体像元的遥感反射率在所述预设波段的光谱角余弦大于或者等于预设值时，确定所述高光谱影像的水体像元所属的目标水体类型为该任一水体类型；当所述高光谱影像的水体像元的遥感反射率与每个水体类型对应的平均遥感反射率光谱在所述预设波段的光谱角余弦均小于预设值时，分别计算所述高光谱影像的水体像元
相对于每个水体类型的权重值。5.根据权利要求4所述的高光谱卫星的水体叶绿素a浓度反演方法，其特征在于，采用以下第三计算公式计算所述高光谱影像的水体像元相对于贫营养水体类型的第一权重值：采用以下第四计算公式计算所述高光谱影像的水体像元相对于中营养水体类型的第二权重值：采用以下第五计算公式计算所述高光谱影像的水体像元相对于富营养水体类型的第二权重值：其中，SA1表示所述高光谱影像的水体像元的遥感反射率与贫营养水体类型对应的平均遥感反射率光谱在所述预设波段的光谱角余弦，SA2表示所述高光谱影像的水体像元的遥感反射率与中营养水体类型对应的平均遥感反射率光谱在所述预设波段的光谱角余弦，SA3表示所述高光谱影...
【专利技术属性】
技术研发人员：刘瑶，刘智伟，肖晨超，
申请(专利权)人：自然资源部国土卫星遥感应用中心，
类型：发明
国别省市：
还没有人留言评论。发表了对其他浏览者有用的留言会获得科技券。