基于地形调节植被指数估算长汀县植被覆盖度

植被覆盖度遥感估算最常用的方法是基于植被指数构建模型,但大部分的植被指数没有考虑地形的影响。以福建省长汀县作为研究区,引入能消除地形影响的地形调节植被指数(topography adjusted vegetation index,TAVI),利用像元二分模型估算植被覆盖度,旨在研究TAVI对植被覆盖度估算结果的影响,并与基于归一化差值植被指数(normalized difference vegetation index,NDVI)估算的结果进行比较。根据目视效果和统计指标的分析表明:基于TAVI估算的植被覆盖度精度高于基于NDVI的估算结果,并能有效降低阴坡阳坡间的差异,提高阴坡区域植被覆盖度的估算精度。     

Estimating fractional vegetation cover and the vegetation index of bare soil and highly dense vegetation with a physically based method

Normalized difference vegetation index (NDVI) of highly dense vegetation (NDVI_v) and bare soil (NDVI_s), identified as the key parameters for Fractional Vegetation Cover (FVC) estimation, are usually obtained with empirical statistical methods However, it is often difficult to obtain reasonable values of NDVI_v and NDVI_s at a coarse resolution (e.g., 1 km), or in arid, semiarid, and evergreen areas. The uncertainty of estimated NDVI_s and NDVI_v can cause substantial errors in FVC estimations when a simple linear mixture model is used. To address this problem, this paper proposes a physically based method. The leaf area index (LAI) and directional NDVI are introduced in a gap fraction model and a linear mixture model for FVC estimation to calculate NDVI_v and NDVI_s. The model incorporates the Moderate Resolution Imaging Spectroradiometer (MODIS) Bidirectional Reflectance Distribution Function (BRDF) model parameters product (MCD43B1) and LAI product, which are convenient to acquire. Two types of evaluation experiments are designed 1) with data simulated by a canopy radiative transfer model and 2) with satellite observations. The root-mean-square deviation (RMSD) for simulated data is less than 0.117, depending on the type of noise added on the data. In the real data experiment, the RMSD for cropland is 0.127, for grassland is 0.075, and for forest is 0.107. The experimental areas respectively lack fully vegetated and non-vegetated pixels at 1 km resolution. Consequently, a relatively large uncertainty is found while using the statistical methods and the RMSD ranges from 0.110 to 0.363 based on the real data. The proposed method is convenient to produce NDVI_v and NDVI_s maps for FVC estimation on regional and global scales.     

Seasonal trend analysis (STA) of MODIS vegetation index time series for the mangrove canopy of the Teacapan-Agua Brava lagoon system,Mexico

Monthly time series, from 2001 to 2016, of the Normalized Difference Vegetation Index (NDVI) and the Enhanced Vegetation Index (EVI) from MOD13Q1 products were analyzed with Seasonal Trend Analysis (STA), assessing seasonal and long-term changes in the mangrove canopy of the Teacapan-Agua Brava lagoon system, the largest mangrove ecosystem in the Mexican Pacific coast. Profiles from both vegetation indices described similar phenological trends, but the EVI was more sensitive in detecting intra-annual changes. We identified a seasonal cycle dominated by Laguncularia racemosa and Rhizophora mangle mixed patches, with the more closed canopy occurring in the early autumn, and the maximum opening in the dry season. Mangrove patches dominated by Avicennia germinans displayed seasonal peaks in the winter. Curves fitted for the seasonal vegetation indices were better correlated with accumulated precipitation and solar radiation among the assessed climate variables (Pearson’s correlation coefficients, estimated for most of the variables, were r ≥ 0.58 p < 0.0001), driving seasonality for tidal basins with mangroves dominated by L. racemosa and R. mangle. For tidal basins dominated by A. germinans, the maximum and minimum temperatures and monthly precipitation fit better seasonally with the vegetation indices (r ≥ 0.58, p < 0.0001). Significant mangrove canopy reductions were identified in all the analyzed tidal basins (z values for the Mann-Kendall test ≤ ?1.96), but positive change trends were recorded in four of the basins, while most of the mangrove canopy (approximately 87%) displayed only seasonal canopy changes or canopy recovery (z > ?1.96). The most resilient mangrove forests were distributed in tidal basins dominated by L. racemosa and R. mangle (Mann-Kendal Tau t ≥ 0.4, p ≤ 0.03), while basins dominated by A. germinans showed the most evidence of disturbance.     

