Remote Sensing — Open Access Journal

An investigation has been carried out to examine the crosstalk contamination in the Visible Infrared Imaging Radiometer Suite (VIIRS) onboard the Suomi National Polar-orbiting Partnership (SNPP) spacecraft. Prior to this study, the cause of the pronounced striping in Earth View (EV) images and obvious discontinuity in the EV brightness temperature (BT) of the thermal emissive bands (TEB) during black body (BB) warm-up cool-down (WUCD) calibration observed since launch has not been identified. Meanwhile, it has been recently demonstrated in the MODerate-resolution Imaging Spectroradiometer (MODIS) long-wave infrared (LWIR) photovoltaic (PV) bands that the crosstalk effect induces the same erroneous features. In this investigation, it is shown that the established lunar imagery analysis indeed verifies the existence of crosstalk contamination in SNPP VIIRS TEB. The crosstalk effect is quantitatively characterized by deriving the crosstalk coefficients from the scheduled lunar observations. The magnitude of the effect is comparatively smaller than that in MODIS LWIR PV bands, but is of a large enough magnitude to induce the aforementioned artificial features. Among all SNPP VIIRS TEB, Band M14 has the largest crosstalk contamination from Band M15, while Bands M13, M15, M16, and I5 have pronounced crosstalk effects as well. One new detail of the crosstalk effect specific to SNPP VIIRS, differing from the MODIS result, is the distinctive two-group pattern of odd and even detectors for each affected band due to the arrangement of the detector on the focal plane assembly (FPA). This is fully consistent with the earlier finding that this odd-even detector arrangement contributes to striping in the sea surface temperature (SST) products. Our analyses additionally suggest an explanation of the large temperature anomalies appearing during the WUCD time periods. The parallel effort examining the potential crosstalk contamination in SNPP VIIRS reflective solar bands, however, reveals no observable effect. Full article

Eutrophication is an increasing problem in coastal waters of the Baltic Sea. Moreover, algal blooms, which occur every summer in the Gulf of Gdansk can deleteriously impact human health, the aquatic environment, and economically important fisheries, tourism, and recreation industries. Traditional laboratory-based techniques for water monitoring are expensive and time consuming, which usually results in limited numbers of observations and discontinuity in space and time. The use of hyperspectral radiometers for coastal water observation provides the potential for more detailed remote optical monitoring. A statistical approach to develop local models for the estimation of optically significant components from in situ measured hyperspectral remote sensing reflectance in case 2 waters is presented in this study. The models, which are based on empirical orthogonal function (EOF) analysis and stepwise multilinear regression, allow for the estimation of parameters strongly correlated with phytoplankton (pigment concentration, absorption coefficient) and coloured detrital matter abundance (absorption coefficient) directly from reflectance spectra measured in situ. Chlorophyll a concentration, which is commonly used as a proxy for phytoplankton biomass, was retrieved with low error (median percent difference, MPD = 17%, root mean square error RMSE = 0.14 in log₁₀ space) and showed a high correlation with chlorophyll a measured in situ (R = 0.84). Furthermore, phycocyanin and phycoerythrin, both characteristic pigments for cyanobacteria species, were also retrieved reliably from reflectance with MPD = 23%, RMSE = 0.23, R² = 0.77 and MPD = 24%, RMSE = 0.15, R² = 0.74, respectively. The EOF technique proved to be accurate in the derivation of the absorption spectra of phytoplankton and coloured detrital matter (CDM), with R² (λ) above 0.83 and RMSE around 0.10. The approach was also applied to satellite multispectral remote sensing reflectance data, thus allowing for improved temporal and spatial resolution compared with the in situ measurements. The EOF method tested on simulated Medium Resolution Imaging Spectrometer (MERIS) or Ocean and Land Colour Instrument (OLCI) data resulted in RMSE = 0.16 for chl-a and RMSE = 0.29 for phycocyanin. The presented methods, applied to both in situ and satellite data, provide a powerful tool for coastal monitoring and management. Full article

