Next Article in Journal
Isolation of Chavibetol and Methyleugenol from Essential Oil of Pimenta pseudocaryophyllus by High Performance Liquid Chromatography
Next Article in Special Issue
Identification of Maize Kernel Vigor under Different Accelerated Aging Times Using Hyperspectral Imaging
Previous Article in Journal
Co-Delivery of Gemcitabine and Paclitaxel in cRGD-Modified Long Circulating Nanoparticles with Asymmetric Lipid Layers for Breast Cancer Treatment
Previous Article in Special Issue
Insight into Rapid DNA-Specific Identification of Animal Origin Based on FTIR Analysis: A Case Study
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Variety Identification of Raisins Using Near-Infrared Hyperspectral Imaging

1
College of Biosystems Engineering and Food Science, Zhejiang University, Hangzhou 310058, China
2
Key Laboratory of Spectroscopy Sensing, Ministry of Agriculture and Rural Affairs, Hangzhou 310058, China
3
College of Information Science and Technology, Shihezi University, Shihezi 832000, China
4
State Key Laboratory of Modern Optical Instrumentation, Zhejiang University, Hangzhou 310058, China
*
Authors to whom correspondence should be addressed.
Molecules 2018, 23(11), 2907; https://doi.org/10.3390/molecules23112907
Submission received: 12 October 2018 / Revised: 6 November 2018 / Accepted: 7 November 2018 / Published: 8 November 2018

Abstract

:
Different varieties of raisins have different nutritional properties and vary in commercial value. An identification method of raisin varieties using hyperspectral imaging was explored. Hyperspectral images of two different varieties of raisins (Wuhebai and Xiangfei) at spectral range of 874–1734 nm were acquired, and each variety contained three grades. Pixel-wise spectra were extracted and preprocessed by wavelet transform and standard normal variate, and object-wise spectra (sample average spectra) were calculated. Principal component analysis (PCA) and independent component analysis (ICA) of object-wise spectra and pixel-wise spectra were conducted to select effective wavelengths. Pixel-wise PCA scores images indicated differences between two varieties and among different grades. SVM (Support Vector Machine), k-NN (k-nearest Neighbors Algorithm), and RBFNN (Radial Basis Function Neural Network) models were built to discriminate two varieties of raisins. Results indicated that both SVM and RBFNN models based on object-wise spectra using optimal wavelengths selected by PCA could be used for raisin variety identification. The visualization maps verified the effectiveness of using hyperspectral imaging to identify raisin varieties.

1. Introduction

Raisins are generally consumed as snacks, and they are also served as popular ingredients in many other food menus. Raisins are dried grapes which are rich in dietary fiber, carbohydrates with a low glycemic index, and minerals like copper and iron, with a low fat content [1,2]. In addition to their nutritional value, they also have medical value, such as regulating blood pressure for individuals with mild increases in blood pressure [2,3,4]. In general, raisins are important commercial products for the grape industry.
The commercial value of raisins differs according to the production area. In China, Xinjiang Uygur Autonomous Region is one of the major producing regions of grape, the perfect producing conditions and climates make it quite suitable for grape planting and deep processing. Variety is another important factor which influences the taste and nutritional compositions of raisins. To satisfy the demands of producers and consumers, different varieties of grapes are developed. Variety is one of the important factors in pricing the raisins. Varieties of raisins can be identified by specialist, experienced famers, and laboratory-based chemical analysis methods. To improve the identification efficiency, advanced non-destructive methods have been introduced, among which computer vision, spectroscopy, and spectral imaging techniques have shown great efficiency and potential for large scale detection at industrial level. Ma et al. achieved rapid non-destructive identification of apple varieties with 96.67% accuracy based on hyperspectral imaging [5]. Zhang et al. identified coffee variety using mid-infrared transmittance spectroscopy combined with pattern recognition algorithm [6]. Yang et al. developed a model for maize seed variety identification based on hyperspectral imaging [7].
Hyperspectral imaging is a technique combining computer vision and spectroscopy. Images of the study objects can be acquired for image analysis, and spectral information can be extracted from each pixel within the image for spectral analysis. A combination of image analysis and spectral analysis can also be explored. Hyperspectral imaging has been widely used in food analysis [8,9], and it has showed great potential in the grape industry. Fernandes et al. estimated grape anthocyanin concentration using hyperspectral imaging data. The squared correlation coefficient value was 0.65 compared to the values measured using conventional laboratory techniques [10]. Rodríguez-Pulido et al. found it was possible to assess the maturation stage in grape seeds based on the near-infrared spectra with prediction models and multivariate analysis methods [11]. Zhao et al. used hyperspectral imaging to identify different varieties of grape seeds. The results indicated that the variety of each single grape seed was accurately identified with 94.3% accuracy of the calibration set and 88.7% accuracy of the prediction set [12].
The general application of hyperspectral images is to conduct data analysis on a predefined region of interest (ROI) [12,13,14]. Spectral information is most widely used in hyperspectral image analysis, due to the advantage that spectral information can be precisely extracted from each pixel within ROIs. In general, pixel-wise spectra are averaged to build calibration models, and some researchers have focused on using pixel-wise spectra to build calibration models [15,16,17,18]. In fact, the size and the shape of raisins, which are key factors for the classification of different varieties, also play an important role in raisin grading within one variety. The raisin size can be influenced by the harvesting procedure of fresh grapes and the air-drying procedure, which can also beget irregular shapes of raisins in addition to storage.
The objective of this study is to explore the feasibility of using near-infrared hyperspectral imaging to identify raisin varieties. The specific objectives are: (1) exploring the influence of fruit size and shape in classification accuracy; (2) exploring spectral preprocessing of standard normal variate (SNV) in classification accuracy; (3) comparing performances of objective-wise analysis and pixel-wise analysis of SVM (Support Vector Machine), k-NN (k-nearest Neighbors Algorithm), and RBFNN (Radial Basis Function Neural Network) models.

