A method for beam shaping based on fitting the power moments of the final beam intensity distribution and independent of the optical system particularities is suggested. It is shown how one can analytically calculate any moment of the final phase space distribution using the moments of the initial distribution and the optical system transfer map. Numerical tests carried out for a final focus system have demonstrated the usefulness of the approach developed here.