基于遥感的区域尺度森林地上生物量估算研究

森林是陆地生态系统最大的碳库,精确估算森林生物量是陆地碳循环研究的关键。首先从机载LiDAR数据中提取高度和密度统计量,采用逐步回归模型进行典型样区生物量估算;然后利用机载LiDAR数据估算的生物量作为样本数据,与多光谱遥感数据Landsat8 OLI的波段反射率及植被指数建立回归模型,实现区域尺度森林地上生物量估算。实验结果显示,机载LiDAR数据估算的鼎湖山样区生物量与地面实测生物量的相关性R2达0.81,生物量RMSE为40.85 t/ha,说明机载LiDAR点云数据的高度和密度统计量与生物量存在较高的相关性。以机载LiDAR数据估算的生物量为样本数据,结合多光谱遥感数据Landsat8 OLI估算粤西北地区的森林地上生物量,精度验证结果为:R2为0.58,RMSE为36.9 t/ha;针叶林、阔叶林和针阔叶混交林等3种不同森林类型生物量的估算结果为:R2分别为0.51(n=251)、0.58(n=235)和0.56(n=241),生物量RMSE分别为24.1 t/ha、31.3 t/ha和29.9 t/ha,估算精度相差不大。总体上看,利用遥感数据可以开展区域尺度的森林地上生物量估算,为森林固碳监测提供有力的参考数据。     

Historical forest biomass dynamics modelled with Landsat spectral trajectories

Estimation of forest aboveground biomass (AGB) is informative of the role of forest ecosystems in local and global carbon budgets. There is a need to retrospectively estimate biomass in order to establish a historical baseline and enable reporting of change. In this research, we used temporal spectral trajectories to inform on forest successional development status in support of modelling and mapping of historic AGB for Mediterranean pines in central Spain. AGB generated with ground plot data from the Spanish National Forest Inventory (NFI), representing two collection periods (1990 and 2000), are linked with static and dynamic spectral data as captured by Landsat Thematic Mapper (TM) and Enhanced Thematic Mapper Plus (ETM+) sensors over a 25 year period (1984–2009). The importance of forest structural complexity on the relationship between AGB and spectral vegetation indices is revealed by the analysis of wavelet transforms. Two-dimensional (2D) wavelet transforms support the identification of spectral trajectory patterns of forest stands that in turn, are associated with traits of individual NFI plots, using a flexible algorithm sensitive to capturing time series similarity. Single-date spectral indices, temporal trajectories, and temporal derivatives associated with succession are used as input variables to non-parametric decision trees for modelling, estimation, and mapping of AGB and carbon sinks over the entire study area. Results indicate that patterns of change found in Normalized Difference Vegetation Index (NDVI) values are associated and relate well to classes of forest AGB. The Tasseled Cap Angle (TCA) index was found to be strongly related with forest density, although the related patterns of change had little relation with variability in historic AGB. By scaling biomass models through small (∼2.5 ha) spatial objects defined by spectral homogeneity, the AGB dynamics in the period 1990–2000 are mapped (70% accuracy when validated with plot values of change), revealing an increase of 18% in AGB irregularly distributed over 814 km² of pines. The accumulation of C calculated in AGB was on average 0.65 t ha⁻¹ y⁻¹, equivalent to a fixation of 2.38 t ha⁻¹ y⁻¹ of carbon dioxide.     

Predicting forest structural parameters using the image texture derived from WorldView-2 multispectral imagery in a dryland forest,Israel