Assessment of RapidEye vegetation indices for estimation of leaf area index and biomass in corn and soybean crops

Leaf area index (LAI) and biomass are important indicators of crop development and the availability of this information during the growing season can support farmer decision making processes. This study demonstrates the applicability of RapidEye multi-spectral data for estimation of LAI and biomass of two crop types (corn and soybean) with different canopy structure, leaf structure and photosynthetic pathways. The advantages of Rapid Eye in terms of increased temporal resolution (∼daily), high spatial resolution (∼5 m) and enhanced spectral information (includes red-edge band) are explored as an individual sensor and as part of a multi-sensor constellation. Seven vegetation indices based on combinations of reflectance in green, red, red-edge and near infrared bands were derived from RapidEye imagery between 2011 and 2013. LAI and biomass data were collected during the same period for calibration and validation of the relationships between vegetation indices and LAI and dry above-ground biomass. Most indices showed sensitivity to LAI from emergence to 8 m²/m². The normalized difference vegetation index (NDVI), the red-edge NDVI and the green NDVI were insensitive to crop type and had coefficients of variations (CV) ranging between 19 and 27%; and coefficients of determination ranging between 86 and 88%. The NDVI performed best for the estimation of dry leaf biomass (CV = 27% and r² = 090) and was also insensitive to crop type. The red-edge indices did not show any significant improvement in LAI and biomass estimation over traditional multispectral indices. Cumulative vegetation indices showed strong performance for estimation of total dry above-ground biomass, especially for corn (CV ≤ 20%). This study demonstrated that continuous crop LAI monitoring over time and space at the field level can be achieved using a combination of RapidEye, Landsat and SPOT data and sensor-dependant best-fit functions. This approach eliminates/reduces the need for reflectance resampling, VIs inter-calibration and spatial resampling.     

南方丘陵区植被覆盖度遥感估算的地形效应评估

植被覆盖变化是生态环境领域的核心研究内容之一,但其估算精度常受到地形效应、土壤背景、大气效应等各种因素影响。以Landsat 8 OLI为遥感数据源,基于像元二分模型,分别利用归一化差值植被指数(NDVI)、经Cosine-C校正的归一化差值植被指数(NDVI)和归一化差值山地植被指数(NDMVI)建立植被覆盖度估算模型,以评估南方丘陵区植被覆盖度的地形效应。结果表明,3种植被覆盖度估算模型均能削弱地形效应,但消除或抑制地形效应影响的能力不同。比较而言,基于NDMVI指数构建的植被覆盖度估算模型的地形效应最小,更适合地形复杂区域的植被覆盖度遥感估算;基于Cosine-C校正的NDVI植被指数构建的植被覆盖度估算模型的地形效应次之,但存在一定的过度校正现象;基于NDVI植被指数构建的植被覆盖度估算模型的地形效应最大,尤其当坡度≥10°时,阴坡植被覆盖度比阳坡明显偏低。     

Angle effects of vegetation indices and the influence on prediction of SPAD values in soybean and maize

