On-site variety discrimination of tomato plant using visible-near infrared reflectance spectroscopy

On-site variety discrimination of tomato plant using visible-near infrared reflectance spectroscopy first-author corresp ‡Corresponding Author ‡Corresponding Author †E-mail:ybying@zju.edu.cn †E-mail:ybying@zju.edu.cn ybying@zju.edu.cn College of Biosystems Engineering and Food Science, Zhejiang University, Hangzhou 310029, China College of Biosystems Engineering and Food Science, Zhejiang University, Hangzhou 310029, China The use of visible-near infrared (NIR) spectroscopy was explored as a tool to discriminate two new tomato plant varieties in China (Zheza205 and Zheza207). In this study, 82 top-canopy leaves of Zheza205 and 86 top-canopy leaves of Zheza207 were measured in visible-NIR reflectance mode. Discriminant models were developed using principal component analysis (PCA), discriminant analysis (DA), and discriminant partial least squares (DPLS) regression methods. After outliers detection, the samples were randomly split into two sets, one used as a calibration set ( n =82) and the remaining samples as a validation set ( n =82). When predicting the variety of the samples in validation set, the classification correctness of the DPLS model after optimizing spectral pretreatment was up to 93%. The DPLS model with raw spectra after multiplicative scatter correction and Savitzky-Golay filter smoothing pretreatments had the best satisfactory calibration and prediction abilities (correlation coefficient of calibration ( R c )=0.920, root mean square errors of calibration=0.196, and root mean square errors of prediction=0.216). The results show that visible-NIR spectroscopy might be a suitable alternative tool to discriminate tomato plant varieties on-site. Visible-NIR spectroscopy Tomato plant variety Discrimination Principal component analysis (PCA) Discriminant analysis (DA) Discriminant partial least squares (DPLS) INTRODUCTION 1 The development of agricultural science and technology is promoting a large range of commercial varieties in many important agricultural species. The identification of varieties of agricultural and horticultural crops is important during all stages of seed production, trading, and agricultural production and processing (Cooke, 1995b 1995b ). The expansion of the international trade creates increasing interest in the descriptive characterization of plant varieties in the context of intellectual property protection, stimulated by the recent agreements within the framework of the World Trade Organization (Cooke, 1999 1999 ). Traditionally, variety identification has been carried out using morphological characters. However, due to modern breeding technology, lack of phenotypic variation causes difficulties for variety identification by the conventional morphological methods. Moreover, modern crop production increasingly demands methods that are more rapid and discriminating. In recent years, a great deal of research has been done on various techniques for realizing a robust identification method, such as molecular markers (Cooke, 1995a 1995a ; Lee et al., 2005 2005 ; de Riek et al., 2007 2007 ) and gene expression profiling (Ancillo et al., 2007 2007 ), but these methods cannot perform on-site measuring. In addition, they are also expensive and time-consuming. In previous studies, it was found that near infrared (NIR) spectroscopy has the potential to determine plant pathological conditions (Riedell and Blackmer; 1999 1999 ; Xu et al., 2007 2007 ) and to authenticate species based on samples’ spectral properties (Chen et al., 1998 1998 ). Kim et al.( 2000 2000 ) used linear discriminant analysis (LDA) and non-linear techniques based on multilayer perceptrons (MLPs) with variations on back-propagation (BP) learning to classify kiwifruit grown under different conditions, using a range of features extracted from visible-NIR spectra. Results indicate that the performance of non-linear models was significantly better in comparison to the linear model, LDA. Noble and Crowe ( 2001 2001 ) developed classifiers using five features extracted from diffuse reflectance spectra in the region of 250~2500 nm to distinguish six plant species. The accuracies were greater than 90%. Perez et al.( 2001 2001 ) used NIR reflectance spectroscopy to authenticate two green asparagus varieties using modified partial least square regression. In their results, the standard error of cross-validation and coefficient of determination were 0.082 and 0.97, respectively. Fu et al.( 2007 2007 ) explored the use of Fourier transform NIR spectrometry to discriminate three pear varieties, applying discriminant partial least squares (DPLS), discriminant analysis (DA) and probabilistic neural networks calibration models developed in the region of 800~2500 nm, and gained a high classification accuracy of 97.5%~100%. Wang and Paliwal ( 2006 2006 ) used NIR spectroscopy to differentiate six classes of wheat by five classifiers. Classification was performed on principal component score (PCS), Fourier coefficients, and wavelet coefficients of the original spectral data. LDA combined with principal component analysis (PCA) as a pre-processing method had the best classification accuracies. Yee et al.( 2006 2006 ) used NIR spectra to identify the potato tuber variety based on LDA pattern recognition techniques. The success rate in separating the spectra into respective classes was 93%. Chen et al.( 2007 2007 ) studied the differentiation of three tea categories using support vector machine (SVM). The top five latent variables extracted by PCA were used as the input of SVM classifiers. The identification rates of the three teas were up to 90%, 100%, and 93.3%. The objective of this study was to apply visible-NIR reflectance spectroscopy for discrimination of two new commercial tomato plant varieties in China (Zheza205 and Zheza207). The specific goals were: (1) to investigate the potential of visible-NIR reflectance spectroscopy in on-site plant varieties discrimination; (2) to compare the modeling performance of the three discriminant models based on PCA, DA, and DPLS. MATERIALS AND METHODS 1 Experimental materials 2 Two new varieties of tomato plants (Zheza205 and Zheza207, bred by Zhejiang Academy of Agricultural Science, China) were cultivated in an environmentally controlled artificial climate chamber at Zhejiang University. Leaf-level reflectance measurements were made at the mid growth stage of tomato plant (about 20 cm height). We selected 82 top-canopy leaves of Zheza205 and 86 top-canopy leaves of Zheza207 to perform on-site measuring on the same day (at 30.4 °C with relatively humidity of 57.8%). NIR spectra acquisition 2 The visible-NIR reflectance spectra were collected by a commercial USB2000 spectrometer (Ocean Optics Inc., USA). The USB2000 spectrometer was equipped with a grating (600 lines blazed at 500 nm), a detector (350~1100 nm) with a resolution of 1.5 nm, and an optical fiber cable. The optical fiber was placed 5 mm away from the leaf and directed to the center of the illuminated area at 45° to the light fiber to collect the dispersed light from the leaf as described elsewhere (Cervantes-Martínez et al., 2002 2002 ). The spectrum of each leaf sample was the average of 5 successive scans. Spectral data acquisition and manipulation were performed in real time with a microcomputer and commercially available software OOIBase32 (Ocean Optics Inc., USA). In this study, only the spectral information between 400~900 nm (Fig. 1 1 ) was used for analyzing because the signal-to-noise ratio was very low beyond 900 nm. Data processing and discriminant analysis 2 PCA, DA, and DPLS methods were performed using TQ Analyst v6.2.1 software (Thermo Electron Corp., USA). After spectral outliers removal using a standardized H value (Workman, 2001 2001 ; Steuer et al., 2001 2001 ), the remaining 164 samples’ spectra were analyzed using PCA and DA on the raw spectra with multiplicative scatter correction (MSC) pretreatment. For DPLS discriminant analysis, the remaining 164 samples were divided into two sets randomly, namely, a calibration set and a validation set, with each comprising 82 top-canopy leaves of Zheza205 or 82 top-canopy leaves of Zheza207. The derivative treatments “0, 5, 3” and “1, 5, 3” were used to reduce baseline variation and spectral smoothing after MSC pretreatment, respectively, where the 1st digit is the number of the derivative, the 2nd one is the length of the Savitzky-Golay filter, and the 3rd one is the polynomial order of the Savitzky-Golay filter (Perez et al., 2001 2001 ). PCA is a method to reduce the dimensionality of a dataset composed of a large number of variables by transforming the original dataset into a new set of variables, the PCS, which are orthogonal and uncorrelated to each other and are ordered so that the first few PCS retain most of the variation present in all of the original variables (Cozzolino et al., 2003 2003 ; Lee et al., 2003 2003 ). PCA was used in further analysis to examine the interrelationships between different variables, and to interpret sample patterns, groupings, similarities, or differences (Cozzolino et al., 2003 2003 ; Cozzolino and Murray, 2004 2004 ; Mouazen et al., 2006 2006 ; Xie et al., 2007b 2007b ). DA is a qualitative tool used to group (or classify) spectra into pre-defined classes. In this study, DA was used to classify the two varieties of tomato plants, Zheza205 and Zheza207, by computing a sample’s distance from each class center in Mahalanobis distance (MD) units. The MD is the parameter that is calculated to determine how close to the center of gravity of its group an individual sample spectrum is. The closer MD is to zero, the better is the match. As such, the smallest MD is used to pick the group that the individual fits best. In DPLS regression technique, analysis was performed using the dummy regression technique as described elsewhere (Perez et al., 2001 2001 ; Cozzolino et al., 2003 2003 ; Xie et al., 2007b 2007b ). Each sample in the calibration and validation sets is assigned a dummy variable (1 or 2) as a reference value, an arbitrary number indicating whether the sample belongs to a particular group or not (Xie et al., 2007b 2007b ). In this case, samples of Zheza205 were assigned a numeric value of 1 and samples of Zheza207 a numeric value of 2. The DPLS model was developed by regression of the spectral data (the first derivative) against the assigned reference value (dummy variable), and then was tested for accuracy by using it to predict the variety of samples in the validation set. A sample was classified as Zheza205 if its predicted value was between 0.5 and 1.5, and it was classified as Zheza207 if the value was between 1.5 and 2.5 as described by Cozzolino et al.( 2003 2003 ) and Xie et al.( 2007b 2007b ). The discriminant accuracy was described by the correlation coefficient of calibration ( R c ) and the root mean square errors of prediction ( RMSEP ). The DPLS regression equation can be simply formulated as described by Alsberg et al.( 1998 1998 ): Y = XB , where X is the data matrix created with the spectra of samples. Each row of the matrix X contains the spectrum of a sample, and each column contains the intensity at each wavelength; the matrix Y constructed with dummy variables contains information about class memberships of samples; the matrix B represents the estimated regression coefficients (Wang et al., 1998 1998 ). And the RMSEP is done by $$RMSEP = \sqrt {\frac{1}{{N_{\text{p}} }}\sum\limits_{i = 1}^{N_{\text{p}} } {(\hat y_i - y_i )^2 } }$$ , where y ˆ i represents the predicted value of the i th observation, y i represents the measured value of the i th observation, and N p is the number of observations in the prediction set. Microscopy analysis 2 Tomato leaf tissues (1 mm×3 mm) of Zheza205 and Zheza207 were fixed with 2.5% (v/v) glutaraldehyde overnight at 4 °C and subsequently with 1% (w/v) osmium acid for 1~2 h. After fixation, the tissues were washed three times in phosphate buffer (pH=7.0) for 15 min each, and dehydrated in a graded ethanol series (50%, 70%, 80%, 90%, 95%, and 100%, v/v) for 15~20 min each and in 100% acetone for 20 min. The tissues were embedded in 1:1 (v/v) mixture of the embedding medium (Spurr epoxy resin) and 100% acetone for 1 h at room temperature, then in 3:1 (v/v) mixture of the embedding medium and 100% acetone for 3 h, and in 100% embedding medium overnight. Finally, the tissues were placed in capsules containing embedding medium at 70 °C overnight. The embedded material was sectioned into slices (70~90 nm) with a diamond knife on an ultra microtome (Reichert Ultracut E, Germany) and was stained by 50% (v/v) uranyl acetate and lead citrate for 15 min, respectively. Stained samples were observed with a JEM-1230 transmission electron microscope (JEOL Co., Ltd., Japan). RESULTS 1 PCA was performed on the raw spectra to examine qualitative differences between the two tomato plant varieties. The first three PCS (95.422% of the variation explained) of the visible-NIR data were plotted in Fig. 2 2 . PC1 explained 66.342% of the total variation in the samples, PC2 26.859%, and PC3 2.222%. The three-dimension score plot showed two major clusters of samples related with the two different tomato plant varieties. This confirmed the assumptions that discrimination between varieties was possible and that different spectral attributes between samples were associated with characteristics of the tomato plants depending on the varieties (Cozzolino and Murray, 2004 2004 ; Xie et al., 2007b 2007b ). As PCA only indicates the visualizing dimension spaces, DA and DPLS were utilized for an improved discrimination as described elsewhere (Xie et al., 2007b 2007b ). A DA model was developed on the first eight PCS resulting from the PCA as they covered most of the variation (99.9% of variation described) contained in the raw spectra. This model was used to classify the two tomato plant varieties. Two Zheza207 samples were misclassified as Zheza205; i.e., the MD to the center of the gravity of Zheza207 group is superior to that to the center of the gravity of Zheza205 group, and 100% of Zheza205 samples were classified correctly by the DA discriminant model (Fig. 3 3 ). The calibration and validation statistics and performances of the dummy classification models using PLS regression method are shown in Table 1 1 . All the models developed had a coefficient of determination in calibration ( R c ) higher than 0.80 and most of the samples in calibration and prediction sets were classified correctly by the classification models ( > 93%). From Table 1 1 , it is clear that models using the 1st derivative data led to high RMSEP values and relatively low root mean square errors of calibration ( RMSEC ) values compared with raw spectral data. These results generally could be due to overfitting of the calibration. Similar discrimination results were reported previously by Wang et al.( 2006 2006 ). The discrimination model by the “1, 5, 3” derivative treatment with three factors resulted in RMSEC =0.184 and RMSEP =0.239, which were better than those of the same derivative treatment with only one factor. Factor increase led to reduction of the RMSEC and improved the classification accuracy in calibration, but not as significant as in prediction. Overall, the results show that the model producing the best prediction is the one that combines MSC and Savitzky-Golay filter smoothing pretreatments in the range of 400~900 nm using raw visible-NIR spectra with three factors. MSC: multiplicative scatter correction; R c : coefficient of determination in calibration; RMSEC : root mean square errors of calibration; RMSEP : root mean square errors of prediction; n c : number classified correctly; P c : percentage classified correctly 2 No. 2 Data pretreatment 2 Derivative treatment 2 Factors 4 Calibration 3 Prediction R c RMSEC n c P c (%) RMSEP n c P c (%) 1 MSC “0, 5, 3” 3 0.920 0.196 81 98.78 0.216 80 97.56 2 MSC “1, 5, 3” 1 0.800 0.300 77 93.90 0.278 77 93.90 3 MSC “1, 5, 3” 3 0.930 0.184 82 100 0.239 77 93.90 Fig. 4 4 shows the visible-NIR predictions of the two tomato plant varieties in the calibration set ( n =164) using the best DPLS model by the “0, 5, 3” derivative treatment. Fig. 4 4 clearly shows that two Zheza207 samples were misclassified as Zheza205 (i.e., the calculated values below 1.5) but 100% of Zheza205 samples were classified correctly. Fig. 5 5 shows the transmission electron microscope (TEM) images of the sections of Zheza205 and Zheza207 tomato plant leaf tissues. As shown in the TEM images, there is a transparent waxy cuticle layer at the top of the Zheza207 tomato leaf (Fig. 5b 5b ), but no waxy cuticle layer at the top of Zheza205 tomato leaf (Fig. 5a 5a ). The cuticular wax plays a critical role in the interaction of light with the leaves (Holmes and Keiller, 2002 2002 ), and is known to increase the reflectance of the leaves (Vogelmann, 1993 1993 ). We also see from Fig. 1 1 that reflectance was significantly higher from the Zheza207 tomato leaf with a waxy cuticle layer than from the Zheza205 tomato leaf. DISCUSSION 1 The ability of visible-NIR model to discriminate or identify plant varieties is based on the leaf spectral properties that are related to the leaf biochemical composition and structure, which depend on many factors such as the plant species, the development or microclimate position of the leaf on the plant (Jacquemoud and Ustin, 2001 2001 ). The visible-NIR optical domain ranges can be divided into two parts: in the visible spectra (400~750 nm), leaf reflectance is low because of absorption by photosynthetic pigments (mainly chlorophylls and carotenoids); in the short-wavelength near infrared (SW-NIR) domain (750~1300 nm), where there are no strong absorption features, and the magnitude of reflectance is governed by structural discontinuities encountered in the leaf. Leaf structure difference and physiological and ecological changes on plants caused by the influence of man or nature can affect the optical properties (Borregaard et al., 2000 2000 ). Several studies on plant variety discrimination based on their optical properties have been carried out, such as authentication of green asparagus varieties (Perez et al., 2001 2001 ), discrimination of different classes of wheat (Wang and Paliwal, 2006 2006 ), identification of green, black and Oolong teas (Chen et al., 2007 2007 ), and classification of nontransgenic and transgenic tomatoes (Xie et al., 2007a 2007a ; 2007c 2007c ). These studies have proven that visible-NIR spectroscopy is an accurate technology for identification and authentication of plant varieties. CONCLUSION 1 This study shows that visible-NIR spectroscopy, together with chemometric techniques such as PCA, DA, and DPLS, could be used for the identification of plant varieties on-site. Our results indicate that an excellent classification can be obtained by DPLS method after optimizing spectral data pretreatment, with accuracy up to 93%. Further studies are needed to improve the calibration accuracy and robustness with more samples and to extend the discrimination to other plant varieties. References BK Alsberg , DB Kell , R Goodacre . Variable selection in discriminant partial least-squares analysis . Anal Chem , 1998 . 70 ( 19 ): 4126 - 4133 . DOI: 10.1021/ac980506o http://doi.org/10.1021/ac980506o . Variable selection in discriminant partial least-squares analysis G Ancillo , J Gadea , J Forment , J Guerri , L Navarro . Class prediction of closely related plant varieties using gene expression profiling . J Exp Bot , 2007 . 58 ( 8 ): 1927 - 1933 . 17452756 DOI: 10.1093/jxb/erm054 http://doi.org/10.1093/jxb/erm054 . Class prediction of closely related plant varieties using gene expression profiling T Borregaard , H Nielsen , L Nørgaard , H Have . Crop-weed discrimination by line imaging spectroscopy . J Agric Eng Res , 2000 . 75 ( 4 ): 389 - 400 . DOI: 10.1006/jaer.1999.0519 http://doi.org/10.1006/jaer.1999.0519 . Crop-weed discrimination by line imaging spectroscopy J Cervantes-Martínez , R Flores-Hernández , B Rodríguez-Garay , F Santacruz-Ruvalcaba . Detection of bacterial infection of agave plants by laser-induced fluorescence . Appl Optics , 2002 . 41 ( 13 ): 2541 - 2545 . DOI: 10.1364/AO.41.002541 http://doi.org/10.1364/AO.41.002541 . Detection of bacterial infection of agave plants by laser-induced fluorescence LM Chen , NC Carpita , WD Reiter , RH Wilson , C Jeffries , MC McCann . A rapid method to screen for cell-wall mutants using discriminant analysis of Fourier transform infrared spectra . Plant J , 1998 . 16 ( 3 ): 385 - 392 . 9881159 DOI: 10.1046/j.1365-313x.1998.00301.x http://doi.org/10.1046/j.1365-313x.1998.00301.x . A rapid method to screen for cell-wall mutants using discriminant analysis of Fourier transform infrared spectra QS Chen , JW Zhao , CH Fang , DM Wang . Feasibility study on identification of green, black and oolong teas using near infrared reflectance spectroscopy based on support vector machine . Spectrochim Acta A , 2007 . 66 ( 3 ): 568 - 574 . DOI: 10.1016/j.saa.2006.03.038 http://doi.org/10.1016/j.saa.2006.03.038 . Feasibility study on identification of green, black and oolong teas using near infrared reflectance spectroscopy based on support vector machine RJ Cooke . Gel electrophoresis for the identification of plant varieties . J Chromatogr A , 1995 . 698 ( 1-2 ): 281 - 299 . DOI: 10.1016/0021-9673(94)00649-T http://doi.org/10.1016/0021-9673(94)00649-T . Gel electrophoresis for the identification of plant varieties RJ Cooke . CW Wrigley . Introduction: The Reasons for Variety Identification . Identification of Food Grain Varieties , 1995 . : USA AACC , 1 - 18 . . Introduction: The Reasons for Variety Identification RJ Cooke . New approaches to potato variety identification . Potato Res , 1999 . 42 ( 3-4 ): 529 - 539 . DOI: 10.1007/BF02358169 http://doi.org/10.1007/BF02358169 . New approaches to potato variety identification D Cozzolino , I Murray . Identification of animal meat muscles by visible and near infrared reflectance spectroscopy . Lebensm-Wiss u-Technol , 2004 . 37 ( 4 ): 447 - 452 . DOI: 10.1016/j.lwt.2003.10.013 http://doi.org/10.1016/j.lwt.2003.10.013 . Identification of animal meat muscles by visible and near infrared reflectance spectroscopy D Cozzolino , HE Smyth , M Gishen . Feasibility study on the use of visible and near infrared spectroscopy together with chemometrics to discriminate between commercial white wines of different varietal origins . J Agric Food Chem , 2003 . 51 ( 26 ): 7703 - 7708 . 14664532 DOI: 10.1021/jf034959s http://doi.org/10.1021/jf034959s . Feasibility study on the use of visible and near infrared spectroscopy together with chemometrics to discriminate between commercial white wines of different varietal origins J de Riek , I Everaert , D Esselink , E Calsyn , MJM Smulders , B Vosman . Assignment tests for variety identification compared to genetic similarity-based methods using experimental datasets from different marker systems in sugar beet . Crop Sci , 2007 . 47 ( 5 ): 1964 - 1974 . DOI: 10.2135/cropsci2006.04.0265 http://doi.org/10.2135/cropsci2006.04.0265 . Assignment tests for variety identification compared to genetic similarity-based methods using experimental datasets from different marker systems in sugar beet XP Fu , Y Zhou , YB Ying , HS Lu , HR Xu . Discrimination of pear varieties using three classification methods based on near infrared spectroscopy . Trans ASABE , 2007 . 50 ( 4 ): 1355 - 1361 . . Discrimination of pear varieties using three classification methods based on near infrared spectroscopy MG Holmes , DR Keiller . Effects of pubescence and waxes on the reflectance of leaves in the ultraviolet and photosynthetic wavebands: a comparison of a range of species . Plant Cell Environ , 2002 . 25 ( 1 ): 85 - 93 . DOI: 10.1046/j.1365-3040.2002.00779.x http://doi.org/10.1046/j.1365-3040.2002.00779.x . Effects of pubescence and waxes on the reflectance of leaves in the ultraviolet and photosynthetic wavebands: a comparison of a range of species S Jacquemoud , SL Ustin . Leaf Optical Properties: A State of the Art . Proceedings of the 8th International Symposium on Physical Measurements and Signatures in Remote Sensing , 2001 . : Aussois, France 223 - 232 . . Leaf Optical Properties: A State of the Art J Kim , A Mowat , P Poole , N Kasabov . Linear and non-linear pattern recognition models for classification of fruit from visible-near infrared spectra . Chemometr Intell Lab Syst , 2000 . 51 ( 2 ): 201 - 216 . DOI: 10.1016/S0169-7439(00)00070-8 http://doi.org/10.1016/S0169-7439(00)00070-8 . Linear and non-linear pattern recognition models for classification of fruit from visible-near infrared spectra D Lee , JC Reeves , RJ Cooke . DNA profiling and plant variety registration: 1. the use of random amplified DNA polymorphisms to discriminate between varieties of oilseed rape . Electrophoresis , 2005 . 17 ( 1 ): 261 - 265 . DOI: 10.1002/elps.1150170145 http://doi.org/10.1002/elps.1150170145 . DNA profiling and plant variety registration: 1. the use of random amplified DNA polymorphisms to discriminate between varieties of oilseed rape KH Lee , NQ Zhang , S Das . Comparing Adaptive Neuro-fuzzy Inference System (ANFIS) to Partial Least-squares (PLS) Method for Simultaneous Prediction of Multiple Soil Properties , 2003 . : Las Vegas, Nevada, USA ASAE Paper No. 033144 . Comparing Adaptive Neuro-fuzzy Inference System (ANFIS) to Partial Least-squares (PLS) Method for Simultaneous Prediction of Multiple Soil Properties AM Mouazen , R Karoui , JD Baerdemaeker , 等 . Classification of Soils into Different Moisture Content Levels Based on VIS-NIR Spectra , 2006 . : Portland, USA ASABE Paper No. 061067 . Classification of Soils into Different Moisture Content Levels Based on VIS-NIR Spectra SD Noble , TG Crowe . Plant Discrimination Based on Leaf Reflectance , 2001 . : Sacramento, California, USA ASABE Paper No. 011150 . Plant Discrimination Based on Leaf Reflectance DP Perez , MT Sanchez , G Cano , A Garrido . Authentication of green asparagus varieties by near infrared reflectance spectroscopy . J Food Sci , 2001 . 66 ( 2 ): 323 - 327 . DOI: 10.1111/j.1365-2621.2001.tb11340.x http://doi.org/10.1111/j.1365-2621.2001.tb11340.x . Authentication of green asparagus varieties by near infrared reflectance spectroscopy WE Riedell , TM Blackmer . Leaf reflectance spectra of cereal aphid-damaged wheat . Crop Sci , 1999 . 39 ( 6 ): 1835 - 1840 . . Leaf reflectance spectra of cereal aphid-damaged wheat B Steuer , H Schulz , E Lager . Classification and analysis of citrus oils by NIR spectroscopy . Food Chem , 2001 . 72 ( 1 ): 113 - 117 . DOI: 10.1016/S0308-8146(00)00209-0 http://doi.org/10.1016/S0308-8146(00)00209-0 . Classification and analysis of citrus oils by NIR spectroscopy TC Vogelmann . Plant tissue optics . Annu Rev Plant Physiol Plant Mol Biol , 1993 . 44 ( 1 ): 231 - 251 . DOI: 10.1146/annurev.pp.44.060193.001311 http://doi.org/10.1146/annurev.pp.44.060193.001311 . Plant tissue optics CY Wang , CT Chen , CP Chiang , ST Young , SN Chow , HK Chiang . Partial least-squares discriminant analysis on autofluorescence spectra of oral carcinogenesis . Appl Spectrosc , 1998 . 52 ( 9 ): 1190 - 1195 . DOI: 10.1366/0003702981945002 http://doi.org/10.1366/0003702981945002 . Partial least-squares discriminant analysis on autofluorescence spectra of oral carcinogenesis L Wang , SC Frank Lee , XR Wang , Y He . Feasibility study of quantifying and discriminating soybean oil adulteration in camellia oils by attenuated total reflectance MIR and fiber optic diffuse reflectance NIR . Food Chem , 2006 . 95 ( 3 ): 529 - 536 . DOI: 10.1016/j.foodchem.2005.04.015 http://doi.org/10.1016/j.foodchem.2005.04.015 . Feasibility study of quantifying and discriminating soybean oil adulteration in camellia oils by attenuated total reflectance MIR and fiber optic diffuse reflectance NIR W Wang , J Paliwal . Spectral data compression and analyses techniques to discriminate wheat classes . Trans ASABE , 2006 . 49 ( 5 ): 1607 - 1612 . . Spectral data compression and analyses techniques to discriminate wheat classes J Workman . DA Burns , EW Ciurczak . NIR Spectroscopy Calibration Basics . Handbook of near Infrared Analysis , 2001 . 2nd Ed : New York Marcfl Dfkker, Inc , 90 - 129 . . NIR Spectroscopy Calibration Basics LJ Xie , YB Ying , TJ Ying . Combination and comparison of chemometrics methods for identification of transgenic tomatoes using visible and near-infrared diffuse transmittance technique . J Food Eng , 2007 . 82 ( 3 ): 395 - 401 . DOI: 10.1016/j.jfoodeng.2007.02.062 http://doi.org/10.1016/j.jfoodeng.2007.02.062 . Combination and comparison of chemometrics methods for identification of transgenic tomatoes using visible and near-infrared diffuse transmittance technique LJ Xie , YB Ying , TJ Ying , HY Yu . Discrimination of transgenic tomatoes based on visible/near-infrared spectra . Anal Chim Acta , 2007 . 584 ( 2 ): 379 - 384 . 17386628 DOI: 10.1016/j.aca.2006.11.071 http://doi.org/10.1016/j.aca.2006.11.071 . Discrimination of transgenic tomatoes based on visible/near-infrared spectra LJ Xie , YB Ying , TJ Ying . Quantification of chlorophyll content and classification of nontransgenic and transgenic tomato leaves using visible/near-infrared diffuse reflectance spectroscopy . J Agric Food Chem , 2007 . 55 ( 12 ): 4645 - 4650 . 17503831 DOI: 10.1021/jf063664m http://doi.org/10.1021/jf063664m . Quantification of chlorophyll content and classification of nontransgenic and transgenic tomato leaves using visible/near-infrared diffuse reflectance spectroscopy HR Xu , YB Ying , XP Fu , SP Zhu . Near infrared spectroscopy in detecting leaf miner damage on tomato leaf . Biosyst Eng , 2007 . 96 ( 4 ): 447 - 454 . DOI: 10.1016/j.biosystemseng.2007.01.008 http://doi.org/10.1016/j.biosystemseng.2007.01.008 . Near infrared spectroscopy in detecting leaf miner damage on tomato leaf N Yee , WT Bussell , GG Coghill . Use of near infrared spectra to identify cultivar in potato ( Solanum tuberosum ) crisps . New Zeal J Crop Hort , 2006 . 34 ( 2 ): 177 - 181 . . Use of near infrared spectra to identify cultivar in potato (Solanum tuberosum) crisps 10.1631/jzus.B0820200 Project (No. 60405003) supported by the National Natural Science Foundation of China 2008-06-24 2009-01-07 2008-12-10 2022-01-05 Copyright © 2009, Journal of Zhejiang University Science 2009 Journal of Zhejiang University-SCIENCE B (Biomedicine and Biotechnology) 02 10