2. Results and Discussion

2.1. Spectral Profiles

In this research, 200 wavelength variables ranging from 975 to 1646 nm of hyperspectral images were studied. Figure 1 presents average spectra of each grade of raisins of Wuhebai (WHB) and Xiangfei (XF) with standard deviation (SD) at peaks and valleys (1123, 1210, 1308, and 1473 nm). The absorbance bands at 1123, 1210, and 1308 nm are largely attributed to the C–H stretching mode and overtone [19]. The wavelength around 1473 nm is a characteristic water wavelength [20]. It was obvious that a large proportion of overlap exists among eight curves, so it was necessary to conduct further study to make a better distinction between the two varieties of raisins.

2.2. PCA Scores Image Visualization

Pixel-wise PCA scores could be used to depict the PCA scores image. The first seven PCs explained over 99% of the total variance. Figure 2 shows visualized hyperspectral images of the first seven principal components (PC1–PC7) of two varieties of raisins. As can be seen from Figure 2, the warm color (yellow-red) accounted for the majority in WHB scores image of PC1 and PC2. In contrast to WHB scores image, the cold color (green-blue) was more obvious in XF scores image of PC1 and PC2, which revealed differences between two varieties. PCA scores image of PC5, PC6 and PC7 of XF exhibited obvious difference in color for Grade1 and two other grades, which showed differences among different grades. Although the PCA scores image could be used to distinguish different varieties and grades of raisins to some extent, it was necessary to conduct further study in order to obtain satisfactory classification results.

2.3. Effective Wavelength Selection

PCA loadings were used to select effective wavelengths for raisin cultivars classification. Since the first seven PCs explained over 99% of total variance, the loadings of these PCs were used. To examine the differences of object-wise analysis (average spectra) and pixel-wise analysis (pixel-wise spectra), PCA was conducted on both object-wise spectra and pixel-wise spectra of two different varieties of raisins. Figure 3 and Table 1 shows 20 or 17 optimal wavelengths selected by PCA based on object-wise analysis or pixel-wise analysis, respectively. PCA loadings plots of object-wise spectra and pixel-wise spectra were quite similar. As shown in Table 1, corresponding optimal wavelengths for object-wise spectra and pixel-wise spectra were nearly the same, with slight differences caused by different varieties and grades.
ICA was also conducted on object-wise spectra and pixel-wise spectra. To compare with PCA, the same numbers of optimal wavelengths are showed in Table 2.

2.4. Raisin Variety Classification Models Based on Different Grades

The PCA analysis above indicated that there were differences between different varieties of raisins, and there were also differences among different grades of raisins. Thus, grade was an important factor which influenced classification results of two varieties of raisins.
To evaluate the influence of different grades on model performances, samples from the same grade of raisins were formed as calibration set, and the remaining samples were used as prediction set. SVM (Support Vector Machine) models were built using effective wavelengths selected by PCA, and the results are shown in Table 3.
For the calibration model built based on Grade1, classification results of the calibration set were good for both varieties and prediction results of three grades of WHB were good, while prediction results of XF were poor. There were no correctly classified samples for Grade3 of XF. The SVM model built based on Grade2 exhibited better performance compared with SVM model based on Grade1. For the calibration model built based on Grade2, both the calibration set and prediction set obtained satisfactory results, and WHB and XF both obtained good classification results. When the calibration set was built according to Grade3, classification results of calibration set were good. The prediction results of the three grades of WHB were good, and the prediction results of Grade2 and Grade3 of XF were also good. However, the prediction result of Grade1 was poor with classification accuracy lower than 20%.
These results revealed that different grades of raisins had influences on variety classification. As shown in Table 3, classification models based on Grade1 and Grade3 obtained poorer classification results compared with classification models based on Grade2. The reason might be that feature differences between Grade1 and Grade3 of XF were significant (for example the difference of sample size of different grades for same raisin variety was obvious as shown in Figure 4). The calibration set using Grade1 or Grade3 might not cover enough sample features used for PCA scores analysis.
Table 4 shows the results of SVM models built based on optimal wavelengths selected by ICA. Compared with Table 3, the calibration and prediction accuracies of SVM models based on Grade1 and Grade3 using optimal wavelengths selected by ICA were close to SVM models using optimal wavelengths selected by PCA. The accuracies based on Grade2 using optimal wavelengths selected by ICA were lower than SVM models based on PCA optimal wavelengths selection for both varieties of raisins.

2.5. Classification Results of Pixel-Wise and Object-Wise Spectra

