Nothing Special   »   [go: up one dir, main page]
Next Article in Journal
Design of a New Acoustic Logging While Drilling Tool
Previous Article in Journal
Combining Inertial Sensors and Machine Learning to Predict vGRF and Knee Biomechanics during a Double Limb Jump Landing Task
You seem to have javascript disabled. Please note that many of the page functionalities won't work as expected without javascript enabled.
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Sensors
Volume 21
Issue 13
10.3390/s21134384
sensors-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Rice Seed Purity Identification Technology Using Hyperspectral Image with LASSO Logistic Regression Model

by
Weihua Liu
1,
Shan Zeng
2,*,
Guiju Wu
3,
Hao Li
2 and
Feifei Chen
2
1
School of Electric & Electronic Engineering, Wuhan Polytechnic University, Wuhan 430023, China
2
School of Mathematics & Computer Science, Wuhan Polytechnic University, Wuhan 430023, China
3
The Key Laboratory of Earthquake Geodesy, Institute of Seismology, China Earthquake Administration, Wuhan 430023, China
*
Author to whom correspondence should be addressed.
Sensors 2021, 21(13), 4384; https://doi.org/10.3390/s21134384
Submission received: 29 May 2021 / Revised: 22 June 2021 / Accepted: 24 June 2021 / Published: 26 June 2021
(This article belongs to the Special Issue Research and Application of Robust Hyperspectral Image)
Browse Figures
Figure 1
<p>Schematic of the hyperspectral imaging system for rice seeds.</p> "> Figure 2
<p>Image of four rice seed samples through the viewfinder of the hyperspectral camera; the first (top) to fourth (bottom) rows are Huguangxiang, Xiangwan Japonica, Huanghuazhan, and Japonica 530, respectively.</p> "> Figure 3
<p>Flowchart of hyperspectral image processing, including data pre-processing, background segmentation, data set preparation, modelling, and identification accuracy output.</p> "> Figure 4
<p>Example of background segmentation using synthesised spectral slope greyscale image. (<b>a</b>) Colour image of a seed; (<b>b</b>) Binarised image using single-band (601.55 nm) information; (<b>c</b>) Segmentation result using single-band (601.55 nm) information; (<b>d</b>) Synthesised spectral slope greyscale image; (<b>e</b>) Binarised image using spectral slope information; (<b>f</b>) Segmentation result using spectral slope information.</p> "> Figure 5
<p>Average spectral features of four rice seeds. (<b>a</b>) Original reflectance spectrum; (<b>b</b>) spectrum after SNV; (<b>c</b>) spectrum after FD’ (<b>d</b>) spectrum after SD.</p> "> Figure 5 Cont.
<p>Average spectral features of four rice seeds. (<b>a</b>) Original reflectance spectrum; (<b>b</b>) spectrum after SNV; (<b>c</b>) spectrum after FD’ (<b>d</b>) spectrum after SD.</p> "> Figure 6
<p>Change in parameters with different <math display="inline"><semantics> <mi>λ</mi> </semantics></math> values. (<b>a</b>) Regression coefficient; (<b>b</b>) fitting deviation; (<b>c</b>) number of feature wavelength bands; (<b>d</b>) prediction accuracy.</p> "> Figure 7
<p>Wavelength band selection using SNV data for different rice varieties under the optimal λ. (<b>a</b>) Huguangxiang; (<b>b</b>) Xiangwan Japonica; (<b>c</b>) Huanghuazhan; (<b>d</b>) Japonica 530.</p> "> Figure 8
<p>Wavelength band selection using FD data for different rice varieties under the optimal λ. (<b>a</b>) Huguangxiang. (<b>b</b>) Xiangwan Japonica. (<b>c</b>) Huanghuazhan. (<b>d</b>) Japonica 530.</p> "> Figure 9
<p>Wavelength band selection using SD data for different rice varieties under the optimal <math display="inline"><semantics> <mi>λ</mi> </semantics></math>. (<b>a</b>) Huguangxiang; (<b>b</b>) Xiangwan Japonica; (<b>c</b>) Huanghuazhan; (<b>d</b>) Japonica 530.</p> ">
Review Reports Versions Notes

Abstract

:
Hyperspectral technology is used to obtain spectral and spatial information of samples simultaneously and demonstrates significant potential for use in seed purity identification. However, it has certain limitations, such as high acquisition cost and massive redundant information. This study integrates the advantages of the sparse feature of the least absolute shrinkage and selection operator (LASSO) algorithm and the classification feature of the logistic regression model (LRM). We propose a hyperspectral rice seed purity identification method based on the LASSO logistic regression model (LLRM). The feasibility of using LLRM for the selection of feature wavelength bands and seed purity identification are discussed using four types of rice seeds as research objects. The results of 13 different adulteration cases revealed that the value of the regularisation parameter was different in each case. The recognition accuracy of LLRM and average recognition accuracy were 91.67–100% and 98.47%, respectively. Furthermore, the recognition accuracy of full-band LRM was 71.60–100%. However, the average recognition accuracy was merely 89.63%. These results indicate that LLRM can select the feature wavelength bands stably and improve the recognition accuracy of rice seeds, demonstrating the feasibility of developing a hyperspectral technology with LLRM for seed purity identification.

1. Introduction

Paddy is one of the most important food sources worldwide, and its large number of varieties is continuing to increase. The phenomenon of seed mixing exists inevitably in each stage of paddy planting, production, and circulation, which affects the yield, nutrient composition, taste, and price of paddy rice. Therefore, the identification of rice seed purity is crucial in the grain industry. However, the identification task is primarily conducted manually with the naked eye [1] and can thus be easily affected by subjective factors such as experience and tiredness. Certain chemical or biological methods, such as phenol staining and protein electrophoresis, can accurately identify seed purity [2,3]; however, these methods are destructive, time-consuming, require large equipment or expertise, and cannot be batch processed. Therefore, it is vital to develop a rapid and accurate automatic identification method for rice seed purity.
Recently, nondestructive testing methods based on optical technology, such as infrared spectroscopy technology [4,5], machine vision, and hyperspectral technology [6,7,8,9,10,11,12,13], have received significant attention. Infrared spectroscopy technology identifies the spectral features of seeds [14] and analyses their internal chemical compositions (such as protein and starch) to classify them. However, this technology can only acquire a small area of spectral information from the seed samples at a time, and many repetitive measurements in different places are needed to represent the entire scenario of adulteration [15]. Machine vision technology captures colour RGB images of seeds and analyses their external physical features to classify them. Several reports on relevant machine vision technology have been published [16,17]. For example, Chaugule et al. (2014) classified four types of rice seeds based on their colour, shape, and texture with a recognition accuracy of 90.30–94.00% [18]. Additionally, they used the angle of seeds as a new identification feature in 2016, improving the identification accuracy to 97.6% [19]. In 2017, Huang et al. analysed the key details of both ends of rice seeds and found that the identification accuracies of three rice seeds with similar appearances were 92.68%, 97.35%, and 96.57%, respectively [20]. Furthermore, in 2021, Ansari et al. adopted machine vision technology combined with a multivariate analysis method to establish a rapid detection method of rice seeds based on different purities, obtaining the highest accuracy of 93.9% [21]. However, the phenomena of overlapping and superposition of external physical features are often observed owing to the similar appearances of many rice seeds. Therefore, it is difficult to obtain a stable classification effect using spatial information alone [22].
Hyperspectral technology can obtain the spectral and spatial information of seed samples, and it has attracted significant attention in the field of seed identification [23,24,25]. There have been many reports on the use of hyperspectral technology to identify corn, wheat, soybean, rice, and other grain particles. For example, in 2013, Silvia et al. used partial least squares discriminant analysis to identify hyperspectral images of groats and oats with an accuracy of approximately 100% [26]. In 2015, Yang et al. used a support vector machine (SVM) to identify waxy corn seeds, achieving the highest identification accuracy of 98.2% [27]. In 2020, Zhu et al. used a convolutional neural network to identify hyperspectral images of six types of soybean seeds with an average recognition rate of 91% [28]. The identification methods of rice can be broadly classified into two categories. The first incorporates the popular deep learning algorithms for classification and recognition. For example, in 2018, Qiu et al. used a convolutional neural network to identify four types of rice; they achieved a recognition rate of 87% using 12,000 training samples [29]. Furthermore, in 2020, Weng et al. used a deep learning network, namely principal component analysis (PCA) network, to train 4320 samples, obtaining the best recognition rate of over 98% for 10 rice varieties [30]. However, although the deep learning method does not require feature extraction, it requires a large number of training samples, leading to relatively high costs of hyperspectral image collection and data processing. The second category uses the traditional machine learning methods for identification. The most effective method is to use PCA dimension reduction and SVM classification. In 2015, Sun et al. used a combination of spectral, image, and textural features to identify four types of rice using SVM with a recognition accuracy of 91.6% [31]. Furthermore, in 2020, Liu et al. used a combination of spectral, image, and textural features to identify four types of rice seeds with an average identification accuracy of 98% [32]. These methods only require a small number of training samples to achieve a satisfactory recognition rate. However, hyperspectral spatial information contains a significant amount of redundant information and unpredictable noise, leading to high costs of hyperspectral image acquisition and processing. The cost of data acquisition and processing can be further reduced if the feature wavelength bands can be selected via some mathematical methods such as PCA and analysis of variance [33].
There are several methods to reduce the dimension of hyperspectral data [34,35], such as correlation analysis [36], beta coefficient of partial least squares regression analyses [37], and PCA [38]. PCA is widely used in the field of grain recognition [39,40,41,42]. However, it is essentially a dimensionality reduction method. Feature wavelength bands selected by PCA methods are not suitable for physical interpretation; it is easy to ignore those wavelength bands with a small contribution but contain more useful information for classification. Moreover, feature band selection and classification are two independent algorithms. However, there exists a complex interaction between feature band selection and classification modelling. The algorithm featuring simultaneous selection and classification is rarely used for the identification of rice seed purity.
In 1996, Tibshirani et al. proposed LASSO regression, a data sparse method [43] often used for hyperspectral data unmixing [44,45,46]. For the first time in 2017, Yang et al. introduced the LASSO method into hyperspectral images for band selection by conducting tests on public datasets. Band selection based on LASSO regression is more explanatory in physics and helps reveal the physical meaning [47]. The logistic regression model (LRM) is a probability classification method [48,49]. As a classifier, it is widely used in hyperspectral image processing, deep learning and other fields, such as medicine, social science, along with others [50,51,52,53]. This study combines the advantages of the sparse feature of the LASSO algorithm and the classification feature of the logistic regression model. Then, an identification method for hyperspectral rice seed purity based on the LASSO logistic regression model (LLRM) is proposed. We selected feature wavelength bands from four types of rice seeds according to different adulteration conditions. Thereafter, we established the prediction model and compared it with the full-band LRM. Finally, we discuss the feasibility of LLRM to select the feature bands and detect the purity of rice seeds.

