一种投影滤波式快速光谱去噪方法技术

本发明专利技术提出一种投影滤波式快速光谱去噪方法。方法包括：(1)将三维噪声光谱数据沿光谱维度加权投影，形成低噪声投影图像；(2)利用二维数据去噪方法进一步抑制投影图像的噪声，得到高信噪比参考图像；(3)以参考图像为先验，对光谱各通道分别进行滤波，得到低噪声光谱数据。本发明专利技术充分利用光谱数据的稀疏特性，将去噪过程应用于二维投影图像，能够同时具备高去噪性能和低计算复杂度。同时，本发明专利技术可根据实际应用效果，灵活选用最佳的二维投影图像去噪算法。

【技术实现步骤摘要】

本专利技术涉及多维信号(图像)处理领域，特别涉及一种投影滤波式光谱去噪方法。
技术介绍
真实世界光线通常具有很宽的光谱范围。尽管传统的彩色图像能够满足人类视觉系统需求，但是光谱图像可以提供更多的光谱维度细节信息，比如说精细农业、遥感成像、安全监控等。近年来出现了一些光谱图像获取方法，比如说扫描式、滤波式、编码光圈式、断层扫描式、棱镜掩膜式等。与传统相机相比，这些多光谱采集系统通过单束光线分光到不同成像单元以得到高维度光谱信息，因此光谱分辨率的提升不可避免地带来成像信噪比的降低。现有光谱去噪方法主要分为三类：1.通过把光谱数据看成每个波段图像的叠加，直接应用二维图像方法进行光谱去噪，比如K-SVD(K-SVD：K-Singular Value Decomposition)和BM3D(Block-Matching and 3D Filtering)。这类方法的不足在于，没有充分利用光谱数据不同波段间的内在相关性。2.通过把光谱数据看成三维张量，利用张量方法进行光谱去噪，比如PARAFAC(Parallel Factor Analysis)和TDL(Decomposable Nonlocal Tensor Dictionary Learning)。3.通过把光谱数据看成三维立体，利用基于三维立体的方法进行光谱去噪，比如BM4D(Block-Matching and 4D Filtering)。然而，随着光谱分辨率的提升，这些基于张量和三维立体方法的计算复杂度呈指数级上涨。
技术实现思路
本专利技术目的是旨在提出一种投影滤波式快速光谱去噪方法，解决上述现有方法无法兼有高去噪性能和低计算复杂度的缺陷。本专利技术采用的技术方案如下：一种投影滤波式快速光谱去噪方法，包括以下步骤：(1)将三维噪声光谱数据沿光谱维度加权投影，形成低噪声投影图像；(2)利用二维数据去噪方法进一步抑制投影图像的噪声，得到高信噪比参考图像；(3)以参考图像为先验，对光谱各通道分别进行滤波，得到低噪声光谱数据。本专利技术充分利用光谱数据的稀疏特性，能够同时具备高去噪性能和低计算复杂度的优点。具体地，步骤(1)的加权投影模型，能够保持每个光谱通道的结构细节，具备O(1)的计算复杂度，并且可根据噪声模型选择投影权数以获得最佳降噪效果；步骤(2)可基于实际需求，灵活选用最佳的二维降噪方法以获得最佳的降噪性能；步骤(3)的滤波算法，对于光谱各通道之间相互独立，支持并行化编程，具备高时间效率。总体来说，本专利技术的计算时间比现有技术低2-3个数量级，同时能够取得高于现有技术的降噪效果。附图说明图1为本专利技术投影滤波式快速光谱去噪方法的流程图；图2为本专利技术一个实施例的投影滤波式快速光谱去噪过程示意图。图3为本专利技术实施例与现有技术去噪方法的去噪数据对比；图4为本专利技术实施例与现有技术去噪方法的去噪图像效果对比，(a)原图像，(b)噪声图像，(c)BM3D方法，(d)BM4D方法，(e)PARAFAC方法，(f)TDL方法，(g)本专利技术方法；图5为本专利技术实施例与现有技术去噪方法的计算时间对比。具体实施方式图1是根据本专利技术投影滤波式快速光谱去噪方法的流程图，主要包括以下步骤：步骤一，将三维噪声光谱数据沿光谱维度加权投影，形成低噪声投影图像，三维噪声光谱数据包括二维空间信息和一维光谱信息。具体地，假设所求低噪声投影图像为Z，Z是M×N维数据；噪声光谱数据为Y＝[Y1，Y2，Y3...，YB]，Y是M×N×B维数据，B为光谱通道数，其中Ym(m＝1，2，3...，B)表示单波段图像；投影权数表示为向量w＝[w1，w2，w3，…，wB]。加权投影公式如下：Z＝w·Y＝w1Y1+w2Y2+w3Y3+…+wBYB该公式的突出优点为：(1)光谱维度投影，能够保持每个光谱通道的结构细节，(2)计算复杂度低，(3)可根据不同的噪声模型，灵活选择权数w，以获得最佳的降噪效果。各光谱通道投影权数由噪声模型决定，本专利技术的一个实施例中，对于泊松噪声模型，例如低照度成像，各通道投影权数相等，推导过程如下：光谱信号X与含噪声光谱数据Y的关系为：Y～Possion(X)峰值信噪比数学定义如下： S N R ( Z ) = Z ‾ V a r ( Z ) ]]>是信号的期望(均值)，Var(Z)是信号的标准差。因为Y服从独立泊松分布，Z的期望和方差具有如下性质： Z ‾ = Σ m = 1 B w m Y m ‾ = Σ m = 1 B w m Y m ‾ ]]> V a r ( Σ m = 1 B z ) = V a r ( Σ m = 1 B w m Y m ) = Σ m = 1 B w m 2 Y m V a r ( Y m ) = Σ m = 1 B w m 2 本文档来自技高网...

【技术保护点】
一种投影滤波式快速光谱去噪方法，其特征在于，包括以下步骤：(1)将三维噪声光谱数据沿光谱维度加权投影，形成低噪声投影图像；(2)利用二维数据去噪方法进一步抑制投影图像的噪声，得到高信噪比参考图像；(3)以参考图像为先验，对光谱各通道分别进行滤波，得到低噪声光谱数据。

【技术特征摘要】
1.一种投影滤波式快速光谱去噪方法，其特征在于，包括以下步骤：(1)将三维噪声光谱数据沿光谱维度加权投影，形成低噪声投影图像；(2)利用二维数据去噪方法进一步抑制投影图像的噪声，得到高信噪比参考图像；(3)以参考图像为先验，对光谱各通道分别进行滤波，得到低噪声光谱数据。2.如权利要求1所述的一种投影滤波式快速光谱去噪方法，其特征在于，所述三维噪声光谱数据包括二维空间信息和一维光谱信息。3.如权利要求1所述的一种投影滤波式快速光谱去噪方法，其特征在于，所述步骤(1)中，加权投影的数学公式为：Z＝w·Y＝w1Y1+w2Y2+w3Y3+…+wBYB其中，Z为...

【专利技术属性】
技术研发人员：曹汛，陈林森，岳涛，
申请(专利权)人：南京大学，
类型：发明
国别省市：江苏;32