Estimation of forest structural parameters by field-based data collection methods is both expensive and time consuming. Satellite remote sensing is a low-cost alternative in modeling and mapping structural parameters in large forest areas. The current study investigates the potential of using WordView-2 multispectral satellite imagery for predicting forest structural parameters in a dryland plantation forest in Israel. The relationships between image texture features and the several structural parameters such as Number of Trees (NT), Basal Area (BA), Stem Volume (SV), Clark-Evans Index (CEI), Diameter Differentiation Index (DDI), Contagion Index (CI), Gini Coefficient (GC), and Standard Deviation of Diameters at Breast Heights (SDDBH) were examined using correlation analyses. These variables were obtained from 30 m × 30 m square-shaped plots. The Standard Deviation of Gray Levels (SDGL) as a first order texture feature and the second order texture variables based on Gray Level Co-occurrence Matrix (GLCM) were calculated for the pixels that corresponds to field plots. The results of the correlation analysis indicate that the forest structural parameters are significantly correlated with the image texture features. The highest correlation coefficients were calculated for the relationships between the SDDBH and the contrast of red band (r = 0.75, p < 0.01), the BA and the entropy of blue band (r = 0.73, p < 0.01), and the GC and the contrast of blue band (r = 0.71, p < 0.01). Each forest structural parameter was modeled as a function of texture measures derived from the satellite image using stepwise multi linear regression analyses. The determination coefficient (R²) and root mean square error (RMSE) values of the best fitting models, respectively, are 0.38 and 109.56 ha⁻¹ for the NT; 0.54 and 1.79 m² ha⁻¹ for the BA; 0.42 and 27.18 m³ ha⁻¹ for the SV; 0.23 and 0.16 for the CEI; 0.32 and 0.05 for the DDI; 0.25 and 0.06 for the CI; 0.50 and 0.05 for the GC; and 0.67 and 0.70 for the SDDBH. The leave-one-out cross-validation technique was applied for validation of the best-fitted models (R² > 0.50). In conclusion, cross-validated statistics confirmed that the structural parameters including the BA, SDDBH, and GC can be predicted and mapped with a reasonable accuracy using the texture features extracted from the spectral bands of WorldView-2 image.     

Change detection and classification of land cover at Hustai National Park in Mongolia

Land cover types of Hustai National Park (HNP) in Mongolia, a hotspot area with rare species, were classified and their temporal changes were evaluated using Landsat MSS TM/ETM data between 1994 and 2000. Maximum-likelihood classification analysis showed an overall accuracy of 88.0% and 85.0% for the 1994 and 2000 images, respectively. Kappa coefficients associated with the classification were resulted to 0.85 for 1994 and 0.82 for 2000 image. Land cover types revealed significant temporal changes in the classification maps between 1994 and 2000. The area has increased considerably by 166.5 km² for mountain steppe and by 12 km² for a sand dune. By contrast, agricultural areas and degraded areas affected by human being activity were decreased by 46.1 km² and 194.8 km² over the 6-year span, respectively. These areas were replaced by mountain steppe area. Specifically, forest area was noticeably fragmented, accompanied by the decrease of ∼400 ha. The forest area revealed a pattern with systematic gain and loss associated with the specific phenomenon called as ‘forest free-south slope’. We discussed the potential environmental conditions responsible for the systematic pattern and addressed other biological impacts by outbreaks of forest pests and ungulates.     

Estimating the parameters of forest inventory using machine learning and the reduction of remote sensing features

Locally computed statistics of image texture and a case-based reasoning (CBR) system were evaluated for mapping of forest attributes. Cluster analysis was preferred to regression models, as a pre-selection method of features. The best stand-based accuracy using satellite sensor images was 74.64 m⁻³ ha⁻¹ (36%) RMSE for stand volume, 1.98 m⁻³ ha⁻¹ a⁻¹ (49%) for annual increase in stand volume, where κ = 0.23 for stand growth classes and κ = 0.41 for dominant tree species in stands. The top pixel-based accuracy using orthophotos was 76.54 m⁻³ ha⁻¹ (41%) RMSE for stand volume, 1.87 m⁻³ ha⁻¹ a⁻¹ (44%) for annual increase in stand volume, where κ = 0.24 for stand growth classes and κ = 0.38 for dominant tree species in stands. Mean saturation in 30 m radius was the most useful feature when orthophotos were used, and standard deviation of Landsat ETM 6.2 values in 80 m radius was the best when satellite sensor images were used. The most valuable feature components (radii, channels and local statistics) for orthophotos were: 30 m kernel radius, lightness and the mean of pixel values; for satellite sensor images: 80 m kernel radius, near-infrared channel (ETM 4) and the mean of pixel values. Locally computed statistics.     

