一种高光谱图像滤波降噪的方法技术

本发明专利技术涉及高光谱图像技术领域，具体公开了一种高光谱图像滤波降噪的方法，其中，包括：获取待滤波的高光谱图像的图像信息，所述图像信息包括图像平面坐标轴和波段；根据所述图像信息选择匹配的三维滤波窗口，并确定与所述三维滤波窗口对应的拟合多项式；根据所述三维滤波窗口遍历所述待滤波的高光谱图像，并通过所述拟合多项式实现对所述待滤波的高光谱图像内的所有数据点进行拟合；根据拟合后的数据生成新的高光谱图像，所述新的高光谱图像的空间大小与所述待滤波的高光谱图像一致。本发明专利技术提供的高光谱图像滤波降噪的方法针对高光谱图像平滑效果更加准确。像平滑效果更加准确。像平滑效果更加准确。

【技术实现步骤摘要】
一种高光谱图像滤波降噪的方法


[0001]本专利技术涉及高光谱图像
，尤其涉及一种高光谱图像滤波降噪的方法。

技术介绍

[0002]高光谱图像是由光学传感器对同一场景利用几十甚至上百个光谱波段获得的图像。高光谱数据在成像、传输的过程中会受到多种噪声源的污染，噪声对高光谱图像的分类识别带来了干扰。光谱图像中的噪声分布复杂，降噪滤波是制约高光谱图像分析精度的难题。高光谱图像噪声滤除的主要方法有直方图匹配、矩匹配、光谱微分、S
‑
G(Savitzky
‑
Golay，平滑去噪)滤波等方法，其中S
‑
G滤波被广泛应用于高光谱数据降噪滤波。
[0003]目前的S
‑
G滤波方法对一维曲线的降噪平滑效果非常出色，但在高光谱数据中，由于空间信息无法利用，遇到个别噪声特别严重的波段，原有的S
‑
G滤波方法仍不能达到较好的处理效果，并且所带来的信息丢失会导致图像失真更为严重。

技术实现思路