To study the anisotropy of vegetation indices (VIs) and explore its influence on the retrieval accuracy of canopy soil-plant analyzer development (SPAD) value, the bidirectional reflectance distribution function (BRDF) models of soybean and maize are calculated from the multi-angle hyperspectral images acquired by UAV, respectively. According to the reflectance extracted from the BRDF model, the dependences of 16 commonly-used VIs on observation angles are analyzed, and the SPAD values of maize and soybean canopy are predicted by using the 16 VI values at different observation angles and their combinations as input parameters. The results show that the 16 VIs have different sensitivity to angle in the principal plane: green ratio vegetation index (GRVI), ratio vegetation index (RVI), red edge chlorophyll index (CIRE), and modified chlorophyll absorption in reflectance index/optimized soil-adjusted vegetation index (MCARI/OSAVI) are very sensitive to angles, among which MCARI/OSAVI of maize fluctuated the most (138.83 %); in contrast, the green optimal soil adjusted vegetation index (GOSAVI), normalized difference vegetation index (NDVI), and green normalized difference vegetation index (GNDVI) hardly change with the observation angles. In terms of SPAD prediction, the accuracy of different VI is different, the mean absolute error (MAE) showed that MCARI1 provided the highest accuracy of retrieval for soybean (MAE=1.617), while for maize it was MCARI/OSAVI (MAE=2.422). However, when using the same VI, there was no significant difference in the accuracy of the predicted results, whether the VI from different angles was used or the combination of multi-angles was used. The present results provide guiding significance and practical value for the retrieval of SPAD value in vegetation canopies and in-depth applications of multi-angular remote sensing.     

Comparative analysis of different uni- and multi-variate methods for estimation of vegetation water content using hyper-spectral measurements

Assessment of vegetation water content is critical for monitoring vegetation condition, detecting plant water stress, assessing the risk of forest fires and evaluating water status for irrigation. The main objective of this study was to investigate the performance of various mono- and multi-variate statistical methods for estimating vegetation water content (VWC) from hyper-spectral data. Hyper-spectral data is influenced by multi-collinearity because of a large number of (independent) spectral bands being modeled by a small number of (dependent) biophysical variables. Therefore, some full spectrum methods that are known to be suitable for analyzing multi-collinear data set were chosen. Canopy spectral reflectance was obtained with a GER 3700 spectro-radiometer (400–2400 nm) in a laboratory setting and VWC was measured by calculating wet/dry weight difference per unit of ground area (g/m²) of each plant canopy (n = 95). Three multivariate statistical methods were applied to estimate VWC: (1) partial least square regression, (2) artificial neural network and (3) principal component regression. They were selected to minimize the problem related to multi-collinearity. For comparison, uni-variate techniques including narrow band ratio water index (RWI), normalized difference water index (NDWI), second soil adjusted vegetation index (SAVI2) and transferred soil adjusted vegetation index (TSAVI) were applied. For each type of vegetation index, all two-band combinations were evaluated to determine the best band combination. Validation of the methods was based on the cross validation procedure and using three statistical indicators: R², RMSE and relative RMSE. The cross-validated results identified PLSR as the regression model providing the most accurate estimates of VWC among the various methods. The result revealed that this model is highly recommended for use with multi-collinear datasets ($R_{\text{CV}}^2=0.94$, RRMSE_CV = 0.23). Principal component regression exhibited the lowest accuracy among the multivariate models ($R_{\text{CV}}^2=0.78$, RRMSE_CV = 0.41).     

叶面积指数遥感尺度效应与尺度纠正

由于地表空间异质性的普遍存在,遥感反演模型的非线性必然会导致不同分辨率观测的遥感结果不一致,从而产生遥感产品尺度效应。本文研究了遥感产品尺度效应概念、模拟方法和定量计算模型,并利用锡林浩特草原研究区的实测数据,对尺度效应模型和方法进行了定量计算与验证分析。首先,基于不同升尺度方法与多尺度遥感成像机理之间的机理联系,通过“先反演再平均”与“先平均再反演”之间的差异,可计算“高”分辨率与“低”分辨率之间的遥感产品尺度差异。其次,分别以红光、近红外两波段反射率和归一化植被指数(NDVI)为自变量,对叶面积指数(LAI)非线性遥感模型进行泰勒展开,研究了模型非线性、遥感数据空间异质性对LAI遥感产品尺度差异的影响,发现高阶项可忽略,利用二阶导数项和遥感数据方差项可定量计算遥感产品尺度差异,经过二阶导数项纠正后的尺度差异相对偏差从5.6%分别降低到0.78%和1.45%。最后,分析了LAI遥感产品尺度效应的特征规律,得出以下结论:随着植被覆盖的增大,同等遥感空间异质性的LAI遥感产品尺度差异越大,且红光波段比近红外波段的尺度差异敏感性高近2个数量级;对于绝大部分陆地植被区域,存在“低分辨率低估”尺度效应,且遥感产品尺度差异的主导要素为LAI模型非线性,NDVI变量自身非线性对尺度效应贡献占23.5%;对于湿地类植被与水体混合情形,NDVI变量非线性的贡献为主导贡献,出现“低分辨率高估”尺度效应,必须利用红光、近红外两波段的二阶导数项非线性尺度差异,才能解释这一类型的LAI遥感产品尺度效应。本文建立了具有一定普适意义的遥感产品尺度效应定量模拟与尺度纠正方法,对推动定量遥感的尺度问题研究有一定参考价值。     