Assessing Performance of Tomo-SAR and Backscattering Coefficient for Hemi-Boreal Forest Aboveground Biomass Estimation

Tomo-SAR technique has been used for hemi-boreal forest height and further forest biomass estimation through allometric equation. Backscattering coefficient especially in longer wavelength (L- or P-band) is thought as a useful parameter for hemi-boreal forest biomass retrieval. The aim of this paper is to assess the performance of vertical backscattering power and backscattering coefficient for hemi-boreal forest aboveground biomass (AGB) estimation with airborne P-band data. The test site locates in southern Sweden called Remningstorp test site, and the in-situ forest AGB ranges from 14 t/ha to 245 t/ha at stand level. Multi-baseline P-band Pol-InSAR data in repeat-path mode collected during March and May in 2007 at Remningstorp test site was used. We found that the correlation coefficient (R) between backscattering coefficient of P-band HH polarization and the in-situ forest biomass reached 0.87. The R for P-band VV backscattering power at 5 m is 0.71 and 10 m is 0.72. Backscattering coefficient in HH polarization and vertical backscattering power at 5 m and 10 m were applied to construct a model for hemi-boreal forest AGB estimation by backward step-wise regression and cross-validation approach. The results showed that the estimated forest AGB ranges from 19 to 240 t/ha, and the constructed model obtained a higher R and smaller RMSE, the value of R is 0.91, RMSE is 30.43 t/ha at Remningstorp test site.     

Accurate mapping of forest types using dense seasonal Landsat time-series

An accurate map of forest types is important for proper usage and management of forestry resources. Medium resolution satellite images (e.g., Landsat) have been widely used for forest type mapping because they are able to cover large areas more efficiently than the traditional forest inventory. However, the results of a detailed forest type classification based on these images are still not satisfactory. To improve forest mapping accuracy, this study proposed an operational method to get detailed forest types from dense Landsat time-series incorporating with or without topographic information provided by DEM. This method integrated a feature selection and a training-sample-adding procedure into a hierarchical classification framework. The proposed method has been tested in Vinton County of southeastern Ohio. The detailed forest types include pine forest, oak forest, and mixed-mesophytic forest. The proposed method was trained and validated using ground samples from field plots. The three forest types were classified with an overall accuracy of 90.52% using dense Landsat time-series, while topographic information can only slightly improve the accuracy to 92.63%. Moreover, the comparison between results of using Landsat time-series and a single image reveals that time-series data can largely improve the accuracy of forest type mapping, indicating the importance of phenological information contained in multi-seasonal images for discriminating different forest types. Thanks to zero cost of all input remotely sensed datasets and ease of implementation, this approach has the potential to be applied to map forest types at regional or global scales.     

Monitoring aboveground forest biomass dynamics over three decades using Landsat time-series and single-date inventory data

