激光雷达回波模型辅助的坡地森林冠层高度反演

大光斑激光雷达数据已广泛应用于森林冠层高度提取,但通常仅限于地形坡度小于20°的平缓地区。在地形坡度大于20°的陡峭山区,地形引起的波形展宽使得地面回波和植被回波信息混合在一起,给森林冠层高度提取带来巨大挑战。本文利用激光雷达回波模型和地形信息,提出了一种模型辅助的坡地森林冠层高度反演算法。该方法以激光雷达回波信号截止点为参考,定义了波形高度指数H50和H75,使用激光雷达回波模型与已知地形信息模拟裸地的激光雷达回波,将裸地回波信号截止点与森林激光雷达回波信号截止点对齐,利用裸地回波计算常用的波形相对高度指数RH50和RH75,对森林冠层高度进行反演。并与高斯波形分解法和波形参数法的反演结果进行了比较。研究结果表明:(1)利用所提取的波形指数RH50和RH75对胸高断面积加权平均高(Lorey’s height)进行了估算,在坡度小于20°时,高斯波形分解法、波形参数法和模型辅助法的估算结果与实测值线性拟合的相关系数(R2)分别为0.70,0.78和0.98,对应的均方根误差(RMSE)分别为2.90 m,2.48 m和0.60 m,模型辅助法略优于其他两种方法;(2)在坡度大于20°时,高斯波形分解法、波形参数法和模型辅助法的R2分别为0.14,0.28和0.97,相应的RMSE分别为4.93 m,4.53 m和0.81 m,模型辅助法明显优于其他两种方法;(3)在0°—40°时,模型辅助法对Lorey’s height估算结果与实测值的R2为0.97,RMSE为0.80 m。本研究提出的模型辅助法具有更好的地形适应性,在0°—40°的坡度范围内具备对坡地森林冠层高度反演的潜力。     

Forest canopy height estimation using ICESat/GLAS data and error factor analysis in Hokkaido,Japan

Spaceborne light detection and ranging (LiDAR) enables us to obtain information about vertical forest structure directly, and it has often been used to measure forest canopy height or above-ground biomass. However, little attention has been given to comparisons of the accuracy of the different estimation methods of canopy height or to the evaluation of the error factors in canopy height estimation. In this study, we tested three methods of estimating canopy height using the Geoscience Laser Altimeter System (GLAS) onboard NASA’s Ice, Cloud, and land Elevation Satellite (ICESat), and evaluated several factors that affected accuracy. Our study areas were Tomakomai and Kushiro, two forested areas on Hokkaido in Japan. The accuracy of the canopy height estimates was verified by ground-based measurements. We also conducted a multivariate analysis using quantification theory type I (multiple-regression analysis of qualitative data) and identified the observation conditions that had a large influence on estimation accuracy. The method using the digital elevation model was the most accurate, with a root-mean-square error (RMSE) of 3.2 m. However, GLAS data with a low signal-to-noise ratio (⩽10.0) and that taken from September to October 2009 had to be excluded from the analysis because the estimation accuracy of canopy height was remarkably low. After these data were excluded, the multivariate analysis showed that surface slope had the greatest effect on estimation accuracy, and the accuracy dropped the most in steeply sloped areas. We developed a second model with two equations to estimate canopy height depending on the surface slope, which improved estimation accuracy (RMSE = 2.8 m). These results should prove useful and provide practical suggestions for estimating forest canopy height using spaceborne LiDAR.     

A novel approach to estimate canopy height using ICESat/GLAS data: A case study in the New Forest National Park,UK

