Evaluating the potential of multispectral imagery to map multiple stages of tree mortality

Insect outbreaks are major forest disturbances, causing tree mortality across millions of ha in North America. Resultant spatial and temporal patterns of tree mortality can profoundly affect ecosystem structure and function. In this study, we evaluated the classification accuracy of multispectral imagery at different spatial resolutions. We used four-band digital aerial imagery (30-cm spatial resolution and aggregated to coarser resolutions) acquired over lodgepole pine-dominated stands in central Colorado recently attacked by mountain pine beetle. Classes of interest included green trees and multiple stages of post-insect attack tree mortality, including dead trees with red needles (“red-attack”), dead trees without needles (“gray-attack”), and non-forest. The 30-cm resolution image facilitated delineation of trees located in the field, which were used in image classification. We employed a maximum likelihood classifier using the green band, Red-Green Index (RGI), and Normalized Difference Vegetation Index (NDVI). Pixel-level classification accuracies using this imagery were good (overall accuracy of 87%, kappa = 0.84), although misclassification occurred between a) sunlit crowns of live (green) trees and herbaceous vegetation, and b) sunlit crowns of gray- and red-attack trees and bare soil. We explored the capability of coarser resolution imagery, aggregated from the 30-cm resolution to 1.2, 2.4, and 4.2 m, to improve classification accuracy. We found the highest accuracy at the 2.4-m resolution, where reduction in omission and commission errors and increases in overall accuracy (90%) and kappa (0.88) were achieved, and visual inspection indicated improved mapping. Pixels at this resolution included more shadow in forested regions than pixels in finer resolution imagery, thereby reducing forest canopy reflectance and allowing improved separation between forest and non-forest classes, yet were fine enough to resolve individual tree crowns better than the 4.2-m imagery. Our results illustrate that a classification of an image with a spatial resolution similar to the area of a tree crown outperforms that of finer and coarser resolution imagery for mapping tree mortality and non-forest classes. We also demonstrate that multispectral imagery can be used to separate multiple postoutbreak attack stages (i.e., red-attack and gray-attack) from other classes in the image.     

Assessing canopy mortality during a mountain pine beetle outbreak using GeoEye-1 high spatial resolution satellite data

Continuing, severe outbreaks of mountain pine beetle (Dendroctonus ponderosae) across western North America have resulted in widespread mortality of lodgepole pine (Pinus contorta). Multiple studies have used high spatial resolution satellite data to map areas of beetle kill; these studies have largely focused on mapping red canopy cover associated with recent tree mortality and have not examined mapping gray canopy cover that occurs after red needles have dropped. The work presented here examines the use of newly available GeoEye-1 data for mapping both red and gray canopy area in southeastern Wyoming lodgepole pine forest. A 0.5 m spatial resolution, pan-sharpened GeoEye-1 image was used to classify areas of green, red, and gray canopy cover. Reference data were collected at twelve 500 m² field plots. Shadow-normalized green, red, and gray canopy area from classified GeoEye-1 data closely agreed with field-estimated green, red, and gray canopy area. Mean absolute error in canopy cover for the twelve sample plots was 8.3% for the green class, 5.4% for the red class, and 7.2% for the gray class. When all twelve plots were aggregated, remotely sensed estimates of green, red, and gray cover were within 1.7% of the field-estimated cover. Our results demonstrate that high spatial resolution spaceborne multispectral data are a promising tool for mapping canopy mortality caused by mountain pine beetle outbreaks.     

Sensitivity of the thematic mapper enhanced wetness difference index to detect mountain pine beetle red-attack damage

