Combined use of LiDAR data and multispectral earth observation imagery for wetland habitat mapping

Although wetlands play a key role in controlling flooding and nonpoint source pollution, sequestering carbon and providing an abundance of ecological services, the inventory and characterization of wetland habitats are most often limited to small areas. This explains why the understanding of their ecological functioning is still insufficient for a reliable functional assessment on areas larger than a few hectares. While LiDAR data and multispectral Earth Observation (EO) images are often used separately to map wetland habitats, their combined use is currently being assessed for different habitat types. The aim of this study is to evaluate the combination of multispectral and multiseasonal imagery and LiDAR data to precisely map the distribution of wetland habitats. The image classification was performed combining an object-based approach and decision-tree modeling. Four multispectral images with high (SPOT-5) and very high spatial resolution (Quickbird, KOMPSAT-2, aerial photographs) were classified separately. Another classification was then applied integrating summer and winter multispectral image data and three layers derived from LiDAR data: vegetation height, microtopography and intensity return. The comparison of classification results shows that some habitats are better identified on the winter image and others on the summer image (overall accuracies = 58.5 and 57.6%). They also point out that classification accuracy is highly improved (overall accuracy = 86.5%) when combining LiDAR data and multispectral images. Moreover, this study highlights the advantage of integrating vegetation height, microtopography and intensity parameters in the classification process. This article demonstrates that information provided by the synergetic use of multispectral images and LiDAR data can help in wetland functional assessment     

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.     

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.     

Texture augmented detection of macrophyte species using decision trees

Image classification using multispectral sensors has shown good performance in detecting macrophytes at the species level. However, species level classification often does not utilize the texture information provided by high resolution images. This study investigated whether image texture provides useful vector(s) for the discrimination of monospecific stands of three floating macrophyte species in Quickbird imagery of the South Nation River. Semivariograms indicated that window sizes of 5 × 5 and 13 × 13 pixels were the most appropriate spatial scales for calculation of the grey level co-occurrence matrix and subsequent texture attributes from the multispectral and panchromatic bands. Of the 214 investigated vectors (13 Haralick texture attributes * 15 bands + 9 spectral bands + 10 transformations/indices), feature selection determined which combination of spectral and textural vectors had the greatest class separability based on the Mann–Whitney U-test and Jefferies–Matusita distance. While multispectral red and near infrared (NIR) performed satisfactorily, the addition of panchromatic-dissimilarity slightly improved class separability and the accuracy of a decision tree classifier (Kappa: red/NIR/panchromatic-dissimilarity – 93.2% versus red/NIR – 90.4%). Class separability improved by incorporating a second texture attribute, but resulted in a decrease in classification accuracy. The results suggest that incorporating image texture may be beneficial for separating stands with high spatial heterogeneity. However, the benefits may be limited and must be weighed against the increased complexity of the classifier.     

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.     

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.     

CBERS-02B卫星遥感影像的区域网平差

针对中巴资源一号卫星(CBERS-02B)卫星遥感影像姿态角误差较大的特点,提出了利用区域网平差方法提高其对地目标定位精度的策略和具体计算过程。首先对参与平差的每景影像选取4个地面控制点进行影像姿态角常差检校,然后采用与地形无关方案解求各自的RPC参数,最后选取带仿射变换项的有理函数模型(RFM)进行多重覆盖影像的区域网平差。对两个地区的0级CBERS-02B单条CCD立体影像对的区域网平差试验表明,对地目标定位的平面和高程精度均达到了±3个像元的水平,且高程精度明显优于平面精度。相对于常规的卫星遥感影像区域网平差方法,平面和高程精度均有明显提升,几乎达到国外同等高分辨率卫星遥感影像的几何定位精度。这说明中国卫星遥感影像亦具有较好的几何定位潜力,在区域网平差之前进行系统误差预改正是必要和可行的。     

Exterior orientation of CBERS-2B imagery using multi-feature control and orbital data

