Detection and mapping the spatial distribution of bracken fern weeds using the Landsat 8 OLI new generation sensor

Bracken fern is an invasive plant that presents serious environmental, ecological and economic problems around the world. An understanding of the spatial distribution of bracken fern weeds is therefore essential for providing appropriate management strategies at both local and regional scales. The aim of this study was to assess the utility of the freely available medium resolution Landsat 8 OLI sensor in the detection and mapping of bracken fern at the Cathedral Peak, South Africa. To achieve this objective, the results obtained from Landsat 8 OLI were compared with those derived using the costly, high spatial resolution WorldView-2 imagery. Since previous studies have already successfully mapped bracken fern using high spatial resolution WorldView-2 image, the comparison was done to investigate the magnitude of difference in accuracy between the two sensors in relation to their acquisition costs. To evaluate the performance of Landsat 8 OLI in discriminating bracken fern compared to that of Worldview-2, we tested the utility of (i) spectral bands; (ii) derived vegetation indices as well as (iii) the combination of spectral bands and vegetation indices based on discriminant analysis classification algorithm. After resampling the training and testing data and reclassifying several times (n = 100) based on the combined data sets, the overall accuracies for both Landsat 8 and WorldView-2 were tested for significant differences based on Mann-Whitney U test. The results showed that the integration of the spectral bands and derived vegetation indices yielded the best overall classification accuracy (80.08% and 87.80% for Landsat 8 OLI and WorldView-2 respectively). Additionally, the use of derived vegetation indices as a standalone data set produced the weakest overall accuracy results of 62.14% and 82.11% for both the Landsat 8 OLI and WorldView-2 images. There were significant differences {U (100) = 569.5, z = −10.8242, p < 0.01} between the classification accuracies derived based on Landsat OLI 8 and those derived using WorldView-2 sensor. Although there were significant differences between Landsat and WorldView-2 accuracies, the magnitude of variation (9%) between the two sensors was within an acceptable range. Therefore, the findings of this study demonstrated that the recently launched Landsat 8 OLI multispectral sensor provides valuable information that could aid in the long term continuous monitoring and formulation of effective bracken fern management with acceptable accuracies that are comparable to those obtained from the high resolution WorldView-2 commercial sensor.     

Mapping trees in high resolution imagery across large areas using locally variable thresholds guided by medium resolution tree maps

Large area tree maps, important for environmental monitoring and natural resource management, are often based on medium resolution satellite imagery. These data have difficulty in detecting trees in fragmented woodlands, and have significant omission errors in modified agricultural areas. High resolution imagery can better detect these trees, however, as most high resolution imagery is not normalised it is difficult to automate a tree classification method over large areas. The method developed here used an existing medium resolution map derived from either Landsat or SPOT5 satellite imagery to guide the classification of the high resolution imagery. It selected a spatially-variable threshold on the green band, calculated based on the spatially-variable percentage of trees in the existing map of tree cover. The green band proved more consistent at classifying trees across different images than several common band combinations. The method was tested on 0.5 m resolution imagery from airborne digital sensor (ADS) imagery across New South Wales (NSW), Australia using both Landsat and SPOT5 derived tree maps to guide the threshold selection. Accuracy was assessed across 6 large image mosaics revealing a more accurate result when the more accurate tree map from SPOT5 imagery was used. The resulting maps achieved an overall accuracy with 95% confidence intervals of 93% (90–95%), while the overall accuracy of the previous SPOT5 tree map was 87% (86–89%). The method reduced omission errors by mapping more scattered trees, although it did increase commission errors caused by dark pixels from water, building shadows, topographic shadows, and some soils and crops. The method allows trees to be automatically mapped at 5 m resolution from high resolution imagery, provided a medium resolution tree map already exists.     

Resolution vs. image quality in pre-tsunami imagery used for tsunami impact models in Aceh,Indonesia