According to Zhang et al. (2018) [21], pixel-wise spectra can be used to build classification models, and can achieve good prediction results on sample average spectra. The results of average spectra showed that sample size might be a factor influencing classification results. The advantage of hyperspectral imaging was to obtain spectral information of each pixel within the sample. Previous results have showed that pixel-wise spectra analysis has great value in hyperspectral image analysis [22,23]. For each variety of raisins, pixel-wise spectra were extracted. In all, there were about 300,000 pixels of each grade within the calibration sets of WHB and XF raisins. Establishing calibration models using such a great number of pixels (over 1,800,000 pixels) requires quite heavy computation. Selecting effective wavelengths could reduce the data volume significantly, but the data volume was still large. As for about 300,000 pixels of each grade of raisins, there might be redundant pixels for modelling, so representative pixels should be selected to reduce the data amount.
To select representative pixel-wise spectra, a calibration set selection procedure was proposed by Kang et al. (2004) [24]. Firstly, for pixel-wise spectra of all grades of a variety, the collected pixel-wise spectra were clustered into 3000 groups using the k-means algorithm. Secondly, the Euclidean distance between sample and group centroid was calculated, and the sample with smallest Euclidean distance was selected into the new calibration set.
SVM, k-NN, and RBFNN models were built using selected pixel-wise spectra or object-wise spectra, and prediction set was also formed based on selected pixel-wise spectra or object-wise spectra. The results are shown in Table 5. The value of sensitivity means the classification accuracy of WHB, and the value of sensitivity means the classification accuracy of XF.
From Table 5, the results of SVM and RBFNN models using pixel-wise spectra to predict pixel-wise spectra were acceptable, with classification accuracy of calibration and prediction about 80%–90%. Meanwhile, the results of SVM and RBFNN models using pixel-wise spectra to predict object-wise spectra also obtained good results for calibration set, with 93.62% and 88.40% accuracy, respectively. Compared with SVM and RBFNN models based on pixel-wise spectral, the results of k-NN were slightly poorer, with accuracies varied from 40%–90%. Three models based on object-wise spectra all obtained acceptable results for the calibration set, with accuracies ranging from 87%–99%. SVM, k-NN and RBFNN models using object-wise spectra to predict object-wise spectra obtained better results for the prediction set compared with models using pixel-wise spectra to predict object-wise spectra, with all accuracies above 80%.
The overall results indicated that SVM and RBFNN models using object-wise spectra to predict object-wise spectra could be used to identify raisin varieties. Selection of representative samples was of significance for stable and accurate models, which should be further studied.
Table 6 shows the classification results for SVM, k-NN, and RBFNN models based on optimal wavelengths selected by ICA. The accuracies based on models using pixel-wise spectra to predict pixel-wise spectra or object-wise spectra were much lower than the same models using optimal wavelengths selected by PCA. The prediction set results of models using object-wise spectra to predict pixel-wise spectra were poor, with accuracies varying from 48%–63%. Models using object-wise spectra to predict object-wise spectra obtained acceptable results, with all accuracies above 90%. However, the calibration results of three models using object-wise spectra as the calibration set using optimal wavelengths selected by ICA were slightly lower than the results of three models using optimal wavelengths selected by PCA.

2.6. Prediction Maps of Raisin Variety Detection

Based on the developed models, prediction maps of different raisins varieties could be formed. According to the results in Table 5, we used the pixel-wise SVM model to form prediction maps. Raisin grades of the corresponding pixel were predicted to form classification maps, and prediction maps are shown in Figure 4. Pixel-wise prediction maps indicate that most of the pixels could be correctly classified. The prediction maps show clear a difference between WHB and XF according to different visualized prediction color.

3. Materials and Methods

3.1. Sample Preparation

Two varieties of raisins, including Wuhebai (WHB) and Xiangfei (XF), were collected from a local market in Shihezi, Xinjiang Uygur Autonomous Region, China. For each variety of raisin, three grades (Grade1-large size, Grade2-medium size, Grade3-small size) of raisins were manually collected according to the raisin size. For each grade, 450 g samples were divided into 30 groups (nearly 15 g per group) to acquire 30 hyperspectral images. Two varieties of raisins were all produced in 2017. RGB images of the two varieties of raisins are shown in Figure 5.

3.2. Hyperspectral Imaging System

The hyperspectral image acquisition was carried out using an assembled near-infrared hyperspectral imaging system with the spectral range of 975–1646 nm. The system consisted of an ImSpector N17E imaging spectrograph (Spectral Imaging Ltd., Oulu, Finland), a Xeva 992 camera (Xenics Infrared Solutions, Leuven, Belgium) installed with an OLES22 lens (Spectral Imaging Ltd., Oulu, Finland), two 150 W tungsten halogen lamps (3900 Lightsource, Illumination Technologies Inc., Elbridge, NY, USA) that were symmetrically placed and served as the light source, and a conveyer belt (Isuzu Optics Corp., Taiwan). The imaging system was controlled by the software (Xenics N17E, Isuzu Optics Corp., Taiwan), which can be used to calibrated and analyze the images as well. The sketch of the hyperspectral imaging system is presented in Figure 6.

3.3. Hyperspectral Image Acquisition and Correction

To acquire hyperspectral images, the distance between sample plane and the camera was set to 16 cm, the moving speed of the plate was set to 13.5 mm/s, and the exposure time of the camera was set to 4 ms. After adjustment, white reference image was collected by a white Teflon bar whose reflectivity is approximately 100%, and dark reference image was acquired by turning off the light source and covering the lens with lens cap whose reflectivity is about 0%. The white and dark reference images were used to calibrate the light intensity and reduce the dark current. For each group, one image was acquired. In all, 30 images were acquired for each grade of raisin.
After raw hyperspectral images acquisition, hyperspectral images were then corrected by the following equation:
I c = I r I d I w I d
where Ic is the corrected image, Ir is the raw image, Iw is the white reference image and Id is the dark reference image.

3.4. Spectral Data Preprocessing and Extraction

