Printer Friendly

An Integrated Field and Hyperspectral Remote Sensing Method for the Estimation of Pigments Content of Stipa Purpurea in Shenzha, Tibet.

1. Introduction

Stipa purpurea, distributed in northern Tibet, is the most important and largest ecological system and the survival and development materials for herdsmen [1, 2]. At the same time, Stipa purpurea grassland known as alpine biological germplasm resources has unique species diversity and high biodiversity conservation value. Due to the severe natural conditions in northern Tibet, the ecological system of Stipa purpurea grassland is extremely fragile and the antidisturbance ability is poor [3,4]. Once the grassland vegetation is destroyed, it is difficult to recover. Photosynthetic pigment content, as an effective plant health indicator for the detection of photosynthesis and disease pollution [5,6], can be estimated and analyzed by remote sensing. This will be useful not only for analyzing the growth and health status of Stipa purpurea grassland, but also for reflecting the succession of Stipa purpurea community. It could play an important role in regulating and improving the environment of Qinghai-Tibet plateau or even the whole Eurasia.

The nutritional status of plants is closely related to pigment content and compositions in vegetation biochemical parameters. Chlorophyll is such kind of vegetation biochemical parameter and closely related to nitrogen, protein, water, and others. When the grassland is damaged, the chlorophyll content will be reduced, which can indirectly reflect the physiological condition of the grassland under external stress. Chlorophyll can absorb light and convert it into electrical energy, while carotenoid has the function of preventing chlorophyll photooxidation besides collecting and transmitting light energy [7]. The traditional quadratic survey of chlorophyll has some shortcomings such as sparse sampling points and measurement difficulties in large scale, which leads to low reliability of analysis results. The hyperspectral remote sensing provides the possibility for estimating photosynthetic pigment contents in grassland [8,9] and is important for the estimating of photosynthetic capacity, species identification, production dynamic observation, and precision fertilization [10,11].

At present, hyperspectral remote sensing has made progress in monitoring vegetation pigment contents, and some hyperspectral chlorophyll indices have been proposed by maximizing the reflection information of vegetation and minimizing the influence of external factors (e.g., aerosol, soil background, and nonphotosynthetic materials) [12-20]. For example, Horler et al. [12] studied the relationship between vegetation spectrum and chlorophyll content and analyzed the role of red edge position (wavelength of vegetation derivative spectrum near 700 nm) for the assessment of chlorophyll concentration. The contents of biochemical components (chlorophyll, lignin, nitrogen, cellulose, etc.) of 12 kinds of leaves were computed by using the first-order derivative reflectance spectroscopy and absorption characteristics method. Blackburn [15] proposed two new indices for estimating the pigment concentrations of four deciduous broad-leaved species at different wilting stages by using a simple ratio of specific pigment simple index and a normalized difference index of specific pigments. Richardson et al. [17] used traditional extraction techniques to measure chlorophyll a (Chla), chlorophyll b (Chlb), and total chlorophyll content (Chls) of 100 paper birch leaves and found that some indices based on reflectance, such as the chlorophyll NDI, were much better indicators of chlorophyll content than some of the more commonly used indices, such as [lambda]RE or NDVI. Dash and Curran [18] proposed a new index called the MERIS terrestrial chlorophyll index (MTCI) using data in three red/NIR wavebands centered at 681.25, 708.75, and 753.75 nm, and the MTCI appeared to be a most suitable index for the estimation of chlorophyll content with Medium Resolution Imaging Spectrometer (MERIS) data. The MERIS index is now a standard level-2 product of European Space Agency.

The above-mentioned indices can be classified into four types [21], single spectral index, normalized spectral index, ratio spectral index, and multiband spectral index, and can obtain very fine vegetation physiological conditions and environmental stress differences by their unique fine spectral characteristics. But, the majority of hyperspectral chlorophyll indices were developed for crops [22] or forests other than grassland [23]. Nichol and Grace [24] found that several previously published indices performed relatively poorly and yielded low coefficients, and some ecologically important species remain to be explored, despite their "applicability" across many species types [25]. Especially in alpine grassland, the hyperspectral estimation model for pigment content of Stipa purpurea community is rarely reported.

The objectives of this paper are (1) to take the sparse grass group of Stipa purpurea community as the target in alpine grassland and analyze its canopy and community spectral characteristics; (2) to extract spectral information and establish the correlation with pigments content; and (3) to analyze the spatial distribution of pigments content of Stipa purpurea in Shenzha. The remainder of this paper describes the field experiment, spectral analysis, and characteristics and the Stipa purpurea mapping results. The specific methods are as follows: (1) the hyperspectral characteristic variables were selected by spectral analysis, and the spectral data were processed by statistical method; (2) correlation analysis was carried out and the sensitive parameters were selected to establish the pigments content and spectral characteristics of the estimation model; (3) in order to reflect the distribution and succession of Stipa purpurea community, the maps are produced and analyzed by GIS tools.

2. Materials and Methods

