CN101908208B - Self-adaptive confirming method of smooth wave-filtering spatial scale facing to picture edge detection - Google Patents

Self-adaptive confirming method of smooth wave-filtering spatial scale facing to picture edge detection Download PDF

Info

Publication number
CN101908208B
CN101908208B CN2010102379141A CN201010237914A CN101908208B CN 101908208 B CN101908208 B CN 101908208B CN 2010102379141 A CN2010102379141 A CN 2010102379141A CN 201010237914 A CN201010237914 A CN 201010237914A CN 101908208 B CN101908208 B CN 101908208B
Authority
CN
China
Prior art keywords
prime
image
sliding window
pixel
color
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.)
Expired - Fee Related
Application number
CN2010102379141A
Other languages
Chinese (zh)
Other versions
CN101908208A (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.)
Zhejiang University ZJU
Original Assignee
Zhejiang University ZJU
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 Zhejiang University ZJU filed Critical Zhejiang University ZJU
Priority to CN2010102379141A priority Critical patent/CN101908208B/en
Publication of CN101908208A publication Critical patent/CN101908208A/en
Application granted granted Critical
Publication of CN101908208B publication Critical patent/CN101908208B/en
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

  • Image Analysis (AREA)
  • Image Processing (AREA)

Abstract

本发明公开了一种面向图像边缘检测的平滑滤波空间尺度的自适应确定方法。本发明可以为图像中每个像素点自适应确定最佳的平滑滤波空间尺度,也就是落在大面积平滑均匀区域中像素点的平滑滤波空间尺度相对比较大,而处于图像边缘邻近像素点的平滑滤波空间尺度相对比较小。该自适应确定图像中各个像素点的平滑滤波空间尺度可以结合到灰度或彩色图像的边缘检测或平滑滤波中,能有效地消除图像中独立分布高斯白噪声或空间相关纹理噪声的干扰,同时保持图像中的边缘不会被过模糊,可以应用于纺织印染行业中。

Figure 201010237914

The invention discloses an adaptive determination method of smoothing filter space scale oriented to image edge detection. The present invention can adaptively determine the best smoothing filter spatial scale for each pixel in the image, that is, the smoothing filter spatial scale of the pixel in the large-area smooth and uniform area is relatively large, while the smoothing filter spatial scale of the pixel in the edge of the image is relatively large. The smoothing filter spatial scale is relatively small. The adaptive determination of the smoothing filter spatial scale of each pixel in the image can be combined with the edge detection or smoothing filter of grayscale or color images, which can effectively eliminate the interference of independently distributed Gaussian white noise or spatially correlated texture noise in the image, and at the same time Keep the edges of the image from being blurred, and can be applied in the textile printing and dyeing industry.

Figure 201010237914

Description