The hyperspectral image at 1119 nm was selected for background segmentation since the reflectance difference between sample and background was more obvious. To differentiate background from foreground, we set the segmentation threshold to 0.122 for hyperspectral image binarization.
After image correction, spectral information was extracted from hyperspectral images. Each raisin kernel was defined as the region of interest (ROI). Pixel-wise spectra within the ROI were firstly extracted, and wavelet transform (WT) was used for smoothing. Wavelet function of Daubechies 7 with a decomposition level of 3 was used to reduce random noises. After WT preprocessing, standard normal variate (SNV) was used to reduce the influence of scattering of pixel-wise spectra [25]. Then, average spectra calculated according to pixel-wise spectra within each ROI were used to represent the sample. In this study, pixel-wise spectra and average spectra were both used for analysis.
To extract spectral information, a binary image was obtained using the gray-scale image at 1199 nm, the background was set as 0 and the kernel regions were set as 1. The binary image was applied to the gray-scale images at each wavelength, and the background information was thus removed. After the background removal, pixel-wise spectra were extracted and preprocessed.

3.5. Sample Set Division

For each grade of raisin, 30 hyperspectral images were acquired. Hyperspectral images were randomly split into the calibration set and prediction set at the ratio of 2:1, with 20 hyperspectral images in the calibration set and 10 hyperspectral images in the prediction set for each grade.

3.6. Data Analysis Methods

3.6.1. Principal Component Analysis

Principal component analysis (PCA) was used to explore the qualitative differences among different varieties of raisins [11,26,27,28]. In hyperspectral images, object-wise analysis and pixel-wise analysis were studied. For object-wise analysis, the average spectrum of each raisin kernel was used to conduct PCA; for pixel-wise analysis, pixel-wise spectra were used to conduct PCA. To explore the differences among raisins, the samples in the calibration set were used to conduct PCA. Then, scores values of each PC were then assigned to each kernel or each pixel to form the PCA scores image.
Hyperspectral imaging suffers from the large volume of data, and effectively reducing the data volume is of significance for data processing. There are also collinearity and redundancy in the spectra, which will affect the data analysis procedure. Variable selection is an effective strategy to reduce the data volume and select informative wavelengths. In this study, PCA loadings were used to select effective wavelengths. Loadings of each principal component (PC) indicate the correlation between the original variables and new feature variables. The higher the loading value is, the more important the variable is. The wavelengths with high absolute loading values can be selected as effective variables.

3.6.2. Independent Component Analysis

Independent component analysis (ICA) is a technique which is widely used in feature selection and feature extraction. It extracts independent source signals which are statistical independent by linear or nonlinear transformation. Independent component (IC) is obtained by a high-order statistic. Given a spectral matrix X, X can be expressed as Equation (2):
X = A s
where s are the independent components (ICs) and A is the mixing matrix. For spectral data matrix X, s is unknown, and the general procedure is to find the estimation of s by the following equation:
s = W X
where s is the estimation of s and W is the weight matrix for unmixing.
The procedure to select optimal wavelengths is as follow [29]. The average absolute weight value of each variable in W is calculated, and the variables with larger average absolute weight values are selected as optimal wavelengths. To compare with PCA, the same number of optimal wavelengths was selected by ICA. Fast ICA proposed by Hyvarinen and Oja was used to perform ICA in this study [30].

3.6.3. Discriminant Models

Support vector machine (SVM) is used to build models to classify different varieties of raisins. SVM is a supervised machine learning method, which is efficient to deal with linear and nonlinear data for classification and regression. For classification issues, SVM maps the original data into new feature spaces [31,32,33]. According to linearly separable data, a simple linear classifier can be constructed. For non-linearly separable data, the original data should be mapped into high-dimensional feature spaces so that the non-linearly separable issue can be transformed to a linearly separable issue. Kernel functions are the key for the mapping. Radial basis function (RBF) is a widely used kernel function with good performances for nonlinear data, and it was used as kernel function in this study. To conduct SVM with RBF kernel function, model penalty coefficient (C) and kernel parameter (γ) were determined by a grid-search procedure. In this study, the range of C and γ was 2−8–28.
The k-nearest neighbors algorithm (k-NN) is a type of instance-based learning method used for classification and regression [34,35]. Both for classification and regression, a useful technique can be used to assign weight to the contributions of the neighbors, so that the nearest neighbors contribute more to the average than the more distant ones. The k-NN algorithm is among the simplest of all machine learning algorithms.
Radial basis function neural network (RBFNN) is an efficient feedforward neural network, which has the best approximation performance and global optimal characteristics that outperforms other feedforward networks, and has a simple structure as well as a fast training speed. On the other hand, it is also a neural network model that is widely used in pattern recognition, nonlinear function approximation, and other fields [36,37].

3.6.4. Software and Model Evaluation

The performance of classification models was evaluated by the classification accuracy, specificity, and sensitivity [38]. Hyperspectral images analysis, spectral data extraction, spectral preprocessing, PCA analysis, SVM, K-NN, and RBFNN were conducted on Matlab R2014b (The MathWorks, Natick, MA, USA).

3.6.5. Visualization of Prediction Maps

One of the advantages of hyperspectral imaging is that prediction maps can be formed to visualize the distribution of physical and chemical features. Object-wise or pixel-wise calibration models using spectra extracted from the hyperspectral images can be used to predict object-wise or pixel-wise features, and prediction maps can be formed with the predicted values [6,13,39].

4. Conclusions