A polarimetric ground-based radar (GB-radar) system operated in the L-band has been developed. The frequency range of GB-radar is 1.215 to 1.3 GHz, which is the same as that of Japanese satellite-borne SAR, PALSAR-2 (Phased Array type L-band Synthetic Aperture Radar), and Japanese L-band air-borne radar, i.e., Pi-SAR-L2 (Polarimetric and Interferometric Airborne Synthetic Aperture Radar L2). Polarimetric calibration was carried out twice in the field to calibrate and validate the GB-radar data. Cross-talk and channel imbalance are improved for the both experiments, and are from −13.3 dB to −30.7 dB, and from 1.06 to 1.00, respectively for the first experiment, after calibration. The calibrated cross-talk and channel imbalance values were comparable to −31.7 dB and 1.013, which were obtained using PALSAR. Radiometric calibration and antenna pattern correction were also carried out in the second experiment. Forest observations were also carried out simultaneously by GB-radar and Pi-SAR-L2 in the second experiment. The range profiles obtained by GB-radar and Pi-SAR-L2 were compared for several polarimetric parameters, namely, radar backscattering coefficient, polarimetric coherence, eigenvalue-decomposition parameters, and four-component-decomposition parameters. Both range profiles matched moderately well and showed good performance that could compensate for the limited possibility of satellite/air-borne SAR observation. Full article

The Tibetan Plateau, the world’s largest orogenic plateau, hosts thousands of lakes that play prominent roles as water resources, environmental archives, and sources of natural hazards such as glacier lake outburst floods. Previous studies have reported that the size of lakes on the Tibetan Plateau has changed rapidly in recent years, possibly because of atmospheric warming. Tracking these changes systematically with remote sensing data is challenging given the different spectral signatures of water, the potential for confusing lakes with glaciers, and difficulties in classifying frozen or partly frozen lakes. Object-based image analysis (OBIA) offers new opportunities for automated classification in this context, and we have explored this method for mapping lakes from LANDSAT images and Shuttle Radar Topography Mission (SRTM) elevation data. We tested our algorithm for most of the Tibetan Plateau, where lakes in tectonic depressions or blocked by glaciers and sediments have different surface colours and seasonal ice cover in images obtained in 1995 and 2015. We combined a modified normalised difference water index (MNDWI) with OBIA and local topographic slope data in order to classify lakes with an area >10 km². Our method derived 323 water bodies, with a total area of 31,258 km², or 2.6% of the study area (in 2015). The same number of lakes had covered only 24,892 km² in 1995; lake area has increased by ~26% in the past two decades. The classification had estimated producer’s and user’s accuracies of 0.98, with a Cohen’s kappa and F-score of 0.98, and may thus be a useful approximation for quantifying regional hydrological budgets. We have shown that our method is flexible and transferable to detecting lakes in diverse physical settings on several continents with similar success rates. Full article

The Moderate Resolution Imaging Spectroradiometer (MODIS) provides widespread Aerosol Optical Depth (AOD) datasets for climatological and environmental health research. Since MODIS AOD clearly lacks coverage in orbit-scanning gaps and cloud obscuration, some applications will benefit from data recovery using multi-temporal AOD. Aimed at qualitatively describing the relationship between multi-temporal AOD, AOD loadings and Normalized Difference Vegetation Index (NDVI) have been considered based on the mechanism of satellite AOD retrieval. Accordingly, the NDVI-based Weighted Linear Regression (NWLR) has been proposed to recover AOD by synthetically weighing AOD similarity, spatial proximity, and NDVI similarity. To evaluate the performance of AOD recovery, simulated experiments applying gap and window masks were conducted in South Asia and Beijing, respectively. The evaluation results demonstrated that the linear regression R2 achieved 0.8 and the absolute relative errors remained steady. Further validation was conducted between the recovered and actual AODs using 56 Aerosol Robotic Network (AERONET) sites in East and South Asia from 2013 to 2015, which demonstrated that over 41% of recovered AODs fell within the expected error (EE) envelope. Additional validation conducted in South Asia and Beijing showed that recovery by NWLR did not expand satellite-derived AOD errors, and the accuracy of recovered AOD was consistent with the accuracy of the original Aqua MODIS Deep Blue (DB) AOD. The recovery results illustrated that AOD coverage was improved in most regions, especially in North China, Mongolia, and South Asia, which could provide better support in aerosol spatio-temporal analysis and aerosol data assimilation. Full article