Red-attack damage caused by mountain pine beetle (Dentroctonus ponderosa Hopkins) infestation in stands of lodgepole pine (Pinus contorta) in the Prince George Forest Region of British Columbia was examined using multitemporal Landsat-7 ETM+ imagery acquired in 1999, 2000, and 2001. The image data were geometrically and atmospherically corrected, and processed using the Tasseled Cap Transformation (TCT) to obtain wetness indices. The final steps included pixel subtraction, enhancement, and thresholding of the wetness index differences. The resulting enhanced wetness difference index (EWDI) was used to interpret spectral patterns in stands with confirmed (through aerial survey) red-attack damage in 2001, and these EWDI patterns were compared to the patterns of reflectance in normal-colour composites. We stratified the aerial survey dataset into two levels and used the EWDI to discriminate classes of 10-29 red-attack trees and 30-50 red-attack trees, and a sample of healthy forest collected from inventory data. Classification accuracy of red-attack damage based on the EWDI ranged from 67% to 78% correct.     

Detection of mountain pine beetle-killed ponderosa pine in a heterogeneous landscape using high-resolution aerial imagery

Previous studies have used remote-sensing images to map tree mortality caused by mountain pine beetle (Dendroctonus ponderosae; MPB) in relatively homogeneous lodgepole pine (Pinus contorta) forests; however, classification methods have not been tested for the patchy landscape of ponderosa pine-dominated (Pinus ponderosae) montane forests characterized by highly variable tree density. This study explores two supervised classification methods to identify MPB-caused mortality (red attack) in heterogeneous montane forests of the Colorado Front Range using 1 m-resolution 2011 imagery of the National Agriculture Imagery Program (NAIP): maximum likelihood using the red, green, and blue bands, and the red-green index (RGI), and a thresholding technique using the RGI. Two variations of the RGI threshold method were also explored: the addition of a green-band threshold and the incorporation of a focal analysis. Evaluation pixels were used to assess the accuracy of the classification methods. The maximum likelihood (97 Percentage Correctly Classified (PCC); 11% error of commission for red attack) and RGI threshold (85 PCC; 46% error of commission for red attack) classification methods overestimated the red attack. The RGI and green band threshold classification reduced the error of commission (5%) and had high overall accuracy (97 PCC). In a comparison of classification methods across tree-density sites, we found the maximum likelihood classification had a very high accuracy in the high-density site (95 PCC), but substantially lower accuracy in the low-density site (85 PCC) due to the presence of more visible cover types. The RGI threshold classification with the green band constraint produced more consistent PCCs across tree densities: high (93.7 PCC), moderate (95.2 PCC), and low (92.0 PCC). Our results indicate forest structure may affect the classification accuracy and should be considered when selecting a classification method for a landscape.     

Assessment of QuickBird high spatial resolution imagery to detect red attack damage due to mountain pine beetle infestation

High spatial resolution remotely sensed data has the potential to complement existing forest health programs for both strategic planning over large areas, as well as for detailed and precise identification of tree crowns subject to stress and infestation. The area impacted by the current mountain pine beetle (Dendroctonus ponderosae Hopkins) outbreak in British Columbia, Canada, has increased 40-fold over the previous 5 years, with approximately 8.5 million ha of forest infested in 2005. As a result of the spatial extent and intensity of the outbreak, new technologies are being assessed to help detect, map, and monitor the damage caused by the beetle, and to inform mitigation of future beetle outbreaks. In this paper, we evaluate the capacity of high spatial resolution QuickBird multi-spectral imagery to detect mountain pine beetle red attack damage. ANOVA testing of individual spectral bands, as well as the Normalized Difference Vegetation Index (NDVI) and a ratio of red to green reflectance (Red-Green Index or RGI), indicated that the RGI was the most successful (p < 0.001) at separating non-attack crowns from red attack crowns. Based on this result, the RGI was subsequently used to develop a binary classification of red attack and non-attack pixels. The total number of QuickBird pixels classified as having red attack damage within a 50 m buffer of a known forest health survey point were compared to the number of red attack trees recorded at the time of the forest health survey. The relationship between the number of red attack pixels and observed red attack crowns was assessed using independent validation data and was found to be significant (r² = 0.48, p < 0.001, standard error = 2.8 crowns). A comparison of the number of QuickBird pixels classified as red attack, and a broader scale index of mountain pine beetle red attack damage (Enhanced Wetness Difference Index, calculated from a time series of Landsat imagery), was significant (r² = 0.61, p < 0.001, standard error = 1.3 crowns). These results suggest that high spatial resolution imagery, in particular QuickBird satellite imagery, has a valuable role to play in identifying tree crowns with red attack damage. This information could subsequently be used to augment existing detailed forest health surveys, calibrate synoptic estimates of red attack damage generated from overview surveys and/or coarse scale remotely sensed data, and facilitate the generation of value-added information products, such as estimates of timber volume impacts at the forest stand level.     

