A method for processing an image using a bilateral filter. The bilateral filter is reformulated at each pixel location in the image into a sum of the original signal value of a central pixel at said pixel location and a bilateral correction term which is a function of local signal differences between the central pixel and its neighbors. The bilateral correction term is calculated using a computationally efficient approximation.