Understanding the spatial–temporal patterns of air pollution is crucial for mitigation strategies, a task fostered nowadays by the generation of continuous concentration maps by remote sensing technologies. We applied spatial modelling to analyze such spatial–temporal patterns in Lombardy, Italy, one of the most polluted regions in Europe. We conducted monthly spatial autocorrelation (global and local) of the daily average concentrations of PM2.5, PM10, O3, NO2, SO2, and CO from 2016 to 2020, using 10 × 10 km satellite data from the Copernicus Atmosphere Monitoring Service (CAMS), aggregated on districts of approximately 100,000 population. Land-use classes were computed on identified clusters, and the significance of the differences was evaluated through the Wilcoxon rank-sum test with Bonferroni correction. The global Moran’s I autocorrelation was overall high (>0.6), indicating a strong clustering. The local autocorrelation revealed high–high clusters of PM2.5 and PM10 in the central urbanized zones in winter (January–December), and in the agrarian southern districts in summer and autumn (May–October). The temporal decomposition showed that values of PMs are particularly high in winter. Low–low clusters emerged in the northern districts for all the pollutants except O3. Seasonal peaks for O3 occurred in the summer months, with high–high clusters mostly in the hilly and mildly urban districts in the northwest. These findings elaborate the spatial patterns of air pollution concentration, providing insights for effective land-use-based pollution management strategies.