Next Article in Journal
A Reliable Data Transmission Model for IEEE 802.15.4e Enabled Wireless Sensor Network under WiFi Interference
Next Article in Special Issue
A Weekly Indicator of Surface Moisture Status from Satellite Data for Operational Monitoring of Crop Conditions
Previous Article in Journal
Faller Classification in Older Adults Using Wearable Sensors Based on Turn and Straight-Walking Accelerometer-Based Features
Previous Article in Special Issue
Using Color, Texture and Object-Based Image Analysis of Multi-Temporal Camera Data to Monitor Soil Aggregate Breakdown
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Adaptive Grouping Distributed Compressive Sensing Reconstruction of Plant Hyperspectral Data

College of Life Information Science & Instrument Engineering, Hangzhou Dianzi University, Hangzhou 310018, China
*
Authors to whom correspondence should be addressed.
Sensors 2017, 17(6), 1322; https://doi.org/10.3390/s17061322
Submission received: 30 March 2017 / Revised: 26 May 2017 / Accepted: 1 June 2017 / Published: 7 June 2017
(This article belongs to the Special Issue Precision Agriculture and Remote Sensing Data Fusion)

Abstract

:
With the development of hyperspectral technology, to establish an effective spectral data compressive reconstruction method that can improve data storage, transmission, and maintaining spectral information is critical for quantitative remote sensing research and application in vegetation. The spectral adaptive grouping distributed compressive sensing (AGDCS) algorithm is proposed, which enables a distributed compressed sensing reconstruction of plant hyperspectral data. The spectral characteristics of hyperspectral data are analyzed and the joint sparse model is constructed. The spectral bands are adaptively grouped and the hyperspectral data are compressed and reconstructed on the basis of grouping. The experimental results showed that, compared with orthogonal matching pursuit (OMP) and gradient projection for sparse reconstruction (GPSR), AGDCS can significantly improve the visual effect of image reconstruction in the spatial domain. The peak signal-to-noise ratio (PSNR) at a low sampling rate (the sampling rate is lower than 0.2) increases by 13.72 dB than OMP and 1.66 dB than GPSR. In the spectral domain, the average normalized root mean square error, the mean absolute percentage error, and the mean absolute error of AGDCS is 35.38%, 31.83%, and 33.33% lower than GPSR, respectively. Additionally, AGDCS can achieve relatively high reconstructed efficiency.

1. Introduction

Hyperspectral technology is a breakthrough technology in agriculture remote sensing which enables the dynamic and precise monitoring of crop types and crop growth. Hyperspectral remote sensing technology has been widely used in estimating the yield of crops, agricultural resources surveying, agricultural disaster monitoring, and precision agriculture [1]. Plant quantitative remote sensing technology is widely used in a variety of applications by mining spectral information and setting up spectral retrieving model. For the estimation of crop yield, Nuarrsa et al. extracted a rice area with an overall accuracy of 87.91% using the normalized difference vegetation index (NDVI), radar vegetation index (RVI), and soil-adjusted vegetation index (SAVI) from MODIS time series data [2]. Tornosa et al. assessed the potential of different spectral indices for monitoring rice agricultural practices and hydroperiod dynamics by combining phenometrics and statistical time series approaches [3]. Kang et al. developed spectral indices that can reduce the effects of varied canopy structure and growth stages for the estimation of leaf Chl [4]. Fei et al. optimized the band combinations further, and identified the optimized central bands and suitable bandwidths of the three-band nitrogen planar domain index (NPDI) for estimating the aerial N uptake, N concentration, and above-ground biomass [5]. Atherton et al. linked spectral measurements of fluorescence and the PRI to photosynthesis dynamics at the leaf scale and over short time-scales [6]. Heli et al. found that the values of crop variables may not be accurately determined when they are based solely on the measurements of leaves, especially only the upper leaves, as the values varied greatly among the different vertical leaf and stem layers and the different modules, such as leaves, stems, and spikes [7]. Qiang et al. proposed an iterative method which integrate MODIS, VEGETATION, and MISR data to improve the estimation of leaf area index (LAI) climatology [8]. Mahlein et al. developed specific spectral disease indices (SDIs) for the detection of diseases in crops [9]. Alicia et al. assessed AS1 and AS2 behavior over a cotton crop growing period, testing whether function-fitting procedures can be used to model MODIS AS1 and AS2 and NDVI time series and derive objective AS1 and AS2 phenological metrics that can be used to monitor cotton phenological stages [10]. Veraverbeke et al. evaluated the discriminatory power of existing VIs and thermally-enhanced indices in burned land applications [11]. Jochem et al. introduced an automated spectral band analysis tool (BAT) based on Gaussian process regression (GPR) for the spectral analysis of vegetation properties [12]. Abderrazak et al. proposed a spatiotemporal monitoring method of soil salinization in the Tadla plain in Central Morocco using spectral indices derived from Thematic Mapper (TM) and Operational Land Imager (OLI) data [13]. Ferner et al. tested whether spatio-temporal information on the quality (metabolizable energy content, ME) and quantity (green biomass, BM) of West African forage resources can be correlated to in situ-measured reflectance data [14]. Oz et al. found informative spectral bands in three types of models—vegetation indices (VI), neural network (NN), and partial least squares (PLS) regression—for estimating leaf chlorophyll (Chl) and carotenoid (Car) contents of three unrelated tree species and to assess the accuracy of the models using a minimal number of bands [15]. Jesús et al. proposed a two-step approach to realize simultaneous LAI mapping over green and senescent croplands [16]. Dibyendu et al. demonstrated that total polyphenols of tea can be precisely estimated from a field spectroradiometer at the leaf level irrespective of age of the bushes and farming practices [17].
With the aid of Internet technology, agricultural informatization and agricultural big data have become inevitable trends. In recent years, with the successive launch of hyperspectral satellites and the development of microhyperspectral imagers of UAVs, the applications of hyperspectral remote sensing has become widely available. However, on the other hand, the increase of data volume brings great challenges with respect to data transmission, analysis, and storage [18]. Candes, Donoho, and Tao et al. proposed a new data acquisition and processing theory called compressive sensing (CS) [19,20,21]. Compressive sensing samples data at far below the Nyquist sampling rate by constructing an uncorrelated observation matrix, and the original data is reconstructed by a reconstruction algorithm. It, thus, provides a new way for compressing and reconstructing data with large volume.
At present, many studies have been conducted on applying compressive sensing in processing high-dimensional data. Kang et al. [22] proposed a method of distributed compressive sensing to grouping the video sequences efficiently by studying the correlation of the video sequences. Ly et al. [23] pointed out that the hyperspectral data should be stochastically separated by spectral and spatial partitioning. Chen et al. [24] proposed a sparse method for hyperspectral image target detection. Wang et al. proposed a pixel-based distributed compressive sensing [25], which divides the hyperspectral data into endmember extraction and abundance estimation through a linear mixture model. However, all of these methods did not concern applications of compressive sensing in agriculture, but mainly focused on the spatial reconstruction.
To promote the application of hyperspectral remote sensing in agriculture, the compressive sensing method provides a new method for the compression and recovery of hyperspectral data. However, at present, the research of compressive sensing is mainly focused on the reconstruction of the spatial image, and the spectral dimension of the hyperspectral data needs to be concerned in information reconstruction. Hyperspectral images have spectral-spatial correlations. The compressive sensing of high-dimensional data using the autocorrelation nature of data to improve the data sparse representation, which is able to reduce the complexity and improve the accuracy of reconstruction. In this study, through the analysis of plant spectral characteristics, a distributed spectral adaptive grouping compressive sensing is proposed and verified.

