[发明专利]叶绿素a浓度插值方法、装置、电子设备及存储介质

说明书

本发明提供了一种叶绿素a浓度插值方法，应用于测绘技术领域，包括：获取待插值的叶绿素a浓度图像，对行数据和列数据进行EMD分解得到行数据和列数据的固有模态函数分量和余量，并进行插值计算，得到行数据插值后的固有模态函数分量和列数据插值后的固有模态函数分量，将行数据插值后的固有模态函数分量和行数据的余量相加，得到行数据插值后的信号数据，以及，将列数据插值后的固有模态函数分量和列数据的余量相加，得到列数据插值后的信号数据，计算行数据和列数据插值后的信号数据的平均值，基于行数据和列数据插值后的信号数据的平均值，得到插值后的叶绿素a浓度图像。本申请还公开了一种叶绿素a浓度插值装置、电子设备及存储介质。

技术领域

本申请涉及测绘技术领域，尤其涉及一种叶绿素a浓度插值方法、装置、电子设备及存储介质。

背景技术

中分辨率成像光谱仪(MODIS，Moderate-resolution ImagingSpectroradiometer)是美国地球观测系统(EOS，Earth Observation System)计划中用于观测全球生物和物理过程的遥感仪器，搭载在Terra和Aqua两颗卫星上，MODIS具有36个光谱波段(0.4μm-14.4μm)，分辨率为250m、500m和1000m，每1-2天对地球表面观测一次，获取陆地表面、海洋水色、云、生物地球化学、气溶胶、大气温度、浮游植物等目标的图像，且能够免费获取。但由于全球叶绿素a浓度产品数据存在大面积缺失情况，数据可利用率很低。

发明内容

本申请的主要目的在于提供一种叶绿素a浓度插值方法、装置、电子设备及存储介质，可有效提高叶绿素a浓度产品利用率。

为实现上述目的，本申请实施例第一方面提供一种叶绿素a浓度插值方法，包括：

获取待插值的叶绿素a浓度图像，所述叶绿素a浓度图像包括行数据以及列数据；

对所述行数据进行EMD分解得到行数据的固有模态函数分量和余量，以及，对所述列数据进行EMD分解得到列数据的固有模态函数分量和余量；

对所述行数据的固有模态函数分量，以及，所述列数据的固有模态函数分量利用三次样条插值法进行插值计算，得到所述行数据插值后的固有模态函数分量，以及，所述列数据插值后的固有模态函数分量；

将所述行数据插值后的固有模态函数分量和所述行数据的余量相加，得到所述行数据插值后的信号数据，以及，将所述列数据插值后的固有模态函数分量和所述列数据的余量相加，得到所述列数据插值后的信号数据；

计算所述行数据插值后的信号数据的平均值，以及，所述列数据插值后的信号数据的平均值；

基于所述行数据插值后的信号数据的平均值，以及，所述列数据插值后的信号数据的平均值，得到插值后的叶绿素a浓度图像。

在一实施例中，所述获取待插值的叶绿素a浓度图像包括：

获取待插值的叶绿素a浓度的MODIS数据，以及，与所述对象同一时间的GOCI遥感影像；

将所述GOCI遥感影像反演为叶绿素a浓度的数据；

将所述叶绿素a浓度的数据的分辨率转化成与所述对象相同的分辨率，得到转化后的叶绿素a浓度的数据；

将所述转化后的叶绿素a浓度的MODIS数据与所述GOCI遥感影像融合成所述待插值的叶绿素a浓度图像。

在一实施例中，对于每一行行数据，采用第一方式对所述行数据进行EMD分解得到行数据的固有模态函数分量和余量；

所述第一方式包括：

S1、查找所述行数据中的极大值点和极小值点，拟合出上包络线和下包络线；

S2、计算所述行数据与所述行数据的上包络线和下包络线的平均值的差值；