Next Article in Journal
Understanding the Development of Heterocyclic Aromatic Amines in Fried Bacon and in the Remaining Oil after Pan-Frying in Five Different Vegetable Oils
Next Article in Special Issue
Sustainable Environmental Assessment of Waste-to-Energy Practices: Co-Pyrolysis of Food Waste and Discarded Meal Boxes
Previous Article in Journal
Typical Characterization of Commercial Camellia Oil Products Using Different Processing Techniques: Triacylglycerol Profile, Bioactive Compounds, Oxidative Stability, Antioxidant Activity and Volatile Compounds
Previous Article in Special Issue
Butterfly Pea Flower as a Novel Ingredient to Produce Antioxidant-Enriched Yellow Pea-Based Breakfast Cereals
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

A Study of the Reliability and Accuracy of the Real-Time Detection of Forage Maize Quality Using a Home-Built Near-Infrared Spectrometer

1
China-Canada Joint Lab of Food Nutrition and Health (Beijing), School of Food and Health, Beijing Technology and Business University, 11 Fucheng Road, Beijing 100048, China
2
Biomass and Bioresource Utilization Laboratory, College of Engineering, China Agricultural University, Beijing 100083, China
*
Author to whom correspondence should be addressed.
Foods 2022, 11(21), 3490; https://doi.org/10.3390/foods11213490
Submission received: 12 August 2022 / Revised: 24 September 2022 / Accepted: 2 October 2022 / Published: 3 November 2022

Abstract

:
The current study was conducted to explore the real-time detection capability of a home-built grating-type near-infrared (NIR) spectroscopy online system to determine forage maize quality. The factor parameters affecting the online NIR spectrum collection were analyzed, and the results indicated that the detection optical path of 12 cm, conveyor speeds of 10 cm s−1, and number of scans of 32 were the optimal parameters. Choosing the crude protein and moisture of forage maize as quality indicators, the reliability of the home-built NIR online spectrometer was confirmed compared with other general research NIR instruments. In addition, an NIR online multivariate analysis model developed using the partial least squares (PLS) method for the prediction of forage maize quality was established, and the reliability, applicability, and stability of the NIR model were further discussed. The results illustrated that the home-built grating-type NIR online system performed satisfying and comparable accuracy and repeatability of the real-time prediction of forage maize quality.

1. Introduction

Maize is an abundant and important crop in the world and is widely consumed as an important animal feed source. The utilization of maize in compound feed is 50–70%, and the total consumption of forage maize will continue to grow [1,2]. In addition, global demand for livestock products has continuously increased in recent decades, as the consumption of meat has increased in many countries [3,4]. It is noteworthy that the quality of forage maize directly affects livestock health and animal products such as meat, eggs, and milk. Therefore, the control of forage maize quality is of importance to ensure animal product safety and even human health [5].
Near-infrared (NIR) is the region in which absorptions corresponding to overtones and combinations of the fundamental vibrational transitions occur [6,7]. In recent decades, with the development of new NIR technologies, NIR online analytical technology has been used in agriculture, food, medicine, tobacco, and other industries as a rapid, efficient, and effective tool for the process of real-time analysis [6,8,9,10,11,12,13,14]. The advantages of NIR spectroscopy over traditional methods are speed, accuracy and precision, no sample destruction, simple sample preparation, the ability to provide simultaneous information on several properties and constituents, relatively reduced operational costs, and no use of toxic reagents [8,15]. Combined with the fiber optic probe, non-invasive NIR online measurements can be carried out from the production of raw materials to the quality control in production processes, as well as the final quality analysis of products [16].
However, most feed quality detection in practical production have mainly relied upon the laboratory NIR spectrometer, which cannot satisfactorily guide the real-time production process. Therefore, for industrial applications, realizing online detection in real-time is essential for the rapid development of the modern feed industry, which provides the ability to detect the raw material, process control, and end product at the same time. Online detection also improves the analysis efficiency and prompt feedback for process control [17].
Actually, the use of the NIR technique for online process monitoring is still limited. The stability of the spectral collection is influenced by the sample state or the transport environment in the actual production process [18]. In addition, the online spectra acquisition compared with the stationary measurement may result in a relatively low signal-to-noise ratio of spectra due to the high noise generated by the mechanical vibrations of the collection system [19]. Consequently, establishing a reliable and stable NIR online analysis model for the whole process of online detection and forage maize quality control is challenging, with important practical significance. Moreover, for industrial applications, it is also essential to choose the adequate operation parameters and chemometric tools to obtain a robust model. Herein, partial least squares (PLS) was chosen to establish a model between component concentrations or properties of samples and their spectral absorbance measured at different wavelengths because of the simplicity of use, speed, and good performance [20].
In the current study, we focused on the real-time detection of forage maize quality based on a home-built grating-type NIR online system. The factor parameters (sample surface temperature, number of scans, detection optical paths, particle sizes, conveyor speeds) of the online NIR spectra acquisition were analyzed and optimized. The other two general research NIR spectrometers were chosen for the study of the reliability of the online NIR spectrometer. In addition, an NIR online multivariate analysis model for the prediction of forage maize quality (crude protein, moisture, crude ash, and energy) was established, and its reliability, applicability, and stability were analyzed and validated.

2. Materials and Methods

2.1. Sample Preparation

The experimental samples (138) of forage maize were collected from 18 provinces (Shandong, Hebei, Shanxi, Shaanxi, Sichuan, Yunnan, etc.) in China. The first batch was 108 forage maize samples for the parameter optimization and the reliability analysis of the established online model (75 samples for calibration and 33 samples for validation). The second batch was 30 samples selected from the other 110 forage maize using the Mahalanobis distance method. The raw samples passed the 8 mm standard sieve (Highway Instrument, Zhejiang, China), and the impurities such as sand and soil were removed. Then the samples were packaged in zip-lock bags and stored in a sample storage room at 20 ± 2 °C for use.

2.2. Experimental Analysis of Forage Maize Quality

