土壤背景对冠层NDVI的影响分析

归一化差值植被指数NDVI是植被遥感中应用最为广泛的指数之一, 但它受土壤背景等因素的干扰比较强烈。结合实测的土壤数据以及公式推导、PROSAIL 模型模拟等方法分析了这种影响。首先, 假定与土壤线性混合且叶片呈水平分布的植被冠层, 根据土壤与植被分别在红光、近红外波段处的反射率值、植被覆盖度等参数, 利用公式推导了土壤背景对不同覆盖度下冠层NDVI的影响。其次, 利用PROSAIL冠层光谱模拟模型, 模拟分析了土壤背景对不同LAI下冠层NDVI的影响。分析的结果表明:LAI 越小, 土壤背景的影响越大; 暗土壤背景下的冠层NDVI值大于亮土壤背景下冠层的NDVI值; 并且,暗土壤条件下,NDVI值对土壤亮度的变化更敏感,而亮土壤下,NDVI值则对LAI或覆盖度的变化更敏感。最后利用实测的不同土壤背景下的冬小麦冠层光谱数据, 验证了公式推导和模型模拟的结果。     

基于无人机高光谱数据的玉米叶面积指数估算

无人机高光谱遥感是低成本、高精度获取精细尺度农作物生物物理参数和生物化学参数的新型手段,以此快速反演叶面积指数(Leaf Area Index,LAI)对作物长势评价、产量预测具有重要意义。以山东禹城市玉米为研究对象,利用PROSAIL辐射传输模型模拟玉米冠层反射率获取LAI特征响应波段结合相关性定量分析获取对LAI变化最为敏感的波段,并以此计算6种植被指数(Vegetation Index,VI),利用6种回归模型分别对单一特征波段和VI进行反演建模,以实测LAI评定模型精度。研究表明,光谱反射率中516、636、702、760和867 nm等波段对LAI变化最为敏感,以此建立的单一特征波段反演模型预测LAI精度R~2为0.44～0.58;RMSE为0.16～0.18,其中636 nm建立的模型(LAI=21.86exp(-29.47R636))相比其他反演模型预测精度较高(R~2=0.58,RMSE=0.16);6种植被指数与LAI高度相关,相关性系数R2为0.85～0.86,以此建立的反演模型相比单一特征波段反演模型精度有所提高,R~2为0.66～0.72,RMSE为0.12～0.14;其中mNDVI构建的LAI估算模型(LAI=exp(2.76～1.77/mNDVI))精度最高(R~2=0.72,RMSE=0.13)。无人机高光谱遥感是快速、无损监测农作物生长信息的有效手段,为指导精细化尺度作物管理提供依据。     

基于无人机高光谱数据的玉米叶面积指数估算

无人机高光谱遥感是低成本、高精度获取精细尺度农作物生物物理参数和生物化学参数的新型手段,以此快速反演叶面积指数（Leaf Area Index, LAI）对作物长势评价、产量预测具有重要意义。以山东禹城市玉米为研究对象,利用PROSAIL辐射传输模型模拟玉米冠层反射率获取LAI特征响应波段结合相关性定量分析获取对LAI变化最为敏感的波段,并以此计算6种植被指数（Vegetation Index,VI）,利用6种回归模型分别对单一特征波段和VI进行反演建模,以实测LAI评定模型精度。研究表明,光谱反射率中516、636、702、760和867 nm等波段对LAI变化最为敏感,以此建立的单一特征波段反演模型预测LAI精度R2为0.44~0.58;RMSE为0.16~0.18,其中636 nm建立的模型(LAI=21.86exp(-29.47R₆₃₆))相比其他反演模型预测精度较高（R2=0.58,RMSE=0.16）;6种植被指数与LAI高度相关,相关性系数R ²为0.85~0.86,以此建立的反演模型相比单一特征波段反演模型精度有所提高,R2为0.66~0.72,RMSE为0.12~0.14;其中mNDVI构建的LAI估算模型（LAI=exp(2.76~1.77/mNDVI))精度最高（R2=0.72,RMSE=0.13）。无人机高光谱遥感是快速、无损监测农作物生长信息的有效手段,为指导精细化尺度作物管理提供依据。     

基于无人机高光谱数据的玉米叶面积指数估算

