CN105894478A - Image denoising method based on statistical local rank characteristics - Google Patents

Image denoising method based on statistical local rank characteristics Download PDF

Info

Publication number
CN105894478A
CN105894478A CN201610435536.5A CN201610435536A CN105894478A CN 105894478 A CN105894478 A CN 105894478A CN 201610435536 A CN201610435536 A CN 201610435536A CN 105894478 A CN105894478 A CN 105894478A
Authority
CN
China
Prior art keywords
image
local order
denoising
lrt
statistics
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
CN201610435536.5A
Other languages
Chinese (zh)
Other versions
CN105894478B (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.)
Shanghai Lisha Technology Co ltd
Original Assignee
Chongqing University
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 Chongqing University filed Critical Chongqing University
Priority to CN201610435536.5A priority Critical patent/CN105894478B/en
Publication of CN105894478A publication Critical patent/CN105894478A/en
Application granted granted Critical
Publication of CN105894478B publication Critical patent/CN105894478B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T5/00Image enhancement or restoration
    • G06T5/70Denoising; Smoothing

Landscapes

  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Engineering & Computer Science (AREA)
  • Theoretical Computer Science (AREA)
  • Image Processing (AREA)
  • Facsimile Image Signal Circuits (AREA)

Abstract

The invention discloses an image denoising method based on statistical local rank characteristics. The image denoising method comprises the following steps: performing local rank transformation on the image under different parameter conditions by utilizing a local rank operator, thereby obtaining positive local rank transformation and negative local rank transformation of the image; adding the local rank transformation and negative local rank transformation to obtain statistical local rank characteristics with continuous parameter changes; taking the statistical local rank characteristics as constraint conditions on the basis of an image denoising method with sparse representation, and performing primary denoising on the image; and finally, performing secondary denoising on the image by controlling the difference of the statistical local rank characteristics between the images before and after denoising, and removing the image noise, thereby obtaining a final clear image. The method has the obvious effects that compared with the traditional denoising method based on sparse representation, the method has a better denoising effect, can acquire a denoised image with high quality, and further can effectively guarantee the reliability of subsequent image processing and analyzing.

Description