Parametrized BRDF for atmospheric and topographic correction and albedo estimation in Jiangxi rugged terrain,China

A method is presented for bi‐directional reflectance distribution function (BRDF) parametrization for topographic correction and surface reflectance estimation from Landsat Thematic Mapper (TM) over rugged terrain. Following this reflectance, albedo is calculated accurately. BRDF is parametrized using a land‐cover map and Landsat TM to build a BRDF factor to remove the variation of relative solar incident angle and relative sensor viewing angle per pixel. Based on the BRDF factor and radiative transfer model, solar direct radiance correction, sky diffuse radiance and adjacent terrain reflected radiance correction were introduced into the atmospheric‐topographic correction method. Solar direct radiance, sky diffuse radiance and adjacent terrain reflected radiance, as well as atmospheric transmittance and path radiance, are analysed in detail and calculated per pixel using a look‐up table (LUT) with a digital elevation model (DEM). The method is applied to Landsat TM imagery that covers a rugged area in Jiangxi province, China. Results show that atmospheric and topographic correction based on BRDF gives better surface reflectance compared with sole atmospheric correction and two other useful atmospheric‐topographic correction methods. Finally, surface albedo is calculated based on this topography‐corrected reflectance and shows a reasonable accuracy in albedo estimation.     

Mapping percentage tree cover from Envisat MERIS data using linear and nonlinear techniques

The aim of this study was to predict percentage tree cover from Envisat Medium Resolution Imaging Spectrometer (MERIS) imagery with a spatial resolution of 300 m by comparing four common models: a multiple linear regression (MLR) model, a linear mixture model (LMM), an artificial neural network (ANN) model and a regression tree (RT) model. The training data set was derived from a fine spatial resolution land cover classification of IKONOS imagery. Specifically, this classification was aggregated to predict percentage tree cover at the MERIS spatial resolution. The predictor variables included the MERIS wavebands plus biophysical variables (the normalized difference vegetation index (NDVI), leaf area index (LAI), fraction of photosynthetically active radiation (fPAR), fraction of green vegetation covering a unit area of horizontal soil (fCover) and MERIS terrestrial chlorophyll index (MTCI)) estimated from the MERIS data. An RT algorithm was the most accurate model to predict percentage tree cover based on the Envisat MERIS bands and vegetation biophysical variables. This study showed that Envisat MERIS data can be used to predict percentage tree cover with considerable spatial detail. Inclusion of the biophysical variables led to greater accuracy in predicting percentage tree cover. This finer-scale depiction should be useful for environmental monitoring purposes at the regional scale.     

The potential for integrating Sentinel 2 MSI with SPOT 5 HRG and Landsat 8 OLI imagery for monitoring semi-arid savannah woody cover

