CN102034227A - Method for de-noising image - Google Patents
Method for de-noising image Download PDFInfo
- 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
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
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:
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:
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)
- 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.
- 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.
- 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.
- 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.
- 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.
- 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.
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)
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)
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 |
-
2010
- 2010-12-29 CN CN2010106108647A patent/CN102034227B/en active Active
Patent Citations (6)
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)
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 |