Empirical climate classification is a process that makes environmental conditions understandable to humans by using climatic elements. Köppen-Geiger (KG) is a popular climate classification method that uses long-term precipitation and temperature data to classify climate into five primary groups. However, long-term continuous meteorological data is heavily exposed to data scarcity, particularly in a national scale. This research study addresses this challenge by leveraging satellite imageries, multilinear regression models and spatial interpolation within the context of entire country of Iran between 2016 and 2019. Accordingly, this study examined statistical relationship between 14 explanatory variables under four main categories of MODIS-LST, MODIS-NDVI, MODIS-TVDI, GPM-precipitation and SRTM-DEM against ground-based precipitation and temperature data (dependent variables). The spatial interpolation model (i.e. Krigging and CO-krigging) was directly developed from weather observation station datasets. A total of 332 synoptic stations were selected, 67% of which were used in modeling and the remaining 33% in testing. Accuracy assessment was performed with Kappa statistics. Overall, this research study developed three KG classification maps. These include a map per precipitation and temperature from regression model and spatial interpolation and a point-based maps from unused climate data in modelling. This study identified three KG main climate groups of arid, warm temperate and snow and eight KG sub-groups of hot desert, cold steppe, cold desert, hot steppe, warm temperate climate with dry hot summer, Snow climate with dry hot summer, Warm temperate climate with dry warm summer and Snow climate with dry warm summer. A comparison between those maps (kappa = 0.75) showed the higher accuracy of regression-based KG maps against spatial interpolation maps. This study contributes to a more detailed monitor of climate change across countries and regions with sparse distribution of weather observation data.