Variance using median filter

Hi,

In calculating variance using definition and implementing median filter why does in the code

variance = medfilt2(double(I_noise).^2, [N N], 'symmetric') - (medfilt2(double(I_noise), [N N],'symmetric')).^2;

produce zero in all points of an image?

How do I correctly calculate the variance of a median filtered image?

Answer 1

[Jon]

In your expression:

variance = medfilt2(double(I_noise).^2, [N N], 'symmetric') - (medfilt2(double(I_noise), [N N],'symmetric')).^2;

both of the terms on the right hand side of the expression are identical. Subtracting one from the other will obviously give zero.

This is the same for example as the expression:

x = 2 - 2

which of course gives an answer of x = 0;

Also I think you are getting confused about median and average.