2. Methods

2.1. Spectral Adaptive Distributed Compressive Sensing

Figure 1 shows the flowchart of the proposed plant hyperspectral compressive sensing reconstruction algorithm. Firstly, the spectral characteristics of hyperspectral data are analyzed and the joint sparse model is constructed. Secondly, the spectral bands are adaptively grouped and the hyperspectral data are compressed and reconstructed on the basis of grouping, so as to determine the optimal grouping threshold. Then, to evaluate the reconstruction effect, AGDCS, OMP [26], and GPSR [27] are used to analyze the PSNR in the space and calculate errors among spectra.

2.2. Analysis of Plant Spectral Characteristics

In this study, a visible and near infrared (VIS-NIR) hyperspectral imaging system covering the spectral wavelengths of 380–1030 nm was used. The system includes a CCD camera (C8484-05, Hamamatsu, Hamamatsu city, Japan), an imaging spectrograph, a lens, two light sources provided by two 150 W quartz tungsten halogen lamps and V10E software (Isuzu Optics Corp., Taiwan) for the computer operating the spectral image system. The spectral resolution is 2.8 nm and the area CCD array detector of the camera has 6,726,512 pixels. The system scans the samples line by line, and the reflected light was dispersed by the spectrograph and captured by the area CCD array detector in spatial-spectral axes.
The data used in the experiment are hyperspectral images of 12 pieces of Camellia sinensis. Data from wavelengths of 450–850 nm and a total of 320 bands of data are chosen to perform the following experiments, and a single pixel is defined by a 12-bit unsigned integer. The images at 661 nm, 553 nm, and 449 nm from the original data are selected as the red, green, and blue channels of the false color composite RGB images, as shown in Figure 2a. Within each leaf, three regions of interest (ROIs) are marked with the size of 3 × 3 pixels. The averaged spectrum for each ROI is demonstrated in Figure 2b.
It can be observed from the spectral curve that despite the averaged spectral curves of different regions are different, the trend of the curves is consistent. The fluctuation of spectral reflectance values in 450–680 nm are relatively small. The spectrum in 680–760 nm rised rapidly. In 760–850 nm, the spectral reflectance is stablized at relatively high values, but the fluctuation within the local adjacent wave bands is relatively large. In order to further study on the distribution of spectral correlation of plant hyperspectral images, the correlation matrix of the hyperspectral bands is shown in Figure 3. It can be observed from the correlation coefficient diagram that the spectral reflectance is high in 450–680 nm and 760–850 nm, and the absolute value of autocorrelation and cross-correlation coefficient between the bands is above 0.9. The correlation between 680–760 nm is small, which is consistent with the variation of the spectral curve. Based on the qualitative analysis of the correlation characteristics of the plant spectrum, the following conclusions are obtained by analyzing the correlation coefficient (r) of each band. 59.41% bands are extremely correlated (|r| > 0.95); 6.93% bands are highly correlated (0.8 < |r| ≤ 0.95); 8.48% bands are moderately correlated (0.5 < |r| ≤ 0.8); whereas 25.18% bands are lowly correlated (0.3 < |r| ≤ 0.5) or not correlated (|r| ≤ 0.3).