Land cover roughness coefficients (LCRs) have been used in multivariate spatial models to test the mitigation potential of coastal vegetation to reduce impacts of the 2004 tsunami in Aceh, Indonesia. Previously, a Landsat 2002 satellite imagery was employed to derive land cover maps, which were then combined with vegetation characteristics, i.e., stand height, stem diameter and planting density to obtain LCRs. The present study tested LCRs extracted from 2003 and 2004 Landsat (30 m) images as well as a combination of 2003 and 2004 higher spatial resolution SPOT (10 m) imagery, while keeping the previous vegetation characteristics. Transects along the coast were used to extract land cover, whenever availability and visibility allowed. These new LCRs applied in previously developed tsunami impact models on wave outreach, casualties and damages confirmed previous findings regarding distance to the shoreline as a main factor reducing tsunami impacts. Nevertheless, the models using the new LCRs did not perform better than the original one. Particularly casualties models using 2002 LCRs performed better (δAIC > 2) than the more recent Landsat and SPOT counterparts. Cloud cover at image acquisition for Landsat and low area coverage for SPOT images decreased statistical predictive power (fewer observations). Due to the large spatial heterogeneity of tsunami characteristics as well as topographic and land-use features, it was more important to cover a larger area. Nevertheless, if more land cover classes would be referenced and high resolution imagery with low cloud cover would be available, the full benefits of higher spatial resolution imagery used to extract more precise land use roughness coefficients could be exploited.     

Updating Landsat-based forest cover maps with MODIS images using multiscale spectral-spatial-temporal superresolution mapping

With the high deforestation rates of global forest covers during the past decades, there is an ever-increasing need to monitor forest covers at both fine spatial and temporal resolutions. Moderate Resolution Imaging Spectroradiometer (MODIS) and Landsat series images have been used commonly for satellite-derived forest cover mapping. However, the spatial resolution of MODIS images and the temporal resolution of Landsat images are too coarse to observe forest cover at both fine spatial and temporal resolutions. In this paper, a novel multiscale spectral-spatial-temporal superresolution mapping (MSSTSRM) approach is proposed to update Landsat-based forest maps by integrating current MODIS images with the previous forest maps generated from Landsat image. Both the 240 m MODIS bands and 480 m MODIS bands were used as inputs of the spectral energy function of the MSSTSRM model. The principle of maximal spatial dependence was used as the spatial energy function to make the updated forest map spatially smooth. The temporal energy function was based on a multiscale spatial-temporal dependence model, and considers the land cover changes between the previous and current time. The novel MSSTSRM model was able to update Landsat-based forest maps more accurately, in terms of both visual and quantitative evaluation, than traditional pixel-based classification and the latest sub-pixel based super-resolution mapping methods The results demonstrate the great efficiency and potential of MSSTSRM for updating fine temporal resolution Landsat-based forest maps using MODIS images.     

Phenology from Landsat when data is scarce: Using MODIS and Dynamic Time-Warping to combine multi-year Landsat imagery to derive annual phenology curves

Green-leaf phenology describes the development of vegetation throughout a growing season and greatly affects the interaction between climate and the biosphere. Remote sensing is a valuable tool to characterize phenology over large areas but doing at fine- to medium resolution (e.g., with Landsat data) is difficult because of low numbers of cloud-free images in a single year. One way to overcome data availability limitations is to merge multi-year imagery into one time series, but this requires accounting for phenological differences among years. Here we present a new approach that employed a time series of a MODIS vegetation index data to quantify interannual differences in phenology, and Dynamic Time Warping (DTW) to re-align multi-year Landsat images to a common phenology that eliminates year-to-year phenological differences. This allowed us to estimate annual phenology curves from Landsat between 2002 and 2012 from which we extracted key phenological dates in a Monte-Carlo simulation design, including green-up (GU), start-of-season (SoS), maturity (Mat), senescence (Sen), end-of-season (EoS) and dormancy (Dorm). We tested our approach in eight locations across the United States that represented forests of different types and without signs of recent forest disturbance. We compared Landsat-based phenological transition dates to those derived from MODIS and ground-based camera data from the PhenoCam-network. The Landsat and MODIS comparison showed strong agreement. Dates of green-up, start-of-season and maturity were highly correlated (r 0.86-0.95), as were senescence and end-of-season dates (r > 0.85) and dormancy (r > 0.75). Agreement between the Landsat and PhenoCam was generally lower, but correlation coefficients still exceeded 0.8 for all dates. In addition, because of the high data density in the new Landsat time series, the confidence intervals of the estimated keydates were substantially lower than in case of MODIS and PhenoCam. Our study thus suggests that by exploiting multi-year Landsat imagery and calibrating it with MODIS data it is possible to describe green-leaf phenology at much finer spatial resolution than previously possible, highlighting the potential for fine scale phenology maps using the rich Landsat data archive over large areas.     

