Landsat-7 ETM+ radiometric normalization comparison for northern mapping applications

Relative radiometric normalization has long been performed to generate consistency among individual Landsat scenes for production of composites containing multiple scenes. Normalization methods have relied on matching identical and assumed invariant features in both images of an overlapping pair, or on invariant targets that are not necessarily the same features. Problems with overlap normalization methods include sensitivity to outliers in overlap data caused by atmospheric or land cover change between scenes, which can lead to radiometric error propagation across a mosaic caused by a normalized scene becoming a reference for the subsequent scene entered into the mosaic. Solutions to such problems include interactive outlier removal to generate a normalization function using a ‘no change’ data set and methods that are robust against outliers to automatically generate normalization functions with minimal user input. This paper compares two normalization methods that use a robust regression technique called Theil-Sen with an established overlap normalization method. The first method uses Theil-Sen regression to generate a normalization function between overlap regions, while the second uses Theil-Sen to normalize to coarse-resolution composite reflectance data from the SPOT VEGETATION (VGT) sensor. The results of the normalizations were evaluated in two ways: (1) using statistics generated between overlap regions; and (2) separately using coarse-resolution data as a reference. Both overlap normalization methods performed almost identically; however, Theil-Sen was faster and easier to implement than its traditional counterpart due to its insensitivity to outliers and capability for full automation. While overlap and coarse-resolution normalizations each outperformed the other when evaluated against its calibration set, error propagation caused by outliers in overlap samples was avoided in the normalization to coarse-resolution imagery. Advantages offered by normalization to coarse-resolution data using robust regression, including full automation, make this method particularly attractive for generation of large area mosaics containing 100 Landsat scenes or more.     

Comparing the spectral signal of burned surfaces between Landsat 7 ETM+ and Landsat 8 OLI sensors

There is a long history of the use of Landsat data in burned land mapping mainly due to certain characteristics of the Landsat imagery including the spatial, spectral, and temporal data resolution, the low cost (Landsat data are now freely available), and the existence of an almost 35-year historical archive (excluding Landsat 1–3). Landsat 8 (Operational Land Imager (OLI) and Thermal Infrared Sensor (TIRS)) was launched on 11 February 2013 and it captures data in three new bands along with two additional thermal bands. However, is the spectral signal of burned surfaces in satellite remote-sensing data of Landsat series consistent and robust enough to allow the successful application of the techniques developed so far for Landsat 8? In this article, we compare the spectral signal of burned surfaces between Landsat 7 Enhanced Thematic Mapper Plus (ETM+) and Landsat 8 OLI sensors using five case studies that correspond to five large fire events in different biophysical environments in Greece, for which both Landsat 7 ETM+ and Landsat 8 OLI data were available. From the comparative analysis using histogram data plots of burned (post-fire image) and vegetated (pre-fire image) areas, spectral signature plots and separability indices of certain land-cover types, estimated using the same sampling areas over both satellite images, a general consistency was observed between the two sensors. Slight differences between the sensors were attributed to differences in the acquisition dates and were related to the type of vegetation rather than the sensors used to record the satellite images. Neither sensor provided improved discrimination over the other.     

Identifying ring geomorphic features in Lebanon using ASTER and Landsat 7 ETM+ images

Satellite images enable us to identify several geomorphologic features on the Earth's surface that could not be easily recognized on ground surface or by using conventional methods. This is mainly attributed to the optical advantages of remote sensing techniques. Thus, suspicious geomorphologic signatures can be observed on the terrain surface. These are mostly geologic‐controlled. Linear aspects are given most attention in many geological studies to reflect subsurface, hidden structures. In this study, using ENVI‐4.3 and ERDAS Imagine‐9.1 software to analyse ASTER and Landsat 7 ETM+ images of Lebanon, we exposed a miscellany of geomorphic ring structures on different geologic formations of the Lebanese terrain. The diameter of the recognized ring structures ranges from a few hundreds of metres to several kilometres. Preliminary field surveys were carried out on some of these structures in order to identify their origin. All the recognized structures are circular and semi‐rounded in their geomorphology; however, they are attributed either to subsurface crust processes or to the impact of falling meteoritic from outer space. The former were indicated through rock deformations and existence of intrusive bodies. Meteorite impact was induced from the glassy and metallic materials prevailing within the ring structure area, as well as from the interference of some rings.     

Assessment of fire severity and species diversity in the southern Appalachians using Landsat TM and ETM+ imagery

