留言板

尊敬的读者、作者、审稿人, 关于本刊的投稿、审稿、编辑和出版的任何问题, 您可以本页添加留言。我们将尽快给您答复。谢谢您的支持!

姓名
邮箱
手机号码
标题
留言内容
验证码

木薯块根灰分和水分近红外光谱预测模型的构建与优化

张逸杰 王思琦 陆小静 宋记明 王睿 曹敏 张瑞 王红刚 吴金山

张逸杰, 王思琦, 陆小静, 宋记明, 王睿, 曹敏, 张瑞, 王红刚, 吴金山. 木薯块根灰分和水分近红外光谱预测模型的构建与优化[J]. 热带生物学报, 2024, 15(3): 259-267. doi: 10.15886/j.cnki.rdswxb.20230031
引用本文: 张逸杰, 王思琦, 陆小静, 宋记明, 王睿, 曹敏, 张瑞, 王红刚, 吴金山. 木薯块根灰分和水分近红外光谱预测模型的构建与优化[J]. 热带生物学报, 2024, 15(3): 259-267. doi: 10.15886/j.cnki.rdswxb.20230031
ZHANG Yijie, WANG Siqi, LU Xiaojing, SONG Jiming, WANG Rui, CAO Min, ZHANG Rui, WANG Honggang, WU Jinshan. Construction and optimization of prediction models for cassava ash and water content by using near-infrared spectroscopy[J]. Journal of Tropical Biology, 2024, 15(3): 259-267. doi: 10.15886/j.cnki.rdswxb.20230031
Citation: ZHANG Yijie, WANG Siqi, LU Xiaojing, SONG Jiming, WANG Rui, CAO Min, ZHANG Rui, WANG Honggang, WU Jinshan. Construction and optimization of prediction models for cassava ash and water content by using near-infrared spectroscopy[J]. Journal of Tropical Biology, 2024, 15(3): 259-267. doi: 10.15886/j.cnki.rdswxb.20230031

木薯块根灰分和水分近红外光谱预测模型的构建与优化

doi: 10.15886/j.cnki.rdswxb.20230031
基金项目: 

国家现代木薯产业技术体系 CARS-11-HNCYH

详细信息
    第一作者:

    张逸杰(1998-)男,海南大学林学院2020级硕士研究生。E-mail:530259881@qq.com

    通信作者:

    吴金山(1984-)男,博士,副教授,研究方向:主要研究方向为植物生理生化。E-mail: wujsh2007@163.com

  • 叶静
  • 中图分类号: S533

