A method of producing a magnetic resonance (MR) image of a region of interest is provided. The method includes the steps of: acquiring an initial MR image of the region of interest, the initial MR image mapping values of an MR-sensitive, physical property at positions over the region; determining a corresponding map of the estimated uncertainties in the values of the MR-sensitive, physical property over the region; and calculating a weighted MR image of the region, the weighted MR image mapping values of a function which combines, at each position of the initial image, the respective value of the MR-sensitive, physical property and the respective estimated uncertainty, the function applying a higher weighting to positions with relatively low estimated uncertainties than to positions with relatively high estimated uncertainties.