A method of deblurring images e.g. from an MRI scan, comprises calculating the distribution of the slope value of the intensities between the pixels in the image, and using a Levy distribution of power factor k of 1 (Cauchy distribution) up to but not including (2) to determine the correction term to be applied.