Near-infrared (NIR) phosphors based on Cr3+ doped garnets present great potential in the next generation of NIR light sources. Nevertheless, the huge searching space for the garnet composition makes the rapid discovery of NIR phosphors with high performance remain a great challenge for the scientific community. Herein, a generalizable machine learning (ML) strategy is designed to accelerate the exploration of innovative NIR phosphors via establishing the relationship between key parameters and emission peak wavelength (EPW). We propose a semi-supervised co-training model based on kernel ridge regression (KRR) and support vector regression (SVR), which successfully establishes an expanded dataset with unlabeled dataset (previously unidentified garnets), addressing the overfitting issue resulted from a small dataset and greatly improving the model generalization capability. The model is then interpreted to extract valuable insights into the contribution originated from different features. And a new type NIR luminescent material of Lu3Y2Ga3O12: Cr3+ (EPW~750 nm) is efficiently screened, which demonstrates a high internal (external) quantum efficiency of 97.1% (38.8%) and good thermal stability, particularly exhibiting promising application in the NIR phosphor-converted LEDs (pc-LED). These results suggest the strategy proposed in this work could provide new viewpoint and direction for developing NIR luminescence materials.