2. Materials and Methods

2.1. Sample Preparation and Data Acquisition

We selected the following four types of rice seeds with plump grains: Huguangxiang, Xiangwan Japonica, Huanghuazhan, and Japonica 530. Each type contained 72 grains, which were naturally air-dried to be considered as experimental samples. Different types of seeds were neatly placed on the sample table with a spacing of approximately 5 mm to ensure that all the seeds had the same collection environment. Hyperspectral image data of all the seeds were collected four times owing to a large number of seeds.
A schematic of the hyperspectral imaging system for rice seeds is illustrated in Figure 1. The rice seed samples were placed directly on a standard whiteboard. Two bilinear light sources, including six halogen bulbs (H8, OSRAM light AG; Munich, Germany), were symmetrically placed on both sides of the whiteboard for illumination. Furthermore, the hyperspectral camera (Specim IQ, Specim, Spectral Imaging Ltd.; Oulu, Finland) was fixed on the top with a bracket [54]. Specim IQ is a handheld hyperspectral camera, which integrates a processor (NVIDIA Tegra K1, NVIDIA Corporation; Santa Clara, CA, USA), spectral camera (CMOS technology), viewfinder camera (5 Mpix), focus camera, and a scanner with the motor for optics movement. Its spectral resolution is 7 nm with 204 spectral bands across the wavelength range 400–1000 nm, and the image resolution is 512 × 512 px. The original data collected by the hyperspectral camera were light intensity information. These data need to be converted into reflectivity data for subsequent processing. The formula for reflection spectrum conversion is as follows:
R = D R a w − D D a r k D W h i t e − D D a r k ,
where D R a w , D D a r k , and D W h i t e denote the light intensity of the received sample, dark current intensity, and light intensity reflected by the whiteboard, respectively.
Figure 2 depicts the visible image of four rice seed samples obtained through the viewfinder of the hyperspectral camera; they are Huguangxiang, Xiangwan Japonica, Huanghuazhan, and Japonica 530 from the first (top) to fourth (bottom) rows, respectively. Huangxiang and Huanghuazhan are indica rice with a long and narrow shape, whereas Xiangwan Japonica and Japonica 530 are japonica rice with a short and wide shape. Although there are apparent differences in appearance between the japonica and indica rice, the shape features of the two japonica (or indica) varieties are significantly similar. Thus, it is difficult to identify the four varieties based on their image features alone owing to their biological diversity.
The original hyperspectral image also contained other objects, such as the reference whiteboard, in addition to the seed samples. We cut out the irrelevant parts to reduce the amount of data and save data processing resources, retaining only the region of interest (ROI). The original hyperspectral image contained 204 bands of information. The first and last few bands contained no useful information owing to significant noise. Therefore, only the data from the 15th to 190th bands (450–950 nm) were selected for subsequent processing.

2.2. Data Processing Flow

Figure 3 depicts the data processing flowchart of the hyperspectral rice seed purity identification technology. The first step was data pre-processing. The reflectivity data of the ROI were processed by standard normal variate (SNV) transformation and derivative methods. The second step was background segmentation, which involved segmenting the rice seed pixels from the background image and obtaining and analysing the spectral features of each rice seed. The third step involved preparing the training and the test set. The spectral features of 72 pure seeds of each type and the other three hybrid seeds were grouped together to simulate the purity of the seed, producing 13 datasets. These datasets were divided into training and test sets. The fourth step was modelling. LLRM was built based on the training set, and the selected spectral features were used to identify seeds. The last step was data output, wherein the accuracy of seed purity detection was determined based on the classification results.

2.3. Spectral Feature Extraction

2.3.1. Data Pre-Processing

Various noises exist in the original hyperspectral data due to the limitations of the acquisition environment and the instrument, requiring that all data need to be pre-processed [55]. There are many methods suitable for pre-processing hyperspectral data, among which SNV and derivative methods are the most used. SNV is often used to eliminate the influence of sample surface scattering and optical path change on the spectrum [56,57]. Derivative processing includes first-order derivatives (FDs) and second-order derivatives (SDs) to eliminate other background interference and improve spectral resolution. Hence, this study uses SNV, FD, and SD methods to pre-process the data.
The SNV transformation formula is as follows:
S i j ( m ) = x i j ( m ) − x ¯ i j 1 N − 1 ∑ m = 1 N ( x i j ( m ) − x ¯ i j ) 2 ,
where S i j of x i j represents the spectral features before and after the transformation of the pixels in row i and column j , respectively. Furthermore, x ¯ i j denotes the average reflectivity of x i j , which is provided by:
x ¯ i j = ∑ m = 1 N x i j ( m ) N .
FDs and SDs are fulfilled by the Savitzky–Golay smoothing and differentiation filter [58,59], and the calculation formula as follows:
d R ( λ ) = 1 110 Δ λ { 5 [ R ( λ + 5 Δ λ ) − R ( λ − 5 Δ λ ) ] + 4 [ R ( λ + 4 Δ λ ) − R ( λ − 4 Δ λ ) ] − 3 [ R ( λ + 3 Δ λ ) − R ( λ − 3 Δ λ ) ] + 2 [ R ( λ + 2 Δ λ ) − R ( λ − 2 Δ λ ) ] + [ R ( λ + Δ λ ) − R ( λ − Δ λ ) ] } ,
d 2 R ( λ ) = 1 429 Δ λ 2 { 15 [ R ( λ + 5 Δ λ ) + R ( λ − 5 Δ λ ) ] + 6 [ R ( λ + 4 Δ λ ) + R ( λ − 4 Δ λ ) ] − [ R ( λ + 3 Δ λ ) + R ( λ − 3 Δ λ ) ] − 6 [ R ( λ + 2 Δ λ ) + R ( λ − 2 Δ λ ) ] − 9 [ R ( λ + Δ λ ) + R ( λ − Δ λ ) ] − 10 R ( λ ) } ,
where the polynomial order is 2 and smoothing point is 11, and R ( λ ) represents the reflectance data for wavelength λ , and Δ λ = 2.941 nm is the wavelength interval.