Evaluating MODIS snow products for modelling snowmelt runoff: Case study of the Rio Grande headwaters

Snow-covered area (SCA) is a key variable in the Snowmelt-Runoff Model (SRM) and in other models for simulating discharge from snowmelt. Landsat Thematic Mapper (TM), Enhanced Thematic Mapper (ETM + ) or Operational Land Imager (OLI) provide remotely sensed data at an appropriate spatial resolution for mapping SCA in small headwater basins, but the temporal resolution of the data is low and may not always provide sufficient cloud-free dates. The coarser spatial resolution Moderate Resolution Imaging Spectroradiometer (MODIS) offers better temporal resolution and in cloudy years, MODIS data offer the best alternative for mapping snow cover when finer spatial resolution data are unavailable. However, MODIS’ coarse spatial resolution (500 m) can obscure fine spatial patterning in snow cover and some MODIS products are not sensitive to end-of-season snow cover. In this study, we aimed to test MODIS snow products for use in simulating snowmelt runoff from smaller headwater basins by a) comparing maps of TM and MODIS-based SCA and b) determining how SRM streamflow simulations are changed by the different estimates of seasonal snow depletion. We compared gridded MODIS snow products (Collection 5 MOD10A1 fractional and binary SCA; SCA derived from Collection 6 MOD10A1 Normalised Difference Snow Index (NDSI) Snow Cover), and the MODIS Snow Covered-Area and Grain size retrieval (MODSCAG) canopy-corrected fractional SCA (SCA_MG), with reference SCA maps (SCA_REF) generated from binary classification of TM imagery. SCA_MG showed strong agreement with SCA_REF; excluding true negatives (where both methods agreed no snow was present) the median percent difference between SCA_REF and SCA_MG ranged between −2.4% and 4.7%. We simulated runoff for each of the four study years using SRM populated with and calibrated for snow depletion curves derived from SCA_REF. We then substituted in each of the MODIS-derived depletion curves. With efficiency coefficients ranging between 0.73 and 0.93, SRM simulation results from the SCA_MG runs yielded the best results of all the MODIS products and only slightly underestimated discharge volume (between 7 and 11% of measured annual discharge). SRM simulations that used SCA derived from Collection 6 NDSI Snow Cover also yielded promising results, with efficiency coefficients ranging between 0.73 and 0.91.In conclusion, we recommend that when simulating snowmelt runoff from small basins (<4000 km²) with SRM, we recommend that users select either canopy-corrected MODSCAG or create their own site-specific products from the Collection 6 MOD10A1 NDSI.     

Very high spatial resolution optical and radar imagery in tracking water level fluctuations of a small inland reservoir

Tracking water level fluctuations in small lakes and reservoirs is important in order to better understand and manage these ecosystems. A geographic object-based image analysis (GEOBIA) method using very high spatial and temporal resolution optical (Pléiades) and radar (COSMO-SkyMed and TerraSAR-X) remote sensing imagery is presented here which (1) tracks water level fluctuations via variations in water surface area and (2) avoids common difficulties found in using single-band radar images for water-land image classification. Results are robust, with over 98% of image surface area correctly classified into land or water, R² = 0.963 and RMSE = 0.42 m for a total water level fluctuation range of 5.94 m. Multispectral optical imagery is found to be more straightforward in producing results than single-band radar imagery, but the latter crucially increase temporal resolution to the point where fluctuations can be satisfactorily tracked in time. Moreover, an analysis suggests that high and medium spatial resolution imagery is sufficient, in at least some cases, in tracking the water level fluctuations of small inland reservoirs. Finally, limitations of the methodology presented here are briefly discussed along with potential solutions to overcome them.     

