CN102034227A - Method for de-noising image - Google Patents

Method for de-noising image Download PDF

Info

Publication number
CN102034227A
CN102034227A CN 201010610864 CN201010610864A CN102034227A CN 102034227 A CN102034227 A CN 102034227A CN 201010610864 CN201010610864 CN 201010610864 CN 201010610864 A CN201010610864 A CN 201010610864A CN 102034227 A CN102034227 A CN 102034227A
Authority
CN
China
Prior art keywords
image
window
gray
denoising
threshold range
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
CN 201010610864
Other languages
Chinese (zh)
Other versions
CN102034227B (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.)
Sichuan Jiuzhou Investment Holding Group Co.,Ltd.
Original Assignee
Sichuan Jiuzhou Electric Group Co Ltd
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 Sichuan Jiuzhou Electric Group Co Ltd filed Critical Sichuan Jiuzhou Electric Group Co Ltd
Priority to CN2010106108647A priority Critical patent/CN102034227B/en
Publication of CN102034227A publication Critical patent/CN102034227A/en
Application granted granted Critical
Publication of CN102034227B publication Critical patent/CN102034227B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Landscapes

  • Image Processing (AREA)

Abstract

The invention relates to a method for processing a video image, and discloses a method for de-noising an image. The method comprises the following steps of: A, calculating the joint probability distribution P(I,J) of the image; B, determining the threshold range [LI,UI] of gray scale according to the P(I,J); C, ensuring that Wi,j is a N*N window, the coordinate of a central point of the window is (i,j), x is used for representing the central pixel point of the window Wi,j, and the window Wi,j moves to the top left corner of the image; and D, sequencing N*N pixel points in the Wi,j according to gray scale values, finding a gray scale value J corresponding to a median (N*N+1)/2 of the sequenced gray scale values, not processing the pixel point if J is in the threshold range [LI,UI] in the step B, and replacing the central point x of the window Wi,j with the gray scale value J if J is not in the threshold range [LI,UI]. The method has the advantage of solving the problems that image detail information is lost and the image blurs when pulse noise is removed by the conventional filtering method.

Description

