A method of processing a pixellated image to retrieve a phase distribution representative of the image. The phase distribution representative of the image is in the Fourier domain. The method includes padding the image pixels with padding, or non-image, pixels to increase the total number of pixels in the pixellated image. The method further includes processing the padding or non-image pixels differently to the image pixels in each iteration of the iterative method.