To improve the usability of Enhanced Thematic Mapper Plus ETM+ scan line corrector SLC-off data, this article proposes using HJ-1A/1B imagery as auxiliary i.e. reference data to recover the SLC-off ETM+ data. The least-median-of-squares LMedS method is newly proposed to recover missing pixels of Landsat 7 by removing the variant or abnormal digital number values. In particular, for the visible and near-infrared bands, using HJ-1A/1B for recovery has three clear advantages: the same spatial resolution, similar spectral resolution, and approximate temporal resolution. The experiments show that a...