CN106780353A - 一种基于时相光谱角度量的多时相云遮挡数据恢复方法 - Google Patents

Abstract

一种基于时相光谱角度量的多时相云遮挡数据恢复方法，本发明涉及基于时相光谱角度量的多时相云遮挡数据恢复方法。本发明的目的是为了解决现有云遮挡数据恢复方法恢复精度低的问题。具体过程为：一、输入存在云遮挡的多时相多光谱遥感图像；二、计算存在云遮挡的多时相多光谱遥感图像地理坐标上各点的云遮挡数据缺失度；三、将缺失度最小的点作为要填充的点；四、提取X和相应缺失时相地理坐标；五、利用时相光谱角函数计算X与除X点外的其他各点相似度；六、找到相似度最大点Y；七、利用Y和缺失数据填补算法对要填充点的多时相多光谱数据X进行填补；八、迭代一，直至全部云遮挡数据被填充。本发明用于云遮挡数据恢复领域。

Description

一种基于时相光谱角度量的多时相云遮挡数据恢复方法

技术领域

本发明涉及基于时相光谱角度量的多时相云遮挡数据恢复方法。

背景技术

遥感图像的可靠性和完整性是通过遥感图像实现对社会服务的基础，然而由于不可避免的云遮挡情况的出现，导致很多情况下遥感图像数据存在了云霞观测数据的缺失，降低了遥感图像的分析的适用范围。实现对遥感图像中云遮挡数据的恢复，可有效扩大遥感图像应用范畴，提升卫星光学图像的应用效率，具有十分重要的意义。

当前遥感图像云遮挡区域的数据恢复，主要分为两大类：一是仅利用当前图像中缺失数据周围的一些未缺失数据进行线性预测和填补，但这种情况在云遮挡面积较大的恢复的精度非常低；第二种情况是利用一幅相同区域相同传感器在其他时刻获得的完全未缺失数据进行填补，但是，由于光谱漂移的存在，预测精度仍然不能满足需求，同时完全没有云遮挡的情况很难满足，因此这种恢复算法的适用性不够。

发明内容

本发明的目的是为了解决现有云遮挡数据恢复方法恢复精度低的问题，而提出一种基于时相光谱角度量的多时相云遮挡数据恢复方法。

一种基于时相光谱角度量的多时相云遮挡数据恢复方法具体过程为：

步骤一、输入存在云遮挡的多时相多光谱遥感图像；

步骤二、计算存在云遮挡的多时相多光谱遥感图像地理坐标上各点的云遮挡数据缺失度；

步骤三、将存在云遮挡的多时相多光谱遥感图像空间上各点的云遮挡数据缺失度最小的点作为要填充的点；

步骤四、提取要填充点的多时相多光谱数据X和相应缺失时相地理坐标；

步骤五、利用时相光谱角函数计算要填充点的多时相多光谱数据X与除X点外的其他各点相似度；

步骤六、找到相似度最大点Y；

步骤七、利用Y和缺失数据填补算法对要填充点的多时相多光谱数据X进行填补；

步骤八、迭代步骤一，直至全部云遮挡数据被填充。

本发明的有益效果为：

本发明提出了一种新的基于多时相信息联合利用的遥感图像云遮挡数据恢复算法，改善了传统遥感光学图像云遮挡数据恢复精度和可恢复区域的范围，多个时相数据联合利用，提高了遥感光学图像的应用范围和应用能力，提高云遮挡数据恢复的适用范围，提升遥感光学图像的数据完整性，提高光学遥感图像的应用深度。本方法相对经典的NSPI方法，其各波段平均的平方根误差由39.99缩小到19.45。

为了验证本发明所提出的方法的性能，针对一组Landsat8多光谱遥感图像进行了测试，该组数据共9个时相，其中各个时相中均不存在云遮挡位置，然后我们通过将其他图像中真实的云位置引入这些图像中，实现了对云遮挡的仿真。各个时相中云量均为30％左右，云的位置随机不固定。图2a为原始的第2个时相中的RGB三波段的原始图像数据，图2b为对图2a进行的云遮挡仿真，图2c为通过经典的云遮挡数据恢复算法(NSPI方法)进行的数据恢复结果，仍然显示的是第二个时相的恢复图，图2d为本文方法进行的云遮挡数据恢复，也为第二时相的恢复结果。通过对比图2c和图2d，我们可以清晰的看到，传统方法仍然有部分数据不能恢复(黑色部分)，而我们的方法能够恢复所有数据且恢复的数据与未缺失图像(图2a)非常近似，恢复效果非常优秀。图3a为利用NSPI方法对第二个时相的近红外波段的恢复结果与原始未缺失恢复结果的散点图。图3b为利用本文方法对第二个时相的近红外波段的恢复结果与原始未缺失恢复结果的散点图。好的恢复算法其恢复结果应该更接近原始的未缺失数据，好的恢复算法的散点图应该更接近于对角线。对比图3a和图3b，可以看到我们方法的散点图要比NSPI方法更接近对角线，即我们的方法的恢复结果要比NSPI方法好。

