Fig. 14. The NTL change maps, the time series plot of a selected pixel, and high-resolution images for Milan, Italy. (a) The accumulated annual NTL change maps from 2013 to 2021 with the latest detected change year presented. (b) The day-of-year NTL change map in 2014 over the region enlarged from the white rectangle from Fig. 14a. (c) The time series plot of a selected pixel (red cross in Fig. 14b) and the corresponding VZA-COLD detection results, in which the red, green, blue, and grey colors indicate the different VZA intervals, lines represent the estimated models, small dots are the DNB observations, and the large magenta dot is the detected changes. (d) The high-resolution Google Earth image in March 2013. The red rectangle represents the location and size of the selected pixel. (e) The high-resolution Google Earth image in July 2015. The red rectangle represents the location and size of the selected pixel. The dark background in Fig. 14a and Fig. 14b is the Esri ArcMap Dark Gray Canvas base map.