Crude protein, moisture, crude ash, and energy were selected as the quality index of forage maize. The crude protein was detected according to GB/T 24901-2010 using an automatic Kay nitrogen analyzer (Kjeltec 2300, FOSS, Hillerød, Denmark) [21]. The moisture was measured according to GB/T 24900-2010 using an electric heating constant temperature drying oven (DHG-9123A, Shanghai Jinghong Laboratory Instrument Co., Ltd., Shanghai, China) [22]. The crude ash was tested based on GB/T 6438-2007 using a program-controlled chamber electric furnace (Shanghai Jinghong Laboratory Instrument Co., Ltd., Shanghai, China) [23]. The energy was measured using an automated oxygen bomb calorimeter (Parr 6300, PARR, Moline, IL, USA).

2.3. NIR Spectroscopy Systems

2.3.1. Different NIR Spectrometers

The main NIR spectrometer was the home-built grating online analysis system (Figure S1A) composed of a light source (two tungsten lamps), an optical splitting system (RockSolid), a detector (InGaAs), a transport component, and a sample attachment, denoted as Instrument A. The other two general research NIR spectrometers (Figure S1B,C) were as follows: Instrument B was a SupNIR-2700 grating NIR spectrometer (Focused Photonics Inc., Hangzhou, China, instrument B) with an integrated design including a light source system, a spectral acquisition component, and a transportation component. The light source system mainly included a fiber-coupled sensor and a light source (single tungsten lamp); the spectral acquisition components were mainly composed of an interferometer (RockSolid), a detector (InGaAs), and a circuit system; the transport component was a rotating sample stage. Instrument C was a Spectrum 400 diffuse reflection Fourier transform NIR spectrometer (Perkin Elmer, Waltham, MA, USA, instrument C) with an integrated design was measured using the integrated sphere diffuse reflection method. It was equipped with an InGaAs detector and a sealed anti-vibration optical platform protecting the optical system from moisture. The transport unit was a rotating sample stage.
In addition, the spectrum acquisition software of instrument A and B was from Focused Photonics Inc., while the spectra acquired software of instrument C was from Perkin Elmer. The experimental forage maize samples were equilibrated at room temperature for 72 h before NIR spectra collection. All three NIR instruments adopted non-contact diffuse reflection.

2.3.2. Spectra Acquisition

(1) Instrument A spectra collection: the sample was placed in the cuboid sample box (100 × 8 × 2 cm) on the conveyor belt, the optical detection path was adjusted to 12 cm, and the conveyor speed was set to 10 cm s−1. The NIR spectra were collected when the sample box was moved under the acquisition window. The sample was repeatedly loaded 3 times, and the average spectra were taken as the sample representative spectra. Spectrometer scanning parameters were as follows: the wavelength range was 1000–2500 nm, the spectral resolution was 10 nm, and the number of scans was 32.
(2) Instrument B spectra acquisition: the transport unit was a rotating sample stage. The sample was repeatedly loaded 3 times, and the average spectra were taken as the sample representative spectra. The spectrometer scanning parameters were as follows: the wavelength range was 1000–2500 nm, the spectral resolution was 10 nm, and the number of scans was 32.
(3) Instrument C spectra acquisition: the transport unit was a rotating sample stage. The sample was repeatedly loaded 3 times, and the average spectra were taken as the sample spectra. Spectrometer scanning parameters were as follows: the wavenumber range was 10,000–4000 cm−1, the spectral resolution was 8 cm−1, and the number of scans was 32.

2.4. Analysis of Influence Factors for NIR Online Spectrometer

2.4.1. Sample Surface Temperature

The halogen tungsten lamp of the grating NIR online spectrometer used in this study generates a large amount of heat at the focus of the light, which rapidly increases the surface temperature of the sample. To study the variation of sample surface temperature and the corresponding NIR spectra during the scanning process, the spectra of forage maize samples were measured by the NIR online spectrometer, and the real-time surface temperatures of samples were detected by the MC infrared thermometer (Zhixianglingyu Technology Co., Ltd., Shenzhen, China). For each forage maize sample, eight spectra were measured at the surface temperatures of 25, 30, 35, 40, 45, 50, 55, and 60 °C. The wavelength range was 1000–2500 nm, the number of scans was 32, the detection optical path was 12 cm, and the conveyor speed was 0 cm s−1. The differences of average spectra of samples collected at different temperatures were compared and analyzed to study the variation of NIR spectra of forage maize with different surface temperatures.

2.4.2. Different Number of Spectra Scans

Different scanning times cause variations of the sample surface temperature, which results in spectral data variation. The scanning time was decided by the different number of scans. Therefore, the number of scans was set to 1, 16, 32, and 64 to study the influence of the number of scans on the NIR spectra of the forage maize. The forage maize samples at different numbers of scans were collected and repeatedly tested 10 times.
The evaluation index of spectral reproducibility is expressed by the standard deviation (SD), and the calculation formula is as follows:
  S D l = i = 1 10 ( A i l A a v l ) 10 1
where A i l is the absorbance of l wavelength of the i spectrum, and A a v l is the average absorbance of l wavelength.
The spectral reproducibility under different numbers of scans was calculated by the average SD of the representative samples. The wavelength range was 1000–2500 nm, the detection optical path was 12 cm, the environment temperature was 25 °C, and the conveyor speed was 0 cm s−1. The effects of different numbers of scans on the spectral reproducibility of forage maize were studied by analyzing the changes of spectral SD.

2.4.3. Different Detection Optical Paths

The detection optical path was defined as the distance between the spectrometer scanning window and the surface of the sample. Five different detection optical paths (8, 10, 12, 14, and 16 cm) were set to evaluate the effects of different detection optical paths on NIR spectra by adding different spacers under the sample cup. The wavelength range was 1000–2500 nm, the scanning times was 32, the environment temperature was 25 °C, and the conveyor speed was 0 cm s−1. Using crude protein as the quality index, the quantitative analysis models of different detection optical paths were established to study the effects of different detection optical paths on the NIR spectra of forage maize.

2.4.4. Different Particle Sizes

NIR spectra of raw forage maize and two other different particle size samples (1 mm and 0.425 mm) were prepared to explore the effects of different particle sizes on NIR spectra. The wavelength range was 1000–2500 nm, the detection optical path was 12 cm, the number of scans was 32, the environment temperature was 25 °C, and the conveyor speed was 0 cm s−1. The quantitative analysis model of crude protein with different particle sizes for forage maize was established to study the effects of different particle sizes on NIR spectra.

