Mastering the evolution of urban land cover is important for urban management and planning. In this paper, a method for analyzing land cover evolution within urban built-up areas based on nighttime light data and Landsat data is proposed. The method solves the problem of inaccurate descriptions of urban built-up area boundaries from the use of single-source diurnal or nocturnal remote sensing data and was able to achieve an effective analysis of land cover evolution within built-up areas. Four main procedures are involved: (1) The neighborhood extremum method and maximum likelihood method are used to extract nighttime light data and the urban built-up area boundaries from the Landsat data, respectively; (2) multisource urban boundaries are obtained using boundary pixel fusion of the nighttime light data and Landsat urban built-up area boundaries; (3) the maximum likelihood method is used to classify Landsat data within multisource urban boundaries into land cover classes, such as impervious surface, vegetation and water, and to calculate landscape indexes, such as overall landscape trends, degree of fragmentation and degree of aggregation; (4) the changes in the multisource urban boundaries and landscape indexes were obtained using the abovementioned methods, which were supported by multitemporal nighttime light data and Landsat data, to model the urban land cover evolution. Using the cities of Shenyang, Changchun and Harbin in northeastern China as experimental areas, the multitemporal landscape index showed that the integration and aggregation of land cover in the urban areas had an increasing trend, the natural environment of Shenyang and Harbin was improving, while Changchun laid more emphasis on the construction of artificial facilities. At the same time, the method proposed in this paper to extract built-up areas from multi-source city data showed that the user accuracy, production accuracy, overall accuracy and Kappa coefficient are at least 3%, 1%, 1% and 0.04 higher than the single-source data method.