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.     

Assessment of RapidEye vegetation indices for estimation of leaf area index and biomass in corn and soybean crops

Leaf area index (LAI) and biomass are important indicators of crop development and the availability of this information during the growing season can support farmer decision making processes. This study demonstrates the applicability of RapidEye multi-spectral data for estimation of LAI and biomass of two crop types (corn and soybean) with different canopy structure, leaf structure and photosynthetic pathways. The advantages of Rapid Eye in terms of increased temporal resolution (∼daily), high spatial resolution (∼5 m) and enhanced spectral information (includes red-edge band) are explored as an individual sensor and as part of a multi-sensor constellation. Seven vegetation indices based on combinations of reflectance in green, red, red-edge and near infrared bands were derived from RapidEye imagery between 2011 and 2013. LAI and biomass data were collected during the same period for calibration and validation of the relationships between vegetation indices and LAI and dry above-ground biomass. Most indices showed sensitivity to LAI from emergence to 8 m²/m². The normalized difference vegetation index (NDVI), the red-edge NDVI and the green NDVI were insensitive to crop type and had coefficients of variations (CV) ranging between 19 and 27%; and coefficients of determination ranging between 86 and 88%. The NDVI performed best for the estimation of dry leaf biomass (CV = 27% and r² = 090) and was also insensitive to crop type. The red-edge indices did not show any significant improvement in LAI and biomass estimation over traditional multispectral indices. Cumulative vegetation indices showed strong performance for estimation of total dry above-ground biomass, especially for corn (CV ≤ 20%). This study demonstrated that continuous crop LAI monitoring over time and space at the field level can be achieved using a combination of RapidEye, Landsat and SPOT data and sensor-dependant best-fit functions. This approach eliminates/reduces the need for reflectance resampling, VIs inter-calibration and spatial resampling.     

Using a stochastic gradient boosting algorithm to analyse the effectiveness of Landsat 8 data for montado land cover mapping: Application in southern Portugal

This study aims to develop and propose a methodological approach for montado ecosystem mapping using Landsat 8 multi-spectral data, vegetation indices, and the Stochastic Gradient Boosting (SGB) algorithm. Two Landsat 8 scenes (images from spring and summer 2014) of the same area in southern Portugal were acquired. Six vegetation indices were calculated for each scene: the Enhanced Vegetation Index (EVI), the Short-Wave Infrared Ratio (SWIR32), the Carotenoid Reflectance Index 1 (CRI1), the Green Chlorophyll Index (CI_green), the Normalised Multi-band Drought Index (NMDI), and the Soil-Adjusted Total Vegetation Index (SATVI). Based on this information, two datasets were prepared: (i) Dataset I only included multi-temporal Landsat 8 spectral bands (LS8), and (ii) Dataset II included the same information as Dataset I plus vegetation indices (LS8 + VIs). The integration of the vegetation indices into the classification scheme resulted in a significant improvement in the accuracy of Dataset II’s classifications when compared to Dataset I (McNemar test: Z-value = 4.50), leading to a difference of 4.90% in overall accuracy and 0.06 in the Kappa value. For the montado ecosystem, adding vegetation indices in the classification process showed a relevant increment in producer and user accuracies of 3.64% and 6.26%, respectively. By using the variable importance function from the SGB algorithm, it was found that the six most prominent variables (from a total of 24 tested variables) were the following: EVI_summer; CRI1_spring; SWIR32_spring; B6_summer; B5_summer; and CI_green_summer.     

Spatially detailed retrievals of spring phenology from single-season high-resolution image time series