Relatively little is known about the disturbance ecology of large wildfires in the southern Appalachians. The occurrence of a 4000-ha wildfire in the Linville Gorge Wilderness area in western North Carolina has provided a rare opportunity to study a large fire with a range of severities. The objectives of this study were to 1) assess the potential for using multi-temporal Landsat imagery to map fire severity in the southern Appalachians, 2) examine the influences of topography and forest community type on the spatial pattern of fire severity; and 3) examine the relationship between predicted fire severity and changes in species richness. A non-linear regression equation predicted a field-based composite burn index (CBI) as a function of change in the Normalized Burn Ratio (dNBR) with an R² of 0.71. Fire severity was highest on drier landforms located on upper hillslopes, ridges, and on southwest aspects, and was higher in pine communities than in other forest types. Predicted CBI was positively correlated with changes in species richness and with the post-fire cover of pine seedlings (Pinus virginiana, P. rigida, and P. pungens), suggesting that burn severity maps can be used to predict community-level fire effects across large landscapes. Despite the relatively large size of this fire for the southern Appalachians, severity was strongly linked to topographic variability and pre-fire vegetation, and spatial variation in fire severity was correlated with changes in species richness. Thus, the Linville Gorge fire appears to have generally reinforced the ecological constraints imposed by underlying environmental gradients.     

A consistency analysis of surface reflectance and leaf area index retrieval from overlapping clear-sky Landsat ETM+ imagery

In this study, the consistency of systematic retrievals of surface reflectance and leaf area index was assessed using overlap regions in adjacent Landsat Enhanced Thematic Mapper-Plus (ETM+) scenes. Adjacent scenes were acquired within 7-25 days apart to minimize variations in the land surface reflectance between acquisition dates. Each Landsat ETM+ scene was independently geo-referenced and atmospherically corrected using a variety of standard approaches. Leaf area index (LAI) models were then applied to the surface reflectance data and the difference in LAI between overlapping scenes was evaluated. The results from this analysis show that systematic LAI retrieval from Landsat ETM+ imagery using a baseline atmospheric correction approach that assumes a constant aerosol optical depth equal to 0.06 is consistent to within ±0.61 LAI units. The average absolute difference in LAI retrieval over all 10 image pairs was 26% for a mean LAI of 2.05 and the maximum absolute difference over any one pair was 61% for a mean LAI of 1.13. When no atmospheric correction was performed on the data, the consistency in LAI retrieval was improved by 1%. When a scene-based dense, dark vegetation atmospheric correction algorithm was used, the LAI retrieval differences increased to 28% for a mean LAI of 2.32. This implies that a scene-based atmospheric correction procedure may improve the absolute accuracy of LAI retrieval without having a major impact on retrieval consistency. Such consistency trials provide insight into the current limits concerning surface reflectance and LAI retrieval from fine spatial resolution remote sensing imagery with respect to the variability in clear-sky atmospheric conditions.     

Recovering missing pixels for Landsat ETM + SLC-off imagery using HJ-1A /1B as auxiliary data

To improve the usability of Enhanced Thematic Mapper Plus (ETM+) scan line corrector (SLC)-off data, this article proposes using HJ-1A/1B imagery as auxiliary (i.e. reference) data to recover the SLC-off ETM+ data. The least-median-of-squares (LMedS) method is newly proposed to recover missing pixels of Landsat 7 by removing the variant or abnormal digital number values. In particular, for the visible and near-infrared bands, using HJ-1A/1B for recovery has three clear advantages: the same spatial resolution, similar spectral resolution, and approximate temporal resolution. The experiments show that all of the reference-recovery methods are better than the non-reference-recovery method. The results of using of auxiliary data in reference-recovery methods, from best to worst, are Landsat 8, HJ-1A/1B, and Landsat 7. However, for recovering missing pixels, HJ-1A/1B is superior to the ETM+ auxiliary data due to the shorter time interval in Landsat 7 (a few hours). Hence, HJ-1A/1B should be considered a useful auxiliary data to recover ETM+ SLC-off imagery data.     

Landsat ETM+ imaging for mineral potential mapping: application to Avaj area,Qazvin, Iran

