Bilateral filter

A bilateral filter is a non-linear, edge-preserving and noise-reducing smoothing filter for images. The intensity value at each pixel in an image is replaced by a weighted average of intensity values from nearby pixels. This weight can be based on a Gaussian distribution. Crucially, the weights depend not only on Euclidean distance of pixels, but also on the radiometric differences (e.g. range differences, such as color intensity, depth distance, etc.). This preserves sharp edges by systematically looping through each pixel and adjusting weights to the adjacent pixels accordingly.

The bilateral filter is defined as


I^\text{filtered}(x) = \frac{1}{W_p} \sum_{x_i \in \Omega} I(x_i)f_r(\|I(x_i)-I(x)\|)g_s(\|x_i-x\|),

where the normalization term


W_p = \sum_{x_i \in \Omega}{f_r(\|I(x_i)-I(x)\|)g_s(\|x_i-x\|)}

ensures that the filter preserves image energy and

As mentioned above, the weight W_p is assigned using the spatial closeness and the intensity difference.[1] Consider a pixel located at (i, j) which needs to be denoised in image using its neighbouring pixels and one of its neighbouring pixels is located at (k, l). Then, the weight assigned for pixel (k, l) to denoise the pixel (i, j) is given by: 
w(i, j, k, l)= e^{(-\frac{(i-k)^2+ (j-l)^2}{2 \sigma_d^2}- \frac{\|I(i, j)- I(k, l)\|^2}{2 \sigma_r^2})}

where σd and σr are smoothing parameters and I(i, j) and I(k, l) are the intensity of pixels (i, j) and  (k, l) respectively. After calculating the weights, normalize them. 
I_D(i, j)=\frac{\sum_{k, l}{I(k, l) * w(i, j, k, l)} } { \sum_{k, l}{w(i, j, k, l)} }

where I_D is the denoised intensity of pixel (i, j).

Parameters

Limitations

The bilateral filter in its direct form can introduce several types of image artifacts:

There exist several extensions to the filter that deal with these artifacts. Alternative filters, like the guided filter , have also been proposed as an efficient alternative without these limitations.

Implementations

Adobe Photoshop implements a bilateral filter in its surface blur tool. GIMP implements a bilateral filter in its Filters-->Blur tools; and it is called Selective Gaussian Blur'.

Related models

The Bilateral filter was shown to be an application of the short time kernel of the Beltrami flow [2] see also [3]

See also

External links

References

  1. Carlo Tomasi and Roberto Manduchi, “Bilateral filtering for gray and color images,” in Computer Vision, 1998. Sixth International Conference on . IEEE, 1998, pp. 839– 846.
  2. R. Kimmel, R. Malladi, and N. Sochen. Images as Embedded Maps and Minimal Surfaces: Movies, Color, Texture, and Volumetric Medical Images. International Journal of Computer Vision, 39(2):111-129, Sept. 2000. some color results http://www.cs.technion.ac.il/~ron/PAPERS/KimMalSoc_IJCV2000.pdf
  3. N. Sochen, R. Kimmel, and A.M. Bruckstein. Diffusions and confusions in signal and image processing, Journal of Mathematical Imaging and Vision, 14(3):195-209, 2001.http://www.cs.technion.ac.il/~ron/PAPERS/SocKimBru_JMIV2001.pdf