Multitemporal archived imagery enables the monitoring of savannah woody cover, for ecological purposes. Compatibility in multitemporal, multiple sensor image data would facilitate the monitoring. The decommissioning of SPOT 5 (Système Pour l’Observation de la Terre 5) left a void in multispectral imagery at the 10 m spatial resolution of its high-resolution geometric (HRG) sensor. The subsequent launch of Sentinel 2 presented an opportunity for data continuity to monitor the savannah woody cover, using equivalent 10 m resolution multispectral instrument (MSI) bands. This study examined the integration potential of Sentinel 2 MSI with the longer archive HRG and Landsat 8 (Land Satellite 8) Operational Land Imager (OLI) imagery, in assessing savannah woody cover. Images of three semi-arid savannah sites acquired on same season dates that excluded herbaceous vegetation from the spectral signature were used: November 2014 (HRG) and December 2015 (MSI, OLI). Using equivalent green (G), red (R), and near infrared (NIR) bands at 10 m (MSI, HRG) and 30 m (OLI) resolution, the woody cover was mapped through subpixel classification. The mapped woody cover was compared for statistical differences using χ² analysis at 10 m resolution (MSI, HRG) and at a degradation of the MSI and HRG images to the 30 m OLI pixel size. Conversion to top-of-atmosphere reflectance values facilitated inter-sensor correlation of G, R, and NIR reflectance for field sampling sites where woody cover was quantified. Inter-sensor regression functions in G, R, and NIR band MSI and HRG images were developed. The 10 m resolution classifications of woody cover were not statistically different. Due to spatial resolution similarity, SPOT 5 HRG multispectral imagery was established as suitable for integration with equivalent band MSI imagery in mapping the woody cover in a multitemporal analysis. For dense woody cover, Landsat 8 OLI imagery was more suitable for integration with MSI than HRG images due to higher radiometric sensitivity, which can permit monitoring physiology-related woody reflectance.     

Integrating remotely sensed and ancillary data sources to characterize a mountain pine beetle infestation

The ongoing mountain pine beetle (Dendroctonus ponderosae Hopkins) outbreak in British Columbia, Canada, has reached epidemic proportions, with the beetle expanding into geographic areas outside its known biological range. In this study, estimates of red attack damage were derived from a logistic regression model using multi-date Landsat imagery, and ancillary information including terrain attributes and solar radiation. The model estimates were found to be approximately 70% accurate using an independent set of beetle survey data as validation. This probability surface of red attack damage, along with forest inventory and terrain attributes, were used as inputs to decision tree analyses, in order to identify which forest attributes were associated with stands that had a greater likelihood of mountain pine beetle red attack damage. Three distinct decision tree models were developed, with each having a different set of input variables. The results of the analyses indicated that site index (an indicator of the quality of a forest site) and slope were the principal discriminators of the current mountain pine beetle attack, followed by basal area of pine dominated stands, and to a lesser extent, crown closure and stem density. The results suggest that indicators of site quality, particularly site index, could be a complementary addition to existing stand susceptibility rating models.     

Using hyperspectral remote sensing to detect and quantify southeastern pine senescence effects in red-cockaded woodpecker (Picoides borealis) habitat

Conservation of threatened and endangered species requires maintenance of critical habitat. The red-cockaded woodpecker Picoides borealis (RCW) is a threatened bird species, endemic to the mixed conifer forests of the southeastern United States. RCW nests and forages preferentially in mature longleaf pine Pinus palustris, but also uses mature loblolly pine Pinus taeda and shortleaf pine Pinus echinata forests. In the last century, the extent of mature pine forests has been greatly reduced by logging. The RCW, in contrast to other woodpeckers, excavates nest cavities in living trees and senescence symptoms (year round leaf chlorosis and leaf mortality) have been observed in mature pine stands across the southeast. Widespread mortality of the mature pine forests would threaten the long-term survival of the RCW. We used airborne hyperspectral data across a portion of Ft. Benning Military Installation, Georgia, U.S.A., to determine if senescent trees can be identified and mapped and assess the likely persistence of mature pines in the RCW habitat. Univariate analysis of variance showed good separation between asymptomatic, senesced and dead physiological conditions with asymptomatic trees having significantly higher reflectance for all bands in the wavelength range between 0.719 and 1.1676 µm, senescent trees having significantly lower reflectance for bands in the range between 1.1927 and 1.3122 µm, and dead trees having significantly higher reflectance for bands in the range between 1.8151 and 1.9471 µm. Classification and Regression Tree (CART) models achieved correct classification rates and kappa statistics above 70%. CART models selected information from wavelength regions similar to those identified from the ANOVA, which likely explains their performance. Our aggregated CART model of tree senescence estimated that 141.4 ha (3%) of the study area is affected. RCW nests occurred in areas with significantly higher tree cover, and trees within foraging and home ranges did not have significantly more senescence than areas without RCW. These results indicate that although tree senescence is widespread, mortality is yet to significantly affect RCW habitat. Results and analysis of critical habitat similar to those exemplified in this study can extend our knowledge about the stressors of these important and imperiled components of biodiversity.     

