冠层反射光谱对植被理化参数的全局敏感性分析

植被理化参数与许多有关植物物质能量交换的生态过程密切相关,定量分析植被反射光谱对理化参数的敏感性是遥感反演理化参数含量的前提。本文采用EFAST(Extended Fourier Amplitude Sensitivity Test)全局敏感性分析方法,利用PROSAIL辐射传输模型分析了冠层疏密程度对叶片生化组分含量、冠层结构以及土壤背景等多种参数敏感性的影响,并对植被理化参数反演所需先验知识的精度问题进行了初步探讨。研究表明:(1)对于较为稠密的冠层,可见光波段的冠层反射率主要受叶绿素含量的影响,近红外和中红外波段的冠层反射率主要受干物质量和含水量的影响;(2)对于稀疏的冠层,LAI是影响400—2500 nm波段范围内冠层反射率的最重要参数,土壤湿度次之,叶片生化参数对冠层反射率的敏感性较低;(3)在已知稀疏冠层LAI的情况下进一步确定土壤的干湿状态,可显著提高冠层反射率对叶绿素含量的敏感度,有助于稀疏冠层叶绿素含量的反演。     

针对SAIL冠层模型的土壤背景反射率修正

提出一种基于SAIL模型的地表反射率修正方案,有效减小地形起伏的影响。通过引入太阳直射光的方向-方向反射与大气散射的半球-方向反射,遵循光路可逆原理对地表反射率进行几何修正,同时考虑地表自身热辐射对入瞳辐射的影响从而修正地表反射率,发展适用于SAIL模型的地表反射率修正模型。利用长常高速部分路段的实测植被理化参数及光谱信息对地形修正后的SAIL模型模拟精度进行对比分析,结果表明地形修正后SAIL模型有效提高SAIL模型模拟的植被冠层光谱精度,修正后SAIL模型可为后续南方地区定量遥感的应用提供更精确的数据支持。     

Estimation of forest LAI by inverting canopy reflectance models and multi-angle imagery

Computer simulation models have seldom been applied for estimating the structural and biophysical variables of forest canopy. In this study, an approach for the estimation of leaf area index (LAI) using the information contained in hyperspectral, multi-angle images and the inversion of a computer simulation model are explored. For this purpose, L-systems combined with forest growth model ZELIG were applied to render 3-D forest architectural scenarios. The Radiosity-graphics combined model (RGM) was used to estimate forest LAI from the Compact High-Resolution Imaging Spectrometer/Project for On-Board Autonomy (CHRIS/PROBA) data. LAI inversion was performed using the look-up table (LUT) method. The estimated LAI was evaluated against in situ LAI measurement and compared against the LAI predictions from CHRIS data obtained using the Li-Strahler geometric-optical canopy reflectance model (GOMS). The results indicated that the method used in this study can be efficient strategy to estimate LAI by RGM model inversion.     

Retrieval of wheat leaf area index using price approach based on inversion of canopy reflectance model

This paper reports estimation of the Leaf Area Index (LAI) of wheat crop from IRS-LISS-III data using Price (1993) approach. Empirical approach for LAI estimation with different NDVI estimation procedures viz. radiance, apparent reflectance and dark object subtraction (DOS) based atmospheric correction were also evaluated. Validation of LAI retrieval and NDVI normalizations were carried out using field level measurements of crop LAI and spectral property using canopy analyzer and spectro-radiometer, respectively over selected fields in Bhopal District, Madhya Pradesh. It was observed that empirical relations are sensitive to the NDVI estimation approach and DOS method performed better as compared to other two approaches. It was also observed that LAI estimation from Price algorithm is sensitive to the crop attenuation coefficients. Crop specific attenuation coefficients reported in literature for Indian cultivars gave higher accuracy. The root mean square (RMS) error of 0.77 for LAI estimation was achieved using above described approach.     

Vegetation chlorophyll estimates in the Amazon from multi-angle MODIS observations and canopy reflectance model

