A method for identifying skin region in an image includes: selecting a pixel in the target image; calculating a first probability of the pixel in a first color space, the first probability being a probability that the pixel represents skin in the first color space; calculating a second probability of the pixel in a second color space distinct from the first color space, the second probability being a probability that the pixel represents skin in the second color space; determining a combined probability that the pixel in the target image represents skin based on the first probability and the second probability; expanding the pixel into a skin region in the target image by repeating said steps of determining a combined probability that another pixel adjacent to the pixel represents skin until the skin region reaches a predefined size threshold; and performing a predefined beautifying process to the pixels.