2.3.2. Image Segmentation

The rice seed pixels must be separated from the background before the extraction of the spectral features. The traditional grains segmentation method selects a hyperspectral image of a certain wavelength band [60]. The adaptive threshold method was used to obtain the binarised image, which was used as a template. Thereafter, the hyperspectral data were multiplied to obtain image segmentation.
The traditional method only uses the spatial information of hyperspectral data, and it is difficult to obtain an accurate binarised image using only single wavelength band information due to the significant amount of noise in the hyperspectral data. However, the difference in spectral slope between the rice seeds and background was noticeable. The reflectance of rice seeds increases gradually in the 400–1000 nm region, and the spectral curve is relatively smooth, but the spectrum of the standard whiteboard is flat at 400–1000 nm. Therefore, we can use the slope information of the spectrum to synthesise a new greyscale image. The formula of a synthesised greyscale image is as follows:
g i j = a r c t a n ( m i n ( x i j ) m a x ( x i j ) ) ,
where g i j represents the synthesised grey value of the pixel at the i th row and j th column, and x i j denotes the spectrum of pixel points in the i th row and j th column. The arctan function ensures that the value remains between 0 and 1. Thereafter, the synthesised greyscale image was binarised by Otsu’s method to obtain a more accurate sample image boundary. This method uses spectral slope as an effective discriminator, which is equivalent to using both the spatial and spectral information of hyperspectral data to fulfil image segmentation, so it has a better segmentation effect.
Figure 4 depicts the effect of the image segmentation after using the adaptive threshold method to binarise the synthesised greyscale image. Figure 4a depicts a colour image of a rice seed with a cast shadow due to an illumination problem, significantly affecting the image segmentation. Figure 4b depicts the binarised image using a hyperspectral image at 601.55 nm. Figure 4c showcases the result of segmentation using 601.55 nm wavelength band information. Some background pixels are incorrectly segmented due to the influence of the shadow, thereby affecting the spectral features of the rice seeds. Figure 4d shows the synthesised greyscale image using spectral slope. The spectral information helps overcome the effect of the shadow. Figure 4e represents the binarised spectral slope greyscale image. Figure 4f depicts the segmentation result using spectral slope information. Thus, the boundary of the rice seeds can be accurately determined using spectral slope information, yielding a considerably better result than that depicted in Figure 4c.

2.3.3. Spectral Feature

The pixels from each rice seed formed a connected region after the rice seeds were segmented from the background. Thereafter, all the connected areas in the image were counted. The pixels in each connected area belonged to the same rice seed sample. The spectral features of each rice seed were obtained by averaging the spectral data of the same connected area. Figure 5 depicts the average spectral features of four types of rice seeds. Figure 5a illustrates the reflectance spectrum after correction. Figure 5b depicts the spectrum after SNV processes. The average spectral differences among different paddy varieties were not significant for both raw and SNV data. Figure 5c,d show some differences in the spectra of different varieties after the FD and SD processes. Therefore, it is necessary to establish a mathematical classification model to distinguish these bands.

2.4. LLRM Foundation and Optimisation

2.4.1. Data Grouping

The spectral features of 72 pure seeds of each seed type and the other three adulterated seeds were grouped to simulate different cases of seed purity detection, producing 13 datasets. Table 1 presents the detailed information. Groups 1–10 considered Huguangxiang as the research object. Groups 1–4 were adulterated with the other three types of seeds in equal amounts. However, their proportions were different. Groups 5–7 comprised a mix of the other three seeds in different amounts. Groups 8–10 were adulterated with one other seed. Groups 11–13 considered Xiangwan Japonica, Huanghuazhan, and Japonica 530 as the research objects, respectively. Each group was adulterated with three other seeds, and each of them had 12 grains.

2.4.2. Partition of Sample Sets

Sample set partitioning based on joint x-y distance (SPXY) algorithm was used to randomly divide the data into training and test sets, increasing the difference and representativeness between the samples and improving the stability of the model sets [61,62]. The formula to calculate the distance between two samples is as follows:
d x y ( p , q ) = d x ( p , q ) m a x ( d x ( p , q ) ) + d y ( p , q ) m a x ( d y ( p , q ) ) ,
where
d x ( p , q ) = ∑ k = 1 M [ x p ( k ) − x q ( k ) ] 2   and
d y ( p , q ) = | y p − y q | .
Presented here, ( p , q ) ∈ [ 1 , N ] , p and q denote the numbers of two different samples, x p and x q represent the spectral features of samples p and q , respectively. Furthermore, d x and d y denote the spectral distance and category distance of two different samples, respectively.

2.4.3. LLRM

The identification of rice seed purity can be considered a binary classification problem because it does not require the identification of non-target seed varieties. Let us assume that N represents the number of samples of rice seeds with a P -dimensional feature. Thus, the input feature, X, is an N × P matrix, and the output Y = [ y 1 , y 2 , ⋯ y N ] is a binary variable. If the sample is of the target seed variety, y i = 1 ; else, y i = 0 . The logarithmic probability of the i th sample belonging to target seed variety is:
l o g { p ( y i = 1 | X i , β ) 1 − p ( y i = 1 | X i , β ) } = η β ( X i ) ,
where X i = [ x i 1 , x i 2 , ⋯ x i P ] represents the feature of the i th sample, β = [ β 0 , β 1 , β 2 , ⋯ β P ] denotes the regression coefficient of LLRM, and p ( y i = 1 | X i , β ) represents the probability of the i th sample, which belongs to the target seed variety when the regression coefficient is β .
η β ( X i ) = β 0 + β 1 x i 1 + β 2 x i 2 + ⋯ β p x i P = ∑ j = 0 P β j x i j .
The logarithmic maximum likelihood function of LLRM is as follows:
l ( β ) = ∑ i = 1 N { y i η β ( X i ) − l o g { 1 + e x p [ η β ( X i ) ] } } .
LASSO regression minimises the sum of squares of residuals when the sum of absolute values of regression coefficients is lesser than a normal constraint. Some regression coefficients were gradually compressed during the parameter estimation; when the regression coefficient was lower than the threshold, it was set to 0 to achieve variable screening. The L1 norm was introduced into LLRM as the punishment term. Furthermore, the objective of variable sparsity was achieved. Its regression cost function is as follows:
S λ ( β ) = − l ( β ) + λ ‖ β ‖ 1 .
The estimated regression coefficient, β ^ ( λ ) , was obtained by minimising Equation (7) in the model.
β ^ ( λ ) = a r g m i n β ∑ i = 1 N { y i η β ( X i ) − l o g { 1 + exp [ η β ( X i ) ] } } + λ ∑ j = 1 P | β j | ,
where λ represents the regularisation parameter. The smaller the value of λ , the smaller the punishment, retaining more variables. However, the larger the value of λ , the larger the punishment, retaining fewer variables. Therefore, the selection of an appropriate λ has a significant impact on the classification accuracy of LLRM. This study uses a ten-fold cross-validation method to select the optimal value of λ and can be calculated as follows:
λ ^ C V = a r g m i n λ > 0 K − 1 C V k ( λ ) ,
where K = 10 , C V k ( λ ) represents an error in fitting the kth subset when the optimal value of λ is used.
The detailed steps of the LLRM algorithm implementation are as follows:
(1)
Data pre-processing: The data collected by hyperspectral technology were pre-processed. Then, they were grouped, as presented in Table 1.
(2)
Dividing the training and test sets: The SPXY algorithm was used to randomly divide the data into training and test sets with a 1:1 ratio.
(3)
Selecting the regularisation parameter: The training samples were divided into 10 parts using the cross-validation method, taking turns using 9 of their train, and then calculating the error of fitting the other one. The mean square error of the 10 prediction results was used to estimate the accuracy of the algorithm. Moreover, the optimal value of λ was selected based on the accuracy of the model prediction and by combining the number of feature bands.
(4)
Selecting the feature wavelength bands: The coordinate descent method was used to calculate the regression coefficient at an optimal value of λ . Then, the feature band was selected based on the regression coefficient.
(5)
Modelling: The feature wavelength bands from the test set were selected. Moreover, the regression model was considered the foundation.
(6)
Calculation of classification accuracy: The function is calculated as follows:
p ( i ) = 1 / ( 1 + e x p ( − y i ) ) .
If the probability is greater than 0.5, the sample represents the target seed variety. However, if the probability is less than 0.5, the sample represents a non-target seed variety. Finally, they were compared with the labels and the classification accuracy was determined.