2.3. Joint Sparse Model

According to the spectral characteristics of the hyperspectral images, there is a general correlation between the bands of the plant hyperspectral images and that more than half of the bands are highly correlated. The higher the correlation between the bands, less difference remained. Based on this, the present study proposed a joint sparse model of bands of plant hyperspectral images [22].
Assuming there are two bands Xi and Xi+1, the correlation between Xi and Xi+1 can be calculated, where Xi represents the current band. As the result of different wavelength reflection of the same object in the different bands of hyperspectral images, Xi and Xi+1 have the same spatial information. Meanwhile, Xi and Xi+1 have their own unique spectral information. The Xi and Xi+1 can be described as follows:
Xi = Xc + Xi_r
Xi+1 = Xc + Xi+1_r
where Xc is the same part of Xi and Xi+1, which refers to the same spatial information, Xi_r and Xi+1_r are their own unique parts, which mean the results of different reflections of different wavelengths. Xi is used as a reference to Xi+1, and the same spectral estimation (Xc) and the different information error coding (Xi+1_r) are used as the predictive value of Xi+1 in the spectral coding. The joint sparse model can be expressed as below:
Xi = ΨSi
Xi+1_r = ΨSi+1_r
where Si and Si+1_r are sparse representations of Xi and Xi+1_r, respectively, and Ψ is a canonical orthogonal matrix.

2.4. Distributed Compressive Sensing Based on Spectral Characteristics

In the distributed compressive sensing, hyperspectral band data are divided into a series of group of pictures (GOPs) bands. Each GOP consists of several bands which contain one key band and several non-key bands. Then, the sampling rates of the key band and non-key bands are realized by controlling the sampling matrix:
Yi = ΦXi = ΦΨSi
where Ψ is a canonical orthogonal matrix, Si is a sparse representation of the original signal in the transform domain, ΨHΨ = ΨΨH = I, and I is a unit matrix. The size of the Φ is M × N, M << N, which is the partial block Hadamard matrix [28]. Xi is a band of hyperspectral data, and Yi is an observed value.
Finally, the GPSR algorithm is used in the key band and the reconstruction of the non-key band is assisted by using the information from the key band. Compressive sensing reconstruction is the process of solving Equation (6). The solution of l0 norm is an NP-hard problem. The minimization problem of l1 norm is equivalent to that of l0 norm under certain conditions [29], so Equation (6) can be transformed to Equation (7). In this paper, GPSR and OMP are chosen to reconstruct Xi on the basis of the Yi:
x = arg   min x 0   st   Φ x = y
x = arg   min x 1   st   Φ x = y

2.5. Spectral Adaptive Grouping and Selection of Key Bands

Based on the above analysis, non-key bands can be reconstructed by the key band with side information assist. Therefore, it is very important to effectively group bands and select more effective key bands. In this study, PSNR is used as a basis for adaptive grouping of hyperspectral data and selecting key bands. The steps of adaptive grouping algorithm are as follows:
  • Step 1: Solve all PSNRs between the first band and each of the rest of the bands, and those in the rest of the bands and those whose PSNRs are greater than the threshold are all selected and classified into the group of the first band.
  • Step 2: Set up a new set from the remaining bands and repeat Step 1 to construct a new group.
  • Step 3: Repeat Step 2 until all bands are assigned to different groups.
In the adaptive band grouping algorithm, we can see that in every grouping, the PSNR values of the first band and other bands are greater than the threshold value. Given that the first band in the group has high similarity with the remaining bands, the first band of each group is determined as the key band and the other bands are non-key bands.

2.6. Error Evaluation Methods

The mean absolute percentage error (MAPE), the mean absolute error (MAE), and the root mean square error (RMSE) are used to evaluate the reconstructed effectiveness.
MAPE, also known as mean absolute percentage deviation (MAPD), is a measure of prediction accuracy of a forecasting method in statistics, and is defined by the formula:
M A P E = t = 1 n | A t F t A t | n
MAE is a quantity used to measure how close forecasts or predictions are to the eventual outcomes. The MAE is given by:
M A E = i = 1 n | y i x i | n
RMSE represents the sample standard deviation of the differences between predicted values and observed values. The RMSE is given by:
R M S E = t = 1 n ( y t y t ) 2 n

3. Results and Discussion

In the experiments, the software platform is MATLAB R2012a and the hardware platform is a Lenovo notebook computer in which the CPU is an Intel I3-2350M clocked at 2.3 GHz with 6 GB of memory.

