CN112950500A - Hyperspectral denoising method based on edge detection low-rank total variation model - Google Patents

Hyperspectral denoising method based on edge detection low-rank total variation model Download PDF

Info

Publication number
CN112950500A
CN112950500A CN202110213967.8A CN202110213967A CN112950500A CN 112950500 A CN112950500 A CN 112950500A CN 202110213967 A CN202110213967 A CN 202110213967A CN 112950500 A CN112950500 A CN 112950500A
Authority
CN
China
Prior art keywords
edge detection
total variation
hyperspectral
subproblem
rank
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.)
Granted
Application number
CN202110213967.8A
Other languages
Chinese (zh)
Other versions
CN112950500B (en
Inventor
蒋俊正
蔡万源
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Guilin University of Electronic Technology
Original Assignee
Guilin University of Electronic Technology
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Guilin University of Electronic Technology filed Critical Guilin University of Electronic Technology
Priority to CN202110213967.8A priority Critical patent/CN112950500B/en
Publication of CN112950500A publication Critical patent/CN112950500A/en
Application granted granted Critical
Publication of CN112950500B publication Critical patent/CN112950500B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

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
    • G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
    • G06F17/10—Complex mathematical operations
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F18/00—Pattern recognition
    • G06F18/20—Analysing
    • G06F18/22—Matching criteria, e.g. proximity measures
    • 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/13—Edge detection
    • 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
    • Y02A90/00—Technologies having an indirect contribution to adaptation to climate change
    • Y02A90/30—Assessment of water resources

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Data Mining & Analysis (AREA)
  • Mathematical Physics (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • General Engineering & Computer Science (AREA)
  • Bioinformatics & Computational Biology (AREA)
  • Mathematical Analysis (AREA)
  • Evolutionary Biology (AREA)
  • Bioinformatics & Cheminformatics (AREA)
  • Artificial Intelligence (AREA)
  • Algebra (AREA)
  • Computational Mathematics (AREA)
  • Evolutionary Computation (AREA)
  • Mathematical Optimization (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Pure & Applied Mathematics (AREA)
  • Databases & Information Systems (AREA)
  • Software Systems (AREA)
  • Image Analysis (AREA)

Abstract

The invention discloses a hyperspectral denoising method based on an edge detection low-rank total variation model, which comprises the steps of constructing an input signal model, and constructing and optimizing the edge detection low-rank total variation model based on the input signal model; solving a first sub-problem divided in the edge detection low-rank fully-variable model by using a singular value shrinkage method; dividing the obtained second subproblems based on the number of wave bands of the hyperspectral image, and solving all the second subproblems by using an iterative gradient-based fast edge detection four-neighborhood total variation algorithm; solving the obtained third subproblem by using a soft threshold shrinkage operator, and iterating the obtained results of all the subproblems; and comparing the current iteration result with a set iteration termination condition until the iteration termination condition is met, calculating a corresponding peak signal-to-noise ratio and a corresponding structural similarity value, and improving the denoising effect of the spectral image through experimental verification analysis.

Description

Hyperspectral denoising method based on edge detection low-rank total variation model
Technical Field
The invention relates to the technical field of hyperspectral image data processing, in particular to a hyperspectral denoising method based on an edge detection low-rank total variation model.
Background
The hyperspectral imaging (HSI) technology can acquire two-dimensional images in a wide electromagnetic spectrum range and a higher spectral resolution, and thus is widely applied to the fields of archaeology and art protection, vegetation and water resource control, food quality and safety control, forensic medicine, surgery and diagnosis, crime scene detection, biomedicine, military and the like.
However, due to its unique physical design, HSI is inevitably contaminated by various noises during the acquisition process, and the common types of the contaminated noises include gaussian noise, impulse noise, dead pixel, streak noise, and the like. In recent years, many scientists have proposed various HSI denoising algorithms. Among them, the method for restoring hyperspectral images based on low-rank matrix restoration (LRMR), the method for restoring hyperspectral images based on total variation regularization low-rank matrix decomposition (LRTV), and the like are good in effect. The LRMR method utilizes the low-rank characteristic of the hyperspectral image to represent the main information of the hyperspectral image by using a low-rank matrix, so that a large amount of redundant information in the hyperspectral image is removed. Most of the noise is also contained in the redundant information, so that the denoising effect can be achieved in the process of restoring the hyperspectral image by using the low-rank matrix. However, the denoising method based on low rank only studies the correlation between spectral bands, ignores the spatial correlation of local neighborhood pixels, and thus cannot achieve the optimal denoising effect. The LRTV method describes the spatial correlation and spatial smoothness of local neighborhood pixels on the basis of low-rank characteristics, the denoising effect is better than that of LRMR, but the method has less utilization of the spatial correlation information of the neighborhood pixels, and neglects the protection of the hyperspectral image edge in the smoothing process of the local neighborhood pixels. Therefore, a design method for improving the denoising effect of the spectral image is to be further introduced.
Disclosure of Invention
The invention aims to provide a hyperspectral denoising method based on an edge detection low-rank total variation model, and the hyperspectral image denoising effect is improved.
In order to achieve the purpose, the invention provides a hyperspectral denoising method based on an edge detection low-rank total variation model, which comprises the following steps:
constructing an input signal model, and constructing and optimizing an edge detection low-rank total variation model based on the input signal model;
solving a first sub-problem divided in the edge detection low-rank fully-variable model by using a singular value shrinkage method;
dividing the obtained second subproblems based on the number of wave bands of the hyperspectral image, and solving all the second subproblems by using an iterative gradient-based fast edge detection four-neighborhood total variation algorithm;
solving the obtained third subproblem by using a soft threshold shrinkage operator, and iterating the obtained results of all the subproblems;
and comparing the current iteration result with a set iteration termination condition until the iteration termination condition is met, and calculating a corresponding peak signal-to-noise ratio and a corresponding structural similarity value.
The method comprises the following steps of constructing an input signal model, constructing and optimizing an edge detection low-rank total variation model based on the input signal model, and comprises the following steps:
sequentially adding the acquired original image with randomly generated sparse noise and Gaussian noise to obtain an input signal model;
and constructing an edge detection low-rank total variation model based on the input signal model, optimizing the edge detection low-rank total variation model by using an augmented Lagrange function method, and dividing main sub-problems, wherein the divided sub-problems comprise a first sub-problem, a second sub-problem and a third sub-problem.
Dividing the obtained second subproblems based on the number of the wave bands of the hyperspectral image, and solving all the second subproblems by using an iterative gradient-based fast edge detection four-neighborhood total variation algorithm, wherein the method comprises the following steps of:
dividing the obtained second sub-problem into a plurality of band sub-problems according to the number of the bands of the hyperspectral image;
rewriting and iterating each wave band subproblem, and assigning values to the detected pixel points by using an edge detection operator in the iteration process;
and carrying out iterative calculation on the difference absolute values of the four adjacent domains of all the pixel points which meet the conditions after assignment until the set iterative conditions are met, and obtaining the solution of the corresponding second subproblem.
After assignment, iterative computation is performed on the four-neighbor domain difference absolute values of all pixel points meeting the condition until a set iteration condition is met, and a solution of a corresponding second subproblem is obtained, wherein the iterative computation comprises the following steps:
turning over the value of the assigned pixel point to obtain a detection value;
and adding the difference absolute values of the four adjacent domains of each pixel point, calculating the difference between the result value of the constraint condition of the current calculation result and the result value of the constraint condition of the last calculation result, dividing the obtained absolute value of the difference by the result value of the constraint condition of the current calculation result, terminating iteration if the obtained calculated value is smaller than the set iteration condition, and overlapping all frequency band result values to obtain the solution of the second subproblem.
Comparing the current iteration result with a set iteration termination condition until the iteration termination condition is met, and calculating a corresponding peak signal-to-noise ratio and a corresponding structural similarity value, wherein the method comprises the following steps:
if the current calculation result does not meet the set iteration termination condition, solving the first subproblem, the second subproblem, the third subproblem and all other subproblems again until the set iteration termination condition is met, wherein the other subproblems are the problems except the first subproblem, the second subproblem and the third subproblem;
and after the set iteration termination condition is met, comparing the obtained hyperspectral denoised image with the original image to obtain a corresponding peak signal-to-noise ratio and a corresponding structural similarity value.
The invention discloses a hyperspectral denoising method based on an edge detection low-rank total variation model, which comprises the steps of constructing an input signal model, and constructing and optimizing the edge detection low-rank total variation model based on the input signal model; solving a first sub-problem divided in the edge detection low-rank fully-variable model by using a singular value shrinkage method; dividing the obtained second subproblems based on the number of wave bands of the hyperspectral image, and solving all the second subproblems by using an iterative gradient-based fast edge detection four-neighborhood total variation algorithm; solving the obtained third subproblem by using a soft threshold shrinkage operator, and iterating the obtained results of all the subproblems; comparing the current iteration result with a set iteration termination condition until the iteration termination condition is met, calculating a corresponding peak signal-to-noise ratio and a corresponding structural similarity value, and denoising the hyperspectral image by utilizing an edge detection four-neighborhood total variation algorithm, thereby not only enhancing the smoothness relation between neighborhoods, but also protecting the edges in the hyperspectral image and avoiding the edges from being smoothed to influence the denoising effect. Through experimental verification and analysis, the method has a better denoising effect.
Drawings
In order to more clearly illustrate the embodiments of the present invention or the technical solutions in the prior art, the drawings used in the description of the embodiments or the prior art will be briefly described below, it is obvious that the drawings in the following description are only some embodiments of the present invention, and for those skilled in the art, other drawings can be obtained according to the drawings without creative efforts.
FIG. 1 is a schematic step diagram of a hyperspectral denoising method based on an edge detection low-rank total variation model provided by the invention.
Fig. 2 is a clean artwork provided by the present invention.
Fig. 3 is a noisy image provided by the present invention.
FIG. 4 is an LRMR denoised image provided by the present invention.
FIG. 5 is a denoised LRTV image according to the present invention.
FIG. 6 is an EDTV denoised image according to the present invention.
Fig. 7 is a comparison diagram of PSNR of each band provided by the present invention.
Fig. 8 is a comparison diagram of SSIM of each band provided by the present invention.
Detailed Description
Reference will now be made in detail to embodiments of the present invention, examples of which are illustrated in the accompanying drawings, wherein like or similar reference numerals refer to the same or similar elements or elements having the same or similar function throughout. The embodiments described below with reference to the drawings are illustrative and intended to be illustrative of the invention and are not to be construed as limiting the invention.
In the description of the present invention, "a plurality" means two or more unless specifically defined otherwise.
Referring to fig. 1, the present invention provides a hyperspectral denoising method based on an edge detection low-rank total variation model, including the following steps:
s101, constructing an input signal model, and constructing and optimizing an edge detection low-rank total variation model based on the input signal model.
Specifically, an input signal model Y is constructed, wherein Y represents an input noise signal; x represents a clean original image; s represents sparse noise and is used for depicting pulse noise, dead pixels, stripe noise and the like; n represents gaussian noise.
According to an input signal model, establishing an edge detection low-rank total variation model:
Figure BDA0002952439790000041
Figure BDA0002952439790000042
wherein min represents the value of X and S when the following expression reaches the minimum value; i | · | purple wind*The kernel norm is the sum of matrix singular values, is used for convex approximation rank constraint and is used for depicting the low rank characteristic of the hyperspectral image; | X | non-conducting phosphorEDTVRepresenting the piecewise smoothness of the hyperspectral image; | S | non-woven phosphor1Representing sparse noise; s.t. objectto, with constraints indicated later; i | · | purple windFIs Frobenius norm, which is the square sum of the absolute values of the matrix elements and then the square of the square; in constraint term
Figure BDA00029524397900000510
The square of the F norm of Gaussian noise is represented, and the term is made to be as small as possible after denoising, so that the denoising effect is achieved; epsilon is a number as small as possible and is used for restricting an optimization term; rank (·) denotes the rank of the matrix; r is the set matrix rank size and is used for constraining optimization items to meet the low-rank property; both τ and λ are regular term parameters.
And (3) according to the established edge detection low-rank total variation model, carrying out optimization solution on the model:
the above problem is solved by using an augmented lagrangian method (hereinafter, referred to as ALM). The above problem was first equivalently rewritten as follows:
Figure BDA0002952439790000051
Figure BDA0002952439790000052
the equivalence of L to X before the equivalent overwrite is convenient for the subsequent use of ALM method.
By adopting the ALM method, the optimized augmented Lagrangian function is as follows:
Figure BDA0002952439790000053
s.t.rank(L)≤r
wherein the content of the first and second substances,
Figure BDA00029524397900000511
expression relating to L, X, S, Lambda1,Λ2A minimum function of the function of (a); lambda1,Λ2Optimizing a coefficient matrix;<·,·>represents the inner product; μ is a penalty factor and the initial value is set to 1 e-2.
Dividing the problem into a plurality of sub-problems, and solving one of the sub-problems in an iterative manner:
Figure BDA0002952439790000054
Figure BDA0002952439790000055
Figure BDA0002952439790000056
Figure BDA0002952439790000057
Figure BDA0002952439790000058
wherein
Figure BDA0002952439790000059
The value of the expression when the minimum value is reached is expressed; k represents the kth iteration; *(k+1)Represents the result of formula after the (k + 1) th iteration; *(k)The results of formula after the k-th iteration are shown.
Thus solving the problem to solve L, X, S the three main first to third sub-problems.
S102, solving a first sub-problem divided in the edge detection low-rank total variation model by using a singular value contraction method.
Specifically, the first sub-problem with respect to L is solved using a singular value contraction method.
Figure BDA0002952439790000061
For a given matrix W, it is decomposed using singular value decomposition, resulting in:
W＝UErV*,Er＝diag({σi}1≤i≤r)
u and V are unitary matrixes obtained by singular value decomposition of W; v*Is a conjugate transpose of V; erDetermining a diagonal matrix for the semi-positive; diag (·) denotes a diagonal matrix whose diagonal elements are constructed as · s; { sigma. }i}1≤i≤rRepresenting a set of the first r diagonal elements.
And then using a singular value contraction operator:
Figure BDA0002952439790000062
wherein D isδ(W) is formula
Figure BDA0002952439790000063
When the minimum value is reached and the L is less than r;
Dδ(W)＝UDδ(Er)V*,Dδ(Er)＝diag{max((σiδ),0), where max (x, 0) represents the maximum value for each comparison of the two compared to 0.
S103, dividing the obtained second subproblems based on the number of the wave bands of the hyperspectral image, and solving all the second subproblems by using an iterative gradient-based fast edge detection four-neighborhood total variation algorithm.
Specifically, the second sub-problem for X is solved:
Figure BDA0002952439790000071
the problem is decomposed into p band sub-problems according to the number p of bands of the hyperspectral image, i.e.
Figure BDA0002952439790000072
Wherein j is an integer from 1 to p;
Figure BDA0002952439790000073
represents the value of the jth band X after the (k + 1) th iteration;
Figure BDA0002952439790000074
after assignment, iteration calculation is carried out on the four-neighborhood difference absolute values of all the pixel points which meet the condition that the pixel point is a non-edge pixel point and the four neighborhoods of the pixel point are non-edge pixel points until set iteration conditions are met, and a solution of a corresponding second subproblem is obtained.
The problem is solved by using an iterative gradient-based fast edge detection four-neighborhood total variation algorithm.
For the sub-problem of each band described above, we can rewrite it as an equivalent problem, as follows:
Figure BDA0002952439790000075
the solution of the above problem is
Figure BDA0002952439790000076
Wherein, PCFor the orthographic projection operator, L (p, q) is a matrix pair operator, wherein
L(p,q)i,j＝pi,j+qi,j-pi-1,j-qi,j-1，i＝1,...,m,j＝1,...,n
pi,j＝xi,j-xi+1,j，i＝1,...,m-1,j＝1,...,n
qi,j＝xi,j-xi,j+1，i＝1,...,m,j＝1,...,n-1
Wherein x*,·The value of a pixel point with the abscissa as x and the ordinate as · is defined;
continuously assigning value to X in the iteration process, continuously approaching to the optimal solution, performing edge detection on X by using Sobel edge operator in each iteration process, and if the image is detectedWhen the prime point is an edge, the value is 1, otherwise, the value is 0; the 0, 1 value is turned over to obtain a detection value deltai,j。
And then, the smoothness of the hyperspectral image is depicted by adding the difference absolute values of four adjacent domains of each pixel point, when an edge is detected, the pixel point does not participate in the smoothness depiction, or when a certain point around the pixel point is an edge, the adjacent domain in the direction does not participate in the smoothness depiction, so that the smoothness depiction is obtained
Figure BDA0002952439790000081
Wherein the content of the first and second substances,
Figure BDA0002952439790000082
when iterated to satisfy this
Figure BDA0002952439790000083
And when the absolute value of the subtraction of the value and the value of the equation in the previous iteration is divided by the value of the current iteration to be less than the preset iteration condition 1e-4, stopping the iteration, and solving to obtain the calculated value X.
Calculating value X obtained from each wave bandjAnd (5) overlapping and reducing to obtain complete X.
And S104, solving the obtained third subproblem by using a soft threshold shrinkage operator, and iterating the obtained results of all the subproblems.
Specifically, a soft threshold shrinkage operator is used to solve the third sub-problem with respect to S.
Figure BDA0002952439790000084
Shrinking operators by soft threshold
Figure BDA0002952439790000085
Wherein x ∈ R, Δ > 0; then the solution to this step sub-problem is:
Figure BDA0002952439790000086
iterate to obtain Λ1
Figure BDA0002952439790000091
Iterate to obtain Λ2
Figure BDA0002952439790000092
And S105, comparing the current iteration result with a set iteration termination condition until the iteration termination condition is met, and calculating a corresponding peak signal-to-noise ratio and a corresponding structural similarity value.
Specifically, the penalty factor is set such that ρ is 1.5 and μ is obtained for each amplification stepmax1e6, iterate μ.
μ(k+1)＝min(ρμ,μmax) Wherein min (·) represents comparing values of · and taking the smaller value of the two.
Judging an iteration termination condition, if the condition is met, stopping iteration, wherein the iteration termination condition is as follows:
Figure BDA0002952439790000093
wherein | · | purple∞The matrix is an infinite norm and represents the maximum value of the sum of absolute values of row vectors of the matrix; epsilon1，ε2Is a number as small as possible, and is a constraint parameter.
And if the iteration termination condition is not met, continuously returning to solve the first subproblem again, and performing a new iteration until the termination condition is met to obtain the solution of the materialized model.
And performing peak signal-to-noise ratio (PSNR) and Structural Similarity (SSIM) parameter calculation on the resolved hyperspectral denoised image X and a clean picture without noise, and evaluating the superiority of the contrast denoising effect.
Example analysis:
the hyperspectral image input for this example is a clean hyperspectral image of indiana, i.e., X, of usa, which has 224 wavebands and has dimensions of 145X 224. Adding noise artificially, wherein the added noise comprises sparse noise (including salt and pepper noise) and Gaussian noise, and outputting Y after adding the noise, wherein tau is 0.015; λ is
Figure BDA0002952439790000094
Where M, N is the hyperspectral image dimension, this example M145, N145; r is 10; epsilon1＝ε21 e-6; the initial value of μ is 1 e-2.
Sending the three-dimensional hyperspectral image after noise addition into the model of the technical scheme, and obtaining a contrast image before and after noise removal as shown in figures 2-6, wherein figures 2-6 are gray level images which represent black depth, and the gray level value of the gray level image is smaller when the image is darker; wherein, fig. 2 shows a clean original image, fig. 3 shows a denoised image, fig. 4 shows a denoised image by LRMR method, fig. 5 shows a denoised image by LRTV method, and fig. 6 shows a denoised image according to the present technical solution; according to the comparison between fig. 2 and fig. 3, it can be observed that the experiment adds more serious noise, and the effectiveness of the technical scheme can be further embodied; according to fig. 4, it can be observed that there are still fuzzy points in the image, and the denoising effect is not good enough compared with fig. 5 and 6; observing fig. 5 and fig. 6, it is observed that the denoising effects of the two are similar globally in fig. 6, but when the detail aspect is carefully observed, the denoising effect of fig. 6 is better than that of fig. 5, that is, it can be observed that the denoising effect of the model in the technical scheme is good. Compared with a method (LRMR) for restoring a hyperspectral image based on low-rank matrix restoration proposed by Hongyan Zhang and a method (LRTV) for restoring a hyperspectral image based on holomorphic regular low-rank matrix decomposition proposed by Wei He, the PSNR and SSIM parameters of each wave band after denoising in the technical scheme are compared with the prior art scheme to obtain a graph 7-8, the graph 7 is observed, a curve (solid line) positioned at the top of the image is a PSNR parameter curve of the technical scheme, a curve (dotted line) positioned in the middle is a PSNR parameter curve of the technical scheme, a curve (dotted line) positioned at the bottom is a PSNR parameter curve of the technical scheme, and a linear type corresponding to the technical scheme is arranged at the upper right of the image; the higher the PSNR value is, the better the denoising effect is; according to the observation of FIG. 7, the denoising effect of each wave band of the PSNR is better than that of the LRMR in the technical scheme; the PSNR basically meets the condition that the denoising effect of all bands is superior to that of an LRTV technical scheme (the existing individual band is slightly lower than that of the LRTV technical scheme, but the slightly poor band is almost the same as that of the LRTV technical scheme), and compared with the LRTV technical scheme, the PSNR avoids the deterioration phenomenon of the denoising effect of the individual band; observing fig. 8, the curve (dotted line) at the top of the image is the SSIM parameter curve of the technical scheme, the curve (dotted line) at the middle is the SSIM parameter curve of the LRTV technical scheme, the curve (solid line) at the bottom is the SSIM parameter curve of the LRMR technical scheme, and there is a linear type corresponding to the technical scheme on the upper right of the image; the closer the SSIM parameter is to 1, the closer the similarity of the SSIM parameter and the original image structure is, namely the better the denoising effect is, and the conclusion consistent with the observation of the graph 7 can be observed; in conclusion, the technical scheme is superior to the prior technical scheme.
Averaging the PSNR and SSIM of each wave band in the technical scheme and the prior scheme, and recording the average as a primary result; repeating the experiment for 30 times, and obtaining results for 30 times; the PSNR and SSIM parameters obtained by averaging the results of 30 times are shown in table 1, and it can be observed that the technical solution is superior to the existing solutions.
Table 130 times comparison table of denoising effect of each scheme
Figure BDA0002952439790000101
The invention discloses a hyperspectral denoising method based on an edge detection low-rank total variation model, which comprises the steps of constructing an input signal model, and constructing and optimizing the edge detection low-rank total variation model based on the input signal model; solving a first sub-problem divided in the edge detection low-rank fully-variable model by using a singular value shrinkage method; dividing the obtained second subproblems based on the number of wave bands of the hyperspectral image, and solving all the second subproblems by using an iterative gradient-based fast edge detection four-neighborhood total variation algorithm; solving the obtained third subproblem by using a soft threshold shrinkage operator, and iterating the obtained results of all the subproblems; comparing the current iteration result with a set iteration termination condition until the iteration termination condition is met, calculating a corresponding peak signal-to-noise ratio and a corresponding structural similarity value, and denoising the hyperspectral image by utilizing an edge detection four-neighborhood total variation algorithm, thereby not only enhancing the smoothness relation between neighborhoods, but also protecting the edges in the hyperspectral image and avoiding the edges from being smoothed to influence the denoising effect. Through experimental verification and analysis, the method has a better denoising effect.
While the invention has been described with reference to a preferred embodiment, it will be understood by those skilled in the art that various changes in form and detail may be made therein without departing from the spirit and scope of the invention as defined by the appended claims.

Claims (5)

1. A hyperspectral denoising method based on an edge detection low-rank total variation model is characterized by comprising the following steps:
constructing an input signal model, and constructing and optimizing an edge detection low-rank total variation model based on the input signal model;
solving a first sub-problem divided in the edge detection low-rank fully-variable model by using a singular value shrinkage method;
dividing the obtained second subproblems based on the number of wave bands of the hyperspectral image, and solving all the second subproblems by using an iterative gradient-based fast edge detection four-neighborhood total variation algorithm;
solving the obtained third subproblem by using a soft threshold shrinkage operator, and iterating the obtained results of all the subproblems;
and comparing the current iteration result with a set iteration termination condition until the iteration termination condition is met, and calculating a corresponding peak signal-to-noise ratio and a corresponding structural similarity value.
2. The hyperspectral denoising method based on edge detection low-rank total variation model according to claim 1, wherein constructing an input signal model and constructing and optimizing an edge detection low-rank total variation model based on the input signal model comprises:
sequentially adding the acquired original image with randomly generated sparse noise and Gaussian noise to obtain an input signal model;
and constructing an edge detection low-rank total variation model based on the input signal model, optimizing the edge detection low-rank total variation model by using an augmented Lagrange function method, and dividing main sub-problems, wherein the divided sub-problems comprise a first sub-problem, a second sub-problem and a third sub-problem.
3. The hyperspectral denoising method based on the edge detection low-rank fully variant model according to claim 1, wherein the obtained second subproblems are divided based on the number of hyperspectral image bands, and all the second subproblems are solved by using an iterative gradient-based fast edge detection four-neighborhood fully variant algorithm, which comprises:
dividing the obtained second sub-problem into a plurality of band sub-problems according to the number of the bands of the hyperspectral image;
rewriting and iterating each wave band subproblem, and assigning values to the detected pixel points by using an edge detection operator in the iteration process;
and carrying out iterative calculation on the difference absolute values of the four adjacent domains of all the pixel points which meet the conditions after assignment until the set iterative conditions are met, and obtaining the solution of the corresponding second subproblem.
4. The hyperspectral denoising method based on the edge detection low-rank total variation model according to claim 3, wherein the iterative computation is performed on the four-neighbor domain difference absolute values of all the pixel points which satisfy the condition after being assigned until a set iteration condition is satisfied, and a solution of a corresponding second subproblem is obtained, and comprises:
turning over the value of the assigned pixel point to obtain a detection value;
and adding the difference absolute values of the four adjacent domains of each pixel point, calculating the difference between the result value of the constraint condition of the current calculation result and the result value of the constraint condition of the last calculation result, dividing the obtained absolute value of the difference by the result value of the constraint condition of the current calculation result, terminating iteration if the obtained calculated value is smaller than the set iteration condition, and overlapping all frequency band result values to obtain the solution of the second subproblem.
5. The hyperspectral denoising method based on the edge detection low-rank total variation model according to claim 1, wherein comparing the current iteration result with a set iteration termination condition until the iteration termination condition is satisfied, and calculating a corresponding peak signal-to-noise ratio and a corresponding structural similarity value comprises:
if the current calculation result does not meet the set iteration termination condition, solving the first subproblem, the second subproblem, the third subproblem and all other subproblems again until the set iteration termination condition is met, wherein the other subproblems are the problems except the first subproblem, the second subproblem and the third subproblem;
and after the set iteration termination condition is met, comparing the obtained hyperspectral denoised image with the original image to obtain a corresponding peak signal-to-noise ratio and a corresponding structural similarity value.
CN202110213967.8A 2021-02-25 2021-02-25 Hyperspectral denoising method based on edge detection low-rank total variation model Active CN112950500B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202110213967.8A CN112950500B (en) 2021-02-25 2021-02-25 Hyperspectral denoising method based on edge detection low-rank total variation model

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202110213967.8A CN112950500B (en) 2021-02-25 2021-02-25 Hyperspectral denoising method based on edge detection low-rank total variation model

Publications (2)

Publication Number Publication Date
CN112950500A true CN112950500A (en) 2021-06-11
CN112950500B CN112950500B (en) 2022-08-26

Family

ID=76246317

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202110213967.8A Active CN112950500B (en) 2021-02-25 2021-02-25 Hyperspectral denoising method based on edge detection low-rank total variation model

Country Status (1)

Country Link
CN (1) CN112950500B (en)

Citations (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20110170796A1 (en) * 2006-01-31 2011-07-14 Shen-En Qian Method And System For Increasing Signal-To-Noise Ratio
CN108399608A (en) * 2018-03-01 2018-08-14 桂林电子科技大学 Dimensional images denoising method based on tensor dictionary and full variation
CN109102477A (en) * 2018-08-31 2018-12-28 哈尔滨工业大学 A kind of high-spectrum remote sensing restoration methods based on the constraint of non-convex low-rank sparse
CN109658351A (en) * 2018-11-26 2019-04-19 浙江工业大学 The high spectrum image denoising method that a kind of combination L0 gradient constraint and local low-rank matrix are restored
CN111369457A (en) * 2020-02-28 2020-07-03 西南电子技术研究所（中国电子科技集团公司第十研究所） Remote sensing image denoising method for sparse discrimination tensor robustness PCA
CN111915518A (en) * 2020-07-23 2020-11-10 南京邮电大学 Hyperspectral image denoising method based on triple low-rank model
CN111951183A (en) * 2020-07-14 2020-11-17 浙江工业大学 Low-rank total-variation hyperspectral image restoration method based on near-end alternative penalty algorithm
CN111951186A (en) * 2020-07-16 2020-11-17 南京邮电大学 Hyperspectral image denoising method based on low-rank and total variation constraint
CN112069919A (en) * 2020-08-17 2020-12-11 浙江工业大学 Hyperspectral image denoising method based on non-convex low-rank matrix approximation and total variation regularization
CN112084847A (en) * 2020-07-30 2020-12-15 浙江工业大学 Hyperspectral image denoising method based on multi-channel truncation kernel norm and total variation regularization

Patent Citations (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20110170796A1 (en) * 2006-01-31 2011-07-14 Shen-En Qian Method And System For Increasing Signal-To-Noise Ratio
CN108399608A (en) * 2018-03-01 2018-08-14 桂林电子科技大学 Dimensional images denoising method based on tensor dictionary and full variation
CN109102477A (en) * 2018-08-31 2018-12-28 哈尔滨工业大学 A kind of high-spectrum remote sensing restoration methods based on the constraint of non-convex low-rank sparse
CN109658351A (en) * 2018-11-26 2019-04-19 浙江工业大学 The high spectrum image denoising method that a kind of combination L0 gradient constraint and local low-rank matrix are restored
CN111369457A (en) * 2020-02-28 2020-07-03 西南电子技术研究所（中国电子科技集团公司第十研究所） Remote sensing image denoising method for sparse discrimination tensor robustness PCA
CN111951183A (en) * 2020-07-14 2020-11-17 浙江工业大学 Low-rank total-variation hyperspectral image restoration method based on near-end alternative penalty algorithm
CN111951186A (en) * 2020-07-16 2020-11-17 南京邮电大学 Hyperspectral image denoising method based on low-rank and total variation constraint
CN111915518A (en) * 2020-07-23 2020-11-10 南京邮电大学 Hyperspectral image denoising method based on triple low-rank model
CN112084847A (en) * 2020-07-30 2020-12-15 浙江工业大学 Hyperspectral image denoising method based on multi-channel truncation kernel norm and total variation regularization
CN112069919A (en) * 2020-08-17 2020-12-11 浙江工业大学 Hyperspectral image denoising method based on non-convex low-rank matrix approximation and total variation regularization

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
欧阳宁等: "基于改进的扩散平滑和RBM的高光谱图像分类", 《电视技术》 *
蔡荣荣等: "一种基于低秩表示的高光谱图像恢复方法", 《复旦学报(自然科学版)》 *

Also Published As

Publication number Publication date
CN112950500B (en) 2022-08-26

Similar Documents

Publication Publication Date Title
Xiao et al. Removing stripe noise from infrared cloud images via deep convolutional networks
Sandić-Stanković et al. DIBR synthesized image quality assessment based on morphological wavelets
Dabov et al. Image denoising by sparse 3-D transform-domain collaborative filtering
Schulte et al. A fuzzy noise reduction method for color images
Chierchia et al. A nonlocal structure tensor-based approach for multicomponent image recovery problems
Harikiran et al. Impulse noise removal in digital images
CN106991661B (en) Non-local mean denoising method fusing KL (karhunen-Loeve) transformation and grey correlation degree
Iwasokun et al. Image enhancement methods: a review
Dar et al. Restoration by compression
Vishwakarma Denoising and inpainting of sonar images using convolutional sparse representation
CN113160069B (en) Hyperspectral image denoising method based on image signal
CN111461999B (en) SAR image speckle suppression method based on super-pixel similarity measurement
CN112950500B (en) Hyperspectral denoising method based on edge detection low-rank total variation model
Rakheja et al. Image denoising using various wavelet transforms: a survey
Banerjee et al. Bacterial foraging-fuzzy synergism based image Dehazing
Jebur et al. Image denoising techniques: An overview
Mohamadi et al. A new hybrid denoising model based on PDEs
Saleem et al. Survey on color image enhancement techniques using spatial filtering
Shivagunde et al. Single image super-resolution based on modified interpolation method using MLP and DWT
Hassanpour Image zooming using non-linear partial differential equation
CN112767263A (en) Image mixed noise removing method based on low-rank function and multi-direction total variation
AKINTOYE et al. COMPOSITE MEDIAN WIENER FILTER BASED TECHNIQUE FOR IMAGE ENHANCEMENT.
Pizurica et al. Wavelet domain denoising of single-band and multiband images adapted to the probability of the presence of features of interest
Azam et al. Remote sensing image resolution enlargement algorithm based on wavelet transformation
Mercy et al. Effective image deblurring based on improved image edge information and blur Kernel estimation

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