A kind of method of image denoising
Technical field
The present invention relates to method of video image processing, relate in particular to a kind of method of image denoising.
Background technology
In intelligent transportation field widespread use video analysis and recognition technology, video analysis and identification (video analyzing and recognition) technology refers to use a computer and carry out computing and analysis by computational algorithm software from video, extract the useful information in the video, finish a technology of this information extraction and understanding, noise in the image can hinder the understanding of people to image, simultaneously also can hinder computer vision technique (CV, Computer Vision) accurately extracts useful information in the video, and the purpose of image denoising is exactly the noise of removing in the image, improve the awareness of people, improve image is done further useful information extraction and understood image; The video source of high-quality is provided for relevant Video Applications.Medium filtering denoising method of the prior art is a kind of nonlinear smoothing technology, and the gray-scale value of its each picture element is set to this intermediate value of putting all the picture element gray-scale values in certain neighborhood window.Taking out the odd number data by certain sample window from image sorts; Replacing data to be processed with the intermediate value after the ordering gets final product.But this method can cause losing of image detail information in removal impulsive noise, thereby image is thickened.And can in filtering noise, keep the detailed information of image effectively based on the impulsive noise filtering method that noise spot detects.The impulsive noise filtering method key that detects based on noise spot is: the one, and the detection of noise spot; The 2nd, to the filtering of noise spot.
Summary of the invention
What the present invention had overcome prior art can cause losing of image detail information removing impulsive noise, and therefore the problem that image thickens is necessary to provide a kind of method of image denoising.
A kind of image denoising method, it comprises following steps: steps A: the joint probability distribution P of computed image (I, J); Step B: (I J) determines the threshold range [LI, UI] of gray scale according to P; Step C:W I, jBe the window of N*N, its center point coordinate is that (i, j), x is used for representing window W I, jThe center pixel point, described window W I, jMove to the upper left corner of image; Step D: to W I, jIn N*N picture element sort according to the size of gray-scale value, find out intermediate value (N*N+1)/2 corresponding gray J of ordering gray-scale value, do not handle if J drops in the threshold range [LI, UI] among the step B this point; Otherwise window W I, jCentral point x replaces with gray-scale value J.
Described method also comprises step e: window W I, jMove right one, if window W I, jThe low order end that moves to image then skips to next line; If window W I, jIf do not exceed the lower right corner of image, then go to step D and continue to carry out, otherwise whole denoising process finishes.
Described window W I, jThe upper left corner that moves to image is meant window W I, jThe left edge of a left side, coboundary and image overlap.Described method also comprises step F: the pictorial data after the output denoising.Described step B sets a threshold range for each grey level, and when whether the picture element of a certain grey level of judgement was noise spot in step C, one group of all corresponding with its grey level threshold range compared.Described window W I, jOutwards continuation, the line number of continuation are (N-1)/2.
Compared with prior art, the invention has the beneficial effects as follows and avoided existing filtering method can cause losing of image detail information, the problem that image thickens removing impulsive noise.The detailed information that can in filtering noise, keep image based on the impulsive noise filtering method of noise spot detection effectively.
Description of drawings
Fig. 1 is for being step synoptic diagram of the present invention.
Embodiment
The present invention is further elaborated below in conjunction with accompanying drawing.The invention provides a kind of method of image denoising.It comprises following steps:
Steps A: the joint probability distribution P of computed image (I, J);
Step B: (I J) determines the threshold range [LI, UI] of gray scale according to P;
Step C:W I, jBe the window of N*N, its center point coordinate is that (i, j), x is used for representing window W I, jThe center pixel point, described window W I, jMove to the upper left corner of image;
Step D: to W I, jIn N*N picture element sort according to the size of gray-scale value, find out intermediate value (N*N+1)/2 corresponding gray J of ordering gray-scale value, do not handle if J drops in the threshold range [LI, UI] among the step B this point; Otherwise window W I, jCentral point x replaces with gray-scale value J.
The method of described image denoising also comprises step e: window W I, jMove right one, if window W I, jThe low order end that moves to image then skips to next line; If window W I, jIf do not exceed the lower right corner of image, then go to step D and continue to carry out, otherwise whole denoising process finishes.
Described window W I, jThe upper left corner that moves to image is meant window W I, jThe left edge of a left side, coboundary and image overlap.The present invention also comprises step F: the pictorial data after the output denoising.
Preferably, all set a threshold range for each grey level among the step B, when whether the picture element of a certain grey level of judgement was noise spot in step C, one group of all corresponding with its grey level threshold range compared.By such mode, can reach the actual value of going back original image better.For example, 8 256 color shade figure have 256 groups of threshold values, and each group threshold value all is made up of lower limit LI and upper limit UI.Suppose that S is 8 256 color shade figure that contain the N*N of impulsive noise, Wi, j are the windows of N*N, its center point coordinate be (i, j), x is used for representing window Wi, the center pixel point of j, its gray-scale value is I.In order to judge whether x is noise spot, at first will be for S determines 256 groups of threshold ranges, for this this item technology defined a joint probability distribution P (I, J), I is window Wi, the gray-scale value of j central pixel point x, J is the value of this medium filtering. it is defined as follows:
Figure 623586DEST_PATH_IMAGE001
Formula 1
Figure 956478DEST_PATH_IMAGE002
Formula 2
I ∈ [0,255] in the formula, J ∈ [0,255], and med(medium filtering formula) (Wi, j) expression window Wi, the median-filtered result of j.In fact (I J) has reflected distribution relation between each gray level of image and its possible intermediate value result to P.According to P (I J) determines the threshold range [LI, UI] of each gray level, LI wherein, the definition of UI:
Figure 957801DEST_PATH_IMAGE003
Formula 3
Figure 717947DEST_PATH_IMAGE004
Formula 4
In the following formula, gray level is the medium filtering value of I among a presentation video S, and on behalf of moving window length, k be, is 3 herein, and t is a threshold value, and through experiment, when t got 2, denoising effect was best.
Further among technical scheme steps C and the step D, determined threshold range [LI, UI] afterwards through steps A and step B, to judge among step C and the step D that whether x is noise spot, at first will calculate x at corresponding window Wi, j, median-filtered result, if this result drops in the threshold value [LI, UI], judge that then this point is an image useful information point, if the result drops on threshold value [LI, UI] outside, judge that then this point is noise spot, then filtering.
Further among the technical scheme steps E, window Wi, j (then jump to next line if moved to the low order end of image) to the right and move one, if do not exceed the lower right corner of image, then go to step D and continue to carry out, otherwise whole filtering finishes.
Below with the gray-scale map of 8 512*512 as experimental subjects, added 10%, 20%, 30% low-light (level) noise respectively, Wi, j are the windows of a N*N, its center point coordinate is that (i, j), x is used for representing window Wi, the center pixel point of j.For the correct frontier point of handling image, at first tackle image and carry out expansion limit, border (this paper adopts image symmetrical to expand the limit method), and the line number on expansion limit depends on Wi, the size of N among the j, when N=3, Wi, j is the window of a 3*3, and the video in window border of 3*3 should be to extending out delegation; When N=5, Wi, j are the windows of a 5*5, and the window of 5*5 should be to extending out two row, and the line number that expands the limit in a word is (N-1)/2.Embodiments of the invention adopt the window of 3*3 should outwards extend out two row.
Expand the limit to set the value of the t in formula 3 and the formula 4 after finishing be 2 and carry out noise filtering according to steps A to the method flow of step F and handle also output image.Experimental results show that filter effect the best when t=2.
Practice shows that this noise filtering method can keep the detailed information of image effectively in filtering noise.More can embody the advantage of this method under the low-light (level) situation, compared with prior art low-light (level) video image details obtains extraordinary details reservation, for computer vision technique provides round-the-clock high-quality video image.