3.1. The Results of Spectral Adaptive Threshold Grouping

3.1.1. The Results of Adaptive Grouping and the Different Sampling Rate of Key and Non-Key Bands

Key bands are sampled at a high rate and all sampling rates are from 0.1 bpp (bits per pixel) to 0.5 bpp. In case the sampling rate of the GPSR algorithm is over 0.4 bpp, the PSNR value changes slowly. The sampling rate of the key band is set to 0.5 bpp. The grouping results of non-key bands at different thresholds and the sampling rates of non-key bands are shown in Table 1.
It can be seen from Table 1, with the increase of the number of groupings, that the sampling rate of the non-key band decreases gradually, especially at a low sampling rate. When the threshold is set to 30 dB and the overall sampling rate is 0.1 bpp, the number of groupings will reach 65 and the sampling rate of the non-key band will be negative.

3.1.2. Analysis of Results of Adaptive Band Grouping Reconstruction

According to the sampling rate of the non-key bands calculated in Table 1, experimental results of the reconstructed PSNR of key bands, non-key bands, and all bands are shown as Figure 4, Figure 5 and Figure 6.
It can be seen from Figure 4 that the PSNR values of different groupings of key band reconstruction vary significantly, as the sampling rate of key bands is set to 0.5 bpp. To explain this phenomenon, more experiments are implemented for each band at a sampling rate of 0.5 bpp, as shown in Figure 7.
In Figure 7, when the sampling rate is 0.5 bpp, the reconstructed PSNR values of different bands are quite different. According to the spectral characteristic curve, the difference of the spectrum with a proximate reconstructed PSNR value is insignificant. The lower the PSNR threshold is, the more groups will emerge and the higher the average PSNR value of key bands will be, and vice versa.
As can be seen from Figure 5 and Figure 6, when the sampling rate is larger, the reconstructed PSNR values of non-key bands and all bands are both higher. More groupings indicate less efficiency of the overall reconstruction. In order to obtain good reconstructed PSNR for all sampling rates from 0.1 bpp to 0.5 bpp, the grouping PNSR threshold of 25 dB is chosen in the following experiments.

3.2. Spatial Domain Reconstruction Analysis

Twenty-five decibels was selected as the overall threshold for further analysis using OMP and GPSR algorithms. The reconstructed results are analyzed using the subjective evaluation and average peak signal to noise ratio in the spatial domain. To facilitate a visual comparison, 661 nm, 553 nm, and 449 nm are selected, respectively, as red, green, and blue channels to form the synthesized RGB image. Figure 8 illustrates the experimental results for OMP, GPSR, and AGDCS under the sampling rate from 0.1 bpp to 0.5 bpp. The fidelity of reconstructed images of all algorithms are significantly related with the sampling rate, especially at low sampling rates. The increasing of the sampling rate can significantly improve the subjective quality of the reconstructed image. Better subjective quality can be obtained using AGDCS compared with the other algorithms. The reconstructed subjective quality of AGDCS is very close to that of GPSR and significantly better than that of OMP at high sampling rates (great than or equal to 0.4 bpp). Experimental results show that the side information assists of key bands can improve the quality of reconstruction at low sampling rates.
It can be seen from Figure 9 that the reconstructed PSNRs of AGDCS and GPSR are significantly higher than that of OMP under sampling rates from 0.1 bpp to 0.5 bpp. At low sampling rates from 0.1 bpp to 0.3 bpp, the reconstructed PSNR of AGDCS is significantly higher than that of GPSR. At high sampling rates (greater than or equal to 0.4 bpp), the reconstructed PSNR of AGDCS approaches that of GPSR. It is shown that the joint sparse model improves the reconstructed PSNR in the spatial domain, especially at low sampling rates. Additionally, from the error bar of Figure 9, the standard deviation of AGDCS is close to that of GPSR.

3.3. Spectral Domain Reconstruction Results Analysis

The spectral analysis on vegetation is an important basis for the monitoring of crop growing status and crop stressors in precision farming. According to regions of interest (ROIs) in Figure 1, the average reconstructed spectral curves of ROIs using OMP, GPSR, and AGDCS are calculated at different sampling rates, from 0.1 bpp to 0.5 bpp, respectively.
As shown in Figure 10, Figure 11 and Figure 12, different regions have different sensitivity to compressive sensing reconstruction. There is a significant correlation between the reconstructed algorithms and the sampling rates, and the reconstructed effect is improved with the increase of the sampling rate. Additionally, the reconstructed curves of AGDCS are significantly better than those of the others for ROI2 and ROI3, and that of OMP is the lowest under different sampling rates for these three ROIs. For ROI1, part of reconstructed curves of AGDCS is worse than that of GPSR. The reconstructed effectiveness of the various methods is quantified by introducing MAPE, MAE, and RMSE as shown in Table 2. MAPE, MAE, and RMSE of AGDCS are better than those of the others when the sampling rate is less than 0.2 bpp. MAPE, MAE, and RMSE of AGDCS is comparable to those of GPSR when the sampling rate of greater than or equal to 0.3 bpp, and those of OMP are lower than those of the others.