2.4.4. Optimum λ Selection

The value of λ has a considerable influence in solving the LLRM. An optimum value of λ can establish a balance between the fitting error and model complexity. If λ is significantly large, the degree of model compression increases, the number of selected feature wavelength bands decreases, and the fitting error increases, consequently affecting classification. If λ is significantly small, the variable compression effect is constrained, and the number of selected feature bands increases significantly, wasting computing resources. Therefore, the value of λ must be optimised carefully.
According to the number of feature wavelengths and recognition accuracy, the optimisation range of λ is set to 0.0001–0.2 in this study. We used cross-validation to determine the relatively small range of λ [63,64]. Thereafter, the number of feature bands and the accuracies for different λ values were analysed, and the optimal λ value was selected. Figure 6a depicts the change in the regression coefficient with λ . The relative value of the regression coefficient gradually decreases, and the coefficient is gradually compressed with increasing λ . Figure 6b depicts variation in the data fitting error of Group 2 with the change in λ . The blue curve represents the fitting error, and the thick red line indicates the standard deviation of different fitting errors. The red vertical line corresponds to the value of λ with the smallest fitting error. In general, the smallest fitting error does not imply the best prediction effect because the data contain unknown noise. Therefore, the optimal value of λ was expected to be found near the red vertical line. Figure 6c,d show that the number of feature wavelength bands is 10 when the fitting error is the smallest with a prediction accuracy of 92.57%. However, the prediction accuracy does not decrease with the decreasing number of feature wavelength bands. This indicates the need to further optimise the model. The highest prediction accuracy of 95.37% is obtained below the red vertical line. There are two positions with the highest prediction rate of λ = 0.0270 and λ = 0.0668, which correspond to 10 and 7 feature wavelength bands, respectively. Moreover, the green vertical line corresponds to λ = 0.0668, which requires the minimum number of feature wavelength bands while exhibiting the same prediction accuracy, thereby effectively reducing the cost of data processing. Other data groups can be optimised similarly to obtain the best LLRM.

3. Results

3.1. Results of Wavelength Bands Selection

Table 2 presents the wavelength bands selected using the SNV data under different adulteration groups. Groups 1–10 indicate that the coincidence degree of the selected wavelength bands is significantly higher when the same variety of rice seeds is considered despite the different number and variety of doped seeds. Among them, wavelength bands of 678.71 nm, 810.86 nm, and 871.60 nm appeared in all the data groups. Furthermore, the appearance frequencies of 548.55 nm, 798.77 nm, and 847.25 nm were also significantly higher, indicating the feasibility of using LLRM for band selection. Groups 11–13 show that the coincidence degree of the selected wavelength bands is significantly smaller, implying that different LLRMs must be established for different rice varieties.
Table 3 and Table 4 list the selected wavelength bands using FD and SD data under different adulteration groups, respectively. The coincidence degree of the selected wavelength bands is significantly higher for the same variety of rice seeds. In contrast, it is significantly smaller for different varieties of rice seeds. We compared the results of band selection using SNV, FD, and SD data for the same type of target seeds. The results indicate that the band coincidence is extremely low, implying that the band selection of LLRM is highly sensitive to the data pre-processing procedure. Under the optimal LLRM, the number of feature wavelength bands selected by the SNV, FD, and SD data were 5–17, 3–11, and 6–12, respectively. Moreover, the FD data required the least number of bands among the three.

3.2. Band Selection

Figure 7 depicts the wavelength band selection using the SNV data for different rice varieties. The red lines indicate the average spectra of the target rice seeds, and the blue lines represent the average spectra of the adulterated seeds. The positions indicated by the green lines represent the best bands based on the data of Group 2. Figure 7a shows that the wavelength bands of 548.55 nm, 714.55 nm, and 810.86 nm correspond to small reflectance peaks of the target seeds. The reflectance of the target seed spectra at wavelengths of 643.01 nm, 678.71 nm, and 871.60 nm are significantly different from those of the adulterated seeds. The difference in spectra between the target rice seeds and adulterated seeds are not noticeable at 762.98 nm, 798.77 nm, and 847.25 nm. Certain wavelength bands, such as 896.01 nm and 920.48 nm, are not selected despite the noticeable difference in their spectra, implying that greater spectral differences do not contain more useful information for classification. The number and location of the selected bands in Figure 7b–d are significantly different, indicating that different LLRMs must be established for different varieties.
Figure 8 depicts the wavelength band selection using FD data for different rice varieties. Furthermore, Figure 8a shows that there are significant differences in the data at 466.77 nm, 583.85 nm, 678.71 nm, and 835.11 nm. The greater spectral differences do not imply a greater variety of samples. Comparing Figure 8b–d, we find that the selected wavelength bands are mainly concentrated in the infrared wavelengths of 750–900 nm. Moreover, the spectral differences cannot be observed directly in the hyperspectral data, indicating that the infrared bands contain more useful information than visible light bands. Thus, we conclude that more attention should be allotted to the infrared bands for wavelength band selection.
Figure 9 depicts the wavelength band selection using SD data for different rice varieties. It is difficult to notice the spectral differences in different wavelength bands because the SD data is significantly rough. Hence, the wavelengths with large peaks are not selected. Comparing Figure 7, Figure 8 and Figure 9, we find that the selected feature wavelength bands are almost different under different pre-processing conditions for the same type of target seeds, indicating that the data pre-processing method significantly influences the LLRM. The pre-processing using the FD data exhibited the fastest convergence speed and required the least number of bands; hence, it is the best pre-processing method.
The coincidence degrees of the feature bands selected for different data groups were significantly higher in the Huguangxiang case. The appearance frequencies of 548.55 nm, 643.01 nm, 678.71 nm, 714.55 nm, 810.86 nm, and 871.60 nm were the highest. The coincidence degrees of the feature bands were significantly smaller for different types of rice seeds, indicating that a specific LLRM must be established for different rice varieties. The band selection results of the LLRM were highly sensitive to data pre-processing. It was found that pre-processing with FD was more suitable for LLRM than with SNV or SD.

3.3. Comparison of Model Accuracy

Table 5, Table 6 and Table 7 show the recognition accuracies between LLRM and LRM using SNV, FD an SD data, respectively. As shown in Table 5, the recognition accuracy of LLRM is 98.77–100%, with an average recognition accuracy of 99.64%. The recognition accuracy of the full variable LRM is 71.60–100%, with an average recognition rate of 93.66%. However, the results of LRM are considerably lower than those of LLRM. The recognition accuracies of Group 3 and 4 are particularly more prominent at less than 80%; this indicates that LLRM selects the key bands steadily and eliminates unnecessary noise, thereby improving the signal-to-noise ratio of the hyperspectral data and obtaining a better rice seed recognition rate.
Table 6 presents the comparison of the recognition accuracies between LLRM and LRM using FD data. The recognition accuracy of LLRM is 95.37–100%, with an average recognition accuracy of 99.18%. However, the recognition accuracy of the full variable LRM is 73.33–95.83%, with an average recognition rate of merely 88.63%. Table 7 shows that the recognition accuracy of the full variable LRM using SD data is 77.78–95.61%, with an average recognition rate of merely 86.60%. However, the recognition accuracy of LLRM is 91.67–100%, with an average recognition rate of 96.59%. The average recognition rate of LLRM is higher than the highest rate of LRM. Comparing Table 5 with Table 6 and Table 7 shows that the recognition accuracy of the full variable LRM is 78.89–95.61%, with an average recognition rate of merely 86.60%. However, the recognition accuracy of LLRM is 91.67–100%, with an average recognition rate of 96.59%. This result indicates that it is feasible to develop a hyperspectral technology with LLRM to identify seed purity.

4. Conclusions

