Despite several decades of observations of ice cover in Polish lakes, researchers have not broadly applied satellite images to date. This paper presents a temporal and spatial analysis of the variability in the occurrence of ice cover on lakes in the Drawskie Lakeland in the hydrological years 1984–2022 based on satellite data from Landsat missions 4, 5, 7, 8, and 9. The range of occurrence of ice cover was determined based on the value of the Normalised Difference Snow Index (NDSI) and blue spectral band (ρλblue). The determination of ice cover extent adopted ρλblue values from 0.033 to 0.120 as the threshold values. The analysis covered 67 lakes with an area from 0.07 to 18.71 km2. A total of 53 images were analysed, 14 and 39 out of which showed full and partial ice cover, respectively. The cluster analysis permitted the designation of two groups of lakes characterised by an approximate range of ice cover. The obtained results were analysed in the context of the morphometric parameters of the lakes. It was evidenced that the range of the ice cover on lakes is determined by the surface area of the lakes; their mean and maximum depth, volume, length, and width; and the height of the location above sea level. The results of analyses of the spatial range of ice cover in subsequent scenes allowed for the preparation of maps of probability of ice cover occurrence that permit the complete determination of its variability within each of the lakes. Monitoring of the spatial variability in ice cover within individual lakes as well as in reference to lakes not subject to traditional observations offers new research possibilities in many scientific disciplines focused on these ecosystems.