王 凡 李永玉 彭彦昆 杨炳南 李 龙 尹学清
(1.中国农业大学工学院,北京 100083; 2.中国农业机械化科学研究院,北京 100083)
马铃薯是我国仅次于小麦、水稻和玉米的第四大主要粮食作物。近些年来,随着我国与周边国家经贸合作关系的不断深化,加工特性优良的鲜马铃薯和薯条、薯片及马铃薯淀粉出口将继续保持增加态势[1]。马铃薯内部成分不仅决定其营养价值,而且还会影响加工品质。马铃薯中水分占块茎总质量的70%左右,主要影响油炸过程中的耗油量、加工中蒸发水所需的能耗及烹饪口感;淀粉占马铃薯块茎总质量的10%~20%,是马铃薯最主要的营养成分之一;还原糖占马铃薯块茎总质量的0.1%~5%,它与氨基酸产生的褐变反应会导致马铃薯产品色泽灰暗,甚至产生有害物质。马铃薯品质参差不齐是制约我国马铃薯产业主食化进程的重要因素之一,马铃薯品质快速无损检测对推进马铃薯产业发展有着重要意义。
目前,国内外不少学者对马铃薯内部品质进行了无损检测研究,WALSH等[2]、KANG等[3]使用透射光谱分别建立了干物质、比重的近红外预测模型,相关系数分别为0.79和0.85,HAASE等[4]、SUBEDI等[5]、张小燕等[6]对马铃薯干物质、含水率等主要参数进行了预测,相关系数均在0.9以上。但这些研究均属于检测方法实验研究,检测系统体积大、移动不便,不利于实际推广和应用。随着近红外光谱在农产品品质检测领域的推广应用,我国学者设计了各种农产品专用便携式检测装置,如生鲜肉品质检测装置[7-8]、作物和水果品质检测装置[9-11]等。但是有关马铃薯品质检测便携式装置的研究在国内未见报道,随着我国马铃薯主食化产业的推进,研究适用于马铃薯产地的马铃薯品质便携式实时检测装置显得尤为重要。
本文基于可见/近红外局部透射光谱,根据马铃薯大小及形状特征,设计便携式马铃薯多品质无损检测装置,编写实时检测分析控制软件,植入马铃薯含水率、淀粉质量分数、还原糖质量分数的预测模型,实现马铃薯内部品质的实时无损检测。
当近红外辐射照射到样品表面时,一部分辐射发生反射,一部分辐射透过样品从另一侧射出,而另一部分辐射被样品散射并返回到表面。反射光谱信号强度较好,但受马铃薯表皮影响较大[12],透射光谱能较好地反映内部信息,但因马铃薯大小引起的光程差异较大[13-14]。考虑到马铃薯整体质地较为均匀,本研究采用局部透射的检测方式。整个装置由光谱采集模块、光源模块、控制与显示模块等构成,如图1所示。光源模块由灯珠与反光杯组成,为整个装置提供光照;光谱采集模块由光谱仪和聚焦透镜组成,该部分将透过样品的光收集、转换为电信号并传输给ARM处理器,处理器将处理后的光谱和预测值显示在屏幕上,整个装置由供电模块供电。
图1 硬件系统结构图Fig.1 Schematic diagram of hardware system
1.2.1光谱采集模块
光谱仪是本设备的核心部件,它的优劣很大程度上决定着装置的精度和稳定性。本研究通过多种光谱仪的体积、原理、结构、分辨率、波段范围等对比分析,根据预实验并参照文献 [15-18],选择了波长范围为650~1 100 nm的海洋光学公司生产的STS光谱仪。它是基于CMOS(互补金属氧化物半导体)传感器的微型光谱仪,在结构上采用交叉式Czerny-Turner 结构[19],拥有较小体积(40 mm×42 mm×24 mm),具有较高的信噪比(1 500以上)和分辨率(1.5 nm)。
图2 耦合透镜示意图Fig.2 Diagrams of coupling lens1.聚焦透镜 2.螺纹扣环 3.螺纹准直筒
在大部分近红外光谱检测系统中,光谱仪与样品之间由光纤连接[20-23],但研究证明光在光纤中传导时会损失一部分能量,这会导致积分时间的增加,而且光纤脆性较大,为保证其最小弯折半径,往往需要更大的操作空间。光谱仪虽然可以不通过光纤直接接收光信号,但光谱仪接口暴露会导致灰尘进入,且光谱仪体积较大,直接作为探头,布置较为困难。基于以上考虑,本研究选择海洋光学74-DA型耦合透镜作为检测探头。耦合透镜结构示意图如图2a、2b所示,聚焦透镜由两个带螺纹的扣环固定在螺纹准直筒内,为了验证耦合透镜的效果,使用3种不同连接方式,以100 ms为积分时间采集白参比的光谱,结果如图2c所示。3种不同方式采集的光谱曲线形状大致相似,但使用聚焦透镜可以明显提升光谱信号强度。
1.2.2光源模块
卤素灯是近红外光谱检测中常用的光源,其发射波长涵盖200~2 500 nm,为了选择合适功率的光源,分别测试了0.5、1.0、1.5、2.2 W的灯珠在100 ms下的样品光谱响应。如图3a所示,功率小于1 W的光源发光较弱,样品光强占光谱仪满量程光强(14 000 cd)不足40%,信噪比差,容易受到光谱仪噪声的干扰,当光源功率为2.2 W时,因光源与样品距离较近,样品表皮出现略微灼伤。本文采用了1.5 W(3 V、0.5 A)998070-1型超高亮卤素灯(美国Welch Allyn公司)加反光杯的集成化设计,如图3b所示。
图3 光源模块示意图Fig.3 Schematic of light source1.光谱仪 2.耦合透镜 3.样品 4.反光杯 5.灯珠
1.2.3控制与显示终端
控制单元不仅要实现对光谱仪的控制,还要对采集的数据进行计算处理并传输到显示终端,因此控制单元需要较快的运行速度和合适的外接接口。本研究选择Raspberry Pi 3b作为控制单元,它是一款基于ARM的微型计算机主板,以SD/MicroSD卡为内存硬盘,卡片主板周围有4个USB、1个100 MB以太网接口和HDMI高清视频输出接口,可以实现大量数据流畅运算,但其尺寸只有85 mm×56 mm×17 mm,满足设备小型化的需求。
根据控制单元,选择了与之匹配的树莓派3.2英寸电阻屏(3.2inch RPi LCD (B),85 mm×56 mm×17 mm),该屏幕分辨率为320像素×240像素,与树莓派之间通过I/O接口相连,除了显示功能,该屏幕附带4颗按钮,在本设备中利用其中3颗按钮分别作为黑参考、白参考和样品光谱的采集外触发。
1.2.4供电模块
为了满足现场检测的需求,需要配备可充电式电池。光谱仪、控制器、显示器的输入电压为5 V,光源的输入电压为3 V,选用输出电压为5 V的可充电锂电池,一路直接为控制器供电,另一路经稳压板(LTC3780为芯片)转换为3 V电压后为卤素灯供电。经计算选择电量为6 A·h的锂电池,可保证装置续航时间达到3 h以上。
1.2.5整机设计
基于对光源单元、光谱采集单元、开发板以及稳压板大小的综合考虑,将开发板沿壳体的长度方向横向设置,然后将光谱采集单元和光源竖向设置,电源和光源之间设置稳压板,在整机右侧设置电池,达到最大限度节省空间的目的。便携式马铃薯多品质无损设备如图4所示,该装置整机尺寸11 cm×6.5 cm×9.5 cm,可单手操作。
图4 便携式马铃薯多品质无损检测设备Fig.4 Images of potato hand-held equipment1.ARM处理器-屏幕 2.充电接口 3.开关 4.光谱仪 5.耦合透镜 6.反光杯 7.灯珠 8.稳压板 9.电池 10.按钮
1.2.6黑白参考盒设计
光谱仪的工作状态及电源性能等均会影响近红外光谱的稳定性,为消除这些外部因素对检测结果的影响,采集光谱时一般先进行黑白参考的校正。因便携式装置无支架、暗箱等,设计了黑白参考校正盒,如图5a所示。校正盒分为盒体和盒盖,盒体顶面设置通孔,通孔的形状恰好与探头外轮廓相匹配形成暗室,盒盖防止灰尘落入并起保护作用。参考盒的内部底面设置材料为聚四氟乙烯的白板,如图5b所示。
图5 参考盒结构图Fig.5 Structure diagrams of reference box
图6 马铃薯局部透射光谱Fig.6 Partial transmission spectrum of potato
选用由中国农业科学院花卉蔬菜研究所(北京)提供的大西洋品种马铃薯样品85个,其中65个用于建模,20个用于验证。采摘后送往实验室4℃条件下冷藏待用,采集光谱前将马铃薯块茎洗净置于室温(20℃)中放置24 h(相对湿度60%)以消除温度对实验的影响。
利用自行设计的便携式马铃薯多品质检测装置采集所有马铃薯样品的局部透射光谱。采集前将装置预热10 min,采集黑、白参考并保存,每个马铃薯样品避开芽眼和缺陷在赤道上间隔90°取4个采集点,4点的平均光谱作为该样品的原始光谱,并换算成马铃薯样品局部透射率,公式为
式中T——马铃薯样品的局部透过率,%
I——马铃薯样品的局部透过光谱强度,cd
D——黑参考光谱强度,cd
R——白参考光谱强度,cd
光谱采集完成后,分别用直接干燥法[23]、酸解法[24]和滴定法[25]测定马铃薯样品的含水率、淀粉质量分数和还原糖质量分数。
马铃薯样品的原始光谱先用五点S-G卷积平滑(Savitzky-Golay smooth, SG)进行平滑处理,因实际检测过程中存在样品表皮散射及光程的变化,对SG平滑后的光谱进行标准正态变量变换(Standard normal variable transformation, SNV),原始光谱曲线、SG平滑预处理及SNV处理后的光谱曲线如图6所示。
图7 测量值和预测值对比Fig.7 Comparison charts of measurements and predictions
基于SG-SNV预处理后的光谱曲线分别建立了马铃薯含水率、淀粉质量分数、还原糖质量分数的偏最小二乘预测模型(Partial least squares regression, PLSR),模型结果如图7所示。含水率、淀粉质量分数和还原糖质量分数的校正集相关系数Rc分别为0.937 9、0.913 7和0.937 1,校正集均方根误差(Root mean square error of calibration, RMSEc)分别为0.316 1%、0.309 0%和0.034 8%;验证集相关系数Rv分别为0.927 8、0.914 6和0.933 8,验证集均方根误差(Root mean square error of validation, RMSEv)分别为0.325 3%、0.344 9%和0.041 6%。基于自行设计的马铃薯多品质无损检测装置所建立的马铃薯含水率、淀粉质量分数和还原糖质量分数3个参数的预测模型,其预测精度与相关的近红外漫透射光谱研究[2-3]有所提升,能满足实际现场检测需求。
图9 装置测试结果散点图Fig.9 Model result of external prediction
基于QT开发工具,采用C/C++语言编写了检测装置实时分析控制软件,软件采用模块化思想,分为光谱仪控制模块、光谱处理模块和结果显示模块,光谱仪控制模块对光谱仪采集参数(如积分时间、平均次数等)进行设置并采集光谱数据;光谱处理模块自动对光谱数据进行透过率的转换和预处理,并通过内置的模型预测样品含水率、淀粉质量分数和还原糖质量分数;结果显示模块将预处理后的光谱和品质预测值显示在屏幕上。软件工作流程如图8所示。采集样品光谱前,先点击“Dark”和“White”按钮分别对光谱仪进行黑白校正并保存校正结果。探头对准样品,采集样品的局部透射光谱,自动进行SG平滑和SNV预处理并分别代入植入模型进行马铃薯各品质参数的预测,样品光谱曲线和检测结果实时显示在LCD屏幕。
图8 便携式装置软件工作流程图Fig.8 Working flow chart of software in hard-held device
选取与建模样本无关的大西洋品种马铃薯样品20个,对便携式马铃薯多品质无损检测装置的稳定性和预测精度进行了验证。利用检测装置预测马铃薯样品含水率、淀粉质量分数和还原糖质量分数,每个样品采集4次,以4次的变异系数衡量仪器的稳定性,再用化学分析方法检测马铃薯3个品质参数的标准理化值,并进行对比分析,装置预测值与标准理化值散点图如图9所示。马铃薯样品的含水率、淀粉质量分数和还原糖质量分数的预测值与标准理化值的相关系数Rp分别为0.914 1、0.912 2和0.914 0,预测均方根误差分别为0.352 7%、0.340 4%和0.040 0%,平均偏差(Average deviation, AD)分别为0.295 1%、0.253 6%和0.031 6%。每个样品重复采集4次光谱,分别预测含水率、淀粉质量分数和还原糖质量分数,其最大变异系数分别为0.006 7、0.012 4和0.123 1。结果表明,便携式马铃薯多品质无损检测装置稳定性和检测精度均能满足现场实时检测要求。
(1)基于可见/近红外局部透射光谱,根据马铃薯大小及形状特征,设计了便携式马铃薯多品质无损检测硬件装置,包括光谱采集模块、光源模块、控制与显示模块、供电模块,并设计了专用黑白参考校正盒。整个装置尺寸为11 cm×6.5 cm×9.5 cm,易于携带。
(2)基于便携式马铃薯多品质无损检测装置,建立了马铃薯多品质预测模型。采集马铃薯650~1 100 nm的可见/近红外局部透射光谱,对原始光谱进行SG平滑和SNV处理后,建立了马铃薯含水率、淀粉质量分数和还原糖质量分数的偏最小二乘预测模型。马铃薯含水率、淀粉质量分数和还原糖质量分数预测模型验证集相关系数分别为0.927 8、0.914 6和0.933 8,验证集均方根误差分别为0.325 3%、0.344 9%和0.041 6%。
(3)基于QT开发工具,采用C/C++语言编写了装置实时分析控制软件。将所建马铃薯含水率、淀粉质量分数和还原糖质量分数预测模型植入到硬件装置中,实现了马铃薯多品质参数实时无损检测一键式操作。
(4)验证了便携式马铃薯多品质无损检测装置的检测精度。马铃薯样品含水率、淀粉质量分数和还原糖质量分数的装置预测值与标准理化值的相关系数分别为0.914 1、0.912 2和0.914 0,均方根误差分别为0.352 7%、0.340 4%和0.040 0%,平均偏差分别为0.295 1%、0.253 6%和0.031 6%,重复采样最大变异系数分别为0.006 7、0.012 4和0.123 1。结果表明,该装置可以实现马铃薯含水率、淀粉质量分数和还原糖质量分数的实时无损检测。