Improving the mapping of crop types in the Midwestern U.S. by fusing Landsat and MODIS satellite data

Mapping crop types is of great importance for assessing agricultural production, land-use patterns, and the environmental effects of agriculture. Indeed, both radiometric and spatial resolution of Landsat’s sensors images are optimized for cropland monitoring. However, accurate mapping of crop types requires frequent cloud-free images during the growing season, which are often not available, and this raises the question of whether Landsat data can be combined with data from other satellites. Here, our goal is to evaluate to what degree fusing Landsat with MODIS Nadir Bidirectional Reflectance Distribution Function (BRDF)-Adjusted Reflectance (NBAR) data can improve crop-type classification. Choosing either one or two images from all cloud-free Landsat observations available for the Arlington Agricultural Research Station area in Wisconsin from 2010 to 2014, we generated 87 combinations of images, and used each combination as input into the Spatial and Temporal Adaptive Reflectance Fusion Model (STARFM) algorithm to predict Landsat-like images at the nominal dates of each 8-day MODIS NBAR product. Both the original Landsat and STARFM-predicted images were then classified with a support vector machine (SVM), and we compared the classification errors of three scenarios: 1) classifying the one or two original Landsat images of each combination only, 2) classifying the one or two original Landsat images plus all STARFM-predicted images, and 3) classifying the one or two original Landsat images together with STARFM-predicted images for key dates. Our results indicated that using two Landsat images as the input of STARFM did not significantly improve the STARFM predictions compared to using only one, and predictions using Landsat images between July and August as input were most accurate. Including all STARFM-predicted images together with the Landsat images significantly increased average classification error by 4% points (from 21% to 25%) compared to using only Landsat images. However, incorporating only STARFM-predicted images for key dates decreased average classification error by 2% points (from 21% to 19%) compared to using only Landsat images. In particular, if only a single Landsat image was available, adding STARFM predictions for key dates significantly decreased the average classification error by 4 percentage points from 30% to 26% (p < 0.05). We conclude that adding STARFM-predicted images can be effective for improving crop-type classification when only limited Landsat observations are available, but carefully selecting images from a full set of STARFM predictions is crucial. We developed an approach to identify the optimal subsets of all STARFM predictions, which gives an alternative method of feature selection for future research.     

Integrating age in the detection and mapping of incongruous patches in coffee (Coffea arabica) plantations using multi-temporal Landsat 8 NDVI anomalies

The development of cost-effective, reliable and easy to implement crop condition monitoring methods is urgently required for perennial tree crops such as coffee (Coffea arabica), as they are grown over large areas and represent long term and higher levels of investment. These monitoring methods are useful in identifying farm areas that experience poor crop growth, pest infestation, diseases outbreaks and/or to monitor response to management interventions. This study compares field level coffee mean NDVI and LSWI anomalies and age-adjusted coffee mean NDVI and LSWI anomalies in identifying and mapping incongruous patches across perennial coffee plantations. To achieve this objective, we first derived deviation of coffee pixels from the global coffee mean NDVI and LSWI values of nine sequential Landsat 8 OLI image scenes. We then evaluated the influence of coffee age class (young, mature and old) on Landsat-scale NDVI and LSWI values using a one-way ANOVA and since results showed significant differences, we adjusted NDVI and LSWI anomalies for age-class. We then used the cumulative inverse distribution function (α ≤ 0.05) to identify fields and within field areas with excessive deviation of NDVI and LSWI from the global and the age-expected mean for each of the Landsat 8 OLI scene dates spanning three seasons. Results from accuracy assessment indicated that it was possible to separate incongruous and healthy patches using these anomalies and that using NDVI performed better than using LSWI for both global and age-adjusted mean anomalies. Using the age-adjusted anomalies performed better in separating incongruous and healthy patches than using the global mean for both NDVI (Overall accuracy = 80.9% and 68.1% respectively) and for LSWI (Overall accuracy = 68.1% and 48.9% respectively). When applied to other Landsat 8 OLI scenes, the results showed that the proportions of coffee fields that were modelled incongruent decreased with time for the young age category and while it increased for the mature and old age classes with time. We concluded that the method could be useful for the identification of anomalous patches using Landsat scale time series data to monitor large coffee plantations and provide an indication of areas requiring particular field attention.     