2.1. Study Area. Shenzha County is located in the middle of Tibet, southwest of Naqu area. The average elevation of this area is 4750 m. The annual average temperature is below 0[degrees]C, and the highest temperature is 10[degrees]C every year in July-August. The minimum temperature is -40[degrees]C. The Stipa purpurea grassland is well developed in Shenzha, and it is widely distributed in slopes, hills, rivers, and lakes terraces and lake plains, which are well drained below 4900 m, forming a large grassland community. This region has vast grassland and abundant water resources, which is one of the traditional pastoral areas in Tibet. As a result of long-term overgrazing, some good grasslands have been significantly degraded, such as Stellera chamaejasme Linn, Morina kokonorica Hao, and other poisonous grass.

2.2. Spectral Data Acquisition. The selection of quadratic sampling location is very important for field spectrum acquisition process. The selected target must be representative and can truly reflect the average characteristics. In this paper, the field experiment locations are shown in Figure 1. Each location is a typical Stipa purpurea community in northern Tibet plateau of Shenzha. A total of 9 sampled plots (0.5 m x 0.5 m) were selected randomly, and 60 test units were measured. The spectral reflectance of Stipa purpurea was measured by SVC HR-1024 portable spectrometer. The spectral range of this spectrometer is from 350 nm to 2500 nm and the number of channels is 1024. In order to show the best growth status of alpine grassland, field experiments are seasoned between August and September 2015, which is the growth season of Stipa purpurea community in Shenza area. The weather is sunny, without wind or with very little wind. The measurement time is controlled between 11:00 and 14:00. In the measurement, a single vegetation type was selected. The sensor head of spectrometer is down vertically, the field angle is 4[degrees], the height of the grass canopy is about 20 cm, the coverage is nearly 100%, and the diameter is about 5.2 cm. Four sampling points were selected randomly in each plot for spectral determination, and the mean value was taken as the repeat spectral reflectance value. In order to reduce the impact of atmospheric changes, plots and reference whiteboard was measured alternately (as shown in Figure 2). At each sample point, GPS, grass height, and coverage were recorded at the same time.

2.3. Remote Sensing Data and Its Preprocessing. The HJ-1A hyperspectral imager (HSI) product was downloaded from July to August 2015 from China Centre for Resources Satellite Data and Application, with ground width of 50 km, spatial resolution of 100 m, and 128 bands ranging from 450 nm to 950 nm. The hyperspectral data was selected according to the growing season of Stipa purpurea in study area. The principle of data selection is to remove the image covering the same area with a larger sampling time or poor quality, thus reducing the error caused by mismatch between the ground sampling and satellite monitoring time. Finally, 46 images were selected; 13 images contain quadratic samples. And then these images were postprocessed such as band fusion, geometric correction, and projection conversion.

2.4. Determination of Photosynthetic Pigments. In order to extract photosynthetic pigments of alpine grassland, the optical density was measured by spectrophotometry at the wavelength of 663 nm, 647 nm, 537 nm, and 470 nm. According to Lambert-Beer law, The Chla, Chlb, Chls, and Carotenoids (Cars) were calculated. Each process was repeated three times, and the average value was taken as the photosynthetic pigment content. A total of 70 samples were determined, and four photosynthetic pigment contents were determined for each sample.

[mathematical expression not reproducible]. (1)

2.5. Feature Parameter Extraction of Hyperspectral Data. In order to explore the characteristics of canopy reflectance spectrum of Stipa purpurea, some useful narrow bands were selected from hyperspectral bands for remote sensing information extraction. The aboveground biomass was modeled by using hyperspectral data and various transformations, such as red, yellow, and blue optical parameters, vegetation index, green reflection peak, and red light absorption valley. The spectral data of the experimental units were smoothed and deaveraged four times. The parameters such as the position, area, vegetation index, and derivative of the spectra were calculated (as shown in Table 1). The linear and nonlinear regression equations were established by comparing these spectral variables with the light and pigment contents, and the optimal estimation model of one or more variables was preliminarily screened.

2.6. Regression Method. A stepwise regression analysis method was used to determine the biomass-related spectral variable equation, and the measured biomass and spectral variable value was then input. The output result was a series of multivariate linear equations containing different spectral variables, as well as the corresponding determination coefficient ([R.sup.2]), and test value (F). The stepwise regression model was as follows:

Y = [a.sub.0] + [a.sub.1][x.sub.1] + [a.sub.2][x.sub.2] + ... + [a.sub.i][x.sub.i] + ..., (2)

where [x.sub.i] is the spectral variable value; Y is the biomass estimation value; [a.sub.0] is a constant item; and [a.sub.i] represents the partial regression coefficient.

3. Results and Analysis

3.1. Hyperspectral Analysis of Stipa Purpurea. Grassland vegetation, atmosphere, water, soil, and other factors form a comprehensive reflection spectrum of grassland. The reflectance is affected by many factors, such as grassland population, type, vegetation water content, vegetation coverage, air condition, and soil condition. Therefore, the spectrum analysis of Stipa purpurea is the basis of remote sensing monitoring of forage resources in Tibetan plateau, and its spectral characteristics are the integrated responses of Stipa purpurea and its habitat conditions.

