Remote sensing data analysis in the cloudy area is still a challenging process. Fortunately, remote sensing technology is fast growing. As a result, multitemporal data could be used to overcome the problem of the cloudy area. Using multitemporal data is a common approach to address the cloud problem. However, most methods only use two data, one as the main data and the other as complementary of the cloudy area. In this paper, a method to harness multitemporal remote sensing data for automatically extracting some indices is proposed. In this method, the process of extracting the indices is done without having to mask the cloud. Those indices could be further used for many applications such as the classification of urban built-up. Landsat-8 data that is acquired during 2019 are stacked, therefore each pixel at the same position creates a list. From each list, indices are extracted. In this study, NDVI, NDBI, and NDWI are used to mapping built-up areas. Furthermore, extracted indices are divided into four categories by their value (maximum, quantile 75, median, and mean). Those indices are then combined into a simple formula to mapping built-up to see which produces better accuracy. The Pleiades as high-resolution remote sensing data is used to assist supervised classification for assessment. In this study, the combination of mean NDBI, maximum NDVI, and mean NDWI result highest Kappa coefficient of 0.771.