Medium spatial resolution biomass is a crucial link from the plot to regional and global scales. Although remote-sensing data-based methods have become a primary approach in estimating forest above ground biomass (AGB), many difficulties remain in data resources and prediction approaches. Each kind of sensor type and prediction method has its own merits and limitations. To select the proper estimation algorithm and remote-sensing data source, several forest AGB models were developed using different remote-sensing data sources (Geoscience Laser Altimeter System (GLAS) data and Thematic Mapper (TM) data) and 108 field measurements. Three modeling methods (stepwise regression (SR), support vector regression (SVR) and random forest (RF)) were used to estimate forest AGB over the Daxing’anling Mountains in northeastern China. The results of models using different datasets and three approaches were compared. The random forest AGB model using Landsat5/TM as input data was shown the acceptable modeling accuracy (R2 = 0.95 RMSE = 17.73 Mg/ha) and it was also shown to estimate AGB reliably by cross validation (R2 = 0.71 RMSE = 39.60 Mg/ha). The results also indicated that adding GLAS data significantly improved AGB predictions for the SVR and SR AGB models. In the case of the RF AGB models, including GLAS data no longer led to significant improvement. Finally, a forest biomass map with spatial resolution of 30 m over the Daxing'anling Mountains was generated using the obtained optimal model. Full article

Due to the side-looking imaging geometry of the current synthetic aperture radar (SAR) sensors, only ground deformation along the radar’s line-of-sight (LOS) and azimuth directions can be potentially obtained from a single amplitude pair (SAP) of SAR using offset tracking (OT) procedures. This significantly hinders the accurate assessment of mining-related hazards and better understanding of the mining subsidence mechanism. In this paper, we propose a method for completely retrieving three-dimensional (3-D) mining-induced displacements with OT-derived observations of LOS deformation from a single amplitude pair of SAR (referred to as OT-SAP hereinafter). The OT-SAP method first constructs two extra constraints at each pixel of the mining area based on the proportional relationship between the horizontal motion of the mining area and the gradients of the vertical subsidence in the east and north directions. The full 3-D mining-induced displacements are then solved by coupling the two constructed extra constraints with the OT-derived observations of the LOS deformation. The Daliuta coal mining area in China was selected to test the proposed OT-SAP method. The results show that the maximum 3-D displacements of this mining area were about 4.3 m, 1.1 m, and 1.3 m in the vertical, east, and north directions, respectively, from 21 November 2012 to 6 February 2013. The accuracies of the retrieved displacements in the vertical and horizontal directions are about 0.201 m and 0.214 m, respectively, which are much smaller than the mining-induced displacements in this mining area and can satisfy the basic requirements of mining deformation monitoring. Full article

In this paper, we present the supervised multi-view canonical correlation analysis ensemble (SMVCCAE) and its semi-supervised version (SSMVCCAE), which are novel techniques designed to address heterogeneous domain adaptation problems, i.e., situations in which the data to be processed and recognized are collected from different heterogeneous domains. Specifically, the multi-view canonical correlation analysis scheme is utilized to extract multiple correlation subspaces that are useful for joint representations for data association across domains. This scheme makes homogeneous domain adaption algorithms suitable for heterogeneous domain adaptation problems. Additionally, inspired by fusion methods such as Ensemble Learning (EL), this work proposes a weighted voting scheme based on canonical correlation coefficients to combine classification results in multiple correlation subspaces. Finally, the semi-supervised MVCCAE extends the original procedure by incorporating multiple speed-up spectral regression kernel discriminant analysis (SRKDA). To validate the performances of the proposed supervised procedure, a single-view canonical analysis (SVCCA) with the same base classifier (Random Forests) is used. Similarly, to evaluate the performance of the semi-supervised approach, a comparison is made with other techniques such as Logistic label propagation (LLP) and the Laplacian support vector machine (LapSVM). All of the approaches are tested on two real hyperspectral images, which are considered the target domain, with a classifier trained from synthetic low-dimensional multispectral images, which are considered the original source domain. The experimental results confirm that multi-view canonical correlation can overcome the limitations of SVCCA. Both of the proposed procedures outperform the ones used in the comparison with respect to not only the classification accuracy but also the computational efficiency. Moreover, this research shows that canonical correlation weighted voting (CCWV) is a valid option with respect to other ensemble schemes and that because of their ability to balance diversity and accuracy, canonical views extracted using partially joint random view generation are more effective than those obtained by exploiting disjoint random view generation. Full article

The combination of remote sensing and sensor network technologies can provide unprecedented earth observation capabilities, and has attracted high R&D interest in recent years. However, the procedures and tools used for deployment, geo-referenciation and collection of logged measurements in the case of traditional environmental monitoring stations are not suitable when dealing with hundreds or thousands of sensor nodes deployed in an environment of tenths of hectares. This paper presents a scheme based on Unmanned Aerial Systems that intends to give a step forward in the use of sensor networks for environment observation. The presented scheme includes methods, tools and technologies to solve sensor node deployment, localization and collection of measurements. The presented scheme is scalable—it is suitable for medium–large environments with a high number of sensor nodes—and highly autonomous—it is operated with very low human intervention. This paper presents the scheme including its main components, techniques and technologies, and describes its implementation and evaluation in field experiments. Full article