A non-parametric,supervised classification of vegetation types on the Kaibab National Forest using decision trees

Traditional land classification techniques for large areas that use Landsat Thematic Mapper (TM) imagery are typically limited to the fixed spatial resolution of the sensors (30 m). However, the study of some ecological processes requires land cover classifications at finer spatial resolutions. We model forest vegetation types on the Kaibab National Forest (KNF) in northern Arizona to a 10-m spatial resolution with field data, using topographical information and Landsat TM imagery as auxiliary variables. Vegetation types were identified by clustering the field variables total basal area and proportion of basal area by species, and then using a decision tree based on auxiliary variables to predict vegetation types. Vegetation types modelled included pinyon-juniper, ponderosa pine, mixed conifer, spruce- and deciduous-dominated mixes, and openings. To independently assess the accuracy of the final vegetation maps using reference data from different sources, we used a post-stratified, multivariate composite estimator. Overall accuracy was 74.5% (Kappa statistic = 49.9%). Sources of error included differentiating between mixed conifer and spruce-dominated types and between openings in the forest and deciduous-dominated mixes. Overall, our non-parametric classification method successfully identified dominant vegetation types on the study area at a finer spatial resolution than can typically be achieved using traditional classification techniques.     

Estimation and vicarious validation of urban vegetation abundance by spectral mixture analysis

Both moderate and high spatial resolution imagery can be used to quantify abundance and distribution of urban vegetation for urban landscape management and to provide inputs to physical process models. Estimation of vegetation fraction from Landsat ETM+ and Quickbird allows for operational monitoring and reconnaissance at moderate resolution with calibration and vicarious validation at higher resolution. Establishing a linear correspondence between ETM-derived vegetation fraction and Quickbird-derived vegetation fraction facilitates the validation task by extending the spatial scale from 30 × 30 m to a more manageable 2.8 × 2.8 m. A comparative analysis indicates that urban reflectance can be accurately represented with a three component linear mixture model for both Landsat ETM+ and Quickbird imagery in the New York metro area. The strong linearity of the Substrate Vegetation Dark surface (SVD) mixture model provides consistent estimates of illuminated vegetation fraction that can be used to constrain physical process models that require biophysical inputs related to vegetation abundance. When Quickbird-derived 2.8 m estimates of vegetation fraction are integrated to 30 m scales and coregistered to Landsat-derived 30 m estimates, median estimates agree with the integrated fractions to within 5% for fractions > 0.2. The resulting Quickbird-ETM+ scatter distribution cannot be explained with estimate error alone but is consistent with a 3% to 6% estimation error combined with a 17 m subpixel registration ambiguity. The 3D endmember fraction space obtained from ETM+ imagery forms a ternary distribution of reflectance properties corresponding to distinct biophysical surface types. The SVD model is a reflectance analog to Ridd's V–I–S land cover model but acknowledges the fact that permeable and impermeable surfaces cannot generally be distinguished on the basis of broadband reflectance alone. We therefore propose that vegetation fraction be used as a proxy for permeable surface distribution to avoid the common erroneous assumption that all nonvegetated surfaces along the gray axis are completely impermeable. Comparison of mean vegetation fractions to street tree counts in New York City shows a consistent relationship between minimum fraction and tree count. However, moderate and high resolution areal estimates of vegetation fraction provide complementary information because they image all illuminated vegetation, including that not counted by the in situ street tree inventory.     

