1. A method of color interpolation of an image acquired by a digital color sensor generating values of grey levels for each image pixel as a function of a filter applied to the sensor, interpolating values of missing colors of each image pixel for generating distinct values of primary colors or complementary hues for each image pixel, the method comprising: a) selecting a (2k+1)×(2k+1) working window of said image centered around a pixel to be filtered, k being an integer number; b) determining a direction of interpolation of the central pixel of said working window; c) calculating statistical parameters defined as a function of moments of an intensity of pixels of a same primary color or complementary hue of said (2k+1)×(2k+1) window and of at least another working window of smaller size centered around the pixel to be filtered; and d) as a function of values of said statistical parameters of the (2k+1)×(2k+1) working window and of said at least one working window of smaller size, deciding of filtering said central pixel according to an interpolation chosen in a group that includes a homogeneous low-pass interpolation carried out on pixels of the (2k+1)×(2k+1) working window, an interpolation carried out on pixels of the working window of smaller size, a directional filtering interpolation carried out on pixels of the working window of smaller size, and adaptive interpolations determined by said calculated statistical moments.