A method of retrieving phase information from input intensity information, representative of a target image, in which a Fourier transform is performed on data and the result used in forming a phase estimate, the phase estimate being inverse Fourier transformed, thereby producing magnitude and phase replay, and wherein not only is the phase reply component but also data derived from the magnitude replay component, iteratively fed back.