A class of fourth-order partial differential equations (PDEs) are proposed
to optimize the trade-off between noise removal and edge preservation. The
time evolution of these PDEs seeks to minimize a cost functional which is a
n increasing function of the absolute value of the Laplacian of the image i
ntensity function. Since the Laplacian of an image at a pixel is zero if th
e image is planar in its neighborhood, these PDEs attempt to remove noise a
nd preserve edges by approximating an observed image with a piecewise plana
r image. Piecewise planar images look more natural than step images which a
nisotropic diffusion (second order PDEs) uses to approximate an observed im
age. So the proposed PDEs are able to avoid the blocky effects widely seen
in images processed by anisotropic diffusion, while achieving the degree of
noise removal and edge preservation comparable to anisotropic diffusion. A
lthough both approaches seem to be comparable in removing speckles in the o
bserved images, speckles are more visible in images processed by the propos
ed PDEs, because piecewise planar images are less likely to mask speckles t
han step images and anisotropic diffusion tends to generate multiple false
edges. Speckles can be easily removed by simple algorithms such as the one
presented in this paper.