Continuous wavelet analysis for the detection of green attack damage due to mountain pine beetle infestation

Mountain pine beetle (Dendroctonus ponderosae Hopkins) is the most destructive insect infesting mature pine forests in North America and has devastated millions of hectares of forest in western Canada. Past studies have demonstrated the use of multispectral imagery for remote identification and mapping of visible or red attack damage in forests. This study aims to detect pre-visual or green attack damage in lodgepole pine needles by means of hyperspectral measurements, particularly via continuous wavelet analysis. Field measurements of lodgepole pine stands were conducted at two sites located northwest of Edmonton, Alberta, Canada. In June and August of 2007, reflectance spectra (350-2500 nm) were collected for 16 pairs of trees. Each of the 16 tree pairs included one control tree (healthy), and one stressed tree (girdled to simulate the effects of beetle infestation). In addition, during the period of June through October 2008, spectra were collected from 15 pairs of control- and beetle-infested trees. Spectra derived from these 31 tree pairs were subjected to a continuous wavelet transform, generating a scalogram that compiles the wavelet power as a function of wavelength location and scale of decomposition. Linear relationships were then explored between the wavelet scalograms and chemical properties or class labels (control and non-control) of the sample populations in order to isolate the most useful distinguishing spectral features that related to infested or girdled trees vs. control trees.A deficit in water content is observed in infested trees while an additional deficit in chlorophyll content is seen for girdled trees. The measurable water deficit of infested and girdled tree samples was detectable from the wavelet analysis of the reflectance spectra providing a novel method for the detection of green attack. The spectral features distinguishing control and infested trees are predominantly located between 950 and 1390 nm from scales 1 to 8. Of those, five features between 1318 to 1322 nm at scale 7 are consistently found in the July and August 2008 datasets. These features are located at longer wavelengths than those investigated in previous studies (below 1100 nm) and provide new insights into the potential remote detection of green attack. Spectral features that distinguish control and girdled trees were mostly observed between 1550 and 2370 nm from scales 1 to 5. The differing response of girdled and infested trees appears to indicate that the girdling process does not provide a perfect simulation of the effects caused by beetle infestation.It remains to be determined if the location of the 1318-1322 nm features, near the edge of a strong atmospheric water absorption band, will be sufficiently separable for use in airborne detection of green attack. A plot comparing needle water content and wavelet power at 1320 nm reveals considerable overlap between data derived from both infested and control samples, though the groups are statistically separable. This obstacle may preclude a high accuracy separation of healthy and infected single individuals, but establishing threshold identification levels may provide an economical, efficient and expeditious method for discriminating between healthy and infested tree populations.     

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.     

Using satellite remote sensing for DEM extraction in complex mountainous terrain: Landscape analysis of the Makalu Barun National Park of eastern Nepal

The design and management of national parks and other protected areas requires a broad base of physiographic and geo-ecological information about the landscape. This paper evaluates the effectiveness of satellite remote sensing for photogrammetric stereo-mapping and digital elevation model (DEM) extraction within remote mountainous terrain. As a case study, a landscape analysis of the Makalu Barun National Park and Conservation Area of east Nepal (27.5° N, 87.0° E) was examined. The study area is a highly complex and rugged mountain landscape, with extreme topographic relief and an elevation gradient spanning more than 8300 m. A DEM extracted from stereo SPOT imagery resulted in a median disagreement of 58 m when compared to a DEM generated from a conventionally digitized GIS dataset of topographic contours (scale=1:250 000). Visual comparison of the two DEMs showed substantial agreement at the landscape scale, while larger scale comparison of 100 m contours revealed some localized differences. The SPOT extracted DEM provided equal or better basis for orthorectification of satellite imagery when compared to the conventional DEM. Derivative landscape analysis outputs, such as hydrological modelling, drainage networks and watershed boundaries, compared well with results based upon the conventional dataset. Intermediate map products useful for field research and mapping included production of an orthorectified satellite base-map image. Additionally, a fused multisensor high resolution image of the study area, combining Landsat Thematic Mapper (TM) and SPOT imagery at 10 m resolution, was orthorectified to produce a false-colour satellite image map highlighting the spectral discrimination between land cover classes.     