无人机高光谱遥感是低成本、高精度获取精细尺度农作物生物物理参数和生物化学参数的新型手段，以此快速反演叶面积指数（Leaf Area Index, LAI）对作物长势评价、产量预测具有重要意义。以山东禹城市玉米为研究对象，利用PROSAIL辐射传输模型模拟玉米冠层反射率获取LAI特征响应波段结合相关性定量分析获取对LAI变化最为敏感的波段，并以此计算6种植被指数（Vegetation Index,VI），利用6种回归模型分别对单一特征波段和VI进行反演建模，以实测LAI评定模型精度。研究表明，光谱反射率中516、636、702、760和867 nm等波段对LAI变化最为敏感，以此建立的单一特征波段反演模型预测LAI精度R2为0.44~0.58；RMSE为0.16~0.18，其中636 nm建立的模型(LAI=21.86exp(-29.47R₆₃₆))相比其他反演模型预测精度较高（R2=0.58，RMSE=0.16）；6种植被指数与LAI高度相关，相关性系数R ²为0.85~0.86，以此建立的反演模型相比单一特征波段反演模型精度有所提高，R2为0.66~0.72，RMSE为0.12~0.14；其中mNDVI构建的LAI估算模型（LAI=exp(2.76~1.77/mNDVI))精度最高（R2=0.72，RMSE=0.13）。无人机高光谱遥感是快速、无损监测农作物生长信息的有效手段，为指导精细化尺度作物管理提供依据。     

全球典型植被叶片光谱特征及其对LAI反演的影响分析

在全球范围长时间序列LAI遥感产品反演算法中,植被冠层反射率模型仅使用少量叶片光谱特征代表全球植被全年的典型植被光谱特征,叶片光谱的不确定性导致LAI遥感产品存在一定的误差。目前全球已经构建了多个典型植被叶片波谱数据集,这些数据集包含多个植被物种、不同空间地域及多时相叶片光谱数据,为定量分析叶片光谱特征提供了数据支持。主要利用LOPEX’93、ANGERS’03、中国典型地物波谱数据库和野外实测的叶片光谱数据,以黄边参数、红边参数和叶片光谱指数作为分析指标,探讨不同植被物种、不同气候区和不同物候期的叶片光谱特征差异,及其对植被冠层反射率、LAI反演的影响,为发展考虑现实叶片光谱差异的LAI反演算法提供研究基础。结果表明:植被叶片光谱存在多样性,叶片光谱特征差异主要影响MODIS传感器近红外波段和绿波段反射率值,其中,绿波段反射率值对叶片光谱变化最为敏感;在LAI反演算法中,如果只考虑植被类型而不考虑物种叶片光谱差异,可能会给LAI反演带来大于3的误差。     

波段位置和宽度对河口湿地4种植被NDVI的影响

研究不同波段位置和宽度对植被NDVI的影响,对于进一步认识NDVI指数具有重要的意义。采用ASD(Analytical Spectral Devices)地物光谱仪测定闽江河口互花米草(Spartina alterniflora)、秋茄(Kandelia candel)、芦苇(Phragmites australis)和短叶茳芏(Cyperus malaccensis)冠层光谱,利用ViewSpecPro和Oragin8.0软件对光谱数据进行分析和处理,探讨不同波段位置和波段宽度对河口湿地4种植被NDVI的影响。结果表明:①当红光波段固定,近红外波段以50 nm宽度移动时,4种湿地植被NDVI随近红外波段中心位置增加而迅速增加,之后趋于平稳,在925~1 050 nm出现一个小的谷值,互花米草和短叶茳芏的谷值要比其他两种植物更为明显;不同波段宽度影响表现为:除红边与970 nm附近区域对NDVI的影响较显著外,其他波段影响不显著;②当近红外波段固定,红光波段以10 nm宽度移动时,4种湿地植被NDVI随红光波段中心位置移动先略有增加或变化不大,然后迅速降低;不同波段宽度影响表现为:在650~700 nm波段宽度越宽,NDVI值越小,600~650 nm范围内波段宽度对NDVI的影响不大;③4种湿地植被红光波段宽度对NDVI的影响要大于近红外波段。     

双层植被结构冠层光谱特性的理论模拟

