We propose a novel symmetrization method for solving the transport of intensity equation (TIE) using fast Fourier transforms for situations where the input images may or may not exhibit spatial periodicity. The method is derived from the analysis of intensity conservation law and the internal symmetry of the TIE, and is illustrated for both a computational and an experimental data set.