Skip to main content


Front. Chem., 07 May 2018
Sec. Analytical Chemistry
Volume 6 - 2018 |

Discovery of the Linear Region of Near Infrared Diffuse Reflectance Spectra Using the Kubelka-Munk Theory

Shengyun Dai Xiaoning Pan Lijuan Ma Xingguo Huang Chenzhao Du Yanjiang Qiao* Zhisheng Wu*
  • Key Laboratory of TCM-Information Engineering of State Administration of TCM, Pharmaceutical Engineering and New Drug Development of Traditional Chinese, Medicine of Ministry of Education, Beijing University of Chinese Medicine, Beijing, China

Particle size is of great importance for the quantitative model of the NIR diffuse reflectance. In this paper, the effect of sample particle size on the measurement of harpagoside in Radix Scrophulariae powder by near infrared diffuse (NIR) reflectance spectroscopy was explored. High-performance liquid chromatography (HPLC) was employed as a reference method to construct the quantitative particle size model. Several spectral preprocessing methods were compared, and particle size models obtained by different preprocessing methods for establishing the partial least-squares (PLS) models of harpagoside. Data showed that the particle size distribution of 125–150 μm for Radix Scrophulariae exhibited the best prediction ability with Rpre2 = 0.9513, RMSEP = 0.1029 mg·g−1, and RPD = 4.78. For the hybrid granularity calibration model, the particle size distribution of 90–180 μm exhibited the best prediction ability with Rpre2 = 0.8919, RMSEP = 0.1632 mg·g−1, and RPD = 3.09. Furthermore, the Kubelka-Munk theory was used to relate the absorption coefficient k (concentration-dependent) and scatter coefficient s (particle size-dependent). The scatter coefficient s was calculated based on the Kubelka-Munk theory to study the changes of s after being mathematically preprocessed. A linear relationship was observed between k/s and absorption A within a certain range and the value for k/s was >4. According to this relationship, the model was more accurately constructed with the particle size distribution of 90–180 μm when s was kept constant or in a small linear region. This region provided a good reference for the linear modeling of diffuse reflectance spectroscopy. To establish a diffuse reflectance NIR model, further accurate assessment should be obtained in advance for a precise linear model.


The implementation of process analytical technology (PAT) in the pharmaceutical industry is intended to enhance the quality of products through the measurement of critical quality and performance parameters (Roggo and Ulmschneider, 2008). Near infrared spectroscopy (NIRS) is regarded as a vital tool for the implementation of PAT, as it is increasingly used in pharmaceutical research and development due to its high analysis speed, low-cost, and non-destructive characteristics (De Beer et al., 2011). NIR spectra of chemical species (consisting of C–H, N–H, O–H, and S–H bonds; Sarraguça et al., 2011) can be used to predict their chemical and physical properties (Prieto et al., 2009).

The NIR technology includes two main parts that are transmission spectroscopy and diffuse reflectance spectroscopy. The selection of spectral form is mainly based on the state of samples (i.e., transmission spectroscopy is suitable for liquid samples such as herbal extracts and liquid preparations, while diffuse reflectance spectroscopy is generally used for solid samples such as pharmaceutical powders or granules). Diffuse reflectance spectroscopy is an analytical technique that measures the diffuse reflection of different wavelengths of light to obtain the surface information of the materials.

Various physical, chemical, and biochemical properties in Mediterranean soils were NIR predicted (Zornoza et al., 2008). Chen et al. employed an NIR model for the analysis of total polyphenol content in green tea (Chen Q. et al., 2008). Classification accuracy of about 100 % was obtained by discriminant and classification tree analyses of 82 honey samples by diffuse reflectance mid-infrared Fourier transform spectroscopy (DRIFTS) (Bertelli et al., 2007). Borin et al. utilized NIR technology for the simultaneous quantification of some common adulterants (starch, whey, or sucrose) found in milk powder samples (Borin et al., 2006). All these investigations have illustrated the trend of using NIR technology to predict physical and chemical information.