Claims (6)

  1. An image denoising method, it comprises following steps:
    Steps A: the joint probability distribution P of computed image (I, J);
    Step B: (I J) determines the threshold range [LI, UI] of gray scale according to P;
    Step C:Wi, j are the windows of N*N, its center point coordinate be (i, j), x is used for representing window Wi, the center pixel point of j, described window Wi, j moves to the upper left corner of image;
    Step D: to Wi, the picture element of the N*N among the j sorts according to the size of gray-scale value, finds out intermediate value (N*N+1)/2 corresponding gray J of ordering gray-scale value, if interior this point of threshold range [LI, UI] that J drops among the step B is not handled; Otherwise window Wi, j central point x replaces with gray-scale value J.
  2. Image denoising as claimed in claim 1 method, it is characterized in that described method also comprises step e: window Wi, j move right one, if window Wi, the low order end that j moves to image then skips to next line; If window Wi, j if do not exceed the lower right corner of image, then go to step D and continue to carry out, otherwise whole denoising process finish.
  3. Image denoising as claimed in claim 2 method, it is characterized in that described window Wi, the upper left corner that j moves to image is meant window Wi, the left side of the left side of j, coboundary and image, coboundary overlap.
  4. Image denoising as claimed in claim 3 method, it is characterized in that described method also comprises step F: the pictorial data of output after the denoising.
  5. Image denoising as claimed in claim 4 method, it is characterized in that described step B sets a threshold range for each grey level, when whether the picture element of a certain grey level of judgement was noise spot in step C, one group of all corresponding with its grey level threshold range compared.
  6. Image denoising as claimed in claim 5 method, it is characterized in that described window Wi, the outside continuation of j, the line number of continuation is (N-1)/2.
CN2010106108647A 2010-12-29 2010-12-29 Method for de-noising image Active CN102034227B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN2010106108647A CN102034227B (en) 2010-12-29 2010-12-29 Method for de-noising image

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN2010106108647A CN102034227B (en) 2010-12-29 2010-12-29 Method for de-noising image

Publications (2)

Publication Number Publication Date
CN102034227A true CN102034227A (en) 2011-04-27
CN102034227B CN102034227B (en) 2012-06-06

Family

ID=43887080

Family Applications (1)

Application Number Title Priority Date Filing Date
CN2010106108647A Active CN102034227B (en) 2010-12-29 2010-12-29 Method for de-noising image

Country Status (1)

Country Link
CN (1) CN102034227B (en)

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102413338A (en) * 2011-11-21 2012-04-11 苏州希图视鼎微电子有限公司 Image processing method and device of spotted noise
CN102890819A (en) * 2012-09-07 2013-01-23 浙江工业大学 Image denoising method based on pixel spatial relativity judgment
CN105976341A (en) * 2016-06-17 2016-09-28 天津工业大学 Image adaptive median filtering method
CN108198146A (en) * 2017-12-29 2018-06-22 努比亚技术有限公司 A kind of noise-reduction method, equipment and computer readable storage medium

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5012524A (en) * 1989-02-27 1991-04-30 Motorola, Inc. Automatic inspection method
CN1640113A (en) * 2002-02-28 2005-07-13 皇家飞利浦电子股份有限公司 Noise filtering in images
WO2006032855A2 (en) * 2004-09-21 2006-03-30 Isis Innovation Limited Improvements in image processing
US20060087703A1 (en) * 2004-10-26 2006-04-27 Yunqiang Chen Mutual information regularized Bayesian framework for multiple image restoration
CN101126810A (en) * 2007-09-21 2008-02-20 北京航空航天大学 Synthetic aperture radar image self-adaptive spot noise suppressing method
CN101650439A (en) * 2009-08-28 2010-02-17 西安电子科技大学 Method for detecting change of remote sensing image based on difference edge and joint probability consistency

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5012524A (en) * 1989-02-27 1991-04-30 Motorola, Inc. Automatic inspection method
CN1640113A (en) * 2002-02-28 2005-07-13 皇家飞利浦电子股份有限公司 Noise filtering in images
WO2006032855A2 (en) * 2004-09-21 2006-03-30 Isis Innovation Limited Improvements in image processing
US20060087703A1 (en) * 2004-10-26 2006-04-27 Yunqiang Chen Mutual information regularized Bayesian framework for multiple image restoration
CN101126810A (en) * 2007-09-21 2008-02-20 北京航空航天大学 Synthetic aperture radar image self-adaptive spot noise suppressing method
CN101650439A (en) * 2009-08-28 2010-02-17 西安电子科技大学 Method for detecting change of remote sensing image based on difference edge and joint probability consistency