In this study, we combined the sparse characteristics of the LASSO algorithm and the classification characteristics of the LRM and proposed a method to identify the hyperspectral rice seed purity based on the LLRM. We tested 13 different adulterated groups of four varieties of rice seeds. The results of different data groups indicated that the value of the regularisation parameter λ was different for different adulteration cases. Moreover, the number of selected feature wavelength bands also varied. The best balance between the number of feature wavelength bands and recognition accuracy was obtained by adjusting λ , thereby establishing the optimal LLRM.
The recognition accuracy of the LLRM was 91.67–100% in 13 rice seed adulteration cases, with an average recognition rate of 98.47%. However, the recognition accuracy of the LRM was 71.60–100%, with an average accuracy of merely 89.63%. This result indicates that LLRM selects the key bands steadily and improves the signal-to-noise ratio of the hyperspectral data, obtaining better accuracy in rice seed recognition. The average number of selected feature wavelength bands is 11.3 for SNV data, 8.5 for FD data, and 9.5 for SD data. SNV data selected relatively more feature wavelength bands, which led to better accuracy than FD and SD data. By comparing all the LLRMs, we found that the accuracy of three pre-processing methods was improved by 5.98%, 10.55%, and 9.99%, respectively, and the probability of achieving a 100% recognition accuracy increased after FD processing, implying that FD processing is more suitable than other approaches for LLRM.
These results reveal that the hyperspectral detection technology of rice seed purity based on LLRM can reduce the collection cost and improve the processing speed while ensuring satisfactory recognition accuracy. Thus, LLRM demonstrates significant potential for application in the field of seed purity identification. In future research, we will further analyse and select these candidate feature wavelength bands and add information on the texture and shape of samples to establish a more robust hyperspectral rice seed purity identification technology.

Author Contributions