利用双层冠层反射率模型ACRM，模拟不同叶面积指数LAI、含水量Cw和结构参数N下，波长是820nm和1600nm冠层反射率的角度分布。结果说明，该模型能准确模拟出“热点”效应 |冠层反射率角度分布对LAI的敏感性小于LAI和Cw以及LAI和N的共同作用。其中LAI和Cw共同作用对波长1 600 nm反射率角度分布的影响非常显著，而LAI和N的共同作用在820 nm略微大于1 600 nm。另外，提取冠层含水量的土壤可调节水分指数SAWI受冠层结构的影响也较大。今后在模型选取中应该更好的考虑冠层结构影响。     

基于光谱归一化的马尾松LAI遥感估算研究

通过对比不同传感器间光谱响应函数的差异,研究基于光谱响应函数的不同传感器相似波段的归一化方法,探讨归一化后植被指数在马尾松叶面积指数（LAI）估算中的应用。以某一传感器为基准,根据波段总辐射率比值关系将其他卫星传感器归一化为基准传感器,然后计算其植被指数,建立LAI反演模型。为验证方法可行性,选取永安地区2008年3月获取的BJ-1CCD、IRS-P6LISS3和MODIS数据作为研究对象,根据三者的光谱响应函数差异,将BJ-1CCD和IRS-P6的LISS3的红光和近红外波段归一化为MODIS的相应波段,并分别计算归一化前后的NDVI值。结果表明归一化后不同传感器的植被指数关系与理想的关系y=x更加接近。利用归一化后的IRS-P6影像的NDVI反演马尾松LAI,并将其应用于MODIS和BJ-1传感器,得到归一化后不同传感器的植被指数值基本相等,表明归一化以后的植被指数应用于LAI的估算具有一定的普适性,能适用于多种传感器。     

玉米叶面积指数与高光谱植被指数关系研究

探讨以不同的植被指数建立的高光谱模型对玉米叶面积指数LAI的反演精度。实测不同水肥耦合作用下，玉米冠层的高光谱反射率与叶面积指数（Leaf Area Index）数据，采用高光谱红光波段（631~760 nm）与近红外波段（760~1 074 nm）逐波段构建NDVI、RVI、DVI、TSAVI、PVI植被指数，分别找出与LAI具有最佳相关性波段组合的植被指数，建立玉米LAI估算模型。结果显示，与LAI具有佳相关性的波段组合分别是NDVI（R760，R990）、RVI（R760，R1001）、DVI（R677，R1070）、TSAVI（R 760，R 975）、PVI（R658，R966），它们反演玉米LAI的确定性系数分别：R²>0.72、R²>0.74、R²＝0.95、R²>0.79、R²>0.95。结果表明，在玉米的整个生长季的47个样本中，通过PVI和DVI方式建立的遥感估算模型能够较为准确地估算玉米LAI，TSAVI次之，NDVI、RVI稍差。     

水稻叶面积指数的多光谱遥感估算模型研究

LAI是生态系统研究中最重要的结构参数之一,它是估计多种植冠功能过程的重要参数。通过两年的水稻田间试验,使用美国ASD背挂式野外光谱辐射仪(ASDFieldSpec),获取1999～2000年两年晚稻整个生育期的光谱数据,采用计算机测算图斑面积法测定LAI;根据已有的卫星传感器通道波段(MSS、RBV、SPOT、TM、CH)和它们的组合(比值植被指数、归一化差植被指数),以及具有物理意义的光谱区域(蓝区、绿区、黄边、红光吸收谷、红边、紫区、可见光区、近红外区、全部波段)等共有27个变量构建多光谱变量组,采用5个单变量线性与非线性拟合模型,用1999年试验数据为训练样本,建立水稻LAI的多光谱遥感估算模型。结果表明:适用于水稻LAI估算的多光谱变量是植被指数变量好于波段变量;RVI与NDVI比较,RVI好于NDVI。用2000年试验数据作为测试样本数据,对其精度进行评价和验证,非线性模型的精度高于线性模型的精度,其中以SPOT3/SPOT2为变量的对数模型,拟合R2与预测R2达到了最大,其RMSE和相对误差(%)为最低,因此,认为它是估算LAI的最佳模型。
     