Cited By (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102413338A (en) * 2011-11-21 2012-04-11 苏州希图视鼎微电子有限公司 Image processing method and device of spotted noise
CN102890819A (en) * 2012-09-07 2013-01-23 浙江工业大学 Image denoising method based on pixel spatial relativity judgment
CN102890819B (en) * 2012-09-07 2015-03-04 浙江工业大学 Image denoising method based on pixel spatial relativity judgment
CN105976341A (en) * 2016-06-17 2016-09-28 天津工业大学 Image adaptive median filtering method
CN105976341B (en) * 2016-06-17 2018-10-09 天津工业大学 Image adaptive median filter method
CN108198146A (en) * 2017-12-29 2018-06-22 努比亚技术有限公司 A kind of noise-reduction method, equipment and computer readable storage medium
CN108198146B (en) * 2017-12-29 2020-11-24 深圳市烨弘数码科技有限公司 Noise reduction method, equipment and computer readable storage medium

Also Published As

Publication number Publication date
CN102034227B (en) 2012-06-06

Similar Documents

Publication Publication Date Title
CN108898610B (en) Object contour extraction method based on mask-RCNN
US10108835B2 (en) Method and system for decoding two-dimensional code using weighted average gray-scale algorithm
CN108280823B (en) Method and system for detecting weak edge flaws on optical cable surface in industrial production
CN102682428B (en) Fingerprint image computer automatic mending method based on direction fields
CN109241973B (en) Full-automatic soft segmentation method for characters under texture background
CN110309765B (en) High-efficiency detection method for video moving target
CN108229247A (en) A kind of mobile vehicle detection method
CN103543277A (en) Blood type result recognition algorithm based on grey level analysis and type recognition
CN101908205B (en) Magic square coding-based median filter method
CN102034227B (en) Method for de-noising image
CN105590301B (en) The Impulsive Noise Mitigation Method of adaptive just oblique diesis window mean filter
CN115311507B (en) Building board classification method based on data processing
CN112085126B (en) Single sample target detection method focusing on classification task
US8000535B2 (en) Methods and systems for refining text segmentation results
CN112634301A (en) Equipment area image extraction method and device
CN109101985A (en) It is a kind of based on adaptive neighborhood test image mismatch point to elimination method
CN110276759A (en) A kind of bad line defect diagnostic method of Mobile phone screen based on machine vision
CN111222514A (en) Local map optimization method based on visual positioning
CN114419006A (en) Method and system for removing watermark of gray level video characters changing along with background
Shi et al. Image enhancement for degraded binary document images
CN113205136A (en) Real-time high-precision detection method for appearance defects of power adapter
CN104102911A (en) Image processing for AOI (automated optical inspection)-based bullet appearance defect detection system
CN114067122B (en) Two-stage binarization image processing method
CN105354823A (en) Tree-ring image edge extraction and segmentation system
CN107403192B (en) Multi-classifier-based rapid target detection method and system

Legal Events

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

Effective date of registration: 20210517

Address after: No. 6, Jiuhua Road, khuchuang Park, Mianyang, Sichuan

Patentee after: Sichuan Jiuzhou Investment Holding Group Co.,Ltd.

Address before: 621000 No. 16 Yuejin Road, Fucheng District, Sichuan, Mianyang

Patentee before: SICHUAN JIUZHOU ELECTRIC GROUP Co.,Ltd.

TR01 Transfer of patent right