A second order partial differential operator is applied to an image function. To this end we consider both the Laplacian and a more general elliptic operator. By using a multigrid operator known from the so-called approximation property, we derive a multiresolution decomposition of the image without blurring of edges at coarser levels. We investigate both a linear and a nonlinear variant and compare to some established methods