The Geoscience Laser Altimeter System (GLAS) aboard Ice, Cloud and land Elevation Satellite (ICESat) is a spaceborne LiDAR sensor. It is the first LiDAR instrument which can digitize the backscattered waveform and offer near global coverage. Among others, scientific objectives of the mission include precise measurement of vegetation canopy heights. Existing approaches of waveform processing for canopy height estimation suggest Gaussian decomposition of the waveform which has the limitation to properly characterize significant peaks and results in discrepant information. Moreover, in most cases, Digital Terrain Models (DTMs) are required for canopy height estimation. This paper presents a new automated method of GLAS waveform processing for extracting vegetation canopy height in the absence of a DTM. Canopy heights retrieved from GLAS waveforms were validated with field measured heights. The newly proposed method was able to explain 79% of variation in canopy heights with an RMSE of 3.18 m, in the study area. The unexplained variation in canopy heights retrieved from GLAS data can be due to errors introduced by footprint eccentricity, decay of energy between emitted and received signals, uncertainty in the field measurements and limited number of sampled footprints.Results achieved with the newly proposed method were encouraging and demonstrated its potential of processing full-waveform LiDAR data for estimating forest canopy height. The study also had implications on future full-waveform spaceborne missions and their utility in vegetation studies.     

An improved dark object method to retrieve 500 m-resolution AOT (Aerosol Optical Thickness) image from MODIS data: A case study in the Pearl River Delta area,China

This paper presents an improved Dark Dense Vegetation (DDV) method for retrieving 500 m-resolution aerosol optical depth (AOT) based on MOD04-C005 arithmetic with the Moderate Resolution Imaging Spectroradiometer (MODIS) from the National Aeronautics and Space Administration (NASA). The improvements include change of the movement pattern of retrieval window, selection of a more suitable aerosol type, and storage of the look-up table. The method is then applied to obtain the AOT over the Pearl River Delta region (PRD). By comparing the results with the co-temporal ground sunphotometer observations in 2010, the correlation coefficient is found to be 0.794 with RMSE 0.139 and their variations remain consistent. Contrasts between model values in 2008 and MODIS AOT products in the same date also reveal a high accuracy of the improved DDV method. We also performed sensitivity tests to analyze the impacts of several parameters on apparent reflectance at different bands, and the results show that apparent reflectance is much more sensitive to surface reflectance and AOT than to elevation.     

A method for groundwater prospect zonation in data poor areas using remote sensing and GIS: a case study in Kalikavu Panchayath of Malappuram district,Kerala, India

Abstract

The present study was an attempt to delineate potential groundwater zones in Kalikavu Panchayat of Malappuram district, Kerala, India. The geo-spatial database on geomorphology, landuse, geology, slope and drainage network was generated in a geographic information system (GIS) environment from satellite data, Survey of India topographic sheets and field observations. To understand the movement and occurrence of groundwater, the geology, geomorphology, structural set-up and recharging conditions have to be well understood. In the present study, the potential recharge areas are delineated in terms of geology, geomorphology, land use, slope, drainage pattern, etc. Various thematic data generated were integrated using a heuristic method in the GIS domain to generate maps showing potential groundwater zones. The composite output map scores were reclassified into different zones using a decision rule. The final output map shows different zones of groundwater prospect, viz., very good (15.57% of the area), good (43.74%), moderate (28.38%) and poor (12.31%). Geomorphic units such as valley plains, valley fills and alluvial terraces were identified as good to excellent prospect zones, while the gently sloping lateritic uplands were identified as good to moderate zones. Steeply sloping hilly terrains underlain by hard rocks were identified as poor groundwater prospect zones.     

Development of an image based integrated method for determining and mapping aerosol optical thickness (AOT) over urban areas using the darkest pixel atmospheric correction method,RT equation and GIS: A case study of the Limassol area in Cyprus

This paper presents the development of an image-based integrated method for determining and mapping aerosol optical thickness (AOT). Using the radiative transfer (RT) equation, a methodology was developed to create a Geographical Information System (GIS) model that can visually display the AOT distribution over urban areas. In this paper, the model was applied to eleven Landsat Thematic Mapper/Enhanced Thematic Mapper Plus (TM/ETM+) satellite images of Limassol, Cyprus during 2010 and 2011 to determine the AOT levels in Limassol Cyprus during satellite overpass. The study is innovative and unique in that the RT equation, satellite images, the darkest pixel (DP) method of atmospheric correction and GIS were integrated to derive AOT from satellite images and display the AOT distribution over an urban area without the input of any meteorological or atmospheric parameters. The accuracy of the algorithm was verified through statistical analysis by the strong agreement between the AOT values derived using the algorithm and the in situ AOT values from the ground-based sensors.