Image de-noising method based on statistics local order feature
Technical field
The present invention relates to technical field of image processing, specifically, be a kind of based on statistics local order feature Image de-noising method.
Background technology
Image is obtaining and in transmitting procedure, is inevitably being affected by various noises, cause figure picture element Amount declines, it is impossible to meet the demand of subsequent treatment.In order to improve picture quality, Image Denoising Technology meet the tendency and Raw.
In recent years, the strong tools that the rarefaction representation of signal becomes higher-dimension signal acquisition, characterizes and compress. Sparse representation model assumes that the non-noise composition in image can be by rarefaction representation, and noise contribution can not be dilute Dredging and represent, researcher utilizes this characteristic of sparse representation method to carry out numerous studies.Beams etc. are strange at K- On the basis of different value decomposition algorithm, integrated structure cluster and dictionary learning, devise a kind of based on non local just Then changing the image de-noising method of rarefaction representation, compare traditional K-singular value decomposition method, the method can be more preferably Ground holding image structure information (see document " Shearlet territory SAR image denoising based on rarefaction representation ", It is published in electronics and information journal, the 9th phase in 2012).Grandsons etc. use a kind of shellfish based on sparse representation model This maximum a-posteriori estimation method of leaf is removed picture noise and (is seen document " A novel image denoising Algorithm using linear Bayesian MAP estimation based on sparse representation ", send out Table is in Signal Processing, the 7th phase in 2014).But, above-mentioned image denoising based on rarefaction representation Method is all to process entire image, does not takes into account the difference of image border and non-edge, and this leads Cause denoising result and often there is the flaw on details and edge, even can lose some texture informations.
As can be seen here, need one badly and can either preferably remove picture noise, and figure can be effectively retained simultaneously As edge and the image de-noising method of detailed information.
Summary of the invention
For the deficiencies in the prior art, it is an object of the invention to provide a kind of figure based on statistics local order feature As denoising method, utilize statistics local order feature as noise constraints condition, it is possible to realize effectively removing image While smooth region noise, it is effectively retained the purpose of image edge details information.Compared with additive method, The method has more preferably denoising effect, it is possible to processing for successive image and analyzing provides the higher denoising of quality Image.
For reaching above-mentioned purpose, the technical solution used in the present invention is as follows:
A kind of image de-noising method based on statistics local order feature, it it is critical only that and follow the steps below:
Step 1: for image I, utilizes local order operator according to formula LRTk(I)={ LRTk(Ii)|Ii∈ I} is not Carry out local order conversion under the conditions of same parameter, obtain the positive local order conversion LRT of imagepk(Ii) and negative local order Conversion LRTnk(Ii), wherein, i is picture numbers, IiIt is the i-th amplitude and noise acoustic image, k=0, ± 0.01, ± 0.03 ...;
Step 2: described positive local order is converted LRTpk(Ii) and negative local order conversion LRTnk(Ii) be added, obtain The statistics local order feature of continuous parameters change
Step 3: based on a kind of image de-noising method, and use described statistics this image that locally order feature draws The constraints of denoising method carries out denoising to image, it is thus achieved that primary picture rich in detail
Step 4: by denoising formula to primary picture rich in detailCarry out second denoising process, and in denoising process The primary picture rich in detail of middle controlWith picture rich in detailBetween statistics local order feature difference, it is achieved noise Remove, it is thus achieved that final picture rich in detail
Further, described image de-noising method is based on crossing complete sparse representation model Denoising Algorithm.
Further, constraints described in step 3 is:
m i n I ~ 0 { Σ k E k ( I i ) - Σ k E k ( I ~ 0 i ) } ,
Wherein,For noise image IiStatistics local order feature,For primary clear after filtering Clear imageStatistics local order feature, i is picture numbers, k=0, ± 0.01, ± 0.03 ....
Further, the denoising formula described in step 4 is:
m i n I ~ { ( Σ k LRT p k ( I ~ i ) - Σ k LRT p k ( I ~ 0 i ) ) + ( Σ k LRT n k ( I ~ i ) - Σ k LRT n k ( I ~ 0 i ) ) }
s . t . I N 0 = D α
Wherein,For picture rich in detailJust statistics local order feature,For clear figure PictureNegative statistics local order feature,For primary picture rich in detailJust statistics local order feature,For primary picture rich in detailNegative statistics local order feature,For primary condition, D was complete Standby dictionary, α is sparse coefficient.
The present invention utilizes local order operator, image carries out under the conditions of different parameters local order conversion, obtains The positive local order conversion of image and the order conversion of negative local;By the conversion of described positive local order and negative local order conversion phase Add the statistics local order feature obtaining continuous parameters change;In traditional image denoising side based on rarefaction representation On the basis of method, using described statistics local order feature as constraints denoising first to image, and remain with figure The edge detail information of picture;By controlling before denoising and described statistics local order feature between image after denoising Difference, carries out second denoising to the image remaining with edge detail information, it is achieved the removal of noise obtains required Picture rich in detail.The present invention utilizes statistics local order feature as noise constraints condition, it is possible to realize removing image While smooth region noise, more effectively retain the purpose of image edge details information.
The remarkable result of the present invention is: utilize statistics local order feature as noise constraints condition, it is possible to realize While to the denoising of image smoothing region, more effectively retain the purpose of image edge details information.With biography The denoising method based on rarefaction representation of system is compared, and the inventive method has more preferably denoising effect, it is possible to obtain Obtain quality higher denoising image, thus ensure the validity and reliability that successive image processes and analyzes.
Accompanying drawing explanation
Fig. 1 is the method flow diagram of the present invention.
Detailed description of the invention
Detailed description of the invention and operation principle to the present invention are described in further detail below in conjunction with the accompanying drawings.
As it is shown in figure 1, a kind of image de-noising method based on statistics local order feature, enter according to following steps OK:
Initially enter step 1: for image I, utilize local order operator according to formula LRTk(I)={ LRTk(Ii)|Ii∈ I} carries out local order conversion under the conditions of different parameters, obtains the positive local of image Order conversion LRTpk(Ii) and negative local order conversion LRTnk(Ii);
In the specific implementation, it is possible to use following local order transform operator pair image carries out local order and converts:
LRTδ(I)={ LRTδ(Ii)|Ii∈ I},
Wherein,
LRT δ ( I i ) = Σ j # δ ( I i - I j ) }
# δ ( · ) = 1 , I i - I j > δ 0 , I i - I j ≤ δ
In formula, δ=kIavg, k is constant constant, IavgFor the average of pixel, threshold value in each local neighborhood δ can be on the occasion of or negative value, the corresponding different transformation results of different values.Now the local order conversion of image is as follows:
LRTk(I)={ LRTk(Ii)|Ii∈ I},
For noise image Ii, its positive local order conversion is as follows:
LRTpk(Ii), k=0 ,+0.01 ,+0.03 ...;
Its negative local order conversion is as follows:
LRTnk(Ii), k=0 ,-0.01 ,-0.03 ...
Enter step 2: in order to image border and non-edge are carried out denoising simultaneously, by described positive local order Conversion LRTpk(Ii) and negative local order conversion LRTnk(Ii) be added, the statistics local order obtaining continuous parameters change is special Levy
Secondly entrance step 3: image de-noising method based on rarefaction representation, and use described statistics local order special The constraints of this image de-noising method obtained out carries out denoising to image, it is thus achieved that primary picture rich in detailSpecific as follows:
First noise image rarefaction representation coefficient under crossing complete dictionary D is tried to achieve, as follows:
α = argminI N - Dα 2 2 + λα 1 ;
Therefore, picture rich in detailAcquisition calculated as below,Wherein, α is rarefaction representation coefficient, and λ uses Control degree of rarefication.Cross complete dictionary D and can use fixing dictionary, it would however also be possible to employ learning method obtains.
Owing to this traditional denoising model based on rarefaction representation is not by edge and the non-edge of image Separately discuss, it is possible that the phenomenon of loss in detail while removing noise, in order to avoid this phenomenon Generation, present invention introduces statistics local order feature constraint item, it is achieved denoising.Described statistics local order feature Bound term is expressed as follows:
m i n I ~ 0 { Σ k E k ( I i ) - Σ k E k ( I ~ 0 i ) }
Wherein,For noise image IiStatistics local order feature,For primary clear after filtering Clear imageStatistics local order feature, i is picture numbers, k=0, ± 0.01, ± 0.03 ....
Described statistics local order feature constraint item is intended to the statistics of the positive and negative local sum of ranks value so that filtered image Feature is infinitely close to the statistical nature of original image positive and negative local sum of ranks value.
Finally enter step 4: by denoising formula to primary picture rich in detailCarry out second denoising process, and Primary picture rich in detail is controlled during denoisingWith picture rich in detailBetween the difference of statistics local order feature, real The removal of existing noise, it is thus achieved that final picture rich in detail
Described denoising formula is:
m i n I ~ { ( Σ k LRT p k ( I ~ i ) - Σ k LRT p k ( I ~ 0 i ) ) + ( Σ k LRT n k ( I ~ i ) - Σ k LRT n k ( I ~ 0 i ) ) }
s . t . I N 0 = D α
Wherein,For picture rich in detailJust statistics local order feature,For clear figure PictureNegative statistics local order feature,For primary picture rich in detailJust statistics local order feature,For primary picture rich in detailNegative statistics local order feature,For primary condition, D was complete Standby dictionary, α is sparse coefficient.
For described image de-noising method based on rarefaction representation, use based on crossing complete sparse herein Represent model denoising method, but be not limited to this kind of method.
This programme, will statistics local order feature on the basis of traditional image de-noising method based on rarefaction representation As constraints, image is carried out first denoising, reach to retain the purpose of image edge details information;Subsequently, By controlling before denoising and the difference of described statistics local order feature between image after denoising, to remaining with edge The image of detailed information carries out second denoising, it is achieved the removal of noise.
Therefore, the present invention utilizes statistics local order feature as the constraints of noise, it is achieved that put down image Slide the removal of noise region and to image border and the reservation of detailed information.With traditional based on rarefaction representation Denoising method compare, the inventive method has more preferably denoising effect, it is possible to obtain the higher denoising of quality Image, thus the reliability that effective guarantee successive image processes and analyzes.