Recently, the application of NIR in studying Chinese herbal medicine (CHM) has dramatically increased such as discrimination analysis and quality control for various samples e.g., raw materials, excipients, and dosage forms. Wu et al. used the NIR and different PLS models to quantify the baicalin contents of Yinhuang oral solution based on a total error concept (Wu et al., 2013). Chen et al. employed NIR to distinguish Ganoderma lucidum samples collected from different geographical origins using principal component analysis (PCA) and discriminant analysis algorithms (Chen Y. et al., 2008).

On the other hand, it is well known that the particle size of sample affects NIR spectra. Several studies have been published on the effect of particle size on the determination of drug content in mixed powder products (Norris and Williams, 1984; Aucott and Garthwaite, 1988; Bull, 1991). Franke et al. (1998) reported the particle size determination of lactose using chemometrics-based NIR spectra. However, they did not mention any basic principle to determine particle size in the experiments. Paskatan et al. (2001) reviewed theoretical and practical particle size analysis of powder by NIR spectroscopy. But they did not show the relationship between the basic light scattering principle and the particle size of main contents.

Kubelka-Munk theory (Otsuka, 2004) is the basic quantitative theory of NIRS. The particle size of sample affects the light scattering, directly influencing model construction. It was shown that an accurate knowledge of the particles is crucial in the product development (Blanco and Peguero, 2008). Meanwhile, the differences in CHM particle size could result in different optical path lengths and multiplicative light scattering effects (Jin et al., 2012). Thus, it is important to establish an expeditious method to determine the particle size of CHM.

However, there were a few NIR studies on the simultaneous determination of particle size and active pharmaceutical ingredients of CHM. Wu Z. S. et al. demonstrated that the particle size affected NIR measurement of saikosaponin A in Bupleurum chinense DC (Wu et al., 2015). Bittner et al. employed a successful application of NIR spectroscopy in combination with multivariate data analysis (MVA) for the simultaneous identification and particle size determination of amoxicillin trihydrate particles (Bittner et al., 2011).