Dynamics of vegetation indices in tropical and subtropical savannas defined by ecoregions and Moderate Resolution Imaging Spectroradiometer (MODIS) land cover

In this study, we explored the capacity of vegetation indices derived from the Moderate Resolution Imaging Spectroradiometer (MODIS) reflectance products to characterize global savannas in Australia, Africa and South America. The savannas were spatially defined and subdivided using the World Wildlife Fund (WWF) global ecoregions and MODIS land cover classes. Average annual profiles of Normalized Difference Vegetation Index, shortwave infrared ratio (SWIR32), White Sky Albedo (WSA) and the Structural Scattering Index (SSI) were created. Metrics derived from average annual profiles of vegetation indices were used to classify savanna ecoregions. The response spaces between vegetation indices were used to examine the potential to derive structural and fractional cover measures. The ecoregions showed distinct temporal profiles and formed groups with similar structural properties, including higher levels of woody vegetation, similar forest–savanna mixtures and similar grassland predominance. The potential benefits from the use of combinations of indices to characterize savannas are discussed.     

Integrating temperature vegetation dryness index (TVDI) and regional water stress index (RWSI) for drought assessment with the aid of LANDSAT TM/ETM+ images

This paper presents a new drought assessment method by spatially and temporally integrating temperature vegetation dryness index (TVDI) with regional water stress index (RWSI) based on a synergistic approach. With the aid of LANDSAT TM/ETM data, we were able to retrieve the land-use and land-cover (LULC), vegetation indices (VIs), and land surface temperature (LST), leading to the derivation of three types of modified TVDI, including TVDI_SAVI, TVDI_ANDVI and TVDI_MSAVI, for drought assessment in a fast growing coastal area, Northern China. The categorical classification of four drought impact levels associated with the RWSI values enables us to refine the spatiotemporal relationship between the LST and the VIs. Holistic drought impact assessment between 1987 and 2000 was carried out by linking RWSI with TVDIs group wise. Research findings indicate that: (1) LST and VIs were negatively correlated in most cases of low, medium, and high vegetation cover except the case of high density vegetation cover in 2000 due to the effect of urban heat island (UHI) effect; (2) the shortage of water in 1987 was more salient than that that in 2000 based on all indices of TVDI and RWSI; and (3) TVDIs are more suitable for monitoring mild drought, normal and wet conditions when RWSI is smaller than 0.752; but they are not suitable for monitoring moderate and severe drought conditions.     

A comparison of geographic datasets and field measurements to model soil carbon using random forests and stepwise regressions (British Columbia,Canada)