High-spectral resolution data for monitoring Scots pine (Pinus sylvestris L.) regeneration

High-resolution field spectroradiometric measurements were used to investigate the spectral properties of naturally regenerating Scots pine ( Pinus sylvestris L.) in relation to sapling cover and season. Ratios of reflectance (R757/ R722) and first-derivative reflectance (D719/ D703 and D730/ D700) in the red edge region correlated most strongly with increasing sapling cover. The main control on these ratios was the contrast between the biomass of heather and pine saplings. At the site investigated, the strength of the correlation of reflectance ratios with sapling cover was not affected by seasonal changes in canopy colour.     

Sub-pixel estimation of urban land cover components with linear mixture model analysis and Landsat Thematic Mapper imagery

We examine the utility of linear mixture modelling in the sub-pixel analysis of Landsat Enhanced Thematic Mapper (ETM) imagery to estimate the three key land cover components in an urban/suburban setting: impervious surface, managed/unmanaged lawn and tree cover. The relative effectiveness of two different endmember sets was also compared. The interior endmember set consisted of the median pixel value of the training pixels of each land cover and the exterior endmember set was the extreme pixel value. As a means of accuracy assessment, the resulting land cover estimates were compared with independent estimates obtained from the visual interpretation of digital orthophotography and classified IKONOS imagery. Impervious surface estimates from the Landsat ETM showed a high degree of similarity (RMS error (RMSE) within approximately ±10 to 15%) to that obtained using high spatial resolution digital orthophotography and IKONOS imagery. The partition of the vegetation component into tree vs grass cover was more problematic due to the greater spectral similarity between these land cover types with RMSE of approximately ±12 to 22%. The interior endmember set appeared to provide better differentiation between grass and urban tree cover than the exterior endmember set. The ability to separate the grass vs tree components in urban vegetation is of major importance to the study of the urban/suburban ecosystems as well as watershed assessment.     

A comparison of multi-spectral, multi-angular, and multi-temporal remote sensing datasets for fractional shrub canopy mapping in Arctic Alaska