Hyperspectral imaging was successfully used to identify different varieties of raisins. Three grades of raisins of Wuhebai and Xiangfei were studied. Object-wise and pixel-wise spectra were extracted. PCA analysis was firstly conducted to form PCA scores images, and scores images of the first seven PCs indicated the differences between different varieties and among different grades. PCA and ICA of object-wise spectra and pixel-wise spectra were conducted to select effective wavelengths. The overall results indicates that SVM models and RBFNN models using object-spectra to predict object-spectra based on optimal wavelengths selected by PCA both obtained acceptable results. The overall results showed that hyperspectral imaging was an effective technique to identify raisin varieties, and that both pixel-wise and object-wise could be used to build classification models. Selection of representative samples was important for building a stable and accurate model, and how to select representative samples should be studied in the future.

Author Contributions

Conceptualization, L.F. and C.Z.; Data curation, Y.H.; Formal analysis, P.G.; Funding acquisition, L.F., Y.B., P.G. and Y.H.; Investigation, S.Z.; Methodology, S.Z. and Y.B.; Project administration, L.F. and P.G.; Resources, P.G.; Software, Y.B., P.G. and Y.H.; Supervision, Y.H.; Validation, L.F., C.Z. and Y.H.; Visualization, L.F.; Writing—original draft, L.F. and S.Z.; Writing—review & editing, C.Z. and Y.H.

Funding

This research was funded by National key R&D program of China, grant number 2018YFD0101002; National Natural Science Foundation of China, grant number 31871526 and Major scientific and technological tackling project of Shihezi University, grant number GXJS2015-ZDGG08.

Acknowledgments