Non-destructive estimation of wheat leaf chlorophyll content from hyperspectral measurements through analytical model inversion

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 indirectly by estimating leaf and canopy chlorophyll content. This study evaluated the ability of the PROSAIL canopy-level reflectance model to predict leaf chlorophyll content of spring wheat (Triticum aestivum L.) during the growth stages between pre-tillering (Zadoks Growth Stage (ZGS 15)) to booting (ZGS50). Spring wheat was grown under different N fertility rates (0–200 kg N ha^?1) in 2002. Canopy reflectance, leaf chlorophyll content, N content and leaf area index (LAI) values were measured. There was a weakly significant trend for the PROSAIL model to over-estimate LAI and under-estimate leaf chlorophyll content. To compensate for this interdependency by the model, a canopy chlorophyll content parameter (the product of leaf chlorophyll content and LAI) was calculated. The estimation accuracy for canopy chlorophyll content was generally low earlier in the growing season. This failure of the PROSAIL model to estimate leaf and canopy variables could be attributed to model sensitivity to canopy architecture. Earlier in the growing season, full canopy closure was not yet achieved, resulting in a non-homogenous canopy and strong soil background interference. The canopy chlorophyll content parameter was predicted more accurately than leaf chlorophyll content alone at booting (ZGS 45). A strong relationship between canopy chlorophyll content and canopy N content at ZGS 45 indicates that the PROSAIL model may be used as a tool to predict wheat N status from canopy reflectance measurements at booting or later.     

Mapping daily leaf area index at 30 m resolution over a meadow steppe area by fusing Landsat,Sentinel-2A and MODIS data

ABSTRACT

The leaf area index (LAI) is a key vegetation canopy structure parameter and is closely associated with vegetation photosynthesis, transpiration, and energy balance. Developing a landscape-scale LAI dataset with a high temporal resolution (daily) is essential for capturing rapidly changing vegetation structure at field scales and supporting regional biophysical modeling efforts. In this study, two daily 30 m LAI time series from 2014 to 2016 over a meadow steppe site in northern China were generated using a spatial and temporal adaptive re?ectance fusion model (STARFM) combined with an LAI retrieval radiative transfer model (PROSAIL). Gap-filled Landsat 7, Landsat 8 and Sentinel-2A surface reflectance (SR) images were used to generate fine-resolution LAI maps with the PROSAIL look-up table method. Two daily 500 m moderate-resolution imaging spectroradiometer (MODIS) LAI product-the existing MCD15A3H LAI product and one was generated from the MCD43A4 SR product and the PROSAIL model, were used to provide temporally continuous LAI variations. The STARFM model was then used to fuse the fine-resolution LAI maps with the two 500 m LAI products separately to generate two daily 30 m LAI time series. Both results were assessed for three types of pasture (mowed pasture, grazing pasture, and fenced pasture) using ground measurements from 2014–2015. The results showed that the PROSAIL-generated LAI maps all exhibited a high accuracy, and the root mean squared errors (RMSEs) for the Landsat 7 LAI and Landsat 8 LAI compared to the ground-measured LAI were 0.33 and 0.28 respectively. The Landsat LAI maps also showed good agreement and similar spatial patterns with the Sentinel-2A LAI with mean differences between ± 0.5. The MCD43A4_PROSPECT LAI product exhibited similar seasonal variability to the ground measurements and to the Landsat and Sentinel-2A LAIs, and these data are also smoother and contain fewer noisy points than the gap-filled MCD15A3H LAI product. Compared to the ground measurements, the daily 30 m LAI time series fused from the fine-resolution LAI maps and PROSPECT generated MODIS LAI product demonstrated better performance with an RMSE of 0.44 and a mean absolute error (MAE) of 0.34, which is an improvement from the LAI time series fused from the fine-resolution LAI maps and the existing MCD15A3H LAI product (RMSE of 0.56 and MAE of 0.42). The latter dataset also exhibited abnormal temporal fluctuations, which may have been caused by the interpolation method. The results also demonstrated the very good performance of the STARFM model in grazing and mowed pasture with homogeneous surfaces compared to fenced pasture with smaller patch sizes. The Sentinel-2A data offers increased landscape vegetation observation frequency and provides temporal information about canopy changes that occur between Landsat overpass dates. The scheme developed in this study can be used as a reference for regional vegetation dynamic studies and can be applied to larger areas to improve grassland modeling efforts.     