2.4.5. Different Conveyor Speeds

The NIR spectra of forage maize samples were collected at different conveyor speeds (0, 5, 8, 10, and 20 cm s−1). The wavelength range was 1000–2500 nm, the detection optical path was 12 cm, the scanning times was 32, and the environment temperature was 25 °C. The quantitative analysis model of crude protein for forage maize was established at different conveyor speeds to analyze the effects of different conveyor speeds on NIR spectra.

2.4.6. Orthogonal Experiment of Influence Factors of NIR Online Analysis

Based on the above-mentioned single factor experiment, the effects of the detection optical path, the conveyor speed, and the number of scans during the dynamic spectral data acquisition were further investigated using three-factor and three-level orthogonal experiments. The detection optical paths were 8, 12, and 16 cm; the conveyor speeds were 10, 20, and 30 cm s−1; and the numbers of scans were 8, 16, and 32 times. The orthogonal experiment of the influence factors of the NIR online analysis are shown in Table S1.
The root mean square error (RMS) was used to evaluate the reproducibility of the spectrum and was calculated as follows:
R M S l = i = 1 n ( A i l A ¯ i ) 2 n
where A i l is the absorbance of i wavelength of the l spectrum, A ¯ i is the average absorbance of i wavelength, and n is the total wavelength point.

2.5. Data Processing and Analysis

An NIR quantitative analysis model was established using the PLS method in the PLS-Toolbox 8.0 based on MATLAB. Spectral preprocessing was performed by standard normal variate (SNV), detrend, and first derivative processing [10]. The NIR spectra of the forage maize were divided into the calibration set and validation set based on the concentration gradient method. Outlier samples were rejected by the studentized residual of the sample chemical value and the Mahalanobis distance [15]. The studentized residual is the ratio of the difference between the predicted value and the measured value of the sample to the standard error of the sample. The equation is calculated as follows.
  Y   S t d n t   R e s i d u a l i = y ˜ i y i S E E
where y ˜ i is the model predicted value of the i sample, y i is the true chemical value of the i sample, and SEE is the predicted standard error of all samples.
The Mahalanobis distance of a spectrum refers to the distance between the sample spectrum and the average spectrum of the entire spectrum set. After normalizing the spectral data, the Mahalanobis distance of a single sample is determined by the following equation:
h i i = t i T ( T i T ) 1 t i
where  h i i represents the impact of a single sample on the entire standard sample set. In the near-infrared spectroscopy analysis, h i i expresses the degree of influence of sample i on the regression model, t i is the score vector of sample i, and T i T is the score matrix of the modeling set. The grating NIR online spectrometer and two general research NIR spectrometers were chosen to compare and investigate the reliability of the NIR online analysis method for the forage maize quality. The correlation coefficient of calibration (Rc), the correlation coefficient of prediction (Rp), the root mean square error of calibration (RMSEC), and the root mean square error of prediction (RMSEP) were chosen to evaluate the obtained model and were expressed as follows [24].
R c = [ i = 1 n c ( y ^ i y ¯ c ) i = 1 n c ( y i y ¯ c ) ]
R p = [ i = 1 n v ( y ^ i y ¯ v ) i = 1 n v ( y i y ¯ v ) ]
R M S E C = 1 n c i = 1 n c ( y i y ^ i ) 2
R M S E P = 1 n v i = 1 n v ( y i y ^ i ) 2
where y i is the measured chemical value of the i sample, y ^ i is the NIR absorbance value of the  i sample, n c is the number of samples in the calibration set, and n v is the number of samples in the validation set; y ¯ c and y ¯ v are the average chemical values of calibration set and the validation set, respectively.
The closer the Rc and Rp values are to 1, the closer the predicted value is to the true value. The smaller the RMSEC and RMSEP and the smaller the difference between them, the better the model performance.
Relative analysis error (RPD) is calculated by:
RPD = SD/RMSEP
where SD is the standard deviation of the chemical values of the validation set.
The smaller the difference in model results, the more reliable the NIR online analysis.

3. Results and Discussion

3.1. Parameter Analysis of NIR Spectra Collection

3.1.1. Variation of Surface Temperature and Spectral Absorption during Sample Scanning

Figure 1a depicts the variation of surface temperature of the forage maize sample against the scanning times (8 min and 1 min). The surface temperature increased continuously with the increase of scanning time. In the first minute, the rate of increase was rapid, and the fit curve of temperature against time was studied, which was Y = 0.27X + 25.86, while as time passed, the rising speed showed a decreased trend.
The NIR spectra of the forage maize samples at different temperatures are presented in Figure 1b. The spectral absorbance was the highest at 25 °C and decreased with increasing temperature. The absorption peaks at 1400–1500 nm and 1850–1950 nm corresponding to the moisture peaks varied in the intensities and positions with the changing temperature [18,25]. The temperature variation resulted in the force change between the molecules, which led to the change of vibration and rotation of the O–H, resulting in the NIR spectral shift [26]. According to a previous report [26], temperature affected the spectra in a nonlinear form. As the temperature increased, the molecular energy increased, leading to the molecular transition to higher energy levels, which contributed to the increasing absorption intensity, and the decreasing reflection intensity [27]. Consequently, the absorptions of NIR light energy of the sample molecules were weakened, and the absorption intensities were reduced. Therefore, the environment temperature and sample scanning time must be strictly controlled when collecting NIR spectral data.

3.1.2. Effect of Different Number of Scans on NIR Spectral

NIR spectra and SD values of the forage maize samples with different numbers of scans are presented in Figure 1c, as well as the NIR spectra (SD values) with SNV and detrend preprocessing (Figure 1d). For the original spectra, the spectral noise decreased as the number of scans increased. Comparing the NIR spectra with different numbers of scans, the SD value was the highest under the number of scans of 1 and decreased until the number of scans increased to 32 and then gradually stabilized. The NIR spectra were further preprocessed by SNV and detrended; however, the spectral noise cancellation was not obvious for the spectra with different numbers of scans.
Based on the above mentioned, increasing the number of scans increased the signal-to-noise ratio; however, the increase of number of scans led to the extended scan time and the pronounced change in the sample surface temperature, which thereby caused the NIR spectra shift [18]. Therefore, the reasonable selection of the number of spectra scans can effectively reduce the influence of temperature change on the spectra, as well as improve the detection efficiency in the online process. According to the results, the number of scans of 32 was considered the reasonable parameter at which the spectral data were relatively stable, and the detection efficiency was high.