Assessing the potential of multi-seasonal WorldView-2 imagery for mapping West African agroforestry tree species

High resolution satellite systems enable efficient and detailed mapping of tree cover, with high potential to support both natural resource monitoring and ecological research. This study investigates the capability of multi-seasonal WorldView-2 imagery to map five dominant tree species at the individual tree crown level in a parkland landscape in central Burkina Faso. The Random Forest algorithm is used for object based tree species classification and for assessing the relative importance of WorldView-2 predictors. The classification accuracies from using wet season, dry season and multi-seasonal datasets are compared to gain insights about the optimal timing for image acquisition. The multi-seasonal dataset produced the most accurate classifications, with an overall accuracy (OA) of 83.4%. For classifications based on single date imagery, the dry season (OA = 78.4%) proved to be more suitable than the wet season (OA = 68.1%). The predictors that contributed most to the classification success were based on the red edge band and visible wavelengths, in particular green and yellow. It was therefore concluded that WorldView-2, with its unique band configuration, represents a suitable data source for tree species mapping in West African parklands. These results are particularly promising when considering the recently launched WorldView-3, which provides data both at higher spatial and spectral resolution, including shortwave infrared bands.     

Forest cover trends from time series Landsat data for the Australian continent

In perennial and natural vegetation systems, monitoring changes in vegetation over time is of fundamental interest for identifying and quantifying impacts of management and natural processes. Subtle changes in vegetation cover can be identified by calculating the trends of a vegetation density index over time. In this paper, we apply such an index-trends approach, which has been developed and applied to time series Landsat imagery in rangeland and woodland environments, to continental-scale monitoring of disturbances within forested regions of Australia. This paper describes the operational methods used for the generation of National Forest Trend (NFT) information, which is a time-series summary providing visual indication of within-forest vegetation changes (disturbance and recovery) over time at 25 m resolution. This result is based on a national archive of calibrated Landsat TM/ETM+ data from 1989 to 2006 produced for Australia's National Carbon Accounting System (NCAS). The NCAS was designed in 1999 initially to provide consistent fine-scale classifications for monitoring forest cover extent and changes (i.e. land use change) over the Australian continent using time series Landsat imagery. NFT information identifies more subtle changes within forested areas and provides a capacity to identify processes affecting forests which are of primary interest to ecologists and land managers. The NFT product relies on the identification of an appropriate Landsat-based vegetation cover index (defined as a linear combination of spectral image bands) that is sensitive to changes in forest density. The time series of index values at a location, derived from calibrated imagery, represents a consistent surrogate to track density changes. To produce the trends summary information, statistical summaries of the index response over time (such as slope and quadratic curvature) are calculated. These calculated index responses of woody vegetation cover are then displayed as maps where the different colours indicate the approximate timing, direction (decline or increase), magnitude and spatial extent of the changes in vegetation cover. These trend images provide a self-contained and easily interpretable summary of vegetation change at scales that are relevant for natural resource management (NRM) and environmental reporting.     

Towards decadal soil salinity mapping using Landsat time series data