Self-adaptation towards the smooth wave-filtering spatial scale of Image Edge-Detection is determined method
Technical field
The present invention relates to a kind of picture signal treatment technology, the self-adaptation that relates in particular to a kind of smooth wave-filtering spatial scale towards Image Edge-Detection is determined method.
Background technology
A problem that had not solved is how to determine the smooth wave-filtering spatial scale of each pixel in Image Edge-Detection or smothing filtering, in the traditional images edge detection method, as Canny operator [1] (J.Canny, " A computational approach to edge detection; " IEEE Trans.On Pattern Analysis andMachine Intelligence, 8 (6): 679-698,1986), entire image adopts same smooth wave-filtering spatial scale, but interference of noise in the very difficult removal of images effectively of this method keeps edge of image can not be filtered ripple again simultaneously.
Though also having the researchist to adopt the whole bag of tricks is that each pixel is determined corresponding smooth wave-filtering spatial scale in the image, as Jeong and Kim[2] (H.Jeong and C.I.Kim, " Adaptive determinationof filter scales for edge detection; " IEEE Trans.On Pattern Analysis and MachineIntelligence, 14 (5): 579-585,1992) propose the employing regularization method and optimize the filtering space scale that the energy function of metric space is determined each pixel in the image; And Lindeberg[3] (T.Lindeberg, " Edgedetection and ridge detection with automatic scale selection; " in Proc.of CVPR, 1996) suggestion is that characteristics of image detects the space scale of determining feature detection etc. by the local maximum of seeking normalization feature amplitude at metric space.But all there is weak point in these methods, the gradient descending method that is adopted in the energy function optimizing process is absorbed in the local optimum point, Jeong and Kim[2] smooth wave-filtering spatial scale of each pixel is initialized as 3 in the method that proposes, and then slowly increase, but such space scale is excessive for the pixel that is in the image border, causes that easily crossing filtering because of the image border blurs; Lindeberg[3] then need the space partial derivative of higher order of each pixel in the computed image in the method, and for being subjected to the noise image, especially towards the textile printing and dyeing image of the texture noise that is subjected to space correlation, because signal to noise ratio (S/N ratio) is very low in the image, the calculating of image partial derivative of higher order is insecure.
Summary of the invention
The object of the present invention is to provide a kind of self-adaptation of the smooth wave-filtering spatial scale towards Image Edge-Detection to determine method, if utilized pixel to drop in the even smooth region of image, each sub regions in this zone then, as the left side, right-hand part, the sample mean of the first half and Lower Half gray scale or color is that whole area grayscale or color mean value one does not have estimation partially, under the situation of given confidence level 1-α, judge the left side that this is regional, right-hand part, whether the gray scale of the first half and Lower Half or the sample mean of color drop on whole area grayscale or the corresponding fiducial interval of color mean value.Can select to select best smooth wave-filtering spatial scale adaptively for each pixel in the set for smooth wave-filtering spatial scale at one group according to these characteristics.
The step of the technical solution used in the present invention is as follows:
(1) set a discrete values set selecting smooth wave-filtering spatial scale for each pixel in the image, the numerical value in this discrete set is to arrange from big to small;
(2) the image smoothing wave filter is a moving window, and according to grating lattice scan mode point by point scanning view picture input picture, input picture is gray level image or coloured image, and described smooth wave-filtering spatial scale is exactly the radius r of this moving window;
(3) when moving window scans a pixel, when promptly the center of moving window rests on this pixel, from gathering, the selective discrete values of above-mentioned smooth wave-filtering spatial scale gets the radius r of a numerical value as moving window according to order;
(4) calculate the mean value that drops on interior all pixel gray scales of this moving window or color, and the variance of gray scale or color:
c ‾ ( x , y ) = 1 ( 2 r + 1 ) 2 Σ x ′ = x - r x + r Σ y ′ = y - r y + r c ( x ′ , y ′ ) - - - ( 1 )
σ 2 ( x , y ) = 1 ( 2 r + 1 ) 2 Σ x ′ = x - r x + r Σ y ′ = y - r y + r ( c ( x ′ , y ′ ) - c ‾ ( x , y ) ) 2 - - - ( 2 )
(x wherein, y) volume coordinate of central pixel point in the expression moving window, r represents the radius of described moving window, r is a unit with pixel interval in the digital picture, (x ', y ') expression drops on that (x y) is the center, and radius is the volume coordinate of the moving window interior pixel point of r with pixel, c (x ', y ') expression is positioned at the gray scale or the color of pixel (x ', y ')
Figure BSA00000206449600023
(x y) is the center, and radius is the mean value of interior all pixel gray scales of the moving window of r or color, σ with pixel in expression 2(x, y) variance of all pixel gray scales or color in the described moving window of expression;
(5) calculating drops on described moving window left side respectively, right-hand part, and the mean value of the first half and Lower Half pixel gray scale or color:
c ‾ L ( x , y ) = 1 ( 2 r + 1 ) r Σ x ′ = x - r x - 1 Σ y ′ = y - r y + r c ( x ′ , y ′ ) - - - ( 3 )
c ‾ R ( x , y ) = 1 ( 2 r + 1 ) r Σ x ′ = x + 1 x + r Σ y ′ = y - r y + r c ( x ′ , y ′ ) - - - ( 4 )
c ‾ U ( x , y ) = 1 ( 2 r + 1 ) r Σ x ′ = x - r x + r Σ y ′ = y - r y + 1 c ( x ′ , y ′ ) - - - ( 5 )
c ‾ D ( x , y ) = 1 ( 2 r + 1 ) r Σ x ′ = x - r x + r Σ y ′ = y + 1 y + r c ( x ′ , y ′ ) - - - ( 6 )
Here With
Figure BSA00000206449600033
Expression drops on described moving window left side respectively, right-hand part, and the mean value of the first half and Lower Half pixel gray scale or color is if input picture is a gray level image, then
Figure BSA00000206449600034
With
Figure BSA00000206449600035
It all is scalar; If input picture is a coloured image, then
Figure BSA00000206449600036
With
Figure BSA00000206449600037
Be vector, each vector all has three components, and is promptly red, green and blue component;
(6) press the definition of α quantile in the t distribution on the statistics, check described moving window left side, right-hand part, the mean value of the first half and Lower Half pixel gray scale or color, promptly With
Figure BSA00000206449600039
The mean value that whether drops on whole slide window gray scale or color is in confidence level is fiducial interval under the 1-α:
P { | c &OverBar; j i ( x , y ) - c &OverBar; i ( x , y ) &sigma; i ( x , y ) / n | < z &alpha; / 2 } = 1 - &alpha; - - - ( 7 )
If input picture is a gray level image, subscript i=1 here is because have only a gray-scale value on each pixel in the image, and when input picture is coloured image, here subscript i represents redness corresponding in color mean value or the color variance, green or blue component, and
Figure BSA000002064496000311
Middle subscript j then represents it is described moving window left side, right-hand part, the first half or Lower Half, n represents to drop on described moving window left side, right-hand part, the number of pixel in the first half or the Lower Half, if described moving window is a square, n numerical value then just equals (2r+1) r, according to given α numerical value, obtains z by the inquiry of t distribution table α/2Numerical value;
(7) if described moving window left side, right-hand part, the mean value of the first half and Lower Half pixel gray scale or color, promptly
Figure BSA000002064496000312
With
Figure BSA000002064496000313
In each component, all satisfy formula (7), judge that then (x y) drops on a length of side and is the homogeneous area of (2r+1), jumps to step (9) for the pixel at place, described moving window center;
(8) if described moving window left side, right-hand part, the mean value of the first half and Lower Half pixel gray scale or color, promptly
Figure BSA000002064496000314
With
Figure BSA000002064496000315
In each component have one not satisfy formula (7), the image-region that current moving window place is described is not a homogeneous area, then next than the radius r of fractional value as moving window according to order selection from the alternative discrete values set of above-mentioned smooth wave-filtering spatial scale from big to small, repeat above-mentioned steps (4) then to step (6);
(9) the current moving window radius r of getting is as the smooth wave-filtering spatial scale of the pixel at place, moving window center;
(10) described moving window moves on to next pixel according to the grating lattice scanning sequency, and repeating step (3) is to step (9), and all pixels are all determined smooth wave-filtering spatial scale separately in image.
The beneficial effect that the present invention has is:
This self-adaptation determines that the smooth wave-filtering spatial scale of each pixel can be incorporated in the rim detection or smothing filtering of gray scale or coloured image in the image, independent distribution white Gaussian noise or space correlation texture interference of noise in the removal of images effectively, keep the edge in the image can not blured excessively simultaneously, can be applied in the textile printing and dyeing industry.
Description of drawings
Accompanying drawing is a process flow diagram of the present invention.
Embodiment
As shown in drawings, workflow of the present invention is as follows:
(1) input one width of cloth numerical imaging.
(2) the given discrete values set that can select smooth wave-filtering spatial scale for each pixel in the image, as 6,5,4,3,2,1}.
(3) input confidence level 1-α is 0.05 as the α value.
(4) the image smoothing wave filter is a square moving window, and the radius r of moving window is equal to smooth wave-filtering spatial scale.
(5) maximal value desirable according to the radius r of moving window, as 6 in this example, image is done symmetry expansion up and down accordingly, if the width of input picture and highly be respectively Width and Height, expand the width of back image so and highly be respectively Width+2r and Height+2r, can be earlier the row of image be carried out the symmetry expansion, and then row are carried out symmetry expand, promptly
I row - ext c ( i , j ) = I c ( i , r + 2 - j ) ; 1 &le; i &le; Height , 1 &le; j &le; r I c ( i , j - r ) ; 1 &le; i &le; Height , r + 1 &le; j &le; r + Width I c ( i , 2 Width + r - j ) ; 1 &le; i &le; Height , r + Width + 1 &le; j &le; 2 r + Width - - - ( 1 )
I ext c ( i , j ) = I row - ext c ( w + 2 - i , j ) ; 1 &le; i &le; w , 1 &le; j &le; 2 w + Width I row - ext c ( i - w , j ) ; w + 1 &le; i &le; w + Height , 1 &le; j &le; 2 w + Width I row - ext c ( 2 Height + w - i , j ) ; w + 1 + Height &le; i &le; 2 w + Height , 1 &le; j &le; 2 w + Width - - - ( 2 )
Here I cBe a certain component of expression input picture,, just have only one-component if input picture is a gray level image; If input picture is a coloured image, then at redness, green or the blue component of rgb color space.Here (i, j) volume coordinate of pixel in the image of back is expanded in expression.
Figure BSA00000206449600043
Be to component I cThe symmetry expansion of line direction, and
Figure BSA00000206449600044
Be
Figure BSA00000206449600045
Again column direction is carried out symmetry expansion back image on the basis.
Perhaps also can be earlier the row of image be carried out the symmetry expansion, and then row is carried out symmetry expand, promptly
I col - ext c ( i , j ) = I c ( r + 2 - i , j ) ; 1 &le; i &le; r , 1 &le; j &le; Width I c ( i - r , j ) ; r + 1 &le; i &le; r + Height , 1 &le; j &le; Width I c ( 2 Height + r - i , j ) ; r + Height + 1 &le; i &le; 2 r + Height , 1 &le; j &le; Width - - - ( 3 )
I ext c ( i , j ) = I col - ext c ( i , r + 2 - j ) ; 1 &le; i &le; 2 r + Height , 1 &le; j &le; r I col - ext c ( i , j - r ) ; 1 &le; i &le; 2 r + Height , 1 + r &le; j &le; r + Width I col - ext c ( i , 2 r + Width - j ) ; 1 &le; i &le; 2 r + Height , r + 1 + Width &le; j &le; 2 r + Width - - - ( 4 )
Wherein
Figure BSA00000206449600053
Be to component I cSymmetry expansion on column direction, and Be
Figure BSA00000206449600055
On the basis line direction is carried out image after the symmetry expansion, the expanded images that two kinds of symmetrical extended modes obtain is the same.
(6) moving window is according to inner the work from top to bottom of grating lattice scan mode image after expansion, zigzag scanning from left to right, promptly remove the image inside of moving window maximum radius size frame, the size of expanded images inside is the same with the original size of input picture.
(7) when moving window scans a pixel, when promptly the center of moving window rests on this pixel, get the current radius r of a numerical value according to order as moving window from the alternative discrete values set of above-mentioned smooth wave-filtering spatial scale;
(8) calculate the mean value that drops on this moving window interior pixel point gray scale or color, and gray scale or color variance
c &OverBar; ( x , y ) = 1 ( 2 r + 1 ) 2 &Sigma; x &prime; = x - r x + r &Sigma; y &prime; = y - r y + r c ( x &prime; , y &prime; ) - - - ( 5 )
&sigma; 2 ( x , y ) = 1 ( 2 r + 1 ) 2 &Sigma; x &prime; = x - r x + r &Sigma; y &prime; = y - r y + r ( c ( x &prime; , y &prime; ) - c &OverBar; ( x , y ) ) 2 - - - ( 6 )
(x wherein, y) volume coordinate of central pixel point in the expression moving window, r represents the current numerical value r that gets radius of described moving window, r is a unit with pixel interval in the digital picture, (x ', y ') expression drops on that (x y) is the center, and radius is the volume coordinate of the moving window interior pixels point of r with pixel, c (x ', y ') expression is positioned at the gray scale or the color of pixel (x ', y ') (x y) is the center, and radius is the mean value of interior all pixel gray scales of the moving window of r or color, σ with pixel in expression 2(x, y) variance of inner all pixel gray scales of the described moving window of expression or color.If input picture is a gray level image, then
Figure BSA00000206449600059
And σ 2(x y) is scalar; If input picture is a coloured image, then
Figure BSA000002064496000510
And σ 2(x y) is vector, and each vector all has three components, and is promptly red, green and blue component.
(9) calculating drops on described moving window left side respectively, right-hand part, and the mean value of the first half and Lower Half pixel gray scale or color:
c &OverBar; L ( x , y ) = 1 ( 2 r + 1 ) r &Sigma; x &prime; = x - r x - 1 &Sigma; y &prime; = y - r y + r c ( x &prime; , y &prime; ) - - - ( 7 )
c &OverBar; R ( x , y ) = 1 ( 2 r + 1 ) r &Sigma; x &prime; = x + 1 x + r &Sigma; y &prime; = y - r y + r c ( x &prime; , y &prime; ) - - - ( 8 )
c &OverBar; U ( x , y ) = 1 ( 2 r + 1 ) r &Sigma; x &prime; = x - r x + r &Sigma; y &prime; = y - r y + 1 c ( x &prime; , y &prime; ) - - - ( 9 )
c &OverBar; D ( x , y ) = 1 ( 2 r + 1 ) r &Sigma; x &prime; = x - r x + r &Sigma; y &prime; = y + 1 y + r c ( x &prime; , y &prime; ) - - - ( 10 )
Here
Figure BSA00000206449600065
With
Figure BSA00000206449600066
Expression drops on described moving window left side respectively, right-hand part, and the mean value of the first half and Lower Half pixel gray scale or color is if input picture is a gray level image, then
Figure BSA00000206449600067
With
Figure BSA00000206449600068
Be scalar, and if input picture is a coloured image, then
Figure BSA00000206449600069
With
Figure BSA000002064496000610
Be vector, each all has three components, and is promptly red, green and blue component.
(10) press the definition of α quantile in the statistics t distribution, check described moving window left side, right-hand part, the mean value of the first half and Lower Half pixel gray scale or color, promptly
Figure BSA000002064496000611
With
Figure BSA000002064496000612
Whether drop in the fiducial interval that confidence level is 1-α:
P { | c &OverBar; j i ( x , y ) - c &OverBar; i ( x , y ) &sigma; i ( x , y ) / n | < z &alpha; / 2 } = 1 - &alpha; - - - ( 11 )
If input picture is a gray level image, subscript i=1 here is because have only a gray-scale value on each pixel in the image, and when input picture is coloured image, here subscript i represents corresponding redness in color mean value or the color variance, green or blue component, and
Figure BSA000002064496000614
Middle subscript j then represents it is described moving window left side, right-hand part, the first half or Lower Half, n represents to drop on described moving window left side, right-hand part, the number of pixel in the first half or the Lower Half, if described moving window is a square, n numerical value then just equals (2r+1) r, according to given α numerical value, can obtain z by the inquiry of statistics t distribution table α/2Numerical value.
(11) if described moving window left side, right-hand part, the mean value of the first half and Lower Half pixel gray scale or color, promptly
Figure BSA000002064496000615
With In each component, all satisfy formula (11), judge that then (x y) drops on a length of side and is the homogeneous area of (2r+1), jumps to step (13) for the pixel at place, described moving window center.
(12) if described moving window left side, right-hand part, the mean value of the first half and Lower Half pixel gray scale or color, promptly
Figure BSA00000206449600071
With
Figure BSA00000206449600072
In each component have one not satisfy formula (11), the image-region that current moving window place is described is not a homogeneous area, then next than the radius r of fractional value as moving window according to order selection from the alternative discrete values set of above-mentioned smoothing filter space scale from big to small, repeat above-mentioned steps (8) then to (10).
(13) the current moving window radius r of getting is as the smooth wave-filtering spatial scale of the pixel at place, moving window center.
(14) described moving window moves on to next pixel according to the grating lattice scanning sequency, and repeating step (7) is to step (13), and all pixels are all determined smooth wave-filtering spatial scale separately in image.
(15), can be incorporated in Image Edge-Detection or the smothing filtering according to the smooth wave-filtering spatial scale of each pixel.If in conjunction with method for detecting image edge, can be earlier image be carried out smothing filtering according to the smooth wave-filtering spatial scale of each pixel the best based on gradient; On the basis of removing picture noise, the space derivative of computed image; Be steps such as local non-maximum value inhibition and hysteresis threshold processing then, just can obtain edge of image.