As a preparatory study for future hyperspectral missions that can measure canopy chemistry, we introduce a novel approach to investigate whether multi-angle Moderate Resolution Imaging Spectroradiometer (MODIS) data can be used to generate a preliminary database with long-term estimates of chlorophyll. MODIS monthly chlorophyll estimates between 2000 and 2015, derived from a fully coupled canopy reflectance model (ProSAIL), were inspected for consistency with eddy covariance fluxes, tower-based hyperspectral images and chlorophyll measurements. MODIS chlorophyll estimates from the inverse model showed strong seasonal variations across two flux-tower sites in central and eastern Amazon. Marked increases in chlorophyll concentrations were observed during the early dry season. Remotely sensed chlorophyll concentrations were correlated to field measurements (r² = 0.73 and r² = 0.98) but the data deviated from the 1:1 line with root mean square errors (RMSE) ranging from 0.355 μg cm⁻² (Tapajós tower) to 0.470 μg cm⁻² (Manaus tower). The chlorophyll estimates were consistent with flux tower measurements of photosynthetically active radiation (PAR) and net ecosystem productivity (NEP). We also applied ProSAIL to mono-angle hyperspectral observations from a camera installed on a tower to scale modeled chlorophyll pigments to MODIS observations (r² = 0.73). Chlorophyll pigment concentrations (Chl_A+B) were correlated to changes in the amount of young and mature leaf area per month (0.59 ≤  r² ≤  0.64). Increases in MODIS observed Chl_A+B were preceded by increased PAR during the dry season (0.61 ≤ r² ≤  0.62) and followed by changes in net carbon uptake. We conclude that, at these two sites, changes in LAI, coupled with changes in leaf chlorophyll, are comparable with seasonality of plant productivity. Our results allowed the preliminary development of a 15-year time series of chlorophyll estimates over the Amazon to support canopy chemistry studies using future hyperspectral sensors.     

Effects of solar and view angles on reflectance for paddy field canopies

It is very important to know the spectral characteristics for the sake of understanding the remote sensing data. The reflectance characteristics of paddy field canopies vary with time or observational conditions (solar zenith angle, solar azimuth angle, and view zenith angle). A number of field studies have clarified the effects of these conditions on grain canopy reflectance. Most of the field data used in these study, however, were conducted only through the growing season in one year or by grains planted in pots. A series of authors’ experiments were initiated in 1982 and continued from the spring to the autumn every year to 1987. In this study we describe that the remotely sensed spectral data measured on the ground are influenced not only by the grain type, observational conditions, and growing season but also by the solar zenith angle, solar azimuth angle and view zenith angle in relation to scene. In this paper we report the results from the investigation of these various fundamental properties.     

基于计算机模拟模型LESS的落叶松林木质元素影响分析

目前大部分植被辐射传输模型在模拟太阳辐射与植被之间的相互作用过程时,将植物结构进行了简化,只保留了叶片的结构和空间分布特征,而忽略了木质元素(枝干等)对冠层反射特性的影响.计算机模拟模型LESS能够充分考虑植被的多种组分光谱和结构特征,精确模拟植被冠层内部的光散射和辐射过程.本文以地面实测数据为基础,发展了以单木为基本...     

Assessing the impact of hydrocarbon leakages on vegetation using reflectance spectroscopy

