There is disclosed a method of enhancing an image of a scene. The method divides the image into a plurality of frequency zones, and performs a plurality of Richardson and Lucy (RL) iterations of each of the frequency zones to obtain an enhanced image of the scene. The frequency zones and the plurality of RL iterations are determined using a Taguchi orthogonal array computation that minimizes radiometric errors in the enhanced image. The image may be Rayleigh sampled image obtained by a Rayleigh sampling system. A truth scene is not needed to use the method.