Claims (1)

1.一种面向图像边缘检测的平滑滤波空间尺度的自适应确定方法,其特征在于该方法的步骤如下:1. an adaptive method for determining the smoothing filter spatial scale facing image edge detection, characterized in that the steps of the method are as follows: (1)设定供图像中每个像素点选择平滑滤波空间尺度的一个离散数值集合,这个离散集合中的数值是从大到小排列;(1) Set a discrete value set for each pixel in the image to select a smoothing filter spatial scale, and the values in this discrete set are arranged from large to small; (2)图像平滑滤波器是一个滑动窗口,按照光栅点阵扫描方式逐点扫描整幅输入图像,输入图像是灰度图像或彩色图像,而所述的平滑滤波空间尺度就是这个滑动窗口的半径r;(2) The image smoothing filter is a sliding window, which scans the entire input image point by point according to the raster lattice scanning method, the input image is a grayscale image or a color image, and the smoothing filter spatial scale is the radius of the sliding window r; (3)当滑动窗口扫描到一个像素点时,即滑动窗口的中心停留在该像素点时,从上述平滑滤波空间尺度供选择的离散数值集合中按照次序取一个数值作为滑动窗口的半径r;(3) When the sliding window scans to a pixel point, that is, when the center of the sliding window stays at the pixel point, a numerical value is taken as the radius r of the sliding window in order from the set of discrete numerical values for selection of the above-mentioned smoothing filter spatial scale; (4)计算落在这个滑动窗口内所有像素点灰度或色彩的平均值,以及灰度或色彩的方差:(4) Calculate the average value of the grayscale or color of all pixels falling within the sliding window, and the variance of the grayscale or color: cc &OverBar;&OverBar; (( xx ,, ythe y )) == 11 (( 22 rr ++ 11 )) 22 &Sigma;&Sigma; xx &prime;&prime; == xx -- rr xx ++ rr &Sigma;&Sigma; ythe y &prime;&prime; == ythe y -- rr ythe y ++ rr cc (( xx &prime;&prime; ,, ythe y &prime;&prime; )) -- -- -- (( 11 )) &sigma;&sigma; 22 (( xx ,, ythe y )) == 11 (( 22 rr ++ 11 )) 22 &Sigma;&Sigma; xx &prime;&prime; == xx -- rr xx ++ rr &Sigma;&Sigma; ythe y &prime;&prime; == ythe y -- rr ythe y ++ rr (( cc (( xx &prime;&prime; ,, ythe y &prime;&prime; )) -- cc &OverBar;&OverBar; (( xx ,, ythe y )) )) 22 -- -- -- (( 22 )) 其中(x,y)表示滑动窗口内中心像素点的空间坐标,r表示所述滑动窗口的半径,r以数字图像中像素点间距为单位,(x′,y′)表示落在以像素点(x,y)为中心,半径为r的滑动窗口内像素点的空间坐标,c(x′,y′)表示位于像素点(x′,y′)的灰度或色彩,
Figure FSB00000580858800013
表示以像素点(x,y)为中心,半径为r的滑动窗口内所有像素点灰度或色彩的平均值,σ2(x,y)表示所述滑动窗口内所有像素点灰度或色彩的方差;
Wherein (x, y) represents the spatial coordinates of the central pixel in the sliding window, r represents the radius of the sliding window, r is in the pixel pitch in the digital image, and (x', y') represents the distance between pixels (x, y) is the center, the spatial coordinates of the pixel in the sliding window with a radius of r, c(x', y') represents the grayscale or color of the pixel (x', y'),
Figure FSB00000580858800013
Indicates the average value of the grayscale or color of all pixels in the sliding window with the pixel (x, y) as the center and radius r, σ 2 (x, y) represents the grayscale or color of all pixels in the sliding window Variance;
(5)分别计算落在所述滑动窗口左半部,右半部,上半部和下半部像素点灰度或色彩的平均值:(5) Calculate the average value of pixel grayscale or color falling on the left half of the sliding window, the right half, the upper half and the lower half respectively: cc &OverBar;&OverBar; LL (( xx ,, ythe y )) == 11 (( 22 rr ++ 11 )) rr &Sigma;&Sigma; xx &prime;&prime; == xx -- rr xx -- 11 &Sigma;&Sigma; ythe y &prime;&prime; == ythe y -- rr ythe y ++ rr cc (( xx &prime;&prime; ,, ythe y &prime;&prime; )) -- -- -- (( 33 )) cc &OverBar;&OverBar; RR (( xx ,, ythe y )) == 11 (( 22 rr ++ 11 )) rr &Sigma;&Sigma; xx &prime;&prime; == xx ++ 11 xx ++ rr &Sigma;&Sigma; ythe y &prime;&prime; == ythe y -- rr ythe y ++ rr cc (( xx &prime;&prime; ,, ythe y &prime;&prime; )) -- -- -- (( 44 )) cc &OverBar;&OverBar; Uu (( xx ,, ythe y )) == 11 (( 22 rr ++ 11 )) rr &Sigma;&Sigma; xx &prime;&prime; == xx -- rr xx ++ rr &Sigma;&Sigma; ythe y &prime;&prime; == ythe y -- rr ythe y -- 11 cc (( xx &prime;&prime; ,, ythe y &prime;&prime; )) -- -- -- (( 55 )) cc &OverBar;&OverBar; DD. (( xx ,, ythe y )) == 11 (( 22 rr ++ 11 )) rr &Sigma;&Sigma; xx &prime;&prime; == xx -- rr xx ++ rr &Sigma;&Sigma; ythe y &prime;&prime; == ythe y ++ 11 ythe y ++ rr cc (( xx &prime;&prime; ,, ythe y &prime;&prime; )) -- -- -- (( 66 )) 这里
Figure FSB00000580858800022
Figure FSB00000580858800023
分别表示落在所述滑动窗口左半部,右半部,上半部和下半部像素点灰度或色彩的平均值,如果输入图像是灰度图像,则
Figure FSB00000580858800024
Figure FSB00000580858800025
都是标量;如果输入图像是彩色图像,则
Figure FSB00000580858800026
Figure FSB00000580858800027
是矢量,每个矢量都有三个分量,即红色,绿色和蓝色分量;
here
Figure FSB00000580858800022
and
Figure FSB00000580858800023
Represents the average value of the grayscale or color of the pixels falling in the left half, right half, upper half and lower half of the sliding window, if the input image is a grayscale image, then
Figure FSB00000580858800024
and
Figure FSB00000580858800025
are scalars; if the input image is a color image, then
Figure FSB00000580858800026
and
Figure FSB00000580858800027
are vectors, and each vector has three components, namely red, green and blue components;
(6)按统计学上t分布上α分位点的定义,检查所述滑动窗口左半部,右半部,上半部和下半部像素点灰度或色彩的平均值,即
Figure FSB00000580858800028
是否落在整个滑动窗口灰度或色彩的平均值在置信水平为1-α下的置信区间内: P { | c &OverBar; j i ( x , y ) - c &OverBar; i ( x , y ) &sigma; i ( x , y ) / n | < z &alpha; / 2 } = 1 - &alpha; - - - ( 7 )
(6) According to the definition of the α quantile point on the t distribution in statistics, check the left half of the sliding window, the right half, the average value of the grayscale or color of the pixel points in the upper half and the lower half, that is
Figure FSB00000580858800028
and Whether it falls within the confidence interval of the mean value of the grayscale or color of the entire sliding window at a confidence level of 1-α: P { | c &OverBar; j i ( x , the y ) - c &OverBar; i ( x , the y ) &sigma; i ( x , the y ) / no | < z &alpha; / 2 } = 1 - &alpha; - - - ( 7 )
如果输入图像是灰度图像,这里上标i=1,因为图像中每个像素点上只有一个灰度值,而当输入图像是彩色图像时,这里上标i表示色彩平均值或色彩方差中对应的红色,绿色或蓝色分量,而
Figure FSB000005808588000211
中下标j则表示是所述滑动窗口左半部,右半部,上半部或下半部,n表示落在所述滑动窗口左半部,右半部,上半部或下半部中像素点的个数,如果所述滑动窗口是正方形,n数值则刚好等于(2r+1)r,根据给定α数值,通过t分布表查询得到zα/2的数值;
If the input image is a grayscale image, the superscript i=1 here, because there is only one grayscale value on each pixel in the image, and when the input image is a color image, the superscript i here represents the color average or color variance corresponding red, green or blue components, while
Figure FSB000005808588000211
The middle subscript j indicates that it is the left half, the right half, the upper half or the lower half of the sliding window, and n indicates that it falls on the left half, the right half, the upper half or the lower half of the sliding window The number of pixels in the middle, if the sliding window is a square, the n value is just equal to (2r+1)r, according to the given α value, the value of z α/2 is obtained by querying the t distribution table;
(7)如果所述滑动窗口左半部,右半部,上半部和下半部像素点灰度或色彩的平均值,即
Figure FSB000005808588000213
中的每个分量,都满足公式(7),则判断所述滑动窗口中心所在的像素点(x,y)落在一个边长为(2r+1)的均匀区域,跳到步骤(9);
(7) If the left half of the sliding window, the right half, the average value of pixel grayscale or color in the upper half and lower half, that is and
Figure FSB000005808588000213
Each component in all satisfies the formula (7), then it is judged that the pixel point (x, y) where the center of the sliding window is located falls in a uniform area with a side length of (2r+1), skip to step (9) ;
(8)如果所述滑动窗口左半部,右半部,上半部和下半部像素点灰度或色彩的平均值,即
Figure FSB000005808588000214
Figure FSB000005808588000215
中的每个分量有一个不满足公式(7),说明当前滑动窗口所在的图像区域不是一个均匀区域,则按照从大到小次序从上述平滑滤波空间尺度可供选择的离散数值集合中选择下一个较小数值作为滑动窗口的半径r,然后回到步骤(4);
(8) If the left half of the sliding window, the right half, the average value of pixel grayscale or color in the upper half and lower half, that is
Figure FSB000005808588000214
and
Figure FSB000005808588000215
If one of the components in α does not satisfy the formula (7), it means that the image area where the current sliding window is located is not a uniform area, then select from the discrete value set of the above smoothing filter spatial scale in order from large to small. A smaller value is used as the radius r of the sliding window, and then returns to step (4);
(9)当前滑动窗口所取的半径r作为滑动窗口中心所在的像素点的平滑滤波空间尺度;(9) The radius r taken by the current sliding window is used as the smoothing filter spatial scale of the pixel at the center of the sliding window; (10)所述滑动窗口按照光栅点阵扫描顺序移到下一个像素点,重复步骤(3)到步骤(9),直至图像中所有像素都确定各自的平滑滤波空间尺度。(10) The sliding window is moved to the next pixel according to the scanning order of the raster lattice, and steps (3) to (9) are repeated until all pixels in the image have determined their respective spatial scales for smoothing and filtering.
CN2010102379141A 2010-07-27 2010-07-27 Self-adaptive confirming method of smooth wave-filtering spatial scale facing to picture edge detection Expired - Fee Related CN101908208B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN2010102379141A CN101908208B (en) 2010-07-27 2010-07-27 Self-adaptive confirming method of smooth wave-filtering spatial scale facing to picture edge detection

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN2010102379141A CN101908208B (en) 2010-07-27 2010-07-27 Self-adaptive confirming method of smooth wave-filtering spatial scale facing to picture edge detection