This paper assesses the capability of hyperspectral remote sensing to detect hydrocarbon leakages in pipelines using vegetation status as an indicator of contamination. A field experiment in real scale and in tropical weather was conducted in which Brachiaria brizantha H.S. pasture plants were grown over soils contaminated with small volumes of liquid hydrocarbons (HCs). The contaminations involved volumes of hydrocarbons that ranged between 2 L and 12.7 L of gasoline and diesel per m³ of soil, which were applied to the crop parcels over the course of 30 days. The leaf and canopy reflectance spectra of contaminated and control plants were acquired within 350–2500 nm wavelengths. The leaf and canopy reflectance spectra were mathematically transformed by means of first derivative (FD) and continuum removal (CR) techniques. Using principal component analysis (PCA), the spectral measurements could be grouped into either two or three contamination groups. Wavelengths in the red edge were found to contain the largest spectral differences between plants at distinct, evolving contamination stages. Wavelengths centred on water absorption bands were also important to differentiating contaminated from healthy plants. The red edge position of contaminated plants, calculated on the basis of FD spectra, shifted substantially to shorter wavelengths with increasing contamination, whereas non-contaminated plants displayed a red shift (in leaf spectra) or small blue shift (in canopy spectra). At leaf scale, contaminated plants were differentiated from healthy plants between 550–750 nm, 1380–1550 nm, 1850–2000 nm and 2006–2196 nm. At canopy scale, differences were substantial between 470–518 nm, 550–750 nm, 910–1081 nm, 1116–1284 nm, 1736–1786 nm, 2006–2196 nm and 2222–2378 nm. The results of this study suggests that remote sensing of B. brizantha H.S. at both leaf and canopy scales can be used as an indicator of gasoline and diesel contaminations for the detection of small leakages in pipelines.     

Modeling three-dimensional forest structures to drive canopy radiative transfer simulations of bidirectional reflectance factor

Three-dimensional (3-D) Monte Carlo-based radiative transfer (MCRT) models are usually used for benchmarking in intercomparisons of the canopy radiative transfer (RT) simulations. However, the 3-D MCRT models are rarely applied to develop remote sensing algorithms to estimate essential climate variables of forests, due mainly to the difficulties in obtaining realistic stand structures for different forest biomes over regional to global scales. Fortunately, some of important tree structure parameters such as canopy height and tree density distribution have been available globally. This enables to run the intermediate complexities of the 3-D MCRT models. We consequently developed a statistical approach to generate forest structures with intermediate complexities depending on the inputs of canopy height and tree density. It aims at facilitating applications of the 3-D MCRT models to develop remote sensing retrieval algorithms. The proposed approach was evaluated using field measurements of two boreal forest stands at Estonia and USA, respectively. Results demonstrated that the simulations of bidirectional reflectance factor (BRF) based on the measured forest structures agreed well with the BRF based on the generated structures from the proposed approach with the root mean square error (RMSE) and relative RMSE (rRMSE) ranging from 0.002 to 0.006 and from 0.7% to 19.8%, respectively. Comparison of the computed BRF with corresponding MODIS reflectance data yielded RMSE and rRMSE lower than 0.03 and 20%, respectively. Although the results from the current study are limited in two boreal forest stands, our approach has the potential to generate stand structures for different forest biomes.     

Retrieval of leaf area index in alpine wetlands using a two-layer canopy reflectance model

In this paper, we focused on the retrieval of the LAI in an alpine wetland located in western part of China in late August and early July 2011. A two-layer canopy reflectance model (ACRM) was used to establish the relationships between the LAI and the reflectance of near-infrared (NIR) and red (RED) wavebands. The reflectance data were derived from Landsat TM L1T product and the Terra and Aqua MODIS 16-day and 8-day composite reflectance products (MOD/MYD09) at 250 m resolution. Due to the lack of the information about some major input parameters for ACRM, which are sensitive to model outputs in the reflectance of NIR and RED wavebands, the inverse problem was ill-posed. To overcome this problem, a method of increasing the sensitivity of the LAI while reducing the influence of other model free parameters based on the study of free parameters’ sensitivity to the ACRM outputs and the region’s features was studied. The area of interest was divided into two parts using the approximately statistic normalized difference vegetation index (NDVI) value around 0.5. One part was sparse vegetation (0.1 < NDVI < 0.5), which is more sensitive to soil background effects and less sensitive to the canopy biophysical and biochemical variables. The other part was dense vegetation (0.5 ≤ NDVI < 1.0), which is less sensitive to soil background effects and more sensitive to plant canopies and leaf parameters. Then, the relationships of ρ_nir–LAI and ρ_red–LAI were established using a look-up table algorithm for the two parts. Furthermore, a regularization technique for fast pixel-wise retrieval was introduced to reduce the elements of LUT sets while maintaining a relatively high accuracy. The results were very promising compared to the field measured LAI values that the correlation (R²) of the measured LAI values and retrieved LAI values reached 0.95, and the root-mean-square deviation (RMSD) was 0.33 for late August, 2011, while the R² reached 0.82 and RMSD was 0.25 for early July 2011.     