3.4. Results of Spectral Indices of Physiological Properties

As it can been observed from Table 3, Table 4, Table 5 and Table 6, RMSE comparisons of the spectral indices of CRI550, CRI700, Dep550–750, and Area550–750, of different methods for three ROIs at different sampling rates are given. For spectral indices of CRI550, CRI700, Dep550–750, and Area550–750 for three ROIs at different sampling rates, the RMSE of OMP is relatively higher than that of the other methods. RMSEs of these four spectral indices of AGDCS are almost all less than those of GPSR for ROI2 and ROI3 at different sampling rates. For ROI1, the situation is different. The RMSEs of these four spectral indices of GPSR are sometimes less than those of AGDCS. These results are just like Figure 10 shows and, although part of the results for AGDCS are worse than GPSR, the overall results of AGDCS are better than GPSR.

3.5. Results of Average Reconstructed Time

The experimental results of average reconstructed time of all bands for different algorithms are shown in Figure 13. The reconstruction efficiency of the OMP algorithm is positively correlated with the sampling rate, while the reconstruction efficiency of the GPSR algorithm is negatively correlated with the sampling rate. The average reconstruction time of AGDCS is relatively stable. The average reconstruction time of GPSR is much higher than that of the others. As a result, AGDCS exhibited relatively high reconstructed efficiency.

3.6. Discussion

The OMP is a typical greedy algorithm whose reconstructed effect is worse than that of GPSR, which is a type of convex optimization algorithm. Thus, GPRS is chosen as the reconstructed algorithm for AGDCS to reconstruct plant hyperspectral data. AGDCS is based on the joint sparse model of plant hyperspectral images which can make full use of high inter-spectral correlation and the key bands to assist reconstruction of hyperspectral data. Unlike the AGDCS, GPSR does not take inter-spectral correlation into account, but only uses intra-spectral correlation to process hyperspectral data band by band. Additionally, the grouping strategy of AGDCS is to put those bands whose PSNR values are close to its key band into a group in which there is only one key band so the same spatial information of one group provides the primary information and different information errors become relatively low. However, GPSR processes each band as a key band. It leads to the reconstructed performance of AGDCS being quite good. This explains the phenomenon that the reconstructed performance of AGDCS is better than that of GPSR. For AGDCS, the reconstruction of the non-key bands plays the main role, whose main computational burden is to process the different information errors between the key band and the non-key band, which is significantly lower than that of the key-band. Thus, the reconstructed efficiency of AGDCE is extraordinarily higher than that of GPSR.

4. Conclusions

There is a high spectral correlation for plant hyperspectral data. The joint sparse model based on spectral characteristics can not only improve the fidelity of reconstructed plant hyperspectral images, but also effectively reduces the reconstructed error in spectral domain more efficiently. Under a relatively low sampling rate (less than 0.2 bpp), the PSNR of AGDCS is 13.72 dB higher than that of the OMP algorithm, and is 1.66 dB higher than that of the GPSR algorithm. For the errors in the spectral domain, the average normalized root mean square error, the mean absolute percentage error, and the mean absolute error of AGDCS decrease by 35.38%, 31.83%, and 33.33% than those of the GPSR algorithm, respectively. Additionally, AGDCS can achieve a relatively high reconstructed efficiency.
More extensive research can be studied in the future. For example, new compressive sensing algorithms can be involved in denoising the plant hyperspectral data. A high-performance reconstructed method should be proposed to improve the reconstructed efficiency.

Acknowledgments

This work was supported by the National Nature Science Foundation of China (NSFC) under grant No. 41671415 and 61205200, the National Key Foundation for Exploring Scientific Instruments of China under grants No. 61427808, and the Zhejiang Public Welfare Technology Application Research Project of China under grant No. 2016C32087. We are grateful to all of the subjects for their participation in this study.

Author Contributions