Vegetation indices derived from satellite image time series have been extensively used to estimate the timing of phenological events like season onset. Medium spatial resolution (≥250 m) satellite sensors with daily revisit capability are typically employed for this purpose. In recent years, phenology is being retrieved at higher resolution (≤30 m) in response to increasing availability of high-resolution satellite data. To overcome the reduced acquisition frequency of such data, previous attempts involved fusion between high- and medium-resolution data, or combinations of multi-year acquisitions in a single phenological reconstruction. The objectives of this study are to demonstrate that phenological parameters can now be retrieved from single-season high-resolution time series, and to compare these retrievals against those derived from multi-year high-resolution and single-season medium-resolution satellite data. The study focuses on the island of Schiermonnikoog, the Netherlands, which comprises a highly-dynamic saltmarsh, dune vegetation, and agricultural land. Combining NDVI series derived from atmospherically-corrected images from RapidEye (5 m-resolution) and the SPOT5 Take5 experiment (10m-resolution) acquired between March and August 2015, phenological parameters were estimated using a function fitting approach. We then compared results with phenology retrieved from four years of 30 m Landsat 8 OLI data, and single-year 100 m Proba-V and 250 m MODIS temporal composites of the same period. Retrieved phenological parameters from combined RapidEye/SPOT5 displayed spatially consistent results and a large spatial variability, providing complementary information to existing vegetation community maps. Retrievals that combined four years of Landsat observations into a single synthetic year were affected by the inclusion of years with warmer spring temperatures, whereas adjustment of the average phenology to 2015 observations was only feasible for a few pixels due to cloud cover around phenological transition dates. The Proba-V and MODIS phenology retrievals scaled poorly relative to their high-resolution equivalents, indicating that medium-resolution phenology retrievals need to be interpreted with care, particularly in landscapes with fine-scale land cover variability.     

Landsat 8 vs. Landsat 5: A comparison based on urban and peri-urban land cover mapping

An image dataset from the Landsat OLI spaceborne sensor is compared with the Landsat TM in order to evaluate the excellence of the new imagery in urban landcover classification. Widely known pixel-based and object-based image analysis methods have been implemented in this work like Maximum Likelihood, Support Vector Machine, k-Nearest Neighbor, Feature Analyst and Sub-pixel. Classification results from Landsat OLI provide more accurate results comparing to the Landsat TM. Object-based classifications produced a more uniform result, but suffer from the absorption of small rare classes into large homogenous areas, as a consequence of the segmentation, merging and the spatial parameters in the spatial resolution (30 m) of Landsat images. Based exclusively on the overall accuracy reports, the SVM pixel-based classification from Landsat 8 proved to be the most accurate for the purpose of mapping urban land cover, using medium spatial resolution imagery.     

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.     

Testing the discrimination and detection limits of WorldView-2 imagery on a challenging invasive plant target

