一种适用于卫星高光谱数据的噪声去除方法与流程

本发明属于高光谱遥感图像处理，具体涉及一种适用于卫星高光谱数据的噪声去除方法。

背景技术：

1、国产卫星高光谱数据波段多，高分5号卫星有330个波段，资源卫星02d有166个波段。在利用高光谱卫星数据提取信息的过程中，容易受到噪声的影响，可能会对提取结果产生误判，需要同时对这些噪声去除其对数据处理的干扰。

2、原来的常用噪声去除方法主要是单波段灰度插值，利用两侧的均值去除填充噪声，这种方法应用与高光谱卫星图像效率低，实用性差，

3、因此，需要开发新的方法快速高效全面去除噪声的影响。

技术实现思路

1、本发明的目的在于提供一种适用于卫星高光谱数据的噪声去除方法，该方法能够弥补原有的技术方法存在的不足，降低噪声的影响，提高数据提取的精度和效率。

2、实现本发明目的的技术方案：

3、一种适用于卫星高光谱数据的噪声去除方法，所述方法包括：

4、步骤1，获取国产高光谱卫星影像数据，对影像做大气校正，获取地面反射率的影像数据；

5、步骤2，对影像数据进行重采样，提取波段在高分5号卫星或者资源02d卫星1956-2462nm波段范围内波段影像；

6、步骤3，对影像数据进行重采样，获得高分5号卫星或资源02d卫星在1956-2462nm波段范围内的三个相邻波段n1、n2、n3的像元值；

7、步骤4，根据波段n1、n2、n3的像元值，获得影像像元值h1和h2；

8、步骤5，根据影像像元值h1或h2，获得影像像元值h。

9、所述步骤3中波段n1、n2、n3为星在1956-2462nm波段范围内任意连续的三个波段。

10、所述步骤4具体为：将波段n1的像元值除以波段n3的像元值，获得影像像元值h1；将波段n1的像元值除以波段n2的像元值，获得影像像元值h2。

11、所述步骤5具体为：

12、当影像像元值h1大于1.3，判断该像元满足条件一；

13、当影像像元值h1小于0.7，判断该像元满足条件二；

14、当影像像元值h2大于1.2，判断该像元满足条件三；

15、当影像像元值h2小于0.8，判断该像元满足条件四；

16、选择满足条件一、条件二、条件三、条件四之一的结果，获得影像像元值h。

17、所述步骤5中获取影像像元值h的方法为：

18、满足条件一、条件二、条件三、条件四四个条件中的一个条件，则该像元为噪声像元或者像元值为1，如果四个条件都不满足，则像元值为0或者为正常像元。

19、所述步骤5中将得到的噪声像元生成掩膜，叠加到矿物信息提取结果图，使噪声区域在最终结果图上不显示。

20、所述步骤5中将正常像元的值变成1，与最终结果相乘，去除噪声的影响。

21、本发明的有益技术效果在于：

22、1、本发明提供的一种适用于卫星高光谱数据的噪声去除方法根据噪声的特点，利用波段比值判别进行计算，可以快速去除噪声，具有速度快、精度高的特点。对国产高光谱卫星数据的快速处理具有较好的作用和意义。

23、2、本发明提供的一种适用于卫星高光谱数据的噪声去除方法对条带噪声效果明显，可以有效去除条带噪声的影响。

24、3、本发明提供的一种适用于卫星高光谱数据的噪声去除方法利用三个连续波段，可以避免两个波段的弊端，减少误差和遗漏。在一定波段范围内选择这三个相邻波段，可以减少水汽等干扰波段的影响，避免数据噪声的扩大化造成有用数据的丢失。

技术特征：

1.一种适用于卫星高光谱数据的噪声去除方法，其特征在于，所述方法包括：

2.根据权利要求1所述的一种适用于卫星高光谱数据的噪声去除方法，其特征在于，所述步骤3中波段n1、n2、n3为星在1956-2462nm波段范围内任意连续的三个波段。

3.根据权利要求1所述的一种适用于卫星高光谱数据的噪声去除方法，其特征在于，所述步骤4具体为：将波段n1的像元值除以波段n3的像元值，获得影像像元值h1；将波段n1的像元值除以波段n2的像元值，获得影像像元值h2。

4.根据权利要求3所述的一种适用于卫星高光谱数据的噪声去除方法，其特征在于，所述步骤5具体为：

5.根据权利要求4所述的一种适用于卫星高光谱数据的噪声去除方法，其特征在于，所述步骤5中获取影像像元值h的方法为：

6.根据权利要求5所述的一种适用于卫星高光谱数据的噪声去除方法，其特征在于，所述步骤5中将得到的噪声像元生成掩膜，叠加到矿物信息提取结果图，使噪声区域在最终结果图上不显示。

7.根据权利要求5所述的一种适用于卫星高光谱数据的噪声去除方法，其特征在于，所述步骤5中将正常像元的值变成1，与最终结果相乘，去除噪声的影响。

技术总结
本发明属于高光谱遥感图像处理技术领域，具体涉及一种适用于卫星高光谱数据的噪声去除方法，包括：步骤1，获取国产高光谱卫星影像数据，对影像做大气校正，获取地面反射率的影像数据；步骤2，对影像数据进行重采样，提取波段在高分5号卫星或者资源02D卫星1956‑2462nm波段范围内波段影像；步骤3，对影像数据进行重采样，获得高分5号卫星或资源02D卫星在1956‑2462nm波段范围内的三个相邻波段N1、N2、N3的像元值；步骤4，根据波段N1、N2、N3的像元值，获得影像像元值H1和H2；步骤5，根据影像像元值H1或H2，获得影像像元值H。本发明方法能够弥补原有的技术方法存在的不足，降低噪声的影响，提高数据提取的精度和效率。

技术研发人员：杨燕杰,叶发旺,李新春,邱骏挺
受保护的技术使用者：核工业北京地质研究院
技术研发日：
技术公布日：2024/1/15