It can be seen from Figure 3 that the overall trend of 5 different coverages of Stipa purpurea samples is consistent, although the spectral reflectance is varying. The blue-band with the largest water penetration is located in the chlorophyll absorption zone of green plants between 490 and 520 nm. It shows that blue-band is sensitive to the chlorophyll concentration and is helpful for the identification of plant stress. This conclusion is consistent with Jin et al.'s research [26]. The green reflectance of healthy grass between 490 and 520 nm is sensitive to grass reflectance, and it could be used to identify grassland type and plant productivity state. At the same time, the green band reflection peaks also increase with the lift of coverage and growth conditions. The sharp rising portion between the red and the near-infrared bands of high reflectance is located in the vicinity of the red edge of 700-750 nm. In the decreasing process of coverage, the red edge of the spectrum will shift to the long wave direction, and the offset is about 10~40 nm. The near-infrared band is located in the meadow high reflection area (760-900 nm). The spectral characteristics are controlled by the cell structure of the grassland and are sensitive to the change of the grassland type, density, and growth force. The higher the vegetation height, the coverage, and the productivity, the higher the reflectance is. The canopy reflectance is usually between 20% and 40%, which is due to the multiple reflection and scattering of the internal tissue structure of plant leaves and is mainly determined by biomass and leaf area index.

3.2. Data Analysis. The hyperspectral data after transformation were correlated with pigments content and the results are shown in Table 2. There are significant correlations between chlorophyll content and spectral characteristic variables; in particular, the correlation coefficient of chlorophyll content with green peak reflectance, red valley reflectance, green peak area, RVI, NDVI, EVI, RDVI, VI1, VI2, and VI3 reached 0.01 significant level. The correlation coefficient of Chlb of Stipa purpurea and RVI is the largest value of 0.728. The correlation coefficient of Chla, Chlb, and Chls to RVI, EVI, RDVI, and VI3 is 0.6. The correlation coefficient between carotenoid content and 16 characteristic variables is generally low, which indicates that the spectral characteristic variables of Stipa purpurea are correlated with chlorophyll, and the chlorophyll content could be estimated by using the characteristic variables with good correlation.

3.3. Model Construction and Verification. In many studies, hyperspectral remote sensing data is used to estimate the chlorophyll content. The main method was to find the bands relation to chlorophyll content by regression method. In high spectral narrow band data, the green band uses 553 nm reflectivity, the red band uses 670 nm reflectivity, the near-infrared band uses 801 nm reflectivity to carry on the inversion parameter operation, and the result is valid. This is due to the fact that 553 nm and 670 nm are used to determine the most responsive narrow band position according to [R.sub.g], [R.sub.r] two parameters; 801 nm is used to screen out the narrow band at near-infrared band according to many experts in the grassland vegetation index parameters [26]. The regression models of chlorophyll and band combination are shown in Table 3.

Based on the analysis of typical hyperspectral variables sensitive to chlorophyll content of Stipa purpurea, 10 variables with significant correlation with chlorophyll were extracted in Table 3, and SPSS22.0 was used to analyze Chla, Chlb, and Chls. We analyzed the linear, quadratic, logarithmic, and exponential models and selected the higher coefficient of determination of Chla and Chlb and total chlorophyll as the basis for the modeling of Stipa purpurea (as shown in Figure 4). The observed data were randomly selected and divided into two groups: one for the training samples (n = 30) is used to establish the regression relationship.

As shown in Table 4, all the regression equations [R.sup.2] are tested by 0.01 significance and varied between 0.623 and 0.708 with an average value of 0.658. For the Chla content estimation, the best fitting model is EVI variable parabola equation, [R.sup.2] is 0.708, F-test value is 31.9, and the estimated standard error is small; for Chlb, the best fitting model is the parabolic equation of RVI, NDVI, and VI2 variables, the F-test value is large, and the standard error is small. For the estimation of Chls, the best fitting model is also the parabolic equation of RVI and NDVI, [R.sup.2] value is 0.644, and F test value reached 12.79.

3.4. Spatial Pattern Analysis. The best models for Chla, Chlb, and Chls estimation were obtained and validated by correlation coefficient model (as shown in Table 3, Figure 4). The hyperspectral raster data can be directly used to calculate the spatial results. In the calculation process, the green band uses 553 nm reflectivity which is 38 band of HJ-1A image, 670 nm reflectivity as red band which is 70 band of HJ-1A image, and 801 nm reflectivity as the near-infrared band. Finally, EVI parameters of Chla, RVI value of Chlb, and EVI parameters of Chls were calculated and substituted into the estimation model, respectively. The photosynthetic pigment contents in Shenzha area were estimated by HJ-1A images.