Shrub cover appears to be increasing across many areas of the Arctic tundra biome, and increasing shrub cover in the Arctic has the potential to significantly impact global carbon budgets and the global climate system. For most of the Arctic, however, there is no existing baseline inventory of shrub canopy cover, as existing maps of Arctic vegetation provide little information about the density of shrub cover at a moderate spatial resolution across the region. Remotely-sensed fractional shrub canopy maps can provide this necessary baseline inventory of shrub cover. In this study, we compare the accuracy of fractional shrub canopy (> 0.5 m tall) maps derived from multi-spectral, multi-angular, and multi-temporal datasets from Landsat imagery at 30 m spatial resolution, Moderate Resolution Imaging SpectroRadiometer (MODIS) imagery at 250 m and 500 m spatial resolution, and MultiAngle Imaging Spectroradiometer (MISR) imagery at 275 m spatial resolution for a 1067 km² study area in Arctic Alaska. The study area is centered at 69 °N, ranges in elevation from 130 to 770 m, is composed primarily of rolling topography with gentle slopes less than 10°, and is free of glaciers and perennial snow cover. Shrubs > 0.5 m in height cover 2.9% of the study area and are primarily confined to patches associated with specific landscape features. Reference fractional shrub canopy is determined from in situ shrub canopy measurements and a high spatial resolution IKONOS image swath. Regression tree models are constructed to estimate fractional canopy cover at 250 m using different combinations of input data from Landsat, MODIS, and MISR. Results indicate that multi-spectral data provide substantially more accurate estimates of fractional shrub canopy cover than multi-angular or multi-temporal data. Higher spatial resolution datasets also provide more accurate estimates of fractional shrub canopy cover (aggregated to moderate spatial resolutions) than lower spatial resolution datasets, an expected result for a study area where most shrub cover is concentrated in narrow patches associated with rivers, drainages, and slopes. Including the middle infrared bands available from Landsat and MODIS in the regression tree models (in addition to the four standard visible and near-infrared spectral bands) typically results in a slight boost in accuracy. Including the multi-angular red band data available from MISR in the regression tree models, however, typically boosts accuracy more substantially, resulting in moderate resolution fractional shrub canopy estimates approaching the accuracy of estimates derived from the much higher spatial resolution Landsat sensor. Given the poor availability of snow and cloud-free Landsat scenes in many areas of the Arctic and the promising results demonstrated here by the MISR sensor, MISR may be the best choice for large area fractional shrub canopy mapping in the Alaskan Arctic for the period 2000-2009.     

Geo-atmospheric processing of airborne imaging spectrometry data. Part 2: Atmospheric/topographic correction

A method for the radiometric correction of wide field-of-view airborne imagery has been developed that accounts for the angular dependence of the path radiance and atmospheric transmittance functions to remove atmospheric and topographic effects. The first part of processing is the parametric geocoding of the scene to obtain a geocoded, orthorectified image and the view geometry (scan and azimuth angles) for each pixel as described in part 1 of this jointly submitted paper. The second part of the processing performs the combined atmospheric/ topographic correction. It uses a database of look-up tables of the atmospheric correction functions (path radiance, atmospheric transmittance, direct and diffuse solar flux) calculated with a radiative transfer code. Additionally, the terrain shape obtained from a digital elevation model is taken into account. The issues of the database size and accuracy requirements are critically discussed. The method supports all common types of imaging airborne optical instruments: panchromatic, multispectral and hyperspectral, including fore/aft tilt sensors covering the wavelength range 0.35-2.55 w m and 8-14 w m. The processor is designed and optimized for imaging spectrometer data. Examples of processing of hyperspectral imagery in flat and rugged terrain are presented. A comparison of ground reflectance measurements with surface reflectance spectra derived from airborne imagery demonstrates that an accuracy of 1-3% reflectance units can be achieved.     

The Comparison of Urban Tree Crown Extraction based on Airborne LiDARElevation Difference and High Resolution Imagery

With the support of airborne Light Detection and Ranging (LiDAR) data and high spatial resolution aerial imagery,this paper presents an individual tree extraction method that takes the region of urban as the study area.The elevation difference model originated from LiDAR data was used to extract regions of interest including trees. Then,masking was applied to the high spatial resolution aerial imagery to get the same regions. Besides,image segmentations,based on the marked watershed algorithm,were processed on the high spatial resolution aerial imagery and the elevation difference model separately to extract individual tree crowns. Finally,we took a visual interpretation to delineate tree crowns manually and this result was regarded as the reference crowns map. The extraction accuracies were assessed by comparing the spatial relationships of the reference crowns and the automated delineated tree crowns based on the elevation difference model and the high resolution imagery. The results show that the LiDAR data is developed to improve the efficiency of obtaining forest region that the canopy height model include 85.25% forest information. In addition,the tree crowns extraction accuracy based on the high resolution aerial imagery is 57.14%,while another extraction accuracy based on the elevation difference model is 42.47%,which indicated that the marked watershed algorithm proposed in this paper is effective and the high resolution imagery is better than the elevation difference model to extract tree crowns.