Construction and optimization of prediction models for cassava ash and water content by using near-infrared spectroscopy

  • 摘要: 为快速检测木薯的灰分和水分含量,以同地块同一时期木薯种质资源为材料进行建模,采用GB 5009.3—2016和GB/T5009.4—2016法对木薯灰分和水分含量进行测定,同时使用近红外光谱分析仪对137份样品进行光谱采集,利用TQ Analyst 9.0分析软件,采用偏最小二乘法(PLS)构建木薯灰分和水分近红外定标模型。实验结果显示, 木薯灰分、水分模型相关系数(R)分别为0.94、0.93,校正均方根误差(RMSEC)分别为0.22、0.48,预测均方根误差(RMSEP)分别为0.21、1.46,交叉验证均方差(RMSECV)分别为0.40、1.54;选用未参与建模的20份木薯种质资源对该模型进行外部验证,预测值与真实值进行差异性分析(P>0.05),P值分别为0.464、0.459说明差异无显著性,表明该模型可适用于木薯灰分和水分检测。
    注释:
    1)  叶静
  • 图  2  木薯样品光谱数据方差分析图

    Fig.  2  Analysis of variance chart

    图  3  木薯样品预测值和参考值的相关图

    Fig.  3  Correlation diagram between predicted value and reference value of cassava sample

    表  1  不同木薯品种的水分和灰分含量

    Table  1  Water and ash contents of cassava Germplasm Resources

    编号No.品种Name灰分/%Ash/%水分/%Water/%编号No.品种Name灰分/%Ash/%水分/%Water/%
    1F20002.78 ± 0.044.24 ± 0.2528SC2011.98 ± 0.062.95 ± 0.33
    2东莞红尾1.94 ± 0.044.19 ± 0.3729F2022.14 ± 0.053.60 ± 0.16
    3Q101.89 ± 0.004.86 ± 0.3130F03P2.34 ± 0.024.50 ± 0.35
    4JG13012.53 ± 0.033.65 ± 0.3931COL7122.06 ± 0.044.71 ± 1.39
    5SC92.14 ± 0.064.81 ± 0.2532SC21.71 ± 0.014.14 ± 0.86
    6MF-元引9831.19 ± 0.023.34 ± 0.8133CH201.39 ± 0.014.87 ± 1.27
    72B1.75 ± 0.023.58 ± 2.1334SC62.04 ± 0.065.13 ± 1.44
    8老板娘3号1.68 ± 0.015.21 ± 0.0735COL7151.93 ± 0.054.95 ± 0.65
    9GR32.24 ± 0.015.05 ± 0.3336T122.00 ± 0.023.39 ± 0.57
    10ME1911.18 ± 0.035.43 ± 0.283782292.07 ± 0.024.29 ± 0.78
    11C31.29 ± 0.033.14 ± 2.5238ITBB011.97 ± 0.033.90 ± 0.30
    12泰国1.15 ± 0.004.73 ± 0.0839KU582.14 ± 0.054.02 ± 0.24
    13云南思茅1.72 ± 0.003.56 ± 0.0440ME1902.26 ± 0.053.31 ± 0.05
    14老板娘2号1.03 ± 0.033.22 ± 0.7841CM39921.97 ± 0.054.18 ± 0.66
    15GR9111.95 ± 0.044.46 ± 0.3042SC2051.72 ± 0.004.31 ± 0.10
    16KU501.83 ± 0.033.90 ± 0.4743C2221.74 ± 0.005.92 ± 3.29
    17SC124多倍体1.84 ± 0.063.27 ± 0.9844SC101.64 ± 0.023.47 ± 1.24
    187751.72 ± 0.002.36 ± 1.1845F2972.06 ± 0.024.74 ± 0.17
    19Dec-461.16 ± 0.034.45 ± 0.7146水果木薯1.72 ± 0.013.87 ± 0.61
    20SC162.03 ± 0.044.47 ± 0.5247SC1241.64 ± 0.013.20 ± 0.11
    21哥伦比亚18R2.10 ± 0.052.43 ± 0.0548SC122.16 ± 0.044.42 ± 2.23
    22SC41.41 ± 0.051.89 ± 1.2149南美9号2.25 ± 0.041.57 ± 0.42
    23利民1.73 ± 0.015.46 ± 2.8650花叶木薯1.93 ± 0.044.21 ± 0.28
    24罗勇801.23 ± 0.003.30 ± 0.4551罗勇901.72 ± 0.017.07 ± 0.27
    25COL7131.05 ± 0.023.54 ± 0.4052宝岛9-51.96 ± 0.015.78 ± 0.01
    26SC51.12 ± 0.034.47 ± 0.2653E13401.52 ± 0.016.08 ± 0.74
    27GR61.91 ± 0.044.16 ± 1.0054ME1891.39 ± 0.014.78 ± 0.55
    续表1 Tab.1 Continued
    下载: 导出CSV
    编号No.品种Name灰分/%Ash/%水分/%Water/%编号No.品种Name灰分/%Ash/%水分/%Water/%
    55NZ1881.84 ± 0.065.76 ± 0.43921963.67 ± 0.334.04 ± 0.11
    56琼中1号1.63 ± 0.044.71 ± 0.21935714.24 ± 0.193.86 ± 0.16
    57泰国种1.57 ± 0.025.82 ± 0.00944173.05 ± 0.078.45 ± 0.21
    58新选0481.65 ± 0.014.22 ± 1.9495ZM82294.12 ± 0.164.38 ± 0.21
    59H6601.54 ± 0.055.3 ± 0.39965213.88 ± 0.004.17 ± 0.08
    60SC8辐射2.41 ± 0.023.55 ± 2.6497SC80133.06 ± 0.083.33 ± 0.15
    61NZ1991.23 ± 0.034.01 ± 1.69982743.27 ± 0.074.52 ± 0.07
    62ME1871.89 ± 0.026.27 ± 0.29994284.20 ± 0.043.51 ± 0.05
    632761.78 ± 0.033.86 ± 1.0010060683.00 ± 0.056.42 ± 0.22
    6418R2.09 ± 0.034.85 ± 0.11101MB3.69 ± 0.143.61 ± 0.14
    65COC1732.11 ± 0.063.49 ± 2.23102T3B12.35 ± 0.066.1 ± 0.04
    66C91.39 ± 0.015.22 ± 0.51103T2B1.46 ± 0.036.24 ± 0.00
    67M-SC1241.77 ± 0.065.64 ± 0.01104T1A12.63 ± 0.045.20 ± 0.21
    68MF5561.25 ± 0.013.1 ± 2.34105T4C12.42 ± 0.035.47 ± 0.14
    69老板娘1号1.16 ± 0.036.39 ± 1.2106T2A11.43 ± 0.036.01 ± 0.05
    70MF811.57 ± 0.044.91 ± 0.04107T3A11.98 ± 0.055.08 ± 0.07
    71F5391.65 ± 0.014.95 ± 0.46108T4B11.32 ± 0.025.59 ± 0.07
    72M-GR9111.65 ± 0.043.35 ± 1.16109T1C11.92 ± 0.055.01 ± 0.17
    73GR8911.97 ± 0.034.24 ± 0.13110T4A11.26 ± 0.035.32 ± 0.35
    74MF6481.41 ± 0.024.19 ± 0.75111T3C13.36 ± 0.045.35 ± 0.14
    75MF6711.29 ± 0.025.10 ± 0.30112T2C12.66 ± 0.054.76 ± 0.28
    76SC72.92 ± 0.044.98 ± 0.50113T1B1.96 ± 0.045.18 ± 0.04
    77M 糯米1.49 ± 0.044.55 ± 0.521141 崖州扁豆左62.36 ± 0.035.46 ± 0.21
    78C5761.84 ± 0.034.34 ± 1.471151 撒播黑豆右33.08 ± 0.024.33 ± 0.07
    79C10301.56 ± 0.057.70 ± 0.901161 柱花草左 72.67 ± 0.013.72 ± 0.14
    80MF5321.99 ± 0.034.29 ± 0.931171 崖州扁豆左52.73 ± 0.094.22 ± 0.28
    81SC81.64 ± 0.034.50 ± 0.631181 花生左62.29 ± 0.064.53 ± 0.32
    82MF3591.19 ± 0.016.13 ± 1.751191 花生右13.26 ± 0.014.36 ± 0.24
    83M-FC101.67 ± 0.013.50 ± 0.691201 对照2右31.99 ± 0.052.41 ± 0.08
    84MSC51.07 ± 0.023.94 ± 0.931211 崖州扁豆右41.99 ± 0.023.61 ± 0.08
    85E4071.77 ± 0.003.61 ± 0.961221 柱花草右31.85 ± 0.033.25 ± 0.13
    86MF1811.12 ± 0.024.02 ± 1.171231 崖州扁豆右61.20 ± 0.023.16 ± 0.11
    87M-SC41.41 ± 0.012.47 ± 0.011241 油菜花3左12.72 ± 0.053.37 ± 0.23
    88OOO1.25 ± 0.013.96 ± 0.421251左撒崖州扁豆1.78 ± 0.034.24 ± 0.09
    89SC111.79 ± 0.005.68 ± 0.021261黑豆撒左42.59 ± 0.083.59 ± 0.35
    90F2742.26 ± 0.024.10 ± 1.681271花生左42.48 ± 0.013.17 ± 0.20
    91504.26 ± 0.074.07 ± 0.071281对照左12.04 ± 0.063.61 ± 0.23
    续表1 Tab.1 Continued
    下载: 导出CSV
    编号No.品种Name灰分/%Ash/%水分/%Water/%编号No.品种Name灰分/%Ash/%水分/%Water/%
    1291崖州扁豆撒右22.61 ± 0.072.56 ± 0.401341 油菜间作撒右12.26 ± 0.066.35 ± 0.13
    1301柱花草2右22.40 ± 0.026.44 ± 0.281351柱花草1左22.20 ± 0.027.63 ± 0.10
    1311柱花草右51.90 ± 0.036.39 ± 0.081361右8对照1.06 ± 0.027.14 ± 0.09
    1321撒播柱花草右42.29 ± 0.016.49 ± 0.071371对照左81.07 ± 0.018.15 ± 0.11
    1331花生右71.17 ± 0.047.56 ± 0.57
    注:崖州扁豆、柱花草、花生、油菜花分别代表与木薯SC9品种间作的不同植物种类。
    Note:Macrotyloma uniflorum ‘Yazhou’,stylosanthes, peanut and rape flower represent different plant species intercropping with cassava SC9 variety respectively
    下载: 导出CSV

    表  2  校正集与验证集样本信息

    Table  2  Sample information of correction set and verification set

    成分Component校正集Correction set验证集Validation set
    样本量sample size范围Range样本量sample size范围Range
    灰分/%Ash/%1171.030~4.260201.150~4.200
    水分/%Water/%1171.570~8.450202.470~7.140
    下载: 导出CSV

    表  3  不同预处理以及建模区间的考察结果

    Table  3  Investigation results of different pretreatment and modeling intervals

    指标Index预处理Pretreatment区间Section相关系数R校正标准差RMSEC预测标准差RMSEP交叉标准差RMSECV
    灰分AshConstant(原始)a0.742 80.457 00.636 00.514 0
    b0.759 80.444 00.591 00.498 0
    MSCa0.769 40.436 00.403 00.509 0
    b0.774 40.432 00.675 00.575 0
    SNVa0.770 80.435 00.404 00.518 0
    b0.729 40.467 00.633 00.579 0
    MSC+1S+SGa0.826 390.384 00.636 00.515 0
    b0.904 00.292 00.521 00.535 0
    SNV+1S+SGa0.828 80.382 00.628 00.513 0
    b0.906 20.289 00.554 00.495 0
    MSC+2S+SGa0.873 50.332 00.820 00.590 0
    b0.944 90.223 00.210 00.398 0
    SNV+2S+SGa0.876 00.329 00.810 00.623 0
    b0.925 70.270 000.247 00.495 0
    水分WaterConstant(原始)a0.788 80.784 01.280 01.440 0
    b0.693 90.919 01.300 01.440 0
    MSCa0.778 00.801 01.340 01.280 0
    b0.685 50.929 01.280 01.420 0
    SNVa0.784 10.792 01.380 01.440 0
    b0.696 20.916 01.400 01.470 0
    MSC+1S+SGa0.841 40.689 01.450 01.500 0
    b0.820 70.729 01.220 01.530 0
    SNV+1S+SGa0.844 80.683 01.450 01.480 0
    b0.884 70.595 01.370 01.640 0
    MSC+2S+SGa0.863 50.643 01.710 01.620 0
    b0.926 40.480 01.460 01.540 0
    SNV+2S+SGa0.862 80.645 01.710 01.620 0
    b0.892 60.575 01.340 01.590 0
    注:区间a:全光谱波段;区间b:光谱波段特征区间(5 800~4 300 cm-1和4 000~3 600 cm-1)。
    Note: Interval a: full spectrum band; Interval b: spectral band characteristic interval (5 800~4 300 cm-1 and 4 000~3 600 cm-1
    下载: 导出CSV
  • [1] SADRAS V O, CALDERINI D F. Preface [M]//Crop Physiology Case Histories for Major Crops. Amsterdam: Elsevier, 2021: xvii-xix.
    [2] SUMARDIONO S, JOS B, PUDJIHASTUTI I, et al. Effect of chemical modification, drying method, and drying temperature on baking expansion and the physicochemical properties of cassava starch [J]. Journal of Food Processing and Preservation, 2022, 46(1): e16111.
    [3] PITAK L, SIRISOMBOON P, SAENGPRACHATANARUG K, et al. Rapid elemental composition measurement of commercial pellets using line-scan hyperspectral imaging analysis [J]. Energy, 2021, 220: 119698.
    [4] CHEN H, TAN C, LIN Z, et al. Quantifying several adulterants of notoginseng powder by near-infrared spectroscopy and multivariate calibration [J]. Spectrochimica Acta Part A: Molecular and Biomolecular Spectroscopy, 2019, 211: 280-286.
    [5] JIANG D, QI G, HU G, et al. A residual neural network based method for the classification of tobacco cultivation regions using near-infrared spectroscopy sensors [J]. Infrared Physics & Technology, 2020, 111: 103494.
    [6] QIU G, LÜ E, WANG N, et al. Cultivar classification of single sweet corn seed using Fourier transform near-infrared spectroscopy combined with discriminant analysis [J]. Applied Sciences, 2019, 9(8): 1530.
    [7] BASATI Z, JAMSHIDI B, RASEKH M, et al. Detection of sunn pest-damaged wheat samples using visible/near-infrared spectroscopy based on pattern recognition [J]. Spectrochimica Acta Part A: Molecular and Biomolecular Spectroscopy, 2018, 203: 308-314.
    [8] CHEN H, TAN C, LIN Z, et al. Quantifying several adulterants of notoginseng powder by near-infrared spectroscopy and multivariate calibration [J]. Spectrochimica Acta Part A: Molecular and Biomolecular Spectroscopy, 2019, 211: 280-286.
    [9] TALLADA J G, PALACIOS-ROJAS N, ARMSTRONG P R. Prediction of maize seed attributes using a rapid single kernel near infrared instrument [J]. Journal of Cereal Science, 2009, 50(3): 381-387.
    [10] SUDARNO, SILALAHI D D, RISMAN T, et al. Rapid determination of oil content in dried-ground oil palm mesocarp and kernel using near infrared spectroscopy [J]. Journal of Near Infrared Spectroscopy, 2017, 25(5): 338-347.
    [11] SAHA U, JACKSON D. Analysis of moisture, oil, and fatty acid composition of olives by near-infrared spectroscopy: development and validation calibration models [J]. Journal of the Science of Food and Agriculture, 2018, 98(5): 1821-1831.
    [12] ALTIERI G, MATERA A, GENOVESE F, et al. Models for the rapid assessment of water and oil content in olive pomace by near-infrared spectrometry [J]. Journal of the Science of Food and Agriculture, 2020, 100(7): 3236-3245.
    [13] ARSLAN M, ZOU X, HU X, et al. Near infrared spectroscopy coupled with chemometric algorithms for predicting chemical components in black goji berries (Lycium ruthenicum Murr.) [J]. Journal of Near Infrared Spectroscopy, 2018, 26(5): 275-286.
    [14] JIN X, SHI C, YU C Y, et al. Determination of leaf water content by visible and near-infrared spectrometry and multivariate calibration in Miscanthus [J]. Frontiers in Plant Science, 2017, 8: 721.
    [15] HEO S, CHOI J Y, KIM J, et al. Prediction of moisture content in steamed and dried purple sweet potato using hyperspectral imaging analysis [J]. Food Science and Biotechnology, 2021, 30(6): 783-791.
    [16] RABANERA J D, GUZMAN J D, YAPTENCO K F. Rapid and Non-destructive measurement of moisture content of peanut (Arachis hypogaea L.) kernel using a near-infrared hyperspectral imaging technique [J]. Journal of Food Measurement and Characterization, 2021, 15(4): 3069-3078.
    [17] RAHMAN A, KANDPAL L, LOHUMI S, et al. Nondestructive estimation of moisture content, pH and soluble solid contents in intact tomatoes using hyperspectral imaging [J]. Applied Sciences, 2017, 7(1): 109.
    [18] ZOU X, ZHAO J, POVEY M J W, et al. Variables selection methods in near-infrared spectroscopy [J]. Analytica Chimica Acta, 2010, 667(1/2): 14-32.
    [19] OZAKI Y, CHRISTY A A, MCCLURE W F. Near-infrared spectroscopy in food science and technology[M]. NewYork: John Wiley & Sons, 2006: 23-25.
    [20] FAZELI BURESTAN N, AFKARI SAYYAH A H, SAFI M. Prediction of amylose content, protein content, breakdown, and setback viscosity of Kadus rice and its flour by near-infrared spectroscopy (NIRS) analysis [J]. Journal of Food Processing and Preservation, 2021, 45(1): e15069.
    [21] OZAKI Y, HUCK C, TSUCHIKAWA S, et al.. Near-Infrared Spectroscopy: Theory, Spectral Analysis, Instrumentation, and Applications [M]. Singapore: Springer Singapore, 2021: 57-63.
    [22] JIAO Y, LI Z, CHEN X, et al. Preprocessing methods for near-infrared spectrum calibration [J]. Journal of Chemometrics, 2020, 34(11): e3306.
    [23] LEBOT V, CHAMPAGNE A, MALAPA R, et al. NIR determination of major constituents in tropical root and Tuber crop flours [J]. Journal of Agricultural and Food Chemistry, 2009, 57(22): 10539-10547.
    [24] LEBOT V, MALAPA R, JUNG M. Use of NIRS for the rapid prediction of total N, minerals, sugars and starch in tropical root and tuber crops [J]. New Zealand Journal of Crop and Horticultural Science, 2013, 41(3): 144-153.
    [25] SÁNCHEZ T, CEBALLOS H, DUFOUR D, et al. Prediction of carotenoids, cyanide and dry matter contents in fresh cassava root using NIRS and Hunter color techniques [J]. Food Chemistry, 2014, 151: 444-451.
    [26] 李荣. 近红外光谱技术对益智质量评价的研究 [D]. 广州: 广东药科大学, 2020.
    [27] MISHRA P, ROGER J M, MARINI F, et al. Parallel pre-processing through orthogonalization (PORTO) and its application to near-infrared spectroscopy [J]. Chemometrics and Intelligent Laboratory Systems, 2021, 212: 104190.
    [28] OLIVEIRA U F, COSTA A M, ROQUE J V, et al. Predicting oil content in ripe Macaw fruits (Acrocomia aculeata) from unripe ones by near infrared spectroscopy and PLS regression [J]. Food Chemistry, 2021, 351: 129314.
    [29] LIU Y, SUN L, DU C, et al. Near-infrared prediction of edible oil frying times based on Bayesian Ridge Regression[J]. Optics, 2020, 218(12): 164950.
    [30] 高美凤, 陶焕明. 改进团队进步算法的近红外光谱波长筛选 [J]. 光谱学与光谱分析, 2021, 41(10): 3032-3038.
    [31] 陶焕明. 光谱波长筛选方法及其在近红外检测中的应用 [D]. 无锡: 江南大学, 2021.
    [32] 洪士军, 黄雯, 张立国, 等. 基于尺度不变特征变换筛选稳定特征波长的近红外光谱模型传递方法 [J]. 分析测试学报, 2020, 39(10): 1260-1266.
    [33] 吕晓菡, 蒋锦琳, 杨静, 等. 基于特征波长建模的近红外光谱技术检测辣椒素含量 [J]. 浙江大学学报(农业与生命科学版), 2019, 45(6): 760-766.
  • [1] 赵金铃, 潘晋龙, 吴宇垚, 韩昊祥, 张怀方, 陈新, 曾长英.  不同氮源处理对4个木薯品种幼苗生长发育的影响 . 热带生物学报, 2024, 15(1): 19-26. doi: 10.15886/j.cnki.rdswxb.20230018
    [2] 冯亚亭, 张逸杰, 林南方, 陈银华, 骆凯.  木薯不同组织内生细菌多样性的比较分析 . 热带生物学报, 2024, 15(2): 141-149. doi: 10.15886/j.cnki.rdswxb.20230077
    [3] 贾素行, 朱寿松, 符仁稳, 李春霞, 陈银华.  木薯MePPD3基因的克隆及功能分析 . 热带生物学报, 2024, 15(1): 10-18. doi: 10.15886/j.cnki.rdswxb.20230035
    [4] 王昕, 韩语, 潘纪汶, 杨诺, 郑立新, 曾纪锋, 郭桂英, 郑继平.  临床分离粪肠球菌的毒力基因和耐药基因检测及其耐药性研究 . 热带生物学报, 2023, 14(2): 138-144. doi: 10.15886/j.cnki.rdswxb.2023.02.001
    [5] 徐春华, 李超萍, 时涛, 王国芬, 蔡吉苗, 李博勋, 黄贵修.  木薯细菌性萎蔫病菌2个候选抗铜基因簇的克隆和分析 . 热带生物学报, 2023, 14(5): 499-505. doi: 10.15886/j.cnki.rdswxb.20220009
    [6] 李昉泽, 詹剑.  木薯渣和甘蔗渣基生物炭对砖红壤的改良效果 . 热带生物学报, 2023, 14(5): 577-584. doi: 10.15886/j.cnki.rdswxb.20230071
    [7] 郭沈沈, 吴志祥, 刘文杰, 赵豪, 李运帅, 王誉博, 王鹏.  热带北缘地区橡胶林水分利用效率的季节差异及归因分析 . 热带生物学报, 2023, 14(4): 424-432. doi: 10.15886/j.cnki.rdswxb.2023.04.010
    [8] 戴云鹏, 潘晋龙, 吴宇垚, 韩昊祥, 张怀方, 陈新, 曾长英.  2种氮素的分根处理对木薯幼苗生长发育和氮效率的影响 . 热带生物学报, 2023, 14(6): 614-621. doi: 10.15886/j.cnki.rdswxb.20230008
    [9] 童俊儒, 安锋, 谢贵水.  基于CROPWAT模型的橡胶树需水量 . 热带生物学报, 2022, 13(1): 48-56. doi: 10.15886/j.cnki.rdswxb.2022.01.008
    [10] 周诗正, 陈琳, 颜洪, 傅鹏程.  基于深度神经网络的芯片上活体虫黄藻检测 . 热带生物学报, 2022, 13(5): 451-456. doi: 10.15886/j.cnki.rdswxb.2022.05.004
    [11] 施晨晓, 羊清雯, 陈珍莉, 刘霄燕, 王小洁.  海南省自动土壤水分站数据质量控制方法研究 . 热带生物学报, 2022, 13(4): 367-375. doi: 10.15886/j.cnki.rdswxb.2022.04.007
    [12] 熊剑, 油珅, 长孙东亭, 李晓丹, 罗素兰.  不同给药时间间隔对吗啡诱导小鼠条件性位置偏爱模型的影响 . 热带生物学报, 2021, 12(3): 363-368. doi: 10.15886/j.cnki.rdswxb.2021.03.013
    [13] 廖立国, 谭正洪, 蒋龙, 符妙, 金艳, 刘应帅, 章杰.  基于过程模型的泰国热带森林碳通量动态模拟解析 . 热带生物学报, 2021, 12(3): 279-288. doi: 10.15886/j.cnki.rdswxb.2021.03.002
    [14] 杜艳楠, 王萌, 马建强, 张宇, 梁晓宇.  植物病原真菌早期检测技术及其在橡胶树炭疽病预测预报中的应用 . 热带生物学报, 2021, 12(1): 124-131. doi: 10.15886/j.cnki.rdswxb.2021.01.018
    [15] 邓雅洁, 刘云豪, 孙琪, 邱越, 宋家明, 赖杭桂.  在干旱胁迫下有性四倍体木薯叶片的蛋白质组学研究 . 热带生物学报, 2021, 12(2): 210-218. doi: 10.15886/j.cnki.rdswxb.2021.02.010
    [16] 李祥燕, 郑妍妍, 侯志文, 张欢, 廖飞雄.  白掌叶片生长定量分析与生长模型的构建 . 热带生物学报, 2020, 11(4): 461-469. doi: 10.15886/j.cnki.rdswxb.2020.04.009
    [17] 赵克学, 许淑娟, 马瑞仙, 耿金静, 刘杨, 李琼毅, 冯若飞.  TRIM21基因SYBR GreenⅠ实时荧光定量PCR检测方法的建立 . 热带生物学报, 2020, 11(1): 111-117. doi: 10.15886/j.cnki.rdswxb.2020.01.017
    [18] 赵锐, 樊若晨, 丁凯旋, 唐枝娟, 王红刚, 黄思源, 蒋凌雁, 陈银华, 耿梦婷.  木薯MeNOSIP与MeRbohD蛋白的互作验证及其对外源激素的响应 . 热带生物学报, 2020, 11(2): 170-176. doi: 10.15886/j.cnki.rdswxb.2020.02.007
    [19] 梅成铭, 王桢昌, 陈海明, 陈文学, 云永欢, 陈卫军, 钟秋平.  冷藏期间鲜切菠萝蜜水分迁移对其品质的影响 . 热带生物学报, 2020, 11(1): 72-78. doi: 10.15886/j.cnki.rdswxb.2020.01.011
    [20] 宋娜, 张刘宁颖, 曹敏, 郭文雅, 洪雨慧, 吴金山, 陈银华, 于晓惠.  木薯ERF转录因子调控的靶基因筛选与表达分析 . 热带生物学报, 2020, 11(2): 177-189. doi: 10.15886/j.cnki.rdswxb.2020.02.008
  • 加载中
图(3) / 表 (5)
计量
  • 文章访问数:  19
  • HTML全文浏览量:  9
  • PDF下载量:  7
  • 被引次数: 0
出版历程
  • 收稿日期:  2023-03-08
  • 修回日期:  2023-05-06
  • 刊出日期:  2024-05-25

木薯块根灰分和水分近红外光谱预测模型的构建与优化

doi: 10.15886/j.cnki.rdswxb.20230031
    基金项目:

    国家现代木薯产业技术体系 CARS-11-HNCYH

    作者简介:

    张逸杰(1998-)男,海南大学林学院2020级硕士研究生。E-mail:530259881@qq.com

    通讯作者: 吴金山(1984-)男,博士,副教授,研究方向:主要研究方向为植物生理生化。E-mail: wujsh2007@163.com
  • 叶静
  • 中图分类号: S533

摘要: 为快速检测木薯的灰分和水分含量,以同地块同一时期木薯种质资源为材料进行建模,采用GB 5009.3—2016和GB/T5009.4—2016法对木薯灰分和水分含量进行测定,同时使用近红外光谱分析仪对137份样品进行光谱采集,利用TQ Analyst 9.0分析软件,采用偏最小二乘法(PLS)构建木薯灰分和水分近红外定标模型。实验结果显示, 木薯灰分、水分模型相关系数(R)分别为0.94、0.93,校正均方根误差(RMSEC)分别为0.22、0.48,预测均方根误差(RMSEP)分别为0.21、1.46,交叉验证均方差(RMSECV)分别为0.40、1.54;选用未参与建模的20份木薯种质资源对该模型进行外部验证,预测值与真实值进行差异性分析(P>0.05),P值分别为0.464、0.459说明差异无显著性,表明该模型可适用于木薯灰分和水分检测。

注释:
1)  叶静

