Fig. 9. The NTL change maps, the time series plot of a selected pixel, and high-resolution images for an urbanized suburban area in Melbourne, Australia. (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 maps over the region enlarged from the white rectangle from Fig. 9a. (c) The time series plot of a selected pixel (red cross in Fig. 9b) 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 change. (d) The high-resolution Google Earth image in January 2015. The red rectangle represents the location and size of the selected pixel. (e) The high-resolution Google Earth image in March 2020. The red rectangle represents the location and size of the selected pixel. The dark background in Fig. 9a and Fig. 9b is the Esri ArcMap Dark Gray Canvas base map.