Object-Based Downscaling Method for Land Surface Temperature with High-Spatial-Resolution Multispectral Data

Land surface temperature (LST) is an important environmental parameter in many fields. However, many studies require high-spatial- and high-temporal-resolution LST products to improve the coarse spatial resolution of moderate-resolution imaging spectroradiometer (MODIS) LSTs. Numerous approaches have downscaled MODIS LST images to a finer spatial resolution using pixel-based image analysis (PBA). Meanwhile, object-based image analysis (OBIA) methods, which have developed rapidly in the analysis of high-spatial-resolution visible and near-infrared (VNIR) band data, have received little attention in the LST downscaling field. In this paper, we propose an object-based downscaling (OBD) method for MODIS LST using high-spatial-resolution multispectral images (e.g., Landsat Thematic Mapper (TM), Enhanced Thematic Mapper Plus (ETM+), and Advanced Spaceborne Thermal Emission and Reflection Radiometer (ASTER)) as auxiliary data. The fundamental principle of this method is to preserve the thermal radiance of the “object”, which is composed of several MODIS LST pixels (partly or entirely) and is unchanged after disaggregation into subpixels in the resulting LST image. The decomposition process consists of two key parts: the thermal radiance (TR) estimation of the object from MODIS LST products and the weight calculation of sub-objects or subpixels. Objects were generated from VNIR data and remote sensing indices (e.g., the normalized difference vegetation index (NDVI), the normalized difference built-up index (NDBI), and fractions of different endmembers) using a multiscale segmentation method. The radiance of subpixels or sub-objects was calculated based on the weights of their parent objects, which were estimated by the relationships between the remote sensing indices and the LST. The accuracy and the efficiency of the OBD method were validated using a pair of ASTER and MODIS datapoints that were acquired at the same time. The decomposed LST results showed that the spatial distribution of the downscaled LST image closely resembled the true LST of the ASTER, with an RMSE of 2.5 K for the entire image. A comparison with PBA methods for pixel downscaling also indicated that the OBD method achieves the lowest root mean square error (RMSE) across different landcovers, including urban areas, water bodies, and natural terrain. Therefore, the proposed OBD method significantly enhances the capability of increasing the spatial resolution of coarse MODIS LST, providing an alternative for improving the spatial resolution of MODIS LST images and expanding their applicability to studies that require high-temporal- and high-spatial-resolution LST products.