Effects of straw mulching on nutrient release and nutrient balance in potato soil
-
摘要:目的
探索马铃薯种植过程中稻草覆盖还田模式下稻草N、P、K释放和土壤N、P、K养分平衡的特征。
方法采用盆栽试验的方法，对稻草覆盖还田后的腐解速率及养分释放状况进行研究。
结果稻草覆盖还田82 d后，其C、P含量变化较小，K含量降低93.96%，而N含量却提高了49.22%；覆盖稻草可为土壤提供大量养分且以K居多，每覆盖还田100 kg稻草可为马铃薯生长提供N 0.42 kg、P 0.14 kg(P2O5)和K 1.70 kg(K2O)；稻草腐解率为70.87%；各养分累积释放率表现为：K>C>P>N，C、N和P累积释放率分别为71.89%、56.85%和68.75%，K累积释放率最高，为98.22%。覆盖稻草和施用钾肥对土壤K平衡影响较大，不施钾肥和稻草土壤的K亏缺达163.80 kg/hm2，覆盖稻草或施用钾肥均可有效降低土壤钾素亏缺量，稻草与钾肥配合施用可使土壤K收支平衡，甚至盈余。
结论在稻草覆盖还田情况下，可适当减少钾肥用量，以降低生产成本，缓解我国钾肥资源匮乏的困境。
Abstract:ObjectiveTo study the characteristics of nitrogen(N), phosphorus(P) and potassium (K) nutrient release and soil N, P, K nutrient balance under straw mulching pattern during potato planting process.
MethodWe studied the decomposition rate and nutrient release status of mulching straw using pot experiment.
ResultStraw carbon (C) and P contents changed slightly after mulching for 82 days, while K content decreased by 93.96%, and N content increased by 49.22%. The mulching straw provided a large amount of nutrients (especially K) with 100 kg mulching straw contributing 0.42 kg N, 0.14 kg P (P2O5), and 1.70 kg K (K2O) to soil used for potato growth. Straw decomposition rate was 70.87%. The nutrient accumulative release rates were performed as: K (98.22%) > C (71.89%) > P (68.75%) > N (56.85%). Straw mulching and K fertilization largely affected soil K balance. The K deficient amount was 163.80 kg/hm 2 without straw mulching and K fertilizer application. K deficient amount was effectively reduced by straw mulching and K fertilization, and even K had surplus with combined application of straw and K fertilizer.
ConclusionThe straw mulching pattern can reduce K fertilizer application amount and production cost, and relieve the dilemma of K fertilizer resource deficiency in China.
-
Keywords:
- potato /
- nutrient release /
- nutrient balance /
- straw mulching /
- decomposition
-
奶牛的采食行为是经常被用于预测奶牛疾病的行为之一[1-2]，奶牛在采食过程中的头部运动变化可以用于评价奶牛的采食状况，进而评价饲草料的优劣。因此在现代化智能养殖场中，监测奶牛采食行为并对奶牛头部的活动进行跟踪记录是非常重要的。研究发现，采食量和产奶量的下降与疾病有关[3]，患有乳腺炎的奶牛的饲料和干物质采食量会下降。跛足牛和非跛足牛在采食时间、采食频率和采食量方面存在显著差异，根据这个差异可以判断奶牛的腿部健康状况[4]。这些研究已经清楚地表明疾病会对采食行为有重大影响。因此，对奶牛采食行为进行监测是保障智能畜牧和提高动物福利的重要手段。
在监测奶牛采食行为的技术中，最基本的问题是准确识别奶牛的采食行为。有研究表明，实时定位系统和超宽带系统能够通过估计奶牛在采食区花费的时间来判断奶牛的采食行为[5]；但不能确定奶牛是否进行有效进食。还有研究将压力传感器安装在奶牛身上监测咀嚼运动，虽然传感器具有良好的监测性能，但使用相对繁琐，且给奶牛佩戴时不方便[6]。利用基于惯性测量单元和射频识别系统的耳挂式传感器可以对奶牛进行定位以及采食行为监测[7-8]，但只能判断奶牛的动作幅度，不能判断奶牛是否有效进食。从以上研究可以看出，在奶牛采食行为监测方面已经取得一些研究成果。但是，上述研究通常成本高或者对采食行为判断不准确，无法应用到实际养殖场饲养环境中。因此，需要进一步寻找有效的解决方法。
近些年来，深度学习在畜牧业中得到了广泛应用[9]。白强等[10]使用改进YOLOv5s网络实现奶牛的站立、喝水、行走、躺卧等行为识别。杨阿庆[11]利用全卷积网络对母猪头部特征进行提取，并将头部与采食区重叠区域作为空间特征来识别母猪采食行为。王政等[12]提出了融合YOLOv5n与通道剪枝算法的轻量化识别方法，实现了奶牛发情行为识别。Bezen等[13]基于CNN模型和RGB-D摄像机设计出一种用于奶牛个体采食量测量的系统。Lao等[14]开发了一种基于预处理深度图像数据的算法，用于研究母猪的摄食行为。Shelley等[15]设计了一种3D图像分析方法，用于测量奶牛饲喂期前后饲喂器中可用食物量的变化。以上研究表明，深度学习算法用于奶牛的采食行为研究是可行的，但是目前尚未有研究对奶牛采食动作进行区分以及对奶牛采食过程进行跟踪，因此对奶牛采食行为识别以及跟踪的研究是必要的。
由于奶牛在采食过程中动作快、采食动作特征复杂，很容易对目标产生错检或者漏检，所以本文从以上问题出发，以YOLOv8模型为基础，通过使用BiFormer模块的多头自适应注意力机制加强对奶牛采食行为特征的学习、提取，使用CoT、EMA模块对奶牛特征进行强化并加强分类，以此来增强对采食行为识别检测的准确性，提出一种TCE-YOLO深度学习算法模型，并与Deep SORT算法结合以实现对奶牛采食时头部的跟踪与轨迹绘制。
1. 材料与方法
1.1 数据采集与预处理
本研究的奶牛采食行为数据集来自山东省泰安市泰安金兰奶牛养殖公司。试验数据利用ZED2双目深度相机(生产商为STEREOLABS公司)进行拍摄，本文使用单目分辨率为
1280 像素×720像素，传输帧率能达到稳定30帧/s。由于牛栏顶端距地面高度为1.35 m、饲料带宽度0.8 m，为了避免干涉奶牛的正常采食行为，本文选择相机拍摄位置为奶牛采食区上方1.75 m处以及奶牛采食区前方1.2 m、高0.8 m处，采集不同拍摄方向下奶牛的采食行为。在图像采集过程中，一共采集视频数据20组，每组包含4个视频，每个视频长度为10~14 min。为了保证动作的连贯性，本文利用ZED API对视频以一秒一帧的方式提取图像，去除重复、模糊、重影以及不包含奶牛采食行为的无效图像，选取有效数据集11 288张图像，17头奶牛的采食动作均包含在数据集中，其中前方拍摄图像数量为5 742张，上方拍摄图像数量为5 546张。
本文使用开源标注工具LabelImg人工对
11288 张奶牛采食行为原始图像进行标注，标签TTJJ为咀嚼行为，DTCS为采食行为，GS为拱草行为。根据动作时间提取连续帧，将这一段连续帧定义一个动作，通过观察和统计得出各行为时长：采食行为持续时长为15~35 s，咀嚼行为持续时长为5~10 s，拱草行为持续时长为1~15 s。将原始数据集分为上方拍摄数据集和前方拍摄数据集，每组数据集选取500张作为验证集，在剩余图像中，将上方拍摄数据集按照6∶1的比例分为4320 张训练数据集图像和726张测试数据集图像，将前方拍摄的数据集按照6∶1的比例分为4484 张训练数据集图像和758张测试数据集图像。在养殖场饲养环境下，奶牛只有在采食区且奶牛头部与饲料接触时才可能进行采食行为，而在采食区抬头时基本处于咀嚼饲料以及准备继续采食的过程，当奶牛摄取足够的饲料后便会离开采食区。因此，在分析奶牛采食行为过程中，将奶牛采食行为分为采食和咀嚼两大部分，其中采食还可分为采食和拱草2种行为。其中，各个行为在数据集中占比如表1所示。
表 1 训练数据集中标注的奶牛不同行为数量Table 1. Number of cows’ various behaviors labeled in the training dataset拍摄方向
Shooting
direction数据集数量
Number of
datasets不同行为数量
Number of various behaviors采食
Feeding咀嚼
Chewing拱草
Arching前方 Ahead 4 484 5 684 792 1 613 上方 Above 4 320 6 958 960 1 946 1.2 改进BCE-YOLO模型
1.2.1 BCE-YOLO模型建立
YOLOv8是目前广泛应用的模型，使用了更为丰富的数据增强技术，更好地增加了训练数据的多样性，有效提高了模型的泛化能力，还具有更好的自适应图像缩放功能，可以确保目标跟踪时不会因为图像长宽比改变而影响检测位置像素判断。YOLOv8模型的检测速度虽然很快，但是存在检测精度不高、对模型的特征提取不全面的问题，因此本文以YOLOv8为基础，使用模型增强模块来提升模型的检测性能。
由于目标跟踪的前提是需要对模型特征进行准确的提取，因此在高效率YOLOv8的基础上增加BiFormer模块，利用BiFormer模块的自适应注意力机制实现对奶牛采食行为特征的准确掌握与判断；在Focus层后增加CoT模块，使得模型的注意力集中于奶牛采食行为的特征，在骨干网络末端融合EMA模块，提升特征图的全局感受野以及模型的分类能力。在模型预测模块方面，引用了CIoU[16]代替GIoU，考虑到尺度、距离、惩罚项和锚框之间的重叠率对损失函数的影响，目标帧回归变得更加稳定，公式如下：
$$ {\mathrm{CIoU}} = {\mathrm{IoU}} - \dfrac{{{\rho ^2}(b,{b^{g_{\mathrm{t}}}})}}{{{c^2}}} - \alpha v \text{，} $$ (1) $$ \alpha = \dfrac{v}{{1 - {\mathrm{IoU}} + v}} \text{，} $$ (2) $$ v = \dfrac{4}{{{{\text{π}} ^2}}}{\Bigg(\arctan \dfrac{{{W^{g_{\mathrm{t}}}}}}{{{H^{g_{\mathrm{t}}}}}} - \arctan \dfrac{W}{H}\Bigg)^2} \text{，} $$ (3) $$ {\mathrm{LOS}}{{\mathrm{S}}_{{\mathrm{CIoU}}}} = 1 - {\mathrm{IoU}} + \dfrac{{{\rho ^2}(b,{b^{g_{\mathrm{t}}}})}}{{{c^2}}} + \alpha v \text{，} $$ (4) 式中，c表示同时包含预测边界框和对应地面真实边界框的最小封闭区域的对角线距离，ρ2(b,
$b^{g_{\mathrm{t}}}$ )表示预测中心点与真实帧之间的欧式距离，b、$b^{g_{\mathrm{t}}} $ 分别表示预测框和真实的中心点，W、H分别表示边界框的宽和长，gt表示真实值，α表示权重参数，v表示度量长宽比的相似性，LOSSCIoU表示考虑重叠面积、中心点距离、纵横比的损失，最后用1−IoU得到相应的损失。BCE-YOLO模型工作流程如图1所示，输入图像进入模型时，首先通过CBS模块进行切片操作，将W、H信息集中到通道空间，之后通过CoT模块，进入C2f网络进行一系列卷积操作，得到高维度特征图后通过第8层CBS结构进入BiFormer模块，利用多头自适应注意力结构对特征进行并行学习，将结果输入至EMA模块来增强模型的全局感受野，然后进入FPN特征金字塔网络[17]和PAN路径[18]增强网络，通过上采样以及张量拼接等方式与骨干网络建立联系，丰富特征图的特征信息，最终将特征图送至YOLO检测模块，得到检测结果。
1.2.2 BiFormer模块
BiFiormer模块[19]是一种新型的视觉网络结构，结构如图2所示，在第1阶段使用重叠块嵌入，在第2~4阶段使用合并模块来降低输入空间分辨率，同时增加通道数(C)，然后采用连续的BiFormer模块做特征变换。该模块通过多头自适应注意力层，扩展了模型专注于不同位置的能力，形成多个子空间，让模型去关注不同方面的信息，可以有效结合奶牛身份特征和奶牛采食时的不同行为，实现对高维全局特征的整合。
![]()
图 2 BiFormer结构及模块细节N1、N2、N3、N4分别表示在阶段1、2、3、4的迭代次数Figure 2. Structure and module details of BiFormerN1, N2, N3 and N4 indicate iterations at stage 1, 2, 3 and 4, respectively1.2.3 CoT模块
CoT模块[20]包括了局部注意力原则，并且增强了信息提取。结构如图3所示，通过将静态特征和动态特征结合、局部与全局结合，节约了参数和计算力，保证了模块的轻量化，可以直接集成到现有的网络模型构架中。引入CoT模块后，对于奶牛采食时的头部行为，网络识别特征覆盖到的奶牛采食时头部特征更多，代表奶牛采食时头部运动的特征已经成为模型关注的重点信息。
1.2.4 EMA模块
EMA模块[21]是一种高效多尺度注意力模块，结构如图4所示，EMA模块通过1×1分支进行图像聚合，实现不同跨通道交互特征。跨纬度学习模块进行激励操作，通过全连接层和激活函数，减少通道从而降低计算量，得到第1张空间注意力图。通过3×3分支，经过激活函数与1×1分支的输出结果点积，获得第2张保留精准信息的空间注意力图，最后再与第1张注意力图聚合。
通过EMA模块后，可以增强通道注意力，使特征图具有全局感受野，增强网络的特征图提取且提升网络对目标的分类情况。
1.2.5 Deep SORT算法
多目标跟踪对采食行为识别能够连续、动态地追踪场景中多只奶牛的行为，从而更准确地捕捉和分析它们的采食模式和频率，这有助于深入研究多头奶牛的采食行为特征，提升对采食行为的识别精度。
Deep SORT算法[22]常用于多目标跟踪，工作流程如图5所示，通过BCE-YOLO将奶牛采食时的个体身份以及对应的采食动作计算出来，然后通过卡尔曼滤波器预测物体下一帧的位置，将预测位置和实际检测位置作特征对比并计算IoU，得到相邻2个目标的相似度，最终通过匈牙利匹配算法得到相邻帧的对应ID (Identity document)，实现对奶牛身份以及采食动作的跟踪。
由于Deep SORT算法使用的是奶牛前后特征对比确认，并且增加了级联匹配和新轨迹确认。如图6所示，目标检测网络得到奶牛的唯一特征后，通过预测、观测、更新等方式对奶牛特征信息进行加强，因此使得Deep SORT算法所检测的奶牛ID不容易改变，从而更加精确地跟踪、检测奶牛的采食行为。
2. 结果与分析
2.1 评价指标
为了准确评价模型的性能，本文采用了目标检测算法常用的4个性能评价指标：精确度(Precision)、召回率(Recall)、平均精确度均值mAP(mean Average precision)和F1分数(F1-score)来检验本文提出的模型的性能。损失曲线拟合速度越快，拟合程度越好，最终损失值越低，代表模型的性能越强。
2.2 试验配置和超参数设置
训练以及测试的硬件平台配置为Intel(R) Xeon(R) W-2135 CPU@3.70GHz处理器，32 GB内存，11 GB GeForce RTX2080 Ti显卡，500 G固态硬盘。平台操作系统为64位Ubuntu18.04操作系统，利用Ubuntu18.04终端实现，程序运行中调用加速环境CUDA、Cudnn和第三方库OpenCV、Pytorch，CUDA版本为10.2.89，Cudnn版本基于CUDA10.2下cuDNN 7.6.5，OpenCV版本为4.4.0，Pytorch版本为1.7.1。在不影响识别精度下，为了增加训练效率，选择输入图像为416像素×416像素。设置初始学习率为0.001，权重衰减率为0.000 5，迭代次数为100次。
2.3 目标跟踪算法检测效果
为了验证奶牛采食行为目标跟踪算法的可靠性，本文分别使用上方和前方拍摄的奶牛采食行为数据集来对算法进行测试。
如图7所示为奶牛采食行为前方数据集的目标跟踪效果，由图7A、7B可以看出，当1号奶牛的采食行为从采食变为咀嚼时，奶牛的ID并不会发生改变，仍然对奶牛的头部运动轨迹进行跟踪绘制；由图7C、7D可以看出，当1号奶牛遮挡住2号奶牛时，2号奶牛的标签消失，但当2号奶牛恢复到检测状态时，2号奶牛的ID没有发生改变，说明目标检测网络对奶牛特征的提取充足，使得在目标跟踪时根据目标的特征进行跟踪不会使ID发生跳动。
![]()
图 7 前方拍摄奶牛采食数据集的目标跟踪Figure 7. Target tracking of dairy cow feed dataset taken in front图8为奶牛采食行为上方数据集的目标跟踪效果，与前方相同，奶牛采食过程中行为变化以及遮挡并不会使得奶牛的ID发生变化。
为了验证所提出的Deep SORT算法与BCE-YOLO模型结合的有效性，基于本文自建数据集与其他主流算法进行对比试验，包括Staple和SiameseRPN，结果如表2所示。Deep SORT的跟踪精度和鲁棒性都占有优势，说明Deep SORT算法与BCE-YOLO模型结合具有更好的跟踪性能。通过以上分析可以看出，BCE-YOLO和Deep SORT结合的目标跟踪算法能精确追踪奶牛采食时的头部运动并绘制轨迹，通过运动总量可以进一步判断采食行为，也为后续对奶牛采食行为其他研究做准备。将奶牛运动轨迹点的数据导出，可以计算奶牛整个采食过程所持续的时间以及奶牛采食过程运动的总量。如图9所示，图9A以奶牛采食行为跟踪视频的左上顶点为原点，通过BBOX计算头部中心点并绘制轨迹。由于逐帧跟踪且视频帧率为30帧/s，可计算采食总时间。计算出当前帧与下一帧中心点位置的距离后，通过视频中像素点之间的距离与实际物体距离的比例换算得到奶牛的最终运动总量，如图9B所示。
表 2 不同算法目标跟踪的跟踪精度Table 2. Tracking accuracy of target tracking by different algorithms算法
Algorithm跟踪精度/%
Tracking accuracy鲁棒性/%
RobustnessDeep SORT 58.09 47.28 Staple 57.98 43.69 SiameseRPN 56.14 37.15 ![]()
图 9 奶牛运动轨迹图(A)和奶牛识别图(B)Figure 9. Trajectory map of cow movement (A) and cow identification map (B)2.4 对比试验结果分析
本文对模型进行了消融试验，为了验证本文提出的BCE-YOLO改进方法提升模型性能的有效性，以YOLOv8模型为基础，通过使用不同的模块来验证模型性能。试验表明，改进后的方法对YOLOv8模型的检测效果有一定程度的影响。接下来，本文分别对前方和上方拍摄奶牛采食行为数据集进行测试，如表3、4所示，对试验数据进行分析，以反映BCE-YOLO模型中不同改进方法对模型检测性能的影响。
表 3 前方拍摄奶牛采食数据集中各模块对模型的影响Table 3. The influence of each module on the model in dairy cow feed dataset taken in frontYOLOv8 BiFormer CoT EMA 精确度/%
Precision召回率/%
Recall平均精确度均值/%
mAP√ 69.76 75.72 75.90 √ √ 70.56 77.56 77.53 √ √ 67.14 81.76 77.82 √ √ 69.75 72.23 73.93 √ √ √ 69.40 85.65 76.64 √ √ √ 67.54 75.80 75.76 √ √ √ 64.32 86.77 77.66 √ √ √ √ 77.73 82.57 83.70 表 4 上方拍摄奶牛采食数据集中各模块对模型的影响Table 4. The impact of each module on the model in dairy cow feed dataset taken aboveYOLOv8 BiFormer CoT EMA 精确度/%
Precision召回率/%
Recall平均精确度均值/%
mAP√ 68.29 77.16 67.23 √ √ 70.93 82.35 72.36 √ √ 73.05 80.52 73.00 √ √ 69.50 80.00 67.63 √ √ √ 73.54 85.92 73.82 √ √ √ 73.24 84.41 72.59 √ √ √ 71.45 82.60 69.47 √ √ √ √ 76.32 86.33 76.81 从表3、4数据可以看出，单模块改进YOLOv8网络时，BiFormer模块对YOLOv8的性能提升最大，无论是前方还是上方数据集，精确度、召回率、平均精确度均值均有所提升。其中前方数据集中，精确度提升0.80个百分点，召回率提升1.84个百分点，平均精确度均值提升1.63个百分点；上方数据集中，精确度提升2.64个百分点，召回率提升5.19个百分点，平均精确度均值提升5.13个百分点。单独使用CoT模块时，模型的召回率和平均精确度均值提升较大，前方数据集的召回率提升6.04个百分点，平均精确度均值提升1.92个百分点；上方数据集的召回率提升3.36个百分点，平均精确度均值提升5.77个百分点。在精确度方面，前方数据集在使用CoT模块后下降了2.62个百分点，上方数据集提升4.76个百分点。本文认为是前方拍摄的奶牛采食行为数据集中奶牛面部特征较为丰富，单一的注意力层对特征的提取并不充分，使得提升效果较差；而上方拍摄的奶牛采食行为数据集中奶牛头部特征相对较少，模型提取特征相对容易。同样的情况也发生在EMA模块中。在前方拍摄数据集中，增加EMA模块后召回率下降3.49个百分点，平均精确度均值下降1.97个百分点；而在上方拍摄数据集中精确度、召回率均有所提升，其中精确度提升1.21个百分点，召回率提升2.84个百分点，平均精确度均值基本保持不变。
在双模块改进的YOLOv8中发现，包含BiFormer模块的双模块改进网络相较于原单模块改进网络均有提升。增加CoT和EMA双模块，前方拍摄数据集的精确度相较于YOLOv8模型降低5.44个百分点，上方拍摄数据集精确度的提升幅度也相对较小。可以看出，增加BiFormer模块可以对模型有较大的提升。
当使用3个模块对YOLOv8进行改进时，精确度、召回率、平均精确度均值均提升明显。前方拍摄数据集中，精确度、召回率和平均精确度均值分别提升7.97、6.85和7.80个百分点；上方拍摄数据集中，精确度、召回率和平均精确度均值分别提升8.03、9.17和9.58个百分点。
从以上分析可以看出，改进模块中，BiFormer模块对模型的贡献度最大，增加3个模块后的BCE-YOLO模型性能相较于YOLOv8模型有较大提升，对奶牛采食行为特征的提取也更加充分。
2.5 图像检测对比分析
为了直观分析BCE-YOLO模型相较于YOLOv8的提升情况，本文使用前方和上方拍摄的奶牛采食行为测试数据集对BCE-YOLO和YOLOv8分别进行测试(图10)。前方数据集中，咀嚼和采食动作的置信度提升6%~10%，提升效果明显；上方数据集中，咀嚼动作置信度提升9%，采食动作提升1%；YOLOv8存在漏检现象。可以看出，前方检测效果优于上方。YOLOv8与BCE-YOLO模型在500张验证集图像上的结果如表5所示，在前方和上方拍摄的数据集中，BCE-YOLO模型精确度、召回率和平均精确度均值都占有优势，并且前方检测效果优于上方。分析原因是因为前方数据集的奶牛头部特征丰富，而且BCE-YOLO提取特征能力有较大提升，使得前方检测效果更强；相对来说，上方数据集奶牛头部特征较少，模型的特征提取作用未发挥完全，使得检测效果较前方差一些。动作类别中，咀嚼类别的置信度要高于采食类别，与图11得到的精确度曲线结果一致。本文分析认为，奶牛进行咀嚼动作时需要抬头，使得奶牛头部的特征更加突出，而采食动作时奶牛与饲料直接接触，使得部分饲料的特征也被提取从而使得置信度降低。
表 5 模型在验证集上的检测结果对比Table 5. Comparison of the detection results of the models on the validation set% 模型
Model精确度
Precision召回率
Recall平均精确度均值
mAP前方
Front上方
Above前方
Front上方
Above前方
Front上方
AboveYOLOv8 71.84 69.97 75.69 76.52 76.33 65.67 BCE-YOLO 79.66 79.23 83.79 82.79 86.94 83.67 2.6 采食行为分类检测分析
本研究将奶牛采食行为分为采食和咀嚼两部分，而采食可以分为采食和拱草两部分。本文从精确度角度分析奶牛采食行为分类对检测结果的影响。由于前方与上方拍摄数据集训练时数据增长趋势基本一致，因此只选择上方拍摄数据集训练数据进行分析。
精确度曲线如图11所示，图11A为将奶牛采食行为分为采食和咀嚼，图11B为将奶牛采食行为分为采食、拱草和咀嚼。可以看出，当分类中增加了拱草类别后，咀嚼类别的精确度没有受到影响，在置信度50%~70%时基本稳定在85%~90%，但采食类别精确度出现了下降，平均下降5个百分点，并且拱草类别的精确度只有20%，相对偏低，使得平均精确度下降了6个百分点。本文分析认为，由于拱草行为与采食行为的特征相似度高，且拱草行为一般在极短时间内完成，模型容易将拱草行为与采食行为判定为同一种行为，因此导致模型检测拱草行为的精确度不高，且使采食行为精确度下降。因此本文决定将不使用拱草类别进行测试。
![]()
图 11 不同分类下模型精确度曲线Figure 11. Precision curves of the model under different classifications3. 结论
本试验以养殖场饲养环境下采食区的奶牛为研究对象，对奶牛的采食行为进行研究。为解决复杂饲养环境下奶牛采食行为识别精度差、特征提取不充分的问题，实现对奶牛采食行为的自动监控，提出一种改进BCE-YOLO模型、并结合Deep SORT算法的识别方法。自建奶牛采食行为数据集进行对比试验，得出如下结论。
1)使用BiFormer、CoT、EMA 3个模型增强模块对YOLOv8模型进行增强，提出改进的BCE-YOLO模型，在前方和上方拍摄的数据集中，精确度分别为77.73%、76.32%，召回率分别为82.57%、86.33%，平均精确度均值分别为83.70%、76.81%，相较于YOLOv8模型，模型整体性能提升6~8个百分点。
2)通过将BCE-YOLO模型与Deep SORT结合，对奶牛采食行为进行了准确跟踪且有效抑制了奶牛ID的变更，可以同步实现多奶牛采食行为的自动监督与记录，为奶牛的福利化养殖发展具有积极的影响。
-
![]()
图 1 稻草钾肥配施对稻草C、N、P和K含量的影响
各图中，相同指标柱子上方的不同小写字母表示差异显著(P<0.05，Duncan’ s法)
Figure 1. Influences of combined application of straw and K fertilizer on C, N, P, K contents of rice straw
In each figure, different lowercase letters on bars of the same indicator indicated significant difference(P<0.05, Duncan’ s test)
![]()
图 2 稻草钾肥配施对稻草腐解率及C、N、P和K累积释放率的影响
相同指标柱子上方的不同小写字母表示差异显著(P<0.05，Duncan’ s法)
Figure 2. Influences of combined application of paddy-rice straw and K fertilizer on paddy-rice decomposition rates and accumulative release rates of C, N, P, K
Different lowercase letters on bars of the same indicator indicated significant difference(P<0.05, Duncan’ s test)
表 1 盆栽试验各处理的施肥和覆盖稻草量
Table 1 Amounts of fertilizer application and straw mulching in pot experiments
kg·hm−2 处理
TreatmentN P2O5 K2O 稻草
Paddy-rice strawK0S0(CK) 160.7 36.9 0 0 K0S2 160.7 36.9 0 6 000 K1S2 160.7 36.9 99.8 6 000 K2S2 160.7 36.9 199.7 6 000 K3S2 160.7 36.9 299.5 6 000 K2S0 160.7 36.9 199.7 0 K2S1 160.7 36.9 199.7 3 000 K2S3 160.7 36.9 199.7 9 000 K1S1 160.7 36.9 99.8 3 000 
下载: 导出CSV
表 2 马铃薯收获时稻草养分累积释放量1)
Table 2 Accumulative release amounts of nutrients at harvesting period of potato
kg·hm−2 处理 Treatment C N P K K0S2 1 528.69±19.86b 25.76±0.14b 3.55±0.02b 102.32±0.01b K1S2 1 508.22±140.85b 26.60±3.96b 3.65±0.31b 102.35±0.39b K2S2 1 573.74±109.86b 27.39±2.30b 3.58±0.19b 102.35±0.23b K3S2 1 456.05±75.99b 24.39±1.96b 3.15±0.28b 101.47±0.48b K2S1 764.81±19.93c 12.23±0.69c 1.70±0.07c 51.07±0.12c K2S3 2 221.43±135.76a 36.58±2.90a 4.94±0.28a 152.74±0.46a K1S1 723.99±17.39c 11.13±0.64c 1.71±0.03c 51.03±0.02c 1)表中数据为3次重复的平均值±标准误，同列数据后的不同小写字母表示差异显著(P<0.05，Duncan’ s法)
1) Data were averages of triplicate samples ± standard deviations (n = 3), and different lowercase letters in the same column indicated significant difference(P<0.05, Duncan’ s test)
下载: 导出CSV
表 3 稻草钾肥配施对土壤N平衡的影响1)
Table 3 Influences of combined application of paddy-rice straw and K fertilizer on soil N balance
处理
Treatment投入量/(kg·hm−2) Input amount 支出量/
(kg·hm−2)
Output amount表观盈亏量/(kg·hm−2)
Apparent surplus and deficient amount表观平衡系数
Apparent balance coefficient实际平衡率/%
Actual equilibrium rate化学N
Chemical N稻草N
Straw NK0S0 160.65 0 89.20a 71.45±1.34e 1.80 80.19 K0S2 160.65 43.35 74.43b 129.57±3.62b 2.76 175.42 K1S2 160.65 43.35 88.63a 115.37±7.56bc 2.33 133.34 K2S2 160.65 43.35 73.63ab 130.37±4.12b 2.78 178.77 K3S2 160.65 43.35 77.41ab 126.59±9.37b 2.72 171.59 K2S0 160.65 0 83.97ab 76.68±0.29e 1.91 91.33 K2S1 160.65 21.68 82.01ab 100.31±1.54d 2.23 122.47 K2S3 160.65 65.03 75.40ab 150.28±3.78a 3.01 200.91 K1S1 160.65 21.68 80.42ab 101.90±2.54cd 2.27 127.15 1)表中数据为3次重复的平均值±标准误，同列数据后的不同小写字母表示差异显著(P<0.05，Duncan’ s法)
1) Data were averages of triplicate samples ± standard deviations (n = 3), and different lowercase letters in the same column indicated significant difference(P<0.05, Duncan’ s test)
下载: 导出CSV
表 4 稻草钾肥配施对土壤P平衡的影响1)
Table 4 Influences of combined application of paddy-rice straw and K fertilizer on soil phosphorus balance
处理
Treatment投入量/(kg·hm−2) Input amount 支出量/
(kg·hm−2)
Output amount表观盈亏量/(kg·hm−2)
Apparent surplus and deficient amount表观平衡系数
Apparent balance coefficient实际平衡率/%
Actual equilibrium rate化学P
Chemical P稻草P
Straw PK0S0 16.11 0 13.53a 2.58±0.07d 1.19 19.09 K0S2 16.11 5.0 12.77a 8.33±1.33ab 1.69 50.60 K1S2 16.11 5.0 15.01a 6.1±1.18bc 1.42 42.32 K2S2 16.11 5.0 11.47a 9.64±1.76a 1.95 59.84 K3S2 16.11 5.0 12.77a 8.34±1.59ab 1.71 49.99 K2S0 16.11 0 12.87a 3.24±0.78cd 1.26 26.08 K2S1 16.11 2.5 13.07a 5.54±0.45bcd 1.43 42.74 K2S3 16.11 7.5 12.98a 10.63±0.48a 1.82 82.34 K1S1 16.11 2.5 12.48a 6.13±0.77bc 1.50 50.21 1)表中数据为3次重复的平均值±标准误，同列数据后的不同小写字母表示差异显著(P<0.05，Duncan’ s法)
1) Data were averages of triplicate samples ± standard deviations (n = 3), and different lowercase letters in the same column indicated significant difference(P<0.05, Duncan’ s test)
下载: 导出CSV
表 5 稻草钾肥配施对土壤K平衡的影响1)
Table 5 Influences of combined application of paddy-rice straw and K fertilizer on soil potassium balance
处理
Treatment投入量/(kg·hm−2) Input amount 支出量/
(kg·hm−2)
Output amount表观盈亏量/(kg·hm−2)
Apparent surplus and deficient amount表观平衡系数
Apparent balance coefficient实际平衡率/%
Actual equilibrium rate化学K
Chemical K稻草K
Straw KK0S0 0 0 163.80b −163.80±2.36f 0 −100.00 K0S2 0 86.60 177.01ab −90.40±13.68e 0.49 −50.64 K1S2 83.19 86.60 214.60a −44.81±16.21d 0.80 −20.26 K2S2 166.38 86.60 178.28ab 74.70±23.28b 1.45 45.39 K3S2 249.56 86.60 179.59ab 156.58±28.92a 1.95 94.73 K2S0 166.38 0 186.66ab −20.28±10.75cd 0.90 −10.47 K2S1 166.38 43.30 207.44ab 2.24±5.09c 1.01 1.16 K2S3 166.38 129.90 202.50ab 93.78±4.32b 1.46 46.40 K1S1 83.19 43.30 183.02ab −56.53±13.26de 0.70 −30.41 1)表中数据为3次重复的平均值±标准误，同列数据后的不同小写字母表示差异显著(P<0.05，Duncan’ s法)
1) Data were averages of triplicate samples ± standard deviations (n = 3), and different lowercase letters in the same column indicated significant difference(P<0.05, Duncan’ s test)
下载: 导出CSV
-
[1] 屈冬玉. 马铃薯产业应在我国农村区域经济发展中发挥更大作用[C]//中国马铃薯作物学会马铃薯专业委员会. 马铃薯产业与农村区域发展. 北京: 中国作物学会, 2013: 1-3. [2] 谢开云, 屈冬玉, 金黎平, 等. 中国马铃薯生产与世界先进国家的比较[J]. 世界农业, 2008, 29(5): 35-38. doi: 10.3969/j.issn.1002-4433.2008.05.012 [3] 米健, 罗其友, 高明杰. 南方冬作区马铃薯发展趋势、区域格局与增产潜力[J]. 中国农业资源与区划, 2012, 33(3): 73-79. [4] 屈冬玉, 金黎平, 谢开云. 中国马铃薯产业10年回顾[M]. 北京: 中国农业科学技术出版社, 2010: 66-69. [5] 高媛, 韦艳萍, 樊明寿. 马铃薯的养分需求[J]. 中国马铃薯, 2011, 25(3): 182-187. doi: 10.3969/j.issn.1672-3635.2011.03.018 [6] MULDER A. Potato diseases-diseases, pests and defects[M]. Holland: Aardappelwereld B V and NIVAP, 2005: 209-221.
[7] GEORGE M S, LU G, ZHOU W. Genotypic variation for potassium uptake and utilization efficiency in sweet potato (Ipomoea batatas L.)[J]. Field Crop Res, 2002, 77(1): 7-15. doi: 10.1016/S0378-4290(02)00043-6
[8] 白由路. 高价格下我国钾肥的应变策略[J]. 中国土壤与肥料, 2009, 45(3): 1-4. doi: 10.3969/j.issn.1673-6257.2009.03.001 [9] 陈杨, 樊明寿, 李斐, 等. 氮素营养诊断技术的发展及其在马铃薯生产中的应用[J]. 中国农学通报, 2009, 25(3): 66-71. [10] 张新明, 张洪秀, 李水源, 等. 恩平市典型种植户冬作马铃薯施肥状况调查分析[J]. 安徽农业科学, 2011, 39(36): 22286-22288. doi: 10.3969/j.issn.0517-6611.2011.36.034 [11] LAL R. World crop residues production and implications of its use as a biofuel[J]. Environ Int, 2005, 31(4): 575-584. doi: 10.1016/j.envint.2004.09.005
[12] 毕于运, 王红彦, 王道龙, 等. 中国稻草资源量估算及其开发利用[J]. 中国农学通报, 2011, 27(15): 137-143. [13] 李国学, 张福锁. 固体废物堆肥化和有机复混肥的生产[M]. 北京: 化学工业出版社, 2000: 1-20. [14] 曹国良, 张小曳, 王亚强, 等. 中国区域农田秸秆露天焚烧排放量的估算[J]. 科学通报, 2007, 52(15): 1826-1831. doi: 10.3321/j.issn:0023-074x.2007.15.017 [15] 谭德水, 金继运, 黄绍文, 等. 施钾和秸秆还田对栗钙土区土壤养分及小麦产量的长期效应研究[J]. 干旱地区农业研究, 2009, 27(2): 194-198. [16] 张静, 温晓霞, 廖允成, 等. 不同玉米秸秆还田量对土壤肥力及冬小麦产量的影响[J]. 植物营养与肥料学报, 2010, 16(3): 612-619. doi: 10.11674/zwyf.2010.0314 [17] CARUSO G, CARPUTO D, CONTI S, et al. Effect of mulching and plant density on out-of-season organic potato growth, yield and quality[J]. Adv Hortic Sci, 2013, 27(3): 115-121.
[18] DVORAK P, KUCHTOVA P, TOMASEK J. Response of surface mulching of potato (Solanum tuberosum) on SPAD value, colorado potato beetle and tuber yield[J]. Int J Agric Biol, 2013, 15(4): 798-800.
[19] 农业部办公厅. 农业部关于印发《测土配方施肥技术规范(2011年修订版)》的通知: 农农发[2011]3号[A].(2011-09-22)[2018-11-02]. http://jiuban.moa.gov.cn/zwllm/tzgg/tz/201109/t20110922_2293389.htm. [20] 鲍士旦. 土壤融化分析[M]. 3版. 北京: 中国农业出版社, 2000. [21] 王志勇, 白由路, 杨俐苹, 等. 低土壤肥力下施钾和秸秆还田对作物产量及土壤钾素平衡的影响[J]. 植物营养与肥料学报, 2012, 18(4): 900-906. [22] TANG Q, ZHANG C. Data processing system (DPS) software with experimental design, statistical analysis and data mining developed for use in entomological research[J]. Insect Sci, 2013, 20(2): 254-260. doi: 10.1111/j.1744-7917.2012.01519.x
[23] 闫超. 水稻秸秆还田腐解规律及对土壤养分和酶活性的影响[D]. 哈尔滨: 东北农业大学, 2012. [24] 李世清, 李生秀. 有机物料在维持土壤微生物体氮库中的作用[J]. 生态学报, 2001, 21(1): 136-142. doi: 10.3321/j.issn:1001-9332.2001.01.043 [25] ABRO S A, TIAN X H, WANG X D, et al. Decomposition characteristics of maize (Zea mays L.) straw with different carbon to nitrogen (C/N) ratios under various moisture regimes[J]. Afr J Biotechnol, 2011, 10(50): 10149. doi: 10.5897/AJB10.2261
[26] SUN B, WANG X, WANG F, et al. Assessing the relative effects of geographic location and soil type on microbial communities associated with straw decomposition[J]. Appl Environ Microbiol, 2013, 79(11): 3327-3335. doi: 10.1128/AEM.00083-13
[27] LU P, LIN Y H, YANG Z Q, et al. Effects of application of corn straw on soil microbial community structure during the maize growing season[J]. J Basic Microbiol, 2015, 55(1): 22-32. doi: 10.1002/jobm.201300744
[28] 李逢雨, 孙锡发, 冯文强, 等. 麦秆、油菜秆还田腐解速率及养分释放规律研究[J]. 植物营养与肥料学报, 2009, 15(2): 374-380. doi: 10.3321/j.issn:1008-505X.2009.02.018 [29] 戴志刚, 鲁剑巍, 李小坤, 等. 不同作物还田秸秆的养分释放特征试验[J]. 农业工程学报, 2010, 26(6): 272-276. doi: 10.3969/j.issn.1002-6819.2010.06.047 [30] 张晋京, 窦森. 玉米秸秆分解期间胡敏酸、富里酸动态变化的研究[J]. 土壤通报, 2005(1): 134-136. doi: 10.3321/j.issn:0564-3945.2005.01.037 [31] 张亮. 关中麦玉轮作区施氮对秸秆还田小麦产量和秸秆养分释放的影响[D]. 杨凌: 西北农林科技大学, 2012. [32] 鲁如坤, 时正元, 施建平. 我国南方6省农田养分平衡现状评价和动态变化研究[J]. 中国农业科学, 2000, 33(2): 63-67. doi: 10.3321/j.issn:0578-1752.2000.02.010 [33] 鲁如坤, 刘鸿翔, 闻大中, 等. 我国典型地区农业生态系统养分循环和平衡研究:Ⅳ:农田养分平衡的评价方法和原则[J]. 土壤通报, 1996, 27(5): 197-199. [34] 谭德水, 金继运, 黄绍文, 等. 不同种植制度下长期施钾与秸秆还田对作物产量和土壤钾素的影响[J]. 中国农业科学, 2007, 40(1): 133-139. doi: 10.3321/j.issn:0578-1752.2007.01.018 -
期刊类型引用(1)
1. 许凯，马英，陈三雄，骆金初，黄伟城，郭微，王龙远. 广东德庆县香山森林公园黄牛木种群特征研究. 安徽农业科学. 2022(13): 125-129 . 
百度学术
其他类型引用(0)
计量
- 文章访问数: 1391
- HTML全文浏览量: 3
- PDF下载量: 1346
- 被引次数: 1