Hyperspectral image (HSI) clustering is generally a challenging task because of the complex spectral-spatial structure. Based on the assumption that all the pixels are sampled from the union of subspaces, recent works have introduced a robust technique—the sparse subspace clustering (SSC) algorithm and its enhanced versions (SSC models incorporating spatial information)—to cluster HSIs, achieving excellent performances. However, these methods are all based on the linear representation model, which conflicts with the well-known nonlinear structure of HSIs and limits their performance to a large degree. In this paper, to overcome these obstacles, we present a new kernel sparse subspace clustering algorithm with a spatial max pooling operation (KSSC-SMP) for hyperspectral remote sensing data interpretation. The proposed approach maps the feature points into a much higher dimensional kernel space to extend the linear sparse subspace clustering model to nonlinear manifolds, which can better fit the complex nonlinear structure of HSIs. With the help of the kernel sparse representation, a more accurate representation coefficient matrix can be obtained. A spatial max pooling operation is utilized for the representation coefficients to generate more discriminant features by integrating the spatial-contextual information, which is essential for the accurate modeling of HSIs. This paper is an extension of our previous conference paper, and a number of enhancements are put forward. The proposed algorithm was evaluated on two well-known hyperspectral data sets—the Salinas image and the University of Pavia image—and the experimental results clearly indicate that the newly developed KSSC-SMP algorithm can obtain very competitive clustering results for HSIs, outperforming the current state-of-the-art clustering methods. Full article

Wind farms (WFs) can affect the local climate, and local climate change may influence underlying vegetation. Some studies have shown that WFs affect certain aspects of the regional climate, such as temperature and rainfall. However, there is still no evidence to demonstrate whether WFs can affect local vegetation growth, a significant part of the overall assessment of WF effects. In this research, based on the moderate-resolution imaging spectroradiometer (MODIS) vegetation index, productivity and other remote-sensing data from 2003 to 2014, the effects of WFs in the Bashang area of Northern China on vegetation growth and productivity in the summer (June–August) were analyzed. The results showed that: (1) WFs had a significant inhibiting effect on vegetation growth, as demonstrated by decreases in the leaf area index (LAI), the enhanced vegetation index (EVI), and the normalized difference vegetation index (NDVI) of approximately 14.5%, 14.8%, and 8.9%, respectively, in the 2003–2014 summers. There was also an inhibiting effect of 8.9% on summer gross primary production (GPP) and 4.0% on annual net primary production (NPP) coupled with WFs; and (2) the major impact factors might be the changes in temperature and soil moisture: WFs suppressed soil moisture and enhanced water stress in the study area. This research provides significant observational evidence that WFs can inhibit the growth and productivity of the underlying vegetation. Full article

A large part of arid areas in tropical and sub-tropical regions are dominated by sparse xerophytic vegetation, which are essential for providing products and services for local populations. While a large number of researches already exist for the derivation of wall-to-wall estimations of above ground biomass (AGB) with remotely sensed data, only a few of them are based on the direct use of non-photogrammetric aerial photography. In this contribution we present an experiment carried out in a study area located in the Santiago Island in the Cape Verde archipelago where a National Forest Inventory (NFI) was recently carried out together with a new acquisition of a visible high-resolution aerial orthophotography. We contrasted two approaches: single-tree, based on the automatic delineation of tree canopies; and area-based, on the basis of an automatic image classification. Using 184 field plots collected for the NFI we created parametric models to predict AGB on the basis of the crown projection area (CPA) estimated from the two approaches. Both the methods produced similar root mean square errors (RMSE) at pixel level 45% for the single-tree and 42% for the area-based. However, the latest was able to better predict the AGB along all the variable range, limiting the saturation problem which is evident when the CPA tends to reach the full coverage of the field plots. These findings demonstrate that in regions dominated by sparse vegetation, a simple aerial orthophoto can be used to successfully create AGB wall-to-wall predictions. The level of these estimations’ uncertainty permits the derivation of small area estimations useful for supporting a more correct implementation of sustainable management practices of wood resources. Full article