The Avaj area is located in the northwestern part of Iran, on the boundary between the Orumieh-Dokhtar and Sanandaj-Sirjan geological zones. It corresponds to two different geological sub-zones:Abe-Garm, to the north, which is related to the volcanic belt of Orumieh-Dokhtar; and the Razan sub-zone, to the south, which is related to the Sanandaj-Sirjan zone. Lithological aspects of different rock units indicate the presence of four distinct volcanic groups of andesite, basalt, tuff, and dacite. The present research utilizes various Landsat Enhanced Thematic Mapper Plus (ETM+) image-processing techniques, including false colour composite images, colour composite ratio images, and standard principal component (PC) image analysis on six bands and colour composite selective PC images, filtering and supervised classification. Nonetheless, colour composite selective images may be the most reliable method for enhancement of areas with hydrothermal alteration. All the techniques used clearly show the pervasive alteration of kaolinite-argillite in the area. These alterations are mainly related to the andesite, and tuffs have the same trend as fault directions.     

A comparison of empirical and neural network approaches for estimating corn and soybean leaf area index from Landsat ETM+ imagery

Plant foliage density expressed as leaf area index (LAI) is used in many ecological, meteorological, and agronomic models, and as a means of quantifying crop spatial variability for precision farming. LAI retrieval using spectral vegetation indices (SVI) from optical remotely sensed data usually requires site-specific calibration values from the surface or the use of within-scene image information without surface calibrations to invert radiative transfer models. An evaluation of LAI retrieval methods was conducted using (1) empirical methods employing the normalized difference vegetation index (NDVI) and a new SVI that uses green wavelength reflectance, (2) a scaled NDVI approach that uses no calibration measurements, and (3) a hybrid approach that uses a neural network (NN) and a radiative transfer model without site-specific calibration measurements. While research has shown that under a variety of conditions NDVI is not optimal for LAI retrieval, its continued use for remote sensing applications and in analysis seeking to develop improved parameter retrieval algorithms based on NDVI suggests its value as a “benchmark” or standard against which other methods can be compared. Landsat-7 ETM+ data for July 1 and July 8 from the Soil Moisture EXperiment 2002 (SMEX02) field campaign in the Walnut Creek watershed south of Ames, IA, were used for the analysis. Sun photometer data collected from a site within the watershed were used to atmospherically correct the imagery to surface reflectance. LAI validation measurements of corn and soybeans were collected close to the dates of the Landsat-7 overpasses. Comparable results were obtained with the empirical SVI methods and the scaled SVI method within each date. The hybrid method, although promising, did not account for as much of the variability as the SVI methods. Higher atmospheric optical depths for July 8 leading to surface reflectance errors are believed to have resulted in overall poorer performance for this date. Use of SVIs employing green wavelengths, improved method for the definition of image minimum and maximum clusters used by the scaled NDVI method, and further development of a soil reflectance index used by the hybrid NN approach are warranted. More importantly, the results demonstrate that reasonable LAI estimates are possible using optical remote sensing methods without in situ, site-specific calibration measurements.     

Comparison of AVIRIS and Landsat ETM+ detection capabilities for burn severity

Our study compares data on burn severity collected from multi-temporal Airborne Visible and Infrared Imaging Spectrometer (AVIRIS) with similar data from the Enhanced Thematic Mapper Plus (ETM+) using the differenced Normalized Burn Ratio (dNBR). Two AVIRIS and ETM+ data acquisitions recorded surface conditions immediately before the Hoover Fire began to spread rapidly and again the following year. Data were validated with 63 field plots using the Composite Burn Index (CBI). The relationship between spectral channels and burn severity was examined by comparing pre- and post-fire datasets. Based on the high burn severity comparison, AVIRIS channels 47 and 60 at wavelengths of 788 and 913 nm showed the greatest negative response to fire. Post-fire reflectance values decreased the most on average at those wavelengths, while channel 210 at 2370 nm showed the greatest positive response on average. Fire increased reflectance the most at that wavelength over the entire measured spectral range. Furthermore, channel 210 at 2370 nm exhibited the greatest variation in spectral response, suggesting potentially high information content for fire severity. Based on general remote sensing principles and the logic of variable spectral responses to fire, dNBR from both sensors should produce useful results in quantifying burn severity. The results verify the band-response relationships to burn severity as seen with ETM+ data and confirm the relationships by way of a distinctly different sensor system.     

Comparison of two types of forest disturbance using multitemporal Landsat TM/ETM+ imagery and field vegetation data