[0004]本专利技术提供了一种高光谱图像滤波降噪的方法，解决相关技术中存在的S
‑
G滤波方法对于高光谱图像滤波降噪效果不好的问题。
[0005]作为本专利技术的一个方面，提供一种高光谱图像滤波降噪的方法，其中，包括：
[0006]获取待滤波的高光谱图像的图像信息，所述图像信息包括图像平面坐标轴和波段；
[0007]根据所述图像信息选择匹配的三维滤波窗口，并确定与所述三维滤波窗口对应的拟合多项式；
[0008]根据所述三维滤波窗口遍历所述待滤波的高光谱图像，并通过所述拟合多项式实现对所述待滤波的高光谱图像内的所有数据点进行拟合；
[0009]根据拟合后的数据生成新的高光谱图像，所述新的高光谱图像的空间大小与所述待滤波的高光谱图像一致。
[0010]进一步地，所述根据所述图像信息选择匹配的三维滤波窗口，并确定与所述三维滤波窗口对应的拟合多项式，包括：
[0011]设定三维滤波窗口的窗口大小m值，以及设定拟合多项式的阶数k值，其中2m＞k，m为大于或者等于1的自然数；
[0012]根据所述m值和k值通过最小二乘法确定拟合多项式的拟合参数；
[0013]根据所述拟合参数计算得到该拟合参数下的三维滤波窗口。
[0014]进一步地，所述图像信息包括图像平面坐标轴(x，y)和波段z，所述三维滤波窗口为一组以(x，y，z)＝(0，0，0)为中心且三个方向上均为2m+1个数据的立方体。
[0015]进一步地，所述根据所述m值和k值通过最小二乘法确定拟合多项式的拟合参数，包括：
[0016]根据所述m值和k值确定拟合多项式的表达式；
[0017]根据所述拟合多项式的表达式构成线性方程组；
[0018]对所述线性方程组进行最小二乘法拟合，确定拟合参数。
[0019]进一步地，所述根据所述拟合参数计算得到该拟合参数下的三维滤波窗口，包括：
[0020]根据所述拟合参数求解最佳系数矩阵；
[0021]对所述最佳系数矩阵求解得到三维滤波窗口的表达式。
[0022]进一步地，所述拟合多项式的表达式为：
[0023][0024]进一步地，所述根据所述三维滤波窗口遍历所述待滤波的高光谱图像，并通过所述拟合多项式实现对所述待滤波的高光谱图像内的所有数据点进行拟合，包括：
[0025]对所述待滤波的高光谱图像进行扩充预处理；
[0026]遍历扩充预处理后的待滤波的高光谱图像，并与所述三维滤波窗口进行运算得到所述待滤波的高光谱图像内的所有数据点的拟合值。
[0027]进一步地，所述对所述待滤波的高光谱图像进行扩充预处理，包括：
[0028]对所述待滤波的高光谱图像的各维度边缘进行扩充，分别各自沿x轴和y轴的正负方向扩充大小为m个像素，z轴正负方向分别扩充m个波段，扩充的数据值以其所在原有长方体所在面为对称轴成对称关系。
[0029]本专利技术提供的高光谱图像滤波降噪的方法，通过选择匹配的三维滤波窗口遍历待滤波的高光谱图像，并通过拟合多项式对待滤波的高光谱图像内的所有数据点进行拟合，得到新的高光谱图像，从而实现对待滤波的高光谱图像的滤波降噪，且该方法由于考虑高光谱图像的空间信息，且考虑同一波段相邻像素点的相似性，平滑效果更加准确。另外，还可以通过设置不同的m值和k值，对高光谱空间滤波进行优化，对于噪声严重的波段，在达到平滑的同时，可以更好的保留图像的细节。
附图说明
[0030]附图是用来提供对本专利技术的进一步理解，并且构成说明书的一部分，与下面的具体实施方式一起用于解释本专利技术，但并不构成对本专利技术的限制。在附图中：
[0031]图1为本专利技术提供的高光谱图像滤波降噪的方法的流程图。
[0032]图2为本专利技术提供的三维滤波窗口的示意图。
具体实施方式
[0033]需要说明的是，在不冲突的情况下，本专利技术中的实施例及实施例中的特征可以相互结合。下面将参考附图并结合实施例来详细说明本专利技术。
[0034]为了使本领域技术人员更好地理解本专利技术方案，下面将结合本专利技术实施例中的附图，对本专利技术实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本专利技术一部分的实施例，而不是全部的实施例。基于本专利技术中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都应当属于本专利技术保护的范围。
[0035]需要说明的是，本专利技术的说明书和权利要求书及上述附图中的术语“第一”、“第二”等是用于区别类似的对象，而不必用于描述特定的顺序或先后次序。应该理解这样使用
的数据在适当情况下可以互换，以便这里描述的本专利技术的实施例。此外，术语“包括”和“具有”以及他们的任何变形，意图在于覆盖不排他的包括，例如，包含了一系列步骤或单元的过程、方法、系统、产品或设备不必限于清楚地列出的那些步骤或单元，而是可包括没有清楚地列出的或对于这些过程、方法、产品或设备固有的其它步骤或单元。
[0036]在本实施例中提供了一种高光谱图像滤波降噪的方法，图1是根据本专利技术实施例提供的高光谱图像滤波降噪的方法的流程图，如图1所示，包括：
[0037]S110、获取待滤波的高光谱图像的图像信息，所述图像信息包括图像平面坐标轴和波段；
[0038]需要说明的是，对于相同物质或者物体，相邻像素点的光谱图基本一致，而相邻波段的DN值(Digital Number，遥感影像像元亮度值)又十分接近，因此可以考虑三个维度的平滑，图像平面坐标轴(x,y)和波段z。
[0039]S120、根据所述图像信息选择匹配的三维滤波窗口，并确定与所述三维滤波窗口对应的拟合多项式；
[0040]在本专利技术实施例中，具体可以包括：
[0041]设定三维滤波窗口的窗口大小m值，以及设定拟合多项式的阶数k值，其中2m＞k，m为大于或者等于1的自然数；
[0042]根本文档来自技高网...

【技术保护点】

【技术特征摘要】
1.一种高光谱图像滤波降噪的方法，其特征在于，包括：获取待滤波的高光谱图像的图像信息，所述图像信息包括图像平面坐标轴和波段；根据所述图像信息选择匹配的三维滤波窗口，并确定与所述三维滤波窗口对应的拟合多项式；根据所述三维滤波窗口遍历所述待滤波的高光谱图像，并通过所述拟合多项式实现对所述待滤波的高光谱图像内的所有数据点进行拟合；根据拟合后的数据生成新的高光谱图像，所述新的高光谱图像的空间大小与所述待滤波的高光谱图像一致。2.根据权利要求1所述的高光谱图像滤波降噪的方法，其特征在于，所述根据所述图像信息选择匹配的三维滤波窗口，并确定与所述三维滤波窗口对应的拟合多项式，包括：设定三维滤波窗口的窗口大小m值，以及设定拟合多项式的阶数k值，其中2m＞k，m为大于或者等于1的自然数；根据所述m值和k值通过最小二乘法确定拟合多项式的拟合参数；根据所述拟合参数计算得到该拟合参数下的三维滤波窗口。3.根据权利要求2所述的高光谱图像滤波降噪的方法，其特征在于，所述图像信息包括图像平面坐标轴(x，y)和波段z，所述三维滤波窗口为一组以(x，y，z)＝(0，0，0)为中心且三个方向上均为2m+1个数据的立方体。4.根据权利要求3所述的高光谱图像滤波降噪的方法，其特征在于，所述根据所述m值和k值通过最小二乘法确定拟合多项式的...

【专利技术属性】
技术研发人员：李志峰，
申请(专利权)人：江苏云脑数据科技有限公司，
类型：发明
国别省市：