Landsat time series commonly contain missing observations, i.e., gaps, due to the orbit and sensing geometry, data acquisition strategy, and cloud contamination. A spectral-angle-mapper (SAM) based spatio-temporal similarity (SAMSTS) gap-filling algorithm is presented that is designed to fill small and large area gaps in Landsat data, using one year or less of data and without using other satellite data. Each gap pixel is filled by an alternative similar pixel that is located in a non-missing region of the image. The alternative similar pixel locations are identified by comparison of reflectance time series using a SAM metric revised to be adaptive to missing observations. A time series segmentation-and-clustering approach is used to increase the search efficiency. The SAMSTS algorithm is demonstrated using six months of Landsat 8 Operational Land Imager (OLI) reflectance time series over three 150 × 150 km (5000 × 5000 30 m pixels) areas in California, Minnesota and Kansas. The three areas contain different land cover types, especially crops that have different phenology and abrupt changes due to agricultural harvesting, which make gap filling challenging. Fillings on simulated gaps, which are equivalent to 36% of 5000 × 5000 images in each test area, are presented. The gap filling accuracy is assessed quantitatively, and the SAMSTS algorithm is shown to perform better than the simple closest temporal pixel substitution gap filling approach and the sinusoidal harmonic model-based gap filling approach. The SAMSTS algorithm provides gap-filled data with five-band reflective-wavelength root-mean-square differences less the 0.02, which is comparable to the OLI reflectance calibration accuracy.Temporal interpolation (TI) gap-filling approaches have been developed that fit time series statistical models to predict reflectance or vegetation index values on a given day. Linear, logistic, or sum of sinusoidal models have been used [32][33][34][35][36][37][38][39]. The model fits are conducted on individual pixel time series. Model fitting is sensitive to the quality, number of available observations, and the seasonality of missing observations. TI methods are typically less reliable for surfaces that have abrupt changes, for example, due to land cover change, flooding, fire, or for surfaces with complex phenology, such as double or triple agricultural cropping [35,40,41]. Time series change detection methods that account for both abrupt and gradual trends in coarse spatial resolution data have been developed [41], but are less well suited for Landsat application as they require a higher observation frequency than provided by Landsat.The ASP gap-filling approach fills a gap pixel with the values of one or more alternative pixels selected from non-missing pixels found usually in the same image. Alternative similar pixel locations are identified from a reference image that may be the same, a previous, or a subsequent image that has no gap at the location to be filled [42][43][44][45][46][47][48][49][50]. The ASP method was ...