附图说明

图1为基于时相光谱角度量的多时相数据恢复算法流程图；

图2a为一幅青岛地区不存在云遮挡的landsat5多光谱图像，landsat5为卫星传感器；

图2b为对图2a中进行仿真云缺失的图像(黑色为缺失部分)；

图2c为通过经典的NSPI方法对图2b进行遮挡数据恢复的结果，NSPI为neighborhood similar pixel interpolation；

图2d为通过本算法对图2b进行遮挡数据恢复的结果；

图3a为利用NSPI方法对第2个时相NIR波段进行恢复的散点图，NIR为近红外，NIRband of origin image为原始光谱近红外波段(第2个时相)，NIR band of filled image为恢复光谱近红外波段(第2个时相)；

图3b为利用本文方法对第2个时相NIR波段进行恢复的散点图，NIR band oforigin image为原始光谱近红外波段(第2个时相)，NIR band of filled image为恢复光谱近红外波段(第2个时相)。

具体实施方式

具体实施方式一：结合图1说明本实施方式，本实施方式的一种基于时相光谱角度量的多时相云遮挡数据恢复方法具体过程为：

步骤一、输入存在云遮挡的多时相多光谱遥感图像；

步骤二、计算存在云遮挡的多时相多光谱遥感图像地理坐标上各点的云遮挡数据缺失度；

步骤三、将存在云遮挡的多时相多光谱遥感图像空间上各点的云遮挡数据缺失度最小的点作为要填充的点；

步骤四、提取要填充点的多时相多光谱数据X和相应缺失时相地理坐标；

步骤五、利用时相光谱角函数计算要填充点的多时相多光谱数据X与除X点外的其他各点相似度；每个点都是一个多时相多光谱的数据；

步骤六、找到相似度最大点Y；

步骤七、利用Y和缺失数据填补算法对要填充点的多时相多光谱数据X进行填补；

步骤八、迭代步骤一，直至全部云遮挡数据被填充。

具体实施方式二：本实施方式与具体实施方式一不同的是：所述步骤一中输入存在云遮挡的多时相多光谱遥感图像；具体为：

输入存在云遮挡的多时相多光谱遥感图像数据空间上任意一点都为一个矩阵，其行为时相，列为光谱。

其它步骤及参数与具体实施方式一相同。

具体实施方式三：本实施方式与具体实施方式一或二不同的是：所述步骤二中计算存在云遮挡的多时相多光谱遥感图像空间上各点的云遮挡数据缺失度；具体过程为：

多时相多光谱遥感图像空间上点a的云遮挡数据缺失度计算公式如下：

其它步骤及参数与具体实施方式一或二相同。

具体实施方式四：本实施方式与具体实施方式一至三之一不同的是：所述步骤五中利用时相光谱相似度计算要填充点的多时相多光谱数据X与除X外其他各点相似度；具体过程为：

要填充点的多时相多光谱数据X与除X外其他任意点Z的时相光谱相似度计算公式如下：

式中，X'为X的转置，Z'为Z的转置。

其它步骤及参数与具体实施方式一至三之一相同。

具体实施方式五：本实施方式与具体实施方式一至四之一不同的是：所述步骤七中利用Y和预测公式对要填充点的多时相多光谱数据X进行填补；具体过程为：

步骤五一、利用X、Y及其相应缺失时相位置进行对要填充点的多时相多光谱数据X的局部缺失数据的预测，即对要填充点的多时相多光谱数据X中缺失但Y未缺失的时相数据的预测，预测公式如下：

其中X_∩和Y_∩表示X和Y中都未出现云遮挡时相的数据；当X点处第1、3、5时相未出现云遮挡，Y点第1、3、6时相未被云遮挡，则此时X_∩和Y_∩分别表示X点的第1、3时相的数据和Y点的第1、3时相的数据；|| ||₁表示矩阵的一范数；

步骤五二、利用步骤五一得到的预测数据P对要填充点X进行填补，其中填补位置仅为X缺失但Y未缺失的部分。

其它步骤及参数与具体实施方式一至四之一相同。

采用以下实施例验证本发明的有益效果：

实施例一：

本实施例一种基于时相光谱角度量的多时相云遮挡数据恢复方法具体是按照以下步骤制备的：