3.1.3. Effect of Different Detection Optical Paths on NIR Spectra

Figure 1e presents the NIR spectra of forage maize with different detection optical paths in which the spectra with detection paths of 8 cm and 10 cm were overlapped, while the other spectra showed obvious differences with different detection paths. The spectral absorbance increased with the increase of the detection optical path. The NIR spectra with different detection optical paths were further analyzed by principal component analysis (PCA) and are shown in Figure 1f. Most samples at the detection optical path of 8 cm and 10 cm were overlapped, and they partly overlapped with the sample of the 12 cm detection optical path, while the samples with the detection optical path of 12 cm and 16 cm presented the obvious difference, which indicated the difference among the forage maize spectra that presented an increasing trend with the increase of the detection optical path.
Table S2 lists the results of the quantitative analysis model for crude protein of forage maize at different detection optical paths. The selected optimal preprocessing method was the SNV, detrend, and first derivative processing [20]. The results revealed that the NIR quantitative model of crude protein at the optical path of 12 cm was ideal (RPD = 2.32), in which the coefficient of calibration and prediction were both larger than those in other detection optical paths. As the detection optical path decreased or increased, the prediction effect of the model was reduced.

3.1.4. Effect of Different Particle Sizes on NIR Spectra

The NIR spectra of forage maize with different particle sizes are depicted in Figure 1g, which revealed there are distinct differences between the raw forage maize and the crushed samples (1 mm and 0.425 mm). The spectral absorbance increased with the increase of the particle size. This is because the change of particle size affects the light penetration and reflection characteristics, which contributes to the light propagation in the sample variation to some extent and influences the scattering and absorption of light by the sample. In addition, the scattering coefficient of the sample decreased as the particle size increased; thus, the larger the particle size, the higher the absorbance of the overall spectra [19]. The PCA results of NIR spectra with different particle sizes (Figure 1h) illustrated that the principal component score of different particle size presented differences. The smaller the difference in the particle size, the smaller the difference in the performance of the principal component score [27].
Table S2 shows the results of the quantitative analysis model for crude protein of forage maize in different particle sizes. The Rc of the crude protein analysis model established by NIR online spectra were all greater than 0.94, and the RPD were higher than 2.32. The prediction of the quantitative model for 0.425 mm crushed maize was better than the predictions for the raw maize and 1 mm crushed maize.

3.1.5. Effect of Different Conveyor Speeds on NIR Spectra

The results of the quantitative analysis model for crude protein of forage maize at different conveyor speeds are displayed in Table S2. The crude protein model at 5 cm s−1 was better than those at 10 and 20 cm s−1, of which Rc and Rp were 0.94 and 0.86, respectively, RMSEC and RMSEP were 0.22 and 0.31, respectively, and RPD was 1.87. With the increase of conveyor speed, the performance of the prediction model for crude protein was decreased but not distinct.
Treating the spectra at 0 cm s−1 as the standard spectra, the spectral residuals between the spectra at 5, 10, and 20 cm s−1 and the standard spectra were considered as noise. Figure 1i–l depicts the residual distribution and residual probability density distribution without and with SNV + detrend + first derivative preprocessing. According to the residual distribution, noises at different wavelengths were different. However, SNV + detrend + first derivative preprocessing significantly reduced the noise [11,20]. Residual probability density distribution after spectral preprocessing also conformed to the Gaussian distribution; therefore, spectral preprocessing could effectively reduce the influence of speed noise.
Table S3 shows the signal-to-noise ratios at different conveyor speeds. The signal-to-noise ratios of water, crude protein, and crude fat peaks gradually decreased with increasing speed. This was consistent with the results in Table S2 showing that the quantitative model for crude protein of forage maize declined as the speed increased.
The change rate of the spectral area is also the evaluation index of spectral stability [14]. Table S3 presents the change rate of the spectral area under different conveyor speeds, in which the change rate at the condition of 5 cm s−1 was the smallest (0.003). As the speed increased, the change rate of the spectral area increased continuously. If there is no noise in the spectra, the spectra collected under the same conditions should be coincident [18]. That is, the smaller the change rate of the spectral area, the more stable the spectral detection system. Therefore, as the conveyor speed increased, the spectral noise increased, and the result of the quantitative analysis model for crude protein of forage maize continued to decline.

3.1.6. Parameter Optimization of Influence Factors of NIR Online Analysis

The parameter combinations of orthogonal experiments for the different influence factors of NIR online analysis are listed in Table S4, which also reveals the different RMS values. The results of multivariate analysis of variance are presented in Table S4. The number of spectra scans and the detection optical path were both significant at the significance level of 0.05 (F > F0.95(2,2) = 19), while the conveyor speed was not significant (F = 0.633 < F0.95(2,2) = 19). Furthermore, the significance of the number of spectra scans was greater than that of the detected optical path (107.985 > 46.848).
Based on the above mentioned, when the detection optical path was 12 cm, the conveyor speed was 10 cm s−1, and the number of spectra scans was 32, and the RMS value was the smallest, that is, the spectral reproducibility was the best, which can be used as the optimal parameter reference for NIR online dynamic spectral data acquisition [14].

3.2. Establishment and Comparison of Forage Maize Model Using Three Different NIR Spectrometers

3.2.1. Component Analysis of Forage Maize

The statistical results of crude protein and moisture content of experimental forage maize are presented in Table 1, the averages of which were 8.01% and 12.18%. The crude protein content ranged from 5.44% to 10.76%, and the variance was 0.79; the moisture content ranged from 7.16% to 15.87%, and the variance was 1.84. According to the distribution of crude protein and moisture content (data not shown), protein and moisture both followed a normal distribution [28].

3.2.2. Spectral Characteristics of Different NIR Spectrometers