Claims (4)

1. an image de-noising method based on statistics local order feature, it is characterised in that enter according to following steps OK:
Step 1: for image I, utilizes local order operator according to formula LRTk(I)={ LRTk(Ii)|Ii∈ I} is not Carry out local order conversion under the conditions of same parameter, obtain the positive local order conversion LRT of imagepk(Ii) and negative local order Conversion LRTnk(Ii), wherein, i is picture numbers, IiIt is the i-th amplitude and noise acoustic image, k=0, ± 0.01, ± 0.03 ...;
Step 2: described positive local order is converted LRTpk(Ii) and negative local order conversion LRTnk(Ii) be added, obtain The statistics local order feature of continuous parameters change
Step 3: based on a kind of image de-noising method, and use described statistics this image that locally order feature draws The constraints of denoising method carries out denoising to image, it is thus achieved that primary picture rich in detail
Step 4: by denoising formula to primary picture rich in detailCarry out second denoising process, and in denoising process The primary picture rich in detail of middle controlWith picture rich in detailBetween statistics local order feature difference, it is achieved noise Remove, it is thus achieved that final picture rich in detail
Image de-noising method based on statistics local order feature the most according to claim 1, its feature exists In: described image de-noising method is based on crossing complete sparse representation model Denoising Algorithm.
Image de-noising method based on statistics local order feature the most according to claim 2, its feature exists In: constraints described in step 3 is:
m i n I ~ 0 { Σ k E k ( I i ) - Σ k E k ( I ~ 0 i ) } ,
Wherein,For noise image IiStatistics local order feature,For primary clear after filtering Clear imageStatistics local order feature, i is picture numbers, k=0, ± 0.01, ± 0.03 ....
Image de-noising method based on statistics local order feature the most according to claim 3, its feature exists In: the denoising formula described in step 4 is:
min I ~ { ( Σ k LRT p k ( I ~ i ) - Σ k LRT p k ( I ~ 0 i ) ) + ( Σ k LRT n k ( I ~ i ) - Σ k LRT n k ( I ~ 0 i ) ) }
s . t . I N 0 = D α
Wherein,For picture rich in detailJust statistics local order feature,For clear figure PictureNegative statistics local order feature,For primary picture rich in detailJust statistics local order feature,For primary picture rich in detailNegative statistics local order feature,For primary condition, D was complete Standby dictionary, α is sparse coefficient.
CN201610435536.5A 2016-06-17 2016-06-17 Image de-noising method based on the local order feature of statistics Active CN105894478B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201610435536.5A CN105894478B (en) 2016-06-17 2016-06-17 Image de-noising method based on the local order feature of statistics

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201610435536.5A CN105894478B (en) 2016-06-17 2016-06-17 Image de-noising method based on the local order feature of statistics