English Abstract

张逸杰, 王思琦, 陆小静, 宋记明, 王睿, 曹敏, 张瑞, 王红刚, 吴金山. 木薯块根灰分和水分近红外光谱预测模型的构建与优化[J]. 热带生物学报, 2024, 15(3): 259-267. doi: 10.15886/j.cnki.rdswxb.20230031
引用本文: 张逸杰, 王思琦, 陆小静, 宋记明, 王睿, 曹敏, 张瑞, 王红刚, 吴金山. 木薯块根灰分和水分近红外光谱预测模型的构建与优化[J]. 热带生物学报, 2024, 15(3): 259-267. doi: 10.15886/j.cnki.rdswxb.20230031
ZHANG Yijie, WANG Siqi, LU Xiaojing, SONG Jiming, WANG Rui, CAO Min, ZHANG Rui, WANG Honggang, WU Jinshan. Construction and optimization of prediction models for cassava ash and water content by using near-infrared spectroscopy[J]. Journal of Tropical Biology, 2024, 15(3): 259-267. doi: 10.15886/j.cnki.rdswxb.20230031
Citation: ZHANG Yijie, WANG Siqi, LU Xiaojing, SONG Jiming, WANG Rui, CAO Min, ZHANG Rui, WANG Honggang, WU Jinshan. Construction and optimization of prediction models for cassava ash and water content by using near-infrared spectroscopy[J]. Journal of Tropical Biology, 2024, 15(3): 259-267. doi: 10.15886/j.cnki.rdswxb.20230031
  • 主持人:徐 冉

    木薯 (Manihot esculenta Crantz)别名木蕃薯、树薯,是世界三大薯类作物之一,在食用、饲用和工业等开发利用上均有着较大的价值[1]。木薯的灰分和水分指标与木薯品质密切相关,但在木薯品种选育过程中,利用马弗炉法和烘干法检测其所含灰分和水分的流程较为繁琐[2],效率较低,因此现阶段迫切需要一种高效率、高精度和低成本的方法来检测批量木薯灰分和水分的含量指标。随着近红外光谱技术[near-infrared (NIR) spectroscopy]与化学计量方法的快速发展,近红外光谱技术被广泛应用于农产品品质分析领域,并已成为不可或缺的快速准确的测定方法[3],是工业生产过程中质量控制的首选技术。近红外光谱技术能满足常规方法不能完成的实时监测需求,同时具有无损样品、操作简单、检测快速等特点,在成品检测和在线检测中都有优势[3-5]。研究人员在玉米[6]、小麦[7]等饲料原料中利用近红外光谱技术检测其营养成分和利用近红外光谱技术对三七[8]、玉米[9]、棕榈[10]和橄榄[11-12]果实的含油量进行定量分析。Arslan等[13]采用近红外光谱法对枸杞总黄酮和花青素的化学成分进行了定量分析。Jin等[14]利用可见光和近红外光谱(VIS-NIR)可在30 770~9 300 cm-1范围内对芒草叶片含水量进行预测。Heo等[15]构建了PLS预测模型(R2P=0.95),利用选定波长分别为961.12、1 065.50、1 083.93、1 173.23、1 233.89 nm的HIS系统预测了蒸干紫薯中的水分含量。Rabanera等[16]对波长范围为900~1 700 nm的区段利用HSI系统,测定了花生仁水分含量,并构建了R2P为0.94、RMSEP为1.95%的PLS预测模型。Rahman等[17]在波段1 000~1 550 nm范围内利用HSI技术,开发出一套无损检测番茄中水分含量和可溶性固形物含量的方法。目前国内外基于近红外光谱技术对木薯的研究多集中在表征特征分析等方面,而对其灰分和水分含量进行快速检测模型的研究较少。本研究以海南省海南大学儋州校区农科基地木薯种质资源圃为依托,拟以137份不同木薯种质块根样品为定标和验证材料,进行样品光谱收集,获取木薯有效的近红外光谱信息,并运用各种光谱数据预处理技术和统计方法进行校正、验证,确定各参数的最佳优化设置,初步建立了木薯灰分和水分含量预测模型,以期为进一步建立其他重要品质性状指标模型提供基础,为木薯育种及资源品质改良研究提供高效准确快速的分析方法。

    • 本研究所用的137份木薯种质种植于海南大学儋州校区农科基地木薯种植资源圃(海拔:380 m,109°49′E、19°51′N)。实验地块每年统一连续耕作,2021年4月种植,2022年3月收获,土质为红壤土。实验用的137份木薯块根取自同一地块,生长期同为6个月龄。取好样品后将木薯块根削皮、切段,并在60 ℃烘箱干燥24 h后,用研磨机研磨成粉状,过80目(0.18 mm)筛。粉碎后每个样品采用四分法分为2个重复,装入自封袋干燥避光处保存。

    • 布鲁克MPA近红外光谱仪(美国布鲁克公司);DHG-9070A电热鼓风干燥箱(上海一恒科学仪器有限公司);EL204电子分析天平[0.001, 梅特勒有限-托利多仪器(上海)有限公司];300Y多功能粉碎机(中国伯欧五金制品有限公司);RE52AA旋转蒸发仪(上海亚荣生化仪器厂)。

    • 采用漫反射扫描(diffuse reflectance scanning)方式扫描样品。扫描样品前将光谱仪开机预热30 min,以保证样品测试的稳定性[15]。在约25 ℃环境下,取适量处理后的木薯样品置于IN311/C旋转台上的IN311-S样品杯(直径50 mm)中,摊平,压实,以空气为空白背景扣除(每扫描1次样品粉末,均要以空气为空白进行背景扣除),用近红外光谱仪采集137份木薯样品NIRS图和CSV.格式的数据矩阵,并对137份样品进行全谱段光谱扫描,光谱范围为12 500~3 600 cm-1,测量通道为Sphere Macrosample Rotating,分辨率为16 cm-1,样品扫描64次,光谱点数为1 154。每个样品扫描3次,取3次的平均光谱。除所需样品信息外,获得的近红外光谱会因混入仪器和测试环境产生的无关信息,导致光谱基线漂移(spectral baseline drift)、散射和噪声等问题[16]。常规是通过多元散射校正(MSC)、一阶导数(1S)、二阶导数(2S)、Savitzky-Golay平滑模式(SG)、标准正态变量(SNV)等方法予以校正 [17-18]。导数处理常用来消除光谱中存在的背景干扰,Savitzky-Golay平滑消除噪声,SNV和MSC主要消除光散射[19]。在近红外光谱中,采用不同的前处理方法会得到不同的模型结果。

    • 采用GB5009.3—2016和GB/T5009.4—2016法对木薯灰分和水分含量进行测定。将木薯样品在105 ℃烘箱中烘干过夜,直至达到恒重,以获得水分含量。灰分含量是在550 ℃下焚烧4 h后所得,每个样品重复测定2次。水分含量和灰分含量计算公式如下,

      X含水量= M0+M-M1M × 100%, (1)

      式(1)中,X含水量:样品含水量;M0:称量瓶恒重后的质量;M:样品的质量;M1:称量瓶与样品恒重后的质量。

      X灰分含量=M1-M0M×(1-W)×100%, (2)

      式(2)中,X灰分含量:样品灰分含量;M1:坩埚和燃烧后残渣的质量;M0:空坩埚的质量;M:样品的质量;W:样品的含水量。

    • 采用偏最小二乘回归(PLS)对定量数据进行校正。PLS法是一种基于“精度矩阵(the concentration matrix)包含无用信息,分解精度矩阵,消除精度矩阵中无用信息”的回归方法。通过最小化误差平方和(sum of squares of errors)来确定目标数据的最佳匹配函数,然后用最简单的最小化误差平方和的方法来获得真值。

    • 模型A采用交叉验证均方根误差(RMSECV)评价模型的整体校正能力,由式(3)计算。

      RMSECV=i=1N(Y-X)2N-P, (3)

      式(3)中,N:样本数量;P:主成分数量;Y:预测值;X:真实值。

      模型B:预测均方根误差(RMSEP)是模型的整体预测能力评价指标,由式(4)计算。

      RMSEP=i=1N(Y-X)2N, (4)

      式(4)中,N:样本数量;Y:预测值;X:真实值。

    • 不同木薯品种的水分和灰分含量(表1)的对比分析结果表明,在不同品种的淀粉灰分含量中,14号(老板娘2号)含量最低,为1.03%,91号(50)含量最高,为4.26%,参试品种的灰分含量平均值为2.00%。通过对木薯不同品种的水分含量分析得知,49号(南美9号)含量最低,为1.57%,94号(417)含量最高,为8.45%,参试品种的水分含量平均值为4.54%。

      表 1  不同木薯品种的水分和灰分含量

      Table 1.  Water and ash contents of cassava Germplasm Resources

      编号No.品种Name灰分/%Ash/%水分/%Water/%编号No.品种Name灰分/%Ash/%水分/%Water/%
      1F20002.78 ± 0.044.24 ± 0.2528SC2011.98 ± 0.062.95 ± 0.33
      2东莞红尾1.94 ± 0.044.19 ± 0.3729F2022.14 ± 0.053.60 ± 0.16
      3Q101.89 ± 0.004.86 ± 0.3130F03P2.34 ± 0.024.50 ± 0.35
      4JG13012.53 ± 0.033.65 ± 0.3931COL7122.06 ± 0.044.71 ± 1.39
      5SC92.14 ± 0.064.81 ± 0.2532SC21.71 ± 0.014.14 ± 0.86
      6MF-元引9831.19 ± 0.023.34 ± 0.8133CH201.39 ± 0.014.87 ± 1.27
      72B1.75 ± 0.023.58 ± 2.1334SC62.04 ± 0.065.13 ± 1.44
      8老板娘3号1.68 ± 0.015.21 ± 0.0735COL7151.93 ± 0.054.95 ± 0.65
      9GR32.24 ± 0.015.05 ± 0.3336T122.00 ± 0.023.39 ± 0.57
      10ME1911.18 ± 0.035.43 ± 0.283782292.07 ± 0.024.29 ± 0.78
      11C31.29 ± 0.033.14 ± 2.5238ITBB011.97 ± 0.033.90 ± 0.30
      12泰国1.15 ± 0.004.73 ± 0.0839KU582.14 ± 0.054.02 ± 0.24
      13云南思茅1.72 ± 0.003.56 ± 0.0440ME1902.26 ± 0.053.31 ± 0.05
      14老板娘2号1.03 ± 0.033.22 ± 0.7841CM39921.97 ± 0.054.18 ± 0.66
      15GR9111.95 ± 0.044.46 ± 0.3042SC2051.72 ± 0.004.31 ± 0.10
      16KU501.83 ± 0.033.90 ± 0.4743C2221.74 ± 0.005.92 ± 3.29
      17SC124多倍体1.84 ± 0.063.27 ± 0.9844SC101.64 ± 0.023.47 ± 1.24
      187751.72 ± 0.002.36 ± 1.1845F2972.06 ± 0.024.74 ± 0.17
      19Dec-461.16 ± 0.034.45 ± 0.7146水果木薯1.72 ± 0.013.87 ± 0.61
      20SC162.03 ± 0.044.47 ± 0.5247SC1241.64 ± 0.013.20 ± 0.11
      21哥伦比亚18R2.10 ± 0.052.43 ± 0.0548SC122.16 ± 0.044.42 ± 2.23
      22SC41.41 ± 0.051.89 ± 1.2149南美9号2.25 ± 0.041.57 ± 0.42
      23利民1.73 ± 0.015.46 ± 2.8650花叶木薯1.93 ± 0.044.21 ± 0.28
      24罗勇801.23 ± 0.003.30 ± 0.4551罗勇901.72 ± 0.017.07 ± 0.27
      25COL7131.05 ± 0.023.54 ± 0.4052宝岛9-51.96 ± 0.015.78 ± 0.01
      26SC51.12 ± 0.034.47 ± 0.2653E13401.52 ± 0.016.08 ± 0.74
      27GR61.91 ± 0.044.16 ± 1.0054ME1891.39 ± 0.014.78 ± 0.55
      续表1 Tab.1 Continued
      编号No.品种Name灰分/%Ash/%水分/%Water/%编号No.品种Name灰分/%Ash/%水分/%Water/%
      55NZ1881.84 ± 0.065.76 ± 0.43921963.67 ± 0.334.04 ± 0.11
      56琼中1号1.63 ± 0.044.71 ± 0.21935714.24 ± 0.193.86 ± 0.16
      57泰国种1.57 ± 0.025.82 ± 0.00944173.05 ± 0.078.45 ± 0.21
      58新选0481.65 ± 0.014.22 ± 1.9495ZM82294.12 ± 0.164.38 ± 0.21
      59H6601.54 ± 0.055.3 ± 0.39965213.88 ± 0.004.17 ± 0.08
      60SC8辐射2.41 ± 0.023.55 ± 2.6497SC80133.06 ± 0.083.33 ± 0.15
      61NZ1991.23 ± 0.034.01 ± 1.69982743.27 ± 0.074.52 ± 0.07
      62ME1871.89 ± 0.026.27 ± 0.29994284.20 ± 0.043.51 ± 0.05
      632761.78 ± 0.033.86 ± 1.0010060683.00 ± 0.056.42 ± 0.22
      6418R2.09 ± 0.034.85 ± 0.11101MB3.69 ± 0.143.61 ± 0.14
      65COC1732.11 ± 0.063.49 ± 2.23102T3B12.35 ± 0.066.1 ± 0.04
      66C91.39 ± 0.015.22 ± 0.51103T2B1.46 ± 0.036.24 ± 0.00
      67M-SC1241.77 ± 0.065.64 ± 0.01104T1A12.63 ± 0.045.20 ± 0.21
      68MF5561.25 ± 0.013.1 ± 2.34105T4C12.42 ± 0.035.47 ± 0.14
      69老板娘1号1.16 ± 0.036.39 ± 1.2106T2A11.43 ± 0.036.01 ± 0.05
      70MF811.57 ± 0.044.91 ± 0.04107T3A11.98 ± 0.055.08 ± 0.07
      71F5391.65 ± 0.014.95 ± 0.46108T4B11.32 ± 0.025.59 ± 0.07
      72M-GR9111.65 ± 0.043.35 ± 1.16109T1C11.92 ± 0.055.01 ± 0.17
      73GR8911.97 ± 0.034.24 ± 0.13110T4A11.26 ± 0.035.32 ± 0.35
      74MF6481.41 ± 0.024.19 ± 0.75111T3C13.36 ± 0.045.35 ± 0.14
      75MF6711.29 ± 0.025.10 ± 0.30112T2C12.66 ± 0.054.76 ± 0.28
      76SC72.92 ± 0.044.98 ± 0.50113T1B1.96 ± 0.045.18 ± 0.04
      77M 糯米1.49 ± 0.044.55 ± 0.521141 崖州扁豆左62.36 ± 0.035.46 ± 0.21
      78C5761.84 ± 0.034.34 ± 1.471151 撒播黑豆右33.08 ± 0.024.33 ± 0.07
      79C10301.56 ± 0.057.70 ± 0.901161 柱花草左 72.67 ± 0.013.72 ± 0.14
      80MF5321.99 ± 0.034.29 ± 0.931171 崖州扁豆左52.73 ± 0.094.22 ± 0.28
      81SC81.64 ± 0.034.50 ± 0.631181 花生左62.29 ± 0.064.53 ± 0.32
      82MF3591.19 ± 0.016.13 ± 1.751191 花生右13.26 ± 0.014.36 ± 0.24
      83M-FC101.67 ± 0.013.50 ± 0.691201 对照2右31.99 ± 0.052.41 ± 0.08
      84MSC51.07 ± 0.023.94 ± 0.931211 崖州扁豆右41.99 ± 0.023.61 ± 0.08
      85E4071.77 ± 0.003.61 ± 0.961221 柱花草右31.85 ± 0.033.25 ± 0.13
      86MF1811.12 ± 0.024.02 ± 1.171231 崖州扁豆右61.20 ± 0.023.16 ± 0.11
      87M-SC41.41 ± 0.012.47 ± 0.011241 油菜花3左12.72 ± 0.053.37 ± 0.23
      88OOO1.25 ± 0.013.96 ± 0.421251左撒崖州扁豆1.78 ± 0.034.24 ± 0.09
      89SC111.79 ± 0.005.68 ± 0.021261黑豆撒左42.59 ± 0.083.59 ± 0.35
      90F2742.26 ± 0.024.10 ± 1.681271花生左42.48 ± 0.013.17 ± 0.20
      91504.26 ± 0.074.07 ± 0.071281对照左12.04 ± 0.063.61 ± 0.23
      续表1 Tab.1 Continued
      编号No.品种Name灰分/%Ash/%水分/%Water/%编号No.品种Name灰分/%Ash/%水分/%Water/%
      1291崖州扁豆撒右22.61 ± 0.072.56 ± 0.401341 油菜间作撒右12.26 ± 0.066.35 ± 0.13
      1301柱花草2右22.40 ± 0.026.44 ± 0.281351柱花草1左22.20 ± 0.027.63 ± 0.10
      1311柱花草右51.90 ± 0.036.39 ± 0.081361右8对照1.06 ± 0.027.14 ± 0.09
      1321撒播柱花草右42.29 ± 0.016.49 ± 0.071371对照左81.07 ± 0.018.15 ± 0.11
      1331花生右71.17 ± 0.047.56 ± 0.57
      注:崖州扁豆、柱花草、花生、油菜花分别代表与木薯SC9品种间作的不同植物种类。
      Note:Macrotyloma uniflorum ‘Yazhou’,stylosanthes, peanut and rape flower represent different plant species intercropping with cassava SC9 variety respectively
    • 在4 000~12 000 cm-1范围内计算了137个样品的平均光谱。由原始光谱结果(图1)可知,137个样品的吸收光谱变化趋势相似。但由于木薯成分复杂,不同碱基对应不同组分(如甲基和亚甲基)或相同碱基在不同化学环境下的近红外吸收波长和强度有显著差异,因此木薯各成分的光谱吸收强度不同。图1分布在9 000~12 000 cm-1波段范围内的谱线明显可分为两类,由此表明根据化学成分和分子结构可将137个样品分为两类。

      图  1  参试木薯种质资源近红外光谱图

      Figure 1.  Near infrared spectrum of cassava Germplasm Resources

    • 在TQ Analyst 9.0分析软件中,采用偏最小二乘法(PLS)对137份木薯样品进行光谱离群值剔除(Spectrum Qutlier)分析,未发现存在影响定量模型建立的异常样本[20],同时根据软件优化选择测定灰分和水分2个指标的校正集与验证集,保证将最小和最大含量值的样本划入校正集中,具体划分见表2

      表 2  校正集与验证集样本信息

      Table 2.  Sample information of correction set and verification set

      成分Component校正集Correction set验证集Validation set
      样本量sample size范围Range样本量sample size范围Range
      灰分/%Ash/%1171.030~4.260201.150~4.200
      水分/%Water/%1171.570~8.450202.470~7.140
    • 对木薯样品光谱数据在测定波长范围内进行总体方差计算,以每一波长点的方差值[D(X)]对波长作图(图2),观察光谱数据变异性大小[21]。由图2可见,在6 184.92~4 309.97 cm-1和3 750.18~3 594.66 cm-1波段内曲线波动明显,数据变异性强,可见在此区间内的信息量丰富。以此波段为特征区间,考察全光谱波段(用a表示)与此特征区间(用b表示)的建模效果。

      图  2  木薯样品光谱数据方差分析图

      Figure 2.  Analysis of variance chart

      为了将全光谱波段与模型优化预测波段进行比较,确定最佳建模波段,分别考察了全光谱波段和光谱波段特征区间(5 800~4 300cm-1和4 000~3 600 cm-1)的建模效果(表3)。检测MSC、1S、2S、SG和SNV 5种预处理方法对定量建模影响,以校正模型标准差(RMSEC)、预测均方根误差(RMSEP)、交叉验证均方差(RMSECV)和相关系数(R)作为模型的评价指标,结果表明,对于水分和灰分,经 MSC+2S+SG预处理后,在选择b波长区间分别建立模型时,R值最大,RMSEC、RMSEP以及RMSECV接近于0。这说明经MSC+2S+SG预处理后,选择b为建模区间得到的模型预测结果更准确[22]

      表 3  不同预处理以及建模区间的考察结果

      Table 3.  Investigation results of different pretreatment and modeling intervals

      指标Index预处理Pretreatment区间Section相关系数R校正标准差RMSEC预测标准差RMSEP交叉标准差RMSECV
      灰分AshConstant(原始)a0.742 80.457 00.636 00.514 0
      b0.759 80.444 00.591 00.498 0
      MSCa0.769 40.436 00.403 00.509 0
      b0.774 40.432 00.675 00.575 0
      SNVa0.770 80.435 00.404 00.518 0
      b0.729 40.467 00.633 00.579 0
      MSC+1S+SGa0.826 390.384 00.636 00.515 0
      b0.904 00.292 00.521 00.535 0
      SNV+1S+SGa0.828 80.382 00.628 00.513 0
      b0.906 20.289 00.554 00.495 0
      MSC+2S+SGa0.873 50.332 00.820 00.590 0
      b0.944 90.223 00.210 00.398 0
      SNV+2S+SGa0.876 00.329 00.810 00.623 0
      b0.925 70.270 000.247 00.495 0
      水分WaterConstant(原始)a0.788 80.784 01.280 01.440 0
      b0.693 90.919 01.300 01.440 0
      MSCa0.778 00.801 01.340 01.280 0
      b0.685 50.929 01.280 01.420 0
      SNVa0.784 10.792 01.380 01.440 0
      b0.696 20.916 01.400 01.470 0
      MSC+1S+SGa0.841 40.689 01.450 01.500 0
      b0.820 70.729 01.220 01.530 0
      SNV+1S+SGa0.844 80.683 01.450 01.480 0
      b0.884 70.595 01.370 01.640 0
      MSC+2S+SGa0.863 50.643 01.710 01.620 0
      b0.926 40.480 01.460 01.540 0
      SNV+2S+SGa0.862 80.645 01.710 01.620 0
      b0.892 60.575 01.340 01.590 0
      注:区间a:全光谱波段;区间b:光谱波段特征区间(5 800~4 300 cm-1和4 000~3 600 cm-1)。
      Note: Interval a: full spectrum band; Interval b: spectral band characteristic interval (5 800~4 300 cm-1 and 4 000~3 600 cm-1
    • 根据上述最佳光谱预处理法和最优建模区间,采用PLS建立灰分和水分的近红外定量分析模型。最优定量模型均为通过MSC+2S+SG处理,以5 800~4 300 cm-1和4 000~3 600 cm-1为建模区间建立而成。通过RMSECV对主因子数作图获得最佳主因子数分别为14、16;模型相关系数(R)分别为0.944 9、0.926 4;RMSEC依次为0.223 0、0.480 0;RMSEP依次为0.210 0、1.460;RMSECV依次为0.398 0和1.540。R值均接近于1,表明预测值与参考值间存在良好的相关性(图3)。

      图  3  木薯样品预测值和参考值的相关图

      Figure 3.  Correlation diagram between predicted value and reference value of cassava sample

    • 为进一步研究灰分与水分2种优化模型的可预测性,分析偏最小二乘法拟合得到验证集指标值与真实值的拟合程度,在95%置信区间内进行配对t检验,分别检验最优模型的拟合值,验证定量模型可靠性,可在实际检测中实时反馈木薯品质。选用未参与建模的20份木薯种质资源对该模型进行外部验证,预测值与真实值进行差异性分析,(P>0.05),P值分别为0.464、0.459,说明差异无显著性。由此表明,实验建立的预测模型可实现对不同产地木薯的灰分和水分的快速无损检测。

    • 由于化学成分在不同波长处吸收强度差异,不同波长段所蕴藏的信息差异明显,选择合理波长范围和预处理方法可有效地提高模型性能。近红外光谱(NIRS)与多变量分析相结合已被用于测定木薯块根中的类胡萝卜素、纤维素、氰化物、干物质、淀粉、总糖和总氮含量等多种营养成分[23-25],但对灰分和水分含量的测定模型构建仍然有待研究,本研究结果可以对其进行补充完善。李荣[26]在对波长5 296~12 000 cm-1内益智水分和灰分含量进行NIRS定量分析模型构建后发现,MSC+2S+SG和SNV+2S分别为益智水分和灰分建模的最佳预处理方法,这与本实验结果相类似[27-29]。定标模型建好后,检测木薯样品中的水分、灰分的含量仅需光谱扫描和软件分析2步即可完成。此外,利用近红外方法进行分析测定,所得结果较为准确,且稳定性较高,可以有效减少样品的损耗,具有较好的应用前景[30-33]

      本研究采集的木薯样品均来自同一产地,所测灰分和水分含量分布范围广,具有一定的代表性,可用于分析模型的建立。考虑到木薯产地分布广泛,种属较多,定量分析模型在实际应用过程中需补充新的样品对模型进行完善,以扩大模型的适应性。

    • 本研究建立了新型的近红外光谱分析新方法,能够成功用于快速检测木薯水分、灰分的含量。结果显示,水分和灰分的最优定量模型均为通过MSC+2S+SG处理,以5 800~4 300 cm-1和4 000~3 600 cm-1为建模区间建立而成。预测值与真实值之间无显著差异,由此表明2种方法之间不存在明显系统误差。

参考文献 (33)

目录

    /

    返回文章
    返回