The spectra of forage maize collected from three different NIR spectrometers are depicted in Figure 2, which reveals there were some differences in the spectral absorbance among the three instruments. The frequency multiplication and frequency combination of the atom vibration in the molecule mostly occurred in the NIR region. As seen in Figure 2, the NIR original spectra of maize had three distinct characteristic peaks, namely, 1196 nm, 1450 nm, and 1930 nm. The peaks at 1196 and 1450 nm corresponded to the second-order vibration of the O–H bond, and the characteristic peak at 1930 nm corresponded to the first-order vibration of the O–H bond, all of which were mostly related to the moisture content of maize [7].
Figure 2d displays the NIR average spectra (second derivative) of forage maize collected by three different spectrometers. The peak around 1700–1735 nm was mainly the C–H first-order vibration, which is the lipid absorption caused by the structure of CH2 and CH3 in maize oil; the peak at 2248 nm was the combined vibration of C–H and N–H, mainly due to the absorption of crude protein; the peaks at about 2330 nm and 2349 nm were C–H stretching vibration and deformation vibration, which are related to the fiber composition of maize [18,29].

3.2.3. Comparative Analysis of Models of Different NIR Spectrometers

The crude protein and water models of raw forage maize for different NIR spectrometers were established using PLS combined with spectral preprocessing (SNV, detrend, and first derivative). The model parameters and the scatter plots are presented in Figure 3. Comparing the optimal quantitative models of crude protein established by three NIR spectrometers, the Rp and RMSEP of instrument A, B, and C were 0.78 and 0.40, 0.86 and 0.31, and 0.86 and 0.36, respectively. The crude protein results illustrated that the instrument B and C performed slightly better than instrument A, but the difference was not significant (data not shown). As for the moisture, the models of the three instruments were similar.
The SupNIR-2700 grating NIR spectrometer (instrument B) and Spectrum 400 diffuse reflection Fourier transform NIR spectrometer (instrument C) are market-oriented products, and the accuracy and stability have been widely recognized in different industries. Comparing the three NIR spectrometers, the quantitative analysis model of forage maize established by the NIR online spectrometer (instrument A) was basically the same as the other two general research NIR spectrometers, indicating that the prediction of the NIR online analysis spectrometer was reliable and meets the needs of practical applications.
As seen in Figure 3, the crude protein model was better than the moisture model, which is because the moisture of maize particle did not appear evenly distributed. From the outer to the inner of the grain, the moisture content was reduced in turn, and the moisture content was also easily affected by the external environment or the temperature compared with the crude protein [30]. Therefore, the crude protein spectral information was superior to the moisture spectral information when acquiring the NIR spectra.

3.3. Validation Analysis of NIR Online System for Forage Maize Quality Detection

3.3.1. Component Statistics of Forage Maize

Table 1 lists the contents of crude protein, moisture, crude ash, and energy of raw and crushed (0.425 mm) forage maize. Some outliers have been removed due to the abnormal values. In raw samples, 5, 6, 4, and 5 samples were rejected for crude protein, moisture, crude ash, and energy, respectively, while in the crushed sample, 0, 7, 4, and 0 samples were respectively removed. As seen in Table 1, the component contents between the raw and crushed forage maize were slightly different, which may be caused by the mechanical milling process.

3.3.2. Establishment of NIR Online Model for Forage Maize Quality

The NIR online multivariate analysis model of raw and crushed forage maize was established by PLS combined with spectral processing (SNV, detrend, and first derivative). The model parameters and scatter plots are depicted in Figure 4. The Rc and Rp of models for the crude protein, moisture, crude ash, and energy of the raw forage maize were all lower than those of the crushed forage maize, while the RMSEC and RMSEP of model for the raw maize were greater than those of the crushed maize, which illustrated that the prediction of the model for the crushed forage maize performed better than that for the raw forage maize, which was also in accordance with the results from the quantitative analysis model of crude protein of forage maize with different particle sizes [31]. Comparing the model results of the crude protein, moisture, crude ash, and energy, the crude protein model was the best, the moisture and energy model were slightly worse, and the crude ash model was not satisfied, mainly because of the low crude ash content and the uneven distribution in the forage maize.

3.4. Verification for Practical Application of NIR Online Model for Forage Maize Quality

3.4.1. Repeatability Analysis of NIR Online Model for Forage Maize

Based on the dynamic online analysis system, the representative 10 forage maize samples were continuously and repeatedly recorded 20 times to obtain the mean values of crude protein, moisture, crude ash, and energy. The plots of prediction results over time for the raw and crushed forage maize were established and are presented in Figure S2. According to the real-time online spectral data within 4 min, the dynamic prediction results of crude protein, moisture, crude ash, and energy of forage maize fluctuated around the true values of 7.98, 12.22, 1.10, and 16.72, respectively, indicating that the prediction ability of the model had good repeatability [28,31].
In addition, the fluctuation of dynamic prediction results for the crude protein, moisture, crude ash, and energy of crushed forage maize were significantly less obvious than those of raw forage maize, which illustrated the prediction reproducibility of the crushed maize model was much better than that of the raw maize model. This might be because the crushed sample was more homogenous compared with the raw sample, which resulted in a more stable spectra being generated by the crushed sample.

3.4.2. Accuracy Analysis of NIR Online Model for Forage Maize Quality

Table 2 shows the comparison between the true values of the forage maizes (raw and crushed sample) and their NIR predicted values. According to the statistical result of the representative raw sample, the relative errors of most crude proteins were around 0.3%; the relative errors of all moisture were below 3%; the relative errors of crude ash were all greater than 15%, mainly due to the low crude ash content and non-uniform distribution in forage maize; furthermore, the relative errors of energy were all lower than 2.5%.
From the statistical results of the representative crushed samples, the relative error of crude protein and moisture were all less than 2.5%; the relative error of energy was lower than 1.1%, all of which were satisfied with the relative error requirement. However, the relative error of crude ash was not as good as the other components. Based on the statistical analysis of repeatability and accuracy of model prediction, NIR provides an effective data reference for the quality control of forage maize in practical production process [28].

4. Conclusions