We used geographic datasets and field measurements to examine the mechanisms that affect soil carbon (SC) storage for 65 grazed and non-grazed pastures in southern interior grasslands of British Columbia, Canada. Stepwise linear regression (SR) modeling was compared with random forest (RF) modeling. Models produced with SR performed better than those produced using RF models (r² = 0.56–0.77 AIC = 0.16–0.30 for SR models; r² = 0.38–0.53 and AIC = 0.18–0.30 for RF models). The factors most significant when predicting SC were elevation, precipitation, and the normalized difference vegetation index (NDVI). NDVI was evaluated at two scales using: (1) the MOD 13Q1 (250 m/16-day resolution) NDVI data product from the moderate resolution imaging spectro-radiometer (MODIS) (NDVI_MODIS), and (2) a handheld multispectral radiometer (MSR, 1 m resolution) (NDVI_MSR) in order to understand the potential for increasing model accuracy by increasing the spatial resolution of the gridded geographic datasets. When NDVI_MSR data were used to predict SC, the percentage of the variance explained by the model was greater than for models that relied on NDVI_MODIS data (r² = 0.68 for SC for non-grazed systems, modeled with SR based on NDVI_MODIS data; r² = 0.77 for SC for non-grazed systems, modeled with SR based on NDVI_MSR data). The outcomes of this study provide the groundwork for effective monitoring of SC using geographic datasets to enable a carbon offset program for the ranching industry.     

Retrieval of effective leaf area index (LAIe) and leaf area density (LAD) profile at individual tree level using high density multi-return airborne LiDAR

As an important canopy structure indicator, leaf area index (LAI) proved to be of considerable implications for forest ecosystem and ecological studies, and efficient techniques for accurate LAI acquisitions have long been highlighted. Airborne light detection and ranging (LiDAR), often termed as airborne laser scanning (ALS), once was extensively investigated for this task but showed limited performance due to its low sampling density. Now, ALS systems exhibit more competing capacities such as high density and multi-return sampling, and hence, people began to ask the questions like—“can ALS now work better on the task of LAI prediction?” As a re-examination, this study investigated the feasibility of LAI retrievals at the individual tree level based on high density and multi-return ALS, by directly considering the vertical distributions of laser points lying within each tree crown instead of by proposing feature variables such as quantiles involving laser point distribution modes at the plot level. The examination was operated in the case of four tree species (i.e. Picea abies, Pinus sylvestris, Populus tremula and Quercus robur) in a mixed forest, with their LAI-related reference data collected by using static terrestrial laser scanning (TLS). In light of the differences between ALS- and TLS-based LAI characterizations, the methods of voxelization of 3D scattered laser points, effective LAI (LAIe) that does not distinguish branches from canopies and unified cumulative LAI (ucLAI) that is often used to characterize the vertical profiles of crown leaf area densities (LADs) was used; then, the relationships between the ALS- and TLS-derived LAIes were determined, and so did ucLAIs. Tests indicated that the tree-level LAIes for the four tree species can be estimated based on the used airborne LiDAR (R² = 0.07, 0.26, 0.43 and 0.21, respectively) and their ucLAIs can also be derived. Overall, this study has validated the usage of the contemporary high density multi-return airborne LiDARs for LAIe and LAD profile retrievals at the individual tree level, and the contribution are of high potential for advancing forest ecosystem modeling and ecological understanding.     

Creating high-resolution bare-earth digital elevation models (DEMs) from stereo imagery in an area of densely vegetated deciduous forest using combinations of procedures designed for lidar point cloud filtering

For areas of the world that do not have access to lidar, fine-scale digital elevation models (DEMs) can be photogrammetrically created using globally available high-spatial resolution stereo satellite imagery. The resultant DEM is best termed a digital surface model (DSM) because it includes heights of surface features. In densely vegetated conditions, this inclusion can limit its usefulness in applications requiring a bare-earth DEM. This study explores the use of techniques designed for filtering lidar point clouds to mitigate the elevation artifacts caused by above ground features, within the context of a case study of Prince William Forest Park, Virginia, USA. The influences of land cover and leaf-on vs. leaf-off conditions are investigated, and the accuracy of the raw photogrammetric DSM extracted from leaf-on imagery was between that of a lidar bare-earth DEM and the Shuttle Radar Topography Mission DEM. Although the filtered leaf-on photogrammetric DEM retains some artifacts of the vegetation canopy and may not be useful for some applications, filtering procedures significantly improved the accuracy of the modeled terrain. The accuracy of the DSM extracted in leaf-off conditions was comparable in most areas to the lidar bare-earth DEM and filtering procedures resulted in accuracy comparable of that to the lidar DEM.     

