An image color correction method, wherein a corrected image is computed from original color components in an input image, the method comprising the steps of - computing color components of the corrected image corresponding to an interpolation between a first and a second corrected color components vectors, wherein - the first corrected color component vector corresponds at least partly to scaling all color components of original color component vectors from the input image by factors that reduce a deviation of a characteristic color component vector for the input image from a target vector, and wherein - the second corrected color component vector corresponds to at least partial replacement of at least one of the color components of the original color component vectors from the input image by values computed from the other color components of the original color component vectors from the input image; - computing (42, 43) an interpolation coefficient from a ratio of an average value of the at least one of the color components of the original color component vectors of the input image or one or more reference images and an average obtained from one or more of the other color components of the original color component vectors from the input image or the one or more reference images; - controlling a position of the interpolation between the first and the second corrected color component vectors by said interpolation coefficient, the position being moved increasingly toward the second corrected color component vector with decreasing values of said ratio. - computing color components of the corrected image corresponding to an interpolation between a first and a second corrected color components vectors, wherein - the first corrected color component vector corresponds at least partly to scaling all color components of original color component vectors from the input image by factors that reduce a deviation of a characteristic color component vector for the input image from a target vector, and wherein - the second corrected color component vector corresponds to at least partial replacement of at least one of the color components of the original color component vectors from the input image by values computed from the other color components of the original color component vectors from the input image; - the first corrected color component vector corresponds at least partly to scaling all color components of original color component vectors from the input image by factors that reduce a deviation of a characteristic color component vector for the input image from a target vector, and wherein - the second corrected color component vector corresponds to at least partial replacement of at least one of the color components of the original color component vectors from the input image by values computed from the other color components of the original color component vectors from the input image; - computing (42, 43) an interpolation coefficient from a ratio of an average value of the at least one of the color components of the original color component vectors of the input image or one or more reference images and an average obtained from one or more of the other color components of the original color component vectors from the input image or the one or more reference images; - controlling a position of the interpolation between the first and the second corrected color component vectors by said interpolation coefficient, the position being moved increasingly toward the second corrected color component vector with decreasing values of said ratio.