Invasive plants pose significant threats to biodiversity and ecosystem function globally, leading to costly monitoring and management effort. While remote sensing promises cost-effective, robust and repeatable monitoring tools to support intervention, it has been largely restricted to airborne platforms that have higher spatial and spectral resolutions, but which lack the coverage and versatility of satellite-based platforms. This study tests the ability of the WorldView-2 (WV2) eight-band satellite sensor for detecting the invasive shrub mesquite (Prosopis spp.) in the north-west Pilbara region of Australia. Detectability was challenged by the target taxa being largely defoliated by a leaf-tying biological control agent (Gelechiidae: Evippe sp. #1) and the presence of other shrubs and trees. Variable importance in the projection (VIP) scores identified bands offering greatest capacity for discrimination were those covering the near-infrared, red, and red-edge wavelengths. Wavelengths between 400 nm and 630 nm (coastal blue, blue, green, yellow) were not useful for species level discrimination in this case. Classification accuracy was tested on three band sets (simulated standard multispectral, all bands, and bands with VIP scores ≥1). Overall accuracies were comparable amongst all band-sets (Kappa = 0.71–0.77). However, mesquite omission rates were unacceptably high (21.3%) when using all eight bands relative to the simulated standard multispectral band-set (9.5%) and the band-set informed by VIP scores (11.9%). An incremental cover evaluation on the latter identified most omissions to be for objects <16 m². Mesquite omissions reduced to 2.6% and overall accuracy significantly improved (Kappa = 0.88) when these objects were left out of the confusion matrix calculations. Very high mapping accuracy of objects >16 m² allows application for mapping mesquite shrubs and coalesced stands, the former not previously possible, even with 3 m resolution hyperspectral imagery. WV2 imagery offers excellent portability potential for detecting other species where spectral/spatial resolution or coverage has been an impediment. New generation satellite sensors are removing barriers previously preventing widespread adoption of remote sensing technologies in natural resource management.     

Comparison of UAV and WorldView-2 imagery for mapping leaf area index of mangrove forest

Unmanned Aerial Vehicle (UAV) remote sensing has opened the door to new sources of data to effectively characterize vegetation metrics at very high spatial resolution and at flexible revisit frequencies. Successful estimation of the leaf area index (LAI) in precision agriculture with a UAV image has been reported in several studies. However, in most forests, the challenges associated with the interference from a complex background and a variety of vegetation species have hindered research using UAV images. To the best of our knowledge, very few studies have mapped the forest LAI with a UAV image. In addition, the drawbacks and advantages of estimating the forest LAI with UAV and satellite images at high spatial resolution remain a knowledge gap in existing literature. Therefore, this paper aims to map LAI in a mangrove forest with a complex background and a variety of vegetation species using a UAV image and compare it with a WorldView-2 image (WV2).In this study, three representative NDVIs, average NDVI (AvNDVI), vegetated specific NDVI (VsNDVI), and scaled NDVI (ScNDVI), were acquired with UAV and WV2 to predict the plot level (10 × 10 m) LAI. The results showed that AvNDVI achieved the highest accuracy for WV2 (R² = 0.778, RMSE = 0.424), whereas ScNDVI obtained the optimal accuracy for UAV (R² = 0.817, RMSE = 0.423). In addition, an overall comparison results of the WV2 and UAV derived LAIs indicated that UAV obtained a better accuracy than WV2 in the plots that were covered with homogeneous mangrove species or in the low LAI plots, which was because UAV can effectively eliminate the influence from the background and the vegetation species owing to its high spatial resolution. However, WV2 obtained a slightly higher accuracy than UAV in the plots covered with a variety of mangrove species, which was because the UAV sensor provides a negative spectral response function(SRF) than WV2 in terms of the mangrove LAI estimation.     

A Spatio-Temporal Enhancement Method for medium resolution LAI (STEM-LAI)

Satellite remote sensing has been used successfully to map leaf area index (LAI) across landscapes, but advances are still needed to exploit multi-scale data streams for producing LAI at both high spatial and temporal resolution. A multi-scale Spatio-Temporal Enhancement Method for medium resolution LAI (STEM-LAI) has been developed to generate 4-day time-series of Landsat-scale LAI from existing medium resolution LAI products. STEM-LAI has been designed to meet the demands of applications requiring frequent and spatially explicit information, such as effectively resolving rapidly evolving vegetation dynamics at sub-field (30 m) scales. In this study, STEM-LAI is applied to Moderate Resolution Imaging Spectroradiometer (MODIS) based LAI data and utilizes a reference-based regression tree approach for producing MODIS-consistent, but Landsat-based, LAI. The Spatial and Temporal Adaptive Reflectance Fusion Model (STARFM) is used to interpolate the downscaled LAI between Landsat acquisition dates, providing a high spatial and temporal resolution improvement over existing LAI products. STARFM predicts high resolution LAI by blending MODIS and Landsat based information from a common acquisition date, with MODIS data from a prediction date. To demonstrate its capacity to reproduce fine-scale spatial features observed in actual Landsat LAI, the STEM-LAI approach is tested over an agricultural region in Nebraska. The implementation of a 250 m resolution LAI product, derived from MODIS 1 km data and using a scale consistent approach based on the Normalized Difference Vegetation Index (NDVI), is found to significantly improve accuracies of spatial pattern prediction, with the coefficient of efficiency (E) ranging from 0.77–0.94 compared to 0.01–0.85 when using 1 km LAI inputs alone. Comparisons against an 11-year record of in-situ measured LAI over maize and soybean highlight the utility of STEM-LAI in reproducing observed LAI dynamics (both characterized by r² = 0.86) over a range of plant development stages. Overall, STEM-LAI represents an effective downscaling and temporal enhancement mechanism that predicts in-situ measured LAI better than estimates derived through linear interpolation between Landsat acquisitions. This is particularly true when the in-situ measurement date is greater than 10 days from the nearest Landsat acquisition, with prediction errors reduced by up to 50%. With a streamlined and completely automated processing interface, STEM-LAI represents a flexible tool for LAI disaggregation in space and time that is adaptable to different land cover types, landscape heterogeneities, and cloud cover conditions.     

Estimating winter wheat biomass by assimilating leaf area index derived from fusion of Landsat-8 and MODIS data

A sufficient number of satellite acquisitions in a growing season are essential for deriving agronomic indicators, such as green leaf area index (GLAI), to be assimilated into crop models for crop productivity estimation. However, for most high resolution orbital optical satellites, it is often difficult to obtain images frequently due to their long revisit cycles and unfavorable weather conditions. Data fusion algorithms, such as the Spatial and Temporal Adaptive Reflectance Fusion Model (STARFM) and the Enhanced STARFM (ESTARFM), have been developed to generate synthetic data with high spatial and temporal resolution to address this issue. In this study, we evaluated the approach of assimilating GLAI into the Simple Algorithm for Yield Estimation model (SAFY) for winter wheat biomass estimation. GLAI was estimated using the two-band Enhanced Vegetation Index (EVI2) derived from data acquired by the Operational Land Imager (OLI) onboard the Landsat-8 and a fusion dataset generated by blending the Moderate-Resolution Imaging Spectroradiometer (MODIS) data and the OLI data using the STARFM and ESTARFM models. The fusion dataset had the temporal resolution of the MODIS data and the spatial resolution of the OLI data. Key parameters of the SAFY model were optimised through assimilation of the estimated GLAI into the crop model using the Shuffled Complex Evolution-University of Arizona (SCE-UA) algorithm. A good agreement was achieved between the estimated and field measured biomass by assimilating the GLAI derived from the OLI data (GLAI_L) alone (R² = 0.77 and RMSE = 231 g m⁻²). Assimilation of GLAI derived from the fusion dataset (GLAI_F) resulted in a R² of 0.71 and RMSE of 193 g m⁻² while assimilating the combination of GLAI_L and GLAI_F led to further improvements (R² = 0.76 and RMSE = 176 g m⁻²). Our results demonstrated the potential of using the fusion algorithms to improve crop growth monitoring and crop productivity estimation when the number of high resolution remote sensing data acquisitions is limited.     

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.     

Estimating soil salinity in Pingluo County of China using QuickBird data and soil reflectance spectra

Soil salinization is a worldwide environmental problem with severe economic and social consequences. In this paper, estimating the soil salinity of Pingluo County, China by a partial least squares regression (PLSR) predictive model was carried out using QuickBird data and soil reflectance spectra. At first, a relationship between the sensitive bands of soil salinity acquired from measured reflectance spectra and the spectral coverage of seven commonly used optical sensors was analyzed. Secondly, the potentiality of QuickBird data in estimating soil salinity by analyzing the correlations between the measured reflectance spectra and reflectance spectra derived from QuickBird data and analyzing the contributions of each band of QuickBird data to soil salinity estimation Finally, a PLSR predictive model of soil salinity was developed using reflectance spectra from QuickBird data and eight spectral indices derived from QuickBird data. The results indicated that the sensitive bands covered several bands of each optical sensor and these sensors can be used for soil salinity estimation. The result of estimation model showed that an accurate prediction of soil salinity can be made based on the PLSR method (R² = 0.992, RMSE = 0.195). The PLSR model's performance was better than that of the stepwise multiple regression (SMR) method. The results also indicated that using spectral indices such as intensity within spectral bands (Int1, Int2), soil salinity indices (SI1, SI2, SI3), the brightness index (BI), the normalized difference vegetation index (NDVI) and the ratio vegetation index (RVI) as independent model variables can help to increase the accuracy of soil salinity mapping. The NDVI and RVI can help to reduce the influences of vegetation cover and soil moisture on prediction accuracy. The method developed in this paper can be applied in other arid and semi-arid areas, such as western China.     

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.     

Modelling high arctic percent vegetation cover using field digital images and high resolution satellite data

In this study, digital images collected at a study site in the Canadian High Arctic were processed and classified to examine the spatial-temporal patterns of percent vegetation cover (PVC). To obtain the PVC of different plant functional groups (i.e., forbs, graminoids/sedges and mosses), field near infrared-green-blue (NGB) digital images were classified using an object-based image analysis (OBIA) approach. The PVC analyses comparing different vegetation types confirmed: (i) the polar semi-desert exhibited the lowest PVC with a large proportion of bare soil/rock cover; (ii) the mesic tundra cover consisted of approximately 60% mosses; and (iii) the wet sedge consisted almost exclusively of graminoids and sedges. As expected, the PVC and green normalized difference vegetation index (GNDVI; (R_NIR − R_Green)/(R_NIR + R_Green)), derived from field NGB digital images, increased during the summer growing season for each vegetation type: i.e., ∼5% (0.01) for polar semi-desert; ∼10% (0.04) for mesic tundra; and ∼12% (0.03) for wet sedge respectively. PVC derived from field images was found to be strongly correlated with WorldView-2 derived normalized difference spectral indices (NDSI; (R_x − R_y)/(R_x + R_y)), where R_x is the reflectance of the red edge (724.1 nm) or near infrared (832.9 nm and 949.3 nm) bands; R_y is the reflectance of the yellow (607.7 nm) or red (658.8 nm) bands with R²’s ranging from 0.74 to 0.81. NDSIs that incorporated the yellow band (607.7 nm) performed slightly better than the NDSIs without, indicating that this band may be more useful for investigating Arctic vegetation that often includes large proportions of senescent vegetation throughout the growing season.     

Detecting forest damage after a low-severity fire using remote sensing at multiple scales

Remote sensing technologies are an ideal platform to examine the extent and impact of fire on the landscape. In this study we assess that capacity of the RapidEye constellation and Landsat (Thematic Mapper and Operational Land Imager to map fine-scale burn attributes for a small, low severity prescribed fire in a dry Western Canadian forest. Estimates of burn severity from field data were collated into a simple burn index and correlated with a selected suite of common spectral vegetation indices. Burn severity classes were then derived to map fire impacts and estimate consumed woody surface fuels (diameter ≥2.6 cm). All correlations between the simple burn index and vegetation indices produced significant results (p < 0.01), but varied substantially in their overall accuracy. Although the Landsat Soil Adjusted Vegetation Index provided the best regression fit (R² = 0.56), results suggested that RapidEye provided much more spatially detailed estimates of tree damage (Soil Adjusted Vegetation Index, R² = 0.51). Consumption estimates of woody surface fuels ranged from 3.38 ± 1.03 Mg ha⁻¹ to 11.73 ± 1.84 Mg ha⁻¹, across four derived severity classes with uncertainties likely a result of changing foliage moisture between the before and after fire images. While not containing spectral information in the short wave infrared, the spatial variability provided by the RapidEye imagery has potential for mapping and monitoring fine scale forest attributes, as well as the potential to resolve fire damage at the individual tree level.     

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.     

Reducing background effects in orchards through spectral vegetation index correction

Satellite remote sensing provides an alternative to time-consuming and labor intensive in situ measurements of biophysical variables in agricultural crops required for precision agriculture applications. In orchards, however, the spatial resolution causes mixtures of canopies and background (i.e. soil, grass and shadow), hampering the estimation of these biophysical variables. Furthermore, variable background mixtures obstruct meaningful comparisons between different orchard blocks, rows or within each row. Current correction methodologies use spectral differences between canopies and background, but struggle with a vegetated orchard floor. This background influence and the lack of a generic solution are addressed in this study.Firstly, the problem was demonstrated in a controlled environment for vegetation indices sensitive to chlorophyll content, water content and leaf area index. Afterwards, traditional background correction methods (i.e. soil-adjusted vegetation indices and signal unmixing) were compared to the proposed vegetation index correction. This correction was based on the mixing degree of each pixel (i.e. tree cover fraction) to rescale the vegetation indices accordingly and was applied to synthetic and WorldView-2 satellite imagery. Through the correction, the effect of background admixture for vegetation indices was reduced, and the estimation of biophysical variables was improved (ΔR² = 0.2–0.31).     

Mapping seagrass coverage and spatial patterns with high spatial resolution IKONOS imagery

Seagrass habitats in subtidal coastal waters provide a variety of ecosystem functions and services and there is an increasing need to acquire information on spatial and temporal dynamics of this resource. Here, we explored the capability of IKONOS (IKO) data of high resolution (4 m) for mapping seagrass cover [submerged aquatic vegetation (%SAV) cover] along the mid-western coast of Florida, USA. We also compared seagrass maps produced with IKO data with that obtained using the Landsat TM sensor with lower resolution (30 m). Both IKO and TM data, collected in October 2009, were preprocessed to calculate water depth invariant bands to normalize the effect of varying depth on bottom spectra recorded by the two satellite sensors and further the textural information was extracted from IKO data. Our results demonstrate that the high resolution IKO sensor produced a higher accuracy than the TM sensor in a three-class % SAV cover classification. Of note is that the OA of %SAV cover mapping at our study area created with IKO data was 5–20% higher than that from other studies published. We also examined the spatial distribution of seagrass over a spatial range of 4–240 m using the Ripley’s K function [L(d)] and IKO data that represented four different grain sizes [4 m (one IKO pixel), 8 m (2 × 2 IKO pixels), 12 m (3 × 3 IKO pixels), and 16 m (4 × 4 IKO pixels)] from moderate-dense seagrass cover along a set of six transects. The Ripley’s K metric repeatedly indicated that seagrass cover representing 4 m × 4 m pixels displayed a dispersed (or slightly dispersed) pattern over distances of <4–8 m, and a random or slightly clustered pattern of cover over 9–240 m. The spatial pattern of seagrass cover created with the three additional grain sizes (i.e., 2 × 24 m IKO pixels, 3 × 34 m IKO pixels, and 4 × 4 m IKO pixels) show a dispersed (or slightly dispersed) pattern across 4–32 m and a random or slightly clustered pattern across 33–240 m. Given the first report on using satellite observations to quantify seagrass spatial patterns at a spatial scale from 4 m to 240 m, our novel analyses of moderate-dense SAV cover utilizing Ripley’s K function illustrate how data obtained from the IKO sensor revealed seagrass spatial information that would be undetected by the TM sensor with a 30 m pixel size. Use of the seagrass classification scheme here, along with data from the IKO sensor with enhanced resolution, offers an opportunity to synoptically record seagrass cover dynamics at both small and large spatial scales.     

Intercomparison and validation of MODIS and GLASS leaf area index (LAI) products over mountain areas: A case study in southwestern China

The validation study of leaf area index (LAI) products over rugged surfaces not only gives additional insights into data quality of LAI products, but deepens understanding of uncertainties regarding land surface process models depended on LAI data over complex terrain. This study evaluated the performance of MODIS and GLASS LAI products using the intercomparison and direct validation methods over southwestern China. The spatio-temporal consistencies, such as the spatial distributions of LAI products and their statistical relationship as a function of topographic indices, time, and vegetation types, respectively, were investigated through intercomparison between MODIS and GLASS products during the period 2011–2013. The accuracies and change ranges of these two products were evaluated against available LAI reference maps over 10 sampling regions which standed for typical vegetation types and topographic gradients in southwestern China.The results show that GLASS LAI exhibits higher percentage of good quality data (i.e. successful retrievals) and smoother temporal profiles than MODIS LAI. The percentage of successful retrievals for MODIS and GLASS is vulnerable to topographic indices, especially to relief amplitude. Besides, the two products do not capture seasonal dynamics of crop, especially in spring over heterogeneously hilly regions. The yearly mean LAI differences between MODIS and GLASS are within ±0.5 for 64.70% of the total retrieval pixels over southwestern China. The spatial distribution of mean differences and temporal profiles of these two products are inclined to be dominated by vegetation types other than topographic indices. The spatial and temporal consistency of these two products is good over most area of grasses/cereal crops; however, it is poor for evergreen broadleaf forest. MODIS presents more reliable change range of LAI than GLASS through comparison with fine resolution reference maps over most of sampling regions. The accuracies of direct validation are obtained for GLASS LAI (r = 0.35, RMSE = 1.72, mean bias = −0.71) and MODIS LAI (r = 0.49, RMSE = 1.75, mean bias = −0.67). GLASS performs similarly to MODIS, but may be marginally inferior to MODIS based on our direct validation results. The validation experience demonstrates the necessity and importance of topographic consideration for LAI estimation over mountain areas. Considerable attention will be paid to the improvements of surface reflectance, retrieval algorithm and land cover types so as to enhance the quality of LAI products in topographically complex terrain.     

Assessing the accuracy of hyperspectral and multispectral satellite imagery for categorical and Quantitative mapping of salinity stress in sugarcane fields

This study evaluates the feasibility of hyperspectral and multispectral satellite imagery for categorical and quantitative mapping of salinity stress in sugarcane fields located in the southwest of Iran. For this purpose a Hyperion image acquired on September 2, 2010 and a Landsat7 ETM+ image acquired on September 7, 2010 were used as hyperspectral and multispectral satellite imagery. Field data including soil salinity in the sugarcane root zone was collected at 191 locations in 25 fields during September 2010. In the first section of the paper, based on the yield potential of sugarcane as influenced by different soil salinity levels provided by FAO, soil salinity was classified into three classes, low salinity (1.7–3.4 dS/m), moderate salinity (3.5–5.9 dS/m) and high salinity (6–9.5) by applying different classification methods including Support Vector Machine (SVM), Spectral Angle Mapper (SAM), Minimum Distance (MD) and Maximum Likelihood (ML) on Hyperion and Landsat images. In the second part of the paper the performance of nine vegetation indices (eight indices from literature and a new developed index in this study) extracted from Hyperion and Landsat data was evaluated for quantitative mapping of salinity stress. The experimental results indicated that for categorical classification of salinity stress, Landsat data resulted in a higher overall accuracy (OA) and Kappa coefficient (KC) than Hyperion, of which the MD classifier using all bands or PCA (1–5) as an input performed best with an overall accuracy and kappa coefficient of 84.84% and 0.77 respectively. Vice versa for the quantitative estimation of salinity stress, Hyperion outperformed Landsat. In this case, the salinity and water stress index (SWSI) has the best prediction of salinity stress with an R² of 0.68 and RMSE of 1.15 dS/m for Hyperion followed by Landsat data with an R² and RMSE of 0.56 and 1.75 dS/m respectively. It was concluded that categorical mapping of salinity stress is the best option for monitoring agricultural fields and for this purpose Landsat data are most suitable.