Scrophularia radix (Xuanshen), the root of Scrophularia ningpoensis Hemsl., was a typical CHM with a history going back over 1000 years (The State Pharmacopoeia Commission of People's Republic of China, 2015). It is originally from Zhejiang province and it is a component of the natural herbal supplement named “Zhe Ba Wei.” The major ingredients of Scrophularia radix are iridoids, and harpagoside is one of the main bioactive components with antioxidant, antimicrobial and antitumor activities (Miyazawa and Okuno, 2003; Jing et al., 2011).

In this study, Scrophularia radix was taken as an example and harpagoside was regarded as an API of Scrophularia radix. HPLC was used as a reference method to determine the harpagoside content. NIR was used to monitor the prediction potential of the models of single particle size and mix particle size simultaneously. To our best knowledge, this paper is the first to study on particle size and harpagoside determination in Scrophularia radix with NIR diffuse reflectance spectroscopy. The differences between single particle size model and mix particle size model from the perspective of the Kubelka-Munk theory were explained.

Materials and Method


Ten batches of S. ningpoensis Hemsl. radix were gifted from Daozhen (Guizhou, China), three representative samples were taken from each batch. All samples were identified by Prof. Chunsheng Liu (Beijing University of Chinese Medicine, China). Harpagoside reference standard (lot: 111730-201307) was purchased from the National Institutes for Food and Drug Control (Beijing, China). Acetonitrile (Fisher Scientific, Pittsburgh, PA) was of HPLC-grade. Acetic acid (Beijing Chemical Works, Beijing, China) was of analytical grade. Deionised water was purchased from Hangzhou Wahaha Co., Ltd (Zhejiang, China).

Preparation of Samples

Scrophularia radix samples were crushed into pieces by a disintegrator after brushing off soil dust from the surface. Thirty samples of Scrophularia radix were then pulverized with a blender and screened through a 10-mesh sieve. Finally, the powders were divided into four parts. One part was used for HPLC determination of the harpagoside content. The remaining parts were then smashed and screened through 24-, 50-, 65-, 80-, 100-, 120-, and 150-mesh sieves.

An amount of each sieved sample of Scrophularia radix powder (1 g) was accurately weighed and placed in a 100-mL Erlenmeyer flask. The sample was extracted with 50 mL of 50% ethanol under ultrasonic vibration (40 kHZ, 220 V) for 45 min. After cooling to room temperature, the solution was filtered through a 0.45-μm membrane filter for HPLC analysis.

NIR Equipment and Measurement

The NIR spectra were recorded by a XDS Rapid Content Analyser and VISION software (Metrohm NIR Systems, Florida, USA). The wavelength range for the spectra was 780–2,500 nm. Each spectrum was an average of 64 scans with air as the background, and the wavelength increment was of 0.5 nm. Unless stated otherwise, each sample was measured in triplicate and its mean value was used in the subsequent analysis.

HPLC Method

A certain amount of harpagoside standard was accurately weighed with an XS205DU electronic balance (Mettler Toledo, Greifensee, Switzerland) and then dissolved in 100 mL of methanol to obtain the concentration of 0.02432 mg·mL−1.

HPLC analysis of Scrophularia radix (according to Chinese Pharmacopoeia, 2010 ed) was carried out using a Waters 2695 HPLC system, Waters 2996 DAD detector and auto-sampler (Waters Technologies, Palo Alto, CA). Ten microliters aliquots of the sample solutions were chromatographically analyzed in gradient elution mode on an octadecylsilyl column [250 × 4.6 mm, 5 μm (Dikma, China)] with the mobile phase consisting of acetonitrile and 0.4% acetic acid (v/v) at a flow rate of 1.0 mL·min−1 (Table 1). The column temperature was kept at 30°C and the detection wavelength set at 280 nm. This chromatographic method exhibited good linearity (Y = 3 × 106X−104747, R2 = 0.9998) over the concentration range 0.04864–0.02432 mg·mL−1.


Table 1. HPLC gradient elution of Scrophularia radix extract.


Data analysis was performed by the Unscrambler version 9.6 software package (CAMO Software AS, Oslo, Norway) and home-made routines programmed in MATLAB code (MATLAB v7.0, Math Works, Natick, MA). Following the Kennard-Stone algorithm, 210 samples were divided into 140 calibration samples and 70 validation samples. The root mean square error of calibration (RMSEC), root mean square error of cross-validation (RMSECV), root mean square error of prediction (RMSEP) and corresponding R2 were used to evaluate the PLS model.

In order to establish a robust harpagoside model, a number of preprocessing methods were selected. For instance, multiplicative scatter correction (MSC) and standard normal variate (SNV) were used to eliminate redundant effects of particle size. Derivative methods including first derivative (1D) and second derivative (2D) were obtained to reduce baseline variations observed in original diffuse reflectance spectra and to enhance spectral features. Meanwhile, a nine-point Savitzky-Golay smoothing filter (SG) was employed to depress the background noise amplified by the derivative. For the particle size model, MSC, SNV, and second derivative were not appropriate for an effect to be modeled, so 1D + SG, normalization and baseline subtraction were used. Leave-one-out cross-validation was used to validate the validity of methods. The lowest predicted residual sum of squares (PRESS) value was used to determine the optimum latent variables.

Quantitative Models of NIR Diffuse Reflectance Using the Kubelka-Munk Theory

Kubelka-Munk theory is the theoretical basis for the establishment of quantitative models of NIR diffuse reflectance and its function is as follows (Otsuka, 2004):


According to the Kubelka-Munk function, reflectance is inversely to proportional to the light-scattering coefficient (s), and the s value is inversely proportional to particle size.

The absorbance of NIR diffuse reflectance is expressed by the Kubelka-Munk equation:


Results and Discussion

Spectral Characteristics of NIR Diffuse Reflectance Spectra of Different Particle Size Samples

The representative raw spectra of Scrophularia radix with different particle sizes are shown in Figure 1 i.e., the spectral profiles were similar in shape. However, the main influences of particle size variation on diffuse reflectance spectra was the baseline offset. The well-known phenomenon that larger particles showed a stronger absorption, illustrates that the particle size is vital to the response. Some weak absorption peaks were demonstrated in the second overtone region (SCOT, 1,000–1,400 nm) of the fundamental C-H stretching bands, while much fluctuations in the region of first combination-overtone (FCOT, 1,400–2,040 cm−1) and combination region (CR, 2,040-2,500 nm) were observed. Those absorption peaks might be caused by the diffuse reflectance on different particle sizes.


Figure 1. NIR diffuse reflectance spectra of Scrophularia radix.

HPLC Determination of Harpagoside Content in Scrophularia Radix

The HPLC chromatograms of the representative sample and standard are shown in Figure 2. The retention time of harpagoside in a sample extract was the same as that for the standard solution. Figure 3 shows the harpagoside concentration of 30 samples. There is a significant difference in harpagoside concentration of samples of different particle sizes. The biggest difference of the particle sizes was located in the range of 180–250 μm, but the overall concentration design was suitable for the modeling.


Figure 2. Representative HPLC chromatograms of Scrophularia radix sample and harpagoside standard.


Figure 3. Harpagoside concentration of 30 samples of different particle sizes.

PLS Models for NIR Diffuse Reflectance Data Using Scrophularia Radix of Each Single Particle Size

Based on different preprocessing methods, the PLS model for each particle size was constructed. Figure 4 showed the relationship between the latent variables and PRESS for different preprocessing methods. In general, the lowest PRESS value means the best latent variables (Pan et al., 2015). The model was validated for prediction by internal sample set. Moreover, the model performance values for each particle size using different preprocessing methods are illustrated in Table 2. Data showed that the raw spectra were the best to construct the particle size model of 355–850 μm and <90 μm. While the best preprocessing method for the particle size model of 250–355 μm, 180–250 μm, 150–180 μm, 125–150 μm, and 90-150 μm was EMSC, SG9, SG9, SNV, and MSC, respectively.


Figure 4. The PRESS values of different preprocessing methods for single particle size model.


Table 2. PLS model using preprocessing methods for different single particle sizes.

In addition, the model evaluation parameters, i.e., RMSEC, RMSECV, RMSEP, and RPD, for the particle size of 355–850 μm was 0.0576, 0.1642, 0.2094, and 2.02, respectively. The parameter values of other particle sizes are summarized in Table 2. The relation map between predicted value and reference value is shown in Figure 5, indicating that the best prediction result was for the particle size of 125–150 μm. Therefore, it could be known that the NIR model was influenced by different particle sizes and its quantitative characteristics was explored according to different particle sizes.


Figure 5. The relation map of the reference value and predicted value using each different particle size.

PLS Models for NIR Diffuse Reflectance Data Using Scrophularia Radix of Mix Particle Size

The comparison of model performance for different types of mix particle size (i.e., seven, six, five, four, and three types of particle size) manifests that the mix particle size model was best constructed for 3-type mix particle size (Table 3). Preprocessing methods were also various, such as MSC, SNV, EMSC and SG9. It can be seen from Figure 6, the optimum preprocessing method for the mixed particle size of 180–850 μm, 150–355 μm, 125–250 μm, 90–180 μm, and 0–150 μm was SG9, untreated original spectra, EMSC, EMSC, and MSC, respectively, as this model has the lowest PRESS value.


Table 3. Preprocessing methods for different mix particle size models (3 particle size ranges).


Figure 6. The PRESS values of different preprocessing methods for mix particle size model.

The best prediction from the mix particle size model was for 90–180 μm with RPD value >3 (Table 3). The RPD values of other mix particle size models were also about 2, meaning that the model performance of the mix particle size models was similar. This result further revealed that particle size was vital to quantitative model performance of diffuse reflectance spectra using NIR sensor. In order to make the relationship clearer, a detailed comparison of the model of the single particle size and mixed particle size was summarized.

Comparison of the Model Performance for Single Particle Size and Mix Particle Size

It can be concluded from the comparison between the single particle size and mixed particle size models that the RPD value of the former was better than the latter. Although the prediction results were good in the prediction performance in a certain particle size range by using a single particle size model, the prediction results of single particle size model were not stable. Most of the applications of NIR diffuse reflectance spectra were for a relatively broad range of particle sizes. As a result, a mix particle size calibration model was used for prediction in subsequent studies.

Moreover, the mix particle size correction model was also used to predict the validation set for each particle size for examining which particle size samples could be more accurately predicted as well as achieving the guideline for subsequent sample preparation. The model for particle size of 90–180 μm was selected to predict the particle size of 150–180 μm, 125–150 μm, and 90–125 μm and the best preprocessing method is MSC (Table 4) and RPD values of the three prediction models are 3.81, 5.78, and 2.81 (Table 5).


Table 4. The prediction model for the single particle size by using the mix particle size model.


Table 5. Predicted results of different samples of single Scrophulariaceae Radix particle size model and calibration particle size model.

On the other hand, the RPD values of the models of single particle size were 3.40, 4.78, and 2.52. Compared with the single particle size model, the RPD value of the mix particle size model was better illustrating that the prediction of the mix particle size correction model was more accurate (Table 5). The relation map between the reference and validation sets was shown in Figure 7. The correlation between reference and prediction values was good, which further demonstrated that the mix particle size model was better than the single particle size model. Why particle size was of great importance to the quantitative model of the NIR diffuse reflectance? It was performed by the Kubelka-Munk theory, which is a critical theory in the NIR diffuse reflectance.


Figure 7. The relation map of the calibration particle size models.

Discovery of the Linear Region of NIR Diffuse Reflectance Spectra Using the Kubelka-Munk Theory

In practice, NIR diffuse reflectance is usually used for solid particle determination and its quantitative evidence is based on the Kubelka-Munk theory (Figure 8).


Figure 8. The NIR diffuse reflectance and Kubelka-Munk theory.

It can be learnt from the equation that the absorbance had relationship with the k/s value. A linear relationship was discovered between k/s value and A within a certain range. As illustrated in Figure 9, the value for k/s was >4 obviously indicating that a linear region existed. This results also explained and guided the modeling performance of NIR diffuse reflectance. It was found that such a linear region provides a reference for the linear modeling of diffuse reflectance spectra. It is important to note that the linear region is beneficial for establishing a NIR diffuse reflectance model. According to our data, when the scatter coefficient s does not change, the absorption coefficient k is proportional to the sample concentration. In this study, the quantitative models for single particle size and mix particle size were both constructed to minimize the limitation that the particle size of samples was only available in a certain range. The model of single particle size was better than the mix particle size owing to a small change in the scattering coefficient s.


Figure 9. The relationship between the absorbance (A) and k/s value.


Particle size is of great importance to the quantitative model of the NIR diffuse reflectance. In this study, the single particle size and mix particle size models of Radix Scrophulariae were constructed using PLS methods. For the single particle size model, it was obvious that the best prediction model was for the particle size distribution of 125–150 μm. This particle size distribution illustrated that small particle size was beneficial to construct the quantitative model of harpagoside in Radix Scrophulariae.

For the mix particle size model, a better prediction was obtained for the particle size distribution of 90–180 μm indicating that the mix particle size model could explain more variation in the sample, and the accuracy and robustness of the mix particle size model would be improved. Meanwhile, the quantitative evidence of NIR diffuse reflectance of different particle sizes was based on the Kubelka-Munk theory. A linear relationship was discovered between k/s value and A within a certain range. Data showed that a narrow range of the scatter coefficients s resulted in a better model. Besides, the value for k/s was >4 clearly indicating that a linear region exited. This linear region helped explain and guide the modeling performance of NIR diffuse reflectance data. Finding such a linear region provided a methodological reference for the linear modeling of NIR diffuse reflectance spectra. Thus, further accurate assessment should be obtained in advance for a precise linear model.

Our study also showed that the quantitative analysis of CHM samples was more accurate when the scattering coefficient s remains unchanged or differs insignificantly at theoretical level.

Author Contributions

ZW and YQ: conceived the research; XP: performed the experiment; SD: wrote the manuscript; CD, LM, and XH: analyzed the data. All the authors prepared the manuscript and discussed the results.

Conflict of Interest Statement

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.


This work was supported by the National Natural Science Foundation of China (81773914), Beijing Nova Program of China (xx2016050), and Science Fund for Distinguished Young Scholars in BUCM (2015-JYB-XYQ-003). The authors thank the Key Laboratory of TCM Information Engineering of State Administration of Traditional Chinese Medicine, (Beijing, China) for the assistance in data processing, Modernization of Traditional Chinese Medicine of Daozhen county of China.


NIR, Near Infrared Diffuse; PAT, Process Analytical Technology; DRIFTS, Diffuse Reflectance Mid-infrared Fourier Transform Spectroscopy; CHM, Chinese Herbal Medicine; NIRS, NIR Spectroscopy; MVA, Multivariate Data Analysis; HPLC, High Performance Liquid Chromatography; RMSEC, Root Mean Square Error of Calibration; RMSECV; Root Mean Square Error of Cross-Validation; RMSEP, Root Mean Square Error of Prediction; MSC, Multiplicative Scatter Correction; SNV, Standard Normal Variate; 1D, First Derivative; 2D, Second Derivative; SG, Savitzky–Golay; PRESS, Predicted Residual Sum of Squares; PLS, Partial Least Squares; SCOT, Second Overtones Region; FCOT, First Combination-Overtone; RPD, Residual Predictive Deviation; API, Active Pharmaceutical Ingredient; EMSC, Extended Multiplicative Scatter Correction.


Aucott, L. S., and Garthwaite, P. H. (1988). Transformations to reduce the effect of particle size in near-infrared spectra. Analyst 113, 1849–1854. doi: 10.1039/an9881301849

CrossRef Full Text | Google Scholar

Bertelli, D., Plessi, M., Sabatini, A. G., Lollia, M., and Grillenzonib, F. (2007). Classification of Italian honeys by mid-infrared diffuse reflectance spectroscopy (DRIFTS). Food Chem. 101, 1565–1570. doi: 10.1016/j.foodchem.2006.04.010

CrossRef Full Text | Google Scholar

Bittner, L. K., Heigl, N., Petter, C. H., Noisternig, M. F., Griesser, U. J., Bonn, G. K., et al. (2011). Near-infrared reflection spectroscopy (NIRS) as a successful tool for simultaneous identification and particle size determination of amoxicillin trihydrate. J. Pharmaceut. Biomed. 54, 1059–1064. doi: 10.1016/j.jpba.2010.12.019

PubMed Abstract | CrossRef Full Text | Google Scholar

Blanco, M., and Peguero, A. (2008). An expeditious method for determining particle size distribution by near infrared spectroscopy: comparison of PLS2 and ANN models. Talanta 77, 647–651. doi: 10.1016/j.talanta.2008.07.015

CrossRef Full Text | Google Scholar

Borin, A., Ferrao, M. F., Mello, C., Maretto, D. A., and Poppi, R. J. (2006). Least-squares support vector machines and near infrared spectroscopy for quantification of common adulterants in powdered milk. Anal. Chim. Acta 579, 25–32. doi: 10.1016/j.aca.2006.07.008

PubMed Abstract | CrossRef Full Text | Google Scholar

Bull, C. R. (1991). Compensation for particle size effects in near infrared reflectance. Analyst 116, 781–786. doi: 10.1039/an9911600781

CrossRef Full Text | Google Scholar

Chen, Q., Zhao, J., Liu, M., Cai, J., and Liu, J. (2008). Determination of total polyphenols content in green tea using FT-NIR spectroscopy and different PLS algorithms. J. Pharm. Biomed. Anal. 46, 568–573. doi: 10.1016/j.jpba.2007.10.031

PubMed Abstract | CrossRef Full Text | Google Scholar

Chen, Y., Xie, M. Y., Yan, Y., Zhu, S. B., Nie, S. P., Li, C., et al. (2008). Discrimination of Ganoderma lucidum according to geographical origin with near infrared diffuse reflectance spectroscopy and pattern recognition techniques. Anal. Chim. Acta 618, 121–130. doi: 10.1016/j.aca.2008.04.055

PubMed Abstract | CrossRef Full Text | Google Scholar

De Beer, T., Burggraeve, A., Fonteyne, M., Saerens, L., Remon, J. P., and Vervaet, C. (2011). Near infrared and Raman spectroscopy for the in-process monitoring ofpharmaceutical production processes. Int. J. Pharm. 417, 32–47. doi: 10.1016/j.ijpharm.2010.12.012

PubMed Abstract | CrossRef Full Text | Google Scholar

Franke, P., Gill, I., Luscombe, C. N., Rudd, D. R., Waterhouse, J., and Jayasooriya, U. A. (1998). Near-infrared mass median particle size determination of lactose monohydrate. Evaluating several chemometric approaches. Analyst 123, 2043–2046. doi: 10.1039/a802532k

CrossRef Full Text | Google Scholar

Jin, J. W., Chen, Z. P., Li, L. M., Steponavicius, R., Thennadil, S. N., Yang, J., et al. (2012). Quantitative spectroscopic analysis of heterogeneous mixtures: the correction of multiplicative effects caused by variations in physical properties of samples. Anal. Chem. 84, 320–326. doi: 10.1021/ac202598f

PubMed Abstract | CrossRef Full Text | Google Scholar

Jing, J., Chan, C., Xu, L., Jin, D., Cao, X., Mok, D. K. W., et al. (2011). Development of an in-line HPLC fingerprint ion-trap mass spectrometric method for identification and quality control of Radix Scrophulariae. J. Pharmaceut. Biomed. 56, 830–835. doi: 10.1016/j.jpba.2011.07.032

PubMed Abstract | CrossRef Full Text | Google Scholar

Miyazawa, M., and Okuno, Y. (2003). Volatile components from the roots of Scrophularia ningpoensis Hemsl. Flavour Frag. J. 18, 398–400. doi: 10.1002/ffj.1232

CrossRef Full Text | Google Scholar

Norris, K. H., and Williams, P. C. (1984). Optimization of mathematical treatments of raw near-infrared signal in the measurement of protein in hard red spring wheat: I. Influence of particle size. Cereal Chem. 61, 158–165.

Google Scholar

Otsuka, M. (2004). Comparative particle size determination of phenacetin bulk powder by using Kubelka–Munk theory and principal component regression analysis based on near-infrared spectroscopy. Powder Technol. 141, 244–250. doi: 10.1016/j.powtec.2004.01.025

CrossRef Full Text | Google Scholar

Pan, X. L., Li, F. Y., Wu, Z. S., Zhang, Q., Lin, Z. Z., Shi, X. Y., et al. (2015). Near infrared spectroscopy model development and variable importance in projection assignment of particle size and lobetyolin content of Codonopsis radix. J. Near Infrared Spec. 23, 327–335. doi: 10.1255/jnirs.1175

CrossRef Full Text | Google Scholar

Paskatan, M. C., Steele, J. L., Spillman, C. K., and Haque, E. (2001). Near infrared reflectance spectroscopy for online particle size analysis of powders and ground materials. J. Near Infrared Spectrosc. 9, 153–164. doi: 10.1255/jnirs.303

CrossRef Full Text | Google Scholar

Chinese Pharmacopoeia Commission (2010). Chinese Pharmacopoeia Vol. 1. Beijing: China Medical Science Press, 108.

Prieto, N., Ross, D. W., Navajas, E. A., Nute, G. R., Richardson, R. I., Hyslop, J. J., et al. (2009). On-line application of visible and near infrared reflectance spectroscopy to predict chemical-physical and sensory characteristics of beef quality. Meat Sci. 83, 96–103. doi: 10.1016/j.meatsci.2009.04.005

PubMed Abstract | CrossRef Full Text | Google Scholar

Roggo, Y., and Ulmschneider, M. (2008). “Chapter 4.3: Chemical imaging and chemometrics: useful tools for process analytical technology,” in Pharmaceutical Manufacturing Handbook: Regulations and Quality, ed S. C. Gad (New York, NY: John Wiley & Sons), 411–431.

Sarraguça, M. C., Cruz, A. V., Amaral, H. R., Costa, P. C., and Lopes, J. A. (2011). Comparison of different chemometric and analytical methods for the prediction of particle size distribution in pharmaceutical powders. Anal. Bioanal. Chem. 399, 2137–2147. doi: 10.1007/s00216-010-4230-6

PubMed Abstract | CrossRef Full Text | Google Scholar

The State Pharmacopoeia Commission of People's Republic of China (2015). Pharmacopoeia of the People's Republic of China Press. China Medical Science and Technology Press. Beijing.

Wu, Z. S., Du, M., Shi, X. Y., Xu, B., and Qiao, Y. (2015). Robust PLS prediction model for Saikosaponin A in Bupleurum chinense DC. coupled with granularity-hybrid calibration set. J. Anal. Methods Chem. 2015:583841. doi: 10.1155/2015/583841

PubMed Abstract | CrossRef Full Text | Google Scholar

Wu, Z., Ma, Q., Lin, Z., Peng, Y., Ai, L., Shi, X., et al. (2013). A novel model selection strategy using total error concept. Talanta 107, 248–254. doi: 10.1016/j.talanta.2012.12.057

PubMed Abstract | CrossRef Full Text | Google Scholar

Zornoza, R., Guerrero, C., Mataix-Solera, J., Scow, K. M., Arcenegui, V., and Mataix-Beneyto, J. (2008). Near infrared spectroscopy for determination of various physical, chemical and biochemical properties in Mediterranean soils. Soil Biol. Biochem. 40, 1923–1930. doi: 10.1016/j.soilbio.2008.04.003

PubMed Abstract | CrossRef Full Text | Google Scholar

Keywords: Kubelka-Munk theory, Near infrared (NIR) diffuse reflectance spectroscopy, particle size, PLS, harpagoside, Radix Scrophulariae

Citation: Dai S, Pan X, Ma L, Huang X, Du C, Qiao Y and Wu Z (2018) Discovery of the Linear Region of Near Infrared Diffuse Reflectance Spectra Using the Kubelka-Munk Theory. Front. Chem. 6:154. doi: 10.3389/fchem.2018.00154

Received: 30 November 2017; Accepted: 19 April 2018;
Published: 07 May 2018.

Edited by:

Hoang Vu Dang, Hanoi University of Pharmacy, Vietnam

Reviewed by:

Eleonora-Mihaela Ungureanu, Politehnica University of Bucharest, Romania
Michalina Kotyczka-Moranska, Institute for Chemical Processing of Coal, Poland

Copyright © 2018 Dai, Pan, Ma, Huang, Du, Qiao and Wu. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.

*Correspondence: Yanjiang Qiao,
Zhisheng Wu,