Publications (2)

Publication Number Publication Date
CN105894478A true CN105894478A (en) 2016-08-24
CN105894478B CN105894478B (en) 2018-08-31

Family

ID=56729758

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201610435536.5A Active CN105894478B (en) 2016-06-17 2016-06-17 Image de-noising method based on the local order feature of statistics

Country Status (1)

Country Link
CN (1) CN105894478B (en)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107301632A (en) * 2017-06-28 2017-10-27 重庆大学 A kind of SAR image method for reducing speckle represented based on sequence joint sparse
WO2023279270A1 (en) * 2021-07-07 2023-01-12 Qualcomm Incorporated Cascade image processing for noise reduction

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102184526A (en) * 2011-04-22 2011-09-14 西安电子科技大学 Natural image denoising method based on dictionary learning and block matching
US20130156340A1 (en) * 2011-12-20 2013-06-20 Fatih Porikli Image Filtering by Sparse Reconstruction on Affinity Net
CN104123707A (en) * 2014-08-07 2014-10-29 重庆大学 Local rank priori based single-image super-resolution reconstruction method

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102184526A (en) * 2011-04-22 2011-09-14 西安电子科技大学 Natural image denoising method based on dictionary learning and block matching
US20130156340A1 (en) * 2011-12-20 2013-06-20 Fatih Porikli Image Filtering by Sparse Reconstruction on Affinity Net
CN104123707A (en) * 2014-08-07 2014-10-29 重庆大学 Local rank priori based single-image super-resolution reconstruction method

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
JIE TANG 等: "Image Deblurring Based on Non-Local Total Variation and Global Non-Zero Local Rank Penalty", 《OPEN JOURNAL OF NATURE SCIENCE》 *
YONG-QIANG ZHAO 等: "Hyperspectral Image Denoising via Sparse Representation and Low-Rank Constraint", 《IEEE TRANSACTIONS ON GEOSCIENCE AND REMOTE SENSING》 *
胡伦庭: "基于稀疏表示和局部秩的单幅图像超分辨率重建方法研究", 《中国优秀硕士学位论文全文数据库 信息科技辑》 *

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107301632A (en) * 2017-06-28 2017-10-27 重庆大学 A kind of SAR image method for reducing speckle represented based on sequence joint sparse
WO2023279270A1 (en) * 2021-07-07 2023-01-12 Qualcomm Incorporated Cascade image processing for noise reduction

