A computer-implemented method for enhancing the contrast of an electronic representation of an original image represented by an array of pixel values by processing said image, said processing comprising the steps of a) decomposing said digital image into a set of detail images at multiple resolution levels, said detail images being the difference between an image at a certain scale and at a smaller scale, and a residual image at a resolution level lower than said multiple resolution levels, b) processing at least one pixel of said detail images, c) computing a processed image by applying a reconstruction algorithm to the residual image and the processed detail images, said reconstruction algorithm being such that if it were applied to the residual image and the detail images without processing, then said digital image would be obtained, characterized in that said processing comprises the steps of: d) performing for said at least one pixel of said detail images, a mathematical operation to calculate a statistical measure var k ( i, j) on the translation difference image pixel value in a neighbourhood of an approximation image pixel, wherein an approximation image pixel represents the original grey value image at the same or a larger scale, or at the same or a lower resolution and wherein a translation difference image pixel value is defined as the result of a difference between an approximation image pixel g l ( ri, rj) and another approximation image pixel g l ( ri + m, rj + n) expressed by g l ri rj − g l ri + m, rj + n e) modifying the value of said detail image pixel value d k ( i, j) to obtain a modified detail image pixel value d k out ( i, j) by applying a function f k of said statistical measure var k ( i, j) and said image detail pixel value d k ( i, j) such that: d k out i j = a k i j d k i j and a k i j = f k var k i j wherein a k ( i, j) is an amplification image pixel value. and f k is a non-linear function with a steeper slope in the region of argument values that correspond to low contrast than in other regions. a) decomposing said digital image into a set of detail images at multiple resolution levels, said detail images being the difference between an image at a certain scale and at a smaller scale, and a residual image at a resolution level lower than said multiple resolution levels, b) processing at least one pixel of said detail images, c) computing a processed image by applying a reconstruction algorithm to the residual image and the processed detail images, said reconstruction algorithm being such that if it were applied to the residual image and the detail images without processing, then said digital image would be obtained, d) performing for said at least one pixel of said detail images, a mathematical operation to calculate a statistical measure var k ( i, j) on the translation difference image pixel value in a neighbourhood of an approximation image pixel, wherein an approximation image pixel represents the original grey value image at the same or a larger scale, or at the same or a lower resolution and wherein a translation difference image pixel value is defined as the result of a difference between an approximation image pixel g l ( ri, rj) and another approximation image pixel g l ( ri + m, rj + n) expressed by g l ri rj − g l ri + m, rj + n e) modifying the value of said detail image pixel value d k ( i, j) to obtain a modified detail image pixel value d k out ( i, j) by applying a function f k of said statistical measure var k ( i, j) and said image detail pixel value d k ( i, j) such that: d k out i j = a k i j d k i j and a k i j = f k var k i j wherein a k ( i, j) is an amplification image pixel value. and f k is a non-linear function with a steeper slope in the region of argument values that correspond to low contrast than in other regions.