CN107644406B - Image denoising method based on improved orthogonal matching pursuit - Google Patents

Image denoising method based on improved orthogonal matching pursuit Download PDF

Info

Publication number
CN107644406B
CN107644406B CN201710868505.3A CN201710868505A CN107644406B CN 107644406 B CN107644406 B CN 107644406B CN 201710868505 A CN201710868505 A CN 201710868505A CN 107644406 B CN107644406 B CN 107644406B
Authority
CN
China
Prior art keywords
image
denoising
residual
orthogonal matching
algorithm
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
Application number
CN201710868505.3A
Other languages
Chinese (zh)
Other versions
CN107644406A (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.)
Nanjing University of Science and Technology
Original Assignee
Nanjing University of Science and 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 Nanjing University of Science and Technology filed Critical Nanjing University of Science and Technology
Priority to CN201710868505.3A priority Critical patent/CN107644406B/en
Publication of CN107644406A publication Critical patent/CN107644406A/en
Application granted granted Critical
Publication of CN107644406B publication Critical patent/CN107644406B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

  • Image Processing (AREA)

Abstract

The invention discloses an image denoising method based on improved orthogonal matching pursuit. The method comprises the following steps: firstly, carrying out digital processing on an image, and regarding each pixel in the image as an atom; then calculating correlation coefficient by the absolute value of inner product between residual r and each atom in measurement matrix phi to select L values, and storing the index values of the L values in index matrix AmAnd updating the measurement matrix phi until the residual error is small enough when the condition is met, and the finally obtained residual error is smaller than the last residual error, and reconstructing an image according to the obtained residual error to realize denoising. The invention carries out sparse coding of the signals while carrying out signal sampling, realizes sparse decomposition and reconstruction of the signals by utilizing the self sparsity of the signals, reserves the detail characteristics of the images and improves the denoising precision.

Description

