A method for Magnetic Resonance Imaging to depict a 3-dimensional object by an image having pixels representing volume elements of the object, comprising: immobilising the object and acquiring a reference image (S R,p) of the object with a receiver coil at a first echo time (TE R) immediately following an excitation by a transmitter coil, wherein said reference image (S R,p) is complex-valued, representing each volume element by a pixel with a reference magnitude value (m R,p) and a reference phase value (θ R,p); keeping the object immobilised and acquiring a target image (S T,p) of the object with said receiver coil at a pre-selected second echo time (TE T), longer than said first echo time (TE R), following the same or another excitation by said transmitter coil, wherein said target image (S T,p) is complex-valued, representing each volume element by a pixel with a target magnitude value (m T,p) and a target phase value (θ T,p); subtracting, pixel by pixel, the reference phase value (θ R,p) from the target phase value (θ T,p) to obtain a corrected phase value (θ K,p) for each pixel; and obtaining said image (S K,p) from said target magnitude values (m T,p) and said corrected phase values (θ K,p), characterised in that the aforementioned steps (4, 7, 10) are applied to each of a plurality of receiver coils arranged around said immobilised object to obtain a respective plurality of said images (S K,p), followed by the step of calculating a combined phase image (Θ S), pixel by pixel, according to θ S = ∠ ∑ p m T, p ⋅ e − i θ T, p − θ R, p with ∠(·) denoting the four-quadrant tangent inverse operator, m T,p being the target magnitude value of a pixel, representing a volume element of the object, in a target image (S T,p) acquired with the p th receiver coil, θ T,p being the target phase value of said pixel, θ R,p being the reference phase value of a pixel, representing said volume element, in a reference image (S R,p) acquired with the p th receiver coil, and θ S being the phase value of a pixel, representing said volume element, in the combined phase image (Θ S). immobilising the object and acquiring a reference image (S R,p) of the object with a receiver coil at a first echo time (TE R) immediately following an excitation by a transmitter coil, wherein said reference image (S R,p) is complex-valued, representing each volume element by a pixel with a reference magnitude value (m R,p) and a reference phase value (θ R,p); keeping the object immobilised and acquiring a target image (S T,p) of the object with said receiver coil at a pre-selected second echo time (TE T), longer than said first echo time (TE R), following the same or another excitation by said transmitter coil, wherein said target image (S T,p) is complex-valued, representing each volume element by a pixel with a target magnitude value (m T,p) and a target phase value (θ T,p); subtracting, pixel by pixel, the reference phase value (θ R,p) from the target phase value (θ T,p) to obtain a corrected phase value (θ K,p) for each pixel; and obtaining said image (S K,p) from said target magnitude values (m T,p) and said corrected phase values (θ K,p), characterised in that the aforementioned steps (4, 7, 10) are applied to each of a plurality of receiver coils arranged around said immobilised object to obtain a respective plurality of said images (S K,p), followed by the step of calculating a combined phase image (Θ S), pixel by pixel, according to θ S = ∠ ∑ p m T, p ⋅ e − i θ T, p − θ R, p with ∠(·) denoting the four-quadrant tangent inverse operator, m T,p being the target magnitude value of a pixel, representing a volume element of the object, in a target image (S T,p) acquired with the p th receiver coil, θ T,p being the target phase value of said pixel, θ R,p being the reference phase value of a pixel, representing said volume element, in a reference image (S R,p) acquired with the p th receiver coil, and θ S being the phase value of a pixel, representing said volume element, in the combined phase image (Θ S). ∠(·) denoting the four-quadrant tangent inverse operator, m T,p being the target magnitude value of a pixel, representing a volume element of the object, in a target image (S T,p) acquired with the p th receiver coil, θ T,p being the target phase value of said pixel, θ R,p being the reference phase value of a pixel, representing said volume element, in a reference image (S R,p) acquired with the p th receiver coil, and θ S being the phase value of a pixel, representing said volume element, in the combined phase image (Θ S).