New moment relationships between the function of wave front aberrations of an optical system and numerical characteristics (moments) of intensity distribution in the image are obtained. These are used as a basis for the method of wave front reconstruction from images of a source in several parallel planes.