Estimation of forest canopy structural parameters using kernel-driven bi-directional reflectance model based multi-angular vegetation indices

Near-surface bi-directional reflectance and high-spatial resolution true-color imagery of several forested canopies were acquired using an unmanned helicopter. The observed reflectance from multiple view-zenith angles were simulated with a kernel-driven bidirectional reflectance model, and the BRDF parameters were retrieved. Based on the retrieved BRDF parameters, kernel-derived multi-angular vegetation indices (KMVIs) were computed. The potential of KMVI for prediction of canopy structural parameters such as canopy fraction and canopy volume was assessed. The performance of each KMVI was tested by comparison to field measured canopy fraction and canopy volume. For the prediction of canopy fraction, the KMVI that included the nadir-based NDVI performed better than other KMVI emphasizing the importance of nadir observation for remote estimation of the canopy fraction. The Nadir BRDF-adjusted NDVI was found to be superior for the prediction of canopy fraction, which could explain 77% variation of the canopy fraction. However, none of the existing KMVI predicted the canopy volume better than Nadir BRDF-adjusted NDVI and Nadir-view NDVI. The Canopy structural index (CSI) was proposed with the combination of normalized difference between dark-spot near infrared reflectance and hot-spot red reflectance. The CSI could establish an improved relationship with the canopy volume over Nadir BRDF-adjusted NDVI and Nadir-view NDVI, explaining 72% variation in canopy volume. In addition, MODIS based KMVI were evaluated for the prediction of canopy fraction and canopy volume. MODIS based KMVI also showed similar results to the helicopter based KMVI. The promising results shown by the CSI suggest that it could be an appropriate candidate for remote estimation of three-dimensional canopy structure.     

Non-destructive estimation of potato leaf chlorophyll from canopy hyperspectral reflectance using the inverted PROSAIL model

Optimizing nitrogen (N) fertilization in crop production by in-season measurements of crop N status may improve fertilizer N use efficiency. Hyperspectral measurements may be used to assess crop N status by estimating leaf chlorophyll content. This study evaluated the ability of the PROSAIL canopy-level reflectance model to predict leaf chlorophyll content. Trials were conducted with two potato cultivars under different N fertility rates (0–300 kg N ha⁻¹). Canopy reflectance, leaf area index (LAI) and leaf chlorophyll and N contents were measured. The PROSAIL model was able to predict leaf chlorophyll content with reasonable accuracy later in the growing season. The low estimation accuracy earlier in the growing season could be due to model sensitivity to non-homogenous canopy architecture and soil background interference before full canopy closure. Canopy chlorophyll content (leaf chlorophyll content × LAI) was predicted less accurately than leaf chlrophyll content due to the low estimation accuracy of LAI for values higher than 4.5.     

双冠层反射率模型分析各向异性平整指数(AFX)对植被参数的敏感性

