This work estimated the land surface emissivities (LSEs) for MODIS thermal infrared channels 29 (8.4-8.7 μm), 31 (10.78-11.28 μm), and 32 (11.77-12.27 μm) using an improved normalized difference vegetation index (NDVI)-based threshold method. The channel LSEs are expressed as functions of atmospherically corrected reflectance from the MODIS visible and near-infrared channels with wavelengths ranging from 0.4 to 2.2 μm for bare soil. To retain the angular information, the vegetation LSEs were explicitly expressed in the NDVI function. The results exhibited a root mean square error (RMSE) among the estimated LSEs using the improved method, and those calculated using spectral data from Johns Hopkins University (JHU) are below 0.01 for channels 31 and 32. The MODIS land surface temperature/emissivity (LST/E) products, MOD11_L2 with LSE derived via the classification-based method with 1 km resolution and MOD11C1 with LSE retrieved via the day/night LST retrieval method at 0.05°resolution, were used to validate the proposed method. The resultant variances and entropies for the LSEs estimated using the proposed method were larger than those extracted from MOD11_L2, which indicates that the proposed method better described the spectral variation for different land covers. In addition, comparing the estimated LSEs to those from MOD11C1 yielded RMSEs of approximately 0.02 for the three channels; however, more than 70% of pixels exhibited LSE differences within 0.01 for channels 31 and 32, which indicates that the proposed method feasibly depicts LSE variation for different land covers.