The contents of Chla, Chlb, and Chls in Shenzha area are shown in Figure 5. The Chla is in the range of 0-30 mg/g. The Chla values in bare soil and alpine steppe region are low, the Siling Co. and Kering Co. are up to 30 mg/g, and the alpine swamp meadow near lake is 18 mg/g. The range of Chlb is 0-15 mg/g. The Chlb values in the alpine swamp meadow and alpine meadow near Kering Co. are 6-10 mg/g and in the alpine grassland are varied from 10 to 15 mg/g. The Chls is between 7 and 50 mg/g, and its spatial distribution is similar to Chlb. The above results indicate that the distribution of Chla is opposite to Chlb and total chlorophyll. Chlb determines the distribution pattern of Chls in the water area, but the Chlb value is lower than the Chla. The Chls in the wetland is between 15 and 25 mg/g, while the chlorophyll in the alpine steppe and alpine desert steppe is lower than that in the dry soil condition.

3.5. Photosynthetic Pigment Analysis. The Chla and Chlb are the most studied parameters in photosynthetic pigments content. The lutein, carotenoids, and other photosynthetic pigments are not only the composition of light-emitting pigment. Under the external stress they are more involved in the leaves of yellow Cyclin and other metabolic pathways. The correlation coefficients of the pigment content of Stipa grandis were correlated with each other by 0.01. In particular, the correlation coefficient between the Chla and Chls is 0.98 (Table 5), indicating that the trend of change between different pigments is consistent, and the response of different pigments in the spectral characteristic parameters are similarity (Table 3), such as Chla, Chla, and Chls of RVI.

In single factor model, [R.sup.2] of RVI is higher than 0.65, and [R.sup.g] and [SD.sub.g] are both lower than 0.42. It is shown that the characteristic parameters are consistent with the simulation of the photosynthetic pigment content. Some parameters are related to all the pigment correlations, and some are independent. However, in the stepwise regression model, not all parameters with higher correlation with the pigment can iterate into the equation but rather enter the equation in the form of complementary effects. Several similar parameters are excluded from the model operation (Table 4). The spectral parameters of the hyperspectral remote sensing image are calculated by comparing the band of the characteristic parameters involved in the regression model with the band of the hyperspectral image of HJ-1A and then are substituted into the regression model. The results showed the spatial distribution of Chla, Chlb, and Chls (Figure 5).

The chlorophyll content of Shenzha County was calculated by the spatial distribution data of photosynthetic pigments (Table 6). The mean estimated value content of Chla in Stipa bungeana is about 4.88 times that of Chlb, while Cars is generally slightly more than Chlb, so chlorophyll a plays a dominant role in absorption. The chlorophyll content of the measured sample was slightly higher than that of the spatial estimation, because the change of photosynthetic pigment in the community of Stipa purpurea occurred in the alpine steppe region with low soil water content, and the change of the pigment in the lake and water was small. Some of the images of the snow-covered area (seasonal snow) affect the actual monitoring of photosynthetic pigments, so the estimate is slightly lower, but the estimated relative error is small, below 7.24%.

4. Conclusion

Based on the correlation between ground measured spectral characteristics and photosynthetic pigments, [R.sub.g], [R.sub.r], [SD.sub.g], RVI, NDVI, EVI, RDVI, VI1, VI2, and VI3 characteristic parameters with the correlation level of 0.01 and ground photosynthetic pigments are used to establish single factor equations and select the appropriate spectral characteristic parameters. In these parameters, [R.sup.2] of [R.sub.g], [R.sub.r], and [SD.sub.g] are less than 0.42, so these three parameters do not participate in the next step budget. The three spectral inversion models of Chla, Chlb, and Chls were established by using RVI, NDVI, EVI, RDVI, VI1, VI2, and VI3. Finally, the spatial distribution pattern of Chla, Chlb, and Chls was obtained by using the scale push method and HJ-1A remote sensing image, and the characteristic parameter band in the model was mapped to the band in the hyperspectral remote sensing image. The relative error is small, and the monitoring and investigation of the grassland growth and health status of large-scale area in high altitude and alpine region are realized by using the spatial estimation value and the ground measured value.

(1) In the fitting analysis of pigments content and hyperspectral parameters, 670 nm is used as the red band, 801 nm as the near-infrared band, and 553 nm as the green band to construct pigments content inversion model. Among the 16 spectral parameters, 10 parameters are correlated with chlorophyll, and the correlation order is Chlb > Chla > Chls. The best correlation index between Chla, Chlb, and Chls is RVI, and the correlation coefficients are 0.707, 0.728, and 0.694, respectively. The results indicate that the estimation of the chlorophyll content of Stipa purpurea is more feasible by using spectral characteristic variables.

(2) The chlorophyll estimation model of Stipa purpurea was established by using 10 extremely significant correlation spectral parameters. The estimation model of Chla is based on vegetation index EVI, estimation of Chlb is based on RVI, NDVI, and VI2, Chls model is based on RVI and NDVI, and the coefficient of determination ([R.sup.2]) of fitted equation is 0.6. The relative error of the model estimation is small, which is below 7.24%, indicating that the model is feasible in the region.