Understanding forest biomass dynamics is crucial for carbon and environmental monitoring, especially in the context of climate change. In this study, we propose a robust approach for monitoring aboveground forest biomass (AGB) dynamics by combining Landsat time-series with single-date inventory data. We developed a Random Forest (RF) based kNN model to produce annual maps of AGB from 1988 to 2017 over 7.2 million ha of forests in Victoria, Australia. The model was internally evaluated using a bootstrapping technique. Predictions of AGB and its change were then independently evaluated using multi-temporal Lidar data (2008 and 2016). To understand how natural and anthropogenic processes impact forest AGB, we analysed trends in relation to the history of disturbance and recovery. Specifically, change metrics (e.g., AGB loss and gain, Years to Recovery - Y2R) were calculated at the pixel level to characterise the patterns of AGB change resulting from forest dynamics. The imputation model achieved a RMSE value of 132.9 Mg ha⁻¹ (RMSE% = 46.3%) and R² value of 0.56. Independent assessments of prediction maps in 2008 and 2016 using Lidar-based AGB data achieved relatively high accuracies, with a RMSE of 108.6 Mg ha⁻¹ and 135.9 Mg ha⁻¹ for 2008 and 2016, respectively. Annual validations of AGB maps using un-changed, homogenous Lidar plots suggest that our model is transferable through time (RMSE ranging from 109.65 Mg ha⁻¹ to 112.27 Mg ha⁻¹ and RMSE% ranging from 25.38% to 25.99%). In addition, changes in AGB values associated with forest disturbance and recovery (decrease and increase, respectively) were captured by predicted maps. AGB change metrics indicate that AGB loss and Y2R varied across bioregions and were highly dependent on levels of disturbance severity (i.e., a greater loss and longer recovery time were associated with a higher severity disturbance). On average, high severity fire burnt from 200 Mg ha⁻¹ to 550 Mg ha⁻¹ of AGB and required up to 15 years to recover while clear-fell logging caused a reduction in 250 Mg ha⁻¹ to 600 Mg ha⁻¹ of AGB and required nearly 20 years to recover. In addition, AGB within un-disturbed forests showed statistically significant but monotonic trends, suggesting a mild gradual drop over time across most bioregions. Our methods are designed to support forest managers and researchers in developing forest monitoring systems, especially in developing regions, where only a single date forestry inventory exists.     

Terrestrial laser scanning to estimate plot-level forest canopy fuel properties

This paper evaluates the potential of a terrestrial laser scanner (TLS) to characterize forest canopy fuel characteristics at plot level. Several canopy properties, namely canopy height, canopy cover, canopy base height and fuel strata gap were estimated. Different approaches were tested to avoid the effect of canopy shadowing on canopy height estimation caused by deployment of the TLS below the canopy. Estimation of canopy height using a grid approach provided a coefficient of determination of R² = 0.81 and an RMSE of 2.47 m. A similar RMSE was obtained using the 99th percentile of the height distribution of the highest points, representing the 1% of the data, although the coefficient of determination was lower (R² = 0.70). Canopy cover (CC) was estimated as a function of the occupied cells of a grid superimposed upon the TLS point clouds. It was found that CC estimates were dependent on the cell size selected, with 3 cm being the optimum resolution for this study. The effect of the zenith view angle on CC estimates was also analyzed. A simple method was developed to estimate canopy base height from the vegetation vertical profiles derived from an occupied/non-occupied voxels approach. Canopy base height was estimated with an RMSE of 3.09 m and an R² = 0.86. Terrestrial laser scanning also provides a unique opportunity to estimate the fuel strata gap (FSG), which has not been previously derived from remotely sensed data. The FSG was also derived from the vegetation vertical profile with an RMSE of 1.53 m and an R² = 0.87.     

Stratified aboveground forest biomass estimation by remote sensing data

Remote sensing-assisted estimates of aboveground forest biomass are essential for modeling carbon budgets. It has been suggested that estimates can be improved by building species- or strata-specific biomass models. However, few studies have attempted a systematic analysis of the benefits of such stratification, especially in combination with other factors such as sensor type, statistical prediction method and sampling design of the reference inventory data. We addressed this topic by analyzing the impact of stratifying forest data into three classes (broadleaved, coniferous and mixed forest). We compare predictive accuracy (a) between the strata (b) to a case without stratification for a set of pre-selected predictors from airborne LiDAR and hyperspectral data obtained in a managed mixed forest site in southwestern Germany. We used 5 commonly applied algorithms for biomass predictions on bootstrapped subsamples of the data to obtain cross validated RMSE and r² diagnostics. Those values were analyzed in a factorial design by an analysis of variance (ANOVA) to rank the relative importance of each factor. Selected models were used for wall-to-wall mapping of biomass estimates and their associated uncertainty. The results revealed marginal advantages for the strata-specific prediction models over the unstratified ones, which were more obvious on the wall-to-wall mapped area-based predictions. Yet further tests are necessary to establish the generality of these results. Input data type and statistical prediction method are concluded to remain the two most crucial factors for the quality of remote sensing-assisted biomass models.     

