An image processing method for saturation enhancement comprising the steps of: setting up a function x0=f 1 (x i)=α·x i +β, where x i ∈{r i, g i, b i} is a component value in the range [0.. 1] of a point A i (r i, g i, b i) in normalized RGB color space, x 0 ∈{r 0, g 0, b 0} is a component value of a point Ao(ro, go, bo), and α and β are scalar parameters which are constant for each point A i; and setting up a function S(r, g, b)=1-3·min[r, g, b]/(r+g+b) defining the saturation of a point A(r, g, b), characterized in comprising the steps of: based on S(r, g, b), setting up a function S 0 =f 2 (r i, g i, b i)=1-3·(α·min[r i, g i, b i]+β)/(r i +g i +b i) defining the saturation of point A 0 (r 0, g 0, b 0) expressed in dependency of the component values of point A i (r i, g i, b i) under the condition that the brightness l i =r i +g i +b i of point A i is identical to the brightness lo=ro+go+bo of point Ao; setting the saturation S 0 to a value of 1, which is a maximum while a value of 0 is a minimum, and for each color point A i subject to saturation enhancement solving for α and β using functions f 1 (x i), f 2 (r i, g i, b i) and the condition (r 0 +g 0 +b 0)=(r i +g i +b i) to obtain the corresponding point Ao(ro, go, bo); obtaining the desired output saturation S 0 desired for input point A i using a predefined saturation enhancement function S 0 desired =f(S i); performing saturation enhancement on point A i using the desired output saturation according to the functions defined above in the RGB color space if max(ro, go, b 0)≤1, otherwise converting the RGB color components r i, g i, b i of point A i into the CMY color space using the conversion c i =1-r i, m i =1-g i, y i =1-b i and performing saturation enhancement using the same value S 0 desired and the same functions f 1 () and f 2 () as for the r, g, b coordinates on the c, m, y coordinates. setting up a function x0=f 1 (x i)=α·x i +β, where x i ∈{r i, g i, b i} is a component value in the range [0.. 1] of a point A i (r i, g i, b i) in normalized RGB color space, x 0 ∈{r 0, g 0, b 0} is a component value of a point Ao(ro, go, bo), and α and β are scalar parameters which are constant for each point A i; and setting up a function S(r, g, b)=1-3·min[r, g, b]/(r+g+b) defining the saturation of a point A(r, g, b), characterized in comprising the steps of: based on S(r, g, b), setting up a function S 0 =f 2 (r i, g i, b i)=1-3·(α·min[r i, g i, b i]+β)/(r i +g i +b i) defining the saturation of point A 0 (r 0, g 0, b 0) expressed in dependency of the component values of point A i (r i, g i, b i) under the condition that the brightness l i =r i +g i +b i of point A i is identical to the brightness lo=ro+go+bo of point Ao; setting the saturation S 0 to a value of 1, which is a maximum while a value of 0 is a minimum, and for each color point A i subject to saturation enhancement solving for α and β using functions f 1 (x i), f 2 (r i, g i, b i) and the condition (r 0 +g 0 +b 0)=(r i +g i +b i) to obtain the corresponding point Ao(ro, go, bo); obtaining the desired output saturation S 0 desired for input point A i using a predefined saturation enhancement function S 0 desired =f(S i); performing saturation enhancement on point A i using the desired output saturation according to the functions defined above in the RGB color space if max(ro, go, b 0)≤1, otherwise converting the RGB color components r i, g i, b i of point A i into the CMY color space using the conversion c i =1-r i, m i =1-g i, y i =1-b i and performing saturation enhancement using the same value S 0 desired and the same functions f 1 () and f 2 () as for the r, g, b coordinates on the c, m, y coordinates. based on S(r, g, b), setting up a function S 0 =f 2 (r i, g i, b i)=1-3·(α·min[r i, g i, b i]+β)/(r i +g i +b i) defining the saturation of point A 0 (r 0, g 0, b 0) expressed in dependency of the component values of point A i (r i, g i, b i) under the condition that the brightness l i =r i +g i +b i of point A i is identical to the brightness lo=ro+go+bo of point Ao; setting the saturation S 0 to a value of 1, which is a maximum while a value of 0 is a minimum, and for each color point A i subject to saturation enhancement solving for α and β using functions f 1 (x i), f 2 (r i, g i, b i) and the condition (r 0 +g 0 +b 0)=(r i +g i +b i) to obtain the corresponding point Ao(ro, go, bo); obtaining the desired output saturation S 0 desired for input point A i using a predefined saturation enhancement function S 0 desired =f(S i); performing saturation enhancement on point A i using the desired output saturation according to the functions defined above in the RGB color space if max(ro, go, b 0)≤1, otherwise converting the RGB color components r i, g i, b i of point A i into the CMY color space using the conversion c i =1-r i, m i =1-g i, y i =1-b i and performing saturation enhancement using the same value S 0 desired and the same functions f 1 () and f 2 () as for the r, g, b coordinates on the c, m, y coordinates.