We thank Hongji Lin in College of electrical engineering, Zhejiang University for editing and improving the manuscript.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Williamson, G.; Carughi, A. Polyphenol content and health benefits of raisins. Nutr. Res. 2010, 30, 511–519. [Google Scholar] [CrossRef] [PubMed]
  2. Margaret, J.S.; Xinyue, W.; Tiffany, H.; James, E.P. A Comprehensive review of raisins and raisin components and their relationship to human health. J. Nutr. Health 2017, 50, 203. [Google Scholar] [CrossRef] [Green Version]
  3. Kanellos, P.T.; Kaliora, A.C.; Gioxari, A.; Christopoulou, G.O.; Kalogeropoulos, N.; Karathanos, V.T. Absorption and bioavailability of antioxidant phytochemicals and increase of serum oxidation resistance in healthy subjects following supplementation with raisins. Plant Food Hum. Nutr. 2013, 68, 411–415. [Google Scholar] [CrossRef] [PubMed]
  4. Bays, H.E.; Schmitz, K.; Christian, A.; Ritchey, M.; Anderson, J. Raisins and Blood Pressure: A Randomized, Controlled Trial. J. Am. Coll. Cardiol. 2012, 59, E1721. [Google Scholar] [CrossRef]
  5. Huiling, M.; Wang, R.; Cheng, C.; Dong, W. Rapid Identification of Apple Varieties Based on Hyperspectral Imaging. Trans. CSAE 2017, 48, 305–312. [Google Scholar] [CrossRef]
  6. Zhang, C.; Wang, C.; Liu, F.; He, Y. Mid-Infrared Spectroscopy for Coffee Variety Identification: Comparison of Pattern Recognition Methods. J. Spectrosc. 2016, 2016, 1–7. [Google Scholar] [CrossRef]
  7. Yang, S.; Zhu, Q.B.; Huang, M.; Qin, J.W. Hyperspectral Image-Based Variety Discrimination of Maize Seeds by Using a Multi-Model Strategy Coupled with Unsupervised Joint Skewness-Based Wavelength Selection Algorithm. Food Anal. Method 2017, 10, 1–10. [Google Scholar] [CrossRef]
  8. Bao, Y.; Liu, F.; Kong, W.; Sun, D.W.; He, Y.; Qiu, Z. Measurement of Soluble Solid Contents and pH of White Vinegars Using VIS/NIR Spectroscopy and Least Squares Support Vector Machine. Food Bioprocess Technol. 2014, 7, 54–61. [Google Scholar] [CrossRef]
  9. Wang, H.; Peng, J.; Xie, C.; Bao, Y.; He, Y. Fruit Quality Evaluation Using Spectroscopy Technology: A Review. Sensors 2015, 15, 11889. [Google Scholar] [CrossRef] [PubMed]
  10. Fernandes, A.M.; Oliveira, P.; Moura, J.P.; Oliveira, A.A.; Falco, V.; Correia, M.J.; Melopinto, P. Determination of anthocyanin concentration in whole grape skins using hyperspectral imaging and adaptive boosting neural networks. J. Food Eng. 2011, 105, 216–226. [Google Scholar] [CrossRef]
  11. Rodríguez-Pulido, F.J.; Barbin, D.F.; Sun, D.W.; Gordillo, B.; González-Miret, M.L.; Heredia, F.J. Grape seed characterization by NIR hyperspectral imaging. Postharvest Biol. Technol. 2013, 76, 74–82. [Google Scholar] [CrossRef]
  12. Zhao, Y.; Zhang, C.; Zhu, S.; Gao, P.; Feng, L.; He, Y. Non-Destructive and Rapid Variety Discrimination and Visualization of Single Grape Seed Using Near-Infrared Hyperspectral Imaging Technique and Multivariate Analysis. Molecules 2018, 23, 1352. [Google Scholar] [CrossRef] [PubMed]
  13. Zhang, C.; Wang, Q.; Liu, F.; He, Y.; Xiao, Y. Rapid and non-destructive measurement of spinach pigments content during storage using hyperspectral imaging with chemometrics. Measurement 2016, 97, 149–155. [Google Scholar] [CrossRef]
  14. Kong, W.; Zhang, C.; Cao, F.; Liu, F.; Luo, S.; Tang, Y.; He, Y. Detection of Sclerotinia Stem Rot on Oilseed Rape (Brassica napus L.) Leaves Using Hyperspectral Imaging. Sensors 2018, 18, 1764. [Google Scholar] [CrossRef]
  15. Xing, J.; Symons, S.; Shahin, M.; Hatcher, D. Detection of sprout damage in Canada Western Red Spring wheat with multiple wavebands using visible/near-infrared hyperspectral imaging. Biosyst. Eng. 2010, 106, 188–194. [Google Scholar] [CrossRef]
  16. Williams, P.J.; Kucheryavskiy, S. Classification of maize kernels using NIR hyperspectral imaging. Food Chem. 2016, 209, 131–138. [Google Scholar] [CrossRef]
  17. Arngren, M.; Hansen, P.W.; Eriksen, B.; Larsen, J.; Larsen, R. Analysis of Pregerminated Barley Using Hyperspectral Image Analysis. J. Agric. Food Chem. 2011, 59, 11385–11394. [Google Scholar] [CrossRef] [PubMed]
  18. Yao, H.; Hruska, Z.; Kincaid, R.; Brown, R.L.; Bhatnagar, D.; Cleveland, T.E. Hyperspectral image classification and development of fluorescence index for single corn kernels infected with Aspergillus flavus. Trans. ASABE 2013, 56, 1977–1988. [Google Scholar] [CrossRef]
  19. Khodabux, K.; L’Omelette, M.S.S.; Jhaumeer-Laulloo, S.; Ramasami, P.; Rondeau, P. Chemical and near-infrared determination of moisture, fat and protein in tuna fishes. Food Chem. 2007, 102, 669–675. [Google Scholar] [CrossRef]
  20. Kinoshita, K.; Miyazaki, M.; Morita, H.; Vassileva, M.; Tang, C.; Li, D.; Ishikawa, O.; Kusunoki, H.; Tsenkova, R. Spectral pattern of urinary water as a biomarker of estrus in the giant panda. Sci. Rep. 2012, 2, 856. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  21. Zhang, C.; Liu, F.; He, Y. Identification of coffee bean varieties using hyperspectral imaging: Influence of preprocessing methods and pixel-wise spectra analysis. Sci. Rep. 2018, 8, 2166. [Google Scholar] [CrossRef] [PubMed]
  22. Diezma, B.; Lleó, L.; Roger, J.M.; Herrero-Langreo, A.; Lunadei, L.; Ruiz-Altisent, M. Examination of the quality of spinach leaves using hyperspectral imaging. Postharvest Biol. Technol. 2013, 85, 8–17. [Google Scholar] [CrossRef] [Green Version]
  23. Lara, M.A.; Lleó, L.; Diezma-Iglesias, B.; Roger, J.M.; Ruiz-Altisent, M. Monitoring spinach shelf-life with hyperspectral image through packaging films. J. Food Eng. 2013, 119, 353–361. [Google Scholar] [CrossRef] [Green Version]
  24. Kang, J.; Ryu, K.R.; Kwon, H.C. Using Cluster-Based Sampling to Select Initial Training Set for Active Learning in Text Classification. In Pacific-Asia Conference on Knowledge Discovery and Data Mining; Springer: Berlin/Heidelberg, Germany, 2004; pp. 384–388. [Google Scholar]
  25. Rinnan, Å.; Berg, F.V.D.; Engelsen, S.B. Review of the most common pre-processing techniques for near-infrared spectra. Trac-Trend Anal. Chem. 2009, 28, 1201–1222. [Google Scholar] [CrossRef]
  26. Choudhary, R.; Mahesh, S.; Paliwal, J.; Jayas, D.S. Identification of wheat classes using wavelet features from near infrared hyperspectral images of bulk samples. Biosyst. Eng. 2009, 102, 115–127. [Google Scholar] [CrossRef]
  27. Sun, J.; Jiang, S.; Mao, H.; Wu, X.; Li, Q. Classification of Black Beans Using Visible and Near Infrared Hyperspectral Imaging. Int. J. Food Sci. Technol. 2016, 19, 1687–1695. [Google Scholar] [CrossRef]
  28. Huang, M.; He, C.; Zhu, Q.; Qin, J. Maize Seed Variety Classification Using the Integration of Spectral and Image Features Combined with Feature Transformation Based on Hyperspectral Imaging. Appl. Sci. 2016, 6, 183. [Google Scholar] [CrossRef]
  29. Du, H.; Qi, H.; Wang, X.; Ramanath, R. Band selection using independent component analysis for hyperspectral image processing. In Proceedings of the Applied Imagery Pattern Recognition Workshop, Washington, DC, USA, 15–17 October 2003; IEEE: New York, NY, USA, 2003; pp. 93–98. [Google Scholar]
  30. Hyvärinen, A.; Oja, E. Independent component analysis: Algorithms and applications. Neural Netw. 2000, 13, 411–430. [Google Scholar] [CrossRef]
  31. Feng, X.; Zhao, Y.; Zhang, C.; Cheng, P.; He, Y. Discrimination of Transgenic Maize Kernel Using NIR Hyperspectral Imaging and Multivariate Data Analysis. Sensors 2017, 17, 1894. [Google Scholar] [CrossRef] [PubMed]
  32. Dumont, J.; Hirvonen, T.; Heikkinen, V.; Mistretta, M.; Granlund, L.; Himanen, K.; Fauch, L.; Porali, I.; Hiltunen, J.; Keski-Saari, S. Thermal and hyperspectral imaging for Norway spruce (Picea abies) seeds screening. Comput. Electron. Agric. 2015, 116, 118–124. [Google Scholar] [CrossRef]
  33. Lee, H.; Kim, M.S.; Song, Y.R.; Oh, C.S.; Lim, H.S.; Lee, W.H.; Kang, J.S.; Cho, B.K. Non-destructive evaluation of bacteria-infected watermelon seeds using visible/near-infrared hyperspectral imaging. J. Sci. Food Agric. 2016, 97, 1084. [Google Scholar] [CrossRef] [PubMed]
  34. Guo, G.; Wang, H.; Bell, D.; Bi, Y.; Greer, K. KNN Model-Based Approach in Classification. In Proceedings of the Otm Confederated International Conferences on the Move to Meaningful Internet Systems, Catania, Sicily, Italy, 3–7 November 2003; Springer: Berlin, Germany, 2003; pp. 986–996. [Google Scholar]
  35. Hall, P.; Park, B.U.; Samworth, R.J. Choice of Neighbor Order in Nearest-Neighbor Classification. Ann. Stat. 2008, 36, 2135–2152. [Google Scholar] [CrossRef]
  36. Shcherbakov, V.; Larsson, E. Radial basis function partition of unity methods for pricing vanilla basket options. Comput. Math. Appl. 2016, 71, 185–200. [Google Scholar] [CrossRef] [Green Version]
  37. Fornberg, B.; Flyer, N.; Hovde, S.; Piret, C. Locality properties of radial basis function expansion coefficients for equispaced interpolation. IMA J. Numer. Anal. 2018, 28, 121–142. [Google Scholar] [CrossRef]
  38. Baratloo, A.; Hosseini, M.; Negida, A.; El Ashal, G. Part 1: Simple Definition and Calculation of Accuracy, Sensitivity and Specificity. Emergency 2015, 3, 48–49. [Google Scholar] [PubMed]
  39. He, J.; Zhang, C.; He, Y. Application of Near-Infrared Hyperspectral Imaging to Detect Sulfur Dioxide Residual in the Fritillaria thunbergii Bulbus Treated by Sulfur Fumigation. Appl. Sci. 2017, 7, 77. [Google Scholar] [CrossRef]