Forest disturbances influence many landscape processes, including changes in microclimate, hydrology, and soil erosion. We analyzed the spectral response and temporal progress of two types of disturbances of spruce forest (bark beetle outbreak and clear-cuts) in the central part of Šumava Mountains at the border between the Czech Republic and Germany, Central Europe. The bark beetle (Ips typographus [L.]) outbreak in this region in the last 20 years resulted in regional-scale spruce forest decay. Clear-cutting was done here to prevent further bark-beetle propagation in the buffer zones.The aim of the study is to identify the differences in spectral response between the two types of forest disturbances and their temporal dynamics. General trends were analyzed throughout the study area, with sampled disturbance areas selected to assess the relationship between field vegetation data and their spectral response. Thirteen Landsat TM/ETM+ scenes from 1985 to 2007 were used for the assessment. The following spectral indices were estimated: NDMI, Tasseled Cap (Brightness, Greenness, Wetness), DI, and DI′. The DI′, Wetness, and Brightness indices show the highest sensitivity to forest disturbance for both disturbance types (clear-cuts and bark beetle outbreak). The multitemporal analysis distinguished three different stages of development. The highest spectral differences between the clear-cuts and the bark beetle disturbances were found in the period between 1996 and 2004 with increased levels of forest disturbance (repeated measures ANOVA, Scheffé post hoc test; p ≤ 0.05). Clear-cut disturbance resulted in significantly higher spectral differences from the original forest and occurred as a more discrete event in comparison to bark beetle outbreak.     

An automated scheme for glacial lake dynamics mapping using Landsat imagery and digital elevation models: a case study in the Himalayas

Glacial lakes in alpine regions are sensitive to climate change. Mapping and monitoring these lakes would improve our understanding of regional climate change and glacier-related hazards. However, glacial lake mapping over large areas using remote sensing remains a challenge because of various disturbing factors in glacial and periglacial environments. This article presents an automated mapping algorithm based on hierarchical image segmentation and terrain analysis to delineate glacial lake extents. In this algorithm, each glacial lake is delineated with a local segmentation value, and the topographic features derived from digital elevation models (DEMs) are also used to separate mountain shadows from glacial lakes. About 100 scenes of Landsat Thematic Mapper/Enhanced Thematic Mapper Plus (TM/ETM+) images from circa 1990, circa 2000 and 2009 were used to map the glacial lakes and their changes over the entire Himalayas. The results show that the algorithm can map the glacial lakes effectively and efficiently. Mountain shadows or melting glaciers can be differentiated from glacial lakes automatically, and those lakes with mountain shadows can also be identified. Area changes of more than 1000 glacial lakes show that the glacial lakes in the Himalayas have experienced mixed directions of change, while the overall lake areas are expanding at an accelerated rate in the past two decades, indicating great changes to the glacial lakes in the Himalayas.     

Using lidar and effective LAI data to evaluate IKONOS and Landsat 7 ETM+ vegetation cover estimates in a ponderosa pine forest

Structural and functional analyses of ecosystems benefit when high accuracy vegetation coverages can be derived over large areas. In this study, we utilize IKONOS, Landsat 7 ETM+, and airborne scanning light detection and ranging (lidar) to quantify coniferous forest and understory grass coverages in a ponderosa pine (Pinus ponderosa) dominated ecosystem in the Black Hills of South Dakota. Linear spectral mixture analyses of IKONOS and ETM+ data were used to isolate spectral endmembers (bare soil, understory grass, and tree/shade) and calculate their subpixel fractional coverages. We then compared these endmember cover estimates to similar cover estimates derived from lidar data and field measures. The IKONOS-derived tree/shade fraction was significantly correlated with the field-measured canopy effective leaf area index (LAI_e) (r²=0.55, p<0.001) and with the lidar-derived estimate of tree occurrence (r²=0.79, p<0.001). The enhanced vegetation index (EVI) calculated from IKONOS imagery showed a negative correlation with the field measured tree canopy effective LAI and lidar tree cover response (r²=0.30, r=−0.55 and r²=0.41, r=−0.64, respectively; p<0.001) and further analyses indicate a strong linear relationship between EVI and the IKONOS-derived grass fraction (r²=0.99, p<0.001). We also found that using EVI resulted in better agreement with the subpixel vegetation fractions in this ecosystem than using normalized difference of vegetation index (NDVI). Coarsening the IKONOS data to 30 m resolution imagery revealed a stronger relationship with lidar tree measures (r²=0.77, p<0.001) than at 4 m resolution (r²=0.58, p<0.001). Unmixed tree/shade fractions derived from 30 m resolution ETM+ imagery also showed a significant correlation with the lidar data (r²=0.66, p<0.001). These results demonstrate the power of using high resolution lidar data to validate spectral unmixing results of satellite imagery, and indicate that IKONOS data and Landsat 7 ETM+ data both can serve to make the important distinction between tree/shade coverage and exposed understory grass coverage during peak summertime greenness in a ponderosa pine forest ecosystem.     