Publications (2)

Publication Number Publication Date
CN101908208A CN101908208A (en) 2010-12-08
CN101908208B true CN101908208B (en) 2011-11-09

Family

ID=43263661

Family Applications (1)

Application Number Title Priority Date Filing Date
CN2010102379141A Expired - Fee Related CN101908208B (en) 2010-07-27 2010-07-27 Self-adaptive confirming method of smooth wave-filtering spatial scale facing to picture edge detection

Country Status (1)

Country Link
CN (1) CN101908208B (en)

Families Citing this family (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103795943B (en) * 2012-11-01 2017-05-17 富士通株式会社 Image processing apparatus and image processing method
CN105809665B (en) * 2014-12-31 2018-06-22 展讯通信(上海)有限公司 The detection method and device of picture noise level
CN107274363B (en) * 2017-06-02 2020-09-22 北京理工大学 An edge-preserving image filtering method with scale-sensitive properties
CN108389162B (en) * 2018-01-09 2021-06-25 浙江大学 An Image Edge Preserving Filtering Method Based on Adaptive Neighborhood Shape
CN109448010B (en) * 2018-08-31 2020-08-04 浙江理工大学 A method for automatic generation of square continuous patterns based on content features
CN109840901B (en) * 2019-01-09 2020-08-25 武汉精立电子技术有限公司 Quick judgment method for split screen Mura

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6707940B1 (en) * 2000-03-31 2004-03-16 Intel Corporation Method and apparatus for image segmentation
CN101242488A (en) * 2007-02-06 2008-08-13 西北工业大学 A Smoothing Method of Extremum of Digital Image
CN101447077A (en) * 2008-12-18 2009-06-03 浙江大学 Edge detection method of color textile texture image oriented to textile industry
CN101464998A (en) * 2009-01-15 2009-06-24 浙江大学 Non-gauss veins noise smooth filtering method for textile industry

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP0440786A1 (en) * 1989-08-28 1991-08-14 Eastman Kodak Company A computer based digital image noise reduction method based on overlapping planar approximation

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6707940B1 (en) * 2000-03-31 2004-03-16 Intel Corporation Method and apparatus for image segmentation
CN101242488A (en) * 2007-02-06 2008-08-13 西北工业大学 A Smoothing Method of Extremum of Digital Image
CN101447077A (en) * 2008-12-18 2009-06-03 浙江大学 Edge detection method of color textile texture image oriented to textile industry
CN101464998A (en) * 2009-01-15 2009-06-24 浙江大学 Non-gauss veins noise smooth filtering method for textile industry

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
JP平4-501628A 1992.03.19

Also Published As

Publication number Publication date
CN101908208A (en) 2010-12-08

Similar Documents

Publication Publication Date Title
CN101908208B (en) Self-adaptive confirming method of smooth wave-filtering spatial scale facing to picture edge detection
CN104978715B (en) Non-local mean image denoising method based on filtering window and parameter self-adaption
CN105654445B (en) A kind of handset image denoising method based on wavelet transformation edge detection
CN101464998B (en) Non-gauss veins noise smooth filtering method for textile industry
Zontak et al. Separating signal from noise using patch recurrence across scales
KR100872253B1 (en) How to remove noise in an image created by the image sensor
CN114926436B (en) Defect detection method for periodic pattern fabric
US8571312B2 (en) Image interpolation method and apparatus using pattern characteristics of color filter array
CN102156996B (en) Image edge detection method
KR101637488B1 (en) Image interpolation method and apparatus using pattern characteristics of color filter array
DE602005001967T2 (en) Image processing apparatus and method, recording medium and program
CN103347190B (en) Edge-related and color-combined demosaicing and amplifying method
TWI393073B (en) Image denoising method
CN116205823A (en) A Method of Ultrasound Image Denoising Based on Spatial Filtering
CN104966285B (en) A kind of detection method of salient region
CN102103751A (en) Foreground image extraction method and device
CN115330795B (en) Cloth burr defect detection method
KR20140141027A (en) A method and apparatus for detecting a repetitive pattern in image
CN104835127B (en) An Adaptive Smoothing Filtering Method
CN101447077A (en) Edge detection method of color textile texture image oriented to textile industry
CN106373131B (en) Edge-based image salient region detection method
CN101431606A (en) Self-adapting denoising processing method based on edge detection
CN104484865A (en) Method for removing raindrops in video image
CN104680548A (en) Omnidirectional M-type kinetocardiogram detection method
JP4994955B2 (en) Mobile object identification device and mobile object identification program

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
CF01 Termination of patent right due to non-payment of annual fee
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20111109

Termination date: 20180727