Conceptualization, S.Z.; methodology, W.L. and F.C.; software, G.W. and H.L.; formal analysis, W.L. and G.W.; writing—original draft preparation, W.L.; writing—review and editing, G.W. and H.L.; supervision, S.Z. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the National Natural Science Foundation of China, grant number U1833119 and 42074172, and funded by National Food and Strategic Reserves Administration Foundation, grant number LQ2018501, and funded by Hubei province Natural Science Foundation for Distinguished Young Scholars, grant number 2020CFA063.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The data that support the findings of this study are available upon request from the authors.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Yang, J.; Sun, J.; Du, L.; Chen, B.; Zhang, Z.; Shi, S.; Gong, W. Monitoring of Paddy Rice Varieties Based on the Combination of the Laser-Induced Fluorescence and Multivariate Analysis. Food Anal. Methods 2017, 10, 2398–2403. [Google Scholar] [CrossRef]
  2. Cui, S.; Ma, X.; Wang, X.; Zhang, T.-A.; Hu, J.; Tsang, Y.F.; Gao, M.-T. Phenolic acids derived from rice straw generate peroxides which reduce the viability of Staphylococcus aureus cells in biofilm. Ind. Crops Prod. 2019, 140, 111561. [Google Scholar] [CrossRef]
  3. Nethra, N.; Prasad, S.R.; Vishwanath, K.; Dhanraj, K.; Gowda, R. Identification of rice hybrids and their parental lines based on seed, seedling characters, chemical tests and gel electrophoresis of total soluble seed proteins. Seed Sci. Technol. 2007, 35, 176–186. [Google Scholar] [CrossRef]
  4. Ellepola, S.W.; Choi, S.M.; Ma, C.Y. Conformational study of globulin from rice (Oryza sativa) seeds by Fourier-transform infrared spectroscopy. Int. J. Biol. Macromol. 2005, 37, 12–20. [Google Scholar] [CrossRef] [PubMed]
  5. Chen, J.; Li, M.; Pan, T.; Pang, L.; Yao, L.; Zhang, J. Rapid and non-destructive analysis for the identification of multi-grain rice seeds with near-infrared spectroscopy. Spectrochim. Acta Mol. Biomol. Spectrosc. 2019, 219, 179–185. [Google Scholar] [CrossRef]
  6. Aznan, A.A.; Rukunudin, I.H.; Shakaff, A.; Ruslan, R.; Saad, F. The use of machine vision technique to classify cultivated rice seed variety and weedy rice seed variants for the seed industry. Int. Food Res. J. 2016, 23, 31–35. [Google Scholar]
  7. Kiratiratanapruk, K.; Temniranrat, P.; Sinthupinyo, W.; Prempree, P.; Chaitavon, K.; Porntheeraphat, S.; Prasertsak, A. Development of Paddy Rice Seed Classification Process Using Machine Learning Techniques for Automatic Grading Machine. J. Sens. 2020, 2020, 7041310. [Google Scholar] [CrossRef]
  8. Liu, Z.-Y.; Cheng, F.; Ying, Y.-B.; Rao, X.-Q. Identification of rice seed varieties using neural network. J. Zhejiang Univ. Sci. B 2005, 6, 1095–1100. [Google Scholar] [CrossRef] [Green Version]
  9. Fotiadou, K.; Tsagkatakis, G.; Tsakalides, P. Deep Convolutional Neural Networks for the Classification of Snapshot Mosaic Hyperspectral Imagery. Electron. Imaging 2017, 2017, 185–190. [Google Scholar] [CrossRef] [Green Version]
  10. Fabiyi, S.D.; Vu, H.; Tachtatzis, C.; Murray, P.; Harle, D.; Dao, T.K.; Andonovic, I.; Ren, J.; Marshall, S. Varietal Classification of Rice Seeds Using RGB and Hyperspectral Images. IEEE Access 2020, 8, 22493–22505. [Google Scholar] [CrossRef]
  11. Zhang, L.; Sun, H.; Rao, Z.; Ji, H. Hyperspectral imaging technology combined with deep forest model to identify frost-damaged rice seeds. Spectrochim. Acta Mol. Biomol. Spectrosc. 2020, 229, 117973. [Google Scholar] [CrossRef]
  12. Zhang, J.; Yang, Y.; Feng, X.; Xu, H.; Chen, J.; He, Y. Identification of Bacterial Blight Resistant Rice Seeds Using Terahertz Imaging and Hyperspectral Imaging Combined with Convolutional Neural Network. Front. Plant Sci. 2020, 11, 821. [Google Scholar] [CrossRef]
  13. Yang, Y.; Chen, J.; He, Y.; Liu, F.; Feng, X.; Zhang, J. Assessment of the vigor of rice seeds by near-infrared hyperspectral imaging combined with transfer learning. RSC Adv. 2020, 10, 44149–44158. [Google Scholar] [CrossRef]
  14. Zhou, L.; Zhang, C.; Taha, M.F.; Wei, X.; He, Y.; Qiu, Z.; Liu, Y. Wheat Kernel Variety Identification Based on a Large Near-Infrared Spectral Dataset and a Novel Deep Learning-Based Feature Selection Method. Front. Plant Sci. 2020, 11, 575810. [Google Scholar] [CrossRef] [PubMed]
  15. Kong, W.; Zhang, C.; Liu, F.; Nie, P.; He, Y. Rice Seed Cultivar Identification Using Near-Infrared Hyperspectral Imaging and Multivariate Data Analysis. Sensors 2013, 13, 8916–8927. [Google Scholar] [CrossRef] [Green Version]
  16. Kuo, T.-Y.; Chung, C.-L.; Chen, S.-Y.; Lin, H.-A.; Kuo, Y.-F. Identifying rice grains using image analysis and sparse-representation-based classification. Comput. Electron. Agric. 2016, 127, 716–725. [Google Scholar] [CrossRef]
  17. Anami, B.S.; Malvade, N.N.; Palaiah, S.; Naveen, N.; Surendra, P. Automated recognition and classification of adulteration levels from bulk paddy grain samples. Inf. Process. Agric. 2019, 6, 47–60. [Google Scholar] [CrossRef]
  18. Chaugule, A.; Mali, S.N. Evaluation of Texture and Shape Features for Classification of Four Paddy Varieties. J. Eng. 2014, 2014, 617263. [Google Scholar] [CrossRef] [Green Version]
  19. Chaugule, A.A.; Mali, S.N. Identification of paddy varieties based on novel seed angle features. Comput. Electron. Agric. 2016, 123, 415–422. [Google Scholar] [CrossRef]
  20. Huang, K.-Y.; Chien, M.-C. A Novel Method of Identifying Paddy Seed Varieties. Sensors 2017, 17, 809. [Google Scholar] [CrossRef]
  21. Ansari, N.; Ratri, S.S.; Jahan, A.; Ashik-E-Rabbani, M.; Rahman, A. Inspection of paddy seed varietal purity using machine vision and multivariate analysis. J. Agric. Food Res. 2021, 3, 100109. [Google Scholar] [CrossRef]
  22. Vithu, P.; Moses, J.A. Machine vision system for food grain quality evaluation: A review. Trends Food Sci. Technol. 2016, 56, 13–20. [Google Scholar] [CrossRef]
  23. Bai, Z.; Hu, X.; Tian, J.; Chen, P.; Luo, H.; Huang, D. Rapid and nondestructive detection of sorghum adulteration using optimization algorithms and hyperspectral imaging. Food Chem. 2020, 331, 127290. [Google Scholar] [CrossRef] [PubMed]
  24. Okada, N.; Maekawa, Y.; Owada, N.; Haga, K.; Shibayama, A.; Kawamura, Y. Automated Identification of Mineral Types and Grain Size Using Hyperspectral Imaging and Deep Learning for Mineral Processing. Minerals 2020, 10, 809. [Google Scholar] [CrossRef]
  25. Femenias, A.; Gatius, F.; Ramos, A.J.; Sanchis, V.; Marín, S. Use of hyperspectral imaging as a tool for Fusarium and deoxynivalenol risk management in cereals: A review. Food Control 2020, 108, 106819. [Google Scholar] [CrossRef]
  26. Serranti, S.; Cesare, D.; Marini, F.; Bonifazi, G. Classification of oat and groat kernels using NIR hyperspectral imaging. Talanta 2013, 103, 276–284. [Google Scholar] [CrossRef]
  27. Yang, X.; Hong, H.; You, Z.; Cheng, F. Spectral and Image Integrated Analysis of Hyperspectral Data for Waxy Corn Seed Variety Classification. Sensors 2015, 15, 15578–15594. [Google Scholar] [CrossRef] [Green Version]
  28. Zhu, S.; Zhang, J.; Chao, M.; Xu, X.; Song, P.; Zhang, J.; Huang, Z. A Rapid and Highly Efficient Method for the Identification of Soybean Seed Varieties: Hyperspectral Images Combined with Transfer Learning. Molecules 2020, 25, 152. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  29. Qiu, Z.; Chen, J.; Zhao, Y.; Zhu, S.; He, Y.; Zhang, C. Variety Identification of Single Rice Seed Using Hyperspectral Imaging Combined with Convolutional Neural Network. Appl. Sci. 2018, 8, 212. [Google Scholar] [CrossRef] [Green Version]
  30. Weng, S.; Tang, P.; Yuan, H.; Guo, B.; Yu, S.; Huang, L.; Xu, C. Hyperspectral imaging for accurate determination of rice variety using a deep learning network with multi-feature fusion. Spectrochim. Acta Mol. Biomol. Spectrosc. 2020, 234, 118237. [Google Scholar] [CrossRef] [PubMed]
  31. Sun, J.; Lu, X.; Mao, H.; Jin, X.; Wu, X. A Method for Rapid Identification of Rice Origin by Hyperspectral Imaging Technology. J. Food Process Eng. 2017, 40, e12297. [Google Scholar] [CrossRef] [Green Version]
  32. Liu, W.; Zeng, S.; Li, H.; Xiao, Z.; Huang, X.; Jiang, H. Paddy variety identification using hyperspectral imagery under non-ideal illumination conditions. Int. J. Wavelets Multiresolut. Inf. Process. 2020, 18, 2050036. [Google Scholar] [CrossRef]
  33. Chin, J.A.; Wang, E.C.; Kibbe, M.R. Evaluation of hyperspectral technology for assessing the presence and severity of peripheral artery disease. J. Vasc. Surg. 2011, 54, 1679–1688. [Google Scholar] [CrossRef] [Green Version]
  34. Dai, Q.; Cheng, J.-H.; Sun, D.-W.; Zeng, X.-A. Advances in Feature Selection Methods for Hyperspectral Image Processing in Food Industry Applications: A Review. Crit. Rev. Food Sci. Nutr. 2015, 55, 1368–1382. [Google Scholar] [CrossRef] [PubMed]
  35. Baek, I.; Kim, M.S.; Cho, B.-K.; Mo, C.; Barnaby, J.Y.; McClung, A.M.; Oh, M. Selection of Optimal Hyperspectral Wavebands for Detection of Discolored, Diseased Rice Seeds. Appl. Sci. 2019, 9, 1027. [Google Scholar] [CrossRef] [Green Version]
  36. Zhang, W.; Li, X.; Zhao, L. A Fast Hyperspectral Feature Selection Method Based on Band Correlation Analysis. IEEE Geosci. Remote Sens. Lett. 2018, 15, 1750–1754. [Google Scholar] [CrossRef]
  37. Mehmood, T.; Liland, K.H.; Snipen, L.; Sæbø, S. A review of variable selection methods in Partial Least Squares Regression. Chemom. Intell. Lab. Syst. 2012, 118, 62–69. [Google Scholar] [CrossRef]
  38. Rodarmel, C.; Shan, J. Principal component analysis for hyperspectral image classification. Surv. Land Inf. Sci. 2002, 62, 115–122. [Google Scholar]
  39. Kang, X.; Duan, P.; Li, S. Hyperspectral image visualization with edge-preserving filtering and principal component analysis. Inf. Fusion 2020, 57, 130–143. [Google Scholar] [CrossRef]
  40. Uddin, M.P.; Mamun, M.A.; Afjal, M.I.; Hossain, M.A. Information-theoretic feature selection with segmentation-based folded principal component analysis (PCA) for hyperspectral image classification. Int. J. Remote Sens. 2021, 42, 286–321. [Google Scholar] [CrossRef]
  41. Machidon, A.L.; Del Frate, F.; Picchiani, M.; Machidon, O.M.; Ogrutan, P.L. Geometrical Approximated Principal Component Analysis for Hyperspectral Image Analysis. Remote Sens. 2020, 12, 1698. [Google Scholar] [CrossRef]
  42. Hashjin, S.S.; Khazai, S. A new method to detect targets in hyperspectral images based on principal component analysis. Geocarto Int. 2020, 1–19. [Google Scholar] [CrossRef]
  43. Tibshirani, R. Regression Shrinkage and Selection Via the Lasso. J. R. Stat. Soc. 1996, 58, 267–288. [Google Scholar] [CrossRef]
  44. Uezato, T.; Fauvel, M.; Dobigeon, N. Hyperspectral Unmixing with Spectral Variability Using Adaptive Bundles and Double Sparsity. IEEE Trans. Geosci. Remote Sens. 2019, 57, 3980–3992. [Google Scholar] [CrossRef] [Green Version]
  45. Salehani, Y.E.; Gazor, S.; Cheriet, M. A new weighted ℓp-norm for sparse hyperspectral unmixing. In Proceedings of the 2017 IEEE Global Conference on Signal and Information Processing, Montreal, QC, Canada, 14–16 November 2017; pp. 1255–1259. [Google Scholar]
  46. Sun, L.; Ge, W.; Chen, Y.; Zhang, J.; Jeon, B. Hyperspectral unmixing employing l1–l2 sparsity and total variation regularization. Int. J. Remote Sens. 2018, 39, 6037–6060. [Google Scholar] [CrossRef]
  47. Yang, D.; Bao, W. Group Lasso-Based Band Selection for Hyperspectral Image Classification. IEEE Geosci. Remote Sens. Lett. 2017, 14, 2438–2442. [Google Scholar] [CrossRef]
  48. Dankmar, B. Multinomial logistic regression algorithm. Ann. Inst. Stat. Math. 1992, 44, 197–200. [Google Scholar]
  49. Krishnapuram, B.; Carin, L.; Figueiredo, M.A.T.; Hartemink, A.J. Sparse multinomial logistic regression: Fast algorithms and generalization bounds. IEEE Trans. Pattern Anal. Mach. Intell. 2005, 27, 957–968. [Google Scholar] [CrossRef] [Green Version]
  50. Haut, J.M.; Paoletti, M.E. Cloud Implementation of Multinomial Logistic Regression for UAV Hyperspectral Images. IEEE J. Miniat. Air Space Syst. 2020, 1, 163–171. [Google Scholar] [CrossRef]
  51. Mohammadi, A.; Shaverizade, A. Ensemble deep learning for aspect-based sentiment analysis. Int. J. Nonlinear Anal. Appl. 2021, 12, 29–38. [Google Scholar]
  52. Bayaga, A. Multinomial logistic regression: Usage and application in risk analysis. J. Appl. Quant. Methods 2010, 5, 288–297. [Google Scholar]
  53. Peng, C.-Y.J.; Lee, K.L.; Ingersoll, G.M. An Introduction to Logistic Regression Analysis and Reporting. J. Educ. Res. 2002, 96, 3–14. [Google Scholar] [CrossRef]
  54. Behmann, J.; Acebron, K.; Emin, D.; Bennertz, S.; Matsubara, S.; Thomas, S.; Bohnenkamp, D.; Kuska, M.T.; Jussila, J.; Salo, H.; et al. Specim IQ: Evaluation of a New, Miniaturized Handheld Hyperspectral Camera and Its Application for Plant Phenotyping and Disease Detection. Sensors 2018, 18, 441. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  55. Qian, Y.; Ye, M. Hyperspectral Imagery Restoration Using Nonlocal Spectral-Spatial Structured Sparse Representation with Noise Estimation. IEEE J. Sel. Top. Appl. Earth Obs. Remote Sens. 2012, 6, 499–515. [Google Scholar] [CrossRef]
  56. Barnes, R.J.; Dhanoa, M.S.; Lister, S.J. Standard Normal Variate Transformation and De-Trending of Near-Infrared Diffuse Reflectance Spectra. Appl. Spectrosc. 1989, 43, 772–777. [Google Scholar] [CrossRef]
  57. Rinnan, Å.; van den Berg, F.; Engelsen, S.B. Review of the most common pre-processing techniques for near-infrared spectra. TrAC Trends Anal. Chem. 2009, 28, 1201–1222. [Google Scholar] [CrossRef]
  58. Savitzky, A.; Golay, M.J.E. Smoothing and Differentiation of Data by Simplified Least Squares Procedures. Anal. Chem. 1964, 36, 1627–1639. [Google Scholar] [CrossRef]
  59. Luo, J.; Ying, K.; Bai, J. Savitzky-Golay smoothing and differentiation filter for even number data. Signal Process. 2005, 85, 1429–1434. [Google Scholar] [CrossRef]
  60. Xia, C.; Yang, S.; Huang, M.; Zhu, Q.; Guo, Y.; Qin, J. Maize seed classification using hyperspectral image coupled with multi-linear discriminant analysis. Infrared Phys. Technol. 2019, 103, 103077. [Google Scholar] [CrossRef]
  61. Galvão, R.K.H.; Araujo, M.C.U.; José, G.E.; Pontes, M.J.C.; Silva, E.C.; Saldanha, T.C.B. A method for calibration and validation subset partitioning. Talanta 2005, 67, 736–740. [Google Scholar] [CrossRef] [PubMed]
  62. Tian, H.; Zhang, L.; Li, M.; Wang, Y.; Sheng, D.; Liu, J.; Wang, C. Weighted SPXY method for calibration set selection for composition analysis based on near-infrared spectroscopy. Infrared Phys. Technol. 2018, 95, 88–92. [Google Scholar] [CrossRef]
  63. Takayama, T.; Iwasaki, A. Optimal wavelength selection on hyperspectral data with fused lasso for biomass estimation of tropical rain forest. ISPRS Ann. Photogramm. Remote Sens. Spat. Inf. Sci. 2016, 3, 101–108. [Google Scholar] [CrossRef] [Green Version]
  64. Sirimongkolkasem, T.; Drikvandi, R. On Regularisation Methods for Analysis of High Dimensional Data. Ann. Data Sci. 2019, 6, 737–763. [Google Scholar] [CrossRef] [Green Version]