二向反射分布函数包含地表反射的方向性特征信息。研究二向反射分布函数BRDF(Bidirectional Reflectance Distribution Function)形状对植被结构参数的敏感性,有助于理解植被的二向性反射规律,进而反演植被参数。本文耦合双冠层反射率模型和核驱动的罗斯厚核-李氏稀疏互易核模型,利用EFAST全局敏感性分析方法,以各向异性平整指数为BRDF形状变化的衡量指标,研究了不同天空光比例(SKYL)下,各向异性平整指数AFX对植被参数敏感度的变化,以及SKYL=0.1时AFX的敏感性。结果表明:(1)在红波段,上、下层叶面积指数、上层叶绿素含量,以及上层叶倾角分布是AFX的敏感参数,在近红外波段,上、下层叶面积指数LAI是AFX的敏感参数。(2)冠层尺度上的参数敏感度总体大于叶片尺度。(3)晴天时(SKYL=0.1),红波段主敏感度较大的参数分别是上层LAI、上层叶倾角分布和下层叶片结构参数,近红外波段主敏感度较大的参数主要是上、下层LAI。     

Estimation of forest aboveground biomass from HJ1B imagery using a canopy reflectance model and a forest growth model

Accurately estimating the spatial distribution of forest aboveground biomass (AGB) is important because of its carbon budget forms part of the global carbon cycle. This paper presented three methods for obtaining forest AGB based on a forest growth model, a Multiple-Forward-Mode (MFM) method and a stochastic gradient boosting (SGB) model. A Li-Strahler geometric-optical canopy reflectance model (GOMS) with the ZELIG forest growth model was run using HJ1B imagery to derive forest AGB. GOMS-ZELIG simulated data were used to train the SGB model and AGB estimation. The GOMS-ZELIG AGB estimation was evaluated for 24 field-measured data and compared against the GOMS-SGB model and GOMS-MFM biomass predictions from multispectral HJ1B data. The results show that the estimation accuracy of the GOMS-MFM model is slightly higher than that of the GOMS-SGB model. The GOMS-ZELIG and GOMS-MFM models are considerably more accurate at estimating forest AGB in arid and semiarid regions.     

Predicting foliar biochemistry of tea (Camellia sinensis) using reflectance spectra measured at powder,leaf and canopy levels

Some biochemical compounds are closely related with the quality of tea (Camellia sinensis (L.)). In this study, the concentration of these compounds including total tea polyphenols, free amino acids and soluble sugars were estimated using reflectance spectroscopy at three different levels: powder, leaf and canopy, with partial least squares regression. The focus of this study is to systematically compare the accuracy of tea quality estimations based on spectroscopy at three different levels. At the powder level, the average r² between predictions and observations was 0.89 for polyphenols, 0.81 for amino acids and 0.78 for sugars, with relative root mean square errors (RMSE/mean) of 5.47%, 5.50% and 2.75%, respectively; at the leaf level, the average r² decreased to 0.46–0.81 and the relative RMSE increased to 4.46–7.09%. Compared to the results yielded at the leaf level, the results from canopy spectra were slightly more accurate, yielding average r² values of 0.83, 0.77 and 0.56 and relative RMSE of 6.79%, 5.73% and 4.03% for polyphenols, amino acids and sugars, respectively. We further identified wavelength channels that influenced the prediction model. For powder and leaves, some bands identified can be linked to the absorption features of chemicals of interest (1648 nm for phenolic, 1510 nm for amino acids, 2080 nm and 2270 nm for sugars), while more indirectly related wavelengths were found to be important at the canopy level for predictions of chemical compounds. Overall, the prediction accuracies achieved at canopy level in this study are encouraging for future study on tea quality estimated at the landscape scale using airborne and space-borne sensors.     

江苏省人口分布的时空动态分析

为了探寻江苏省人口分布的时空动态特征以提供人口动态管理的决策依据,该文以2000年第五次人口普查和2010年第六次人口普查数据为基础,应用GIS空间分析功能结合人口地理学方法,对江苏省两次普查的人口分布特征进行分析。结果表明,江苏省人口空间分布不均,存在着苏南、苏中、苏北以及市县之间的差异;人口空间分布不均衡的趋势不断增强,市域与县域人口分布特征存在差异,人口分布呈现出由苏北向苏南转移的特点。     

Use of field reflectance data for crop mapping using airborne hyperspectral image