(3) The distribution of Chla is opposite to Chlb and Chls in the water area. Chlb in water area determines the distribution pattern of Chls, but Chlb value is lower than Chla; thus, the Chls in wetlands near lakes is 15-25 mg/g. The Chls of alpine steppe and alpine desert steppe is mainly affected by Chla in dry soil condition.

(4) The proposed method in this paper can effectively estimate the pigments content of Stipa purpurea in Shenzha, Tibet. The results could reflect the distribution and succession pattern of Stipa purpurea community. However, there are some errors in estimation model for pigment content based on hyperspectral vegetation physiological parameters. This is due to the fact that the vegetation type, growth stage, growth condition, and the signal-to-noise ratio of the spectrometer itself will influence the establishment of statistical model.

Conflicts of Interest

The authors declare that they have no conflicts of interest.


This research is supported by the National Natural Science Funds of China (41301094); the Hundred Young Talents Program of the Institute of Mountain Hazards and Environment (SDSQB-2015-02); the Lead Strategic Project of the Chinese Academy of Sciences (XDB03030507); the 973 National Projects (2015CB452706); the Open Fund for Spatial Information Technology Key Lab Bases of the Ministry of Land and Resources of China (KLGSIT2014-06).


[1] X. M. Zhou, Q. J. Wang, Y. Q. Zhang, X. Q. Zhao, and Y. P. Lin, "Quantitative analysis of vegetation succession in alpine meadow under different grazing intensity," Journal of Plant Ecology and Botany, vol. 11, no. 4, pp. 276-285, 1987.

[2] M. J. Duan, Q. Z. Gao, Y. F. Wan et al., "Effect of grazing on community characteristics and species diversity of stipa purpurea alpine grassland in northern tibet," Acta Ecologica Sinica, vol. 30, no. 14, pp. 3892-3900, 2010.

[3] D. Lu, Y. Zhao, R. Han, L. Wang, and P. Qin, "The complete chloroplast genome sequence of the Purple Feathergrass Stipa purpurea (Poales: Poaceae)," Conservation Genetics Resources, vol. 8, no. 2, pp. 101-104, 2016.

[4] P. P. Yue, X. F. Lu, R. R. Ye et al., "Distribution of Stipa purpurea steppe in the northeastern qinghai-xizang plateau (China)," Russian Journal of Ecology, vol. 42, no. 1, pp. 50-56, 2011.

[5] D. Filella and J. Penuelas, "The red edge position and shape as indicators of plant chlorophyll content, biomass and hydric status," International Journal of Remote Sensing, vol. 15, no. 7, pp. 1459-1470, 1994.

[6] A. C. Madeira, A. Mendonca, and M. E. Ferreira, "Relationship between spectroradiometric and chlorophyll measurements in green beans," Communications in Soil Science and Plant Analysis, vol. 31, no. 5-6, pp. 631-643, 2000.

[7] A. A. Gitelson, G. P. Keydan, and M. N. Merzlyak, "Three-band model for noninvasive estimation of chlorophyll, carotenoids, and anthocyanin contents in higher plant leaves," Geophysical Research Letters, vol. 33, no. 11, Article ID L11402, 2006.

[8] M. J. Gallardo, J. P. Staforelli, P. Meza, I. Bordeu, and S. Torres, "Characterization of Chromobacterium violaceum pigment through a hyperspectral imaging system," AMB Express, vol. 4, no. 4, pp. 2-9, 2014.

[9] M. Gupta, P. K. Srivastava, S. Mukherjee, and G. Sandhya Kiran, "Chlorophyll retrieval using ground based hyperspectral data from a tropical area of india using regression algorithms," in Remote Sensing Applications in Environmental Research, vol. 4, pp. 177-194, 2014.

[10] F. Flores-de-Santiago, J. M. Kovacs, and F. Flores-Verdugo, "The influence of seasonality in estimating mangrove leaf chlorophyll-a content from hyperspectral data," Wetlands Ecology and Management, vol. 21, no. 3, pp. 193-207, 2013.

[11] D. Haboudane, J. R. Miller, N. Tremblay, P. J. Zarco-Tejada, and L. Dextraze, "Integrated narrow-band vegetation indices for prediction of crop chlorophyll content for application to precision agriculture," Remote Sensing of Environment, vol. 81, no. 2-3, pp. 416-426, 2002.

[12] D. N. H. Horler, J. Barber, and A. R. Barringer, "Effects of heavy metals on the absorbance and reflectance spectra of plants," International Journal of Remote Sensing, vol. 1, no. 2, pp. 121-136, 1980.

[13] A. A. Gitelson, C. Buschmann, and H. K. Lichtenthaler, "The chlorophyll fluorescence ratio F735F700 as an accurate measure of the chlorophyll content in plants," Remote Sensing of Environment, vol. 69, no. 3, pp. 296-302, 1999.

[14] D. A. Sims and J. A. Gamon, "Relationships between leaf pigment content and spectral reflectance across a wide range of species, leaf structures and developmental stages," Remote Sensing of Environment, vol. 81, no. 2-3, pp. 337-354, 2002.