In the current study, a home-built grating online analysis system was compared with the other two general research NIR spectrometers. The number of spectra scans, particle size, detection optical path, and conveyor speed were analyzed to investigate the influence factors and parameter optimization. The results indicated that a detection optical path of 12 cm, a conveyor speed of 10 cm s−1, and number of scans of 32 were the optimal parameters. Moreover, the decrease of particle size effectively improved the model performance. Choosing the crude protein and moisture of forage maize as quality indicators, the home-built grating online analysis system (instrument A) and two general research NIR spectrometers (grating instrument B, Fourier transform instrument C) were to establish the quantitative analysis model of forage maize. According to the characteristics of spectral data and the results of the quantitative analysis model, the reliability of the home-built NIR online analysis system was confirmed. Based on the home-built grating NIR online spectrometer, NIR online multivariate analysis models for the raw and crushed forage maize were established. The correlation coefficients of validation for crude protein, moisture, crude ash, and energy of the crushed forage maize were slightly greater than those of the raw maize, both of which were proven to be suitable for online practical application. Online dynamic verification illustrated that the model prediction of the home-built NIR online analysis system had good accuracy and repeatability.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/foods11213490/s1, Figure S1: Three different types of near infrared spectrometers; Figure S2: Repeatability of NIR online model of raw and crushed forage maize for crude protein, moisture, crude ash, and energy (a–d, raw forage maize; e–h, crushed forage maize); Table S1: Orthogonal experimental results of influence factors of NIR online analysis; Table S2: Quantitative model results of crude protein of forage maize at different detection optical path, particle size and conveyor speed; Table S3: Signal-to-noise ratio and change rate of peak area at different conveyor speed; Table S4: Result of multivariate analysis of variance.

Author Contributions

Conceptualization, F.G. and X.L.; methodology, F.G.; software, Y.Z.; validation, Y.Z.; formal analysis, Y.Z.; investigation, Y.Z.; resources, X.L.; data curation, Y.Z.; writing—original draft preparation, F.G.; writing—review and editing, F.G.; visualization, F.G.; supervision, X.L.; project administration, F.G.; funding acquisition, X.L. 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 (32101635), the Research Foundation for Youth Scholars of Beijing Technology and Business University, grant number (19008022138), and the National Key R&D Program Intergovernmental/Hong Kong, Macao, and Taiwan key projects, grant number (2019YFE0103800).

Data Availability Statement

