This paper describes a simplified method of mapping of the ur-ban environment surface to obtain a map of the thermal anomalies distribu-tion and study the structure of the urban heat island. Those maps allow evaluating or planning the urban microclimate optimization methods and studying the effect of land cover type on the site temperature. The article discusses the processing of five satellite images for the summer and winter from 2002 to 2019. We propose a simpler and more automated processing of thermal images for Landsat 7 and Landsat 8. The stages of automatic atmospheric correction according to the DOS1 method and calculation of the emissivity with surface classification are considered. Image processing was carried out in the QGIS software package using the Semiautomatic Classification Plugin extension. As a result, thermal anomalies in Chelya-binsk were localized and a comparison of the thermal map for the specific region before and after urbanization was made.