基于RS和GIS的农作物估产方法研究进展

实时获取农作物长势及产量等信息对于现代农业的发展具有重要意义。近年来,随着遥感技术(remote sensing,RS)和地理信息系统(geographic information system,GIS)广泛应用于农作物估产领域,相继出现了一些较为实用的估产方法,主要有结合辅助数据的估产方法、基于植被指数的估产方法、基于特定模型的估产方法和基于农作物估产平台(软件)的开发等。其中,基于植被指数的估产方法又分为单一和多植被指数估产2类方法。在对近年来该领域大量文献深入研究的基础上,着重就几类热点方法展开论述,并对每类方法的优势和缺陷进行了评述,最后对该领域需要进一步研究的方向进行了探讨和展望,以期为后续研究提供参考。     

Remote sensing and information value (IV) model for regional mapping of fluvial channels and topographic wetness in the Saudi Arabia

This study aims to map regions of near surface fluvial channels, mega-basins and topographic wetness in Saudi Arabia using remote sensing data and an information value (IV) model, which is a modified approach of weight of evidence. We used the new version of the Shuttle Radar Topographic Mission (SRTM) to delineate the fluvial channels, mega-basin, and slope. These hydrological parameters were used to index the topographic wetness of each mega-basin in the region based on IV in a Geographic Information System. We validated our method using the Space Imaging Radar-C and Landsat 8 images and compared the textural features (fluvial channels) evident from SRTM digital elevation model and to determine whether these patterns were different. Our results revealed that the region is drained by nine tributaries and that the Err Rub Al Khali and Sahba mega-basins have the highest value of the IV and topographic wetness values; the Arran and coastal mega-basins have the lowest value of the IV and topographic wetness values. An integrated approach is timely and economically effective and can be applied throughout the arid and semi-arid regions to help hydrologists and urban developers.     

资源环境与城乡规划管理专业地理信息系统课程试验设计

根据地理信息系统课程的特点,针对资源环境与城乡规划管理专业,提出基于Arcview软件的地理信息系统课程实验设计,包括实验要求、实验内容与具体课时安排等。学生通过学习实践,可以加深理解地理信息系统的理论、技术和方法,为本专业相关领域应用问题的分析和解决打下基础。     

Bathymetrical influences on spatial and temporal characteristics of chlorophyll-a concentrations in the Southern Ocean from 2002 to 2012 (October to March) using MODIS

Phytoplankton blooms, particularly in the Southern Ocean, can have significant impact on global biogeochemistry cycling. To investigate the accuracy of chlorophyll-a distribution, and to better understand the spatial and temporal dynamics of phytoplankton biomass, we examine chlorophyll-a estimates (October–March from 2002 to 2012) derived from Moderate Resolution Imaging Spectrometer (MODIS) data following the ocean chlorophyll-a 3 model (OC3M) algorithm. Noticeable seasonality occurs in the temporal distribution of chlorophyll-a concentrations, which shows the highest value in December and January and an increasing tendency during the 2002–2012 period. The spatial distribution of chlorophyll-a varies greatly with latitude, as higher latitudes experience more phytoplankton blooms (chlorophyll-a concentration larger than 1 mg/m³) and marginal seas (Ross Sea and Amundsen Sea) show different bloom anomalies caused by two dominant algae species. Areas at higher latitudes and shallow water (<500 m) experience the shorter ice-free periods with greater seasonality. A noticeable bathymetry gradient exists at 2500-m isobaths, while water at the 500–2500-m depth experiences quite long ice-free periods with a stable water environment. Blooms generally occur near topographic features where currents have strong interactions when the water depth is more than 2500 m. Based on these findings, we can classify the Southern Ocean into two bloom subregions, 0–500 m as an enhanced bloom zone (EBZ), and 500–2500 m as a moderate bloom zone (MBZ). The EBZ has a quite high-bloom probability of about 30%, while the MBZ has only 10%.