Research on peanut hyperspectral image classification method based on SPA-PSO-BP
-
摘要: 为了提高可见-近红外(VNIR)高光谱花生图像分类的准确率和减少分类检测的运算时间，提出了基于连续投影算法(SPA)融合粒子群算法优化后向传播神经网络(PSO-BP)的分类检测模型。利用高光谱成像系统采集了7个花生品种样本的VNIR光谱数据，并进行了背景分割和光谱信息的提取，去除受噪声和杂散光影响大的波段后，运用Savitzky-Golay卷积平滑对400 nm~900 nm范围的波长进行预处理；采用SPA降维及均方根误差值选择了25个特征波长，同时利用PSO-BP神经网络的初始权重和阈值，构建PSO-BP模型作为分类器进行了实验，取得了测试集识别准确率为98.7%、kappa系数为0.98及遗漏误差为3的数据。结果表明，相较4个对比方法构建的分类模型，该模型的准确率分别提高了2.1%、8.6%、3.9%和4.3%。该方法在基于高光谱成像的花生品种分类技术中具有很好的应用前景，为花生品种的高精度、快速无损分类提供了新思路。Abstract: In order to improve the accuracy of visible-near infrared (VNIR) hyperspectral peanut image classification and to reduce the computing time of classification detection, a classification detection model based on successive projection algorithm (SPA) fused with particle swarm optimization back propagation(PSO-BP) neural network was proposed. A hyperspectral imaging system was used to acquire VNIR spectral data of seven peanut species samples and conducts background segmentation and extraction of spectral information. After removing the wavelengths that were highly affected by noise and stray light, the wavelengths in the range of 400 nm~900 nm were preprocessed by using Savitzky-Golay convolutional smoothing. The SPA was used to reduce the dimensionality, and 25 characteristic wavelengths were selected by virtue of the root mean square error values. The PSO was also used to optimize the initial weights and thresholds of the BP neural network, and the PSO-BP model was constructed as a classifier for the experiments, and a recognition accuracy of 98.7%, a kappa coefficient of 0.98, and a miss error of 3 for the test set were obtained. The results demonstrate that the accuracy of the model is improved by 2.1%, 8.6%, 3.9%, and 4.3%, respectively, compared with the classification models constructed by the four comparison methods. The proposed method has good application prospects in peanut variety classification based on hyperspectral imaging technology, and provides a new idea for high accuracy and fast nondestructive classification of peanut varieties.
-
0. 引言
激光辐照气流环境中典型金属材料的热力学响应是激光加工等应用领域关心的热点和关键问题[1-4]，铝合金是工业领域常用材料，研究获得其热响应和熔穿特性，对激光切割、打孔等实际应用具有重要指导意义和价值[5-7]。
在激光辐照材料过程中，材料响应与激光体制参数、材料热物理性质、激光与材料耦合特性、辐照区域内材料所处气流等环境息息相关[8-11]。当发生激光烧蚀材料相变后，金属的相变熔融物会阻碍激光向金属材料内部辐照，将熔融物从熔池中排出需要消耗部分激光能量，导致激光烧蚀清除效率低[12-14]。国外从20世纪70年代已开始从理论上研究切向气流对激光作用下材料清除时间的影响，近年来，国内外在实验和理论方面已有较多气流速度、气流组分等作用金属、复合材料靶的研究，结果表明，当切向气流作用于靶表面时，由于对流换热和切向气流对熔蚀物的清除作用，材料热响应和清除特性与自然对流环境下相比存在显著差异[8, 12, 15-16]。甚至当激光加热使靶材发生软化时，气流引起的前后表面压力差导致的剪切应力可在远低于熔点的情况下实现材料清除[17-18]。虽然理论、实验已有较多关于气流环境下激光辐照研究，但仍以定性的结果和规律认识为主。应用中如何高效、定量地获得材料清除效果，仍需要发展满足物理需求的数值模拟方法。
本文作者构建了考虑主要物理机制的气流和激光联合作用材料热学响应物理模型，同时为量化比较对流换热和清除熔蚀物二者在激光辐照材料中的作用，以铝合金为例，计算获得侧向气流条件不同激光功率密度下铝合金板的烧蚀清除规律，相关工作对铝合金板加工参数选择具有一定参考价值。
1. 数值计算模型
1.1 物理模型
气流环境下激光辐照材料是一个复杂的多物理场强耦合物理过程，存在气流换热、热-力耦合等多种物理效应，激光加载的能量、气流流动引起对流换热导致的能量损失以及材料烧蚀导致能量交换将共同决定靶材内部温度的分布，需要建立相应的物理模型[19-21]。对于激光加载和能量耦合过程，鉴于不透明金属材料对激光吸收深度小，通常约10 nm量级，计算中激光加载能量作为面热源处理，在材料区域内，热传导方程的数学描述为：
ρc∂T∂t=∇⋅(κ∇T)+Qr (1) 式中：ρ为壳材料密度；c为壳材料比热容；T为材料温度；t为时间；∇为矢量微分算符；κ为材料导热系数；Qr为材料相变和氧化反应热源项。
激光加载区域边界条件为：
(κ∂T∂z)+qc+qr=q1 (2) 式中：z为靶厚度；qc为对流换热热流；qr为辐射换热热流；ql为材料表面吸收的激光热流。
加载对流边界换热热流可表示为：
qc=hc(T−T0) (3) 式中：hc对流换热系数；T0为表面气流温度。
辐射换热热流为：
qr=σε(T4−T04) (4) 式中：σ为Stefan-Boltzmann常数；ε为表面发射率。
假定靶温度1000 K，辐射换热热流功率密度约1 W/cm2量级，相对于对流和加载激光功率密度较小，计算中可忽略。激光辐照下金属表面高温物质脱落，激光束直接辐照下一层材料，使得辐照热效应增强。为了模拟这种熔蚀效应，将受激光辐照的最外层未剥离表面定义为熔迹面。热软化物质剥离采用温度准则判据，即达到某高温状态热软化材料被气流剥离[19, 21]。采用“单元生死法”计算边界的移动和热传导，当某一离散单元的温度超过熔化温度或气化温度，定义该单元不再参与计算，对应加载边界施加到新的单元上，整个过程不可逆。
1.2 模型参数及实验验证
通过计算结果与实验结果比较校验计算模型，模型校验用算例参数和数据编号见表 1。表中自然对流下实验结果case 1和case 2参数引自参考文献[16]。自然对流下激光辐照实验靶板温升模拟结果见图 1a，亚声速气流条件下实验的模拟结果见图 1b。图中，R为到靶板中心点距离。模拟计算的case 1和shot 1靶背面不同位置温升与实验结果符合较好，计算结果显示靶板未发生熔穿，靶板升温和辐照结束后自然冷却过程均与实验结果一致；实验测得的case 2和shot 2靶板熔穿时间Δt分别为6.5 s和6.0 s，采用计算铝板靶背面中心点温升获得靶板熔穿时间，模型计算的case 2和shot 2的熔穿时间分别为7.0 s和5.6 s，计算结果与实验结基本一致。通过温升曲线和熔穿时间比较，新建模型可较好地模拟激光辐照靶板的温度和熔穿过程。
表 1 模型校验用算例参数Table 1. Parameters of simulation case for model verificationnumber power density/(W·cm-2) irradiation time/s material material dimension/mm flow velocity/(m·s-1) case 1 1190 10 LY12 ∅30×2 0 case 2 1190 10 LY12 ∅30×1 0 shot 1 800 10 LY12 100×100×3 120 shot 2 1400 10 LY12 100×100×3 120 2. 结果与分析
为充分分析气流环境对清除或气化规律的影响，分别开展了气流条件和无气流条件下2种典型工况的计算分析：有气流环境采用切向速率240 m/s、静温20 ℃的无穷远处均匀空气来流，根据经验公式对流换热系数取0.06 W/(cm2·K)，以熔化温度作为穿孔(或单元删除)判据；无气流条件下以气化温度作为穿孔(或单元删除)判据计算，利用新建数值计算模型对激光辐照薄板厚度d =3 mm和厚板d =7 mm工况分别开展了熔穿规律数值模拟研究。计算中激光参数采用峰值功率密度I0分别取500 W/cm2、1000 W/cm2、1500 W/cm2、2000 W/cm2，光斑采用高斯分布；靶板采用LY12铝，几何模型厚度分别为3 mm和7 mm圆板，激光辐照在中心区域。
气流条件下厚度为3 mm靶板温度演化见图 2。加载激光峰值功率密度为500 W/cm2时，背面温度沿径向方向逐渐减低，背面靶板最高温度450 ℃，尚未达到铝的融化温度；当激光停止后，靶板在气体对流换热作用下，温度逐渐降低。计算获得不同峰值功率密度激光辐照靶板背面中心点的温度演化见图 2b。随着激光功率密度提高，靶温升率逐渐提高，其中在峰值功率密度1000 W/cm2时，辐照6.2 s后靶板熔穿；在峰值功率密度1500 W/cm2时，辐照3.4 s后靶板熔穿；在峰值功率密度2000 W/cm2时，辐照2.4 s后靶板熔穿。
![图 2 激光辐照3 mm厚度靶板背面的温度]()
图 2 激光辐照3 mm厚度靶板背面的温度Figure 2. Temperature at back of 3 mm target plate after laser irradiation气流条件下厚度为7 mm靶板温度演化见图 3。当激光峰值功率密度为2000 W/cm2时，靶板背面不同位置温度变化过程见图 3a，由于结构较厚，时间较长时，在距中心15 cm远位置也出现温升。当激光停止后，靶板在气体对流换热作用下，温度逐渐降低。激光功率密度为2000 W/cm2时板熔穿时间6.8 s。计算获得不同峰值功率密度激光辐照靶板背面中心点的温度演化见图 3b。随着激光功率密度提高，靶温升率逐渐提高，只有在峰值功率密度为2000 W/cm2时，辐照6.8 s后靶板熔穿；由于靶板较厚，在辐照10 s情况下其它较小加载功率密度均未造成熔穿。
![图 3 激光辐照7 mm厚度靶板背面的温度]()
图 3 激光辐照7 mm厚度靶板背面的温度Figure 3. Temperature at back of 7 mm target plate after laser irradiation在无气流条件、无重力等移除熔蚀物的作用情况下，计算获得的不同功率密度时3 mm厚度铝板背面中心点温度见图 4a，激光辐照的能量全部用于加热靶板，无质量迁移力，此时只有在铝板温度达到铝的气化温度后才能穿孔，计算结果显示，即使在功率密度为2000 W/cm2的高斯光束、辐照时间10 s，仍无法烧穿3 mm厚铝板；同样加载激光参数、靶参数情况，当施加质量迁移力，假定熔化后即将熔化层移除，此时3 mm厚度铝板背面中心点温升见图 4b，在功率密度1000 W/cm2~2000 W/cm2情况均出现熔穿。在不考虑气动加热等额外热源情况下，此时计算的熔穿时间即为该工况条件下的最快熔穿时间。
![图 4 熔蚀物清除力对靶背面温升的影响]()
图 4 熔蚀物清除力对靶背面温升的影响Figure 4. Effect of molten material removal force on temperture at back of target plate计算的气流和理想清除力情况下熔穿时间、功率密度参数的关系见图 5。图中，Δt为熔穿时间，A为材料表面吸收率。随着功率密度增加，熔穿时间急剧减小。与参考文献[17]中热平衡积分方法计算获得的气化烧蚀穿透模型结果比较，穿孔时间变化规律符合较好；但定量上由于二者物理模型存在差异，无法直接比较。本文中计算模型中光斑采用高斯分布模型，在考虑了侧向气流对流换热影响，与参考文献[17]中假定均匀光场定态气化假设不同，所以气流工况在低功率密度加载或厚板情况下与文献近似计算结果存在一定偏离；对于薄板、大功率密度加载情况下，对流换热作用时间短，与热平衡积分方法假设的条件更接近，并且本文中算例采用的光斑尺寸较大，在辐照区域内接近1维热扩散，所以参考文献中的理论与数值计算结果更接近。从穿孔时间可以看出，在功率密度较高区域各工况之间的差异小，功率密度较低区域各条件之间差异较大。分析认为，计算中声速固定时，对流换热系数则采用固定值，随着功率密度增加，对流换热的影响比例逐渐较降低。所以，在激光功率密度较小(500 W/cm2附近)时，近似计算或评估中建议不要忽略气流对流换热的影响；在激光功率密度较高(1500 W/cm2以上)时，气流对流换热的影响较小。
![图 5 熔穿时间与功率密度的关系]()
图 5 熔穿时间与功率密度的关系Figure 5. Relationship between melted perforation time and laser power density3. 结论
构建了考虑热传导、对流换热、熔化烧蚀等主要物理机制的气流和激光联合作用材料热学响应物理模型，获得气流环境下不同激光功率密度时铝合金板的热响应和烧蚀熔穿规律。高斯分布光束的出光峰值功率密度为1000 W/cm2、1500 W/cm2和2000 W/cm2时，3 mm厚铝板熔穿时间分别为6.2 s、3.4 s和2.4 s。随着功率密度增加，熔穿时间急剧减小，并且计算结果与参考文献中热平衡积分方法气化烧蚀穿透模型结果规律一致。与文献中假定均匀光场定态气化假设不同，本文中计算模型考虑了侧向气流对流换热影响，可实时更新温度扩散，数值模型与实际物理过程更接近。从穿孔时间可以看出，在功率密度较高区域各工况之间的差异小，功率密度较低区域各条件之间差异较大，在激光功率密度较小(500 W/cm2附近)时，近似计算或评估中建议不要忽略气流对流换热的影响；在激光功率密度较高(1500 W/cm2以上)时，气流对流换热的影响较小。由于气流环境下激光与金属相互作用包含复杂的传热、流动、化学反应以及金属材料的软化、脱落等力学过程，其综合作用机制还有待进一步研究。
-
![]()
图 6 黑花生典型样本预处理前后的光谱曲线
Figure 6. Spectral curves of Black peanut samples before and after pretreatment
![]()
图 7 彩衣仙子典型样本预处理前后的光谱曲线
Figure 7. Spectral curves of Colorful cloth peanut samples before and after pretreatment
![]()
图 8 冀花5号典型样本预处理前后的光谱曲线
Figure 8. Spectral curves of Jihua No.5 peanut samples before and after pretreatment
![]()
图 9 华一罗汉典型样本预处理前后的光谱曲线
Figure 9. Spectral curves of Huayiluohan peanut samples before and after pretreatment
![]()
图 10 双冠白沙典型样本预处理前后的光谱曲线
Figure 10. Spectral curves of Double crown white sand peanut samples before and after pretreatment
![]()
图 11 优品红沙典型样本预处理前后的光谱曲线
Figure 11. Spectral curves of Excellent red peanut samples before and after pretreatment
![]()
图 12 鲁花1号典型样本预处理前后的光谱曲线
Figure 12. Spectral curves of Luhua No.1 peanut samples before and after pretreatment
表 1 4种对比判别算法的参数设置、预测结果及kappa系数
Table 1 Parameter settings, prediction results and kappa coefficients of four comparative discriminative algorithms
recognition algorithm parameter omission error prediction accuracy/% kappa coefficient BP hidden layer: 10 8 96.6 0.96 RF trees: 60 leaf: 3 23 90.1 0.89 RBF spread: 100 12 94.8 0.94 SVM c: 10 kernel function: linear 13 94.4 0.94 
下载: 导出CSV
-
[1] 张智猛, 戴良香, 李美, 等. 花生种子产业现状与发展对策[J]. 中国农业科技导报, 2013, 15(1): 30-37. https://www.cnki.com.cn/Article/CJFDTOTAL-NKDB201301007.htm ZHANG Zh M, DAI L X, LI M, et al. Current situation and development countermeasures of peanut seed industry[J]. China Agricultural Science and Technology Bulletin, 2013, 15(1): 30-37(in Ch-inese). https://www.cnki.com.cn/Article/CJFDTOTAL-NKDB201301007.htm
[2] 李万伦, 甘甫平. 矿山环境高光谱遥感监测研究进展[J]. 国土资源遥感, 2016, 28(2): 1-7. https://www.cnki.com.cn/Article/CJFDTOTAL-GTYG201602001.htm LI W L, GAN F P. Research progress of hyperspectral remote sensing monitoring of mine environment[J]. Remote Sensing of Land and Resources, 2016, 28(2): 1-7(in Chinese). https://www.cnki.com.cn/Article/CJFDTOTAL-GTYG201602001.htm
[3] FENG Y Z, SUN D W. Application ofhyperspectral imaging in food safety inspection and control: A review[J]. Critical Reviews in Food Science & Nutrition, 2012, 52(11): 1039-1058.
[4] ZHU M, HUANG D, HU X J, et al. Application of hyperspectral technology in detection of agricultural products and food: A review[J]. Food Science & Nutrition, 2020, 8(10): 5206-5214.
[5] ZHOU Q, HUANG W Q, TIAN X, et al. Identification of the variety of maize seeds based on hyperspectral images coupled with convolutional neural networks and subregional voting[J]. Journal of the Science of Food and Agriculture, 2021, 101(11): 4532-4542. DOI: 10.1002/jsfa.11095
[6] JIN B Ch, ZHANG Ch, JIA L Q, et al. Identification of rice seed varieties based on near-infrared hyperspectral imaging technology combined with deep learning[J]. ACS Omega, 2022, 7(6): 4735-4749. DOI: 10.1021/acsomega.1c04102
[7] 刘璐, 邵慧, 孙龙, 等. 利用高光谱激光雷达检测木材的霉变与含水量[J]. 激光技术, 2023, 47(5): 620-626. DOI: 10.7510/jgjs.issn.1001-3806.2023.05.007 LIU L, SHAO H, SUN L, et al. Detection of mildew and moisture content in timber by hyperspectral LiDAR[J]. Laser Technology, 2023, 47(5): 620-626(in Chinese). DOI: 10.7510/jgjs.issn.1001-3806.2023.05.007
[8] 邵慧, 靳培龙, 王程, 等. 基于高光谱指数分割的草莓硬度预测与研究[J]. 激光技术, 2024, 48(3): 365-372. DOI: 10.7510/jgjs.issn.1001-3806.2024.03.011 SHAO H, JIN P L, WANG Ch, et al. Prediction and study of strawberry hardness based on hyperspectral index segmentation[J]. Laser Technology, 2024, 48(3): 365-372(in Chinese). DOI: 10.7510/jgjs.issn.1001-3806.2024.03.011
[9] SEKI H, MA T, MURAKAMI H, et al. Visualization of sugar content distribution of white strawberry by near-infrared hyperspectral imaging[J]. Foods, 2023, 12(5): 931. DOI: 10.3390/foods12050931
[10] 梁栋, 刘娜, 张东彦, 等. 利用成像高光谱区分冬小麦白粉病与条锈病[J]. 红外与激光工程, 2017, 46(1): 50-58. https://www.cnki.com.cn/Article/CJFDTOTAL-HWYJ201701008.htm LIANG D, LIU N, ZHANG D Y, et al. Distinguishing winter wheat powdery mildew from stripe rust using imaging hyperspectral[J]. Infrared and Laser Engineering, 2017, 46(1): 50-58(in Chinese). https://www.cnki.com.cn/Article/CJFDTOTAL-HWYJ201701008.htm
[11] XUE H, YANG Y, XU X P, et al. Application of near infrared hyperspectral imaging technology in purity detection of hybrid maize[J]. Applied Sciences, 2023, 13(6): 3507. DOI: 10.3390/app13063507
[12] 朱潘雨, 黄敏, 赵鑫. 基于SMOTE-UVE-SVM的小麦种子纯度高光谱图像检测[J]. 激光技术, 2024, 48(2): 281-287. DOI: 10.7510/jgjs.issn.1001-3806.2024.02.021 ZHU P Y, HUANG M, ZHAO X. Hyperspectral image detection of wheat seed purity based on SMOTE-UVE-SVM[J]. Laser technology, 2024, 48(2): 281-287(in Chinese). DOI: 10.7510/jgjs.issn.1001-3806.2024.02.021
[13] 陈彬, 刘阁, 张贤明. 连续投影算法的润滑油中含水量的近红外光谱分析[J]. 红外与激光工程, 2013, 42(12): 3168-3174. https://www.cnki.com.cn/Article/CJFDTOTAL-HWYJ201312002.htm CHEN B, LIU G, ZHANG X M. NIR spectral analysis of water content in lubricating oil by continuous projection algorithm[J]. Infrared and Laser Engineering, 2013, 42(12): 3168-3174(in Ch-inese). https://www.cnki.com.cn/Article/CJFDTOTAL-HWYJ201312002.htm
[14] 白丽敏, 李粉玲, 常庆瑞, 等. 结合SPA和PLS法提高冬小麦冠层全氮高光谱估算的精确度[J]. 植物营养与肥料学报, 2018, 24(5): 1178-1184. https://www.cnki.com.cn/Article/CJFDTOTAL-ZWYF201805006.htm BAI L M, LI F L, CHANG Q R, et al. Increasing accuracy of hyperspectral remote sensing for total nitrogen of winter wheat canopy by use of SPA and PLS methods[J]. Journal of Plant Nutrition and Fertilizers, 2018, 24(5): 1178-1184 (in Chinese). https://www.cnki.com.cn/Article/CJFDTOTAL-ZWYF201805006.htm
[15] 葛艳, 刘杏杏, 谢俊标. 改进粒子群优化BP神经网络的轮毂分类算法[J]. 中国科技论文, 2019, 14(7): 773-777. https://www.cnki.com.cn/Article/CJFDTOTAL-ZKZX201907012.htm GE Y, LIU X X, XIE J B. Improved particle swarm optimization BP neural network for wheel classification algorithm[J]. China Science and Technology Paper. 2019, 14(7): 773-777 (in Chinese). https://www.cnki.com.cn/Article/CJFDTOTAL-ZKZX201907012.htm
[16] 孟庆龙, 张艳, 尚静. 高光谱成像结合BP网络无损检测李子的硬度[J]. 激光与红外, 2019, 49(8): 968-973. https://www.cnki.com.cn/Article/CJFDTOTAL-JGHW201908011.htm MENG Q L, ZHANG Y, SHANG J. Hyperspectral imaging combined with BP neural network for non-destructive detection of plum hardness[J]. Laser and Infrared, 2019, 49(8): 968-973(in Ch-inese). https://www.cnki.com.cn/Article/CJFDTOTAL-JGHW201908011.htm
[17] 许兆美, 刘永志, 杨刚, 等. 粒子群优化BP神经网络的激光铣削质量预测模型[J]. 红外与激光工程, 2013, 42(9): 2370-2374. https://www.cnki.com.cn/Article/CJFDTOTAL-HWYJ201309013.htm XU Zh M, LIU Y Zh, YANG G, et al. Laser milling quality prediction model with particle swarm optimized BP neural network[J]. Infrared and Laser Engineering, 2013, 42(9): 2370-2374(in Chinese). https://www.cnki.com.cn/Article/CJFDTOTAL-HWYJ201309013.htm
[18] 李东年, 周以齐. 采用改进粒子群优化粒子滤波的三维人手跟踪[J]. 光学精密工程, 2014, 22(10): 2870-2878. https://www.cnki.com.cn/Article/CJFDTOTAL-GXJM201410039.htm LI D N, ZHOU Y Q. 3-D manual tracking using improved particle swarm optimized particle filtering[J]. Optics and Precision Engineering, 2014, 22(10): 2870-2878(in Chinese). https://www.cnki.com.cn/Article/CJFDTOTAL-GXJM201410039.htm
[19] 刘青, 李阳. 基于多尺度径向基函数的时变系统辨识[J]. 北京航空航天大学学报, 2015, 41(9): 1722-1728. https://www.cnki.com.cn/Article/CJFDTOTAL-BJHK201509023.htm LIU Q, LI Y. Identification of time-varying systems using multi-scale radial basis function[J]. Journal of Beijing University of Aeronautics and Astronautics, 2015, 41(9): 1722-1728(in Chinese). https://www.cnki.com.cn/Article/CJFDTOTAL-BJHK201509023.htm
[20] 朱晓琳, 李光辉, 张萌. 基于CARS-MIV-SVR的库尔勒香梨可溶性固体含量预测方法[J]. 光谱学与光谱分析, 2019, 39(11): 3547-3552. https://www.cnki.com.cn/Article/CJFDTOTAL-GUAN201911045.htm ZHU X L, LI G H, ZHANG M. Prediction method of soluble solid content of Korla Pear based on CARS-MIV-SVR[J]. Spectroscopy and Spectral Analysis, 2019, 39(11): 3547-3552(in Chinese). https://www.cnki.com.cn/Article/CJFDTOTAL-GUAN201911045.htm
-
期刊类型引用(0)
其他类型引用(1)
计量
- 文章访问数: 12
- HTML全文浏览量: 1
- PDF下载量: 7
- 被引次数: 1