Global sensitivity analysis of PROSAIL model parameters when simulating Moso bamboo forest canopy reflectance

PROSAIL is a combination of the leaf optical properties spectra (PROSPECT) model and the scattering by arbitrarily inclined leaves (SAIL) canopy bidirectional reflectance model. When modelling forest canopy reflectance using the PROSAIL radiative transfer model, the sensitivities of parameters can affect the modelling accuracy. Traditionally, sensitivities have been assessed using local sensitivity analysis (LSA); however, drawbacks to this approach include a lack of consideration for coupled effects between different parameters. In this study, parameter sensitivities in the PROSAIL model were calculated using two global sensitivity analysis (GSA) methods (the Extended Fourier Amplitude Sensitivity Test (EFAST) method and the Morris method), field measurements, and Landsat 5 Thematic Mapper (TM) data for a Moso bamboo forest. The results of GSA were compared with those of LSA in order to identify the key parameters impacting the Moso bamboo forest canopy reflectance, and to provide a reference for model optimization and vegetation canopy inversion improvement. The results showed that: (1) the sensitivities of six major input parameters of the PROSAIL model were generally consistent with the sorting orders of the two GSA methods, but were not in accordance with those from the LSA method, especially in the mid-infrared band; (2) coupled effects among parameters acting on reflectance simulation in visible light bands were greater than those in infrared bands; (3) the simulated canopy reflectance was evaluated using Landsat 5 TM data, and the results simulated based on LSA analysis showed higher error than those based on GSA analysis, because the LSA method ignored the influence of some parameters on canopy reflectance, e.g. leaf mesophyll structure (N), average leaf angle (ALA), leaf water content (C_w), and leaf dry matter content (C_m). However, GSA was able to fully consider the coupled effects among parameters, and thus identified the sensitive parameters impacting on reflectance more accurately.     

Inversion of a radiative transfer model for estimating vegetation LAI and chlorophyll in a heterogeneous grassland

Radiative transfer models have seldom been applied for studying heterogeneous grassland canopies. Here, the potential of radiative transfer modeling to predict LAI and leaf and canopy chlorophyll contents in a heterogeneous Mediterranean grassland is investigated. The widely used PROSAIL model was inverted with canopy spectral reflectance measurements by means of a look-up table (LUT). Canopy spectral measurements were acquired in the field using a GER 3700 spectroradiometer, along with simultaneous in situ measurements of LAI and leaf chlorophyll content. We tested the impact of using multiple solutions, stratification (according to species richness), and spectral subsetting on parameter retrieval. To assess the performance of the model inversion, the normalized RMSE and R² between independent in situ measurements and estimated parameters were used. Of the three investigated plant characteristics, canopy chlorophyll content was estimated with the highest accuracy (R² = 0.70, NRMSE = 0.18). Leaf chlorophyll content, on the other hand, could not be estimated with acceptable accuracy, while LAI was estimated with intermediate accuracy (R² = 0.59, NRMSE = 0.18). When only sample plots with up to two species were considered (n = 107), the estimation accuracy for all investigated variables (LAI, canopy chlorophyll content and leaf chlorophyll content) increased (NRMSE = 0.14, 0.16, 0.19, respectively). This shows the limits of the PROSAIL radiative transfer model in the case of very heterogeneous conditions. We also found that a carefully selected spectral subset contains sufficient information for a successful model inversion. Our results confirm the potential of model inversion for estimating vegetation biophysical parameters at the canopy scale in (moderately) heterogeneous grasslands using hyperspectral measurements.     

基于宽波段和窄波段植被指数的草地LAI反演对比研究

叶面积指数是一个重要的植被生理生态参数,为探讨不同植被指数反演叶面积指数的可行性,基于同空间分辨率不同光谱分辨率的HJ\|1B CCD1和Hyperion遥感影像数据,以内蒙古自治区赤峰市克斯克腾旗贡格尔草原为研究对象,选取几种常见宽波段植被指数和高光谱窄波段植被指数并结合4种常用回归模型,比较分析了不同植被指数反演叶面积指数的精度。结果表明:对于全部植被指数而言,PVI、MSAVI等综合考虑了土壤、环境等因素的植被指数较传统植被指数NDVI、RVI反演草地LAI精度更高。通过对比发现,在反演草地LAI方面,窄波段植被指数比宽波段植被指数表现出明显的优势。其中,窄波段垂直植被指数PVI验证模型的确定性系数R²为0.65,均方根误差RMSE为0.15,说明实测LAI和模拟LAI值之间具有较好的变化一致性。最后基于Hyperion影像和窄波段垂直植被指数PVI的估算模型生成研究区叶面积指数空间分布图。     