Recent developments in hyperspectral remote sensing technologies enable acquisition of image with high spectral resolution, which is typical to the laboratory or in situ reflectance measurements. There has been an increasing interest in the utilization of in situ reference reflectance spectra for rapid and repeated mapping of various surface features. Here we examined the prospect of classifying airborne hyperspectral image using field reflectance spectra as the training data for crop mapping. Canopy level field reflectance measurements of some important agricultural crops, i.e. alfalfa, winter barley, winter rape, winter rye, and winter wheat collected during four consecutive growing seasons are used for the classification of a HyMAP image acquired for a separate location by (1) mixture tuned matched filtering (MTMF), (2) spectral feature fitting (SFF), and (3) spectral angle mapper (SAM) methods. In order to answer a general research question “what is the prospect of using independent reference reflectance spectra for image classification”, while focussing on the crop classification, the results indicate distinct aspects. On the one hand, field reflectance spectra of winter rape and alfalfa demonstrate excellent crop discrimination and spectral matching with the image across the growing seasons. On the other hand, significant spectral confusion detected among the winter barley, winter rye, and winter wheat rule out the possibility of existence of a meaningful spectral matching between field reflectance spectra and image. While supporting the current notion of “non-existence of characteristic reflectance spectral signatures for vegetation”, results indicate that there exist some crops whose spectral signatures are similar to characteristic spectral signatures with possibility of using them in image classification.     

Spectral reflectance patterns and seasonal dynamics of common understory types in three mature hemi-boreal forests

Due to the growing demand on more accurate prediction of biophysical properties (e.g., leaf area index) or carbon balance models based on remotely sensed data, the understory effect needs to be separated from the overstory. Reflectance models can provide possibility to model and retrieve understory reflectance over large scales, but ground truth data is needed to validate such models and algorithms. In this study, we documented the seasonal variation (April–September) and spectral changes occurring in understory layers of a typical European hemi-boreal forest. The understory composition was recorded and its spectra measured with an ASD FieldSpec Hand-Held UV/VNIR Spectroradiometer eight times at four site types during the growing period (from May to September) in 2013. The collected dataset presented within this study would be of much use to improve and validate algorithms or models for extracting spectral properties of understory from remote sensing data. It can be also further used as a valuable input in radiative transfer simulations that are used to quantify the roles of forest tree layer and understory components in forming a seasonal reflectance course of a hemi-boreal forest, and the upcoming phases of the RAdiation Model Intercomparison (RAMI) experiment.     

MODIS二向反射模型方差函数的研究与应用

地表二向性反射分布函数(BRDF)是表征地物反射随太阳和观测方向变化的物理量。在统计意义上,BRDF表示均值统计量,BRVF(Bidirectional Reflectance Variance Function)表示方差统计量,它们对研究地表各向异性反射特征有着重要意义。本文首先采用误差传播理论,推导出基于MODIS BRDF模型的BRVF表达形式。研究结果表明,BRVF的空间分布模式主要由几何光学核Kgeo和体散射核Kvol的一次项和二次项权重和决定。然后利用EOS地面验证核心站点(EOS Land Validation Core Sites)的MODIS BRDF产品,对BRVF空间分布模式随地表类型、光谱波段和观测角度范围进行验证。验证结果表明,基于MODIS BRDF产品的验证与理论推导有较好的一致性。BRVF空间分布模式和地表类型有关,通常在热点处有一个峰值。在大观测天顶角(60°)下,BRVF随着角度的增大而增大。BRVF在近红外波段整体上大于红波段,表明其波段依赖性。最后,将上述理论成果初步应用于69组地表测量数据的模拟中。模拟结果表明,当大角度缺少观测数据时,模型外延所引起的方向反射方差显著增大,对地表反照率的反演精度和不确定性有较大影响。其中,红波段的白天空反照率的相对误差最大可达38.26%。本研究对利用小角度观测数据进行地表反照率反演的不确定性分析有指导意义;对大角度观测数据缺失情况下,先验知识在地表反照率的反演应用可提供有意义的理论支撑。