Salinization is one of the major soil problems around the world. However, decadal variation in soil salinization has not yet been extensively reported. This study exploited thirty years (1985–2015) of Landsat sensor data, including Landsat-4/5 TM (Thematic Mapper), Landsat-7 ETM+ (Enhanced Thematic Mapper Plus) and Landsat-8 OLI (Operational Land Imager), for monitoring soil salinity of the Yellow River Delta, China. The data were initially corrected for atmospheric effects, and then matched the spectral bands of EO-1 (Earth Observing One) ALI (Advanced Land Imager). Subsequently, soil salinity maps were derived with a previously developed PLSR (Partial Least Square Regression) model. On intra-annual scale, the retrievals showed that soil salinity increased in February, stabilized in March, and decreased in April. On inter-annual scale, soil salinity decreased within 1985–2000 (−0.74 g kg⁻¹/10a, p < 0.001), and increased within 2000–2015 (0.79 g kg⁻¹/10a, p < 0.001). Our study presents a new perspective for use of multiple Landsat data in soil salinity retrieval, and further the understanding of soil salinization development over the Yellow River Delta.     

Monitoring land surface albedo and vegetation dynamics using high spatial and temporal resolution synthetic time series from Landsat and the MODIS BRDF/NBAR/albedo product

Seasonal vegetation phenology can significantly alter surface albedo which in turn affects the global energy balance and the albedo warming/cooling feedbacks that impact climate change. To monitor and quantify the surface dynamics of heterogeneous landscapes, high temporal and spatial resolution synthetic time series of albedo and the enhanced vegetation index (EVI) were generated from the 500 m Moderate Resolution Imaging Spectroradiometer (MODIS) operational Collection V006 daily BRDF/NBAR/albedo products and 30 m Landsat 5 albedo and near-nadir reflectance data through the use of the Spatial and Temporal Adaptive Reflectance Fusion Model (STARFM). The traditional Landsat Albedo (Shuai et al., 2011) makes use of the MODIS BRDF/Albedo products (MCD43) by assigning appropriate BRDFs from coincident MODIS products to each Landsat image to generate a 30 m Landsat albedo product for that acquisition date. The available cloud free Landsat 5 albedos (due to clouds, generated every 16 days at best) were used in conjunction with the daily MODIS albedos to determine the appropriate 30 m albedos for the intervening daily time steps in this study. These enhanced daily 30 m spatial resolution synthetic time series were then used to track albedo and vegetation phenology dynamics over three Ameriflux tower sites (Harvard Forest in 2007, Santa Rita in 2011 and Walker Branch in 2005). These Ameriflux sites were chosen as they are all quite nearby new towers coming on line for the National Ecological Observatory Network (NEON), and thus represent locations which will be served by spatially paired albedo measures in the near future. The availability of data from the NEON towers will greatly expand the sources of tower albedometer data available for evaluation of satellite products. At these three Ameriflux tower sites the synthetic time series of broadband shortwave albedos were evaluated using the tower albedo measurements with a Root Mean Square Error (RMSE) less than 0.013 and a bias within the range of ±0.006. These synthetic time series provide much greater spatial detail than the 500 m gridded MODIS data, especially over more heterogeneous surfaces, which improves the efforts to characterize and monitor the spatial variation across species and communities. The mean of the difference between maximum and minimum synthetic time series of albedo within the MODIS pixels over a subset of satellite data of Harvard Forest (16 km by 14 km) was as high as 0.2 during the snow-covered period and reduced to around 0.1 during the snow-free period. Similarly, we have used STARFM to also couple MODIS Nadir BRDF Adjusted Reflectances (NBAR) values with Landsat 5 reflectances to generate daily synthetic times series of NBAR and thus Enhanced Vegetation Index (NBAR-EVI) at a 30 m resolution. While normally STARFM is used with directional reflectances, the use of the view angle corrected daily MODIS NBAR values will provide more consistent time series. These synthetic times series of EVI are shown to capture seasonal vegetation dynamics with finer spatial and temporal details, especially over heterogeneous land surfaces.     

Texture classification of Mediterranean land cover

Maximum likelihood (ML) and artificial neural network (ANN) classifiers were applied to three Landsat Thematic Mapper (TM) image sub-scenes (termed urban, agricultural and semi-natural) of Cukurova, Turkey. Inputs to the classifications comprised (i) spectral data and (ii) spectral data in combination with texture measures derived on a per-pixel basis. The texture measures used were: the standard deviation and variance and statistics derived from the co-occurrence matrix and the variogram. The addition of texture measures increased classification accuracy for the urban sub-scene but decreased classification accuracy for agricultural and semi-natural sub-scenes. Classification accuracy was dependent on the nature of the spatial variation in the image sub-scene and, in particular, the relation between the frequency of spatial variation and the spatial resolution of the imagery. For Mediterranean land, texture classification applied to Landsat TM imagery may be appropriate for the classification of urban areas only.     

Evaluating metrics derived from Landsat 8 OLI imagery to map crop cover

Developing techniques are required to generate agricultural land cover maps to monitor agricultural fields. Landsat 8 Operational Land Imager (OLI) offers reflectance data over the visible to shortwave-infrared range. OLI offers several advantages, such as adequate spatial and spectral resolution, and 16 day repeat coverage, furthermore, spectral indices derived from Landsat 8 OLI possess great potential for evaluating the status of vegetation. Additionally, classification algorithms are essential for generating accurate maps. Recently, multi-Grained Cascade Forest, which is also called deep forest, was proposed, and it was shown to give highly competitive performance for classification. However, the ability of this algorithm to generate crop maps with satellite data had not yet been evaluated. In this study, the reflectance at 7 bands and 57 spectral indices calculated from Landsat 8 OLI data were evaluated for its potential for crop type identification.     

Improved coastal wetland mapping using very-high 2-meter spatial resolution imagery

Accurate wetland maps are a fundamental requirement for land use management and for wetland restoration planning. Several wetland map products are available today; most of them based on remote sensing images, but their different data sources and mapping methods lead to substantially different estimations of wetland location and extent. We used two very high-resolution (2 m) WorldView-2 satellite images and one (30 m) Landsat 8 Operational Land Imager (OLI) image to assess wetland coverage in two coastal areas of Tampa Bay (Florida): Fort De Soto State Park and Weedon Island Preserve. An initial unsupervised classification derived from WorldView-2 was more accurate at identifying wetlands based on ground truth data collected in the field than the classification derived from Landsat 8 OLI (82% vs. 46% accuracy). The WorldView-2 data was then used to define the parameters of a simple and efficient decision tree with four nodes for a more exacting classification. The criteria for the decision tree were derived by extracting radiance spectra at 1500 separate pixels from the WorldView-2 data within field-validated regions. Results for both study areas showed high accuracy in both wetland (82% at Fort De Soto State Park, and 94% at Weedon Island Preserve) and non-wetland vegetation classes (90% and 83%, respectively). Historical, published land-use maps overestimate wetland surface cover by factors of 2–10 in the study areas. The proposed methods improve speed and efficiency of wetland map production, allow semi-annual monitoring through repeat satellite passes, and improve the accuracy and precision with which wetlands are identified.     

Assessing the potential of integrated Landsat 8 thermal bands,with the traditional reflective bands and derived vegetation indices in classifying urban landscapes

Reliable and up-to-date urban land cover information is valuable in urban planning and policy development. Due to the increasing demand for reliable land cover information there has been a growing need for robust methods and datasets to improve the classification accuracy from remotely sensed imagery. This study sought to assess the potential of the newly launched Landsat 8 sensor’s thermal bands and derived vegetation indices in improving land cover classification in a complex urban landscape using the support vector machine classifier. This study compared the individual and combined performance of Landsat 8’s reflective, thermal bands and vegetation indices in classifying urban land use-land cover. The integration of Landsat 8 reflective bands, derived vegetation indices and thermal bands overall produced significantly higher accuracy classification results than using traditional bands as standalone (i.e. overall, user and producer accuracies). An overall accuracy above 89.33% and a kappa index of 0.86, significantly higher than the one obtained with the use of the traditional reflective bands as a standalone data-set and other analysis stages. On average, the results also indicate high producer and user accuracies (i.e. above 80%) for most of the classes with a McNemar’s Z score of 9.00 at 95% confidence interval showing significant improvement compared with classification using reflective bands as standalone. Overall, the results of this study indicate that the integration of the Landsat 8’s OLI and TIR data presents an invaluable potential for accurate and robust land cover classification in a complex urban landscape, especially in areas where the availability of high resolution datasets remains a challenge.     