为了验证本发明所提出的方法的性能，针对一组Landsat8多光谱遥感图像进行了测试，该组数据共9个时相，其中各个时相中均不存在云遮挡位置，然后我们通过将其他图像中真实的云位置引入这些图像中，实现了对云遮挡的仿真。各个时相中云量均为30％左右，云的位置随机不固定。图2a为原始的第2个时相中的RGB三波段的原始图像数据，图2b为对图2a进行的云遮挡仿真，图2c为通过经典的云遮挡数据恢复算法(NSPI方法)进行的数据恢复结果，仍然显示的是第二个时相的恢复图，图2d为本文方法进行的云遮挡数据恢复，也为第二时相的恢复结果。通过对比图2c和图2d，我们可以清晰的看到，传统方法仍然有部分数据不能恢复(黑色部分)，而我们的方法能够恢复所有数据且恢复的数据与未缺失图像(图2a)非常近似，恢复效果非常优秀。图3a为利用NSPI方法对第二个时相的近红外波段的恢复结果与原始未缺失恢复结果的散点图。图3b为利用本文方法对第二个时相的近红外波段的恢复结果与原始未缺失恢复结果的散点图。好的恢复算法其恢复结果应该更接近原始的未缺失数据，好的恢复算法的散点图应该更接近于对角线。对比图3a和图3b，可以看到我们方法的散点图要比NSPI方法更接近对角线，即我们的方法的恢复结果要比NSPI方法好。

本发明还可有其它多种实施例，在不背离本发明精神及其实质的情况下，本领域技术人员当可根据本发明作出各种相应的改变和变形，但这些相应的改变和变形都应属于本发明所附的权利要求的保护范围。

Claims (5)

1.一种基于时相光谱角度量的多时相云遮挡数据恢复方法，其特征在于：一种基于时相光谱角度量的多时相云遮挡数据恢复方法具体过程为：

步骤一、输入存在云遮挡的多时相多光谱遥感图像；

步骤二、计算存在云遮挡的多时相多光谱遥感图像地理坐标上各点的云遮挡数据缺失度；

步骤三、将存在云遮挡的多时相多光谱遥感图像空间上各点的云遮挡数据缺失度最小的点作为要填充的点；

步骤四、提取要填充点的多时相多光谱数据X和相应缺失时相地理坐标；

步骤五、利用时相光谱角函数计算要填充点的多时相多光谱数据X与除X点外的其他各点相似度；

步骤六、找到相似度最大点Y；

步骤七、利用Y和缺失数据填补算法对要填充点的多时相多光谱数据X进行填补；

步骤八、迭代步骤一，直至全部云遮挡数据被填充。

2.根据权利要求1所述一种基于时相光谱角度量的多时相云遮挡数据恢复方法，其特征在于：所述步骤一中输入存在云遮挡的多时相多光谱遥感图像；具体过程为：

输入存在云遮挡的多时相多光谱遥感图像数据空间上任意一点都为一个矩阵，其行为时相，列为光谱。

3.根据权利要求2所述一种基于时相光谱角度量的多时相云遮挡数据恢复方法，其特征在于：所述步骤二中计算存在云遮挡的多时相多光谱遥感图像空间上各点的云遮挡数据缺失度；具体过程为：

多时相多光谱遥感图像空间上点a的云遮挡数据缺失度计算公式如下：

4.根据权利要求3所述一种基于时相光谱角度量的多时相云遮挡数据恢复方法，其特征在于：所述步骤五中利用时相光谱相似度计算要填充点的多时相多光谱数据X与除X外其他各点相似度；具体过程为：

要填充点的多时相多光谱数据X与除X外其他任意点Z的时相光谱相似度计算公式如下：

式中，X'为X的转置，Z'为Z的转置。

5.根据权利要求4所述一种基于时相光谱角度量的多时相云遮挡数据恢复方法，其特征在于：所述步骤七中利用Y和预测公式对要填充点的多时相多光谱数据X进行填补；具体过程为：

步骤五一、利用X、Y及其相应缺失时相位置进行对要填充点的多时相多光谱数据X的局部缺失数据的预测，即对要填充点的多时相多光谱数据X中缺失但Y未缺失的时相数据的预测，预测公式如下：

$P=Y\×\frac{\left|\right|X_{\∩}|{|}_1}{\left|\right|Y_{\∩}|{|}_1}$

其中X_∩和Y_∩表示X和Y中都未出现云遮挡时相的数据，当X点处第1、3、5时相未出现云遮挡，Y点第1、3、6时相未被云遮挡，则此时X_∩和Y_∩分别表示X点的第1、3时相的数据和Y点的第1、3时相的数据；|| ||₁表示矩阵的一范数；

步骤五二、利用步骤五一得到的预测数据P对要填充点X进行填补，其中填补位置仅为X缺失但Y未缺失的部分。