[15] G. A. Blackburn, "Quantifying chlorophylls and carotenoids at leaf and canopy scales: An evaluation of some hyperspectral approaches," Remote Sensing of Environment, vol. 66, no. 3, pp. 273-285, 1998.

[16] J. G. P. W. Clevers and L. Kooistra, "Using hyperspectral remote sensing data for retrieving total canopy chlorophyll and nitrogen content," in Proceedings of 3rd Workshop on Hyperspectral Image and Signal Processing: Evolution in Remote Sensing, WHISPERS 2011, June 2011.

[17] A. D. Richardson, S. P. Duigan, and G. P. Berlyn, "An evaluation of noninvasive methods to estimate foliar chlorophyll content," New Phytologist, vol. 153, no. 1, pp. 185-194, 2002.

[18] J. Dash and P. J. Curran, "Evaluation of the MERIS terrestrial chlorophyll index (MTCI)," Advances in Space Research, vol. 39, no. 1, pp. 100-104, 2007.

[19] J. Delegido, L. Alonso, G. Gonzalez, and J. Moreno, "Estimating chlorophyll content of crops from hyperspectral data using a normalized area over reflectance curve (NAOC)," International Journal of Applied Earth Observation and Geoinformation, vol. 12, no. 3, pp. 165-174, 2010.

[20] P. Chen, D. Haboudane, N. Tremblay, J. Wang, P. Vigneault, and B. Li, "New spectral indicator assessing the efficiency of crop nitrogen treatment in corn and wheat," Remote Sensing of Environment, vol. 114, no. 9, pp. 1987-1997, 2010.

[21] W. M. Yong and X. M. Wang, "Progress on grassland chlorophyll content estimation by hyperspectral analysis," Progress in Geography, vol. 35, no. 1, pp. 25-34, 2016.

[22] A. Bannari, K. S. Khurshid, K. Staenz, and J. W. Schwarz, "A comparison of hyperspectral chlorophyll indices for wheat crop chlorophyll content estimation using laboratory reflectance measurements," IEEE Transactions on Geoscience and Remote Sensing, vol. 45, no. 10, pp. 3063-3074, 2007.

[23] G. Le Maire, C. Francois, and E. Dufrene, "Towards universal broad leaf chlorophyll indices using PROSPECT simulated database and hyperspectral reflectance measurements," Remote Sensing of Environment, vol. 89, no. 1, pp. 1-28, 2004.

[24] C. J. Nichol and J. Grace, "Determination of leaf pigment content in Calluna vulgaris shoots from spectral reflectance," International Journal of Remote Sensing, vol. 31, no. 20, pp. 5409-5422, 2010.

[25] G. L. Jin, J. Z. Zhu, H. L. Liu, S. M. Tang, and C. H. Wang, "Study on physiology/ecology adaptation of main plant in degraded seriphidium transiliense desert rangeland," Acta Agrectir Sinica, vol. 19, no. 1, pp. 26-30, 2011.

[26] H. A. Jin, D. W. Liu, K. S. Song, Z. M. Wang, F. Li, and H. J. Liu, "Comparing the performance of broad-band and narrowband vegetation indices for estimation of soybean LAI," System Sciences and Comprehensive Studies in Agriculture, vol. 23, no. 4, pp. 503-508, 2007.

Bo Kong, (1) Bing He, (1) Huan Yu, (2) and Yu Liu (2)

(1) Institute of Mountain Hazards and Environment, Chinese Academy of Sciences, Chengdu 610041, China

(2) College of Earth Sciences, Chengdu University of Technology, Chengdu 610059, China

Correspondence should be addressed to Bing He;

Received 28 November 2016; Revised 20 March 2017; Accepted 7 May 2017; Published 22 June 2017

Academic Editor: Huade Guan

Caption: Figure 1: The study area and schematic plot of sample point.

Caption: Figure 2: Field experiment.

Caption: Figure 3: Comparison of different coverage of Stipa purpurea spectral reflectance curve.

Caption: Figure 4: Results of the optimal estimation model of the photosynthetic pigment of Stipa purpurea.

Caption: Figure 5: The map of photosynthetic pigment content.
Table 1: High spectral characteristic parameter selection.

Types               Name               Abbreviation

                    Blue margin        [D.sub.b]

                    Blue margin        [SD.sub.b]

Spectral area and   Yellow margin      [D.sub.y]
location            value

                    Yellow margin      [SD.sub.y]

                    Red margin         [D.sub.r]

                    Red margin         [SD.sub.r]

                    Green peak         [R.sub.g]

                    Green peak         [SD.sub.g]

                    Red valley         [R.sub.r]

                    Ratio between      VI1
                    green peak
                    reflectance and
                    red valley

Ratio parameters    Normalization      VI2
of spectral area    value of green
and location        peak reflectance
                    and red valley

                    Ratio between      VI3
                    red margin area
                    and blue margin

                    Ratio vegetation   RVI

Vegetation          Normalized         NDVI
indices             difference
                    vegetation index

                    Enhanced           EVI
                    vegetation index

                    Again normalized   RDVI
                    vegetation index