Ping Xu conceived of and designed the study. Junfeng Liu performed the experiments and wrote the paper. Bo Qiu helped to revise the paper. Lingyun Xue and Jingcheng Zhang supervised the work.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Tang, Y.; Huang, J. Study on Hyperspectral Remote Sensing in Agriculture. Remote Sens. Technol. Appl. 2001, 16, 248–251. [Google Scholar]
  2. Nuarsa, I.W.; Nishio, F.; Hongo, C.; Mahardika, I.G. Using variance analysis of multitemporal MODIS images for rice field mapping in Bali Province, Indonesia. Int. J. Remote Sens. 2012, 33, 5402–5417. [Google Scholar] [CrossRef]
  3. Tornosa, L.; Huescab, M.; Dominguezc, J.A.; Moyanoa, M.C.; Cicuendezd, V.; Recuerod, L.; Oruetad, A.P. Assessment of MODIS spectral indices for determining rice paddy agricultural practices and hydroperiod. J. Photogramm. Remote Sens. 2015, 101, 110–124. [Google Scholar] [CrossRef]
  4. Kang, Y.; Victoria, L.W.; Xinping, C.; Georg, B. Estimating leaf chlorophyll of barley at different growth stages using spectral indices to reduce soil background and canopy structure effects. J. Photogramm. Remote Sens. 2014, 97, 58–77. [Google Scholar]
  5. Fei, L.; Bodo, M.; Yuncai, H.; Xinping, C.; Urs, S. Optimising three-band spectral indices to assess aerial N concentration, N uptake and aboveground biomass of winter wheat remotely in China and Germany. J. Photogramm. Remote Sens. 2014, 92, 112–123. [Google Scholar]
  6. Atherton, J.; Nichol, C.J.; Porcar-Castell, A. Using spectral chlorophyll fluorescence and the photochemical reflectance index to predict physiological dynamics. Remote Sens. Environ. 2016, 176, 17–30. [Google Scholar] [CrossRef]
  7. Heli, L.; Chunjiang, Z.; Guijun, Y.; Haikuan, F. Variations in crop variables within wheat canopies and responses of canopy spectral characteristics and derived vegetation indices to different vertical leaf layers and spikes. Remote Sens. Environ. 2015, 169, 358–374. [Google Scholar]
  8. Qiang, L.; Shunlin, L.; Zhiqiang, X.; Hongliang, F. Retrieval of leaf area index using temporal, spectral, and angular information from multiple satellite data. Remote Sens. Environ. 2014, 145, 25–37. [Google Scholar]
  9. Mahlein, A.K.; Rumpf, T.; Welkeb, P.; Dehnea, H.W.; Plümerb, L.; Steinera, U.; Oerkea, E.C. Development of spectral indices for detecting and identifying plant diseases. Remote Sens. Environ. 2013, 128, 21–30. [Google Scholar] [CrossRef]
  10. Alicia, P.O.; Margarita, H.; Michael, L.W.; Javier, L.; Shruti, K.; Mónica, G.; Susan, L.U. Derivation of phenological metrics by function fitting to time-series of Spectral Shape Indexes AS1 and AS2: Mapping cotton phenological stages using MODIS time series. Remote Sens. Environ. 2012, 126, 148–159. [Google Scholar]
  11. Veraverbeke, S.; Harris, S.; Hooka, S. Evaluating spectral indices for burned area discrimination using MODIS/ASTER (MASTER) airborne simulator data. Remote Sens. Environ. 2011, 115, 2702–2709. [Google Scholar] [CrossRef]
  12. Jochem, V.; Juan, P.R.; Anatoly, G.; Jesus, D.; José, M.; Gustau, C.V. Spectral band selection for vegetation properties retrieval using Gaussian processes regression. Int. J. Appl. Earth Obs. Geoinf. 2016, 52, 554–567. [Google Scholar]
  13. Abderrazak, E.; Rachid, L.; Karem, C.; Jamal-eddine, O.; Mohamed, H.; El, M.B.; Abderrahmene, E.G. Spatiotemporal monitoring of soil salinization in irrigated Tadla Plain (Morocco) using satellite spectral indices. Int. J. Appl. Earth Obs. Geoinf. 2016, 50, 64–73. [Google Scholar]
  14. Ferner, J.; Linstädter, A.; Südekum, K.H.; Schmidtlein, S. Spectral indicators of forage quality in West Africa’s tropical savannas. Int. J. Appl. Earth Obs. Geoinf. 2015, 41, 99–106. [Google Scholar] [CrossRef]
  15. Oz, K.; Raphael, L.; Anatoly, G. Non-destructive estimation of foliar chlorophyll and carotenoid contents: Focus on informative spectral bands. Int. J. Appl. Earth Obs. Geoinf. 2015, 38, 251–260. [Google Scholar]
  16. Jesús, D.; Jochem, V.; Juan, P.R.; Antonio, R.V.; José, M. Brown and green LAI mapping through spectral indices. Int. J. Appl. Earth Obs. Geoinf. 2015, 35, 350–358. [Google Scholar]
  17. Dibyendu, D.; Prabir, K.D.; Uttam, K.B.; Upasana, S.; Shalini, S.; Jaswant, R.S.; Vinay, K.D. Retrieval of tea polyphenol at leaf level using spectral transformation and multi-variate statistical approach. Int. J. Appl. Earth Obs. Geoinf. 2015, 36, 22–29. [Google Scholar]
  18. Zhang, H.; Li, Z.; Zou, T.; Wei, X.; Yang, G. Overview of Agriculture Big Data Research. Comput. Sci. 2014, 41, 378–392. [Google Scholar]
  19. Candes, E.J.; Romberg, J.; Tao, T. Robust uncertainty principles: Exact signal reconstruction from highly incomplete frequency information. IEEE Trans. Inf. Theory 2006, 52, 489–509. [Google Scholar] [CrossRef]
  20. Donoho, D.L. Compressed sensing. IEEE Trans. Inf. Theory 2006, 52, 1289–1306. [Google Scholar] [CrossRef]
  21. Candes, E.J.; Tao, T. Near-optimal signal recovery from random projections: Universal encoding strategies. IEEE Trans. Inf. Theory 2006, 52, 5406–5425. [Google Scholar] [CrossRef]
  22. Kang, L.W.; Lu, C.S. Distributed compressive video sensing. In Proceedings of the IEEE International Conference on Acoustics, Speech and Signal Processing, Taipei, Taiwan, 19–24 April 2009; pp. 1169–1172. [Google Scholar]
  23. Ly, N.H.; Du, Q.; Fowler, J.E. Reconstruction from Random Projections of Hyperspectral Imagery with Spectral and Spatial Partitioning. IEEE J. Sel. Top. Appl. Earth Obs. Remote Sens. 2013, 6, 466–472. [Google Scholar] [CrossRef]
  24. Chen, Y.; Nasrabadi, N.M.; Tran, T.D. Sparse Representation for Target Detection in Hyperspectral Imagery. IEEE J. Sel. Top. Signal Process. 2011, 5, 629–640. [Google Scholar] [CrossRef]
  25. Wang, Z.; Feng, Y.; Xiao, H.; Wang, L. Distributed compressive sensing imaging and reconstruction of hyperspectral imagery. Opt. Precis. Eng. 2015, 23, 1131–1137. [Google Scholar] [CrossRef]
  26. Tropp, J.A.; Gilbert, A.C. Signal Recovery from Random Measurements via Orthogonal Matching Pursuit. IEEE Trans. Inf. Theory 2007, 53, 4655–4666. [Google Scholar] [CrossRef]
  27. Figueiredo, M.A.T.; Nowak, R.D.; Wright, S.J. Gradient Projection for Sparse Reconstruction: Application to Compressed Sensing and Other Inverse Problems. IEEE J. Sel. Top. Signal Process. 2007, 1, 586–597. [Google Scholar] [CrossRef]
  28. Gan, L.; Do, T.T.; Tran, T.D. Fast compressive imaging using scrambled block Hadamard ensemble. In Proceedings of the 16th European Signal Processing Conference, Lausanne, Switzerland, 25–29 August 2015; pp. 1–5. [Google Scholar]
  29. Tsaig, Y.; Donoho, D.L. Extensions of compressed sensing. Signal Process. 2006, 86, 549–571. [Google Scholar] [CrossRef]