Landslide mapping (LM) has recently become an important research topic in remote sensing and geohazards. The area near the Three Gorges Reservoir (TGR) along the Yangtze River in China is one of the most landslide-prone regions in the world, and the area has suffered widespread and significant landslide events in recent years. In our study, an object-oriented landslide mapping (OOLM) framework was proposed for reliable and accurate LM from ‘ZY-3’ high spatial resolution (HSR) satellite images. The framework was based on random forests (RF) and mathematical morphology (MM). RF was first applied as an object feature information reduction tool to identify the significant features for describing landslides, and it was then combined with MM to map the landslides. Three object-feature domains were extracted from the ‘ZY-3’ HSR data: layer information, texture, and geometric features. A total group of 124 features and 24 landslides were used as inputs to determine the landslide boundaries and evaluate the landslide classification accuracy. The results showed that: (1) the feature selection (FS) method had a positive influence on effective landslide mapping; (2) by dividing the data into two sets, training sets which consisted of 20% of the landslide objects (O_LS) and non-landslide objects (O_NLS), and test sets which consisted of the remaining 80% of the O_LS and O_NLS, the selected feature subsets were combined for training to obtain an overall classification accuracy of 93.3% ± 0.12% of the test sets; (3) four MM operations based on closing and opening were used to improve the performance of the RF classification. Seven accuracy evaluation indices were used to compare the accuracies of these landslide mapping methods. Finally, the landslide inventory maps were obtained. Based on its efficiency and accuracy, the proposed approach can be employed for rapid response to natural hazards in the Three Gorges area. Full article

The segmentation of urban scene mobile laser scanning (MLS) data into meaningful street objects is a great challenge due to the scene complexity of street environments, especially in the vicinity of street objects such as poles and trees. This paper proposes a three-stage method for the segmentation of urban MLS data at the object level. The original unorganized point cloud is first voxelized, and all information needed is stored in the voxels. These voxels are then classified as ground and non-ground voxels. In the second stage, the whole scene is segmented into clusters by applying a density-based clustering method based on two key parameters: local density and minimum distance. In the third stage, a merging step and a re-assignment processing step are applied to address the over-segmentation problem and noise points, respectively. We tested the effectiveness of the proposed methods on two urban MLS datasets. The overall accuracies of the segmentation results for the two test sites are 98.3% and 97%, thereby validating the effectiveness of the proposed method. Full article

Dengue fever (DF), a vector-borne flavivirus, is endemic to the tropical countries of the world with nearly 400 million people becoming infected each year and roughly one-third of the world’s population living in areas of risk. The main vector for DF is the Aedes aegypti mosquito, which is also the same vector of yellow fever, chikungunya, and Zika viruses. To gain an understanding of the spatial aspects that can affect the epidemiological processes across the disease’s geographical range, and the spatial interactions involved, we created and compared Bernoulli and Poisson family Boosted Regression Tree (BRT) models to quantify the overall annual risk of DF incidence by municipality, using the Magdalena River watershed of Colombia as a study site during the time period between 2012 and 2014. A wide range of environmental conditions make this site ideal to develop models that, with minor adjustments, could be applied in many other geographical areas. Our results show that these BRT methods can be successfully used to identify areas at risk and presents great potential for implementation in surveillance programs. Full article

Accurate and reliable information about the situation in refugee or internally displaced person camps is very important for planning any kind of help like health care, infrastructure, or vaccination campaigns. The number and spatial distribution of single dwellings extracted semi-automatically from very high-resolution (VHR) satellite imagery as an indicator for population estimations can provide such important information. The accuracy of the extracted dwellings can vary quite a lot depending on various factors. To enhance established single dwelling extraction approaches, we have tested the integration of stratified template matching methods in object-based image analysis (OBIA) workflows. A template library for various dwelling types (template samples are taken from ten different sites using 16 satellite images), incorporating the shadow effect of dwellings, was established. Altogether, 18 template classes were created covering typically occurring dwellings and their cast shadows. The created template library aims to be generally applicable in similar conditions. Compared to pre-existing OBIA classifications, the approach could increase the producer’s accuracy by 11.7 percentage points on average and slightly increase the user’s accuracy. These results show that the stratified integration of template matching approaches in OBIA workflows is a possibility to further improve the results of semi-automated dwelling extraction, especially in complex situations. Full article