Image denoising method based on improved orthogonal matching pursuit
Technical Field
The invention belongs to the technical field of image processing, and particularly relates to an image denoising method based on improved orthogonal matching pursuit.
Background
Image noise is useless information in an image, the image quality can be reduced, and the noise is often a main obstacle factor in the process of image acquisition and transmission, so that the integrity of the image is damaged, and the use value of the image is greatly influenced. Therefore, extracting useful signals to inhibit noise, improving image quality, preparing for subsequent higher-level processing, and performing denoising processing on images is an indispensable important link. The image denoising is to retain useful information in an image and eliminate or reduce interference and noise in the image.
At present, the more mature image denoising methods can be classified into the following categories:
(1) spatial domain and transform domain processing method
The spatial domain processing is to directly perform data operation on the original image, for example, directly perform data operation on the gray value of the image, and median filtering, wiener filtering and the like belong to the processing methods. The transform domain image denoising, also called frequency domain image denoising, is a very effective image denoising method based on image transform domain processing. The core idea of the method is that after certain transformation, the transformation of the domain range of the processed image from the space domain to the transformation domain is realized, and the operation of a related formula is further applied to the transformation coefficient in the transformation domain, so that the transformation coefficient is restored from the current transformation domain to the original space domain again, and the purpose of denoising is achieved.
For example, in 2013, the 'comparison and research of three spatial domain image denoising methods' of eastern Lin, studied the spatial domain image denoising method; in 2013, the research on a transform domain image denoising method is carried out by the unstable image denoising algorithm research based on the transform domain. They have in common that each point of the sampled signal, which is uniform based on the nyquist sampling theorem, is processed for each pixel in the image, resulting in redundancy of adjacent samples, with a significant waste of resources. And the same thinking and method are used, the data operation is directly carried out on the original image, and the characteristics of the image pixels are ignored. It is not difficult to understand why such methods usually accompany the more severe blurring phenomenon when denoising.
(2) Image denoising method based on orthogonal matching pursuit
The image is sparsely represented, namely, an image is considered to be an image or the image can be transformed into only a limited number of non-zero elements, most of the elements are close to zero, and the limited number of non-zero elements represent most of information of the image, so that the image can reconstruct a reconstructed image with small difference from the original image through the limited number of non-zero elements. The orthogonal matching pursuit algorithm is to use the sparse representation theory to orthogonalize the measurement matrix by an orthogonalization method, then to perform sparse approximation on the projection of the signal on the space formed by the orthogonal atoms to obtain the component and residual component of the signal on each selected atom, then to decompose the residual component by the same method, and to re-describe the signal by the finally obtained elements.
For example, in 2012, the "sparsity-based image denoising overview" of guo de pan and yanghu rain studies the orthogonal matching pursuit algorithm based on sparse representation and its early algorithm; in 2014, the fast implementation of the orthogonal matching pursuit algorithm was researched by the 'segmented regularized orthogonal matching pursuit algorithm' of Wudi, Quinumin and the like. The method realizes the orthogonal matching tracking algorithm, but does not strictly judge the matching tracking process, and excessive phenomenon can be generated in the iterative updating process of the residual error r, so that the edge in the image is unclear.
Disclosure of Invention
The invention aims to provide an image denoising method based on improved orthogonal matching pursuit, so that noise in an image is effectively removed.
The technical solution for realizing the purpose of the invention is as follows: an image denoising method based on improved orthogonal matching pursuit comprises the following steps:
step 1, collecting an image by using an infrared detector or a visible light camera, wherein the image contains noise;
step 2, if the RGB image is collected in the step 1, graying the RGB image; otherwise, performing step 3;
step 3, inputting a measurement matrix phi0Original signal y, sparsity M; initialization residual r0Initializing a support index set as y
Figure GDA0002791419920000021
Initial iteration m is 1, s is not equal to 0;
step 4, residual error r is obtainedm-1And the measurement matrix phim-1Calculating the correlation coefficient by the absolute value of the inner product between each atom, selecting L values, and storing the corresponding index value into AmIn (A)mSelecting an index set for m times;
the pass residual rm-1And the measurement matrix phim-1The absolute value of the inner product between each atom in the equation (a) to calculate the correlation coefficient u, the equation is as follows:
Figure GDA0002791419920000022
in the above formula, u is a correlation coefficient, ukFor the correlation coefficient of the residual with the k-th column of the measurement matrix,
Figure GDA0002791419920000023
is the kth column of the measurement matrix;
calculating a correlation coefficient u through an expression (1), searching index values corresponding to L maximum values from u, and storing the index values into AmPerforming the following steps;
step 5, recording the matrix corresponding to the m index value found
Figure GDA0002791419920000024
Updating a set of measurement matrices
Figure GDA0002791419920000025
Figure GDA0002791419920000026
Wherein phimFor m updated measurement matrices, phim-1The measurement matrix is updated m-1 times;
step 6, obtaining from least squares
Figure GDA0002791419920000031
Updating the mth residual rmAnd approximation of target signal x m
Figure GDA0002791419920000032
The formula is as follows:
Figure GDA0002791419920000033
Figure GDA0002791419920000034
in the above formula, the first and second carbon atoms are,
Figure GDA0002791419920000035
the mth approximation of the target signal x is carried out;
and 7, judging whether the new residual meets the conditions: if | rm2If the value is greater than epsilon, wherein epsilon is a threshold value, and m is equal to m +1, turning to the step 4;
if | rm2| < epsilon > and | < r >m2>‖rm-12Wherein r ism-1Residual error of m-1 time; then the step length is changed to s ═ s/2, and the support set A is formedmIs increased to L ═ L + s, and m ═ m +1, go to step 4;
if | rm2| < epsilon > and | < r >m2≤‖rm-12(ii) a And stopping iteration, and performing final signal reconstruction by using the obtained atoms to obtain a denoised image.
Compared with the prior art, the invention has the following remarkable advantages: (1) based on sparse representation of the image, namely, an image is considered to be changed into only a limited number of non-zero elements after being transformed, most of the elements are close to zero, and the limited number of non-zero elements represent most of information of the image, so that the image can reconstruct a reconstructed image with small difference from an original image through the limited number of non-zero elements, information needing to be reserved is greatly reduced, and resources can be effectively utilized; (2) sparse coding is carried out on the signals while signal sampling is carried out, sparse decomposition and reconstruction of the signals are realized by using the self sparsity of the signals, the neglect of the self characteristics of each pixel is avoided, and the detail characteristics of the image are reserved; (3) the method improves the phenomenon that the residual error r of the existing orthogonal matching pursuit algorithm possibly generates transition in the iterative updating process, and firstly, a support set is increased when the optimal residual error is close to the residual error, so that the precision is improved; and secondly, residual error iteration ending judgment conditions are added to avoid excess.
The present invention is described in further detail below with reference to the attached drawing figures.
Drawings
FIG. 1 is a flowchart of an image denoising method based on improved orthogonal matching pursuit according to the present invention.
Fig. 2 is an original diagram adopted by the present invention, in which (a) is an original infrared image 1, (b) is an original infrared image 2, (c) is an original visible light image 1, and (d) is an original visible light image 2.
FIG. 3 is a graph showing comparison of denoising effects of a conventional mean filtering and orthogonal matching algorithm and the algorithm of the present invention performed on an infrared image 1 having a noise variance of 30, where (a) is the noise infrared image 1, (b) is a mean filtering denoising result graph, (c) is an orthogonal matching algorithm result graph, and (d) is a denoising result graph of the method of the present invention.
FIG. 4 is a graph showing comparison of denoising effects of a conventional mean filtering and orthogonal matching algorithm and the algorithm of the present invention performed on an infrared image 1 having a noise variance of 45, where (a) is the noise infrared image 1, (b) is the b-mean filtering denoising result, (c) is the orthogonal matching algorithm result graph, and (d) is the denoising result graph of the method of the present invention.
Fig. 5 is a graph showing comparison of the denoising effects of the conventional mean filtering and orthogonal matching algorithm and the algorithm of the present invention performed on an infrared image 1 having a noise variance of 60, where (a) is the noise infrared image 1, (b) is a mean filtering denoising result graph, (c) is an orthogonal matching algorithm result graph, and (d) is a denoising result graph of the method of the present invention.
Fig. 6 is a graph showing comparison of the denoising effects of the conventional mean filtering and orthogonal matching algorithm and the algorithm of the present invention performed on an infrared image 2 having a noise variance of 30, where (a) is the noise infrared image 2, (b) is a mean filtering denoising result graph, (c) is an orthogonal matching algorithm result graph, and (d) is a denoising result graph of the method of the present invention.
Fig. 7 is a graph showing comparison of the denoising effects of the conventional mean filtering and orthogonal matching algorithm and the algorithm of the present invention performed on an infrared image 2 having a noise variance of 45, where (a) is the noise infrared image 2, (b) is a mean filtering denoising result graph, (c) is an orthogonal matching algorithm result graph, and (d) is a denoising result graph of the method of the present invention.
Fig. 8 is a graph showing comparison of the denoising effects of the conventional mean filtering and orthogonal matching algorithm and the algorithm of the present invention performed on an infrared image 2 having a noise variance of 60, where (a) is the noise infrared image 2, (b) is a mean filtering denoising result graph, (c) is an orthogonal matching algorithm result graph, and (d) is a denoising result graph of the method of the present invention.
Fig. 9 is a graph showing comparison of denoising effects of the conventional mean filtering and orthogonal matching algorithm and the algorithm of the present invention performed on a visible light image 1 having a noise variance of 30, where (a) is the noise visible light image 1, (b) is a mean filtering denoising result graph, (c) is an orthogonal matching algorithm result graph, and (d) is a denoising result graph of the method of the present invention.
Fig. 10 is a graph showing comparison of the denoising effects of the conventional mean filtering and orthogonal matching algorithm and the algorithm of the present invention performed on a visible light image 1 having a noise variance of 45, where (a) is the noise visible light image 1, (b) is a mean filtering denoising result graph, (c) is an orthogonal matching algorithm result graph, and (d) is a denoising result graph of the method of the present invention.
Fig. 11 is a graph showing comparison of the denoising effect of the conventional mean filtering and orthogonal matching algorithm and the algorithm of the present invention performed on a visible light image 1 having a noise variance of 60, where (a) is the noise visible light image 1, (b) is a mean filtering denoising result graph, (c) is an orthogonal matching algorithm result graph, and (d) is a denoising result graph of the method of the present invention.
Fig. 12 is a graph showing comparison of the denoising effect of the conventional mean filtering and orthogonal matching algorithm and the algorithm of the present invention performed on the visible light image 2 having a noise variance of 30, where (a) is the noise visible light image 2, (b) is the mean filtering denoising result, (c) is the orthogonal matching algorithm result graph, and (d) is the denoising result graph of the method of the present invention.
Fig. 13 is a graph showing comparison of the denoising effect of the conventional mean filtering and orthogonal matching algorithm and the algorithm of the present invention performed on a visible light image 2 having a noise variance of 45, where (a) is the noise visible light image 2, (b) is a mean filtering denoising result graph, (c) is an orthogonal matching algorithm result graph, and (d) is a denoising result graph of the method of the present invention.
Fig. 14 is a graph showing comparison of denoising effects of the conventional mean filtering and orthogonal matching algorithm and the algorithm of the present invention performed on a visible light image 2 having a noise variance of 60, where (a) is the noise visible light image 2, (b) is a mean filtering denoising result graph, (c) is an orthogonal matching algorithm result graph, and (d) is a denoising result graph of the method of the present invention.
Detailed Description
With reference to fig. 1, the image denoising method based on the improved orthogonal matching pursuit of the present invention includes the following steps:
step 1, collecting an image by using an infrared detector or a visible light camera, wherein the image contains noise;
step 2, if the RGB image is collected in the step 1, graying the RGB image; otherwise, performing step 3;
step 3, inputting a measurement matrix phi0Original signal y, sparsity M; initialization residual r0Initializing a support index set as y
Figure GDA0002791419920000051
Initial iteration m is 1, s is not equal to 0;
step 4, residual error r is obtainedm-1And the measurement matrix phim-1Calculating the correlation coefficient by the absolute value of the inner product between each atom, selecting L values, and storing the corresponding index value into AmIn (A)mSelecting an index set for m times;
the pass residual rm-1And the measurement matrix phim-1The absolute value of the inner product between each atom in the equation (a) to calculate the correlation coefficient u, the equation is as follows:
Figure GDA0002791419920000052
in the above formula, u is a correlation coefficient, ukFor the correlation coefficient of the residual with the k-th column of the measurement matrix,
Figure GDA0002791419920000053
is the kth column of the measurement matrix;
calculating a correlation coefficient u through an expression (1), searching index values corresponding to L maximum values from u, and storing the index values into AmPerforming the following steps;
step 5, recording the foundmatrix corresponding to index values of m times
Figure GDA0002791419920000054
Updating a set of measurement matrices
Figure GDA0002791419920000055
Figure GDA0002791419920000056
Wherein phimFor m updated measurement matrices, phim-1The measurement matrix is updated m-1 times;
step 6, obtaining from least squares
Figure GDA0002791419920000057
Updating the mth residual rmAnd approximation of target signal x m
Figure GDA0002791419920000058
The formula is as follows:
Figure GDA0002791419920000059
Figure GDA00027914199200000510
in the above formula, the first and second carbon atoms are,
Figure GDA00027914199200000511
the mth approximation of the target signal x is carried out;
and 7, judging whether the new residual meets the conditions: if | rm2If the value is greater than epsilon, wherein epsilon is a threshold value, and m is equal to m +1, turning to the step 4;
if | rm2| < epsilon > and | < r >m2>‖rm-12Wherein r ism-1Residual error of m-1 time; then the step length is changed to s ═ s/2, and the support set A is formedmIs increased to L ═ L + s, and m ═ m +1, go to step 4;
if | rm2| < epsilon > and | < r >m2≤‖rm-12(ii) a And stopping iteration, and performing final signal reconstruction by using the obtained atoms to obtain a denoised image.
1) If | rm2If the residual error is larger than epsilon, the residual error is still larger, and the image cannot be well reconstructed.
2) If | rm2| < epsilon > and | < r >m2>‖rm-12Then change the step length to s ═ s/2 support set AmIs increased to L-L + s, and m-m +1 goes to step 1. Condition | rm2Epsilon is less than or equal to ensure that the finally obtained residual error is small enough; condition | rm2≥‖rm-12And (4) explaining that the obtained residual error is worse than the previous effect, increasing the support set, and returning to the step (4) to obtain the residual error again.
3) If | rm2| < epsilon > and | < r >m2≤‖rm-12Stopping iteration, and performing final signal reconstruction by using the obtained atoms, or performing the next step. Condition | rm2Epsilon is less than or equal to ensure that the finally obtained residual error is small enough, and epsilon is a small enough quantity; condition | rm2≤‖rm-12And better values are ensured to be obtained continuously in the calculation process, and excessive calculation caused by overlarge step length and other reasons is avoided.
The invention is further described below in connection with simulated embodiments of the invention.
Example 1
Firstly, collecting infrared and visible light images by using an infrared focal plane and a control module thereof or a visible light CCD (charge coupled device), and inputting the images into a computer to obtain the infrared and visible light images; to evaluate the effectiveness of the improved orthogonal matching pursuit algorithm, simulation experiments were performed. The experiment is based on MATLAB R2014b software, and the PC is configured to be CPU 2.4GHz and memory 4G.
To denoise the infrared image and the visible image, 2 sets of experiments were performed. The original images are fig. 2(a) - (d), and the pixel sizes are all 320 × 240.
FIGS. 3(a) - (d) are graphs comparing denoising effects of a conventional mean filtering and orthogonal matching algorithm and an algorithm of the present invention for an infrared image 1 with a noise variance of 30; FIGS. 4(a) - (d) are graphs comparing denoising effects of an infrared image 1 containing a noise variance of 45, a conventional mean filtering and orthogonal matching algorithm, and an algorithm of the present invention; fig. 5(a) - (d) are graphs comparing denoising effects of the infrared image 1 containing the noise variance of 60, the conventional mean filtering and orthogonal matching algorithm, and the algorithm of the present invention.
FIGS. 6(a) - (d) are graphs comparing the denoising effect of the conventional mean filtering and orthogonal matching algorithm and the algorithm of the present invention for the infrared image 2 with the noise variance of 30; FIGS. 7(a) - (d) are graphs comparing the denoising effect of the conventional mean filtering and orthogonal matching algorithm and the algorithm of the present invention for the infrared image 2 with a noise variance of 45; fig. 8(a) - (d) are graphs comparing the denoising effect of the infrared image 2 containing the noise variance of 60, the conventional mean filtering and orthogonal matching algorithm, and the algorithm of the present invention.
FIGS. 9(a) - (d) are graphs comparing denoising effects of a visible light image 1 containing a noise variance of 30, a conventional mean filtering and orthogonal matching algorithm, and an algorithm of the present invention; FIGS. 10(a) - (d) are graphs comparing denoising effects of a visible light image 1 containing a noise variance of 45, a conventional mean filtering and orthogonal matching algorithm, and an algorithm of the present invention; FIGS. 11(a) - (d) are graphs comparing the noise removal effect of a visible light image 1 with a noise variance of 60, a conventional mean filtering and orthogonal matching algorithm, and an algorithm of the present invention
FIGS. 12(a) - (d) are graphs comparing the noise removal effect of the conventional mean filtering and orthogonal matching algorithm and the algorithm of the present invention for the visible light image 2 with the noise variance of 30; fig. 13(a) - (d) are graphs comparing the denoising effect of the conventional mean filtering and orthogonal matching algorithm and the algorithm of the present invention for the visible light image 2 with a noise variance of 45; fig. 14(a) to (d) are graphs showing comparison of the noise removal effect of the conventional mean filtering and orthogonal matching algorithm and the algorithm of the present invention, including the visible light image 2 with the noise variance of 60.
From the denoising results of fig. 3 to fig. 14, it can be found that: the 3 image denoising methods all achieve the purpose of image noise removal, except that an image processed by using the traditional spatial domain mean value filtering still contains more noise and can generate a fuzzy phenomenon, the image processed by using the orthogonal matching pursuit algorithm and the image processed by the algorithm have better effects than the mean value filtering, and the difference is that the image scene edge of the denoising effect image is clearer, such as the boundary between a person and a background, the boundary between the scene edge in the image processed by the algorithm is obvious and the definition is high.

Claims (1)

1. An image denoising method based on improved orthogonal matching pursuit is characterized by comprising the following steps:
step 1, collecting an image by using an infrared detector or a visible light camera, wherein the image contains noise;
step 2, if the RGB image is collected in the step 1, graying the RGB image; otherwise, performing step 3;
step 3, inputting a measurement matrix phi0Original signal y, sparsity M; initialization residual r0Initializing a support index set as y
Figure FDA0002791419910000011
Initial iteration m is 1, s is not equal to 0;
step 4, residual error r is obtainedm-1And the measurement matrix phim-1Calculating the correlation coefficient by the absolute value of the inner product between each atom, selecting L values, and storing the corresponding index value into AmIn (A)mSelecting an index set for m times;
the pass residual rm-1And the measurement matrix phim-1The absolute value of the inner product between each atom in the equation (a) to calculate the correlation coefficient u, the equation is as follows:
Figure FDA0002791419910000012
in the above formula, u is a correlation coefficient, ukFor the correlation of the residual with the kth column of the measurement matrixThe number of the first and second groups is,
Figure FDA0002791419910000013
is the kth column of the measurement matrix;
calculating a correlation coefficient u through an expression (1), searching index values corresponding to L maximum values from u, and storing the index values into AmPerforming the following steps;
step 5, recording the matrix corresponding to the m index value found
Figure FDA0002791419910000014
Updating a set of measurement matrices
Figure FDA0002791419910000015
Figure FDA0002791419910000016
Wherein phimFor m updated measurement matrices, phim-1The measurement matrix is updated m-1 times;
step 6, obtaining from least squares
Figure FDA0002791419910000017
Updating the mth residual rmAnd approximation of target signal x m
Figure FDA0002791419910000018
The formula is as follows:
Figure FDA0002791419910000019
Figure FDA00027914199100000110
in the above formula, the first and second carbon atoms are,
Figure FDA00027914199100000111
the mth approximation of the target signal x is carried out;
and 7, judging whether the new residual meets the conditions: if | rm2If the value is greater than epsilon, wherein epsilon is a threshold value, and m is equal to m +1, turning to the step 4;
if | rm2| < epsilon > and | < r >m2>‖rm-12Wherein r ism-1Residual error of m-1 time; then the step length is changed to s ═ s/2, and the support set A is formedmIs increased to L ═ L + s, and m ═ m +1, go to step 4;
if | rm2| < epsilon > and | < r >m2≤‖rm-12(ii) a And stopping iteration, and performing final signal reconstruction by using the obtained atoms to obtain a denoised image.
CN201710868505.3A 2017-09-22 2017-09-22 Image denoising method based on improved orthogonal matching pursuit Active CN107644406B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201710868505.3A CN107644406B (en) 2017-09-22 2017-09-22 Image denoising method based on improved orthogonal matching pursuit

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201710868505.3A CN107644406B (en) 2017-09-22 2017-09-22 Image denoising method based on improved orthogonal matching pursuit

Publications (2)

Publication Number Publication Date
CN107644406A CN107644406A (en) 2018-01-30
CN107644406B true CN107644406B (en) 2021-05-07

Family

ID=61114072

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201710868505.3A Active CN107644406B (en) 2017-09-22 2017-09-22 Image denoising method based on improved orthogonal matching pursuit

Country Status (1)

Country Link
CN (1) CN107644406B (en)

Families Citing this family (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112083508B (en) * 2020-07-22 2021-08-13 山东大学 Artificial source electromagnetic exploration signal denoising method and system with noise reference channel
CN115564680B (en) * 2022-10-17 2023-03-28 陕西师范大学 Image denoising method based on two-dimensional multipath matching pursuit algorithm

Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103218791A (en) * 2013-05-05 2013-07-24 西安电子科技大学 Image de-noising method based on sparse self-adapted dictionary
CN104601176A (en) * 2014-11-28 2015-05-06 江苏大学 Regularized adaptive matching pursuit based electric system data reconstruction method
CN105281780A (en) * 2015-11-20 2016-01-27 重庆大学 Variable step size regularized adaptive compressed sampling matching pursuit method
CN106656199A (en) * 2016-11-22 2017-05-10 南开大学 Binary inner product orthogonal matching pursuit algorithm based on compressed sensing
CN106707243A (en) * 2017-01-18 2017-05-24 哈尔滨工业大学 Generalized ROMP (Regularized Orthogonal Matching Pursuit) method for reconstructing radar signals
CN106960420A (en) * 2017-02-20 2017-07-18 南京邮电大学 A kind of image reconstructing method of segment iteration matching pursuit algorithm
CN107330946A (en) * 2017-06-05 2017-11-07 中国农业大学 A kind of image processing method and device based on compressed sensing

Patent Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103218791A (en) * 2013-05-05 2013-07-24 西安电子科技大学 Image de-noising method based on sparse self-adapted dictionary
CN104601176A (en) * 2014-11-28 2015-05-06 江苏大学 Regularized adaptive matching pursuit based electric system data reconstruction method
CN105281780A (en) * 2015-11-20 2016-01-27 重庆大学 Variable step size regularized adaptive compressed sampling matching pursuit method
CN106656199A (en) * 2016-11-22 2017-05-10 南开大学 Binary inner product orthogonal matching pursuit algorithm based on compressed sensing
CN106707243A (en) * 2017-01-18 2017-05-24 哈尔滨工业大学 Generalized ROMP (Regularized Orthogonal Matching Pursuit) method for reconstructing radar signals
CN106960420A (en) * 2017-02-20 2017-07-18 南京邮电大学 A kind of image reconstructing method of segment iteration matching pursuit algorithm
CN107330946A (en) * 2017-06-05 2017-11-07 中国农业大学 A kind of image processing method and device based on compressed sensing

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
"Sparse Solution of Underdetermined Systems of Linear Equations by Stagewise Orthogonal Matching Pursuit";David L. 等;《IEEE TRANSACTIONS ON INFORMATION THEORY》;20120228;第58卷(第2期);第1094-1121页 *
"基于压缩感知的变步长自适应匹配追踪重建算法";高睿 等;《光学学报》;20100630;第30卷(第6期);第1639-1644页 *

Also Published As

Publication number Publication date
CN107644406A (en) 2018-01-30

Similar Documents

Publication Publication Date Title
CN109102477B (en) Hyperspectral remote sensing image recovery method based on non-convex low-rank sparse constraint
CN110675347B (en) Image blind restoration method based on group sparse representation
CN108133465B (en) Non-convex low-rank relaxation hyperspectral image recovery method based on spatial spectrum weighted TV
CN111369457B (en) Remote sensing image denoising method for sparse discrimination tensor robustness PCA
CN101944230B (en) Multi-scale-based natural image non-local mean noise reduction method
CN109360157B (en) TV and wavelet regularization-based spatial variation blurred image restoration method
CN103473755B (en) Based on the sparse denoising method of SAR image that change detects
CN109636722B (en) Method for reconstructing super-resolution of online dictionary learning based on sparse representation
Deivalakshmi et al. Removal of high density salt and pepper noise through improved tolerance based selective arithmetic mean filtering with wavelet thresholding
CN107644406B (en) Image denoising method based on improved orthogonal matching pursuit
CN105590296B (en) A kind of single-frame images Super-Resolution method based on doubledictionary study
CN106981046B (en) Single image super resolution ratio reconstruction method based on multi-gradient constrained regression
CN112837220B (en) Method for improving resolution of infrared image and application thereof
CN111640080B (en) CS image denoising reconstruction method based on hyperspectral total variation
Banothu et al. High-order total bounded variation approach for gaussian noise and blur removal
CN108765350B (en) Aerospace-oriented optical remote sensing image quantization filtering method
CN105321156B (en) A kind of image recovery method based on multi-factor structure
CN110335225B (en) Method for fusing infrared light image and visible light image
Jeyaprakash et al. Linearly uncorrelated principal component and deep convolutional image deblurring for natural images
CN108198139B (en) Image denoising method based on conduction total variation regularization
Liu et al. The validity analysis of the non-local mean filter and a derived novel denoising method
Gopi et al. Undecimated double density dual tree wavelet transform based image denoising using a subband adaptive threshold
Agrawal et al. K-means clustering for adaptive wavelet based image denoising
CN110264429B (en) Image enhancement method based on sparsity and similarity priori
Jia et al. Seismic section image detail enhancement method based on bilateral texture filtering and adaptive enhancement of texture details

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