Types              Calculation formula           Literature

                   Maximum of first              Elvidge et al., 1995
                   derivative spectrum in
                   490 to 520 nm

                   Integration of first          Elvidge et al., 1995
                   derivative in
                   490 to 520 nm

Spectral area and  Maximum of first              Elvidge et al., 1995
location           derivative spectrum in
parameters         560 to 640 nm

                   Integration of first          Elvidge et al., 1995
                   derivative in
                   560 to 640 nm

                   Maximum of first              Elvidge et al., 1995
                   derivative spectrum in
                   670 to 760 nm

                   Integration of first          Pu et al., 1993
                   derivative in
                   670 to 760 nm

                   Maximum of band               Elvidge et al., 1995
                   reflectance in
                   520 to 560 nm

                   Area of original              Elvidge et al., 1995
                   spectrum in 520 to 560
                   nm wavelength

                   Minimum of band               Elvidge et al., 1995
                   reflectance in 650 to
                   690 nm wavelength

                   [[rho].sub.g]/                Gitelson et al., 1996

Ratio parameters   [[rho].sub.g] -               Gitelson et al., 1996
of spectral area   [[rho].sub.r]/
and location       [[rho].sub.g] +

                   [SD.sub.r]/                   Gitelson et al., 2003

                   [[rho].sub.NIR]/              Jordan, 1969

Vegetation         [[rho].sub.NIR] -             Rouse, 1973
indices            [[rho].sub.Red]/
                   [[rho].sub.NIR] +

                   [[rho].sub.NIR] -             Liu, 1995
                   [[rho].sub.NIR] +
                   [C.sub.1][[rho].sub.Red] -
                   [C.sub.2][[rho].sub.B] +
                   L (1 + L)

                   [square root of               Richardson, 1992
                   (NDVI x ([[rho].sub.NIR] -

Table 2: The correlation coefficient of Stipa pigment content and
spectral variables.

                Characteristic    Chla      Chlb      Chls      Cars

                  [D.sub.b]      -.409*    -.407*    -.377*    -0.007
Spectral          [D.sub.y]      -0.016    -0.007    -0.038    -0.157
location          [D.dub.r]      -.424*    -.398*    -.417*    0.015
variable          [R.sub.g]      -.579**   -.578**   -.556**   -0.053
                  [R.sub.r]      -.589**   -.592**   -.571**   -0.035

                  [SD.sub.b]     -0.414*   -.404*    -.395*    -0.098
Spectral area     [SD.sub.y]     -0.315    -0.322     -0.32     0.038
variable          [SD.sub.r]      0.123     0.124     0.094    -0.09
                  [SD.sub.g]     -.583**   -.584**   -.561**   -0.06

                     RVI         .707**    .728**    .694**     0.14
Vegetation           NDVI        .594**    .607**    .571**     0.026
variable             EVI         .606**    .623**    .584**     0.075
                     RDVI        .661**    .679**    .644**     0.09

                     VI1         .597**    .627**    .609**     0.131
VI variable          VI2         .569**    .594**    .578**     0.104
                     VI3         .626**    .630**    .585**     0.045

Note. * and ** indicate significant levels at 0.05 and 0.01,

Table 3: Estimation model and evaluation index of the photosynthetic
pigment content of Stipa purpurea.

Variable     Model   Model function

             Chla    y = 18.941[e.sup.-0.101x]
[R.sub.g]    Chlb    y = 4.9395[e.sup.-0.134x]
             Chls    y = 0.0483[x.sup.2] - 2.031x + 24.361

             Chla    y = 0.0498[x.sup.2] - 1.8608x + 20.157
[R.sub.r]    Chlb    y = 0.0103[x.sup.2] - 0.4021x + 4.3275
             Chls    y = 0.0685[x.sup.2] - 2.3681x + 23.633

             Chla    y = 3E - 05[x.sup.2] - 0.0524x + 22.071
[SD.sub.g]   Chlb    y = 19.434[e.sup.-0.003x]
             Chls    y = 4E - 05[x.sup.2] - 0.0554z + 24.53

             Chla    y = 0.0087[x.sup.2] + 2.4105x - 1.2277
RVI          Chlb    y = 0.3113[x.sup.2] - 1.4038x + 2.2996
             Chls    y = 0.1029[x.sup.2] - 1.1332x + 11.11

             Chla    y = 80.457[x.sup.2] - 55.9x + 13.24
NDVI         Chlb    y = 29.211[x.sup.2] - 21.795x + 4.5795
             Chls    y = 140.84[x.sup.2] - 110.27x + 24.82

             Chla    y = 9.7705[x.sup.2] - 13.358x + 8.1238
EVI          Chlb    y = 3.5428[x.sup.2] - 5.4181x + 2.6148
             Chls    y = 20.154[x.sup.2] - 44.399x + 32.094

             Chla    y = 3.2447[x.sup.2] - 1.5165x + 2.8664
RDVI         Chlb    y = 3.114[x.sup.2] - 5.8428x + 3.3762
             Chls    y = 6E - 06[x.sup.2] - 0.0098x + 11.83

             Chla    y = 7.6192[x.sup.2] - 2.3606x - 1.567
VI1          Chlb    y = 13.379[x.sup.2] - 27.295x + 14.61
             Chls    y = 47.726[x.sup.2] - 96.932x + 53.56

             Chla    y = 127.15[x.sup.2] + 17.461x + 3.7539
VI2          Chlb    y = 127.15[x.sup.2] + 17.461x + 3.75
             Chls    y = 337.33[x.sup.2] - 10.643x + 4.124

             Chla    y = 436.67[x.sup.2] - 617.11x + 220.8
VI3          Chlb    y = 0.0378[x.sup.2] - 0.3662x + 1.507
             Chls    y = 0.2864[x.sup.2] - 3.8448x + 17.188

Variable     Model  Training sample          Validation
                     (n = 30)               sample (n = 20)
                     [R.sup.2]     F        [R.sup.2]   SEE

             Chla     0.360      15.755     0.337      0.318
[R.sub.g]    Chlb     0.421      20.358     0.400      0.734
             Chls     0.402      9.086      0.358      3.226

             Chla     0.522      14.737     0.486      3.415
[R.sub.r]    Chlb     0.430      21.154     0.410      0.950
             Chls     0.484      12.686     0.446      2.996

             Chla     0.384      8.409      0.338      3.877
[SD.sub.g]   Chlb     0.425      20.706     0.405      0.731
             Chls     0.401      9.053      0.357      3.229

             Chla     0.655      25.575     0.629      2.903
RVI          Chlb     0.712      33.423     0.691      0.687
             Chls     0.701      29.367     0.662      1.713

             Chla     0.686      29.554     0.663      2.766
NDVI         Chlb     0.653      25.415     0.627      0.755
             Chls     0.612      21.281     0.583      2.600

             Chla     0.710      32.995     0.688      0.864
EVI          Chlb     0.702      31.812     0.680      0.700
             Chls     0.556      16.913     0.523      2.034

             Chla     0.659      26.097     0.634      2.884
RDVI         Chlb     0.699      31.412     0.677      0.703
             Chls     0.655      25.682     0.630      1.792

             Chla     0.579      18.593     0.525      3.285
VI1          Chlb     0.583      18.839     0.552      0.828
             Chls     0.436      10.454     0.395      3.133

             Chla     0.589      19.320     0.558      3.167
VI2          Chlb     0.589      19.349     0.559      0.822
             Chls     0.444      10.784     0.403      3.112

             Chla     0.668      27.196     0.644      2.877
VI3          Chlb     0.552      16.652     0.519      0.858
             Chls     0.549      16.415     0.515      2.804

Table 4: Stepwise regression models of photosynthetic pigment in
Stipa purpurea with enhanced spectral variables (with VIF = 35,
[R.sup.2] is the coefficient of determination, F is the value of F
distribution, and P is the significant level).

Soil properties                Regression models

Chla                          y = 0.466 - 9.17EVI
Chlb               y = 12.69 + 11.22RVI - 6.57NDVI + 0.027VI2
Chls                    y = 0.902 + 0.741RVI + 1.866NDVI

Soil properties    [R.sup.2]   F       P

Chla                0.708     31.9    0.000
Chlb                0.623     48.36   0.000
Chls                0.644     12.79   0.000

Table 5: The correlation coefficient between different pigments.

pigment     Chla     Chlb     Chls     Cars

Chla         1      0.8597   0.9886   0.7540
Chlb       0.8597     1      0.9266   0.5126
Chls       0.9886   0.9266     1      0.7060
Cars       0.7540   0.5126   0.7060     1

Table 6: The chlorophyll content of Shenzha County (mg/g).

Pigment   Maximum      Maximum     Minimum     Minimum     Average
          measured    estimated    measured    estimated   measured
            value       value        value       value       value

Chla        31.5         29.5        1.22        0.65        10.4
Chlb        17.2         14.9        0.40        0.01        2.50
Chls        54.6         48.5        0.96        0.29        12.8

Pigment     Mean      Estimated   Relative
          estimated   standard    error
            value     deviation     (%)

Chla       9.12        7.86       1.57
Chlb       1.87        2.07       7.24
Chls       9.23        10.7       6.95
COPYRIGHT 2017 Hindawi Limited
No portion of this article can be reproduced without the express written permission from the copyright holder.
Copyright 2017 Gale, Cengage Learning. All rights reserved.

Article Details
Printer friendly Cite/link Email Feedback
Title Annotation:Research Article
Author:Kong, Bo; He, Bing; Yu, Huan; Liu, Yu
Publication:Mathematical Problems in Engineering
Article Type:Report
Geographic Code:9CHIN
Date:Jan 1, 2017
Previous Article:Subspace Clustering with Sparsity and Grouping Effect.
Next Article:A General Framework for Splitting Concepts for Cocycles over Generalized Nonautonomous Dynamical Systems.

Terms of use | Privacy policy | Copyright © 2020 Farlex, Inc. | Feedback | For webmasters