Also Published As

Publication number Publication date
CN105894478B (en) 2018-08-31

Similar Documents

Publication Publication Date Title
CN101980284B (en) Two-scale sparse representation-based color image noise reduction method
CN100550978C (en) A kind of self-adapting method for filtering image that keeps the edge
CN103093441B (en) Based on the non-local mean of transform domain and the image de-noising method of two-varaible model
CN104715461A (en) Image noise reduction method
CN103854264A (en) Improved threshold function-based wavelet transformation image denoising method
CN101877123A (en) Image enhancement method and device
CN102663702A (en) Natural image denoising method based on regional division
CN101944230A (en) Multi-scale-based natural image non-local mean noise reduction method
CN102930508B (en) Image residual signal based non-local mean value image de-noising method
CN104835127B (en) A kind of self-adaptive smooth filtering method
CN106570843A (en) Adaptive wavelet threshold function image noise suppression method
CN103761719A (en) Self-adaptive wavelet threshold de-noising method based on neighborhood correlation
CN101504769B (en) Self-adaptive noise intensity estimation method based on encoder frame work
CN104103041A (en) Ultrasonic image mixed noise self-adaptive suppression method
Biswas et al. A new soft-thresholding image denoising method
CN105894478A (en) Image denoising method based on statistical local rank characteristics
CN103886558B (en) Improved adaptive threshold wavelet denoising algorithm based on LoG operator
CN104616259A (en) Non-local mean image de-noising method with noise intensity self-adaptation function
CN106023103A (en) Adaptive orthogonal wavelet image denoising method based on accurate local variance priori modeling
CN102314675A (en) Wavelet high-frequency-based Bayesian denoising method
CN104809714A (en) Image fusion method based on multi-morphological sparse representation
CN102184530B (en) Image denoising method based on gray relation threshold value
CN102622731B (en) Contourlet domain Wiener filtering image denoising method based on two-dimensional Otsu
CN103854258A (en) Image denoising method based on Contourlet transformation self-adaptation direction threshold value
CN108961290B (en) Adaptive image edge detection method based on Ratio operator of Otsu

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant
TR01 Transfer of patent right
TR01 Transfer of patent right

Effective date of registration: 20210513

Address after: 200333 room 5081, 5th floor, No.5 Lane 600, Yunling West Road, Putuo District, Shanghai

Patentee after: Shanghai Lisha Technology Co.,Ltd.

Address before: 400 000 No. 174 Zhengjie, Shapingba District, Chongqing

Patentee before: Chongqing University