CN111369457B - 稀疏判别张量鲁棒性pca的遥感图像去噪方法 - Google Patents
稀疏判别张量鲁棒性pca的遥感图像去噪方法 Download PDFInfo
- Publication number
- CN111369457B CN111369457B CN202010129712.9A CN202010129712A CN111369457B CN 111369457 B CN111369457 B CN 111369457B CN 202010129712 A CN202010129712 A CN 202010129712A CN 111369457 B CN111369457 B CN 111369457B
- Authority
- CN
- China
- Prior art keywords
- tensor
- sparse
- denoising
- image
- remote sensing
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
- 238000000034 method Methods 0.000 title claims abstract description 91
- 238000000513 principal component analysis Methods 0.000 claims abstract description 53
- 238000012360 testing method Methods 0.000 claims abstract description 30
- 238000012549 training Methods 0.000 claims abstract description 17
- 230000000694 effects Effects 0.000 claims abstract description 15
- 238000011156 evaluation Methods 0.000 claims abstract description 6
- 238000003709 image segmentation Methods 0.000 claims abstract description 5
- 238000000354 decomposition reaction Methods 0.000 claims description 18
- 239000011159 matrix material Substances 0.000 claims description 11
- 238000012706 support-vector machine Methods 0.000 claims description 5
- 241000764238 Isis Species 0.000 claims description 2
- 230000008859 change Effects 0.000 claims description 2
- 238000005457 optimization Methods 0.000 claims description 2
- 238000012545 processing Methods 0.000 description 14
- 238000004422 calculation algorithm Methods 0.000 description 11
- 230000008569 process Effects 0.000 description 10
- 230000006870 function Effects 0.000 description 9
- 230000003595 spectral effect Effects 0.000 description 7
- 238000001228 spectrum Methods 0.000 description 6
- 230000009466 transformation Effects 0.000 description 5
- 230000009467 reduction Effects 0.000 description 4
- 230000003044 adaptive effect Effects 0.000 description 3
- 230000005540 biological transmission Effects 0.000 description 3
- 230000000875 corresponding effect Effects 0.000 description 3
- 230000007547 defect Effects 0.000 description 3
- 238000005516 engineering process Methods 0.000 description 3
- 238000001914 filtration Methods 0.000 description 3
- 238000004458 analytical method Methods 0.000 description 2
- 238000007635 classification algorithm Methods 0.000 description 2
- 238000010276 construction Methods 0.000 description 2
- 230000002596 correlated effect Effects 0.000 description 2
- 238000013461 design Methods 0.000 description 2
- 238000005315 distribution function Methods 0.000 description 2
- 238000000605 extraction Methods 0.000 description 2
- 238000003384 imaging method Methods 0.000 description 2
- 230000000007 visual effect Effects 0.000 description 2
- OAICVXFJPJFONN-UHFFFAOYSA-N Phosphorus Chemical compound [P] OAICVXFJPJFONN-UHFFFAOYSA-N 0.000 description 1
- 230000002411 adverse Effects 0.000 description 1
- 238000013459 approach Methods 0.000 description 1
- 230000000903 blocking effect Effects 0.000 description 1
- 230000015556 catabolic process Effects 0.000 description 1
- 230000006835 compression Effects 0.000 description 1
- 238000007906 compression Methods 0.000 description 1
- 238000005314 correlation function Methods 0.000 description 1
- 238000006731 degradation reaction Methods 0.000 description 1
- 238000001514 detection method Methods 0.000 description 1
- 238000009499 grossing Methods 0.000 description 1
- 238000004377 microelectronic Methods 0.000 description 1
- 238000005065 mining Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 230000008092 positive effect Effects 0.000 description 1
- 238000007781 pre-processing Methods 0.000 description 1
- 238000012847 principal component analysis method Methods 0.000 description 1
- 238000003672 processing method Methods 0.000 description 1
- 238000011084 recovery Methods 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 230000011218 segmentation Effects 0.000 description 1
- 238000000926 separation method Methods 0.000 description 1
- 238000007619 statistical method Methods 0.000 description 1
- 238000011426 transformation method Methods 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T5/00—Image enhancement or restoration
- G06T5/70—Denoising; Smoothing
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F18/00—Pattern recognition
- G06F18/20—Analysing
- G06F18/21—Design or setup of recognition systems or techniques; Extraction of features in feature space; Blind source separation
- G06F18/214—Generating training patterns; Bootstrap methods, e.g. bagging or boosting
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/10—Segmentation; Edge detection
- G06T7/11—Region-based segmentation
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/10—Image acquisition modality
- G06T2207/10032—Satellite or aerial image; Remote sensing
- G06T2207/10036—Multispectral image; Hyperspectral image
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/30—Subject of image; Context of image processing
- G06T2207/30204—Marker
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
- Y02A40/00—Adaptation technologies in agriculture, forestry, livestock or agroalimentary production
- Y02A40/10—Adaptation technologies in agriculture, forestry, livestock or agroalimentary production in agriculture
Landscapes
- Engineering & Computer Science (AREA)
- Theoretical Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Physics & Mathematics (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Data Mining & Analysis (AREA)
- Bioinformatics & Computational Biology (AREA)
- Bioinformatics & Cheminformatics (AREA)
- Evolutionary Biology (AREA)
- Evolutionary Computation (AREA)
- Artificial Intelligence (AREA)
- General Engineering & Computer Science (AREA)
- Life Sciences & Earth Sciences (AREA)
- Image Processing (AREA)
Abstract
本发明公开的一种稀疏判别张量鲁棒性PCA的遥感图像去噪方法，旨在提供一种能更有效除去高光谱图像中噪声的去噪方法，本发明通过下述技术方案实现：对输入高光谱遥感图像进行图像分割，按固定尺寸，将高光谱遥感图像分割成不重叠的张量数据块；分别对每个张量数据块进行稀疏判别张量鲁棒性主成分分析：先将张量数据块分解为低秩分量、稀疏判别分量和稀疏噪声分量，构建稀疏判别张量鲁棒性主成分分析模型，迭代求解该模型得到去噪后的张量数据块，获得去噪图像；对重构去噪后的无噪声数据按比例划分训练集和测试集，最后将训练集和测试集输入到分类器，通过分类器输出所有测试样本的类别标记，实现对高光谱遥感图像去噪效果的评估。
Description
技术领域
本发明涉及图像处理领域，尤其是遥感图像处理技术领域的一种基于稀疏判别张量鲁棒性主成分分析(PCA)的高光谱遥感图像去噪方法。
背景技术
在当代高度信息化的社会中，图形和图像在信息传播中所起的作用越来越大，图像是一种重要的信息源，通过图像处理可以帮助人们了解信息的内涵。但是图像在生成和传输过程中常常因受到各种噪声的干扰和影响而使图像降质，这对后续图像的处理(如分割、压缩和图像理解等)将产生不利影响。噪声种类很多，如：电噪声、机械噪声、信道噪声和其他噪声。一般图像的能量主要集中在低频，而图像的细节部分集中在高频区域。由于在图像的存取、数字化和传输中常伴有噪声出现，而这部分干扰信息主要集中在高频区域内，所以消除噪声的一般方法是利用低通滤波来衰减高频分量。但与之同时带来的负面影响使图像的细节也有一定的衰减，从视觉效果上来看图像比处理前模糊。如何从这些被噪声污染的高光谱影像中获取更加有效可靠的信息是亟待解决的重要问题。为了抑制噪声，改善图像质量，便于更高层次的处理，必须对图像进行去噪预处理。消除图像噪声的工作称之为图像滤波或平滑。数字图像噪声去除涉及光学***、微电子技术、计算机科学、数学分析等领域，是一门综合性很强的边缘科学。在理论上噪声可以定义为“不可预测，只能用概率统计方法来认识的随机误差”，因此图像噪声看是多维随机过程，描述噪声的方法可以借用随机过程的描述，即用其概率分布函数和概率密度分布函数。但在很多情况下，这样描述方法是很复杂，甚至不可能的。实际应用中，通常使用其数值特征，即均值方差、相关函数等。因为这些数值特征都可以从某些方面反映出噪声的特征。噪声对图像信号幅度和相位的影响十分复杂，有些噪声和图像信号相互独立不相关，有些是相关的，噪声本身之间也可能相关。遥感图像在成像、获取和传输的过程中,往往会受到很多噪声的影响,其中最为常见的噪声为高斯噪声、云噪声和雾噪声等。遥感图像中,噪声具有随机性与不确定性,这些噪声的存在,会降低图像的质量，抑制有用的信息，影响信息的精度，甚至会导致错误的判断。图像的边缘纹理等细节的模糊,给遥感图像的识别和分析带来一定困难。为获得清晰的、高质量的遥感图像必须进行降噪预处理。遥感图像去噪的目标在于在保护图像细节信息的前提下,最大限度地去除噪声,提高数据的可读性与有效性。遥感图像去噪方法其实质是在去除噪声和维持图像细节上保持平衡。高光谱图像是由二维空间信息和一维光谱信息组成的三维数据，被认为是一个3D张量，对应三阶张量的表现形式。通常原始高光谱图像的尺寸较大，不宜直接进行张量运算。高光谱图像质量受多种复杂因素影响，其获取和传输过程中引入了大量噪声，如热噪声、光子噪声等。这些噪声不仅影响图像视觉效果还影响其后续处理和应用(如高光谱图像特征提取、分类、混合像元分解、目标检测等)。
近年来,能够获取的不同空间分辨率的高光谱卫星遥感数据愈发丰富。然而,高光谱图像在获取和传播过程中难以避免会受到噪声污染。为了获得高质量遥感图像,需对退化图像复原以提升遥感图像质量。现有大多数图像复原算法都是建立在模糊图像点扩散函数和图像噪声已知的前提下进行的。但实际成像过程中,遥感图像退化因素复杂,无法精确知道图像模糊核函数和遥感图像噪声分布。
由于高光谱图像包含数百个波段的光谱，其具有丰富的光谱信息，除此之外，高光谱图像是由两个空间维和一个光谱维构成，因此，其具备良好的空间结构信息。也就是说，一个像素的空间邻域像素有很大概率是同源的，代表同一种地物。Y.Chen等人假设测试像素与其邻域像素具有相同的稀疏结构，提出了联合稀疏表示分类算法。利用邻域像素的相似性在某种程度上也能够抑制噪声的影响。Y.J.Deng等人提出了张量局部保持投影的方法，利用邻域像素块来辅助提升中心像素低维特征的可分性，取得了较好的效果。G.Camps-Valls等人提出了组合核支持向量机分类算法，通过构建光谱核和空-谱核来充分利用高光谱数据的空间特性，将高光谱数据看作一个三阶张量。然而，上述方法均是有针对性的从某种单一任务(分类、特征提取等)的角度考虑空间邻域信息，从而降低噪声对性能的影响。从广泛应用的角度来看，这种方式具有一定的局限性；从对噪声的处理上来看，这种方式仅仅是抑制噪声的负面影响，并没有消除噪声。另外一种处理高光谱数据中噪声的方法是去噪算法，高光谱图像去噪技术也是遥感图像处理技术领域的研究热点。H.Y.Zhang等人提出了基于低秩矩阵恢复的高光谱图像复原算法，能够去除高斯噪声、脉冲噪声、条纹噪声等多种噪声。为了利用高光谱数据的空间结构信息，随后提出了变分正则化的低秩矩阵分解算法。同样，以光谱信息和空间信息的有效利用为目的，J.Z.Xue等人提出空谱联合的低秩正则化去噪算法，同时考虑了光谱域的低秩特性和空间域的非局部低秩特性。在图像去噪领域，还存在一种最为经典的方法—鲁棒性主成分分析，通过将原始数据分解为低秩分量和稀疏噪声分量，构建鲁棒性主成分分析模型，将模型求解后的低秩分量作为无噪声数据。S.H.Mei等人以鲁棒性主成分分析模型为基础，设计了光谱域去噪方法和空间域去噪方法，从分类性能的角度评估，去噪效果较好。然而，一方面，鲁棒性主成分分析模型仅能处理二维数据，三维的高光谱数据只能转换为矩阵的形式进行处理，丢失了原始数据的子空间结构；另一方面，鲁棒性主成分分析模型将原始数据分解为低秩分量和稀疏噪声分量，忽略了每个子空间固有的判别信息，降低了去噪后数据的可分性。图像中的有用信息部分和噪声往往在频带上存在重迭，这是造成基于有用信息和噪声频率特性差别的传统去噪方法存在缺陷的根本原因。
由于高光谱图像各波段噪声强度分布的复杂性，传统的特征降维类去噪方法并不能满足高光谱遥感数据的去噪要求。传统的高光谱遥感图像去噪方法基本可以分为两大类：1)在空间域或频率域进行噪声去除。如Yuan等提出的光谱-空间域自适应全变分去噪模型，Rasti等提出了基于稀疏降秩回归小波的高光谱图像去噪方法,此类方法充分考虑了高光谱数据的二维空间信息，但并没能充分利用高光谱数据各波段间的光谱信息。2)用特征降维方法进行去噪。如Stephan等利用主成分分析(PCA)进行高光谱图像的消噪处理,Green等提出最小噪声分离(Minimum Noise Fraction,MNF)理论,该类方法充分考虑了高光谱图像各波段间的相关性，但由于缺乏对单张图像本身细节(即二维空间信息)的考虑，导致去噪效果也不甚理想。上述方法均在PCA的基础上结合了多尺度几何变换技术应用于高光谱图像去噪，但应用多尺度几何变换对图像进行稀疏表示时得到的字典具有特定的解析公式，能够隐含地快速实现；但其对应的稀疏字典结构是固定的，没有充分考虑图像本身的属性，这也导致其在对图像内容特殊的高光谱遥感图像去噪时，自适应性受到局限。
稀疏表示目前广泛应用于图像去噪领域，而如何构造合适的字典是进行稀疏表示的重要环节。字典构造方法主要可以分为两大类：一类是固定字典构造方法，常用的如离散余弦变换(DCT)字典、小波字典等。此类方法的缺点在于，使用同一种稀疏基对细节特征差别较大的图像进行稀疏表示时效果差距较大。第二类方法是设计自适应性的冗余字典。此类方法弥补了第一类方法的不足，能够根据图像自身的特点，自适应性地选择稀疏基，从而实现对各类图像的精确稀疏逼近。其中，Aharon等提出的K-奇异值分解(K-Singular ValueDecomposition,K-SVD)算法是目前最具代表性、应用最广泛的自适应字典学习算法。以K-SVD算法为基础，Elad等又提出了利用K-SVD算法进行图像去噪的去噪模型。但该算法的使用有个前提，即图像的噪声强度是已知的。传统的PCA去噪方法是选取几个较大的主成分分量用于PCA变换去噪，由于经过PCA变换产生的含噪图像组的图像特征比较复杂，固定的稀疏基难以准确地对图像进行稀疏表示，会导致图像损失较多的细节信息，造成去噪过程中图像失真较为严重，以致融合后的图像细节信息丢失。在实际应用中，需要先对未知图像的噪声强度进行估计，才能继续进行去噪操作。相比PCA方法，PCA-Bish和PCA-Contourlet方法去噪后图像的边缘信息得到了较好保留，但图像的细节信息，尤其是对比度方面损失较为严重，这主要是由于Bish小波阈值收缩方法和Contourlet方法去噪时使用固定的稀疏基，导致融合后的图像细节信息丢失。
目前针对高光谱图像去噪问题提出了多种方法。这些方法主要分为3类：逐波段处理、联合光谱-空间域变换方法和基于张量分解方法。
逐波段处理方法将高光谱图像每个波段的二维图像分别进行处理，但这类方法忽略了光谱相关性。逐波段的BwK-SVD方法，没有考虑每个波段的噪声强度和不同像素空间信息的差异，去噪后的图像过平滑，且在平滑区域仍有部分噪声。
联合光谱-空间域变换方法，对高光谱图像的空间域和光谱域进行特定变换来进行去噪。主要包括混合空间-光谱小波收缩去噪(HSSNR)方法、使用小波收缩和PCA方法、光谱-空间自适应总变分(SSAHTV)模型进行高光谱图像去噪方法、非局部光谱-空间结构稀疏表示的高光谱图像去噪方法等对高光谱图像进行降维。该方法虽然能有效对高光谱图像去噪，但运算时间较长。
基于张量分解方法，通过对高光谱数据进行TUCKER3分解后，执行维纳滤波实现去噪。比如低秩张量近似(LRTA)方法，该方法的特点是通过Tucker分解获得高光谱图像低秩近似进行去噪，但在重构过程会产生伪影。基于张量方法虽然将整个高光谱图像当作一个张量考虑了高光谱图像各波段的相关性，但忽略了高光谱图像的空间信息的非局部相似性，没有反映图像结构信息，导致部分细节信息丢失且去噪结果过平滑。
发明内容
本发明的目的是针对高光谱图像各波段图像噪声分布复杂，传统去噪方法难以达到理想效果而存在的不足之处，提供一种能更有效除去高光谱图像中噪声的去噪方法，尤其是适用于高光谱遥感图像的稀疏判别张量鲁棒性PCA去噪方法。
本发明的上述目的可以通过以下措施来达到，一种稀疏判别张量鲁棒性PCA的遥感图像去噪方法，其特征在于，包括如下步骤：对输入高光谱遥感图像进行图像分割，按固定尺寸，将高光谱遥感图像分割成不重叠的张量数据块；分别对每个张量数据块进行稀疏判别张量鲁棒性主成分分析：首先将张量数据块分解为低秩分量、稀疏判别分量和稀疏噪声分量，对切割的张量数据块逐一进行成分分解，再根据张量数据块组成的成分矩阵，构建稀疏判别的张量鲁棒性主成分分析模型，迭代求解该模型得到去噪后的张量数据块，重构组合图像分块，对每个张量数据块分别进行稀疏判别张量鲁棒性主成分分析去噪，获得去噪图像；对重构去噪后的无噪声数据，按比例划分训练集和测试集，最后将训练集和测试集输入到分类器，通过分类器输出所有测试样本的类别标记，实现对高光谱遥感图像去噪效果的评估。
本发明的有效增益在于：
(1)本发明利用高光谱数据良好的空间结构信息，以张量表示的形式对每一个数据块进行空间结构信息的挖掘，固定尺寸，将高光谱遥感图像分割成不重叠的张量数据块；分别对每个张量数据块进行稀疏判别张量鲁棒性主成分分析：能够更好地表征数据的本征结构，对数据中噪声的去除提供有力支撑。
(2)本发明将原始高光谱数据分解为低秩分量、稀疏判别分量和稀疏噪声分量，有效地保留了各个子空间的判别信息。与鲁棒性主成分分析模型的低秩分量和稀疏噪声分量的二分量分解方式相比，本发明采用的低秩分量、稀疏判别分量、稀疏噪声分量的三分量数据分解方式，尽可能地保留原始高光谱数据中的有价值信息，能够更好地将有用信息与噪声分离，显著提高了去噪后高光谱图像的数据质量。
(3)本发明通过张量表示，利用空间邻域信息挖掘数据的本征结构，进而利用三分量的数据分解方式，保留了原始数据中低秩分量和稀疏判别分量组成的有用信息，去除稀疏噪声分量的干扰信息，增强了数据的判别能力。本发明在公开高光谱数据集上的实验结果表明，所提出的方法在性能上优于其它的图像去噪算法，比如经典的鲁棒性主成分分析方法、空谱联合表示的方法等。
由于本发明的核心是数据成分的多重分解和空间邻域信息的有效利用，因此，只要涉及图像噪声去除和图像复原的，本发明都是有效的。本发明利用高光谱图像是由代表不同地物类别的像素和噪声组成的，相同地物类别的像素表现出低秩特性，噪声表现出稀疏特性，而不同类别像素又具有各自的判别信息，运用判别信息表现出的稀疏特性，将数据分解为低秩分量、稀疏判别分量、稀疏噪声分量，构建稀疏判别张量鲁棒性主成分分析模型，能够更准确地保留数据有价值的成分，去除噪声等干扰成分。稀疏判别张量鲁棒性PCA模型以张量表示的形式来充分挖掘数据的本征结构信息，通过组合图像分块进行重构，获得去噪图像；对重构去噪后的无噪声数据按比例划分训练集和测试集，最后将训练集和测试集输入到分类器，通过分类器输出所有测试样本的类别标记，实现对高光谱遥感图像去噪效果的评估。结果表明上述方法取得了积极效果，从分类性能的角度评价，去噪效果优于经典的鲁棒性主成分分析及其扩展的方法。
本发明所建立的稀疏判别张量鲁棒性PCA模型，适用于高光谱遥感图像，并且能够推广应用于一般的图像去噪。
附图说明
图1是本发明基于稀疏判别张量鲁棒性主成分分析去噪方法的流程图。
图2是图1稀疏判别张量鲁棒性主成分分析模型的求解流程图。
为使本发明的目的、技术方案和优点更加清楚明白，以下结合具体实施例，并参照附图，对本发明进一步详细说明。
具体实施方式
参阅图1。根据本发明，对输入高光谱遥感图像进行图像分割，按固定尺寸，将高光谱遥感图像分割成不重叠的张量数据块；分别对每个张量数据块进行稀疏判别张量鲁棒性主成分分析：首先将张量数据块分解为低秩分量、稀疏判别分量和稀疏噪声分量，对切割的张量数据块逐一进行成分分解，再根据张量数据块组成的成分矩阵，构建稀疏判别的张量鲁棒性主成分分析模型，迭代求解该模型得到去噪后的张量数据块，重构组合图像分块，对每个张量数据块分别进行稀疏判别张量鲁棒性主成分分析去噪，获得去噪图像；对重构去噪后的无噪声数据，按比例划分训练集和测试集，最后将训练集和测试集输入到分类器，通过分类器输出所有测试样本的类别标记，实现对高光谱遥感图像去噪效果的评估。
具体包括以下步骤：
步骤1，由于原始高光谱图像是一个较大尺寸的三维数据矩阵，直接对原始数据进行处理会带来较高的复杂度。需要对输入高光谱遥感图像，进行图像分割；本实施例，按一定的尺度r×r将原始数据H分割成不重叠的子张量数据块X。张量数据块X包含相同类别的共有信息L可以由低秩特性刻画，称为低秩分量；每个子空间特有的判别信息S可以由稀疏特性刻画，称为稀疏判别分量；数据中的噪声分量E一般表现出稀疏特性，称为稀疏噪声分量。
步骤2，对每个张量数据块，分别进行稀疏判别张量鲁棒性主成分分析去噪。
所述步骤2如图2所示，进一步包括以下步骤：在稀疏判别张量鲁棒性主成分分析模型求解中，根据数据中的相同类别像素的共有信息L、各个类别像素中独有的判别信息S、噪声干扰E将张量数据块X按如下形式分解：
X＝L+S+E (1)
然后根据张量数据块成分分解模式X＝L+S+E求解图像块对应的各个分量，对切割的张量数据块X逐一进行成分分解，再根据组成的成分矩阵，构建如下稀疏判别张量鲁棒性主成分分析模型：
式中，||·||*表示张量核范数，||·||1表示l1范数，β和λ表示正则化参数，s.t.表示subject to后接目标函数的约束条件。
在可选的实施例中，根据代数理论，在求解稀疏判别张量鲁棒性主成分分析模型(2)中引入拉格朗日乘子D，构建稀疏判别张量鲁棒性主成分分析模型(2)的拉格朗日函数，然后利用交替方向法求解拉格朗日函数l，得到低秩分量L、稀疏判别分量S和稀疏噪声分量E，更新拉格朗日乘子和惩罚系数μ，写出如下稀疏判别张量鲁棒性主成分分析模型(2)的拉格朗日函数：
其中,l表示拉格朗日函数，μ表示惩罚因子，F表示Frobenius范数。
采用交替迭代的策略优化稀疏判别张量鲁棒性主成分分析模型(2)，求解拉格朗日函数(3)；本实施例采用交替方向法(Alternating Direction Method,ADM)优化求解拉格朗日函数(3)，更新拉格朗日乘子D和惩罚因子μ，核心思想是每次保持其它变量不变，只更新低秩分量L、稀疏判别分量S和稀疏噪声分量E中的一个变量，交替迭代进行求解。各个变量的更新表达式如下：
拉格朗日乘子的更新表达式为：D (l)＝D (l-1)+μ(l-1)(X-L (l)-S (l)-E (l))，
惩罚因子μ的更新表达式如下：μ(l)＝min(ρμ(l-1),106)，
其中，ρ＝1.1，
表示阈值为
的软阈值运算。
交替迭代求解后，判断是否满足前后两次迭代变量的改变量充分小或者达到最大迭代次数的迭代终止条件，否，则返回交替方向法优化求解；是，则终止迭代，输出去噪后的无噪张量数据块(L+S)，然后判断是否完成所有张量数据块的去噪运算：否，则继续执行张量数据块成分分解模式X＝L+S+E(1)与交替迭代求解稀疏判别张量鲁棒性主成分分析模型(2)；是，则终止循环，输出与原始数据尺度对应的无噪数据
步骤3，以分类性能评估稀疏判别张量鲁棒性主成分分析模型(2)去噪效果，将重构的无噪数据
划分为训练集和测试集。
训练集也称为字典，是由一定数量的带类别标签的数据构成的数据矩阵。设训练集为
由N个维度为D的样本组成，其中第i个样本为
是一个一维向量，同样地，测试集为
第i个测试样本为
并且测试集Y包含M个未知类别标签的样本，其中，
表示实数空间。
步骤4，通过分类器输出所有测试样本的类别标记。分类器可以是最近邻分类器、支持向量机分类器、组合核支持向量机分类器等。将训练集和测试集输入到指定的分类器，输出测试样本的类别标签，通过比对测试样本的标签得到分类器的分类精度，并以此来评估所提出去噪模型的去噪效果。
以上所述的具体实施例，对本发明的目的、技术方案和有益效果进行了进一步详细说明，所应理解的是，以上所述仅为本发明的具体实施例而已，并不用于限制本发明，凡在本发明的精神和原则之内，所做的任何修改、等同替换、改进等，均应包含在本发明的保护范围之内。
Claims (10)
1.一种稀疏判别张量鲁棒性PCA的遥感图像去噪方法，其特征在于，包括如下步骤：对输入高光谱遥感图像进行图像分割，按固定尺寸，将高光谱遥感图像分割成不重叠的张量数据块；分别对每个张量数据块进行稀疏判别张量鲁棒性主成分分析：首先将张量数据块分解为低秩分量、稀疏判别分量和稀疏噪声分量，对切割的张量数据块逐一进行成分分解，再根据张量数据块组成的成分矩阵，构建稀疏判别的张量鲁棒性主成分分析模型，迭代求解该模型得到去噪后的张量数据块，重构组合图像分块，对每个张量数据块分别进行稀疏判别张量鲁棒性主成分分析去噪，获得去噪图像；对重构去噪后的无噪声数据，按比例划分训练集和测试集，最后将训练集和测试集输入到分类器，通过分类器输出所有测试样本的类别标记，实现对高光谱遥感图像去噪效果的评估。
2.如权利要求1所述的稀疏判别张量鲁棒性PCA的遥感图像去噪方法，其特征在于：张量数据块X包含相同类别的共有信息L由低秩特性刻画为低秩分量，每个子空间特有的判别信息S由稀疏特性刻画为稀疏判别分量，数据中的噪声分量E表现出稀疏特性，称为稀疏噪声分量。
3.如权利要求1或2所述的稀疏判别张量鲁棒性PCA的遥感图像去噪方法，其特征在于：在稀疏判别张量鲁棒性主成分分析模型求解中，根据数据中的相同类别像素的共有信息L、各个类别像素中独有的判别信息S、噪声干扰E，将张量数据块X按如下形式分解：
X＝L+S+E (1)
然后根据张量数据块成分分解模式X＝L+S+E求解图像块对应的各个分量，对切割的张量数据块X逐一进行成分分解，再根据张量数据块组成的成分矩阵，构建如下稀疏判别张量鲁棒性主成分分析模型：
其中，||·||*表示张量核范数，||·||1表示l1范数，β和λ表示正则化参数，s.t.表示subject to后接目标函数的约束条件。
4.如权利要求3所述的稀疏判别张量鲁棒性PCA的遥感图像去噪方法，其特征在于：根据代数理论，在求解稀疏判别张量鲁棒性主成分分析模型(2)中引入拉格朗日乘子D，构建稀疏判别张量鲁棒性主成分分析模型(2)的拉格朗日函数，然后利用交替方向法求解拉格朗日函数
得到低秩分量L、稀疏判别分量S和稀疏噪声分量E，更新拉格朗日乘子和惩罚系数μ，写出如下稀疏判别张量鲁棒性主成分分析模型(2)的拉格朗日函数：
其中,F表示Frobenius范数。
5.如权利要求4所述的稀疏判别张量鲁棒性PCA的遥感图像去噪方法，其特征在于：采用交替迭代的策略优化稀疏判别张量鲁棒性主成分分析模型求解拉格朗日函数；然后采用交替方向法ADM优化求解拉格朗日函数，更新拉格朗日乘子D和惩罚因子μ，每次保持其它变量不变，只更新低秩分量L、稀疏判别分量S和稀疏噪声分量E中的一个变量，交替迭代进行求解；拉格朗日乘子D的更新表达式为：D (l)＝D (l-1)+μ(l-1)(X-L (l)-S (l)-E (l))，
惩罚因子μ的更新表达式如下：μ(l)＝min(ρμ(l-1),106)，
其中，ρ＝1.1，
表示阈值为
的软阈值运算。
6.如权利要求5所述的稀疏判别张量鲁棒性PCA的遥感图像去噪方法，其特征在于：交替迭代求解后，判断是否满足前后两次迭代变量的改变量充分小或者达到最大迭代次数的迭代终止条件，否，则返回交替方向法优化求解；是，则终止迭代，输出去噪后的无噪张量数据块(L+S)，然后判断是否完成所有张量数据块的去噪运算：否，则继续执行张量数据块成分分解模式X＝L+S+E与交替迭代求解稀疏判别张量鲁棒性主成分分析模型；是，则终止循环，输出与原始数据尺度对应的无噪数据。
7.如权利要求1所述的稀疏判别张量鲁棒性PCA的遥感图像去噪方法，其特征在于：训练集为字典，是由一定数量的带类别标签的数据构成的数据矩阵。
8.如权利要求1所述的稀疏判别张量鲁棒性PCA的遥感图像去噪方法，其特征在于：训练集
由N个维度为D的样本组成，其中，第i个样本为
是一个一维向量，同样地，测试集为
第i个测试样本为
并且测试集Y包含M个未知类别标签的样本，其中，
表示实数空间。
9.如权 利要求1所述的稀疏判别张量鲁棒性PCA的遥感图像去噪方法，其特征在于：分类器是最近邻分类器、支持向量机分类器、组合核支持向量机分类器的一种。
10.如权利要求1所述的稀疏判别张量鲁棒性PCA的遥感图像去噪方法，其特征在于：将训练集和测试集输入到指定的分类器，通过分类器输出测试样本的类别标签，通过比对测试样本的标签得到分类器的分类精度，并以此来评估所提出去噪模型的去噪效果。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202010129712.9A CN111369457B (zh) | 2020-02-28 | 2020-02-28 | 稀疏判别张量鲁棒性pca的遥感图像去噪方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202010129712.9A CN111369457B (zh) | 2020-02-28 | 2020-02-28 | 稀疏判别张量鲁棒性pca的遥感图像去噪方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN111369457A CN111369457A (zh) | 2020-07-03 |
| CN111369457B true CN111369457B (zh) | 2022-05-17 |
Family
ID=71208315
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202010129712.9A Active CN111369457B (zh) | 2020-02-28 | 2020-02-28 | 稀疏判别张量鲁棒性pca的遥感图像去噪方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN111369457B (zh) |
Families Citing this family (9)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN112101381B (zh) * | 2020-08-30 | 2022-10-28 | 西南电子技术研究所（中国电子科技集团公司第十研究所） | 张量协作图判别分析遥感图像特征提取方法 |
| CN112232250A (zh) * | 2020-10-23 | 2021-01-15 | 镇江昭远智能科技有限公司 | 一种基于判别性低秩分解与稀疏表示的人脸识别方法 |
| CN112070181B (zh) * | 2020-11-16 | 2021-02-19 | 深圳市华汉伟业科技有限公司 | 一种基于图像流的协同检测方法及装置、存储介质 |
| CN112633202B (zh) * | 2020-12-29 | 2022-09-16 | 河南大学 | 一种基于双重去噪联合多尺度超像素降维的高光谱图像分类算法 |
| CN112767261B (zh) * | 2021-01-06 | 2024-02-13 | 温州大学 | 一种基于广义非凸张量鲁棒主成分分析模型的针对彩色图像和视频的非局部去噪框架 |
| CN112950500B (zh) * | 2021-02-25 | 2022-08-26 | 桂林电子科技大学 | 一种基于边缘检测低秩全变分模型的高光谱去噪方法 |
| CN113435243B (zh) * | 2021-05-14 | 2024-06-14 | 西安电子科技大学 | 一种高光谱真实下采样模糊核估计方法 |
| CN113408635B (zh) * | 2021-06-29 | 2022-03-08 | 哈尔滨工业大学 | 基于数字表面模型辅助的高光谱图像本征分解方法 |
| CN115131226B (zh) * | 2022-04-01 | 2024-06-14 | 重庆大学 | 一种基于小波张量低秩正则化的图像复原方法 |
Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103971364A (zh) * | 2014-04-04 | 2014-08-06 | 西南交通大学 | 基于加权Gabor小波特征和两级聚类的遥感图像变化检测方法 |
| CN105069758A (zh) * | 2015-08-21 | 2015-11-18 | 武汉大学 | 一种基于鲁棒低秩张量的高光谱图像去噪方法 |
| CN105657436A (zh) * | 2015-12-31 | 2016-06-08 | 中国人民解放军国防科学技术大学 | 一种图像处理方法 |
| CN107491591A (zh) * | 2017-07-23 | 2017-12-19 | 西南电子技术研究所（中国电子科技集团公司第十研究所） | 快速生成敏捷成像卫星对地观测需求筹划的方法 |
| CN109446473A (zh) * | 2018-07-02 | 2019-03-08 | 电子科技大学 | 基于分块的稳健张量主成分分析方法 |
| CN110400276A (zh) * | 2019-08-28 | 2019-11-01 | 上海海洋大学 | 高光谱图像去噪方法、装置 |
-
2020
- 2020-02-28 CN CN202010129712.9A patent/CN111369457B/zh active Active
Patent Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103971364A (zh) * | 2014-04-04 | 2014-08-06 | 西南交通大学 | 基于加权Gabor小波特征和两级聚类的遥感图像变化检测方法 |
| CN105069758A (zh) * | 2015-08-21 | 2015-11-18 | 武汉大学 | 一种基于鲁棒低秩张量的高光谱图像去噪方法 |
| CN105657436A (zh) * | 2015-12-31 | 2016-06-08 | 中国人民解放军国防科学技术大学 | 一种图像处理方法 |
| CN107491591A (zh) * | 2017-07-23 | 2017-12-19 | 西南电子技术研究所（中国电子科技集团公司第十研究所） | 快速生成敏捷成像卫星对地观测需求筹划的方法 |
| CN109446473A (zh) * | 2018-07-02 | 2019-03-08 | 电子科技大学 | 基于分块的稳健张量主成分分析方法 |
| CN110400276A (zh) * | 2019-08-28 | 2019-11-01 | 上海海洋大学 | 高光谱图像去噪方法、装置 |
Non-Patent Citations (1)
| Title |
|---|
| 基于低秩张量分解和统计流形的数据表示及其应用;范海燕;《中国博士学位论文全文数据库 (基础科学辑)》;20200215;A002-83 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN111369457A (zh) | 2020-07-03 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN111369457B (zh) | 稀疏判别张量鲁棒性pca的遥感图像去噪方法 | |
| Borsoi et al. | A fast multiscale spatial regularization for sparse hyperspectral unmixing | |
| Zhong et al. | Multiple-spectral-band CRFs for denoising junk bands of hyperspectral imagery | |
| Lebrun et al. | Secrets of image denoising cuisine | |
| Saxena et al. | Noises and image denoising techniques: A brief survey | |
| Zhang et al. | Joint image denoising using adaptive principal component analysis and self-similarity | |
| CN110400276B (zh) | 高光谱图像去噪方法、装置 | |
| CN104123705B (zh) | 一种超分辨率重建图像质量Contourlet域评价方法 | |
| CN114677300A (zh) | 一种基于双阶段学习框架的高光谱图像深度降噪的方法及*** | |
| Letexier et al. | Nonorthogonal tensor matricization for hyperspectral image filtering | |
| CN104657951A (zh) | 图像乘性噪声移除方法 | |
| CN112767271B (zh) | 一种基于三维变分网络的高光谱图像深度降噪的方法 | |
| Shahdoosti et al. | Combined ripplet and total variation image denoising methods using twin support vector machines | |
| Malladi et al. | Image denoising using superpixel-based PCA | |
| Naveed et al. | Multiscale image denoising using goodness-of-fit test based on EDF statistics | |
| Moudden et al. | Hyperspectral BSS using GMCA with spatio-spectral sparsity constraints | |
| Lu et al. | SAR image despeckling via structural sparse representation | |
| Golla et al. | A novel approach of k-SVD-based algorithm for image denoising | |
| Banothu et al. | High-order total bounded variation approach for gaussian noise and blur removal | |
| Pan et al. | Nonlocal low rank regularization method for fractal image coding under salt-and-pepper noise | |
| CN112862729B (zh) | 一种基于特征曲线指导的遥感图像去噪方法 | |
| CN113160069B (zh) | 一种基于图信号的高光谱图像去噪方法 | |
| Song et al. | Unsupervised denoising for satellite imagery using wavelet subband cyclegan | |
| Ferzo et al. | Image Denoising Techniques Using Unsupervised Machine Learning and Deep Learning Algorithms: A Review | |
| Bhargava et al. | An Effective Method for Image Denoising Using Non-local Means and Statistics based Guided Filter in Nonsubsampled Contourlet Domain. |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |