Due to the nonlinear property of large uniaxial anisotropic scatterers, many iterative optimization methods have a high risk of being trapped in local minima. In this paper, a frequency-hopping subspace-based optimization method (SOM) is proposed to reconstruct the relative permittivity distribution of 2-D large uniaxial anisotropic scatterers with transverse electrical (TE) illumination. This hybrid method utilizes the results obtained at lower frequency to provide good initial guesses for higher frequency reconstruction, which reduces the occurrence of local minima for the inversion at the higher frequency. For lower frequency, it can only obtain coarse resolution image although it is unable to show the details of the scatterers. However, this coarse image provides a priori information for the reconstruction at higher frequencies to get finer resolution. Numerical examples demonstrate that the proposed hybrid method can effectively rebuild large uniaxial anisotropic scatterers (six wavelengths) with higher stability compared with conventional SOM that uses only single-frequency data.