基于NSCT和LSSVM的Landsat ETM+图像土地覆盖分类

提出一种新的基于非下采样Contourlet变换（NSCT）和最小二乘支持向量机（LSSVM）的遥感图像土地覆盖分类方法。该方法动态选择最优的多光谱图像的波段特征进行组合,基于NSCT和IHS对多光谱图像和全色图像进行融合,增强多光谱图像的空间分辨率,基于LSSVM对融合图像进行分类。实验结果表明,提出的方法在保留多光谱图像光谱信息的同时,增强了图像的空间细节表现能力,提供更加可靠的地物分类特征,提高了分类精度,并且优于传统的基于最小距离法、最大似然法、贝叶斯分类法和BPNN分类法的遥感图像分类方法,该方法是有效可行的。     

Comparison of data gap-filling methods for Landsat ETM+ SLC-off imagery for monitoring forest degradation in a semi-deciduous tropical forest in Mexico

A number of methods to overcome the 2003 failure of the Landsat 7 Enhanced Thematic Mapper (ETM+) scan-line corrector (SLC) are compared in this article in a forest-monitoring application in the Yucatan Peninsula, Mexico. The objective of this comparison is to determine the best approach to accomplish SLC-off image gap-filling for the particular landscape in this region, and thereby provide continuity in the Landsat data sensor archive for forest-monitoring purposes. Four methods were tested: (1) local linear histogram matching (LLHM); (2) neighbourhood similar pixel interpolator (NSPI); (3) geostatistical neighbourhood similar pixel interpolator (GNSPI); and (4) weighted linear regression (WLR). All methods generated reasonable SLC-off gap-filling data that were visually consistent and could be employed in subsequent digital image analysis. Overall accuracy, kappa coefficients (κ), and quantity and allocation disagreement indices were used to evaluate unsupervised Iterative Self-Organizing Data Analysis (ISODATA) land-cover classification maps. In addition, Pearson correlation coefficients (r) and root mean squares of the error (RMSEs) were employed for estimates agreement with fractional land cover. The best results visually (overall accuracy > 85%, κ < 9%, quantity disagreement index < 5.5%, and allocation disagreement index < 12.5%) and statistically (r > 0.84 and RMSE < 7%) were obtained from the GNSPI method. These results suggest that the GNSPI method is suitable for routine use in reconstructing the imagery stack of Landsat ETM+ SLC-off gap-filled data for use in forest-monitoring applications in this type of heterogeneous landscape.     

An improved strategy for regression of biophysical variables and Landsat ETM+ data

Empirical models are important tools for relating field-measured biophysical variables to remote sensing data. Regression analysis has been a popular empirical method of linking these two types of data to provide continuous estimates for variables such as biomass, percent woody canopy cover, and leaf area index (LAI). Traditional methods of regression are not sufficient when resulting biophysical surfaces derived from remote sensing are subsequently used to drive ecosystem process models. Most regression analyses in remote sensing rely on a single spectral vegetation index (SVI) based on red and near-infrared reflectance from a single date of imagery. There are compelling reasons for utilizing greater spectral dimensionality, and for including SVIs from multiple dates in a regression analysis. Moreover, when including multiple SVIs and/or dates, it is useful to integrate these into a single index for regression modeling. Selection of an appropriate regression model, use of multiple SVIs from multiple dates of imagery as predictor variables, and employment of canonical correlation analysis (CCA) to integrate these multiple indices into a single index represent a significant strategic improvement over existing uses of regression analysis in remote sensing.To demonstrate this improved strategy, we compared three different types of regression models to predict LAI for an agro-ecosystem and live tree canopy cover for a needleleaf evergreen boreal forest: traditional (Y on X) ordinary least squares (OLS) regression, inverse (X on Y) OLS regression, and an orthogonal regression method called reduced major axis (RMA). Each model incorporated multiple SVIs from multiple dates and CCA was used to integrate these. For a given dataset, the three regression-modeling approaches produced identical coefficients of determination and intercepts, but different slopes, giving rise to divergent predictive characteristics. The traditional approach yielded the lowest root mean square error (RMSE), but the variance in the predictions was lower than the variance in the observed dataset. The inverse method had the highest RMSE and the variance was inflated relative to the variance of the observed dataset. RMA provided an intermediate set of predictions in terms of the RMSE, and the variance in the observations was preserved in the predictions. These results are predictable from regression theory, but that theory has been essentially ignored within the discipline of remote sensing.     

