Neighborhood filters and PDE’s

Denoising images can be achieved by a spatial averaging of nearby pixels. However, although this method removes noise it creates blur. Hence, neighborhood filters are usually preferred. These filters perform an average of neighboring pixels, but only under the condition that their grey level is close enough to the one of the pixel in restoration. This very popular method unfortunately creates shocks and staircasing effects. In this paper, we perform an asymptotic analysis of neighborhood filters as the size of the neighborhood shrinks to zero. We prove that these filters are asymptotically equivalent to the Perona–Malik equation, one of the first nonlinear PDE’s proposed for image restoration. As a solution, we propose an extremely simple variant of the neighborhood filter using a linear regression instead of an average. By analyzing its subjacent PDE, we prove that this variant does not create shocks: it is actually related to the mean curvature motion. We extend the study to more general local polynomial estimates of the image in a grey level neighborhood and introduce two new fourth order evolution equations.This work has been partially financed by the Centre National d’Etudes Spatiales (CNES), the Office of Naval Research under grant N00014-97-1-0839, the Ministerio de Ciencia y Tecnologia under grant MTM2005-08567. During this work, the first author had a fellowship of the Govern de les Illes Balears for the realization of his PhD.