The present invention is a method of processing a digital image that is initially represented by digital data indexed to represent positions on a display. The digital data is indicative of an intensity value Ii(x,y) for each position (x,y) in each i-th spectral band. The intensity value for each position in each i-th spectral band is adjusted to generate an adjusted intensity value for each position in each i-th spectral band in accordance with ∑n=1NWn(log Ii(x,y)-log [Ii(x,y)*Fn(x,y)]),i=1,… ,Swhere Wn is a weighting factor, “*” is the convolution operator and S is the total number of unique spectral bands. For each n, the function Fn(x,y) is a unique surround function applied to each position (x,y) and N is the total number of unique surround functions. Each unique surround function is scaled to improve some aspect of the digital image, e.g., dynamic range compression, color constancy, and lightness rendition. The adjusted intensity value for each position in each i-th spectral band of the image is then filtered with a filter function to generate a filtered intensity value Ri(x,y). To prevent graying of white zones in the image, the maximum of the original intensity value Ii(x,y) and filtered intensity value Ri(x,y) is selected for display.