We explore the rich and unique magnetic quantization of Si-doped graphene defect systems with various concentrations and configurations using the generalized tight-binding model. This model takes into account simultaneously the non-uniform bond lengths, site energies and hopping integrals, and a uniform perpendicular magnetic field (B z ẑ). The magnetic quantized Landau levels (LLs) are classified into four different kinds based on the probability distributions and oscillation modes. The main characteristics of LLs are clearly reflected in the magneto-optical selection rules which cover the dominating ∆ n = |n v − n c | = 0, the coexistent ∆ n = 0 & ∆ n = 1, and the specific ∆ n = 1. These rules for inter-LLs excitations come from the non-equivalence or equivalence of the A i and B i sublattices in a supercell.