Mapping beech (Fagus sylvatica L.) forest structure with airborne hyperspectral imagery

Estimating forest structural attributes using multispectral remote sensing is challenging because of the saturation of multispectral indices at high canopy cover. The objective of this study was to assess the utility of hyperspectral data in estimating and mapping forest structural parameters including mean diameter-at-breast height (DBH), mean tree height and tree density of a closed canopy beech forest (Fagus sylvatica L.). Airborne HyMap images and data on forest structural attributes were collected from the Majella National Park, Italy in July 2004. The predictive performances of vegetation indices (VI) derived from all possible two-band combinations (VI_(i,j) = (R_i − R_j)/(R_i + R_j), where R_i and R_j = reflectance in any two bands) were evaluated using calibration (n = 33) and test (n = 20) data sets. The potential of partial least squares (PLS) regression, a multivariate technique involving several bands was also assessed. New VIs based on the contrast between reflectance in the red-edge shoulder (756–820 nm) and the water absorption feature centred at 1200 nm (1172–1320 nm) were found to show higher correlations with the forest structural parameters than standard VIs derived from NIR and visible reflectance (i.e. the normalised difference vegetation index, NDVI). PLS regression showed a slight improvement in estimating the beech forest structural attributes (prediction errors of 27.6%, 32.6% and 46.4% for mean DBH, height and tree density, respectively) compared to VIs using linear regression models (prediction errors of 27.8%, 35.8% and 48.3% for mean DBH, height and tree density, respectively). Mean DBH was the best predicted variable among the stand parameters (calibration R² = 0.62 for an exponential model fit and standard error of prediction = 5.12 cm, i.e. 25% of the mean). The predicted map of mean DBH revealed high heterogeneity in the beech forest structure in the study area. The spatial variability of mean DBH occurs at less than 450 m. The DBH map could be useful to forest management in many ways, e.g. thinning of coppice to promote diameter growth, to assess the effects of management on forest structure or to detect changes in the forest structure caused by anthropogenic and natural factors.     

Mapping forest biomass from space – Fusion of hyperspectral EO1-hyperion data and Tandem-X and WorldView-2 canopy height models

Spaceborne sensors allow for wide-scale assessments of forest ecosystems. Combining the products of multiple sensors is hypothesized to improve the estimation of forest biomass. We applied interferometric (Tandem-X) and photogrammetric (WorldView-2) based predictors, e.g. canopy height models, in combination with hyperspectral predictors (EO1-Hyperion) by using 4 different machine learning algorithms for biomass estimation in temperate forest stands near Karlsruhe, Germany. An iterative model selection procedure was used to identify the optimal combination of predictors. The most accurate model (Random Forest) reached a r² of 0.73 with a RMSE of 14.9% (29.4 t/ha). Further results revealed that the predictive accuracy depended highly on the statistical model and the area size of the field samples. We conclude that a fusion of canopy height and spectral information allows for accurate estimations of forest biomass from space.     

Estimating aboveground and organ biomass of plant canopies across the entire season of rice growth with terrestrial laser scanning

Non-destructive and accurate estimation of crop biomass is crucial for the quantitative diagnosis of growth status and timely prediction of grain yield. As an active remote sensing technique, terrestrial laser scanning (TLS) has become increasingly available in crop monitoring for its advantages in recording structural properties. Some researchers have attempted to use TLS data in the estimation of crop aboveground biomass, but only for part of the growing season. Previous studies rarely investigated the estimation of biomass for individual organs, such as the panicles in rice canopies, which led to the poor understanding of TLS technology in monitoring biomass partitioning among organs. The objective of this study was to investigate the potential of TLS in estimating the biomass for individual organs and aboveground biomass of rice and to examine the feasibility of developing universal models for the entire growing season. The field plots experiments were conducted in 2017 and 2018 and involved different nitrogen (N) rates, planting techniques and rice varieties. Three regression approaches, stepwise multiple linear regression (SMLR), random forest regression (RF) and linear mixed-effects (LME) modeling, were evaluated in estimating biomass with extensive TLS and biomass data collected at multiple phenological stages of rice growth across the entire season. The models were calibrated with the 2017 dataset and validated independently with the 2018 dataset.The results demonstrated that growth stage in LME modeling was selected as the most significant random effect on rice growth among the three candidates, which were rice variety, growth stage and planting technique. The LME models grouped by growth stage exhibited higher validation accuracies for all biomass variables over the entire season to varying degrees than SMLR models and RF models. The most pronounced improvement with a LME model was obtained for panicle biomass, with an increase of 0.74 in R² (LME: R² = 0.90, SMLR: R² = 0.16) and a decrease of 1.15 t/ha in RMSE (LME: RMSE =0.79 t/ha, SMLR: RMSE =2.94 t/ha). Compared to SMLR and RF, LME modeling yielded similar estimation accuracies of aboveground biomass for pre-heading stages, but significantly higher accuracies for post-heading stages (LME: R² = 0.63, RMSE =2.27 t/ha; SMLR: R² = 0.42, RMSE =2.42 t/ha; RF: R² = 0.57, RMSE =2.80 t/ha). These findings implied that SMLR was only suitable for the estimation of biomass at pre-heading stages and LME modeling performed remarkably well across all growth stages, especially for post-heading. The results suggest coupling TLS with LME modeling is a promising approach to monitoring rice biomass at post-heading stages at high accuracy and to overcoming the saturation of canopy reflectance signals encountered in optical remote sensing. It also has great potential in the monitoring of other crops in cloud-cover conditions and the instantaneous prediction of grain yield any time before harvest.     

Estimation of aboveground biomass and carbon in a tropical rain forest in Gabon using remote sensing and GPS data

The knowledge of biomass stocks in tropical forests is critical for climate change and ecosystem services studies. This research was conducted in a tropical rain forest located near the city of Libreville (the capital of Gabon), in the Akanda Peninsula. The forest cover was stratified in terms of mature, secondary and mangrove forests using Landsat-ETM data. A field inventory was conducted to measure the required basic forest parameters and estimate the aboveground biomass (AGB) and carbon over the different forest classes. The Shuttle Radar Topography Mission (SRTM) data were used in combination with ground-based GPS measurements to derive forest heights. Finally, the relationships between the estimated heights and AGB were established and validated. Highest biomass stocks were found in the mature stands (223 ± 37 MgC/ha), followed by the secondary forests (116 ± 17 MgC/ha) and finally the mangrove forests (36 ± 19 MgC/ha). Strong relationships were found between AGB and forest heights (R² > 0.85).     

Time-series analysis of multi-resolution optical imagery for quantifying forest cover loss in Sumatra and Kalimantan,Indonesia

Monitoring loss of humid tropical forests via remotely sensed imagery is critical for a number of environmental monitoring objectives, including carbon accounting, biodiversity, and climate modeling science applications. Landsat imagery, provided free of charge by the U.S. Geological Survey Center for Earth Resources Observation and Science (USGS/EROS), enables consistent and timely forest cover loss updates from regional to biome scales. The Indonesian islands of Sumatra and Kalimantan are a center of significant forest cover change within the humid tropics with implications for carbon dynamics, biodiversity maintenance and local livelihoods. Sumatra and Kalimantan feature poor observational coverage compared to other centers of humid tropical forest change, such as Mato Grosso, Brazil, due to the lack of ongoing acquisitions from nearby ground stations and the persistence of cloud cover obscuring the land surface. At the same time, forest change in Indonesia is transient and does not always result in deforestation, as cleared forests are rapidly replaced by timber plantations and oil palm estates. Epochal composites, where single best observations are selected over a given time interval and used to quantify change, are one option for monitoring forest change in cloudy regions. However, the frequency of forest cover change in Indonesia confounds the ability of image composite pairs to quantify all change. Transient change occurring between composite periods is often missed and the length of time required for creating a cloud-free composite often obscures change occurring within the composite period itself. In this paper, we analyzed all Landsat 7 imagery with <50% cloud cover and data and products from the Moderate Resolution Imaging Spectroradiometer (MODIS) to quantify forest cover loss for Sumatra and Kalimantan from 2000 to 2005. We demonstrated that time-series approaches examining all good land observations are more accurate in mapping forest cover change in Indonesia than change maps based on image composites. Unlike other time-series analyses employing observations with a consistent periodicity, our study area was characterized by highly unequal observation counts and frequencies due to persistent cloud cover, scan line corrector off (SLC-off) gaps, and the absence of a complete archive. Our method accounts for this variation by generating a generic variable space. We evaluated our results against an independent probability sample-based estimate of gross forest cover loss and expert mapped gross forest cover loss at 64 sample sites. The mapped gross forest cover loss for Sumatra and Kalimantan was 2.86% of the land area, or 2.86 Mha from 2000 to 2005, with the highest concentration having occurred in Riau and Kalimantan Tengah provinces.     

Use of ETM+ images to extend stem volume estimates obtained from LiDAR data

Airborne LiDAR techniques can provide accurate measurements of tree height, from which estimates of stem volume and forest woody biomass can be obtained. These techniques, however, are still expensive to apply repeatedly over large areas. The current paper presents a methodology which first transforms mean stand heights obtained from LiDAR over small strips into relevant stem volume estimates. These are then extended over an entire forest by applying two estimation methods (k-NN and locally calibrated regression) to Landsat ETM+ images. The methodology is tested over a coastal area covered by pine forest in the Regional Park of San Rossore (Central Italy). The results are evaluated by comparison with the ground stem volumes of a recent forest inventory, taking into consideration the effect of stand size. In general, the accuracies of two estimation methods are dependent on the size of the forest stands and are satisfactory only when considering stands larger than 5-10 ha. The outputs of the parametric regression procedure are slightly more stable than those of k-NN and more faithfully reproduce the spatial patterns of the ground data.     

GLAS星载激光雷达和Landsat/ETM+数据的森林生物量估算

基于大脚印激光雷达数据和野外观测数据,该文提出一种获取脚印点内森林生物量的新思路,并结合陆地卫星数据应用于长白山地区森林地上生物量估算。首先,基于3种森林类型(针叶林、阔叶林和针阔混交林),采用多元逐步回归方法建立激光雷达波形指数与脚印点内实测平均树高的回归模型,估算全部脚印点内的平均树高;然后根据脚印点内样方的野外观测数据(平均树高和平均胸径)以及它们与样方生物量的拟合方程估算没有野外调查数据对应的脚印点的生物量;最后对3种森林类型的脚印点森林生物量在各森林覆盖度条件下进行分层分区统计得到生物量等级图。验证比较遥感估算的生物量与野外调查数据推算的生物量,总体误差在0~30(t·hm~(-2))之间,均方根误差为14.66(t·hm~(-2))。     

Evaluation of Landsat TM data for forest cover type and landuse classification in subtropical forests of Kumaon Himalaya (U.P.)

LANDSAT-TM has been evaluated for forest cover type and landuse classification in subtropical forests of Kumaon Himalaya (U.P.) Comparative evaluation of false colour composite generated by using various band combinations has been made. Digital image processing of Landsat-TM data on VIPS-32 RRSSC computer system has been carried out to stratify vegetation types. Conventional band combination in false colour composite is Bands 2, 3 and 4 in Red/Green/Blue sequence of Landsat TM for landuse classification. The present study however suggests that false colour combination using Landsat TM bands viz., 4, 5 and 3 in Red/Green/Blue sequence is the most suitable for visual interpretation of various forest cover types and landuse classes. It is felt that to extract full information from increased spatial and spectral resolution of Landsat TM, it is necessary to process the data digitally to classify land cover features like vegetation. Supervised classification using maximum likelihood algorithm has been attemped to stratify the forest vegetation. Only four bands are sufficient enough to classify vegetaton types. These bands are 2,3,4 and 5. The classification results were smoothed digitaly to increase the readiability of the map. Finally, the classification carred out using digital technique were evaluated using systematic sampling design. It is observed that forest cover type mapping can be achieved upto 80% overall mapping accuracy. Monospecies stand Chirpine can be mapped in two density classes viz., dense pine (<40%) with more than 90% accuracy. Poor accuracy (66%) was observed while mapping pine medium dense areas. The digital smoothening reduced the overall mapping accuracy. Conclusively, Landsat-TM can be used as operatonal sensor for forest cover type mapping even in complex landuse-terrain of Kumaon Himalaya (U.P.)     

Normalized algorithm for mapping and dating forest disturbances and regrowth for the United States

Forest disturbances such as harvesting, wildfire and insect infestation are critical ecosystem processes affecting the carbon cycle. Because carbon dynamics are related to time since disturbance, forest stand age that can be used as a surrogate for major clear-cut/fire disturbance information has recently been recognized as an important input to forest carbon cycle models for improving prediction accuracy. In this study, forest disturbances in the USA for the period of ∼1990–2000 were mapped using 400+ pairs of re-sampled Landsat TM/ETM scenes in 500m resolution, which were provided by the Landsat Ecosystem Disturbance Adaptive Processing System project. The detected disturbances were then separated into two five-year age groups, facilitated by Forest Inventory and Analysis (FIA) data, which was used to calculate the area of forest regeneration for each county in the USA.     

Effects of national forest inventory plot location error on forest carbon stock estimation using k-nearest neighbor algorithm

This paper suggested simulation approaches for quantifying and reducing the effects of National Forest Inventory (NFI) plot location error on aboveground forest biomass and carbon stock estimation using the k-Nearest Neighbor (kNN) algorithm. Additionally, the effects of plot location error in pre-GPS and GPS NFI plots were compared. Two South Korean cities, Sejong and Daejeon, were chosen to represent the study area, for which four Landsat TM images were collected together with two NFI datasets established in both the pre-GPS and GPS eras. The effects of plot location error were investigated in two ways: systematic error simulation, and random error simulation. Systematic error simulation was conducted to determine the effect of plot location error due to mis-registration. All of the NFI plots were successively moved against the satellite image in 360° directions, and the systematic error patterns were analyzed on the basis of the changes of the Root Mean Square Error (RMSE) of kNN estimation. In the random error simulation, the inherent random location errors in NFI plots were quantified by Monte Carlo simulation. After removal of both the estimated systematic and random location errors from the NFI plots, the RMSE% were reduced by 11.7% and 17.7% for the two pre-GPS-era datasets, and by 5.5% and 8.0% for the two GPS-era datasets. The experimental results showed that the pre-GPS NFI plots were more subject to plot location error than were the GPS NFI plots. This study’s findings demonstrate a potential remedy for reducing NFI plot location errors which may improve the accuracy of carbon stock estimation in a practical manner, particularly in the case of pre-GPS NFI data.