The major contribution of this paper relates to the practical advantages of combining Ground Control Points (GCPs), Ground Control Lines (GCLs) and orbital data to estimate the exterior orientation parameters of images collected by CBERS-2B (China–Brazil Earth Resources Satellite) HRC (High-resolution Camera) and CCD (High-resolution CCD Camera) sensors. Although the CBERS-2B is no longer operational, its images are still being used in Brazil, and the next generations of the CBERS satellite will have sensors with similar technical features, which motivates the study presented in this paper. The mathematical models that relate the object and image spaces are based on collinearity (for points) and coplanarity (for lines) conditions. These models were created in an in-house developed software package called TMS (Triangulation with Multiple Sensors) with multi-feature control (GCPs and GCLs). Experiments on a block of four CBERS-2B HRC images and on one CBERS-2B CCD image were performed using both models. It was observed that the combination of GCPs and GCLs provided better bundle block adjustment results than conventional bundle adjustment using only GCPs. The results also demonstrate the advantages of using primarily orbital data when the number of control entities is reduced.     

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.     

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.     

Testing estimation of water surface in Italian rice district from MODIS satellite data

Recent changes in rice crop management within Northern Italy rice district led to a reduction of seeding in flooding condition, which may have an impact on reservoir water management and on the animal and plant communities that depend on the flooded paddies. Therefore, monitoring and quantifying the spatial and temporal variability of water presence in paddy fields is becoming important. In this study we present a method to estimate dynamics of presence of standing water (i.e. fraction of flooded area) in rice fields using MODIS data. First, we produced high resolution water presence maps from Landsat by thresholding the Normalised Difference Flood Index (NDFI) made: we made it by comparing five Landsat 8 images with field-obtained information about rice field status and water presence. Using these data we developed an empirical model to estimate the flooding fraction of each MODIS cell. Finally we validated the MODIS-based flooding maps with both Landsat and ground information. Results showed a good predictability of water surface from Landsat (OA = 92%) and a robust usability of MODIS data to predict water fraction (R² = 0.73, EF = 0.57, RMSE = 0.13 at 1 × 1 km resolution). Analysis showed that the predictive ability of the model decreases with the greening up of rice, so we used NDVI to automatically discriminate estimations for inaccurate cells in order to provide the water maps with a reliability flag. Results demonstrate that it is possible to monitor water dynamics in rice paddies using moderate resolution multispectral satellite data. The achievement is a proof of concept for the analysis of MODIS archives to investigate irrigation dynamics in the last 15 years to retrieve information for ecological and hydrological studies.     

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.     

CBERS-02星图像几何纠正方法试验研究

为了研究中巴地球资源卫星CCD多光谱数据规模化应用的可行性,扩大基础地理信息更新的影像资料来源,以CBERS-02星（简称02星）CCD多光谱数据的标准景为例,对2种几何纠正模型进行了对比研究,分析了3种不同基础控制资料对纠正精度的影响,并对精度评价方法进行了试验研究。在此基础上,总结形成了针对02星CCD多光谱数据1∶10万比例尺影像最为省时有效、精度最佳的几何纠正技术方法。     

Augmenting forest inventory attributes with geometric optical modelling in support of regional susceptibility assessments to bark beetle infestations

Assessment of the susceptibility of forests to mountain pine beetle (Dendroctonus ponderosae Hopkins) infestation is based upon an understanding of the characteristics that predispose the stands to attack. These assessments are typically derived from conventional forest inventory data; however, this information often represents only managed forest areas. It does not cover areas such as forest parks or conservation regions and is often not regularly updated resulting in an inability to assess forest susceptibility. To address these shortcomings, we demonstrate how a geometric optical model (GOM) can be applied to Landsat-5 Thematic Mapper (TM) imagery (30 m spatial resolution) to estimate stand-level susceptibility to mountain pine beetle attack. Spectral mixture analysis was used to determine the proportion of sunlit canopy and background, and shadow of each Landsat pixel enabling per pixel estimates of attributes required for model inversion. Stand structural attributes were then derived from inversion of the geometric optical model and used as basis for susceptibility mapping. Mean stand density estimated by the geometric optical model was 2753 (standard deviation ± 308) stems per hectare and mean horizontal crown radius was 2.09 (standard deviation ± 0.11) metres. When compared to equivalent forest inventory attributes, model predictions of stems per hectare and crown radius were shown to be reasonably estimated using a Kruskal–Wallis ANOVA (p < 0.001). These predictions were then used to create a large area map that provided an assessment of the forest area susceptible to mountain pine beetle damage.     

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.     

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.     