The data presented in this study are available on request from the corresponding author.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Rudel, T.K.; Paul, B.; White, D.; Rao, I.M.; Van Der Hoek, R.; Castro, A.; Boval, M.; Lerner, A.; Schneider, L.; Peters, M. LivestockPlus: Forages, sustainable intensification, and food security in the tropics. Ambio 2015, 44, 685–693. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  2. Williams, P.J.; Geladi, P.; Britz, T.J.; Manley, M. Investigation of fungal development in maize kernels using NIR hyperspectral imaging and multivariate data analysis. J. Cereal Sci. 2012, 55, 272–278. [Google Scholar] [CrossRef]
  3. Saha, U.; Endale, D.; Tillman, P.G.; Johnson, W.C.; Gaskin, J.; Sonon, L.; Schomberg, H.; Yang, Y. Analysis of various quality attributes of sunflower and soybean plants by near infrared reflectance spectroscopy: Development and validation calibration models. Am. J. Anal. Chem. 2017, 8, 462. [Google Scholar] [CrossRef] [Green Version]
  4. Godfray, H.C.J.; Aveyard, P.; Garnett, T.; Hall, J.W.; Key, T.J.; Lorimer, J.; Pierrehumbert, R.T.; Scarborough, P.; Springmann, M.; Jebb, S.A. Meat consumption, health, and the environment. Science 2018, 361, eaam5324. [Google Scholar] [CrossRef] [Green Version]
  5. Gao, F.; Zhou, S.; Han, L.; Yang, Z.; Liu, X. A novel FT-IR spectroscopic method based on lipid characteristics for qualitative and quantitative analysis of animal-derived feedstuff adulterated with ruminant ingredients. Food Chem. 2017, 237, 342–349. [Google Scholar] [CrossRef]
  6. Pullanagari, R.R.; Yule, I.J.; Agnew, M. On-line prediction of lamb fatty acid composition by visible near infrared spectroscopy. Meat Sci. 2015, 100, 156–163. [Google Scholar] [CrossRef]
  7. Ventura, M.; Silva, J.; Andrade, L.; Júnior, R.S.; Lima, S. Near-near-infrared thermal lens spectroscopy to assess overtones and combination bands of sulfentrazone pesticide. Spectrochim. Acta Part A Mol. Biomol. Spectrosc. 2018, 188, 32–36. [Google Scholar] [CrossRef]
  8. Salguero-Chaparro, L.; Baeten, V.; Abbas, O.; Peña-Rodríguez, F. On-line analysis of intact olive fruits by vis–NIR spectroscopy: Optimisation of the acquisition parameters. J. Food Eng. 2012, 112, 152–157. [Google Scholar] [CrossRef]
  9. de Lima, S.M.; Silva, B.F.A.; Pontes, D.V.; Pereira, C.F.; Stragevitch, L.; Pimentel, M.F. In-line monitoring of the transesterification reactions for biodiesel production using NIR spectroscopy. Fuel 2014, 115, 46–53. [Google Scholar] [CrossRef] [Green Version]
  10. Grassi, S.; Alamprese, C. Advances in NIR spectroscopy applied to process analytical technology in food industries. Curr. Opin. Food Sci. 2018, 22, 17–21. [Google Scholar] [CrossRef]
  11. Omar, J.; Slowikowski, B.; Boix, A. Chemometric approach for discriminating tobacco trademarks by near infrared spectroscopy. Forensic Sci. Int. 2019, 294, 15–20. [Google Scholar] [CrossRef] [PubMed]
  12. Liu, R.; Sun, Q.; Hu, T.; Li, L.; Nie, L.; Wang, J.; Zhou, W.; Zang, H. Multi-parameters monitoring during traditional Chinese medicine concentration process with near infrared spectroscopy and chemometrics. Spectrochim. Acta Part A Mol. Biomol. Spectrosc. 2018, 192, 75–81. [Google Scholar] [CrossRef] [PubMed]
  13. Stockl, A.; Lichti, F. Near-infrared spectroscopy (NIRS) for a real time monitoring of the biogas process. Bioresour. Technol. 2018, 247, 1249–1252. [Google Scholar] [CrossRef] [PubMed]
  14. Zhu, M.T.; Long, Y.; Chen, Y.; Huang, Y.S.; Tang, L.J.; Gan, B.; Yu, Q.; Xie, J.H. Fast determination of lipid and protein content in green coffee beans from different origins using NIR spectroscopy and chemometrics. J. Food Compos. Anal. 2021, 102, 104055. [Google Scholar] [CrossRef]
  15. Porep, J.U.; Kammerer, D.R.; Carle, R. On-line application of near infrared (NIR) spectroscopy in food production. Trends Food Sci. Technol. 2015, 46, 211–230. [Google Scholar] [CrossRef]
  16. Alcalà, M.; Blanco, M.; Bautista, M.; González, J.M. On-line monitoring of a granulation process by NIR spectroscopy. J. Pharm. Sci. 2010, 99, 336–345. [Google Scholar] [CrossRef]
  17. Xue, J.; Yang, Z.; Han, L.; Liu, Y.; Liu, Y.; Zhou, C. On-line measurement of proximates and lignocellulose components of corn stover using NIRS. Appl. Energy 2015, 137, 18–25. [Google Scholar] [CrossRef]
  18. Sampaio, P.S.; Soares, A.; Castanho, A.; Almeida, A.S.; Oliveira, J.; Brites, C. Optimization of rice amylose determination by NIR-spectroscopy using PLS chemometrics algorithms. Food Chem. 2018, 242, 196–204. [Google Scholar] [CrossRef]
  19. Salguero-Chaparro, L.; Peña-Rodríguez, F. On-line versus off-line NIRS analysis of intact olives. LWT-Food Sci. Technol. 2014, 56, 363–369. [Google Scholar] [CrossRef]
  20. Thissen, U.; Pepers, M.; Üstün, B.; Melssen, W.; Buydens, L. Comparing support vector machines to PLS for spectral regression applications. Chemom. Intell. Lab. Syst. 2004, 73, 169–179. [Google Scholar] [CrossRef]
  21. GB/T 24901-2010; National Grain and Oil Standardization Technical Committee. Inspection of Grain and Oils-Etermination of Crude Protein in Maize-Near-Infrared Method. Standards Press of China: Beijing, China, 2010.
  22. GB/T 24900-2010; National Grain and Oil Standardization Technical Committee. Inspection of Grain and Oils-Determination of Moisture Content in Maize-Near-Infrared Method. Standards Press of China: Beijing, China, 2010.
  23. GB/T 6438-2007; National Feed Industry Standardization Technical Committee. Animal feeding stuffs-Determination of Crude Ash. Standards Press of China: Beijing, China, 2007.
  24. Morellos, A.; Pantazi, X.-E.; Moshou, D.; Alexandridis, T.; Whetton, R.; Tziotzios, G.; Wiebensohn, J.; Bill, R.; Mouazen, A.M. Machine learning based prediction of soil total nitrogen, organic carbon and moisture content by using VIS-NIR spectroscopy. Biosyst. Eng. 2016, 152, 104–116. [Google Scholar] [CrossRef] [Green Version]
  25. Pandiselvam, R.; Thirupathi, V.; Vennila, P. Fourier transform–Near infrared spectroscopy for rapid and nondestructive measurement of amylose content of paddy. Sci. J. Agric. Eng. 2016, 2, 93–100. [Google Scholar]
  26. Peirs, A.; Scheerlinck, N.; Nicolaï, B.M. Temperature compensation for near infrared reflectance measurement of apple fruit soluble solids contents. Postharvest Biol. Technol. 2003, 30, 233–248. [Google Scholar] [CrossRef]
  27. Huber, K.-P. Molecular Spectra and Molecular Structure: IV. Constants of Diatomic Molecules; Springer Science & Business Media: Berlin/Heidelberg, Germany, 2013. [Google Scholar]
  28. Hetta, M.; Mussadiq, Z.; Wallsten, J.; Halling, M.; Swensson, C.; Geladi, P. Prediction of nutritive values, morphology and agronomic characteristics in forage maize using two applications of NIRS spectrometry. Acta Agric. Scand. Sect. B Soil Plant Sci. 2017, 67, 326–333. [Google Scholar] [CrossRef]
  29. Basri, K.N.; Hussain, M.N.; Bakar, J.; Sharif, Z.; Khir, M.F.A.; Zoolfakar, A.S. Classification and quantification of palm oil adulteration via portable NIR spectroscopy. Spectrochim. Acta Part A Mol. Biomol. Spectrosc. 2017, 173, 335–342. [Google Scholar] [CrossRef] [PubMed]
  30. Hell, J.; Prückler, M.; Danner, L.; Henniges, U.; Apprich, S.; Rosenau, T.; Kneifel, W.; Böhmdorfer, S. A comparison between near-infrared (NIR) and mid-infrared (ATR-FTIR) spectroscopy for the multivariate determination of compositional properties in wheat bran samples. Food Control 2016, 60, 365–369. [Google Scholar] [CrossRef]
  31. Wang, Y.; Guo, W.; Zhu, X.; Liu, Q. Effect of homogenisation on detection of milk protein content based on NIR diffuse reflectance spectroscopy. Int. J. Food Sci. Technol. 2019, 54, 387–395. [Google Scholar] [CrossRef]