Multitemporal fusion of Landsat/TM and ENVISAT/MERIS for crop monitoring

Monitoring Earth dynamics using current and future satellites is one of the most important objectives of the remote sensing community. The exploitation of image time series from sensors with different characteristics provides new opportunities to increase the knowledge about environmental changes and to support many operational applications. This paper presents an image fusion approach based on multiresolution and multisensor regularized spatial unmixing. The approach yields a composite image with the spatial resolution of the high spatial resolution image while retaining the spectral and temporal characteristics of the medium spatial resolution image. The approach is tested using images from Landsat/TM and ENVISAT/MERIS instruments, but is general enough to be applied to other sensor pairs. The potential of the proposed spatial unmixing approach is illustrated in an agricultural monitoring application where Landsat temporal profiles from images acquired over Albacete, Spain, in 2004 and 2009 are complemented with MERIS fused images. The resulting spatial resolution from Landsat allows monitoring small and medium size crops at the required scale while the fine spectral and temporal resolution from MERIS allow a more accurate determination of the crop type and phenology as well as capturing rapidly varying land-cover changes.     

Pixel-based and object-based classifications using high- and medium-spatial-resolution imageries in the urban and suburban landscapes

With the increasing availability of high-spatial-resolution remote sensing imageries and with the observed limitations of pixel-based techniques, the development and testing of geographic object-based image analysis (GEOBIA) techniques for image classification have become one of the main research areas in geospatial science. This paper examines and compares the classification performance of a pixel-based method and an object-based method as applied to high- (QuickBird satellite image) and medium- (Landsat TM image) spatial-resolution imageries in the context of urban and suburban landscapes. For the pixel-based classification, the maximum-likelihood supervised classification approach was employed. And for the object-based classification, the pixel-based classified maps were integrated with a set of image segments produced using various calibrations. The results show evidence that the object-based method can produce classifications that are more accurate for both high- and medium-spatial- resolution imageries in the context of urban and suburban landscapes.     

Empirical models for estimating the suspended sediment concentration in Amazonian white water rivers using Landsat 5/TM

Suspended sediment yield is a very important environmental indicator within Amazonian fluvial systems, especially for rivers dominated by inorganic particles, referred to as white water rivers. For vast portions of Amazonian rivers, suspended sediment concentration (SSC) is measured infrequently or not at all. However, remote sensing techniques have been used to estimate water quality parameters worldwide, from which data for suspended matter is the most successfully retrieved. This paper presents empirical models for SSC retrieval in Amazonian white water rivers using reflectance data derived from Landsat 5/TM. The models use multiple regression for both the entire dataset (global model, N = 504) and for five segmented datasets (regional models) defined by general geological features of drainage basins. The models use VNIR bands, band ratios, and the SWIR band 5 as input. For the global model, the adjusted R² is 0.76, while the adjusted R² values for regional models vary from 0.77 to 0.89, all significant (p-value < 0.0001). The regional models are subject to the leave-one-out cross validation technique, which presents robust results. The findings show that both the average error of estimation and the standard deviation increase as the SSC range increases. Regional models were more accurate when compared with the global model, suggesting changes in optical proprieties of water sampled at different sampling stations. Results confirm the potential for the estimation of SSC from Landsat/TM historical series data for the 1980s and 1990s, for which the in situ database is scarce. Such estimates supplement the SSC temporal series, providing a more comprehensive SSC temporal series which may show environmental dynamics yet unknown.