1. A method for processing an input image, comprising: employing a processor for providing digital input image pixels indexed to represent positions on a display, each pixel being indicative of an intensity value for each position; and determining pixel illumination correction for a pixel based on control parameters computed based on content and the local brightness of said input image, and adjusting the intensity value for the pixel based on the illumination correction to generate an enhanced value, wherein determining pixel illumination correction further includes: estimating the local illumination of the input image by performing a weighted low pass filter (WLPF) operation on the image pixels; separating the filtered image into illumination and reflectance images; adjusting the pixel intensity value of the illumination images using a non-linear mapping function as illumination correction, based on the estimated local illumination at each pixel location, wherein adjusting said intensity value using said illumination correction further includes: separating said input image Y IN (x,y) into illumination image L(x,y) and reflectance image R(x,y) as Y IN (x,y)=L(x,y)·R(x,y), correcting the illumination image L(x,y) by calculating Y IN ′ ⁡ ( x , y ) = L max · ( L ⁡ ( x , y ) L max ) 1 ⁢ / ⁢ γ · R ⁡ ( x , y ) , where Y′ IN (x,y) is the corrected illumination image, L max is the maximum value of the input image, and γ is the correction value; and combining the adjusted illumination image with the reflectance image to generate an output image, wherein said illumination image is estimated by using a weighted low pass filter of the size m×n as: L ⁡ ( x , y ) = ∑ i = 1 m ⁢ ∑ j = 1 n ⁢ g ⁡ (  Y IN ⁡ ( i , j ) - Y IN ⁡ ( x , y )  ) · Y IN ⁡ ( i , j ) ∑ i = 1 m ⁢ ∑ j = 1 n ⁢ g ⁡ (  Y IN ⁡ ( i , j ) - Y IN ⁡ ( x , y )  ) , wherein the weighting function g(X) is determined as: g ⁡ ( X ) = { 0 , if ⁢ ⁢ X > σ · τ 2 , - 1 ( τ 2 - τ 1 ) · σ · ( X - σ · τ 2 ) , else ⁢ ⁢ if ⁢ ⁢ X > σ · τ 1 , 1 , otherwise wherein X is an input intensity difference value, σ is an predefined variance, and τ 1 and τ 2 are first and second thresholds, respectively.