Figure 1. NIR results at different sample surface temperatures, number of scans, detection optical paths, particle sizes, and conveyor speeds ((a) the variation of surface temperature of the forage maize against the scanning time; (b) NIR spectra of the forage maize at different temperatures; (c,d) NIR spectra (SD values) of the forage maize with different number of scans without and with SNV and detrend preprocessing; (e,f) NIR spectra and PCA result of the forage maize with different detection optical path; (g,h) NIR spectra and PCA result of the forage maize with different particle sizes; (il) the residual distribution and residual probability density distribution without and with SNV + detrend + first derivative preprocessing).
Figure 1. NIR results at different sample surface temperatures, number of scans, detection optical paths, particle sizes, and conveyor speeds ((a) the variation of surface temperature of the forage maize against the scanning time; (b) NIR spectra of the forage maize at different temperatures; (c,d) NIR spectra (SD values) of the forage maize with different number of scans without and with SNV and detrend preprocessing; (e,f) NIR spectra and PCA result of the forage maize with different detection optical path; (g,h) NIR spectra and PCA result of the forage maize with different particle sizes; (il) the residual distribution and residual probability density distribution without and with SNV + detrend + first derivative preprocessing).
Foods 11 03490 g001
Figure 2. Spectra of forage maize collected from three different NIR spectrometers and their representative spectra of the second derivative ((a) instrument A; (b) instrument B; (c) instrument C; (d) NIR average spectra (second derivative) of forage maize).
Figure 2. Spectra of forage maize collected from three different NIR spectrometers and their representative spectra of the second derivative ((a) instrument A; (b) instrument B; (c) instrument C; (d) NIR average spectra (second derivative) of forage maize).
Foods 11 03490 g002
Figure 3. Scatter plots of crude protein and moisture models of forage maize for different NIR spectrometers ((a,c,e) crude protein model of instrument A, B and C; (b,d,f) moisture model of instrument A, B, and C).
Figure 3. Scatter plots of crude protein and moisture models of forage maize for different NIR spectrometers ((a,c,e) crude protein model of instrument A, B and C; (b,d,f) moisture model of instrument A, B, and C).
Foods 11 03490 g003
Figure 4. Scatter plots of NIR online analysis model for raw and crushed forage maize by the home-built spectrometer ((a,c,e,g) crude protein, moisture, crude ash, and energy models of raw forage maize; (b,d,f,h) crude protein, moisture, crude ash, and energy models of crushed forage maize).
Figure 4. Scatter plots of NIR online analysis model for raw and crushed forage maize by the home-built spectrometer ((a,c,e,g) crude protein, moisture, crude ash, and energy models of raw forage maize; (b,d,f,h) crude protein, moisture, crude ash, and energy models of crushed forage maize).
Foods 11 03490 g004
Table 1. Component statistics of raw and crushed forage maize samples.
Table 1. Component statistics of raw and crushed forage maize samples.
Sample CalibrationValidation
Maximum
Value (%)
Minimum
Value (%)
Average (%)VarianceMaximum
Value (%)
Minimum
Value (%)
Average (%)Variance
Raw sampleCrude
Protein
10.765.448.010.798.525.447.450.60
moisture15.877.1612.181.8416.4210.8813.241.86
Crude ash1.610.741.110.191.400.771.050.21
Energy17.4215.9416.610.3117.0815.9816.440.30
Crushed sampleCrude
Protein
10.765.447.950.789.255.447.670.94
moisture16.427.1612.251.8315.659.2712.901.63
Crude ash1.610.771.120.191.480.741.030.19
Energy17.4215.9416.590.3117.0115.9816.500.26
Table 2. Comparison of experimental values and NIR predictions of raw and crushed forage maize samples.
Table 2. Comparison of experimental values and NIR predictions of raw and crushed forage maize samples.
Crude Protein (%)Moisture (%)Crude Ash (%)Energy (MJ kg−1)
True ValuePredicted ValueRelative ErrorTrue ValuePredicted ValueRelative ErrorTrue ValuePredicted ValueRelative ErrorTrue ValuePredicted ValueRelative Error
Raw
Sample
18.438.331.189.039.302.991.040.8617.317.0216.860.94
28.097.882.5912.3712.652.261.180.9618.616.7416.382.15
37.887.643.0511.7511.502.131.180.9816.916.6616.530.72
47.877.662.6712.6212.381.901.080.8422.2216.7716.402.21
57.597.843.1612.9412.801.081.020.8615.6816.5916.570.12
67.777.552.8313.2612.942.411.010.7921.7816.5316.361.02
77.987.723.2613.2112.922.201.080.8819.4416.6416.361.68
87.898.123.0412.8112.542.031.100.9315.4516.7416.332.45
97.747.542.5813.1213.170.381.140.8525.4316.6316.401.38
108.588.342.8011.0410.901.271.190.9718.516.9016.850.35
Crushed sample18.438.400.359.039.000.331.041.003.8517.0216.900.71
28.098.272.2212.3712.632.101.181.234.2416.7416.630.65
37.887.961.0111.7511.502.131.181.134.2416.6616.560.60
47.878.042.0312.6212.791.431.081.2414.816.7716.591.07
57.597.762.2412.9412.910.231.021.1513.716.5916.480.66
67.777.962.4513.2613.051.581.011.2220.716.5316.420.72
77.988.152.1313.2113.421.591.081.134.616.6416.480.90
87.898.102.6612.8112.542.031.101.200.0916.7416.600.83
97.747.902.0713.1213.170.381.141.094.3816.6316.480.90
108.588.782.3311.0411.090.451.191.211.6816.9016.810.53
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Gao, F.; Zhang, Y.; Liu, X. A Study of the Reliability and Accuracy of the Real-Time Detection of Forage Maize Quality Using a Home-Built Near-Infrared Spectrometer. Foods 2022, 11, 3490. https://doi.org/10.3390/foods11213490

AMA Style

Gao F, Zhang Y, Liu X. A Study of the Reliability and Accuracy of the Real-Time Detection of Forage Maize Quality Using a Home-Built Near-Infrared Spectrometer. Foods. 2022; 11(21):3490. https://doi.org/10.3390/foods11213490

Chicago/Turabian Style

Gao, Fei, Yuejing Zhang, and Xian Liu. 2022. "A Study of the Reliability and Accuracy of the Real-Time Detection of Forage Maize Quality Using a Home-Built Near-Infrared Spectrometer" Foods 11, no. 21: 3490. https://doi.org/10.3390/foods11213490

APA Style

Gao, F., Zhang, Y., & Liu, X. (2022). A Study of the Reliability and Accuracy of the Real-Time Detection of Forage Maize Quality Using a Home-Built Near-Infrared Spectrometer. Foods, 11(21), 3490. https://doi.org/10.3390/foods11213490

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