Geographic Object-Based Image Analysis (GEOBIA) techniques have become increasingly popular in remote sensing. GEOBIA has been claimed to represent a paradigm shift in remote sensing interpretation. Still, GEOBIA—similar to other emerging paradigms—lacks formal expressions and objective modelling structures and in particular semantic classification methods using ontologies. This study has put forward an object-based semantic classification method for high resolution satellite imagery using an ontology that aims to fully exploit the advantages of ontology to GEOBIA. A three-step workflow has been introduced: ontology modelling, initial classification based on a data-driven machine learning method, and semantic classification based on knowledge-driven semantic rules. The classification part is based on data-driven machine learning, segmentation, feature selection, sample collection and an initial classification. Then, image objects are re-classified based on the ontological model whereby the semantic relations are expressed in the formal languages OWL and SWRL. The results show that the method with ontology—as compared to the decision tree classification without using the ontology—yielded minor statistical improvements in terms of accuracy for this particular image. However, this framework enhances existing GEOBIA methodologies: ontologies express and organize the whole structure of GEOBIA and allow establishing relations, particularly spatially explicit relations between objects as well as multi-scale/hierarchical relations. Full article

The Soil Moisture Active Passive (SMAP) mission was designed to provide global mapping of soil moisture (SM) on nested 3, 9, and 36 km earth grids measured by L-band passive and active microwave sensors. The validation of SMAP SM products is crucial for the application of the products and improvement of the retrieval algorithm. Since the SMAP SM products were released, much effort has been invested in the evaluation of the SMAP radiometer SM product (SMAP_P). However, there has been little validation of SMAP radar (SMAP_A) and active/passive combined (SMAP_AP) SM products. This paper presents an evaluation of SMAP_P, SMAP_A and SMAP_AP SM products by using distributed ground observations networks in different landscapes in the Heihe River Basin of northwestern China. The standard error metrics of SMAP products and relative error are applied to measure the products’ performances. The results show that the SMAP SM products exhibit consistent spatial-temporal variation with the ground measurements and typical precipitation events. Three products show various types of performance capability (e.g., active, passive and combined), surface coverage (e.g., bare, vegetated) and climatic region (e.g., cold, arid). Relatively, the SMAP_P shows the best performance, while the SMAP_A performs the worst. The best performances are observed over bare soils but with overestimation and the largest relative error, and unsatisfactory accuracies are observed over cold regions and woody vegetated surfaces with underestimation. The vegetation effect and the freezing-thawing cycle may be major factors that led to an unsatisfactory performance. Efforts on resolving the influence of these factors are expected to improve the accuracy and to promote the application of SMAP SM products over these regions. Overall, this evaluation provides an understanding of SMAP SM products over cold and arid regions, and suggestions for the further refinement of the SMAP SM retrieval algorithms. Full article

On 24 August 2016, the M 6.2 Amatrice earthquake struck central Italy, well-known as a seismically active region, causing considerable damage to buildings in the town of Amatrice and the surrounding area. Damage from this earthquake was assessed quantitatively by means of multitemporal synthetic aperture radar (SAR) coherence and SAR intensity methods using dual-polarized SAR data obtained from the Sentinel-1 (VV, VH) and ALOS-2 (HH, HV) satellites. We developed linear discriminant functions based on three items: (1) the differential coherence values; (2) the differential backscattering intensity values of pre- and post-event images; and (3) a binary damage map of the optical pre- and post-event imagery. The accuracy of the proposed model was 84% for the Sentinel-1 data and 76% for the ALOS-2 data. The damage proxy maps deduced from the linear discriminant functions can be useful in the parcel-by-parcel assessment of building damage and development of spatial models for the allocation of urban search and rescue operations. Full article

We performed a multi-date composite change detection technique using a dense-time stack of Landsat data to map land-use and land-cover change (LCLUC) in Mainland Southeast Asia (MSEA) with a focus on the expansion of boom crops, primarily tree crops. The supervised classification was performed using Support Vector Machines (SVM), which are supervised non-parametric statistical learning techniques. To select the most suitable SMV classifier and the related parameter settings, we used the training data and performed a two-dimensional grid search with a three-fold internal cross-validation. We worked in seven Landsat footprints and found the linear kernel to be the most suitable for all footprints, but the most suitable regularization parameter C varied across the footprints. We distinguished a total of 41 LCLUCs (13 to 31 classes per footprint) in very dynamic and heterogeneous landscapes. The approach proved useful for distinguishing subtle changes over time and to map a variety of land covers, tree crops, and transformations as long as sufficient training points could be collected for each class. While to date, this approach has only been applied to mapping urban extent and expansion, this study shows that it is also useful for mapping change in rural settings, especially when images from phenologically relevant acquisition dates are included. Full article