In land monitoring applications, clouds and shadows are considered noise that should be removed as automatically and quickly as possible, before further analysis. This paper presents a method to detect clouds and shadows in Suomi NPP satellite’s VIIRS (Visible Infrared Imaging Radiometer Suite) satellite images. The proposed cloud and shadow detection method has two distinct features when compared to many other methods. First, the method does not use the thermal bands and can thus be applied to other sensors which do not contain thermal channels, such as Sentinel-2 data. Secondly, the method uses the ratio between blue and green reflectance to detect shadows. Seven hundred and forty-seven VIIRS images over Fennoscandia from August 2014 to April 2016 were processed to train and develop the method. Twenty four points from every tenth of the images were used in accuracy assessment. These 1752 points were interpreted visually to cloud, cloud shadow and clear classes, then compared to the output of the cloud and shadow detection. The comparison on VIIRS images showed 94.2% correct detection rates and 11.1% false alarms for clouds, and respectively 36.1% and 82.7% for shadows. The results on cloud detection were similar to state-of-the-art methods. Shadows showed correctly on the northern edge of the clouds, but many shadows were wrongly assigned to other classes in some cases (e.g., to water class on lake and forest boundary, or with shadows over cloud). This may be due to the low spatial resolution of VIIRS images, where shadows are only a few pixels wide and contain lots of mixed pixels.

A method to detect abrupt land cover changes using hierarchical clustering of multi-temporal satellite imagery was developed. The Autochange method outputs the pre-change land cover class, the change magnitude, and the change type. Pre-change land cover information is transferred to post-change imagery based on classes derived by unsupervised clustering, enabling using data from different instruments for pre- and post-change. The change magnitude and change types are computed by unsupervised clustering of the post-change image within each cluster, and by comparing the mean intensity values of the lower level clusters with their parent cluster means. A computational approach to determine the change magnitude threshold for the abrupt change was developed. The method was demonstrated with three summer image pairs Sentinel-2/Sentinel-2, Landsat 8/Sentinel-2, and Sentinel-2/ALOS 2 PALSAR in a study area of 12,372 km² in southern Finland for the detection of forest clear cuts and tested with independent data. The Sentinel-2 classification produced an omission error of 5.6% for the cut class and 0.4% for the uncut class. Commission errors were 4.9% for the cut class and 0.4% for the uncut class. For the Landsat 8/Sentinel-2 classifications the equivalent figures were 20.8%, 0.2%, 3.4%, and 1.6% and for the Sentinel-2/ALOS PALSAR classification 16.7%, 1.4%, 17.8%, and 1.3%, respectively. The Autochange algorithm and its software implementation was considered applicable for the mapping of abrupt land cover changes using multi-temporal satellite data. It allowed mixing of images even from the optical and synthetic aperture radar (SAR) sensors in the same change analysis.