Sensors 21 04384 g001 550
Figure 1. Schematic of the hyperspectral imaging system for rice seeds.
Figure 1. Schematic of the hyperspectral imaging system for rice seeds.
Sensors 21 04384 g001
Sensors 21 04384 g002 550
Figure 2. Image of four rice seed samples through the viewfinder of the hyperspectral camera; the first (top) to fourth (bottom) rows are Huguangxiang, Xiangwan Japonica, Huanghuazhan, and Japonica 530, respectively.
Figure 2. Image of four rice seed samples through the viewfinder of the hyperspectral camera; the first (top) to fourth (bottom) rows are Huguangxiang, Xiangwan Japonica, Huanghuazhan, and Japonica 530, respectively.
Sensors 21 04384 g002
Sensors 21 04384 g003 550
Figure 3. Flowchart of hyperspectral image processing, including data pre-processing, background segmentation, data set preparation, modelling, and identification accuracy output.
Figure 3. Flowchart of hyperspectral image processing, including data pre-processing, background segmentation, data set preparation, modelling, and identification accuracy output.
Sensors 21 04384 g003
Sensors 21 04384 g004 550
Figure 4. Example of background segmentation using synthesised spectral slope greyscale image. (a) Colour image of a seed; (b) Binarised image using single-band (601.55 nm) information; (c) Segmentation result using single-band (601.55 nm) information; (d) Synthesised spectral slope greyscale image; (e) Binarised image using spectral slope information; (f) Segmentation result using spectral slope information.
Figure 4. Example of background segmentation using synthesised spectral slope greyscale image. (a) Colour image of a seed; (b) Binarised image using single-band (601.55 nm) information; (c) Segmentation result using single-band (601.55 nm) information; (d) Synthesised spectral slope greyscale image; (e) Binarised image using spectral slope information; (f) Segmentation result using spectral slope information.
Sensors 21 04384 g004
Sensors 21 04384 g005a 550Sensors 21 04384 g005b 550
Figure 5. Average spectral features of four rice seeds. (a) Original reflectance spectrum; (b) spectrum after SNV; (c) spectrum after FD’ (d) spectrum after SD.
Figure 5. Average spectral features of four rice seeds. (a) Original reflectance spectrum; (b) spectrum after SNV; (c) spectrum after FD’ (d) spectrum after SD.
Sensors 21 04384 g005aSensors 21 04384 g005b
Sensors 21 04384 g006 550
Figure 6. Change in parameters with different λ values. (a) Regression coefficient; (b) fitting deviation; (c) number of feature wavelength bands; (d) prediction accuracy.
Figure 6. Change in parameters with different λ values. (a) Regression coefficient; (b) fitting deviation; (c) number of feature wavelength bands; (d) prediction accuracy.
Sensors 21 04384 g006
Sensors 21 04384 g007 550
Figure 7. Wavelength band selection using SNV data for different rice varieties under the optimal λ. (a) Huguangxiang; (b) Xiangwan Japonica; (c) Huanghuazhan; (d) Japonica 530.
Figure 7. Wavelength band selection using SNV data for different rice varieties under the optimal λ. (a) Huguangxiang; (b) Xiangwan Japonica; (c) Huanghuazhan; (d) Japonica 530.
Sensors 21 04384 g007
Sensors 21 04384 g008 550
Figure 8. Wavelength band selection using FD data for different rice varieties under the optimal λ. (a) Huguangxiang. (b) Xiangwan Japonica. (c) Huanghuazhan. (d) Japonica 530.
Figure 8. Wavelength band selection using FD data for different rice varieties under the optimal λ. (a) Huguangxiang. (b) Xiangwan Japonica. (c) Huanghuazhan. (d) Japonica 530.
Sensors 21 04384 g008
Sensors 21 04384 g009 550
Figure 9. Wavelength band selection using SD data for different rice varieties under the optimal λ . (a) Huguangxiang; (b) Xiangwan Japonica; (c) Huanghuazhan; (d) Japonica 530.
Figure 9. Wavelength band selection using SD data for different rice varieties under the optimal λ . (a) Huguangxiang; (b) Xiangwan Japonica; (c) Huanghuazhan; (d) Japonica 530.
Sensors 21 04384 g009
Table
Table 1. Data grouping details to simulate different cases of seed purity detection.
Table 1. Data grouping details to simulate different cases of seed purity detection.
Unit/GrainType 1Type 2Type 3Type 4Sum
Group 172242424144
Group 272121212108
Group 37266690
Group 47233381
Group 57224126114
Group 67212246114
Group 77261224114
Group 8723600108
Group 9720360108
Group 10720036108
Group 1112721212108
Group 1212127212108
Group 1312121272108
Table
Table 2. The results of feature wavelength bands selected using SNV data under different adulteration groups.
Table 2. The results of feature wavelength bands selected using SNV data under different adulteration groups.
GroupFeature Wavelength Bands [nm]Number
Group 1 548.55, 572.07, 643.01, 678.71, 714.55, 762.57, 798.77, 810.86, 822.98, 847.25, 871.60, 896.01, 908.24, 920.4814
Group 2 548.55, 643.01, 678.71, 714.55, 762.57, 798.77, 810.86, 847.25, 871.609
Group 3 548.55, 631.15, 643.01, 678.71, 762.57, 798.77, 810.86, 847.25, 871.60, 920.4810
Group 4 548.55, 619.30, 678.71, 714.55, 798.77, 810.86, 871.607
Group 5 548.55, 572.07, 583.85, 643.01, 678.71, 714.55, 762.57, 786.68, 810.86, 822.98, 847.25, 871.60, 896.01, 920.4814
Group 6 536.82, 548.55, 619.30, 678.71, 714.55, 762.57, 798.77, 810.86, 847.25 871.60, 896.0111
Group 7 548.55, 643.01, 678.71, 714.55, 750.54, 762.57, 798.77, 810.86, 847.25, 871.60, 957.3211
Group 8 548.55, 560.30, 572.07, 631.15, 643.01, 678.71, 762.57, 810.86, 847.25, 871.60, 896.01, 908.24, 920.4813
Group 9 525.10, 536.82, 619.30, 678.71, 726.53, 786.68, 798.77, 810.86, 847.25, 871.60, 896.0111
Group 10 560.30, 643.01, 678.71, 714.55, 798.77, 810.86, 871.60, 908.24, 920.48, 957.3210
Group 11 455.16, 501.72, 525.10, 548.55, 572.07, 583.85, 631.15, 702.58, 714.55, 726.53, 738.53, 762.57, 786.68, 798.77, 810.86, 871.60, 957.3217
Group 12 443.56, 513.40, 560.30, 643.01, 678.71, 702.58, 726.53, 774.62, 786.68, 810.86, 822.98, 835.11, 859.42, 908.24, 932.7415
Group 13 443.56, 513.40, 738.53, 810.86, 822.985
Table
Table 3. The results of feature wavelength bands selected using FD data under different adulteration groups.
Table 3. The results of feature wavelength bands selected using FD data under different adulteration groups.
GroupFeature Wavelength Bands [nm]Number
Group 1 466.77, 583.85, 607.46, 678.71, 750.54, 822.98, 835.11, 847.25, 859.42, 883.79, 932.7411
Group 2 466.77, 583.85, 607.46, 678.71, 786.68, 810.86, 835.11, 859.42, 883.79, 932.7410
Group 3 455.16, 583.85, 666.79, 678.71, 810.86, 835.11, 859.42, 871.60, 883.79, 932.7410
Group 4 466.77, 583.85, 678.71, 810.86, 835.11, 859.42, 883.797
Group 5 466.77, 583.85, 678.71, 835.11, 859.42, 883.796
Group 6 466.77, 654.89, 714.55, 822.98, 859.42, 883.79, 932.747
Group 7 455.16, 466.77, 525.10, 548.55, 607.46, 678.71, 786.68, 835.11, 859.42, 883.79, 932.7411
Group 8 466.77, 583.85, 678.71, 835.11, 859.42, 883.796
Group 9 466.77, 501.72, 619.30, 702.58, 714.55, 822.98, 835.11, 847.25, 859.42, 883.7910
Group 10 455.16, 525.10, 607.46, 786.68, 835.11, 847.25, 859.42, 883.79, 932.749
Group 11 443.56, 501.72, 548.55, 619.30, 762.57, 798.77, 835.11, 847.25, 859.42, 871.60, 883.79, 945.0212
Group 12 443.56, 607.46, 619.30, 726.53, 774.62, 786.68, 822.98, 932.74, 957.329
Group 13 762.57, 810.86, 883.793
Table
Table 4. The results of feature wavelength bands selected using SD data under different adulteration groups.
Table 4. The results of feature wavelength bands selected using SD data under different adulteration groups.
GroupFeature Wavelength Bands [nm]Number
Group 1 478.41, 560.30, 619.30, 643.01, 654.89, 690.64, 726.53, 786.68, 822.98, 847.25, 896.0111
Group 2 595.65, 643.01, 654.89, 738.53, 835.11, 896.01, 945.027
Group 3 595.65, 643.01, 654.89, 726.53, 738.53, 835.11, 859.42, 945.028
Group 4 466.77, 643.01, 654.89, 726.53, 738.53, 774.62, 835.11, 847.25, 859.42, 945.0210
Group 5 631.35, 643.01, 654.89, 690.64, 726.53, 738.53, 786.68, 847.25, 896.019
Group 6 643.01, 654.89, 690.64, 738.53, 786.68, 847.25, 896.01, 945.028
Group 7 478.41, 560.30, 643.01, 654.89, 690.64, 702.58, 738.53, 786.68, 822.98, 847.25, 896.01, 945.0212
Group 8 643.01, 654.89, 690.64, 726.53, 738.53, 786.68, 822.98, 835.11, 847.25, 896.0110
Group 9 455.16, 560.30, 643.01, 654.89, 690.64, 786.68, 896.01, 945.028
Group 10 654.89, 690.64, 702.58, 822.98, 883.79, 896.016
Group 11 455.16, 490.06, 525.10, 607.46, 690.64, 750.54, 774.62, 798.77, 822.98, 847.25, 920.48, 932.7412
Group 12 443.56, 455.16, 583.85, 595.65, 643.01, 762.57, 774.62, 786.68, 859.42, 871.60, 920.48, 957.3212
Group 13 466.77, 490.06, 560.30, 643.01, 690.64, 714.55, 738.53, 750.54, 798.77, 810.8610
Table
Table 5. Comparison of recognition accuracy between LLRM and LRM using SNV data.
Table 5. Comparison of recognition accuracy between LLRM and LRM using SNV data.
Group 1Group 2Group 3Group 4Group 5Group 6Group 7
LLRM99.31%100%100%98.77%100%100%100%
LRM97.22%97.22%78.89%71.60%100%97.37%95.61%
Group 8Group 9Group 10Group 11Group 12Group 13Average
LLRM99.07%100%100%99.07%99.07%100%99.64%
LRM96.30%96.30%100%94.44%97.22%95.37%93.66%
Table
Table 6. Comparison of recognition accuracy between LLRM and LRM using FD data.
Table 6. Comparison of recognition accuracy between LLRM and LRM using FD data.
Group 1Group 2Group 3Group 4Group 5Group 6Group 7
LLRM98.61%100%100%100%100%100%99.12%
LRM95.83%86.11%73.33%74.07%87.72%90.35%94.74%
Group 8Group 9Group 10Group 11Group 12Group 13Average
LLRM100%100%99.07%95.37%97.22%100%99.18%
LRM92.59%87.96%97.22%88.89%84.26%99.07%88.63%
Table
Table 7. Comparison of recognition accuracy between LLRM and LRM using SD data.
Table 7. Comparison of recognition accuracy between LLRM and LRM using SD data.
Group 1Group 2Group 3Group 4Group 5Group 6Group 7
LLRM91.67%95.37%96.67%98.77%93.86%92.98%94.74%
LRM88.19%81.48%78.89%77.78%80.70%83.33%95.61%
Group 8Group 9Group 10Group 11Group 12Group 13Average
LLRM100%97.22%100%98.15%96.30%100%96.59%
LRM86.11%87.96%94.40%91.67%84.26%95.37%86.60%
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Liu, W.; Zeng, S.; Wu, G.; Li, H.; Chen, F. Rice Seed Purity Identification Technology Using Hyperspectral Image with LASSO Logistic Regression Model. Sensors 2021, 21, 4384. https://doi.org/10.3390/s21134384

AMA Style

Liu W, Zeng S, Wu G, Li H, Chen F. Rice Seed Purity Identification Technology Using Hyperspectral Image with LASSO Logistic Regression Model. Sensors. 2021; 21(13):4384. https://doi.org/10.3390/s21134384

Chicago/Turabian Style

Liu, Weihua, Shan Zeng, Guiju Wu, Hao Li, and Feifei Chen. 2021. "Rice Seed Purity Identification Technology Using Hyperspectral Image with LASSO Logistic Regression Model" Sensors 21, no. 13: 4384. https://doi.org/10.3390/s21134384

APA Style

Liu, W., Zeng, S., Wu, G., Li, H., & Chen, F. (2021). Rice Seed Purity Identification Technology Using Hyperspectral Image with LASSO Logistic Regression Model. Sensors, 21(13), 4384. https://doi.org/10.3390/s21134384

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