Sample Availability: Samples of the compounds are available from the authors.
Figure 1. Average spectra with standard deviation (SD) of Wuhebai (WHB) and Xiangfei (XF).
Figure 1. Average spectra with standard deviation (SD) of Wuhebai (WHB) and Xiangfei (XF).
Molecules 23 02907 g001
Figure 2. Scores image for the first seven principal components.
Figure 2. Scores image for the first seven principal components.
Molecules 23 02907 g002
Figure 3. Corresponding optimal wavelengths selected by principal component analysis (PCA): (a) Object-wise analysis. (b) Pixel-wise analysis.
Figure 3. Corresponding optimal wavelengths selected by principal component analysis (PCA): (a) Object-wise analysis. (b) Pixel-wise analysis.
Molecules 23 02907 g003
Figure 4. Classification results using pixel-wise spectra: (a) WHB; (b) XF.
Figure 4. Classification results using pixel-wise spectra: (a) WHB; (b) XF.
Molecules 23 02907 g004
Figure 5. RGB images of the two varieties of raisins: (a) WHB; (b) XF.
Figure 5. RGB images of the two varieties of raisins: (a) WHB; (b) XF.
Molecules 23 02907 g005
Figure 6. Hyperspectral imaging system.
Figure 6. Hyperspectral imaging system.
Molecules 23 02907 g006
Table 1. Corresponding optimal wavelengths selected by PCA.
Table 1. Corresponding optimal wavelengths selected by PCA.
Type of AnalysisNo.Optimal Wavelengths (nm)
Object-wise201005, 1032, 1049, 1086, 1119, 1160, 1173, 1187, 1200, 1220, 1244, 1254, 1278, 1305, 1328, 1352, 1379, 1406, 1433, 1473
Pixel-wise171005, 1029, 1103, 1119, 1164, 1200, 1214, 1251, 1261,
1315, 1328, 1355, 1375, 1406, 1426, 1436, 1473
Table 2. Corresponding optimal wavelengths selected by independent component analysis (ICA).
Table 2. Corresponding optimal wavelengths selected by independent component analysis (ICA).
Type of AnalysisNo.Optimal Wavelengths (nm)
Object-wise20982, 985, 995, 999, 1002, 1009, 1012, 1015, 1019, 1022, 1025, 1029, 1032, 1035, 1039, 1042, 1046, 1049, 1052, 1056
Pixel-wise171139, 1143, 1146, 1150, 1153, 1156, 1207, 1210, 1230, 1521, 1527, 1531, 1548, 1554, 1561, 1575, 1582
Table 3. Classification models based on different grade using optimal wavelengths selected by PCA.
Table 3. Classification models based on different grade using optimal wavelengths selected by PCA.
WHBXFC 4γ 4Cal. ResultPre. Results
WHBXFPre. setWHBXF
Grade1 1Grade113.0665/665245/246Grade31382/13820/602
Grade2930/93122/453
Grade1380/38099/116
Grade2 2Grade225616622/622304/305Grade31371/1382559/602
Grade2305/309146/148
Grade11040/1045323/362
Grade3 3Grade348.59.1950/950405/405Grade3419/432197/197
Grade2658/931434/453
Grade11033/104551/362
1 Grade1 represents large size; 2 Grade2 represents medium size; 3 Grade3 represents small size; 4 C and γ are parameters of SVM model.
Table 4. Classification models based on different grade using optimal wavelengths selected by ICA.
Table 4. Classification models based on different grade using optimal wavelengths selected by ICA.
WHBXFCγCal. ResultPre. Results
WHBXFPre. setWHBXF
Grade1Grade1147.00.3664/665242/246Grade31380/13820/602
Grade2931/93117/453
Grade1379/380100/116
Grade2Grade2147.048.5606/622255/305Grade31360/1382267/602
Grade2296/309119/148
Grade11014/1045306/362
Grade3Grade384.43.0944/950385/405Grade3409/432197/197
Grade2487/931393/453
Grade1899/104515/362
Table 5. Classification results for SVM, k-NN, and RBFNN models based on optimal wavelengths selected by PCA.
Table 5. Classification results for SVM, k-NN, and RBFNN models based on optimal wavelengths selected by PCA.
ModelParameter 5Calibration SetPrediction Set
Acc. 6 (%)Sen. 7Spe. 8Acc. (%)Sen.Spe.
Pixel to pixel 1SVM(256, 5.28)91.830.8980.93980.100.8000.802
k-NN378.480.7000.87078.180.6420.895
RBFNN788.400.8420.92680.890.7970.819
Pixel to object 2SVM(256, 5.28)91.830.8980.93993.620.7850.998
k-NN378.480.7000.87083.820.4640.992
RBFNN788.400.8420.92691.400.7110.997
Object to pixel 3SVM(147, 9.12)99.720.9940.99871.100.8170.626
k-NN595.460.8700.99176.860.7270.803
RBFNN399.780.9940.99954.140.8190.317
Object to object 4SVM(147, 9.12)99.720.9940.99899.120.9870.993
k-NN595.460.8700.99194.060.8390.982
RBFNN399.780.9940.99999.300.9830.997
1 Pixel to pixel means to use models using pixel-wise spectra to predict pixel-wise spectra; 2 Pixel to object means models using pixel-wise spectra to predict object-wise spectra; 3 Object to pixel means to use models using object-wise spectra to predict pixel-wise spectra; 4 Object to object means to use models using object-wise spectra to predict object-wise spectra; 5 Parameters for SVM models are C and γ, parameter for k-NN is number of neighbors (k) and parameter for RBFNN is spread value; 6 Acc. means accuracy; 7 Sen. means sensitivity; 8 Spe. means specificity.
Table 6. Classification results for SVM, k-NN, and RBFNN models based on optimal wavelengths selected by ICA.
Table 6. Classification results for SVM, k-NN, and RBFNN models based on optimal wavelengths selected by ICA.
ModelParameter 5Calibration SetPrediction Set
Acc. 6 (%)Sen. 7Spe. 8Acc. (%)Sen.Spe.
Pixel to pixel 1SVM(256, 16)82.150.7390.90374.90.7080.784
k-NN385.600.7910.89671.130.6180.789
RBFNN678.920.6950.88476.740.7970.819
Pixel to object 2SVM(256, 9.19)82.150.7390.90378.630.2710.998
k-NN385.600.7910.89679.580.3930.962
RBFNN678.920.6950.88480.470.3410.996
Object to pixel 3SVM(147, 84.45)94.680.8790.97654.630.8700.285
k-NN593.640.8490.97462.170.7090.551
RBFNN393.960.8510.97748.340.5650.417
Object to object 4SVM(147, 84.45)94.680.8790.97693.810.8630.969
k-NN593.640.8490.97490.580.8050.947
RBFNN393.960.8510.97793.300.8440.970
1 Pixel to pixel means to use models using pixel-wise spectra to predict pixel-wise spectra; 2 Pixel to object means models using pixel-wise spectra to predict object-wise spectra; 3 Object to pixel means to use models using object-wise spectra to predict pixel-wise spectra; 4 Object to object means to use models using object-wise spectra to predict object-wise spectra; 5 Parameters for SVM models are C and γ, parameter for k-NN is number of neighbors (k) and parameter for RBFNN is spread value; 6 Acc. means accuracy; 7 Sen. means sensitivity; 8 Spe. means specificity.

Share and Cite

MDPI and ACS Style

Feng, L.; Zhu, S.; Zhang, C.; Bao, Y.; Gao, P.; He, Y. Variety Identification of Raisins Using Near-Infrared Hyperspectral Imaging. Molecules 2018, 23, 2907. https://doi.org/10.3390/molecules23112907

AMA Style

Feng L, Zhu S, Zhang C, Bao Y, Gao P, He Y. Variety Identification of Raisins Using Near-Infrared Hyperspectral Imaging. Molecules. 2018; 23(11):2907. https://doi.org/10.3390/molecules23112907

Chicago/Turabian Style

Feng, Lei, Susu Zhu, Chu Zhang, Yidan Bao, Pan Gao, and Yong He. 2018. "Variety Identification of Raisins Using Near-Infrared Hyperspectral Imaging" Molecules 23, no. 11: 2907. https://doi.org/10.3390/molecules23112907

APA Style

Feng, L., Zhu, S., Zhang, C., Bao, Y., Gao, P., & He, Y. (2018). Variety Identification of Raisins Using Near-Infrared Hyperspectral Imaging. Molecules, 23(11), 2907. https://doi.org/10.3390/molecules23112907

Article Metrics

Back to TopTop