1. A method of processing a digital image, comprising the steps of: providing digital data indexed to represent positions of an image having S spectral bands for simultaneous output on a display, said digital data being indicative of an intensity value I i (x,y) for each position (x,y) in each i-th spectral band; defining a classification of said image based on dynamic range of said image in each of said S spectral bands; adjusting said intensity value for said each position in each i-th spectral band to generate an adjusted intensity value for said each position in each i-th spectral band in accordance with 9 &Sum; n = 1 N &it; W n &af; ( log &it; &it; I i &af; ( x, y ) - log &af; [ I i &af; ( x, y ) * F n &af; ( x, y ) ] ), i = 1, …, S where S is the number of unique spectral bands included in said digital data and, for each n, W n is a weighting factor and F n (x,y) is a unique surround function applied to said each position (x,y) and N is the total number of unique surround functions; and filtering said adjusted intensity value for said each position of said image in each of said S spectral bands using a filter function based on said classification of said image wherein a filtered intensity value R i (x,y) is defined.