长江三峡库区Landsat7 ETM+数据的处理方法探讨

Landsat7 ETM+是陆地卫星系列的最新数据,其突出优势在于提高了多光谱波段的空间分辨率,增加了一个分辨率为15 m的全色波段,丰富了数据的信息量。选取长江三峡库区的两景ETM+原始数据,基于不同的地学遥感应用要求,试验性地研究了Landsat7 ETM+新型数据的数字处理方法,并比较了这些处理方法的应用效果。结果表明,选择适当的数字方法对新型Landsat7ETM+遥感数据进行处理,能够快速、有效地获得同时具备较高的光谱分辨率和空间分辨率的遥感数字图像,显示出较好的应用优势。     

A simple and effective method for filling gaps in Landsat ETM+ SLC-off images

The scan-line corrector (SLC) of the Landsat 7 Enhanced Thematic Mapper Plus (ETM+) sensor failed in 2003, resulting in about 22% of the pixels per scene not being scanned. The SLC failure has seriously limited the scientific applications of ETM+ data. While there have been a number of methods developed to fill in the data gaps, each method has shortcomings, especially for heterogeneous landscapes. Based on the assumption that the same-class neighboring pixels around the un-scanned pixels have similar spectral characteristics, and that these neighboring and un-scanned pixels exhibit similar patterns of spectral differences between dates, we developed a simple and effective method to interpolate the values of the pixels within the gaps. We refer to this method as the Neighborhood Similar Pixel Interpolator (NSPI). Simulated and actual SLC-off ETM+ images were used to assess the performance of the NSPI. Results indicate that NSPI can restore the value of un-scanned pixels very accurately, and that it works especially well in heterogeneous regions. In addition, it can work well even if there is a relatively long time interval or significant spectral changes between the input and target image. The filled images appear reasonably spatially continuous without obvious striping patterns. Supervised classification using the maximum likelihood algorithm was done on both gap-filled simulated SLC-off data and the original “gap free” data set, and it was found that classification results, including accuracies, were very comparable. This indicates that gap-filled products generated by NSPI will have relevance to the user community for various land cover applications. In addition, the simple principle and high computational efficiency of NSPI will enable processing large volumes of SLC-off ETM+ data.     

Finer resolution observation and monitoring of global land cover: first mapping results with Landsat TM and ETM+ data

We have produced the first 30 m resolution global land-cover maps using Landsat Thematic Mapper (TM) and Enhanced Thematic Mapper Plus (ETM+) data. We have classified over 6600 scenes of Landsat TM data after 2006, and over 2300 scenes of Landsat TM and ETM+ data before 2006, all selected from the green season. These images cover most of the world's land surface except Antarctica and Greenland. Most of these images came from the United States Geological Survey in level L1T (orthorectified). Four classifiers that were freely available were employed, including the conventional maximum likelihood classifier (MLC), J4.8 decision tree classifier, Random Forest (RF) classifier and support vector machine (SVM) classifier. A total of 91,433 training samples were collected by traversing each scene and finding the most representative and homogeneous samples. A total of 38,664 test samples were collected at preset, fixed locations based on a globally systematic unaligned sampling strategy. Two software tools, Global Analyst and Global Mapper developed by extending the functionality of Google Earth, were used in developing the training and test sample databases by referencing the Moderate Resolution Imaging Spectroradiometer enhanced vegetation index (MODIS EVI) time series for 2010 and high resolution images from Google Earth. A unique land-cover classification system was developed that can be crosswalked to the existing United Nations Food and Agriculture Organization (FAO) land-cover classification system as well as the International Geosphere-Biosphere Programme (IGBP) system. Using the four classification algorithms, we obtained the initial set of global land-cover maps. The SVM produced the highest overall classification accuracy (OCA) of 64.9% assessed with our test samples, with RF (59.8%), J4.8 (57.9%), and MLC (53.9%) ranked from the second to the fourth. We also estimated the OCAs using a subset of our test samples (8629) each of which represented a homogeneous area greater than 500 m?×?500 m. Using this subset, we found the OCA for the SVM to be 71.5%. As a consistent source for estimating the coverage of global land-cover types in the world, estimation from the test samples shows that only 6.90% of the world is planted for agricultural production. The total area of cropland is 11.51% if unplanted croplands are included. The forests, grasslands, and shrublands cover 28.35%, 13.37%, and 11.49% of the world, respectively. The impervious surface covers only 0.66% of the world. Inland waterbodies, barren lands, and snow and ice cover 3.56%, 16.51%, and 12.81% of the world, respectively.     