The effect of atmospheric and topographic correction methods on land cover classification accuracy

Mapping of vegetation in mountain areas based on remote sensing is obstructed by atmospheric and topographic distortions. A variety of atmospheric and topographic correction methods has been proposed to minimize atmospheric and topographic effects and should in principle lead to a better land cover classification. Only a limited number of atmospheric and topographic combinations has been tested and the effect on class accuracy and on different illumination conditions is not yet researched extensively. The purpose of this study was to evaluate the effect of coupled correction methods on land cover classification accuracy. Therefore, all combinations of three atmospheric (no atmospheric correction, dark object subtraction and correction based on transmittance functions) and five topographic corrections (no topographic correction, band ratioing, cosine correction, pixel-based Minnaert and pixel-based C-correction) were applied on two acquisitions (2009 and 2010) of a Landsat image in the Romanian Carpathian mountains. The accuracies of the fifteen resulting land cover maps were evaluated statistically based on two validation sets: a random validation set and a validation subset containing pixels present in the difference area between the uncorrected classification and one of the fourteen corrected classifications. New insights into the differences in classification accuracy were obtained. First, results showed that all corrected images resulted in higher overall classification accuracies than the uncorrected images. The highest accuracy for the full validation set was achieved after combination of an atmospheric correction based on transmittance functions and a pixel-based Minnaert topographic correction. Secondly, class accuracies of especially the coniferous and mixed forest classes were enhanced after correction. There was only a minor improvement for the other land cover classes (broadleaved forest, bare soil, grass and water). This was explained by the position of different land cover types in the landscape. Finally, coupled correction methods showed most efficient on weakly illuminated slopes. After correction, accuracies in the low illumination zone (cos β ≤ 0.65) were improved more than in the moderate and high illumination zones. Considering all results, best overall classification results were achieved after combination of the transmittance function correction with pixel-based Minnaert or pixel-based C-topographic correction. Furthermore, results of this bi-temporal study indicated that the topographic component had a higher influence on classification accuracy than the atmospheric component and that it is worthwhile to invest in both atmospheric and topographic corrections in a multi-temporal study.     

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.     

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.     

Determination of the crop row orientations from Formosat-2 multi-temporal and panchromatic images

This paper presents a technique developed for the retrieval of the orientation of crop rows, over anthropic lands dedicated to agriculture in order to further improve estimate of crop production and soil erosion management. Five crop types are considered: wheat, barley, rapeseed, sunflower, corn and hemp. The study is part of the multi-sensor crop-monitoring experiment, conducted in 2010 throughout the agricultural season (MCM’10) over an area located in southwestern France, near Toulouse. The proposed methodology is based on the use of satellite images acquired by Formosat-2, at high spatial resolution in panchromatic and multispectral modes (with spatial resolution of 2 and 8 m, respectively). Orientations are derived and evaluated for each image and for each plot, using directional spatial filters (45° and 135°) and mathematical morphology algorithms. “Single-date” and “multi-temporal” approaches are considered. The single-date analyses confirm the good performances of the proposed method, but emphasize the limitation of the approach for estimating the crop row orientation over the whole landscape with only one date. The multi-date analyses allow (1) determining the most suitable agricultural period for the detection of the row orientations, and (2) extending the estimation to the entire footprint of the study area. For the winter crops (wheat, barley and rapeseed), best results are obtained with images acquired just after harvest, when surfaces are covered by stubbles or during the period of deep tillage (0.27 > R² > 0.99 and 7.15° > RMSE > 43.02°). For the summer crops (sunflower, corn and hemp), results are strongly crop and date dependents (0 > R² > 0.96, 10.22° > RMSE > 80°), with a well-marked impact of flowering, irrigation equipment and/or maximum crop development. Last, the extent of the method to the whole studied zone allows mapping 90% of the crop row orientations (more than 45,000 ha) with an error inferior to 40°, associated to a confidence index ranging from 1 to 5 for each agricultural plot.