Estimating chlorophyll content of crops from hyperspectral data using a normalized area over reflectance curve (NAOC)

The Normalized Area Over reflectance Curve (NAOC) is proposed as a new index for remote sensing estimation of the leaf chlorophyll content of heterogeneous areas with different crops, different canopies and different types of bare soil. This index is based on the calculation of the area over the reflectance curve obtained by high spectral resolution reflectance measurements, determined, from the integral of the red–near-infrared interval, divided by the maximum reflectance in that spectral region. For this, use has been made of the experimental data of the SPARC campaigns, where in situ measurements were made of leaf chlorophyll content, LAI and fCOVER of 9 different crops – thus, yielding 300 different values with broad variability of these biophysical parameters. In addition, Proba/CHRIS hyperspectral images were obtained simultaneously to the ground measurements. By comparing the spectra of each pixel with its experimental leaf chlorophyll value, the NAOC was proven to exhibit a linear correlation to chlorophyll content. Calculating the correlation between these variables in the 600–800 nm interval, the best correlation was obtained by computing the integral of the spectral reflectance curve between 643 and 795 nm, which practically covers the spectral range of maximum chlorophyll absorption (at around 670 nm) and maximum leaf reflectance in the infrared (750–800 nm). Based on a Proba/CHRIS image, a chlorophyll map was generated using NAOC and compared with the land-use (crops classification) map. The method yielded a leaf chlorophyll content map of the study area, comprising a large heterogeneous zone. An analysis was made to determine whether the method also serves to estimate the total chlorophyll content of a canopy, multiplying the leaf chlorophyll content by the LAI. To validate the method, use was made of the data from another campaign ((SEN2FLEX), in which measurements were made of different biophysical parameters of 7 crops, and hyperspectral images were obtained with the CASI imaging radiometer from an aircraft. Applying the method to a CASI image, a map of leaf chlorophyll content was obtained, which on, establishing comparisons with the experimental data allowed us to estimate chlorophyll with a root mean square error of 4.2 μg/cm², similar or smaller than other methods but with the improvement of applicability to a large set of different crop types.     

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

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

Evaluating the capabilities of Sentinel-2 for quantitative estimation of biophysical variables in vegetation

The red edge position (REP) in the vegetation spectral reflectance is a surrogate measure of vegetation chlorophyll content, and hence can be used to monitor the health and function of vegetation. The Multi-Spectral Instrument (MSI) aboard the future ESA Sentinel-2 (S-2) satellite will provide the opportunity for estimation of the REP at much higher spatial resolution (20 m) than has been previously possible with spaceborne sensors such as Medium Resolution Imaging Spectrometer (MERIS) aboard ENVISAT. This study aims to evaluate the potential of S-2 MSI sensor for estimation of canopy chlorophyll content, leaf area index (LAI) and leaf chlorophyll concentration (LCC) using data from multiple field campaigns. Included in the assessed field campaigns are results from SEN3Exp in Barrax, Spain composed of 35 elementary sampling units (ESUs) of LCC and LAI which have been assessed for correlation with simulated MSI data using a CASI airborne imaging spectrometer. Analysis also presents results from SicilyS2EVAL, a campaign consisting of 25 ESUs in Sicily, Italy supported by a simultaneous Specim Aisa-Eagle data acquisition. In addition, these results were compared to outputs from the PROSAIL model for similar values of biophysical variables in the ESUs. The paper in turn assessed the scope of S-2 for retrieval of biophysical variables using these combined datasets through investigating the performance of the relevant Vegetation Indices (VIs) as well as presenting the novel Inverted Red-Edge Chlorophyll Index (IRECI) and Sentinel-2 Red-Edge Position (S2REP). Results indicated significant relationships between both canopy chlorophyll content and LAI for simulated MSI data using IRECI or the Normalised Difference Vegetation Index (NDVI) while S2REP and the MERIS Terrestrial Chlorophyll Index (MTCI) were found to have the strongest correlation for retrieval of LCC.     

不同传感器的模拟植被指数对水稻叶面积指数的估测精度和敏感性分析

本研究利用水稻冠层高光谱数据,模拟NOAA-AVHRR,Terra-MODIS和Landsat-TM的可见光波段反射率数据,计算各传感器的多种植被指数(NDVI,RVI,EVI,GNDVI,GRVI和Red-edge RVI),比较植被指数模型对水稻LAI的估测精度,分析不同植被指数对LAI变化的敏感性.相对于红波段植被指数,红边比值植被指数(Red-edge RVI)和绿波段指数GRVI与LAI有更好的线性相关关系,而GNDVI和LAI呈现更好的对数相关关系.MODIS的Red-edge RVI指数不仅模型拟合的精度最高,还有独立数据验证的估测精度也最高,而且它的验证精度较拟合精度下降幅度最小;其次是绿波段构建的GNDVI和GRVI植被指数的估测精度,再次是NDVI和EVI的估测精度,而RVI的估测精度最差.敏感性分析发现,13个植被指数对水稻LAI的估测能力都随着LAI的增加而下降,但归一化类植被指数和比值类植被指数对LAI变化反应的差异明显,归一化类植被指数在LAI较低时(LAI<1.5)对LAI变化的反应开始非常敏感,但迅速下降,而比值类植被指数在LAI较低时,明显小于归一化类植被指数,之后随着LAI的增大(LAI>1.5)比值类植被指数对LAI的变化敏感性,则明显高于归一化类植被指数.Red-edge RVI和绿波段指数GRVI和LAI不仅表现了很好的线性相关关系,而且在LAI大于2.9左右保持较高的敏感性.     

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.     

Principal component inversion technique for the retrieval of leaf area index

A simulation study has been carried out to investigate the Principal Component Inversion (PCI) technique for the retrieval of leaf area index (LAI). The PROSAIL model has been used for the forward analysis, i.e., estimation of reflectance for various combinations of LAI, soil reflectance, leaf angle distribution (ø₁), chlorophyll a+b concentration (C_ab), etc. Independent test on sample with LAI range 0.1-7.0 indicated that the retrieved LAI from PCI has higher accuracy (RMSE=0.137) than the classical NDVI-LAI empirical approach (RMSE=1.139). The study needs to be extended to cover retrieval from different types of soil and simultaneous retrieval of different biophysical parameters viz., LAI, C_ab, ø₁ to test the wider applicability of the PCI technique.     

Estimating salinity stress in sugarcane fields with spaceborne hyperspectral vegetation indices

The presence of salt in the soil profile negatively affects the growth and development of vegetation. As a result, the spectral reflectance of vegetation canopies varies for different salinity levels. This research was conducted to (1) investigate the capability of satellite-based hyperspectral vegetation indices (VIs) for estimating soil salinity in agricultural fields, (2) evaluate the performance of 21 existing VIs and (3) develop new VIs based on a combination of wavelengths sensitive for multiple stresses and find the best one for estimating soil salinity. For this purpose a Hyperion image of September 2, 2010, and data on soil salinity at 108 locations in sugarcane (Saccharum officina L.) fields were used. Results show that soil salinity could well be estimated by some of these VIs. Indices related to chlorophyll absorption bands or based on a combination of chlorophyll and water absorption bands had the highest correlation with soil salinity. In contrast, indices that are only based on water absorption bands had low to medium correlations, while indices that use only visible bands did not perform well. From the investigated indices the optimized soil-adjusted vegetation index (OSAVI) had the strongest relationship (R² = 0.69) with soil salinity for the training data, but it did not perform well in the validation phase. The validation procedure showed that the new salinity and water stress indices (SWSI) implemented in this study (SWSI-1, SWSI-2, SWSI-3) and the Vogelmann red edge index yielded the best results for estimating soil salinity for independent fields with root mean square errors of 1.14, 1.15, 1.17 and 1.15 dS/m, respectively. Our results show that soil salinity could be estimated by satellite-based hyperspectral VIs, but validation of obtained models for independent data is essential for selecting the best model.     

多光谱多角度遥感数据综合反演叶面积指数方法研究

叶面积指数是陆地生态系统的一个十分重要的结构参数。用遥感数据求取叶面积指数可以利用光谱的信息,比如通过植被指数来拟合一个经验关系,但很多植被指数明显受土壤背景的影响,对于有明显行结构的农作物,土壤的影响很难消除,植被指数的方法误差较大。多角度遥感包含了大量的地面目标的立体结构信息,具备求解植被特征参数的潜力,但通常多角度遥感反演对光谱信息的利用不足。与以往的反演方法相区别,该文利用行播作物二向反射模型,将多角度与多光谱数据结合进行行播作物LAI反演实验,并对反演算法进行了详细的敏感性分析实验,结果表明采用多角度、多光谱遥感数据相结合的方法可以有效反演行播作物的叶面积指数。     

Spatially Explicit Estimation of Leaf Area Index Using EO-1 Hyperion and Landsat ETM+ Data: Implications of Spectral Bandwidth and Shortwave Infrared Data on Prediction Accuracy in a Tropical Montane Environment

This study evaluated the utility of narrowband (EO-1 Hyperion) and broadband (Landsat ETM+) remote sensing data for the estimation of leaf area index (LAI) in a tropical environment in Sulawesi, Indonesia. LAI was inferred from canopy gap fraction measurements taken in natural tropical forest and cocoa plantations. Single and multiple spectral bands and spectral indices were used as predictor variables in reduced major axis (RMA) and ordinary least squares (OLS) regression models. The predictive power of most regression models was notably higher when employing narrowband data instead of broadband data. Highly significant relationships between LAI and spectral reflectance were observed near the red-edge region and in most shortwave infrared (SWIR) bands. In contrast to most near-infrared (NIR) narrow bands, the correlation between SWIR reflectance and LAI was not confounded when including both vegetation types and did not suffer from saturation. The results demonstrate that leaf area index of a challenging tropical environment can be estimated with satisfactory accuracy from hyperspectral remote sensing data.     

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.     

Inversion of PROSAIL Model for Retrieval of Plant Biophysical Parameters

The current development of satellite technology particularly in the sensors like POLDER and MISR, has emphasized more on directional reflectance measurements (i.e. spectral reflectance of the target measured from different view zenith and azimuth angles) of the earth surface features mainly the vegetation for retrieval of biophysical parameters at regional scale using radiative transfer models. This approach being physical process based and uses directional reflectance measurement has been found to better and more reliable compared to the conventional statistical approach used till date and takes care of anisotropic nature (i.e. reflectance from the target is different if measured from different view angles) of the target. Keeping this in view a field experiment was conducted in mustard crop to evaluate the radiative transfer model for biophysical parameter retrieval through its inversion with the objectives set as (i) to relate canopy biophysical parameters and geometry to its bidirectional reflectance, (ii) to evaluate a canopy reflectance model to best represent the radiative transfer within the canopy for its inversion and (iii) to retrieve crop biophysical parameters through inversion of the model. Two varieties of the mustard crop (Brassica juncea L) were grown with two nitrogen treatments. The bidirectional reflectance data obtained at 5 nm interval for a range of 400–1100 nm were integrated to IRS LISS–II sensor’s four band values using Newton Cotes Integration technique. Biophysical parameters like leaf area index, leaf chlorophyll content, leaf length, plant height and average leaf inclination angle, biomass etc were estimated synchronizing with the bi-directional reflectance measurements. Radiative transfer model PROSAIL model was validated and its inversion was done to retrieve LAI and ALA. Look Up Table (LUT) of Bidirectional reflectance distribution function (BRDF) was prepared simulating through PROSAIL model varying only LAI (0.2 interval from 1.2 to 5.4 ) and ALA (5° interval from 40° to 55°) parameters and inversion was done using a merit function and numerical optimization technique given by Press et al. (1986). The derived LAI and ALA values from inversion were well matched with observed one with RMSE 0.521 and 5.57, respectively.     

Estimation of gross primary production in wheat from in situ measurements

Gross primary production (GPP) is a parameter of significant importance for carbon cycle and climate change research. Remote sensing combined with other climate and meteorological data offers a convenient tool for large scale GPP estimation. This paper presents a study of GPP estimation using three methods with in situ measurements of canopy reflectance, LAI, and the photosynthetically active radiation (PAR). First, because LAI is considered as an indicator of the factor of absorbed PAR (f_APAR), it provides reasonable estimates of GPP for all types of wheat with coefficient of determination R² of 0.7353. The second method uses four kinds of vegetation indices (VIs) to estimate GPP because these indices are suggested to be reliable candidates in the estimation of light use efficiency (LUE). Good determination coefficients were acquired in estimating GPP with R² ranging from the lowest of 0.7604 for NDVI to the highest of 0.8505 for EVI. A new method was proposed for the estimation of GPP following the Monteith logic, which considering GPP as a product of VI × VI × PAR. Results indicated that this method can provide the best estimates of GPP as determination coefficient R² increased largely compared to the other two methods. EVI × EVI × PAR was demonstrated to be the most suitable for the estimation of GPP with the highest R² of 0.9207, which was about 10% larger as compared to GPP estimated from the single EVI. These results will be helpful for the development of new models of GPP estimation with all remote sensing inputs.     

Is it possible to discern Striga weed (Striga hermonthica) infestation levels in maize agro-ecological systems using in-situ spectroscopy?

The invasion by Striga in most cereal crop fields in Africa has posed a significant threat to food security and has caused substantial socioeconomic losses. Hyperspectral remote sensing is an effective means to discriminate plant species, providing possibilities to track such weed invasions and improve precision agriculture. However, essential baseline information using remotely sensed data is missing, specifically for the Striga weed in Africa. In this study, we investigated the spectral uniqueness of Striga compared to other co-occurring maize crops and weeds. We used the in-situ FieldSpec® Handheld 2™ analytical spectral device (ASD), hyperspectral data and their respective narrow-band indices in the visible and near infrared (VNIR) region of the electromagnetic spectrum (EMS) and four machine learning discriminant algorithms (i.e. random forest: RF, linear discriminant analysis: LDA, gradient boosting: GB and support vector machines: SVM) to discriminate among different levels of Striga (Striga hermonthica) infestations in maize fields in western Kenya. We also tested the utility of Sentinel-2 waveband configurations to map and discriminate Striga infestation in heterogenous cereal crop fields. The in-situ hyperspectral reflectance data were resampled to the spectral waveband configurations of Sentinel-2 using published spectral response functions. We sampled and detected seven Striga infestation classes based on three flowering Striga classes (low, moderate and high) against two background endmembers (soil and a mixture of maize and other co-occurring weeds). A guided regularized random forest (GRRF) algorithm was used to select the most relevant hyperspectral wavebands and vegetation indices (VIs) as well as for the resampled Sentinel-2 multispectral wavebands for Striga infestation discrimination. The performance of the four discriminant algorithms was compared using classification accuracy assessment metrics. We were able to positively discriminate Striga from the two background endmembers i.e. soil and co-occurring vegetation (maize and co-occurring weeds) based on the few GRRF selected hyperspectral vegetation indices and the GRRF selected resampled Sentinel-2 multispectral bands. RF outperformed all the other discriminant methods and produced the highest overall accuracy of 91% and 85%, using the hyperspectral and resampled Sentinel-2 multispectral wavebands, respectively, across the four different discriminant models tested in this study. The class with the highest detection accuracy across all the four discriminant algorithms, was the “exclusively maize and other co-occurring weeds” (>70%). The GRRF reduced the dimensionality of the hyperspectral data and selected only 9 most relevant wavebands out of 750 wavebands, 6 VIs out of 15 and 6 out of 10 resampled Sentinel-2 multispectral wavebands for discriminating among the Striga and co-occurring classes. Resampled Sentinel-2 multispectral wavebands 3 (green) and 4 (red) were the most crucial for Striga detection. The use of the most relevant hyperspectral features (i.e. wavebands and VIs) significantly (p ≤ 0.05) increased the overall classification accuracy and Kappa scores (±5% and ±0.2, respectively) in all the machine learning discriminant models. Our results show the potential of hyperspectral, resampled Sentinel-2 multispectral datasets and machine learning discriminant algorithms as a tool to accurately discern Striga in heterogenous maize agro-ecological systems.     

Leaf area index retrieval using gap fractions obtained from high resolution satellite data: Comparisons of approaches,scales and atmospheric effects

This study is aimed at demonstrating the feasibility of the large scale LAI inversion algorithms using red and near infrared reflectance obtained from high resolution satellite imagery. Radiances in digital counts were obtained in 10 m resolution acquired on cloud free day of August 23, 2007, by the SPOT 5 high resolution geometric (HRG) instrument on mostly temperate hardwood forest located in the Great Lakes – St. Lawrence forest in Southern Quebec. Normalized difference vegetation index (NDVI), scaled difference vegetation index (SDVI) and modified soil-adjusted vegetation index (MSAVI) were applied to calculate gap fractions. LAI was inverted from the gap fraction using the common Beer–Lambert's law of light extinction under forest canopy. The robustness of the algorithm was evaluated using the ground-based LAI measurements and by applying the methods for the independently simulated reflectance data using PROSPECT + SAIL coupled radiative transfer models. Furthermore, the high resolution LAI was compared with MODIS LAI product. The effects of atmospheric corrections and scales were investigated for all of the LAI retrieval methods. NDVI was found to be not suitable index for large scale LAI inversion due to the sensitivity to scale and atmospheric effects. SDVI was virtually scale and atmospheric correction invariant. MSAVI was also scale invariant. Considering all sensitivity analysis, MSAVI performed best followed by SDVI for robust LAI inversion from high resolution imagery.     

Analysis of Spectral Vegetation Signal Characteristics as a Function of Soil Moisture Conditions Using Hyperspectral Remote Sensing

The assessment and quantification of spatio-temporal soil characteristics and moisture patterns are important parameters in the monitoring and modeling of soil landscapes. Remote-sensing techniques can be applied to characterize and quantify soil moisture patterns, but only when dealing with bare soil. For soils with vegetation, it is only possible to quantify soil-moisture characteristics through indirect vegetation indicators, i.e. the “vitality” of plants. The “vitality” of vegetation is a sum of many indicators, whereby different stress factors can induce similar changes to the biochemical and physiological characteristics of plants. Analysis of the cause and effect of soil-moisture properties, patterns and stress factors can therefore only be carried out using an experimental approach that specifically separates the causes. The study describes an experimental approach and the results from using an imaging hyperspectral sensor AISA-EAGLE (400–970 nm) and a non-imaging spectral sensor ASD (400–2,500 nm) under controlled and comparable conditions in a laboratory to study the spectral response compared to biochemical and biophysical vegetation parameters (“vitality”) as a function of soil moisture characteristics over the entire blooming period of Ash trees. At the same time that measurements were taken from the hyperspectral sensors, the following vegetation variables were also recorded: leaf area index (LAI), chlorophyll meter value — SPAD-205, vegetation height, C/N content and leaf water content as indicators of the “vitality” and the state of the vegetation. The spectrum of each hyperspectral image was used to calculate a range of vegetation indices (VI’s) with relationships for soil moisture characteristics and stress factors. The relationship between vegetation indices and plant “vitality” indicators was analysed using a Generalized Additive Model (GAM). The results show that leaf water content is the most appropriate vegetation indicator for assessing the “vitality” of vegetation. With the Water Index (WI) it was possible to differentiate between the moisture treatments of the control, moisture drought stress and the moisture flooding treatment over the entire growing season of the plants (R ²?=?0.94). There is a correlation between the “vitality” vegetation parameters (LAI, C/N content and vegetation height) and the indicators NDVI, WI, PRI and Vog2. In our study with Ash trees the vegetation parameter chlorophyll was found not to be a suitable indicator for detecting the “vitality” of plants using the spectral indicators. There is a possibility that the sensitivity of the indicators selected was too low compared to changes in the chlorophyll content of Ash trees. Adding the co-variable ‘time’ strengthens the correlation, whereas incorporating time and moisture treatment only improves the model very slightly. This shows that changes to the biochemical and biophysical characteristics caused by phenology, overlay a differentiation of the moisture treatments.     

Combining UAV-based plant height from crop surface models,visible, and near infrared vegetation indices for biomass monitoring in barley

In this study we combined selected vegetation indices (VIs) and plant height information to estimate biomass in a summer barley experiment. The VIs were calculated from ground-based hyperspectral data and unmanned aerial vehicle (UAV)-based red green blue (RGB) imaging. In addition, the plant height information was obtained from UAV-based multi-temporal crop surface models (CSMs). The test site is a summer barley experiment comprising 18 cultivars and two nitrogen treatments located in Western Germany. We calculated five VIs from hyperspectral data. The normalised ratio index (NRI)-based index GnyLi (Gnyp et al., 2014) showed the highest correlation (R² = 0.83) with dry biomass. In addition, we calculated three visible band VIs: the green red vegetation index (GRVI), the modified GRVI (MGRVI) and the red green blue VI (RGBVI), where the MGRVI and the RGBVI are newly developed VI. We found that the visible band VIs have potential for biomass prediction prior to heading stage. A robust estimate for biomass was obtained from the plant height models (R² = 0.80–0.82). In a cross validation test, we compared plant height, selected VIs and their combination with plant height information. Combining VIs and plant height information by using multiple linear regression or multiple non-linear regression models performed better than the VIs alone. The visible band GRVI and the newly developed RGBVI are promising but need further investigation. However, the relationship between plant height and biomass produced the most robust results. In summary, the results indicate that plant height is competitive with VIs for biomass estimation in summer barley. Moreover, visible band VIs might be a useful addition to biomass estimation. The main limitation is that the visible band VIs work for early growing stages only.     

Reducing background effects in orchards through spectral vegetation index correction

Satellite remote sensing provides an alternative to time-consuming and labor intensive in situ measurements of biophysical variables in agricultural crops required for precision agriculture applications. In orchards, however, the spatial resolution causes mixtures of canopies and background (i.e. soil, grass and shadow), hampering the estimation of these biophysical variables. Furthermore, variable background mixtures obstruct meaningful comparisons between different orchard blocks, rows or within each row. Current correction methodologies use spectral differences between canopies and background, but struggle with a vegetated orchard floor. This background influence and the lack of a generic solution are addressed in this study.Firstly, the problem was demonstrated in a controlled environment for vegetation indices sensitive to chlorophyll content, water content and leaf area index. Afterwards, traditional background correction methods (i.e. soil-adjusted vegetation indices and signal unmixing) were compared to the proposed vegetation index correction. This correction was based on the mixing degree of each pixel (i.e. tree cover fraction) to rescale the vegetation indices accordingly and was applied to synthetic and WorldView-2 satellite imagery. Through the correction, the effect of background admixture for vegetation indices was reduced, and the estimation of biophysical variables was improved (ΔR² = 0.2–0.31).     

PROSAIL模型模拟下的植被指数土壤调节能力比较与适用环境分析

为了消除土壤背景信息对植被指数的影响,近几十年发展了土壤调节植被指数系列(SAVI family)。在不同环境条件下,不同指数抗土壤影响的能力不同。在总结了以消除土壤影响为目的的植被指数建立过程的基础上,利用PROSAIL辐射传输模型模拟的两组数据集,比较分析了NDVI、SAVI、TSAVI、MSAVI、OSAVI和GESAVI在不同叶面积指数(LAI)对应不同土壤背景的情况下抗土壤干扰、表达植被信息的能力,指出了不同植被指数应用的最适环境条件。结合植被指数—信噪比图,将这6种植被指数分成3类:在中低LAI值下,若植被覆盖度均匀,OSAVI和TSAVI有较强的消除土壤影响、表达植被信息的能力;当区域LAI分布不均、植被类型混杂时,MSAVI在表达植被信息时具有较好的稳定性。根据每类植被指数的特征,利用MODIS-VI和MODIS-LAI产品初步验证了上述结论的有效性。     

Leaf to canopy upscaling approach affects the estimation of canopy traits

In remote sensing applications, leaf traits are often upscaled to canopy level using sunlit leaf samples collected from the upper canopy. The implicit assumption is that the top of canopy foliage material dominates canopy reflectance and the variability in leaf traits across the canopy is very small. However, the effect of different approaches of upscaling leaf traits to canopy level on model performance and estimation accuracy remains poorly understood. This is especially important in short or sparse canopies where foliage material from the lower canopy potentially contributes to the canopy reflectance. The principal aim of this study is to examine the effect of different approaches when upscaling leaf traits to canopy level on model performance and estimation accuracy using spectral measurements (in-situ canopy hyperspectral and simulated Sentinel-2 data) in short woody vegetation. To achieve this, we measured foliar nitrogen (N), leaf mass per area (LMA), foliar chlorophyll and carbon together with leaf area index (LAI) at three vertical canopy layers (lower, middle and upper) along the plant stem in a controlled laboratory environment. We then upscaled the leaf traits to canopy level by multiplying leaf traits by LAI based on different combinations of the three canopy layers. Concurrently, in-situ canopy reflectance was measured using an ASD FieldSpec-3 Pro FR spectrometer, and the canopy traits were related to in-situ spectral measurements using partial least square regression (PLSR). The PLSR models were cross-validated based on repeated k-fold, and the normalized root mean square errors (nRMSE_cv) obtained from each upscaling approach were compared using one-way analysis of variance (ANOVA) followed by Tukey’s post hoc test. Results of the study showed that leaf-to-canopy upscaling approaches that consider the contribution of leaf traits from the exposed upper canopy layer together with the shaded middle canopy layer yield significantly (p < 0.05) lower error (nRMSE_cv < 0.2 for canopy N, LMA and carbon) as well as high explained variance (R² > 0.71) for both in-situ hyperspectral and simulated Sentinel-2 data. The widely-used upscaling approach that considers only leaf traits from the upper illuminated canopy layer yielded a relatively high error (nRMSE_cv>0.2) and lower explained variance (R² < 0.71) for canopy N, LMA and carbon. In contrast, canopy chlorophyll upscaled based on leaf samples collected from the upper canopy and total canopy LAI exhibited a more accurate relationship with spectral measurements compared with other upscaling approaches. Results of this study demonstrate that leaf to canopy upscaling approaches have a profound effect on canopy traits estimation for both in-situ hyperspectral measurements and simulated Sentinel-2 data in short woody vegetation. These findings have implications for field sampling protocols of leaf traits measurement as well as upscaling leaf traits to canopy level especially in short and less foliated vegetation where leaves from the lower canopy contribute to the canopy reflectance.     

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.