Design and experiment of seed tube of seedling planter for Panax notoginseng
-
摘要:目的
为提高三七Panax notoginseng育苗机械化播种精确性，针对三七小行株距的特点，研究一种新型导种管。
方法分别建立有、无导种管情况下种子运动的力学模型。在土槽上开展了以投种高度、播种机前进速度、投种角度为试验因素，以理论播种点与实际播种点距离偏差的标准差和弹跳率为试验指标的三因素五水平二次回归正交旋转组合试验；借助Design-Expert 10.0.3软件建立了试验因素与指标间的回归方程及多目标优化数学模型，获取最佳投种参数组合；在最优参数组合下，使用高速摄像技术获取种子运动轨迹，结合运用图像处理技术和曲线拟合法进行种子运动轨迹曲线方程研究；借助EDEM软件进行导种管截面尺寸仿真分析，确定截面尺寸；用3D打印技术试制导种管并开展播种验证试验。
结果在有、无导种管情况下，影响种子离开排种器后运动的共同因素为排种器的周围速度、投种高度和投种角度。通过正交旋转组合试验得出最优投种参数组合为投种高度20 cm、播种机前进速度7.8 m/min、投种角度42°。在未安装导种管条件下，理论播种点与实际播种点距离偏差的标准差、弹跳率分别为51.66 mm和72.31%。研制的新型导种管出口和入口截面尺寸分别为26 mm × 30 mm和26 mm × 60 mm；安装新型导种管后验证试验结果显示，理论播种点与实际播种点距离偏差标准差为26.90，弹跳率为45.20%。
结论新型导种管满足了三七育苗机械化播种要求，播种精确性明显提高，研究结果可为三七育苗播种机导种管设计与田间播种应用提供了参考依据。
Abstract:ObjectiveIn order to improve the accuracy of mechanized seeding of Panax notoginseng, according to the characteristics of small row spacing of P. notoginseng, a new type of seed guide tube was studied.
MethodThe mechanical model was established with and without seed guide tube, respectively. The experiment was carried out in soil trough through the three-factor-five-level quadratic regression orthogonal rotation combination experiment with seed release height, seeder forward speed and seed release angle as the experimental factors, and the standard deviation of distance deviation between theoretical seeding point and actual seeding point, and bounce rate as the experimental indexes. With Design-Expert 10.0.3 software, the regression equation between test factors and indexes and the mathematical model of multi-objective optimization were established to get the best combination of seeding parameters. Under the optimal combination of parameters, high-speed camera technology was used to obtain the seed motion trajectory, and the curve equation of seed motion trajectory was studied by image processing technology and curve fitting method. With EDEM software, the cross-section size of seed guide tube was simulated and analyzed, and the cross-section size was determined. Seed guide tube was manufactured using 3D printing technology and seeding verification test was carried out.
ResultThe common factors affecting seed motion with and without the seed guide tube were circumferential speed of seed merering device, seed release height and seed release angle. Through orthogonal rotation combination test, the optimal combination of seeding parameters was as follows: The seed release height was 20 cm, the seeder forward speed was 7.8 m/min and the seed release angle was 42°. Without seed guide tube, the standard deviation of distance deviation between theoretical seeding point and actual seeding point was 51.66 mm, and the bounce rate was 72.31%. The exit and entrance cross-section sizes of the developed seed guide tube were 26 mm × 30 mm and 26 mm × 60 mm, respectively. The verification test results showed that after installing the new seed guide tube, the standard deviation of distance deviation between theoretical seeding point and actual seeding point was 26.90 mm and the bounce rate was 45.20%.
ConclusionThe new seed guide tube satisfies the requirement of the mechanized seeding of P. notoginseng, the seeding accuracy after installing the new seed guide tube is obviously improved. The research results provide a reference basis for the design and field seeding application of the seed guide tube of P. notoginseng seedling seeder.
-
Keywords:
- Panax notoginseng /
- Mechanized seeding /
- Seed guide tube /
- Motion trajectory /
- Image processing /
- EDEM simulation
-
野外环境下断蕾机器人需对多特征多目标进行识别，特别是农业机器人在切断香蕉花蕾和估产时，要同时识别香蕉果串、果轴和花蕾3种不同特征的目标，其单一目标的大小是一种变量，每次检测的目标尺寸和位置都不一样。另外，目标还受到树叶、遮挡及光照影响，研究复杂环境下的视觉检测方法意义重大。因此，针对香蕉、花蕾及果轴的识别问题，本研究提出多特征目标的快速识别方法。
断蕾机器人的关键视觉技术是对目标进行自动定位，视觉检测一直被认为是最具挑战性的问题之一[1-6]。Tang等[7]用多目视觉跟踪钢管的动态变形，探索其破坏机理；Dai等[8]把深度学习算法用于交通管理；研究还将视觉技术应用于水果采摘机器人识别[9-13]、估产[14]和病虫害检测[15]等。
在植物目标的视觉定位检测中，常用基于颜色、阈值分割、边缘检测等经典的图像处理算法[16]。罗陆锋等[17]结合颜色和聚类方法识别葡萄，精确度达87.63%。Hočevar等[18]基于HSL(色调、饱和度、亮度)识别估算苹果花。然而，颜色与阈值等识别方法受光照变化、枝叶或遮挡等因素影响。与经典模型和手工特征提取方法相比，深度卷积神经网络在识别方面取得了较好的性能[11]，在野外环境中检测目标具有优势。深度学习用于香蕉及果轴的多特征图像检测，可让模型自我学习目标特征，实现自动分类识别。
深度学习算法主要有2种类型的对象检测器，第1类为两阶段目标检测算法，它将检测分为2个阶段。以香蕉为例，第1阶段使含有香蕉的候选区域(Region proposal networks，RPN)生成1个包含待检感兴趣的候选框，选定后再作后续处理；第2阶段是把所感兴趣的香蕉区域提议发送到通道中，进行分类和边界框回归，用卷积神经网络(Convolutional neural network，CNN)完成对候选框中候选香蕉位置与类别的预测，例如R-CNN、Fast R-CNN、Faster R-CNN或Mask R-CNN。CNN是一种前馈神经网络，其卷积计算具有很深的结构，Mask R-CNN增加了图像分割功能，准确率较高，但速度较慢[19]。第2类为单阶段检测器，省去了第1阶段的RPN，直接用卷积神经网络模型预测目标边界，通过获取输入图像将对象检测视为回归问题，同时学习边界框坐标和类标签概率。与2阶段的探测器相比，该算法速度较快，但准确率稍低。常见的单阶段目标检测算法有SSD(Single shot multibox detector)、YOLO(You only look once)及其改进的系列算法YOLOv2、YOLOv3等[1, 20]。
对于第1类方法的研究，Alexe等[21]提出了一种线索综合性的方法来达到更好、更有效的预测效果。Lin等[22]提出一种具有重叠影响的草莓花识别系统，试验表明R-CNN、Fast R-CNN、Faster R-CNN的目标检出率分别约为63%、77%和86%。Tian等[23]提出改进Mask-R-CNN的图像分割算法，用于苹果花的多特征检测，准确率为96%。对于第2类方法的研究，Koirala等[24]用YOLO和R2CNN检测芒果花穗，认为YOLOv3模型用于计数时较优，R2CNN模型用于分类时较准确。陈燕等[25]改进了YOLOv3算法，用于检测荔枝串，平均检测精度为94%，结果可用于采摘机器人的预定位。
视觉检测及深度学习算法主要用于香蕉果实(果串)及树干(假径)、荔枝果串等单一目标检测[26-28]。但是机器人作业时需了解香蕉果实、花蕾及果轴情况，为切下花蕾提供定位信息。不同香蕉果树的香蕉果实、花蕾和果轴大小不等，且部分目标被遮挡，增加了检测复杂度，多特征目标的视觉精确识别方法仍然有待完善。
多种算法在检测精度和速度方面各有特点，目标图像分类识别的结果最终应用于视觉机器人作业控制。因此，研究者需要在目标检测的准确性和速度之间取得平衡[29]。机器视觉对香蕉多目标的快速精确识别成为关键问题。为了能够实现对果轴的精确切割，在图像采集设备将待采摘区域的图像数据输入视觉识别系统后，视觉识别系统应迅速查找画面中的果实、果轴和花蕾，使回归框的各边界准确地与目标边缘相切，并将回归框区域中的画面提取出来，以便进行进一步形态学处理和计算。对此，本文提出一种探索多特征多目标检测方法，优选模型参数，融合聚类算法优化YOLOv3模型的网络结构参数，并建立损失函数模型，形成YOLO-Banana模型。分别选用Faster R-CNN和YOLO-Banana进行多目标分类，对香蕉、花蕾及果轴进行试验验证算法的有效性，为复杂场景多目标多特征精确快速识别提供支持。
1. 试验设备和数据采集
1.1 试验设备
试验设备由硬件设备与软件组成。图像处理使用的计算机主要配置为：i7-7700K 处理器，内存为16 G、2400 MHz；显卡为GTX1080Ti、11G。样本采集时，用普通相机和高像素的手机，图像分辨率是1600万像素。
1.2 图像与数据采集
香蕉树通常只有香蕉果实串、花蕾及连接花蕾与香蕉的花序轴，简称果轴，花蕾及果轴朝下。为了提高产量，当香蕉果实串生长到一定时期，需要从果轴部位切下花蕾，称断蕾。本次采集样品包括香蕉果实、花蕾和果轴(图1)，以便为机器人断蕾提供定位信息。
![]()
图 1 顺、逆光下的香蕉果实、花蕾和果轴Figure 1. Banana fruit，flower bud and fruit axis target detection under frontlight and backlight野外机器人作业时无法通过增加标记或靶标等辅助手段来提高精度[3, 7]，而是采用无标记的目标定位。为了适应实际作业场景，视觉检测时不在目标上做定位靶标。在采样时，我们的策略是，只要有香蕉和花蕾就采样，以便机器人学习判别香蕉开花，做出决策。
香蕉多目标样本试验图像采集拍摄分2批进行。2020年7月25日，晴天，拍摄地点为广东省广州市水果世界岭南鲜果基地。2020年8月4、5日，晴天转阴天，拍摄地点为广东江门果园。采集图像时，包括香蕉果实、果轴、花蕾。图像采集过程中，将相机与花蕾间的距离保持在600 mm左右，在不同背景和光照条件下共采集1285幅图像。通过Labelimg标注软件对图片数据进行批量标注，通过编写脚本按比例自动随机选取训练样本集合、验证样本集合以及测试样本集，其中训练集1029张，占总样本80%；验证集64张，占总样本5%；测试集192张，占总样本15%。
2. 算法描述
2.1 图像的多目标分类
单个采样点会稍微变动拍摄角度来采集多幅图像，因此图片按照默认序列进行训练与检测时可能会出现测试集、验证集等均为少数采样点的若干相似图像的情况，使得训练、检测效果与样本本身有较强的关联性。为了避免这种关联性、减少随机性，在筛去无效图像以后，我们将剩余的图片随机打乱顺序，并按照1个标准格式进行命名。命名完成后，采用Labelimg软件对图像中的目标进行手工标注，并生成1个和图像同名的标签数据文档。图像数据与标签数据将按照VOC2007数据集格式进行分组与整理，按照6∶2∶2的比例划分为训练集、验证集和测试集。采集样图的像素为721×960，在进行图像处理需要把它缩放至像素为416×416。通过输入一张任意大小图片，缩放至像素为416×416，作为网络的输入，即网络的输入是1张像素为416×416、3通道的RGB图。对于比例不为1∶1的图片，通过在图片的上、下或者左、右分别添加2块区域，并将补上的区域以全黑(0，0，0)或全白(255，255，255)像素填充，将其补齐为一张1∶1的3通道图片，由于补齐的区域不包含任何与待检测目标相似的内容，所以对检测结果没有影响。补齐以后通过升或降采样的方式将图片统一为416×416×3的图片，即可输入网络。
2.2 多尺度特征融合
采用多尺度对香蕉多目标进行检测。分辨率对预测的影响主要取决于分辨率信息，也就是像素的数量(图2)。在第2类方法的YOLOv3中，我们采用(类似FPN)上采样(Upsample)和数据融合算法，融合了香蕉图像中的3个尺度(13×13、26×26和52×52)，在多尺度融合特征图上对目标进行独立检测，以优化效果。本文设计了一种新的分类器，通过设计多标签分类的Logistic分类器，用多个独立的Logistic分类器替代传统的Softmax层解决多标签分类问题，并保持准确性。传统Softmax分类器用于处理目标互斥的多分类问题，传统的Logistic分类器用于处理目标互斥的二分类问题。在采样得到的图片中，果轴、果实和花蕾并不存在明确的分界线，三者在图像上均有一定的重叠区域。若采用传统的Softmax分类器，由于各标签是互斥关系，所以这种区域重叠的情况会出现比较明显的分类错误或回归偏差；若采用传统的Logistic分类器，则无法处理这种多目标的检测问题。本文采用由多个独立的Logistic并联组成的分类器来解决这个问题，既能实现多目标检测，又能实现非互斥目标的准确分类。
![]()
图 2 基于香蕉多特征目标的YOLOv3网络结构Figure 2. YOLOv3 network structure based on banana multi-feature target2.3 基于香蕉多特征目标的YOLO v3网络结构
因算法采用多尺度融合特征的图像对多特征香蕉(含轴和花蕾)多目标进行单独检测，有效地增强了对大小不同目标和被遮挡目标的检测效果，并引入跃层连接强化收敛效果。因香蕉树生长于野外，已有的YOLOv3结构仍需要改进(图2)。本文在损失函数的部分进行改进，采用交叉熵损失函数代替方差损失函数，二者均为损失函数，损失函数可以衡量预测模型与真实模型的偏差，其偏差越小(即损失函数的值越小)，则预测模型与真实模型越接近。对于单样本的一次梯度下降过程中，有如下关系：
$$ \mathrm{线}\mathrm{性}\mathrm{部}\mathrm{分}：z=wx+b \text{，} $$ (1) $$ \mathrm{非}\mathrm{线}\mathrm{性}\mathrm{部}\mathrm{分}：\widehat{y}=a=\sigma \left(z\right) \text{，} $$ (2) 对于均方差损失函数，其损失函数定义为：
$$ {\rm{{Loss}_{MSE}}}(y,a)=\frac{1}{2}{\left(y-a\right)}^{2} \text{，} $$ (3) 对于交叉熵损失函数，其损失函数定义为：
$$ {\rm{{Loss}_{CE}}}(y,a)=-\left[y{\rm{lg}}\left(a\right)+\left(1-y\right){\rm{lg}}\left(1-a\right)\right] \text{，} $$ (4) 对于一次梯度下降中权重（w）和偏置（b）的更新，有
$$ w=w-\alpha \frac{{\partial} {\rm{Loss}}(y,a)}{{\partial} w} \text{，} $$ (5) $$ b=b-\alpha \frac{{\rm{\partial}} {\rm{Loss}}(y,a)}{{\rm{\partial}} b} \text{，} $$ (6) 分别将均方差损失函数与交叉熵损失函数对参数求偏导：
$$ 令： z=wx+b \text{，} $$ (7) 则对于均方差损失函数，有：
$$ \frac{\partial {{\rm{Loss}_{MSE}}}\left(y,a\right)}{\partial w}=-\left|y-\sigma \left(z\right)\right|{\sigma }^{\text{'}}\left(z\right)x \text{，} $$ (8) $$ \frac{\partial {{\rm{Loss}_{MSE}}}\left(y,a\right)}{\partial b}=-\left|y-\sigma \left(z\right)\right|{\sigma }^{\text{'}}\left(z\right) \text{，} $$ (9) 对于交叉熵损失函数，有：
$$ \frac{\partial {{\rm{Loss}_{CE}}}\left(y,a\right)}{\partial w}=x\left[\sigma \left(z\right)-y\right] \text{，} $$ (10) $$ \frac{\partial {{\rm{Loss}_{CE}}}\left(y,a\right)}{\partial b}=\sigma \left(z\right)-y \text{，} $$ (11) 可见，对于均方差损失函数，其对参数的偏导中均含有Sigmoid函数的导数，对于Sigmoid函数，其表达式为：
$$ \sigma \left(z\right)=\frac{1}{1+{{\rm{e}}}^{-z}} \text{，} $$ (12) 可见，当输入值很大或很小时，其导数趋近于0，这将导致权重和偏置变化速度急剧下降，即梯度消失。而交叉熵损失函数的偏导中不含有Sigmoid的导数项，可以有效避免梯度消失的问题，从而改善网络训练速度。
本文所采用的网络结构中该算法使用Darknet53作为骨架网络进行特征提取，再对特征图进行坐标回归和分类，结构中的最小组件由Conv+BN+Leaky_relu激活函数组成。
2.4 YOLO-Banana模型与评估
在网络模型参数设计与训练中，通常将epoch参数设置为100进行迭代，batch_size为16~32，花费时间和内存消耗相对增加。
在建立的YOLOv3网络中，优化了Iteration、batch_size和epoch这3个参数，使损失函数Loss曲线趋于稳定，再进一步地对YOLOv3模型的目标候选框维度进行聚类优化，以提高分类识别的精度。尽管YOLOv3自带默认的目标候选框数量、高宽比维度，它具有一定的通用性，但应用于野外复杂背景及变化光照环境下的香蕉多目标还需要优化改进。因此，本文融合聚类算法来优化YOLOv3模型，优化后的模型称为YOLO-Banana模型，利用k均值聚类算法，基于训练数据集，重新获得适应于预测的目标候选框个数、新的高与宽，并进行聚类分析，优化参数见表1。
表 1 香蕉多目标候选框参数Table 1. Parameters of banana multi target candidate box for banana候选框参数 Parameters of candidate box 1 2 3 4 5 6 7 8 9 新候选框高 Height of new candidate box 23 25 28 31 43 51 66 169 172 新候选框宽 Width of new candidate box 31 53 81 148 67 87 111 233 158 原候选框高 Height of original candidate box 10 16 33 30 62 59 116 156 373 原候选框宽 Width of original candidate box 13 30 23 61 45 119 90 198 326 为了评估模型的泛化能力并逐步优化模型，评估参数包括精度(Precision)、召回率(Recall)和二分类模型精确度的指标F1。计算公式如下：
$$ \mathrm{P}\mathrm{r}\mathrm{e}\mathrm{c}\mathrm{i}\mathrm{s}\mathrm{i}\mathrm{o}\mathrm{n}=\frac{\mathrm{T}\mathrm{P}}{\mathrm{T}\mathrm{P}+\mathrm{F}\mathrm{P}}\times 100{\text{%}} \text{，} $$ (13) $$ \mathrm{R}\mathrm{e}\mathrm{c}\mathrm{a}\mathrm{l}\mathrm{l}=\frac{\mathrm{T}\mathrm{P}}{\mathrm{T}\mathrm{P}+\mathrm{F}\mathrm{P}}\times 100{\text{%}} \text{，} $$ (14) $$ {F}_{1}=2\times \frac{\mathrm{P}\mathrm{r}\mathrm{e}\mathrm{c}\mathrm{i}\mathrm{s}\mathrm{i}\mathrm{o}\mathrm{n}\times \mathrm{R}\mathrm{e}\mathrm{c}\mathrm{a}\mathrm{l}\mathrm{l}}{\mathrm{P}\mathrm{r}\mathrm{e}\mathrm{c}\mathrm{i}\mathrm{s}\mathrm{i}\mathrm{o}\mathrm{n}+\mathrm{R}\mathrm{e}\mathrm{c}\mathrm{a}\mathrm{l}\mathrm{l}}\times 100{\text{%}} \text{，} $$ (15) 式中，TP表示正类判定为正类；FP表示负类判定为正类，“存伪”；FN表示正类判定为负类，“去真”。
训练中，batch size设计为8，即每次Iteration训练8张图片，故Iteration为128，对每个香蕉目标的全数据集进行1次运行，每次运行可得到1组Precision和Recall。当对模型算法设定不同的阈值时，会得到多组Precision和Recall，进而画出1条PR曲线，该曲线的面积为平均精度(Average precision, AP)。
3. 结果与分析
用上文提出的方法，在野外环境中完成香蕉多目标样本采集和训练模型。通过采样建立各类图像数据库，训练模型。
用视觉平台在香蕉果园采集香蕉、花蕾和果轴样本。软件系统主要以 OpenCV函数库与 YOLOv3为基础编写而成，硬件设备为有限通用水果采摘机器人及其视觉平台[28]，为末端执行器切断花蕾提供定位(图3)。
将本文数据集在YOLO-Banana模型的神经网络中进行训练和检测，选择训练的最优模型在验证集中进行验证，YOLO-Banana模型网络结构中的Neck部分采用的是FPN结构，通过多尺度反复提取主干层和检测层的特征，以提升网络检测微小目标效率。
用YOLO-Banana模型和Faster R-CNN模型进行比较。分别获得香蕉花蕾、香蕉果实和香蕉果轴的精度、召回率等评价参数。效果见图4、图5。
![]()
图 4 顺光下香蕉测试结果对比图a1, b1, c1表示香蕉原图；a2(1.00, 0.93, 0.99), b2(0.99, 0.92, 0.98), c2(1.00, 0.98, 1.00)表示YOLO-Banana模型处理后的香蕉图像；a3(1.00, 1.00, 1.00), b3(1.00, 1.00, 0.999), c3(1.00, 1.00, 0.999)表示 Faster R-CNN模型处理后的香蕉图像；括号内数字分别代表该算法对香蕉果实, 果轴和花蕾的置信度Figure 4. Comparison of banana test results under frontlighta1, b1 and c1 are original images of banana; a2(1.00, 0.93, 0.99), b2(0.99, 0.92, 0.98), c2(1.00, 0.98, 1.00) are banana images treated by YOLO-Banana model; a3(1.00, 1.00, 1.00), b3(1.00, 1.00, 0.999), c3(1.00, 1.00, 0.999) are banana images treated by Faster R-CNN model; The numbers in parentheses represent the confidence of the algorithm for banana fruit, fruit axis and flower bud, respectively![]()
图 5 逆光下香蕉测试结果对比a1表示香蕉原图；a2(0.97, 0.35, 0.99)表示YOLO-Banana模型处理后的香蕉图像；a3(1.00, 0.999, 0.999)表示 Faster R-CNN模型处理后的香蕉图像；括号内数字分别代表该算法对香蕉果实, 果轴和花蕾的置信度Figure 5. Comparison of banana test results under backlighta1 is original image of banana, a2 (0.97, 0.35, 0.99)is banana image treated by YOLO-Banana model; a3 (1.00, 0.999, 0.999) is banana image treated by Faster R-CNN model; The numbers in parentheses represent the confidence of the algorithm for banana fruit, fruit axis and flower bud, respectively从图4、图5可看出，香蕉和花蕾的检测精度较高，YOLO-Banana对较长的果轴检测精度较高，与Faster R-CNN检测结果接近；当果轴很短为小目标时，部分被向上长的花瓣遮挡，背景复杂且逆光，用YOLO-Banana检测目标的置信度很低，仅为0.35，Faster R-CNN则为0.999(图5)。在这种极端条件下目标的检出置信度会有所下降，但对回归位置没有显著影响。此外实际应用中这种极端条件出现频率极低，此处仅展示在极端恶劣检测环境下的检出效果。
YOLO-Banana与Faster R-CNN的香蕉多目标特征检测结果见表2。表2表明，Faster R-CNN对香蕉果实、花蕾、果轴的多特征目标的总平均识别精度为95.16%，平均每张图像处理的时间为0.434 s。YOLO-Banana对香蕉果实、花蕾、果轴的多特征目标的平均识别精度为91.03%,平均每张图像处理的时间为0.237 s。YOLO-Banana对香蕉多目标精度精度比Faster R-CNN对香蕉多目标精度要低4.13%，但是其识别速度要比Faster R-CNN的快。
表 2 YOLO-Banana与Faster R-CNN的香蕉多目标特征检测结果的对比Table 2. Comparison of multi-target feature detection results for banana between YOLO-Banana and Faster R-CNN模型
Model花蕾 Flower bud 果实 Fruit 果轴 Fruit axis 总平均精度/%
Mean of average
precisiont/s 平均精度/%
Average precision召回率/%
Recall平均精度/%
Average precision召回率/%
Recall平均精度/%
Average precision召回率/%
RecallFaster R-CNN 99 99.48 98 98.05 89 89.51 95.16 0.434 YOLO-Banana 97 97.38 95 95.12 81 85.19 91.03 0.237 为了测试野外光照的影响，选用了144个图片组，每个图片组内图片均为对同一个样本进行采样，包含香蕉、花蕾、果轴这3个部位的顺光与逆光图片进行识别试验(图4和图5)，分别用本文提出的YOLO-Banana模型及Faster R-CNN进行试验，结果见表3。表3表明，本文提出的YOLO-Banana模型在顺光时的平均识别精度与Faster R-CNN精度相近，误差仅为0.06%。
表 3 YOLO-Banana与Faster R-CNN的顺光逆光下的香蕉多目标检测结果的对比Table 3. Comparison of multi-target detection results for banana between YOLO-Banana and Faster R-CNN under frontlight and backlight conditions% 模型
Model花蕾精度
Flower bud average precision果实精度
Fruit average precision果轴精度
Fruit axis average precision总平均精度
Mean of average precision顺光
Frontlight逆光
Backlight顺光
Frontlight逆光
Backlight顺光
Frontlight逆光
Backlight顺光
Frontlight逆光
BacklightFaster R-CNN 100 100 100 100 93 92 97.64 97.47 YOLO-Banana 100 98 100 93 93 83 97.58 91.24 分析原因可知，因特征明显，且形状和色彩相似度高，花蕾识别精度最高；香蕉果实特征明显，但大小形状有一些差异，有时会被树叶遮挡，识别精度比花蕾低；果轴因其形状为近似小圆柱形，常被花瓣部分遮挡，识别率相对较小，特别是在逆光环境下，精度低。本文采用了多尺度检测算法，图像会以不同的分辨率进行检测识别，对不同大小的目标均有较好的检测效果。
断蕾主要考虑的技术指标是识别率，关键在于识别断蕾切割点所在的果轴，因此研究致力于提升果轴的识别率，比如通过聚类算法对各类别求出大小和比例不一的锚框进行网络训练。一般而言，果轴部分在整个生长周期中一直呈绿色，与未成熟的果实和香蕉叶片颜色相近，且四周均被香蕉果实包围，光照条件较差，边界纹理不清晰，且有些果轴短小易被遮挡，所以果轴的识别相对较难，效果不能达到果实和花蕾同时识别的要求。回归框位置的准确性会影响后续断蕾，针对特定目标框进行图像处理，不同类别进行的图像处理方式不一致(纹理色彩等不一致)，如对果轴而言，检测框越准确越能更好地进行后续边缘分割，从而获得果轴主方向，便于果轴切割。
由于深度学习无法确定断蕾位置，作者采用经典图像算法融合几何计算来确定断蕾点位置。在对多目标进行分类识别后，计算花蕾的边缘对称中心及其质心，建立花序轴求解方程，获取花序轴上的切断点。用边缘算法求得花蕾和花序轴图像边缘，利用花蕾的对称性求解花蕾质心与几何对称中心；然后，用几何方法计算花序轴上的切断点，该方法先求解检测到所有直线到质心之间的距离，再通过点线距离最小约束选取距离最短的直线作为切断点所在花序轴的位置，取线段中点坐标作为切断点。
在断蕾时，对切断果轴的轴向尺寸的精度属于自由精度范围，考虑到机器人在夹持切割机构及控制上具有一定的容错性，末端执行器对香蕉果轴检测角度具有一定的鲁棒性，当检测所得的香蕉果轴空间指向与实际香蕉果轴空间指向在一定锥角内，末端执行器均能自适应地修剪，并能获得较为一致的修剪效果。因此，野外断蕾机器人作业时，可以选择本文提出的YOLO-Banana算法。
4. 结论
为了提高野外作业机器人综合工作效率，考虑了检测精度和速度之间的平衡关系，研究了果园环境下用YOLOv3对香蕉多目标的分类识别方法，提出对香蕉目标进行多尺度特征提取及分类模型，对YOLOv3模型的一些参数进行了修改，并在此基础上，提出基于k均值聚类算法改进的YOLO-Banana模型，阐述了损失函数的误差原理，在置信度和分类误差模型中采用交叉熵的损失函数。研究了光照对香蕉果实、花蕾及果轴识别精度的影响。试验结果表明，用YOLO-Banana模型对香蕉果实、花蕾、果轴的检测精度分别为97%、95%和81%，多特征的目标总平均精度为91.03%，平均每张所需时间为0.237 s。Faster R-CNN对香蕉果实、花蕾、果轴的检测精度分别为99%、98%和89%，多特征的目标总平均精度为95.16%，平均每张所需时间为0.434 s。对比表明，前者速度快1.83倍，后者精度高。用YOLO-Banana模型与Faster R-CNN在顺光和逆光环境下进行了对比，两者顺光识别精度的误差很小，仅为0.06%，YOLO-Banana的顺光识别精度高于逆光识别精度，误差为8.26%。研究数据为机器人视觉识别与定位提供了参考。
-
![]()
图 1 三七育苗槽及播种示意图
1：槽肩，2：三七种子，3：土槽基质，4：育苗槽
Figure 1. Structure sketch of Panax notoginseng seedling trough and seeding
1: Slotted shoulder, 2: Panax notoginseng seeds, 3: Soil trough matrix, 4: Nursery trough
![]()
图 2 三七育苗播种机结构示意图
1：种箱，2：窝眼排种滚筒，3：导种装置，4：开沟器，5：镇压弹簧，6：限位拉伸手柄，7：限位轮，8：被动地轮，9：电动机，10：控制柜，11：传动装置，12：机架，13：升降螺杆，14：主动轮，15：万向轮
Figure 2. Structure sketch of Panax notoginseng seedling planter
1: Seed box, 2: Hole roller seed metering device, 3: Seed guide device, 4: Ditching device, 5: Crack spring, 6: Limit stretch handle, 7: Limit wheel, 8: Passive ground wheel, 9: Electric motor, 10: Control cabinet, 11: Transmission, 12: Frame, 13: Lifting screw, 14: Driving wheel, 15: Universal wheel
![]()
图 3 种子与导种管管壁接触过程的受力及运动分析图
建立一个以刚进入导种管种子中心点为坐标原点、以播种机前进速度的反方向为x0轴正方向、以竖直向下为y0轴正方向的坐标系O0；G为种子重力，N；N为导种管直线部分对种子的支持力，N；f为种子与导种管之间的摩擦力，N；β为导种管与水平线的夹角，(°)；H为种子离开排种器时离沟底的高度，m；V0为种子脱离排种器后初始运动速度，m/s；θ为投种角度，(°)；R1为排种器上型孔分布半径，m；n1为排种器转速，r/min；V2为育苗播种机的前进速度，m/s；V1为种子离开导种管时的初速度，m/s
Figure 3. Force and motion analysis graph of seeds which contact with seed guide tube wall
Establish a coordinate system O0 with the center of the seed just entering the seed guide tube as the coordinate origin, the reverse direction of the seeder’s forward speed as the positive direction of x0 axis, the vertical downward direction as positive direction of y0 axis; G is seed gravity, N; N is the support force of the straight part of the seed guide tube to the seed, N; f is the friction between seed and seed guide tube, N; β is the angle between the guide tube and the horizontal line, (°); H is the height of the seed from the seed bottom of the ditch when it leaves the seed metering device, m; V0 is the initial velocity of the seed after it is removed from the seed metering device, m/s; θ is the seed release angle, (°); R1 is the upper hole distribution radius on the seed metering device, m; n1 is the rotation speed of the seed metering device, r/min; V1 is the initial velocity of the seed when it leaves the seed guide tube, m/s; V2 is the forward speed of seedling planter, m/s
![]()
图 4 无导种装置时种子运动结构简图
建立一个以刚脱离排种器种子中心点为坐标原点、以播种机前进速度的反方向为x1轴正方向、以竖直向下为y1轴正方向的坐标系O1；H'为种子离开排种器时离沟底的高度，m；V'0为种子脱离排种器后初始运动速度，m/s；θ'为投种角度，(°)；R'1为排种器上型孔分布半径，m；n'1为排种器转速，r/min；V'2为育苗播种机的前进速度，m/s；V'1为种子离开导种管时的速度，m/s；V'1x为V'1的水平分速度；V'1y为V'1的竖直分速度
Figure 4. Structural sketch of seed motion without seed guiding device
Establish a coordinate system O1 with the center of the seed just separated from seed metering device as the coordinate origin, the reverse direction of the seeder’s forward speed as the positive direction of x1 axis; the vertical downward direction as positive direction of y1 axis; H' is the height of the seed from the bottom of the ditch when it leaves the seed metering device, m; V'0 is the initial movement speed of seeds after they are separated from the seed metering device, m/s; θ' is the seed release angle, (°); R'1 is the distribution radius of the upper hole on the seed metering device, m; n'1 is the rotation speed of the seed metering device, r/min; V'1 is the speed of seeds after leaving the seed guide tube, m/s; V'1x is the horizontal speed of V'1; V'1y is the vertical partial velocity of V'1; V'2 is the forward speed of seedling planter, m/s
![]()
图 5 播种精确性试验台
1：种箱，2：排种器，3：投种器，4：摄像机支架，5：机架，6：电控柜，7：传动链条，8：电动机，9：地轮
Figure 5. Seeding accuracy test bench
1: Seed box, 2: Seed metering device, 3: Seed pitching device, 4: Camera bracket, 5: Frame, 6: Control cabinet, 7: Transmission chain, 8: Electric motor, 9: Ground wheel
![]()
图 6 投种角度与播种精确性关系曲线
Figure 6. The curve of the relation between seed release angle and seeding accuracy
![]()
图 7 图像试验仪器安装图
1：高速摄像仪；2：手持云台；建立一个以排种器中心点为坐标原点、以水平正方向为x2轴正方向、以竖直向下为y2轴正方向的世界坐标系O2；$H^{\prime \prime} $为种子离开排种器时离沟底的高度，m；$V_0^{\prime \prime} $为种子脱离排种器后初始运动速度，m/s；$\theta ^{\prime \prime} $为投种角度，(°)；$R_1^{\prime \prime} $为排种器上型孔分布半径，m；$n_1^{\prime \prime} $为排种器转速，r/min
Figure 7. Installation picture of image testing instrument
1: High speed camera; 2: Handle stabilizer; Establish a world coordinate system O2 with the center of the seed metering device as the coordinate origin, the horizontal positive direction as x2 axis positive direction, the vertical downward direction as positive direction of y2 axis; $ H^{\prime \prime}$ is the height of the seed from the bottom of the ditch when it leaves the seed metering device, m; $V_0^{\prime \prime} $ is the initial movement speed of seeds after they are separated from the seed metering device, m/s; $\theta ^{\prime \prime} $ is the seed release angle, (°); R"1 is the distribution radius of the upper hole on the seed metering device, m; $n_1^{\prime \prime} $ is the rotation speed of the seed metering device, r/min
![]()
图 8 种子下落轨迹序列图像处理过程
Figure 8. Sequential image processing of seed falling trajectory sequence
![]()
图 12 拟合曲线型导种管
R2为导种管圆弧半径，mm；H2为拟合曲线设计高度，mm；H3为过渡放样曲线高度，mm
Figure 12. Fitting curve-shaped seed guide tube
R2 is the arc radius of seed guide tube, mm; H2 is the design height of fitted curve, mm; H3 is the height of transition lofting curve, mm
表 1 播种精确性试验因素水平表
Table 1 Seeding accuracy test factor level table
水平
level投种高度(Z1)/cm
Seed release height播种机前进速度(Z2)/(m·min−1)
Seeder forward
speed投种角度(Z3)/(°)
Seed release
angle1.682 40 10.0 45 1 36 9.2 41 0 30 8.0 35 −1 24 6.8 29 −1.682 20 6.0 25 Δj 6 1.2 6 
下载: 导出CSV
表 2 播种精确性试验方案与结果
Table 2 Seeding accuracy test scheme and result
序
号
No.投种高度 (Z1)
Seed release height播种机前进速度 (Z2)
Seeder forward
speed投种角度 (Z3)
Seed release
angle标准差(P1)/mm
Standard
deviation弹跳率(P2)/%
Bounce
rate1 −1 −1 −1 82.34 79.27 2 1 −1 −1 82.79 81.67 3 −1 1 −1 96.48 80.06 4 1 1 −1 101.52 82.53 5 −1 −1 1 51.08 67.98 6 1 −1 1 56.32 79.84 7 −1 1 1 56.80 75.95 8 1 1 1 79.59 73.25 9 −1.682 0 0 57.24 75.81 10 1.682 0 0 62.87 80.43 11 0 −1.682 0 60.04 78.79 12 0 1.682 0 104.95 83.70 13 0 0 −1.682 60.17 79.86 14 0 0 1.682 63.07 74.74 15 0 0 0 54.43 77.07 16 0 0 0 54.95 81.07 17 0 0 0 51.66 78.80 18 0 0 0 75.10 82.80 19 0 0 0 61.60 79.33 20 0 0 0 61.05 80.93 21 0 0 0 60.23 82.13 22 0 0 0 59.80 83.65 23 0 0 0 75.80 82.53
下载: 导出CSV
表 3 播种精确性试验方差分析1）
Table 3 Variance analysis of seeding accuracy test
来源
Source标准差
Standard deviation弹跳率
Bounce rate平方和
Sum of square均方
Mean squareF P 平方和
Sum of square均方
Mean squareF P 模型
Model3 914.67/
(3 600.15)434.96/
(1 200.05)3.27/
(11.14)0.026 3**/
(0.000 2**)224.35/
(161.81)24.93/
(53.94)3.89/
(7.03)0.013 5*/
(0.002 3**)Z1 135.32 135.32 1.02 0.331 8 34.80/
(34.80)34.80/
(34.80)5.44/
(4.54)0.036 5*/
(0.046 5*)Z2 1382.15/
(1 382.15)1 382.15/
(1 382.15)10.38/
(12.84)0.006 7**/
(0.002**)9.33 9.33 1.46 0.248 9 Z3 959.35/
(959.35)959.35/
(959.35)7.20/
(8.91)0.018 8*/
(0.007 6**)90.32/
(90.32)90.32/
(90.32)14.11/
(11.77)0.002 4**/
(0.002 8**)Z1Z2 61.27 61.27 0.46 0.509 5 26.25 26.25 4.10 0.063 9 Z1Z3 63.51 63.51 0.48 0.502 0 2.30 2.30 0.36 0.559 2 Z2Z3 1.88 1.88 0.014 0.907 2 0.01 0.01 0.00142 0.970 5 $ {Z}_{\text{1}}^{\text{2}} $ 15.36 15.36 0.12 0.739 5 24.40 24.40 3.81 0.072 8 $ {Z}_{\text{2}}^{\text{2}} $ 1 263.42/
(1 258.64)1 263.42/
(1 258.64)9.49/
(11.69)0.008 8**/
(0.002 9**)0.29 0.29 0.045 0.835 7 ${Z}_{\text{3}}^{\text{2}} $ 37.52 37.52 0.28 0.604 5 37.15/
(36.70)37.15/
(36.70)5.80/
(4.78)0.031 6*/
(0.041 4*)残差
Residual1 731.32/
(2 045.85)133.18/
(107.68)83.23/
(145.76)6.40/
(7.67)失拟项
Lack of fit1 147.58/
(1 462.11)229.52/
(132.92)3.15/
(1.82)0.072 9*/
(0.202 1)46.31/
(108.85)9.26/
(9.90)2.01/
(2.14)0.181 8/
(0.144 1)误差
Pure error583.74/
(583.74)72.97/
(72.97)36.91/
(36.91)4.61/
(4.61)总和
Total5645.99/
(5645.99)307.57/
(307.57)1)Z1：投种高度，Z2：播种机前进速度，Z3：投种角度；“*”表示在P≤0.05水平差异显著，“**”表示在P≤0.01水平差异显著；“/”下“( )”内数字为剔除不显著项后重新进行方差分析的结果
1) Z1: Seed release height, Z2: Seeder forward speed, Z3: Seed release angle;“*” indicates significant difference at P≤0.05 level, “**” indicates significant difference at P≤0.01 level; Data in “( )” under “/” are variance analysis result after rejecting unsignificant items
下载: 导出CSV
表 4 新型导种管播种验证试验结果
Table 4 Verification experiment results of seeding with new seed guide tube
试验编号
Number of test标准差/mm
Standard deviation弹跳率/%
Bounce rate1 23.10 41.80 2 23.90 45.30 3 33.70 48.50 均值 Mean 26.90 45.20
下载: 导出CSV
-
[1] 赖庆辉, 于庆旭, 苏微, 等. 三七超窄行气吸式精密排种器设计与试验[J]. 农业机械学报, 2019, 50(4): 102-112. doi: 10.6041/j.issn.1000-1298.2019.04.012 [2] 刘素, 赖庆辉, 董家宇, 等. 气吹式三七精密排种器充填性能的仿真与试验[J]. 华南农业大学学报, 2019, 40(3): 125-132. doi: 10.7671/j.issn.1001-411X.201806023 [3] 奚琪, 杨文彩, 杜迁, 等. 三七槽式育苗配套工程技术指标分析[J]. 广东农业科学, 2013, 40(16): 29-31. doi: 10.3969/j.issn.1004-874X.2013.16.010 [4] 杨文彩, 朱有勇, 杜迁, 等. 基于农机农艺融合的三七机械化精密播种系统研究[J]. 广东农业科学, 2014, 41(2): 175-180. doi: 10.3969/j.issn.1004-874X.2014.02.041 [5] 李宝筏. 农业机械学[M]. 北京: 中国农业大学出版社, 2003. [6] 雷小龙, 廖宜涛, 丛锦玲, 等. 油菜小麦兼用气送式直播机集排器参数优化与试验[J]. 农业工程学报, 2018, 34(12): 16-26. doi: 10.11975/j.issn.1002-6819.2018.12.003 [7] 李欢, 杨文彩, 阮解琼, 等. 三七精密排种器的参数优化设计[J]. 农机化研究, 2015, 37(7): 131-134. doi: 10.3969/j.issn.1003-188X.2015.07.030 [8] 杨航, 杨文彩, 李超群, 等. 2BQ—27型三七精密播种机排种性能试验研究[J]. 中国农机化学报, 2016, 37(2): 33-38. [9] 中华人民共和国国家质量监督检验检疫总局, 中国国家标准化管理委员会. 单粒(精密)播种机试验方法: GB/T 6973—2005[S]. 北京: 中国标准出版社, 2006. [10] 贾洪雷, 郑嘉鑫, 袁洪方, 等. 大豆播种机双V型筑沟器设计与试验[J]. 吉林大学学报(工学版), 2017, 47(1): 323-331. [11] 廖庆喜, 舒彩霞. 精播机输种管对排种均匀性影响的试验研究[J]. 湖北农业科学, 2004(4): 126-128. doi: 10.3969/j.issn.0439-8114.2004.04.045 [12] 王金武, 唐汉, 王金峰, 等. 指夹式玉米精量排种器导种投送运移机理分析与试验[J]. 农业机械学报, 2017, 48(1): 29-37. doi: 10.6041/j.issn.1000-1298.2017.01.005 [13] 马旭, 马成林, 张守勤, 等. 精密播种机接输种管后田间植株分布的研究[J]. 农业工程学报, 1997, 3(3): 119-123. doi: 10.3321/j.issn:1002-6819.1997.03.025 [14] YAZGI A. Effect of seed tubes on corn planter performance[J]. Applied Engineering in Agriculture, 2016, 32(6): 783-790. doi: 10.13031/aea.32.11274
[15] 刘立晶, 杨慧. 基于Geomagic Design软件的导种管三维逆向工程设计[J]. 农业工程学报, 2015, 31(11): 40-45. doi: 10.11975/j.issn.1002-6819.2015.11.006 [16] 杨丽, 颜丙新, 张东兴, 等. 玉米精密播种技术研究进展[J]. 农业机械学报, 2016, 47(11): 38-48. doi: 10.6041/j.issn.1000-1298.2016.11.006 [17] 丁中凯. 2BJZ-10型播种机设计及导种管结构参数研究[D]. 哈尔滨: 东北农业大学, 2012. [18] 贺俊林, 杨涛, 裘祖荣. 精密播种机导种管曲线的设计研究:2BQYF—6A硬茬精密播种机导种管的设计研究[J]. 山西农业大学学报, 2000(4): 389-391. [19] KOCHER M F, COLEMAN J M, SMITH J A, et al. Corn seed spacing uniformity as affected by seed tube condition[J]. Applied Engineering in Agriculture, 2011, 27(2): 177-183. doi: 10.13031/2013.36484
[20] GARNER E B, THIEMKE D B, RYLANDER D J, et al. Seeding machine with seed delivery system: US9807900B2[P]. 2017-10-07.
[21] 陈学庚, 钟陆明. 气吸式排种器带式导种装置的设计与试验[J]. 农业工程学报, 2012, 28(22): 8-15. doi: 10.3969/j.issn.1002-6819.2012.22.002 [22] 康建明, 温浩军, 王士国, 等. 带式导种装置对排种均匀性影响的试验研究[J]. 中国农机化学报, 2015, 36(5): 42-45. [23] 赵淑红, 陈君执, 王加一, 等. 精量播种机V型凹槽拨轮式导种部件设计与试验[J]. 农业机械学报, 2018, 49(6): 146-158. doi: 10.6041/j.issn.1000-1298.2018.06.017 [24] 雷小龙, 廖宜涛, 张闻宇, 等. 油麦兼用气送式集排器输种管道气固两相流仿真与试验[J]. 农业机械学报, 2017, 48(3): 57-68. doi: 10.6041/j.issn.1000-1298.2017.03.007 [25] YATSKUL A, LEMIERE J P, COINTAULT F. Influence of the divider head functioning conditions and geometry on the seed’s distribution accuracy of the air-seeder[J]. Biosystems Engineering, 2017, 161: 120-134. doi: 10.1016/j.biosystemseng.2017.06.015
[26] 戴亿政, 罗锡文, 王在满, 等. 气力集排式水稻分种器设计与试验[J]. 农业工程学报, 2016, 32(24): 36-42. doi: 10.11975/j.issn.1002-6819.2016.24.005 [27] 陶桂香, 衣淑娟, 毛欣, 等. 水稻植质钵盘精量播种装置投种过程的动力学分析[J]. 农业工程学报, 2013, 29(21): 33-39. doi: 10.3969/j.issn.1002-6819.2013.21.005 [28] 余佳佳, 丁幼春, 廖宜涛, 等. 基于高速摄像的气力式油菜精量排种器投种轨迹分析[J]. 华中农业大学学报, 2014, 33(3): 103-108. doi: 10.3969/j.issn.1000-2421.2014.03.018 [29] 刘宏新, 徐晓萌, 刘俊孝, 等. 利用高速摄像及仿真分析立式浅盆型排种器工作特性[J]. 农业工程学报, 2016, 32(2): 13-19. doi: 10.11975/j.issn.1002-6819.2016.02.003 [30] 廖宜涛, 李成良, 廖庆喜, 等. 播种机导种技术与装置研究进展分析[J]. 农业机械学报, 2020, 51(12): 1-14. doi: 10.6041/j.issn.1000-1298.2020.12.001 [31] 郭慧, 王刚, 赵佳乐, 等. 种子纵向分布均匀性指标及空间分布均匀性评价方法[J]. 吉林大学学报(工学版), 2020, 50(3): 1120-1130. [32] 杨航. 2BQ-27型三七精密播种机关键部件设计与试验研究[D]. 昆明: 云南农业大学, 2015. [33] 辛尚龙, 赵武云, 戴飞, 等. 旱区全膜双垄沟播履带式玉米联合收获机的设计[J]. 农业工程学报, 2019, 35(14): 1-11. doi: 10.11975/j.issn.1002-6819.2019.14.001 [34] 沈明霞, 李嘉位, 陆明洲, 等. 基于动态多特征变量的黄羽肉鸡跛行状态定量评价方法[J]. 农业机械学报, 2018, 49(9): 35-44. doi: 10.6041/j.issn.1000-1298.2018.09.004 [35] 屈晶晶, 辛云宏. 连续帧间差分与背景差分相融合的运动目标检测方法[J]. 光子学报, 2014, 43(7): 219-226. [36] 张萌, 钟南, 刘莹莹. 基于生猪外形特征图像的瘦肉率估测方法[J]. 农业工程学报, 2017, 33(12): 308-314. doi: 10.11975/j.issn.1002-6819.2017.12.040 [37] PETROU M, PETROU C. 图像处理基础: 第2版[M]. 章毓晋, 译. 北京: 清华大学出版社, 2013. [38] 李超群. 2BQ-28型三七精密播种机关键部件的研究与设计[D]. 昆明: 云南农业大学, 2016. [39] 于庆旭, 刘燕, 陈小兵, 等. 基于离散元的三七种子仿真参数标定与试验[J]. 农业机械报, 2020, 51(2): 123-132. [40] 严荣俊. 精密播种机导种管结构设计与试验研究[D]. 北京: 中国农业大学, 2012. [41] 杨文彩, 徐路路, 杜一帆, 等. 三七育苗播种压轮仿形开沟装置的设计与试验[J]. 农业工程学报, 2020, 36(7): 53-62. doi: 10.11975/j.issn.1002-6819.2020.07.006 -
期刊类型引用(1)
1. 杨帆，李优佳，SHAMEER K S，阎伟，吕宝乾，蒋方一丁，章雨璐，涂艳，齐可欣. 原产地与入侵地椰子织蛾遗传分化特征. 福建农林大学学报(自然科学版). 2021(02): 178-183 . 
百度学术
其他类型引用(0)
计量
- 文章访问数: 500
- HTML全文浏览量: 18
- PDF下载量: 1067
- 被引次数: 1
