Low-curvature image simplifiers: Global regularity of smooth solutions and Laplacian limiting schemes

We consider a class of fourth-order nonlinear diffusion equations motivated by Tumblin and Turk's "low-curvature image simplifiers" for image denoising and segmentation. The PDE for the image intensity it is of the form u(t) = -del . (g(Deltau)delDeltau) + lambda(f - u), where g(s) = k(2)/(k(2) + s(2)) is a "curvature" threshold and lambda denotes a fidelity-matching parameter. We derive a priori bounds for Deltau that allow us to prove global regularity of smooth solutions in one space dimension, and a geometric constraint for finite-time singularities from smooth initial data ill two space dimensions. This is in sharp contrast to the second-order Perona-Malik equation (an ill-posed problem), on which the original LCIS method is modeled. The estimates also allow us to design a finite difference scheme that satisfies discrete versions of the estimates, in particular, a priori bounds oil the smoothness estimator in both one and two space dimensions. We present computational results that show the effectiveness of such algorithms. Our results are connected to recent results for fourth-order lubrication-type equations and the design of positivity-preserving schemes for such equations. This connection also has relevance for other related fourth-order imaging equations. (C) 2004 Wiley Periodicals, Inc.