A method for providing an expanded dynamic range image from a digital image defined by pixels associated with colors represented in a color space separating luminance from chrominance, comprising: determining a plurality of luminance zones in which luminance values representing said colors of pixels are distributed, for each of said plurality of luminance zones, defining an expansion exponent quadratic function of low-pass filtered luminance determining an expansion exponent value E(p) for each pixel (p) of said luminance zone, inverse tone mapping luminance value of each pixel by applying an expansion exponent value E(p) to said luminance value, wherein said expansion exponent value E(p) is obtained by applying to a low-pass filtered luminance value of said luminance value the expansion exponent quadratic function of the luminance zone to which said luminance value belongs, and providing the expanded dynamic range image based on the inverse tone mapped luminance values, wherein said plurality of luminance zones comprises at least three luminance zones, a low zone for luminances of low values, a high zone for luminances of high values, and a middle zone for luminances comprised between luminances of low values and luminances of high values, wherein said expansion exponent quadratic functions are defined by one set of three parameters a', b', c' and by the position of at least five expansion exponent control points (C1, C2, C3, C4, C5) such that: - said expansion exponent quadratic function determining expansion exponent values in middle zone passes through the second, third, fourth of these five expansion exponent control points, - said expansion exponent quadratic function determining expansion exponent values in low zone passes through the first and second of these five expansion exponent control points, and - said expansion exponent quadratic function determining expansion exponent values in high zone passes through the fourth and fifth of these five expansion exponent control points. determining a plurality of luminance zones in which luminance values representing said colors of pixels are distributed, for each of said plurality of luminance zones, defining an expansion exponent quadratic function of low-pass filtered luminance determining an expansion exponent value E(p) for each pixel (p) of said luminance zone, inverse tone mapping luminance value of each pixel by applying an expansion exponent value E(p) to said luminance value, wherein said expansion exponent value E(p) is obtained by applying to a low-pass filtered luminance value of said luminance value the expansion exponent quadratic function of the luminance zone to which said luminance value belongs, and providing the expanded dynamic range image based on the inverse tone mapped luminance values, wherein said plurality of luminance zones comprises at least three luminance zones, a low zone for luminances of low values, a high zone for luminances of high values, and a middle zone for luminances comprised between luminances of low values and luminances of high values, wherein said expansion exponent quadratic functions are defined by one set of three parameters a', b', c' and by the position of at least five expansion exponent control points (C1, C2, C3, C4, C5) such that: - said expansion exponent quadratic function determining expansion exponent values in middle zone passes through the second, third, fourth of these five expansion exponent control points, - said expansion exponent quadratic function determining expansion exponent values in low zone passes through the first and second of these five expansion exponent control points, and - said expansion exponent quadratic function determining expansion exponent values in high zone passes through the fourth and fifth of these five expansion exponent control points. - said expansion exponent quadratic function determining expansion exponent values in middle zone passes through the second, third, fourth of these five expansion exponent control points, - said expansion exponent quadratic function determining expansion exponent values in low zone passes through the first and second of these five expansion exponent control points, and - said expansion exponent quadratic function determining expansion exponent values in high zone passes through the fourth and fifth of these five expansion exponent control points.