Figure 1. Flowchart of plant hyperspectral compressive sensing reconstruction.
Figure 1. Flowchart of plant hyperspectral compressive sensing reconstruction.
Sensors 17 01322 g001
Figure 2. Spatial image of hyperspectral raw data and three regional spectral curves. (a) Original data space image; (b) The averaged spectrum for each ROI.
Figure 2. Spatial image of hyperspectral raw data and three regional spectral curves. (a) Original data space image; (b) The averaged spectrum for each ROI.
Sensors 17 01322 g002
Figure 3. Correlation coefficient diagram.
Figure 3. Correlation coefficient diagram.
Sensors 17 01322 g003
Figure 4. Average reconstructed PSNR of key bands for different groups.
Figure 4. Average reconstructed PSNR of key bands for different groups.
Sensors 17 01322 g004
Figure 5. Average reconstructed PSNR of non-key bands for different groups.
Figure 5. Average reconstructed PSNR of non-key bands for different groups.
Sensors 17 01322 g005
Figure 6. Average reconstructed PSNR of all bands for different groups.
Figure 6. Average reconstructed PSNR of all bands for different groups.
Sensors 17 01322 g006
Figure 7. Reconstructed PSNR of each band at a sampling rate of 0.5 bpp.
Figure 7. Reconstructed PSNR of each band at a sampling rate of 0.5 bpp.
Sensors 17 01322 g007
Figure 8. Comparison of reconstruction results for different algorithms at different sampling rates.
Figure 8. Comparison of reconstruction results for different algorithms at different sampling rates.
Sensors 17 01322 g008
Figure 9. Comparison of the reconstructed PSNR with different reconstruction methods.
Figure 9. Comparison of the reconstructed PSNR with different reconstruction methods.
Sensors 17 01322 g009
Figure 10. Comparison of reconstructed spectral curves of ROI1. (a) OMP; (b) GPSR; (c) AGDCS.
Figure 10. Comparison of reconstructed spectral curves of ROI1. (a) OMP; (b) GPSR; (c) AGDCS.
Sensors 17 01322 g010
Figure 11. Comparison of reconstructed spectral curves of ROI2. (a) OMP; (b) GPSR; (c) AGDCS.
Figure 11. Comparison of reconstructed spectral curves of ROI2. (a) OMP; (b) GPSR; (c) AGDCS.
Sensors 17 01322 g011
Figure 12. Comparison of reconstructed spectral curves of ROI3. (a) OMP; (b) GPSR; (c) AGDCS.
Figure 12. Comparison of reconstructed spectral curves of ROI3. (a) OMP; (b) GPSR; (c) AGDCS.
Sensors 17 01322 g012
Figure 13. Comparison of the average reconstructed time of all bands for different algorithms.
Figure 13. Comparison of the average reconstructed time of all bands for different algorithms.
Sensors 17 01322 g013
Table 1. Sampling rates of non-key bands in different groups.
Table 1. Sampling rates of non-key bands in different groups.
Threshold/dBGroupsSampling Rate/bpp
0.10.20.30.40.5
2060.0920.1940.2960.3980.500
2170.0910.1930.2960.3980.500
2270.0910.1930.2960.3980.500
2380.0900.1920.2950.3970.500
24150.0800.1850.2900.3950.500
25210.0720.1790.2860.3930.500
26290.0600.1700.2800.3900.500
27360.0490.1620.2750.3870.500
28450.0350.1510.2670.3840.500
29560.0150.1360.2580.3790.500
Table 2. Error analysis.
Table 2. Error analysis.
Error Analysis of Different AlgorithmsSampling Rate/bpp
0.10.20.30.40.5
MAPEOMP0.9683 0.5158 0.2861 0.2105 0.1598
GPSR0.1560 0.1006 0.0796 0.0689 0.0617
AGDCS0.1005 0.0814 0.0728 0.0680 0.0639
MAEOMP0.4599 0.2408 0.1434 0.1035 0.0787
GPSR0.0798 0.0529 0.0424 0.0369 0.0329
AGDCS0.0544 0.0451 0.0401 0.0370 0.0345
RMSEOMP0.5722 0.3003 0.1803 0.1303 0.0994
GPSR0.1065 0.0698 0.0556 0.0481 0.0428
AGDCS0.0710 0.0580 0.0519 0.0479 0.0447
Table 3. RMSE comparison of spectral indices of CRI550.
Table 3. RMSE comparison of spectral indices of CRI550.
Index Analysis of Different AlgorithmsSampling Rate/bpp
0.10.20.30.40.5
ROI1OMP15.1961.21624.31910.7681.027
GPSR0.4990.2090.2190.1870.267
AGDCS0.2620.6230.5140.4950.377
ROI2OMP1.6837.5040.77381.9520.778
GPSR0.6280.4810.3720.3590.292
AGDCS0.5350.4100.3450.3750.432
ROI3OMP0.6231.1161.4000.7070.538
GPSR0.5360.2250.2060.2090.200
AGDCS0.2590.2250.2210.2090.232
Table 4. RMSE comparison of spectral indices of CRI700.
Table 4. RMSE comparison of spectral indices of CRI700.
Index Analysis of Different AlgorithmsSampling Rate/bpp
0.10.20.30.40.5
ROI1OMP103.6981.24628.73313.5301.468
GPSR0.5430.2660.2500.2200.279
AGDCS0.23670.7270.5840.5380.395
ROI2OMP2.054164.5930.95696.0271.340
GPSR0.6010.4240.3150.3120.266
AGDCS0.5130.4050.3330.3500.403
ROI3OMP0.6910.9502.1810.7280.557
GPSR0.6670.3140.2450.2700.273
AGDCS0.3250.3170.2850.3140.327
Table 5. RMSE comparison of spectral indices of Dep550–750.
Table 5. RMSE comparison of spectral indices of Dep550–750.
Index Analysis of Different AlgorithmsSampling Rate/bpp
0.10.20.30.40.5
ROI1OMP0.4030.0840.2100.1160.065
GPSR0.0750.0170.0220.0150.009
AGDCS0.0320.0080.0060.0110.009
ROI2OMP0.0920.2910.0850.0950.043
GPSR0.0800.0490.0370.0320.025
AGDCS0.0570.0370.0290.0270.024
ROI3OMP0.1030.0830.1660.0520.044
GPSR0.0640.0350.0240.0240.024
AGDCS0.0320.0220.0270.0270.022
Table 6. RMSE comparison of spectral indices of Area550–750.
Table 6. RMSE comparison of spectral indices of Area550–750.
Index Analysis of Different AlgorithmsSampling Rate/bpp
0.10.20.30.40.5
ROI1OMP0.3620.0610.1060.0780.059
GPSR0.0930.0130.0220.0190.018
AGDCS0.0140.0310.0260.0300.026
ROI2OMP0.2440.2580.0840.0830.125
GPSR0.1140.0730.0530.0440.036
AGDCS0.0560.0370.0370.0330.033
ROI3OMP0.4560.1710.0950.1320.099
GPSR0.1120.0640.0390.0310.035
AGDCS0.0510.0360.0330.0270.028

Share and Cite

MDPI and ACS Style

Xu, P.; Liu, J.; Xue, L.; Zhang, J.; Qiu, B. Adaptive Grouping Distributed Compressive Sensing Reconstruction of Plant Hyperspectral Data. Sensors 2017, 17, 1322. https://doi.org/10.3390/s17061322

AMA Style

Xu P, Liu J, Xue L, Zhang J, Qiu B. Adaptive Grouping Distributed Compressive Sensing Reconstruction of Plant Hyperspectral Data. Sensors. 2017; 17(6):1322. https://doi.org/10.3390/s17061322

Chicago/Turabian Style

Xu, Ping, Junfeng Liu, Lingyun Xue, Jingcheng Zhang, and Bo Qiu. 2017. "Adaptive Grouping Distributed Compressive Sensing Reconstruction of Plant Hyperspectral Data" Sensors 17, no. 6: 1322. https://doi.org/10.3390/s17061322

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop