Zengyuan Li

Reprint of: Estimation of forest above-ground biomass using multi-parameter remote sensing data over a cold and arid area

Remote sensing is a valuable tool for estimating forest biomass in remote areas. This study explores retrieval of forest above-ground biomass (AGB) over a cold and arid region in Northwest China, using two different methods (non-parametric and parametric), field data, and three different remote sensing data: a SPOT-5 HRG image, multi-temporal dual-polarization ALOS PALSAR and airborne LiDAR data. The non-parametric method was applied in 300 different configurations, varying both the mathematical formulation and the data input (SPOT-5 and ALOS PALSAR), and the quality of the performance of each configuration was evaluated by Leave One Out (LOO) cross-validation against ground measurements. For the parametric method (the multivariate linear regression), the same remote sensing data were used, but in one additional configuration the airborne LiDAR data were used for stepwise multiple regression. The result of the best performing non-parametric configuration was satisfactory (R = 0.69 and RMSE = 20.7 tons/ha). The results for the parametric method were notoriously inaccurate, except for the case where airborne LiDAR data were included. The regression method with airborne low density LiDAR point cloud data was the best of all tested methods (R = 0.84 and RMSE = 15.2 tons/ha). A cross comparison of the two best results showed that the non-parametric method performs nearly as well as the parametric method with LiDAR data, except for some areas where forests have a very heterogeneous structure. It is concluded that the non-parametric method with SPOT data is able to map forest AGB operatively over the cold and arid region as an alternative to the more expensive airborne LiDAR data.

Improving measurement of forest structural parameters by co-registering of high resolution aerial imagery and low density LiDAR data.

Forest structural parameters, such as tree height and crown width, are indispensable for evaluating forest biomass or forest volume. LiDAR is a revolutionary technology for measurement of forest structural parameters, however, the accuracy of crown width extraction is not satisfactory when using a low density LiDAR, especially in high canopy cover forest. We used high resolution aerial imagery with a low density LiDAR system to overcome this shortcoming. A morphological filtering was used to generate a DEM (Digital Elevation Model) and a CHM (Canopy Height Model) from LiDAR data. The LiDAR camera image is matched to the aerial image with an automated keypoints search algorithm. As a result, a high registration accuracy of 0.5 pixels was obtained. A local maximum filter, watershed segmentation, and object-oriented image segmentation are used to obtain tree height and crown width. Results indicate that the camera data collected by the integrated LiDAR system plays an important role in registration with aerial imagery. The synthesis with aerial imagery increases the accuracy of forest structural parameter extraction when compared to only using the low density LiDAR data.

Forest stand biomass estimation using ALOS PALSAR data based on LiDAR-derived prior knowledge in the Qilian Mountain, western China

Studies are needed to evaluate the ability of Advanced Land Observing Satellite (ALOS) Phased Array type L-band Synthetic Aperture Radar (PALSAR) for forest aboveground biomass (AGB) extraction in mountainous areas. In this article, forest biomass was estimated at plot and stand levels, and different biomass grades, respectively. Light detection and ranging (LiDAR) data with about one hit per m2 were first used for forest biomass estimation at the plot level, with R 2 of 0.77. Then the LiDAR-derived biomass, as prior knowledge, was used to investigate the relationship between ALOS PALSAR data and biomass. The results showed that at each biomass level, the range of the back-scatter coefficient in HH and HV polarization (where H and V represent horizontal and vertical polarizations, respectively, and the first of the two letters refers to the transmission polarization and the second to the received polarization) was very large and there was no obvious relationship between the synthetic aperture radar (SAR) back-scatter coefficient and biomass at plot level. At stand level and in different biomass grades, the back-scatter coefficient increased with the increase of forest biomass, and a logarithm equation can be used to describe the relationship. The main reason may be that forest structure is complex at the plot level, while the average value could partly decrease the influence of forest structure at stand level. Meanwhile, terrain radiometric correction (TRC) was investigated and found effective for forest biomass estimation.

Estimating montane forest above-ground biomass in the upper reaches of the Heihe River Basin using Landsat-TM data

In this work, the results of above-ground biomass (AGB) estimates from Landsat Thematic Mapper 5 (TM) images and field data from the fragmented landscape of the upper reaches of the Heihe River Basin (HRB), located in the Qilian Mountains of Gansu province in northwest China, are presented. Estimates of AGB are relevant for sustainable forest management, monitoring global change, and carbon accounting. This is particularly true for the Qilian Mountains, which are a water resource protection zone. We combined forest inventory data from 133 plots with TM images and Advanced Spaceborne Thermal Emission and Reflection Radiometer (ASTER) global digital elevation model (GDEM) V2 products (GDEM) in order to analyse the influence of the sun-canopy-sensor plus C (SCS+C) topographic correction on estimations of forest AGB using the stepwise multiple linear regression (SMLR) and k-nearest neighbour (k-NN) methods. For both methods, our results indicated that the SCS+C correction was necessary for getting more reliable forest AGB estimates within this complex terrain. Remotely sensed AGB estimates were validated against forest inventory data using the leave-one-out (LOO) method. An optimized k-NN method was designed by varying both mathematical formulation of the algorithm and remote-sensing data input, which resulted in 3000 different model configurations. Following topographic correction, performance of the optimized k-NN method was compared to that of the regression method. The optimized k-NN method (R2 = 0.59, root mean square error (RMSE) = 24.92 tonnes ha–1) was found to perform much better than the regression method (R2 = 0.42, RMSE = 29.74 tonnes ha–1) for forest AGB retrieval over this montane area. Our results indicated that the optimized k-NN method is capable of operational application to forest AGB estimates in regions where few inventory data are available.

Wheat lodging monitoring using polarimetric index from RADARSAT-2 data

Abstract The feasibility of monitoring lodging of wheat fields by exploiting fully polarimetric C-band radar images has been investigated in this paper. A set of backscattering intensity features and polarimetric features, derived by target decomposition techniques, was extracted from 5 consecutive Radarsat-2 images. The temporal evolutions of these features of lodging wheat fields were investigated as a function of DAS (day after sowing) during the entire growing season. The temporal behavior was compared between typical lodging fields and normal fields in different growing stages. It was found that polarimetric feature from synthetic aperture radar (SAR) data was very sensitive to wheat lodging. Then a method called polarimetric index, availing the sensitivity of polarimetry to the structure, was put forward to monitor wheat lodging. The method was validated by two sets of in situ data collected in Shangkuli Farmland area, Inner Mongolia, China, at heading and ripe stages of spring wheat. Almost all the lodging fields were successfully distinguished from normal fields. Furthermore, the result revealed that the polarimetric index can reflect the intrinsic feature of lodging wheat with good anti-inference ability such as wheat growth difference. While optical sensors relied on its spectral features to monitor crop lodging, the proposed method based on radar data utilized polarimetric features to monitor crop lodging.

Filling invalid values in a lidar-derived canopy height model with morphological crown control

The light detection and ranging (lidar) technique has rapidly developed worldwide in numerous fields. The canopy height model (CHM), which can be generated from lidar data, is a useful model in forestry research. The CHM shows the canopy height above ground, and it indicates vertical elevation changes and the horizontal distribution of the canopy’s upper surface. Many vegetation parameters, which are important in forest inventory, can be extracted from the CHM. However, some abnormal or sudden changes of the height values (i.e. invalid values), which appear as unnatural holes in an image, exist in CHMs. This article proposes an approach to fill the invalid values in lidar-derived CHMs with morphological crown control. First, the Laplacian operator is applied to an original CHM to determine possible invalid values. Then, the morphological closing operator is applied to recover the crown coverage. By combining the two results, the possible invalid values in the CHM can be confirmed and replaced by corresponding values in the median-filtered CHM. The filling results from this new method are compared with those from other methods and with charge-coupled device images for evaluation. Finally, a CHM with random noise is used to test the filling correctness of the algorithm. The experiments show that this approach can fill the most invalid values well while refraining from overfilling.

Temporal Polarimetric Behavior of Oilseed Rape (Brassica napus L.) at C-Band for Early Season Sowing Date Monitoring

Spatial monitoring of the sowing date plays an important role in crop yield estimation at the regional scale. The feasibility of using polarimetric synthetic aperture radar (SAR) data for early season monitoring of the sowing dates of oilseed rape (Brassica napus L.) fields is explored in this paper. Polarimetric SAR responses of six parameters, relying on polarization decomposition methods, were investigated as a function of days after sowing (DAS) during the entire growing season, by means of five consecutive Radarsat-2 images. A near-continuous temporal evolution of these parameters was observed, based on 88 oilseed rape fields. It provided a solid basis for determining the suitable temporal window and the best polarimetric parameters for sowing date monitoring. A high sensitivity of all polarimetric parameters to the DAS at different growing stages was shown. Simple linear models could be calibrated to estimate sowing dates at early growth stages and were validated on an independent data set. Although Volume and Span parameters provided the highest sowing date estimation accuracy at the earlier growth stages, the other four parameters (Volume/Total, Odd/Total, Entropy and Alpha) were more accurate for a wider temporal window. These results demonstrate the capability and high potential of polarimetric SAR data for monitoring the sowing date of crops in the early season.

Quantitative Evaluation of Polarimetric Classification for Agricultural Crop Mapping

Agricultural crops classification capability of single band full polarization SAR data with different classification methods was evaluated using AIRSAR L-band polarimetric SAR data. It has been found that if only maximum likelihood (ML) classifiers, such as Wishart-maximum likelihood (WML) and normal distribution probability density functions (PDF)based Maximum Likelihood (NML) classifier can be utilized, it is better to choose WML directly applied to complex coherency or covariance matrix. NML cannot achieve acceptable classification result if intensity and phase images derived from coherency matrix are directly used for training the classifier. But if these images were supplied to the spatial-spectral based classifier, Extraction and Classification of Homogenous Objects (ECHO), higher classification accuracy can be obtained. Very low crop types discrimination accuracy has been observed when only H-Alpha polarimetric decomposition resultant images, such as entropy, alpha and anisotropy, were supplied to NML or a spatialspectral based classifier such as ECHO.

LiCHy: The CAF’s LiDAR, CCD and Hyperspectral Integrated Airborne Observation System

We describe the design, implementation and performance of a novel airborne system, which integrates commercial waveform LiDAR, CCD (Charge-Coupled Device) camera and hyperspectral sensors into a common platform system. CAF’s (The Chinese Academy of Forestry) LiCHy (LiDAR, CCD and Hyperspectral) Airborne Observation System is a unique system that permits simultaneous measurements of vegetation vertical structure, horizontal pattern, and foliar spectra from different view angles at very high spatial resolution (~1 m) on a wide range of airborne platforms. The horizontal geo-location accuracy of LiDAR and CCD is about 0.5 m, with LiDAR vertical resolution and accuracy 0.15 m and 0.3 m, respectively. The geo-location accuracy of hyperspectral image is within 2 pixels for nadir view observations and 5–7 pixels for large off-nadir observations of 55° with multi-angle modular when comparing to LiDAR product. The complementary nature of LiCHy’s sensors makes it an effective and comprehensive system for forest inventory, change detection, biodiversity monitoring, carbon accounting and ecosystem service evaluation. The LiCHy system has acquired more than 8000 km2 of data over typical forests across China. These data are being used to investigate potential LiDAR and optical remote sensing applications in forest management, forest carbon accounting, biodiversity evaluation, and to aid in the development of similar satellite configurations. This paper describes the integration of the LiCHy system, the instrument performance and data processing workflow. We also demonstrate LiCHy’s data characteristics, current coverage, and potential vegetation applications.

Isolating individual trees in a closed coniferous forest using small footprint lidar data

Numerous individual tree parameters can be isolated from a canopy height model derived from light detection and ranging. However, closed canopy forests make isolating individual trees difficult. A watershed method based on morphological crown control is studied in this research to isolate individual trees in a closed canopy forest and find how much the tree density influences the individual tree parameters. Morphological crown control is introduced to ensure that the watershed results locate in the crown area. The local maxima algorithm is used to identify potential tree positions in the crown area. Double watershed transformations, in which a simple reconstruction operation is inserted into the two transformations, are then applied to delineate the tree crowns. Finally, the individual trees are isolated, and their parameters are extracted. An experiment is conducted in a closed coniferous forest dominated by Picea crassifolia Kom. in northwest China. Results show that the proposed method can isolate the most dominant and subdominant trees, more than half of the intermediate trees and some suppressed trees.