A computer-implemented method of processing an input digital image is provided. The method includes forming a representation of the input digital image as wavelet components; and for each of a plurality of facial regions in the input digital image: forming a reconstructed facial region based on a subset of the wavelet components; and identifying pixels associated with facial wrinkles in the facial region as pixels In the reconstructed facial region that differ from corresponding pixels in the facial region by a first threshold amount.