Monitoring on-orbit calibration stability of the Terra MODIS and Landsat 7 ETM+ sensors using pseudo-invariant test sites

The ability to detect and quantify changes in the Earth's environment depends on sensors that can provide calibrated, consistent measurements of the Earth's surface features through time. A critical step in this process is to put image data from different sensors onto a common radiometric scale. This work focuses on monitoring the long-term on-orbit calibration stability of the Terra Moderate Resolution Imaging Spectroradiometer (MODIS) and the Landsat 7 (L7) Enhanced Thematic Mapper Plus (ETM+) sensors using the Committee on Earth Observation Satellites (CEOS) reference standard pseudo-invariant test sites (Libya 4, Mauritania 1/2, Algeria 3, Libya 1, and Algeria 5). These sites have been frequently used as radiometric targets because of their relatively stable surface conditions temporally. This study was performed using all cloud-free calibrated images from the Terra MODIS and the L7 ETM+ sensors, acquired from launch to December 2008. Homogeneous regions of interest (ROI) were selected in the calibrated images and the mean target statistics were derived from sensor measurements in terms of top-of-atmosphere (TOA) reflectance. For each band pair, a set of fitted coefficients (slope and offset) is provided to monitor the long-term stability over very stable pseudo-invariant test sites. The average percent differences in intercept from the long-term trends obtained from the ETM + TOA reflectance estimates relative to the MODIS for all the CEOS reference standard test sites range from 2.5% to 15%. This gives an estimate of the collective differences due to the Relative Spectral Response (RSR) characteristics of each sensor, bi-directional reflectance distribution function (BRDF), spectral signature of the ground target, and atmospheric composition. The lifetime TOA reflectance trends from both sensors over 10 years are extremely stable, changing by no more than 0.4% per year in its TOA reflectance over the CEOS reference standard test sites.     

Production of Landsat ETM+ reference imagery of burned areas within Southern African savannahs: comparison of methods and application to MODIS

Accurate production of regional burned area maps are necessary to reduce uncertainty in emission estimates from African savannah fires. Numerous methods have been developed that map burned and unburned surfaces. These methods are typically applied to coarse spatial resolution (1 km) data to produce regional estimates of the area burned, while higher spatial resolution (<30 m) data are used to assess their accuracy with little regard to the accuracy of the higher spatial resolution reference data. In this study we aimed to investigate whether Landsat Enhanced Thematic Mapper (ETM+)‐derived reference imagery can be more accurately produced using such spectrally informed methods. The efficacy of several spectral index methods to discriminate between burned and unburned surfaces over a series of spatial scales (ground, IKONOS, Landsat ETM+ and data from the MOderate Resolution Imaging Spectrometer, MODIS) were evaluated. The optimal Landsat ETM+ reference image of burned area was achieved using a charcoal fraction map derived by linear spectral unmixing (k = 1.00, a = 99.5%), where pixels were defined as burnt if the charcoal fraction per pixel exceeded 50%. Comparison of coincident Landsat ETM+ and IKONOS burned area maps of a neighbouring region in Mongu (Zambia) indicated that the charcoal fraction map method overestimated the area burned by 1.6%. This method was, however, unstable, with the optimal fixed threshold occurring at >65% at the MODIS scale, presumably because of the decrease in signal‐to‐noise ratio as compared to the Landsat scale. At the MODIS scale the Mid‐Infrared Bispectral Index (MIRBI) using a fixed threshold of >1.75 was determined to be the optimal regional burned area mapping index (slope = 0.99, r ² = 0.95, SE = 61.40, y = Landsat burned area, x = MODIS burned area). Application of MIRBI to the entire MODIS temporal series measured the burned area as 10 267 km² during the 2001 fire season. The char fraction map and the MIRBI methodologies, which both produced reasonable burned area maps within southern African savannah environments, should also be evaluated in woodland and forested environments.