Estimating winter wheat (Triticum aestivum) LAI and leaf chlorophyll content from canopy reflectance data by integrating agronomic prior knowledge with the PROSAIL model

Leaf area index (LAI) and leaf chlorophyll content (LCC) are major considerations in management decisions, agricultural planning, and policy-making. When a radiative transfer model (RTM) was used to retrieve these biophysical variables from remote-sensing data, the ill-posed problem was unavoidable. In this study, we focused on the use of agronomic prior knowledge (APK), constructing the relationship between LAI and LCC, to restrict and mitigate the ill-posed inversion results. For this purpose, the inversion results obtained using the SAILH+PROSPECT (PROSAIL) canopy reflectance model alone (no agronomic prior knowledge, NAPK) and those linked with APK were compared. The results showed that LAI inversion had high accuracy. The validation results of the root mean square error (RMSE) between measured and estimated LAI were 0.74 and 0.69 for NAPK and APK, respectively. Compared with NAPK, APK improved LCC estimation; the corresponding RMSE values of NAPK and APK were 13.36 µg cm^–2 and 9.35 µg cm^–2, respectively. Our analysis confirms the operational potential of PROSAIL model inversion for the retrieval of biophysical variables by integrating APK.     

Normalized sensitivity measures for leaf area index estimation using three-band spectral vegetation indices

A great number of spectral vegetation indices (SVIs) have been developed to estimate key biophysical parameters such as leaf area index (LAI). Considerable interest is often given to the local optimization, performance analysis and sensitivity of each spectral band and SVI for LAI estimation given that several confounding factors are present. In this regard, inclusion of shortwave infrared (SWIR) reflectance in traditionally near-infrared (NIR)-red (R)-based SVIs has played a great role for local optimization and increased sensitivity of SVIs to LAI. This study presents the enhanced and normalized sensitivity functions for evaluating (1) the sensitivity of each spectral band and SVI to LAI and (2) the generic performance analysis of empirical model to estimate LAI based on the SVIs. Several alternatives for three-band (NIR-R-SWIR) SVI modifications have been recommended and proven to be simplistic and unbiased way of local optimization.     

A Bayesian network algorithm for retrieving the characterization of land surface vegetation

A hybrid inversion technique based on Bayesian network is proposed for estimating the biochemical and biophysical parameters of land surface vegetation from remotely sensed data. A Bayesian network is a unified knowledge-inferring process that can incorporate information derived from multiple sources including remote sensing and information derived from a priori knowledge. Using this inversion approach, content of chlorophyll a and chlorophyll b (Cab) and leaf area index (LAI) of winter wheat were estimated from data derived from simulations as well as field measurements. Estimations from the simulated data proved accurate, with root mean square errors (RMSEs) of 0.54 m²/m² in LAI and 4.5 μg/cm² in Cab. In validating the estimates against field measurements, it was found that prior knowledge of target parameters improved the accuracy of estimates, in terms of RMSEs from 0.73 to 0.22 m²/m² in LAI and 9.6 to 4.0 μg/cm² in Cab. Bayesian inference in this hybrid inversion scheme produces a posterior probability distribution, which can reveal such properties of the inferred results as updated information contained in the inversion result. Using entropy, the revision of posterior information about the parameters of interest was calculated. Including more data may allow more information to be retrieved about parameters in general. Exceptions were also observed where data from some viewing angles slightly reduced the information on the parameters of interest. It was also found that data from these viewing angles were less sensitive to the parameters. The method proposed here was also validated using LandSat ETM+ imagery provided by the BigFoot project. When used for mapping LAI with ETM+ imagery, the proposed method with an RMSE of 0.70 and a correlation of 0.67 produced a slightly better result than that from empirical regression.