We investigate the ground state properties of anti-ferromagnetic spin-1 Bose gases in one dimensional harmonic potential from the weak repulsion regime to the strong repulsion regime. By diagonalizing the Hamiltonian in the Hilbert space composed of the lowest eigenstates of single particle and spin components, the ground state wavefunction and therefore the density distributions, magnetization distribution, one body density matrix, and momentum distribution for each components are obtained. It is shown that the spinor Bose gases of different magnetization exhibit the same total density profiles in the full interaction regime, which evolve from the single peak structure embodying the properties of Bose gases to the fermionized shell structure of spin-polarized fermions. But each components display different density profiles, and magnetic domains emerge in the strong interaction limit for $M=0.25$. In the strong interaction limit, one body density matrix and the momentum distributions exhibit the same behaviours as those of spin-polarized fermions. The fermionization of momentum distribution takes place, in contrast to the $delta$-function-like distribution of single component Bose gases in the full interaction region.