Monday, October 23, 2017

image processing - Are 2nd Order Edge Detectors More Susceptible to Noise?


I am reading about edge detection and I read that 2nd order detectors are more susceptible to noise. Is there a mathematical proof to this ?



Answer



Suppose that the noise is a random vector X with normal zero-mean components of variance σi, mutually independent, then for the linear combination (the gi being for instance coefficients of a FIR filter): Y=igiXi, the variance of Y will be: V(Y)=ig2iσ2i, which boils down to for a classical stationnary Gaussian noise.


Now start from g^1=[1,-1], and assume you built higher order derivative by self-convolution, for instance g^2=g^1 \ast g^1 = [1 ,-2 ,1], g^3 = g^1 \ast g^1 \ast g^1 =[1,-3,3,1] , etc.


You easily get that the coefficients of g^n are g^n_k=(-1)^k \binom{n}{k}\,.


Hence \|g^n\|_2^2 = \sum_k \binom{n}{k} ^2 = \binom{2n}{n} (a famous identity). For n=1, the amplification is 2, for n=2 it will be 6, and 20 for the third derivative.


So the noise amplification roughly grows with the derivative order (based on Stirling's formula) as \frac{4^n}{\sqrt{\pi n}}\,.



So all in all, noise power tends to explode with derivatives.


This does not happen at order 0, if you consider that traditional low-pass filters having positive coefficients have finite sum, often equal to 1. Hence, for those \|g\|_2 \le \|g\|_1=1. For instance, if one takes standard smoothing "discrete Gaussian" filters of order m with coefficients b^m_k=\frac{ \binom{m}{k}}{2^m}\,, then their \ell_2 squared norm is: \frac{\binom{2m}{m}}{4^m} and one sees that the 4^m term cancels in the asymptotics, and you get that, when m grows: \|h^m\|_2^2 \approx\frac{1}{\sqrt{\pi m}}\,, and thus tend to reduce the noise power.


Finally, this is not the whole story, since we have not talked about edges, which traditionally are wide-band, and especially hard to detect in high frequencies. A study in the Fourier domains may lead to others aspects, that exceed my time today, but is addressed in @Olli Niemitalo's answer.


No comments:

Post a Comment

periodic trends - Comparing radii in lithium, beryllium, magnesium, aluminium and sodium ions

Apparently the of last four, \ce{Mg^2+} is closest in radius to \ce{Li+}. Is this true, and if so, why would a whole larger shell ($\ce{...