• 《中国科学引文数据库(CSCD)》来源期刊
  • 中国科技期刊引证报告(核心版)期刊
  • 《中文核心期刊要目总览》核心期刊
  • RCCSE中国核心学术期刊

不同磷浓度对土壤理化性质及防风生长和药材品质的影响

郝佳, 刘宇航, 殷洁, 刘晶, 田新, 韩忠明, 王云贺, 韩梅, 杨利民

郝佳, 刘宇航, 殷洁, 等. 不同磷浓度对土壤理化性质及防风生长和药材品质的影响[J]. 华南农业大学学报, 2022, 43(3): 59-67. DOI: 10.7671/j.issn.1001-411X.202104024
引用本文: 郝佳, 刘宇航, 殷洁, 等. 不同磷浓度对土壤理化性质及防风生长和药材品质的影响[J]. 华南农业大学学报, 2022, 43(3): 59-67. DOI: 10.7671/j.issn.1001-411X.202104024
HAO Jia, LIU Yuhang, YIN Jie, et al. Effects of different phosphorus concentrations on soil physicochemical property, growth and medicinal materials quality of Saposhnikovia divaricata[J]. Journal of South China Agricultural University, 2022, 43(3): 59-67. DOI: 10.7671/j.issn.1001-411X.202104024
Citation: HAO Jia, LIU Yuhang, YIN Jie, et al. Effects of different phosphorus concentrations on soil physicochemical property, growth and medicinal materials quality of Saposhnikovia divaricata[J]. Journal of South China Agricultural University, 2022, 43(3): 59-67. DOI: 10.7671/j.issn.1001-411X.202104024

不同磷浓度对土壤理化性质及防风生长和药材品质的影响

基金项目: 国家重点研发计划(2019YFC1710700, 2019YFC1710702);吉林省科技发展计划(20200404010YY, 20210204011YY);中央财政林业科技推广示范项目(JLT2020-22, JLT2021-03); 财政部和农业农村部:国家现代农业产业技术体系资助(CARS-21); 吉林省教育厅“十三五”科学技术研究(JJKH20190936KJ)
详细信息
    作者简介:

    郝佳,硕士研究生,主要从事野生动植物保护与利用研究,E-mail: 1719809814@qq.com

    通讯作者:

    韩忠明,教授,博士,主要从事中药资源生态学研究,E-mail: hanzm2008@126.com

  • 中图分类号: R282.2

Effects of different phosphorus concentrations on soil physicochemical property, growth and medicinal materials quality of Saposhnikovia divaricata

  • 摘要:
    目的 

    探讨不同磷浓度处理对土壤特性及防风Saposhnikovia divaricata生长状况、色原酮含量的影响,为防风人工栽培的技术策略制定以及防风对低磷胁迫响应机制的研究提供理论依据。

    方法 

    以2年生防风为材料,设置营养液磷 (NH4H2PO4) 浓度分别为1.0、0.1和0 (无磷) mmol/L,分别在30、60、90 d时测定土壤理化性质及防风生长特性和4种色原酮总含量,分析不同磷浓度处理防风各指标与4种色原酮总含量的关系。

    结果 

    无磷处理下,2年生防风根粗、根生物量在整个处理期间均低于1.0 mmol·L−1磷处理,根长在90 d 时显著高于1.0和0.1 mmol·L−1磷处理;土壤pH 在60 d 后表现为1.0 mmol·L−1磷处理显著高于0.1 mmol·L−1和无磷处理;土壤有机质含量在整个试验期表现为1.0和0.1 mmol·L−1磷处理显著低于无磷处理;90 d 时, 无磷处理的土壤碱解氮含量显著高于1.0和0.1 mmol·L−1磷处理,而土壤有效磷含量0.1 mmol·L−1处理显著高于无磷处理,4 种色原酮总含量0.1 mmol·L−1和无磷处理显著高于1.0 mmol·L−1磷处理。防风根系土壤中酸性磷酸酶活性在30 d 时,3种磷处理间存在显著差异。1.0 mmol·L−1磷处理下,根长、根鲜生物量与4 种色原酮总含量呈正相关;0.1 mmol·L−1磷处理下,碱解氮、有效磷含量和根鲜生物量均与4 种色原酮总含量呈正相关;无磷处理下pH、碱解氮含量、酸性磷酸酶活性均与4种色原酮总含量呈正相关。

    结论 

    防风面对磷胁迫有一系列响应机制,有效成分4种色原酮总含量显著增加,pH、有机质、碱解氮等与酸性磷酸酶协同作用,保证防风对养分的吸收利用。

    Abstract:
    Objective 

    To discuss the effects of different phosphorus concentrations on soil characteristics, growth status and chromone content of Saposhnikovia divaricata, and provide theoretical basis for technical strategy formulation of S. divaricata artificial cultivation and responding mechanism research of S. divaricata to low phosphorus stress.

    Method 

    Two-year-old S. divaricata was used as material, and three phosphorus treatments with NH4H2PO4 concentration of 1.0, 0.1 and 0 (no phosphorus) mmol/L respectively were set up. The soil physical and chemical properties, and the growth characteristics and total contents of four kinds of chromones inS. divaricata were determined on the 30th, 60th, and 90th day to analyze the correlativity between the various indicators and total chromone content under different phosphorus concentrations.

    Result 

    The root diameter and weight of 2-year-old S. divaricata in no phosphorus treatment were lower than those of 1.0 mmol·L−1 phosphorus treatment during the whole treatment period, and the root length after 90 d of treatment was significantly higher than those of 1.0 and 0.1 mmol·L−1 phosphorus treatments. After 60 d of treatment, the soil pH of 1.0 mmol·L−1 phosphorus treatment was significantly higher than those of 0.1 mmol·L−1 and no phosphorus treatments. The contents of soil organic matter in 1.0 and 0.1 mmol·L−1 phosphorus treatments were significantly lower than that in no phosphorus treatment during the whole treatment period. After 90 d of treatment, the alkali-hydrolyzable nitrogen content in soil of no phosphorus treatment was significantly higher than those of 1.0 and 0.1 mmol·L−1 phosphorus treatments, whereas the available phosphorus content in soil of 0.1 mmol·L−1 phosphorus treatment was significantly higher than that of no phosphorus treatment, the total content of four chromones of 0.1 mmol·L−1 and no phosphorus treatments were significantly higher than that of 1.0 mmol·L−1 phosphorus treatment. On the 30 th day of treatment, the activities of acid phosphatase in soil among three treatments were significantly different. There were positive correlations between root diameter/root fresh biomass and chromone content in 1.0 mmol·L−1 phosphorus treatment, between alkali-hydrolyzable nitrogen content/available phosphorus content/root fresh biomass and chromone content in 0.1 mmol·L−1 phosphorus treatment, between pH/alkali-hydrolyzable nitrogen content /acid phosphate activity and chromone content in no phosphorus treatment, respectively.

    Conclusion 

    S. divaricata has a series of response mechanisms to phosphorus stress, including chromone content significantly increasing, synergistic interaction of pH, organic matter, alkali-hydrolyzable nitrogen and acid phosphatase, which ensures the absorption and utilization of nutrients by S. divaricata.

  • 中国是一个农业大国,作为国家发展的根本,农业在国家经济总量中占有很大的比重。如今全球粮食供给偏紧,农业是国家战略性、基础性核心产业,在农作物的生长发育过程中,气候环境、病虫迫害、病菌侵染等严重影响其产量和品质。目前常见的农作物病害有上千种,传统的依靠人工识别农作物病害,不但耗费时间精力,而且效率低。精确迅速地辨识出农作物病害有助于及时救助,对降低农作物病害带来的产量和品质影响具有重大意义。

    随着深度学习技术在物体分类方面的深入研究与计算机图像处理能力的增强,计算机视觉技术也应用在农作物病害识别上,并且取得许多成果。刘翱宇等[1]提出一种基于深度残差网络的植物病害检测网络,引入Focal Loss损失函数,对植物病害图像的检测平均准确率高达97. 96%。Srdjan等[2]将深度卷积神经网络用于植物病害识别,能够区分植物叶片与周围环境,并能够区分健康叶片与13种不同的植物病害,在开发模型测试中,单类病害的识别平均准确率高达96.3%。Albattah等[3]提出了一种改进的CenterNet网络,以DenseNet-77为基础网络对病害样本图片特征进行提取,在PlantVillage数据集上训练,然后分别对番茄、葡萄等在内的14种植物26类病害及健康叶片进行识别,检测精度比当时精确率最高的EfficientNet网络更高。刘洋等[4]针对PlantVillage数据集14种植物中的26种病害对比研究了2个轻量级网络系统MobileNet和Inceptionv3,结果表明MobileNet模型的程序数量更小,执行速率也更快,可部署于安卓操作系统中,在手机端也实现了图像识别。Saleem等[5]对Xception网络精修后,选择性能提升最明显的Adam优化器,对PlantVillage数据集中的26种植物病害进行检测,准确率达到了97.81%。孙俊等[6]使用AlexNet作为基础网络,然后通过批归一化方法、引入全局池化层和减少特征图量等方式,共获得了8种改进模式,对PlantVillage数据集14种植物的26种病害特征进行了鉴别,效果最好的检测模式平均准确度超过了98.1%。李书琴等[7]提出了一种基于轻量级残差网络,通过缩减网络卷积核数目和轻量级残差模块减少网络参数,网络参数大幅减少的同时,在PlantVillage数据集上的检测精度达到98.45%。

    上述深度学习方法不仅在检测准确度方面有了较高的提升,也有向低内存、少参数量改进的趋势,从而达到高精度识别、迅速检测、易部署移动端的目的。然而,检测精度提升的同时,网络参数量也在不断增多,网络深度和宽度也不断加深,浮点运算量加大从而带来较高的延迟,实时检测的效果较差。人工智能公司旷视科技提出的YOLOX-Nano[8],有着参数量小、浮点运算快、延迟低、易部署移动端的优势且检测精度保持着中高水平。本文在原始YOLOX-Nano网络结构的基础上进行改进,在不显著增加网络内存和运算参数量的基础上,提出一种混合卷积注意力模块CBAM[9]、Mixup数据增强策略[10]、Focal Loss分类损失函数[11]和本文设计的基于CenterIOU回归损失函数的改进YOLOX-Nano模型。

    试验所需要的所有病害图片数据来自于公开数据集PlantVillage网站(http://plantvillage.psu.edu),数据集记录包含54309张图像,这些图像涵盖了14种作物:苹果、蓝莓、樱桃、玉米、葡萄、柑橘、桃、甜椒、马铃薯、覆盆子、大豆、南瓜、草莓和番茄。考虑到应尽可能研究某一类农作物的多种病害(3种及以上),因此剔除数据集内仅含健康叶片的农作物数据和仅有2类病害的农作物数据,最终在公开数据集上精心挑选苹果叶片图像4645张,其中,苹果黑星病1000张、苹果黑腐病1000张、苹果松锈病1000张、健康叶片1645张;玉米叶片4354张,其中,玉米灰斑病1000张、玉米锈病1192张、玉米大斑病1000张、健康叶片1162张;葡萄叶片4369张,其中,葡萄黑腐病1180张、葡萄轮斑病1076张、葡萄褐斑病1383张、健康叶片1000张;马铃薯叶片3000张,其中,马铃薯早疫病、马铃薯晚疫病、健康叶片各1000张;草莓叶片4218张,其中,草莓叶焦病2218张,健康叶片为2000张;番茄叶片18835张,其中,番茄疮痂病1404张、番茄早疫病1000张、番茄晚疫病1909张、番茄叶霉病1000张、番茄斑枯病1771张、番茄红蜘蛛损伤1676张、番茄斑点病2127张、番茄黄叶曲叶病5357张、番茄花叶病毒病1000张和健康叶片1591张,共计6类农作物27种病害,保证各类样本间数量均衡。此外,考虑到真实场景背景复杂的情况,数据集添加了无叶片背景图像1143张。试验样本图像如图1所示。

    图  1  PlantVillage部分数据集
    1~4 依次为苹果黑星病、黑腐病、松锈病和健康叶;5~8依次为玉米灰斑病、锈病、健康叶和大斑病;9~12依次为葡萄黑腐病、轮斑病、褐斑病和健康叶;13~15依次为马铃薯早疫病、晚疫病和健康叶;16~17 依次为草莓叶焦病和健康叶;18~27 依次为番茄疮痂病、早疫病、晚疫病、叶霉病、斑枯病、红蜘蛛损伤、斑点病、黄叶曲叶病、花叶病毒病和健康叶;28为无叶片背景图
    Figure  1.  Partial dataset of PlantVillage
    1−4 are apple scab, apple black rot, cedar apple rust and apple healthy leaves in turn; 5−8 are corn Cercospora zeaemaydis Tehon and Daniels, corn Puccinia polysora, corn healthy leaves and corn Curvularia leaf spot fungus; 9−12 are grape black rot, grape black measles, grape leaf blight fungus and grape healthy leaves in turn; 13−15 are potato early blight, potato late blight and potato healthy leaves in turn; 16−17 are strawberry leaf scorch and healthy leaves in turn; 18−27 are tomato bacterial spot bacteria, tomato early blight, tomato late blight water mold, tomato leaf mold, tomato septorial leaf spot, tomato spider mite damage, tomato target spot bacteria, tomato YLCV virus, tomato ToMV disease and tomato healthy leaves in turn; 28 is bladeless background image

    首先将所有样本图像的分辨率转变成416×416,使用标注工具LabelImg将所有样本按Pascal VOC数据集的标注,生成.xml类型的标注文件。此外,从试验数据集图像中随机选择60%样本(约23235张)用于训练模型,20%的样本(约7675张)用来验证,20%样本(约7675张)用于测试。

    试验运用迁移学习,让YOLOX-Nano在大规模公开数据集Image Net上进行训练,获得一个收敛的预训练权重,然后利用这个预训练权重迁移到改进的YOLOX-Nano网络模型中进行参数调优。训练过程中超参数设置为每批量样本数为32,1次遍历完全部训练集数据称为1次迭代,批量设置为100。采用随机梯度下降(Stochastic gradient descent,SGD)优化算法调优网络参数,初始学习率设置为0.01,动量因子为0.9,每经过10次迭代训练,将学习率降低1/10,模型迭代1次保存1次权重。

    YOLOX-Nano作为YOLOX系列的最轻量化版本,对比YOLOX-X、YOLOX-L、YOLOX-M、YOLOX-S、YOLOX-Darknet和YOLOX-Tiny高性能版本,其通过降低网络的宽度和深度、减少模型运算的参数量、取消训练时Mixup数据增强方式等改进,使网络结构更简单,速度更快。

    YOLOX-Nano的结构如图2a所示,整个YOLOX-Nano可以分为3个部分,分别是CSPDarkNet主干提取网络[12]、特征金字塔网络(Feature pyramid network,FPN)[13]和路径聚合网络(Path aggregation network,PAN)[14]中间层以及YoloHead输出层。

    图  2  YOLOX-Nano网络结构和改进后的YOLOX-Nano网络结构
    Figure  2.  Network structure of YOLOX-Nano and improved YOLOX-Nano

    输入的图片首先在CSPDarkNet进行特征提取,提取到3个特征层进行下一步网络的构建,FPN是YOLOX-Nano的加强特征提取网络,在主干部分获得的3个特征层会在这一部分进行特征融合,获取不同尺度的特征信息。在FPN部分,已经获得的有效特征层被用于继续提取特征,在PAN层中,网络会对特征进行阶段上采样和下采样操作实现特征融合。YoloHead是YOLOX-Nano的分类器与回归器,输入图片通过CSPDarkNet、FPN和PAN层已经获得了3个加强特征层,每一个特征图看作每一个特征点的集合,而YoloHead会对特征点进行判断,判断特征点是否有物体与其对应。YOLOX-Tiny的YoloHead被分为了2部分,分别实现检测和分类操作,最后在预测阶段将两者整合输出结果。

    由于YOLOX-Nano宽度和深度较浅,没有充分学习输入特征向量通道与空间之间的重要程度关系,对近似的难分辨目标预测效果较差,其次数据集内多种农作物病害图像样本复杂多样,二分类交叉熵损失(Binary cross entropy loss,BCELoss)作为分类损失函数将困难样本和容易样本赋予相同的权重进行训练,这样造成一些困难样本分类效果不理想,GIOU损失函数对2个候选框没有相交时,IOU=0,不能反映两者的距离大小,也无法精确反映两者的重合度大小,不利于梯度回传和学习训练,致使密集病害定位缺失。

    综合以上问题,本文在YOLOX-Nano模型混合了卷积注意力CBAM模块,引入了Mixup数据增强方式来充实样本集,让网络训练地更充分,最后分类损失函数将BCELoss替换为Focal Loss,回归损失函数将GIOU Loss替换为CenterIOU Loss,改进后的YOLOX-Nano网络结构图如图2b所示。

    在农作物病害检测与识别任务中,农作物叶片病害特征大体相似,只有细节部分略有差异。YOLOX-Nano作为浅层模型,检测精度不够高,这是需要注意的地方。

    注意力机制是运用在深度学习模型中常见的小技巧,是实现网络自适应注意的一个方式,核心重点是让网络关注到它需要关注的特征,忽略不重要的特征,主要分为通道注意力机制和空间注意力机制,卷积注意力机制结构如图3所示。卷积注意力机制的前半部分为通道注意力机制,通道注意力机制先针对输入进来的单个特征层,依次执行全局平均池化方法和全局最大池化方法,之后再针对平均池化和最大池化的结果,利用共享的全连接层进行处理,然后再对处理后的2个结果进行相加,之后再通过Sigmoid激活函数,得到输入特征层中每一条通道的加权值(0~1)。将这个权值乘上原输入特征层即可得到经过通道注意力机制处理后的特征图。卷积注意力机制的后半部分为空间注意力机制,首先对通过通道注意力机制处理的特征层每一个特征点的通道取最大值和平均值,之后将这2个结果进行堆叠,利用通道数为1的卷积核卷积调整通道数,再经过Sigmoid激活函数,获得输入特征层每一个特征点的权值(0~1),然后将这个权值乘上原输入特征层即可得到最终处理的特征层。

    图  3  CBAM卷积注意力机制结构
    Figure  3.  Structure of CBAM convolutional attention mechanism

    本文引入CBAM注意力机制模块,结合通道注意力机制和空间注意力机制,对输入的特征层分别进行通道注意力机制和空间注意力机制的处理,卷积注意力机制对不同通道的特征图进行权重的重新分配,加强网络结构深层信息提取,让网络关注到它需要关注的特征,忽略不重要的特征。

    输入图像经过主干特征提取网络后获得3个有效特征层,在构建FPN层上采样前首先将经由加强特征提取网络的3个有效特征层送入CBAM模块,自动学习特征通道之间的相关性,输出3个尺度的检测结果,再经过卷积层提取特征。随后进行PAN的构建,在PAN融合的下采样前,把主干特征网络获得的3个有效特征层,再经过CBAM注意力机制模块进行融合,然后经过另一个卷积模块再一次提取特征,将最后一个有效特征层卷积,最终获得3个特征。

    Mixup是一种混类增强的数据增强策略,它可以将不同类的图像进行混合,从而扩充训练数据集。 Mixup数据增强策略原理是假设${\text{batch}}_{x1}$是一个批次的图像样本,${\text{batch}}_{y1}$是该批次图像样本对应的标签;${\text{batch}}_{x2}$是另一个批次图像样本,${\text{batch}}_{y2}$是另一个批次图像样本对应的标签,λ是由参数为αβ的Beta分布计算出来的混合系数。

    由此可以得到Mixup原理公式为

    $$ { \lambda =}{\rm{Beta}}(\alpha ,\beta ) \text{,} $$ (1)
    $$ {\mathrm{m}\mathrm{i}\mathrm{x}\mathrm{e}\mathrm{d}\_\mathrm{b}\mathrm{a}\mathrm{t}\mathrm{c}\mathrm{h}}_{{x}}=\mathrm{\lambda }\times {{\rm{batch}}}_{x1}+(1-\mathrm{\lambda })\times {{\rm{batch}}}_{x2} \text{,} $$ (2)
    $$ {\mathrm{m}\mathrm{i}\mathrm{x}\mathrm{e}\mathrm{d}\_\mathrm{b}\mathrm{a}\mathrm{t}\mathrm{c}\mathrm{h}}_{{y}}=\mathrm{\lambda }\times {{\rm{batch}}}_{y1}+(1-\mathrm{\lambda })\times {{\rm{batch}}}_{y2} 。 $$ (3)

    以玉米病害图像为例,在训练时先读取1张样本图像,图像两侧填充,放大图片到416×416(图4a),再随机选取1张样本图像,也放大到416×416(图4b)。然后设置1个融合系数,如0.5,将2张图像加权融合,最终得到完整的图像(图4c)。

    图  4  Mixup数据增强策略示意图(以玉米灰斑病为例)
    Figure  4.  Schematic diagram of mixup data enhancement strategy(Take corn Cercospora zeaemaydis as an example )

    在YOLOX中以Mosaic为基础,Mixup策略作为一种额外的数据增强策略,在浅层模型用Mixup数据增强策略会降低平均精确率(Average precision,AP),所以YOLOX-Nano没有使用Mixup数据增强策略。本文使用Mixup数据增强策略,加大样本数据量防止模型欠拟合。根据图5中的Beta分布概率密度曲线,当α=β= 1时,等于(0,1)均匀分布;当α=β < 1时,表现为两头的概率大,中间的概率小,当α=β→0时,相当于{0,1}二项分布,要么取0,要么取1,等于原始数据没有增强,本试验随机设置Mixup参数α=β∈[0.1,0.9],就可以得到多样化的λ∈(0,1)区间内的概率分布,使数据增强更具有随机性,模型更具有鲁棒性。

    图  5  αβ相同时的Beta分布概率密度曲线
    Figure  5.  Probability density curve of beta distribution while α and β were equal

    一个好的损失函数对于网络参数调优达到事半功倍的效果。由于多种农作物叶片病害纹理特征大致形似,模型对一些病害样本计算的分类置信度产生偏差,进而影响样本分类效果,造成错误分类病害的情况。

    基于以上计算样本分类置信度偏差的问题,YOLOX-Nano将分类损失函数由BCE Loss损失函数替换为Focal Loss损失函数,其能控制容易分类和难分类样本的权重,通过减少易分类样本的权重,使得模型在训练时更专注于难分类的样本,让容易样本对损失产生贡献小,难样本对损失产生贡献大。

    Focal Loss分类损失函数公式定义如式(4)所示,$ \hat {p} $[$ \hat {p} $∈(0,1) ]表示预测值,$ y $表示真实值,$ \gamma $($ \gamma $≥0)表示调节因子。设$p_{\rm{t}}$表示预测某真实标签的概率,定义如式(5)所示,范围为[0,1]。由式(4)、(5)归纳后Focal Loss (FL)损失函数公式如式(6)所示:

    $$ {\rm{f}}\left(x\right)=\left\{\begin{array}{l}-{(1-\hat {p})}^{\gamma }\mathrm{ln}\left(\hat {p}\right), {\rm{if}}\;y=1\\ -{\hat {p}}^{\gamma }\mathrm{ln}\hat {p},{\rm{if}}\;y=0\end{array}\right. \text{,} $$ (4)
    $$ {\rm{f}}\left(x\right)=\left\{\begin{array}{l}\hat {p},{\rm{if}}\;y=1\\ 1-\hat {p},其他\end{array}\right. \text{,} $$ (5)
    $$ {\rm{FL}}\left({p}_{{\rm{t}}}\right)=-{(1-{p}_{{\rm{t}}})}^{\gamma }\mathrm{ln}{p}_{{\rm{t}}} 。 $$ (6)

    与交叉熵函数损失函数相比,公式(6)多了一个${(1-{p}_{{\rm{t}}})}^{\gamma }$系数,实质是对交叉熵损失函数的优化。如图6所示,可以通过改变$ \gamma $调节因子的值,从而控制容易分类和难分类样本对损失的贡献。由公式(6)可知,若某样本类别概率${ p_{\rm{t}} }$趋于1,说明预测值$\hat {p}$接近于真实值$ y $,这类样本就属于容易区分的样本,此时系数${(1-{p}_{{\rm{t}}})}^{\gamma }$趋于0,对损失贡献极少。

    图  6  Focal Loss函数曲线
    Figure  6.  Focal Loss function curve

    YOLOX-Nano在回归候选框时仍使用GIOU损失函数,本文根据数据集图片样本特点,设计出利于回归候选框的CenterIOU损失函数。

    图7所示,已知候选框中心点($ {x}_{\mathrm{p}\mathrm{r}\mathrm{e}} $$ {y}_{\mathrm{p}\mathrm{r}\mathrm{e}} $)、真实框中心点($ {x}_{\mathrm{t}\mathrm{a}\mathrm{g}} $$ {y}_{\mathrm{t}\mathrm{a}\mathrm{g}} $)与最小外接矩形中心点(${{x}}_{\mathrm{c}\mathrm{t}}$$ {y}_{\mathrm{c}\mathrm{t}} $)。根据公式(7)计算候选框和真实框的IOU,式中,${{\mathrm{B}\mathrm{O}\mathrm{X}}_{\mathrm{交}\mathrm{集}}}$${{\mathrm{B}\mathrm{O}\mathrm{X}}_{\mathrm{并}\mathrm{集}}} $分别表示候选框和真实框的交集和并集;构造候选框中心点和真实框中心点两点的一般直线方程$ Y = Ax + By + C $;根据公式(8)计算最小外接矩形中心点($ {x}_{\mathrm{c}\mathrm{t}} $$ {y}_{\mathrm{c}\mathrm{t}} $)到两框中心点的直线距离(Distance);最后代入公式(9),得到候选框回归损失。

    图  7  CenterIOU损失函数图解
    Figure  7.  Diagram of CenterIOU loss function
    $$ \mathrm{I}\mathrm{O}\mathrm{U}=\dfrac{{\mathrm{B}\mathrm{O}\mathrm{X}}_{\mathrm{交}\mathrm{集}}}{{\mathrm{B}\mathrm{O}\mathrm{X}}_{\mathrm{并}\mathrm{集}}} \text{,} $$ (7)
    $$ \mathrm{D}\mathrm{i}\mathrm{s}\mathrm{t}\mathrm{a}\mathrm{n}\mathrm{c}\mathrm{e}=\dfrac{|A{x_0}+B{y_0}+C|}{\sqrt{({A}^{2}+{B}^{2})}} \text{,} $$ (8)
    $$ \mathrm{C}\mathrm{e}\mathrm{n}\mathrm{t}\mathrm{e}\mathrm{r}\mathrm{I}\mathrm{O}\mathrm{U}\_\mathrm{L}\mathrm{o}\mathrm{s}\mathrm{s}= 1-(\mathrm{I}\mathrm{O}\mathrm{U}-\mathrm{D}\mathrm{i}\mathrm{s}\mathrm{t}\mathrm{a}\mathrm{n}\mathrm{c}\mathrm{e}) \text{。} $$ (9)

    由于定位损失函数更换为CenterIOU损失函数,较原GIOU函数让网络学习到预测框和目标框的相对位置,网络在反向传播时能更好地调优回归参数。

    为获取实际状态下病害检测所需时间,对比网络改进前后的效果,以训练好的模型为基础,开发手机端病害识别系统。利用IOS官方的Core ML Tools模型转换库将训练好的pth格式模型转换为IOS平台使用的格式。前端界面提供图库、识别、保存3个按钮,对应调用相册上传识别图片、快速识别作物病害、保存识别记录3种功能。编译后生成安装文件部署在iPhone 8手机,显示如图8所示。

    图  8  病害检测系统界面
    Figure  8.  Interface of disease detection system

    对于农作物叶片病害检测,需要网络具有高精确率和实时性。本文采用 COCO mAP@[0.50:0.05:0.95][15-16] 作为农作物病害检测评价指标来综合评估模型对苹果、玉米、苹果、玉米、葡萄、马铃薯、草莓、番茄6类农作物的27种病害和无叶片背景图像的检测性能,即分别在IOU阈值设置以0.05的步距从0.50开始遍历至0.95,取这些阈值条件下评估所得到的平均准确率(Average precision,AP)及其均值(Mean average precision,mAP)。mAP与准确率(Precision,P)、召回率(Recall, R)有关,其计算公式如下所示:

    $$ {P}=\dfrac{\mathrm{T}\mathrm{P}}{\mathrm{T}\mathrm{P}+\mathrm{A}\mathrm{P}}\times 100{\text{%}} \text{,} $$ (10)
    $$ {R}=\dfrac{\mathrm{T}\mathrm{P}}{\mathrm{T}\mathrm{P}+\mathrm{F}\mathrm{N}}\times 100{\text{%}} \text{,} $$ (11)
    $$ \mathrm{A}\mathrm{P}={\displaystyle\int }_{0}^{1}\mathrm{P}\left(R\right)\mathrm{d}R \text{,} $$ (12)
    $$ \mathrm{m}\mathrm{A}\mathrm{P}=\dfrac{1}{{M}}{\displaystyle\sum} _{k=1}^{\mathrm{M}}\mathrm{A}\mathrm{P}\left(k\right)\times 100\mathrm{{\text{%}}} \text{,} $$ (13)

    式中,TP表示被模型预测正确的正样本数,FP为被模型预测错误的正样本数,FN为被模型预测错误的负样本数,M为类别总数,AP(k)为第k类平均准确率。

    将训练好的模型部署移动端系统后,检测单张样本图像所消耗的平均时间(包含程序加载图片所耗时间和检测时间)作为实时性指标,单位为 s。参数量指标指深度网络模型中可训练参数的个数,以百万为单位,表示为1×106

    使用YOLOX-Nano模型为基础网络,分别引入CBAM卷积注意力机制、Mixup数据增强策略、使用Focal Loss替换原分类损失函数BCELoss,本文设计CenterIOU损失函数替换原回归损失函数GIOU作为改进网络。

    YOLOX-Nano对农作物病害叶片检测与识别的精度较低,单张样本图片检测精准率仅为97.97%,其原因是6种农作物27类病害的叶片病态特征表现比较相似。通道注意力可对特征图的不同通道赋予不同特征,选择性增大不同病害斑点特征通道的权重值;空间注意力对同一特征图不同位置特征点给予不同权重,区别对待特征图内部像素点。

    引入CBAM卷积注意力机制模块后,虽然模型参数量由0.92×106增加到0.96×106、单张样本推断时间增加18 ms,但平均精确率由97.97%提升到98.89%,提升了0.92个百分点,在很大程度上提升了对两者叶片病害识别的准确率。两者结合可总体提升识别准确率,表明了CBAM模块区别相似物体的有效性。

    利用Grad-CAM[17]工具可视化进入FPN层前第一个CBAM卷积注意力机制[18],明确对特征提取的效果。由于玉米病害的样本特征更为接近,可视化效果更佳,本文以玉米灰斑病样本、玉米锈病样本、玉米大斑病样本和健康玉米样本为例进行展示,如图9所示,引入卷积注意力机制前,网络对样本特征提取较随机,对叶片病害区域特征点关注不够。引入CBAM卷积注意力机制后,网络在前向传播过程中,重要特征通道逐渐占有更大的比重,能让网络学习到需重点关注的部分,使改进的YOLOX-Nano模型可以更高效地提取图像难区分的特征信息。

    图  9  Grad-CAM可视化添加CBAM前后玉米叶片特征提取的热力图
    图片的颜色由黑到蓝再到黄直到红,依次表示网络对某区域特征关注逐渐增多
    Figure  9.  Grad-CAM visual characteristic extraction thermal diagram of corn leaf before and after adding CBAM
    The color of the image changes from black to blue, then yellow to red, indicating that the network gradually pays more attention to the characteristics of a certain area

    使用Mixup数据增强策略训练的网络模型单张样本图片检测时间较改进前平均减少了4 ms,其检测准确率却增加了1.02个百分点,一方面减少了检测推断所用时间,另一方面可以提高对农作物病害检测精度,可见引入Mixup数据增强策略扩充数据集可以提升农作物病害检测的效率,模型可以训练更充分,减少推断时间。

    由式(6)可知,调节因子γ可以控制样本对损失的贡献,从而影响到病害分类的性能。试验中通过调节γ值来检测病害分类的效果,试验结果表明,γ取0时mAP为97.97%,γ取1时mAP为98.69%,γ取2时mAP为99.02%,γ取3时mAP为98.42%,γ取4时mAP为97.32%,mAP随γ值的增大呈现先升高后降低的趋势,当γ设置为2时,mAP最大,检测效果最好。

    图10a10b是选取γ为2时Focal Loss损失函数和原BCELoss损失函数在训练和测试阶段的分类损失的变化曲线,可以发现,改进后的YOLOX-Nano分类损失总体趋势都是低于YOLOX-Nano的。

    图  10  训练阶段(a)和测试阶段(b)分类损失变化曲线图
    Figure  10.  Classification loss curve during training(a) and test(b) phases

    改进后的YOLOX-Nano网络训练到100个迭代结束,训练和测试阶段均呈现出持续下降的趋势,且还存在下降空间,说明网络仍在学习,分类精度仍有上升的可能性,而曲线逐渐趋于稳定,说明网络过拟合已经没有学习必要。

    不同的回归框损失函数对检测精度带来不同程度的影响。试验通过比较GIOU和本文设计的CenterIOU损失函数的性能,来确定最适合的回归损失函数,最终试验表明使用GIOU的YOLOX-Nano在mAP50、mAP50:95的检测精度分别为98.23%和97.97%,而使用CenterIOU损失函数的YOLOX-Nano在mAP50、mAP50:95的检测精度分别为98.53%和98.20%,较改进前分别提高了0.30和1.23个百分点,本文设计的CenterIOU损失函数与YOLOX-Nano网络组合后,在数据集上对不同阈值表现得更好,能够明显地提高农作物病害的检测精度。

    改进后的YOLOX-Nano模型参数量虽然增加了0.06×106,单张图片检测时间增加了11 ms,但是改进后的YOLOX-Nano网络对所有农作物病害的平均检测精度达到了99.56%,较改进前精度(97.97%)提升了1.59个百分点。改进后的模型对苹果松锈病、葡萄黑腐病、马铃薯早疫病、马铃薯晚疫病、番茄叶霉病、番茄斑枯病和番茄黄叶曲叶病的识别精度都达到了100%,对玉米灰斑病识别精度(98.66%)稍低,主要原因是玉米大斑病与灰斑病病害症状相似,不易判断,导致出现错分类现象。

    同时,无叶片背景图这一类别存在背景中包含叶片的情况,模型能检测出叶片的存在,进而导致2个模型对无叶片背景图识别精度都表现较低,这也证明了本模型没有出现过拟合的现象,改进的网络有很好的农作物病害检测与识别能力。

    图11是改进前后网络训练平均精度和训练迭代的效果图,可见在训练10个迭代时,改进的YOLOX-Nano就能比改进前产生整体较高的mAP,预训练权重很快适应改进网络,趋于平稳,直到100个迭代训练结束,最终改进后的YOLOX-Nano精度更高。

    图  11  模型改进前后网络精度随迭代次数变化曲线图
    Figure  11.  Curve of network accuracy changing with number of iterations before and after model improvement

    本文还与现阶段较成熟的轻量型模型ResNet-18、MobileNet-v2、YOLOv4-Tiny、YOLOX-Tiny[19-24]进行对比,试验结果如表1所示,改进后的YOLOX-Nano检测准确率较其他轻量型网络高出1.0~2.2个百分点。在真机测试的检测时间对比中,本文提出方法的优势更加明显,单张样本图片检测时间仅需0.187 s,较其他各类模型节省大约20~240 ms的计算时间。不仅如此,在参数量方面,本文方法较ResNet-18减少了91.3%、较YOLOv4-Tiny减少了83.9%、较YOLOX-Tiny减少了80.7%、较MobileNet-v2减少了72.0%,更小的参数量意味着更少的部署成本,因此改进后模型的综合性能非常优秀,对于农业应用具有明显的优势。

    表  1  与主流轻量型网络性能对比
    Table  1.  Performance of the model versus mainstream lightweight network
    模型
    Model
    参数量(×106)
    No. of parameters
    平均精确率/%
    Mean average precision
    单张图片检测时间/s
    Detection time of single image
    YOLOX-Nano 0.92 97.97 0.176
    ResNet-18 11.24 98.60 0.429
    MobileNet-v2 3.40 97.33 0.189
    YOLOv4-Tiny 6.06 97.42 0.382
    YOLOX-Tiny 5.06 98.58 0.286
    改进YOLOX-Nano
    Improved YOLOX-Nano
    0.98 99.56 0.187
    下载: 导出CSV 
    | 显示表格

    为验证本网络移动端部署可行性和所提出网络的实际分类效果,从试验数据集中,随机取出50张图片,录入手机端系统进行识别,手机端检测部分结果如图12a、12b所示。经检测,所有图片病害全部识别准确,置信度保持在99.4%以上,最高置信度表现在苹果、葡萄和马铃薯3种作物的病害检测,均保持在99.8%以上,实际运行最多耗时0.324 s,最少耗时0.169 s,平均耗时为0.187 s。此外,从苹果、玉米、葡萄、马铃薯、草莓、番茄6类农作物病害样本中各筛选5张田间拍摄照片进行检测,部分检测结果如图12c、12d所示,所有病害依然识别准确,仍保持着99%以上的置信度,识别速度并未减弱,也无漏检错检的情况。可见本改进模型的检测效果和运行速度都保持着较高的水平,同时达到了精度和速度的平衡。

    图  12  病害检测结果展示
    Figure  12.  Display of disease detection results

    相对于人工诊断农作物病害,利用计算机视觉方法检测农作物病害具有成本低、准确率高、时延短的优点。本文通过在YOLOX-Nano分类模型中引入空间和通道注意力机制和Mixup数据增强策略,将分类损失函数更换为Focal Loss,提出一种改进的YOLOX-Nano优化网络,利用此网络对农作物病害数据集进行训练与测试,得出以下结论。

    1) 改进的YOLOX-Nano网络对PlantVillige数据集中6种农作物27种病害的检测总体mAP达到了99.56%,较YOLO-Nano基础网络提升了1.59个百分点,单张图片平均检测时间为0.187 s,仅仅较基础网络增加11 ms的推断延迟,达到了精度与速度的平衡,而且网络参数量为0.98×106,利于植保无人机设备和移动端设备的部署。

    2) 卷积注意力机制可以对不同通道的特征图进行权重的重新分配,加强网络结构深层信息的提取,对于很多细粒度特征的提取具有很好的效果,此外Mixup数据增强策略对一些小型网络也能起到促进作用;Focal Loss能控制容易分类和困难分类样本的权重,通过减少易分类样本的权重,在训练时更专注于难分类的样本,让容易样本对损失贡献小,难样本对损失贡献大。

    3) 本文提出的CenterIOU损失函数对本数据集训练具有促进作用,也可以推广应用到其他数据集。

  • 图  1   不同磷浓度对防风土壤理化性质的影响

    各图中,相同处理时间柱子上方的不同小写字母表示不同磷浓度间差异显著(P<0.05,Duncan’s法)

    Figure  1.   Effects of different phosphorus concentrations on physicochemical properties of soil in which Saposhnikovia divaricata was grown

    In each figure, different lowercase letters on the columns of the same treatment time indicate significant differences among different phosphorus concentrations (P<0.05, Duncan’s method)

    图  2   不同磷浓度对防风土壤酸性磷酸酶活性的影响

    相同处理时间柱子上方的不同小写字母表示不同磷浓度间差异显著(P<0.05,Duncan’s法)

    Figure  2.   Effects of different phosphorus concentrations on acid phosphatase activity of soil in which Saposhnikovia divaricata was grown

    Different lowercase letters on the columns of the same treatment time indicate significant differences among different phosphorus concentrations (P<0.05, Duncan’s method)

    图  3   不同磷浓度对防风生长的影响

    各图中,相同处理时间柱子上方的不同小写字母表示不同磷浓度间差异显著(P<0.05,Duncan’s法)

    Figure  3.   Effects of different phosphorus concentrations on growth of Saposhnikovia divaricata

    In each figure, different lowercase letters on the columns of the same treatment time indicate significant differences among different phosphorus concentrations (P<0.05, Duncan’s method)

    图  4   不同磷浓度下防风4种色原酮总含量变化

    相同处理时间柱子上方的不同小写字母表示不同磷浓度间差异显著(P<0.05,Duncan’s法)

    Figure  4.   Changes in total content of four chromones of Saposhnikovia divaricata in different phosphorus concentrations

    Different lowercase letters on the columns of the same treatment indicate significant differences among different phosphorus concentrations (P<0.05, Duncan’s method)

    图  5   不同磷浓度下的土壤理化性质、防风生长指标与4种色原酮总含量的相关关系

    1:根长,2:根粗,3:根鲜生物量,4:根干生物量,5:pH,6:电导率,7:有机质含量,8:碱解氮含量,9:有效磷含量,10:速效钾含量,11:酸性磷酸酶活性,12:4种色原酮总含量;“*”和“**”分别表示2个指标间显著(P<0.05)和极显著(P<0.01)相关(Pearson检验)

    Figure  5.   Correlation between soil physicochemical property, growth index with total content of four chromones of Saposhnikovia divaricataunder different phosphorus concentrations

    1: Root length, 2: Root diameter, 3: Root fresh biomass, 4: Root dry biomass, 5: pH, 6: Conductivity, 7: Organic matter content, 8: Alkali-hydrolyzale nitrogen content, 9: Available phosphorus content, 10: Available potassium content, 11: Acid phosphatase activity, 12: Total content of four chromones; “*” and “**” indicate significant correlations at 0.05 and 0.01 levels, respectively(Pearson test)

    表  1   不同磷浓度下的土壤因子、防风生长指标与4种色原酮总含量的灰色关联度分析

    Table  1   Grey correlation analysis between soil factors, growth index with total content of four chromones of Saposhnikovia divaricataunder different phosphorus concentrations

    因子 Factor 关联系数 Correlation coefficient
    1.0 mmol·L−1 0.1 mmol·L−1 0 mmol·L−1
    根长Root length 0.496 0.447 0.500
    根粗 Root diameter 0.386 0.447 0.444
    根鲜生物量 Root fresh biomass 0.406 0.423 0.396
    根干生物量 Root dry biomass 0.413 0.351 0.423
    pH 0.445 0.448 0.524
    电导率 Conductivity 0.334 0.332 0.320
    有机质含量 Organic matter content 0.295 0.334 0.486
    碱解氮含量 Alkali-hydrolyzale nitrogen content 0.274 0.451 0.442
    速效磷含量 Available phosphorus content 0.323 0.297 0.292
    速效钾含量 Available potassium content 0.270 0.355 0.328
    酸性磷酸酶活性 Acid phosphatase activity 0.394 0.459 0.460
    下载: 导出CSV

    表  2   不同磷浓度下的土壤理化性质、防风生长指标与4种色原酮总含量的通径系数1)

    Table  2   Path coefficient between soil factors, growth index with total content of four chromones of Saposhnikovia divaricata under different phosphorus concentrations

    c(P)/
    (mmol·L−1)
    变量
    Variate
    直接系数
    Direct coefficient
    间接系数 Indirect coefficient
    x1 x2 x3 x4 x5 x6 x7 x8 x9 x10 x11
    1.0 x1 −0.085 −0.012 0.037 0.003 0.082 −0.079 0.001 0.367 0.270 −0.017 −0.090
    x2 0.087 0.012 0.111 −0.014 −0.024 0.251 −0.023 −0.235 −0.116 −0.271 −0.035
    x3 0.173 −0.018 0.056 −0.018 0.012 −0.002 0.025 −0.002 0.029 −0.091 0.141
    x4 −0.041 0.005 0.029 0.076 0.010 0.210 0.007 −0.209 −0.044 −0.115 0.058
    x5 −0.174 0.040 0.012 −0.012 0.002 −0.114 0.049 −0.317 −0.391 0.230 0.249
    x6 0.585 0.012 0.037 −0.001 −0.015 0.034 −0.057 −0.102 0.031 −0.308 −0.313
    x7 −0.099 0.000 0.020 −0.044 0.003 0.086 0.333 0.052 0.143 −0.306 −0.397
    x8 −0.595 0.053 0.035 0.001 −0.014 −0.093 0.100 0.009 −0.365 0.052 0.099
    x9 −0.489 0.047 0.021 −0.010 −0.004 −0.139 −0.037 0.029 −0.444 0.178 0.237
    x10 −0.497 −0.003 0.048 0.032 −0.009 0.080 0.362 −0.061 0.062 0.175 −0.202
    x11 −0.523 −0.015 0.006 −0.047 0.005 0.083 0.350 −0.075 0.112 0.222 −0.192
    0.1 x1 −0.229 0.032 −0.081 0.007 −0.226 0.150 0.002 0.000 0.004 0.057 0.007
    x2 −0.103 0.071 0.140 −0.042 0.256 −0.195 −0.005 −0.027 −0.040 −0.073 −0.008
    x3 0.233 0.079 −0.062 −0.048 0.203 −0.144 −0.002 −0.006 −0.032 −0.043 −0.017
    x4 −0.061 0.027 −0.070 0.181 0.149 −0.144 −0.002 −0.013 −0.037 −0.059 −0.013
    x5 −0.505 −0.102 0.052 −0.094 0.018 0.390 0.003 0.005 −0.030 0.132 0.018
    x6 −0.534 0.064 −0.037 0.063 −0.017 0.369 −0.005 0.012 0.056 −0.162 −0.032
    x7 −0.008 0.042 −0.063 0.044 −0.016 0.150 −0.315 −0.008 −0.006 −0.091 −0.020
    x8 −0.060 0.000 −0.046 0.023 −0.013 0.041 0.106 −0.001 −0.059 0.021 0.013
    x9 −0.113 0.007 −0.037 0.067 −0.020 −0.134 0.266 0.000 −0.032 0.081 0.017
    x10 −0.171 0.077 −0.044 0.059 −0.021 0.391 −0.507 −0.004 0.008 0.053 −0.025
    x11 −0.045 0.037 −0.018 0.087 −0.018 0.206 −0.380 −0.004 0.017 0.042 −0.093
    0 x1 0.157 −0.509 0.277 −0.013 0.647 −0.151 0.003 0.051 −0.544 −0.043 0.036
    x2 −1.214 0.066 0.298 −0.023 1.011 −0.237 0.008 0.050 −0.177 −0.077 0.059
    x3 0.616 0.071 −0.588 −0.030 0.339 −0.136 −0.009 −0.048 −0.146 −0.022 −0.122
    x4 −0.038 0.054 −0.710 0.487 0.397 −0.104 −0.008 −0.007 −0.026 −0.025 −0.105
    x5 −1.234 −0.082 0.994 −0.169 0.012 0.346 −0.013 −0.110 0.412 0.112 −0.048
    x6 −0.508 0.047 −0.567 0.165 −0.008 0.840 0.013 −0.250 0.540 −0.166 −0.175
    x7 0.026 0.020 −0.366 −0.225 0.012 0.642 −0.260 −0.016 0.131 −0.088 0.068
    x8 0.914 0.009 −0.066 −0.033 0.000 0.148 0.139 −0.001 −1.225 0.076 0.160
    x9 −1.664 0.051 −0.129 0.054 −0.001 0.306 0.165 −0.002 0.673 0.078 0.200
    x10 −0.192 0.035 −0.484 0.070 −0.005 0.717 −0.438 0.012 −0.360 0.673 −0.191
    x11 −0.354 −0.016 0.203 0.213 −0.011 −0.168 −0.250 −0.005 −0.413 0.942 −0.104
     1) x1:根长;x2:根粗;x3:根鲜生物量;x4:根干生物量;x5:pH;x6:电导率;x7:有机质含量;x8:碱解氮含量;x9:有效磷含量;x10:速效钾含量;x11:酸性磷酸酶活性
     1) x1: Root length; x2: Root diameter; x3: Root fresh biomass; x4: Root dry biomass; x5: pH; x6: Conductivity; x7: Organic matter content; x8: Alkali-hydrolyzale nitrogen content; x9: Available phosphorus content; x10: Available potassium content; x11: Acid phosphatase activity
    下载: 导出CSV
  • [1] 陈敏旺. 吉林省农田耕层土壤速效氮磷钾养分的时空变化特征[D]. 长春: 吉林农业大学, 2018.
    [2] 张智猛, 戴良香, 张电学, 等. 施磷对土壤有效磷含量、吸磷特性及小麦产量的影响[J]. 河北农业技术师范学院学报, 1999, 1(1): 12-16.
    [3] 裴洲洋, 陈油鸿, 江春, 等. 不同施肥结构对土壤理化性质及烤烟氮、磷、钾含量的影响[J]. 安徽农学通报, 2016, 22(21): 35-37. doi: 10.3969/j.issn.1007-7731.2016.21.014
    [4] 张烨. 低磷胁迫诱导植物酸性磷酸酶的分子调控机制[D]. 北京: 清华大学, 2014.
    [5] 国家药典委员会. 中华人民共和国药典: Ⅰ部[M]. 北京: 中国医药科技出版社, 2020: 156.
    [6] 韩忠明, 王云贺, 胥苗苗, 等. 干旱胁迫对防风生理特性及品质的影响[J]. 西北农林科技大学学报(自然科学版), 2017, 45(11): 100-106.
    [7] 孙晖, 曹玲, 王喜军. 防风药效物质基础季节积累规律的研究[J]. 中国中医药科技, 2003, 10(6): 355-356. doi: 10.3969/j.issn.1005-7072.2003.06.019
    [8] 韩忠明, 胥苗苗, 王云贺, 等. 干旱胁迫对防风叶片保护酶活性、渗透调节物质含量及药材品质的影响[J]. 华南农业大学学报, 2016, 37(6): 91-97. doi: 10.7671/j.issn.1001-411X.2016.06.014
    [9] 胥苗苗. 水分胁迫对防风质量、生理生态特性及关键酶活性影响[D]. 长春: 吉林农业大学, 2016.
    [10] 黄艳茹, 任晓蕾, 霍金海, 等. 影响防风质量因素的研究概况[J]. 中医药信息, 2017, 34(6): 134-138. doi: 10.3969/j.issn.1002-2406.2017.06.038
    [11] 郭旭, 田新, 郝佳, 等. 不同产地防风色原酮含量与土壤因子的关系[J]. 华南农业大学学报, 2020, 41(4): 30-37. doi: 10.7671/j.issn.1001-411X.201912004
    [12] 邵畅畅, 韦小丽, 周紫晶, 等. 供磷水平对棕榈幼苗生长及根系形态的影响[J]. 东北林业大学学报, 2021, 49(1): 12-15.
    [13] 张文献, 李增强, 胡亚丽, 等. 不同浓度磷胁迫对大豆幼苗生长及根系DNA甲基化水平的影响[J]. 中国农业大学学报, 2020, 25(12): 9-18. doi: 10.11841/j.issn.1007-4333.2020.12.02
    [14] 鲍士旦. 土壤农化分析[M]. 3版. 北京: 中国农业出版社, 2000.
    [15] 关松荫. 土壤酶及其研究法[M]. 北京: 农业出版社, 1986.
    [16] 韩忠明, 王云贺, 李轶雯, 等. 防风色原酮微波辅助提取及含量动态变化研究[J]. 中药材, 2011, 34(3): 465-468.
    [17] 梁霞, 刘爱琴, 马祥庆, 等. 磷胁迫对不同杉木无性系酸性磷酸酶活性的影响[J]. 植物生态学报, 2005, 29(1): 54-59. doi: 10.3321/j.issn:1005-264X.2005.01.008
    [18] 唐荣莉, 王春萍, 王红娟, 等. 低磷胁迫对辣椒苗期生长和生理特性的影响[J]. 西南农业学报, 2020, 33(9): 1933-1942.
    [19] 王斐. 有机结合态磷对土壤微生物、玉米生长及磷素利用率的影响[D]. 哈尔滨: 东北农业大学, 2013.
    [20] 余健. 磷胁迫下林木分泌的有机酸及对土壤磷的活化[D]. 南京: 南京林业大学, 2005.
    [21] 王庆仁, 李继云, 李振声. 植物高效利用土壤难溶态磷研究动态及展望[J]. 植物营养与肥料学报, 1998, 4(2): 107-116. doi: 10.3321/j.issn:1008-505X.1998.02.003
    [22] 刘晓霞. 苹果园土壤磷库特征及土壤磷含量和磷源对磷吸收的影响[D]. 泰安: 山东农业大学, 2018.
    [23] 谢安强, 洪伟, 吴承祯, 等. 内生真菌对低磷胁迫下尾巨桉生理及土壤特性的影响[J]. 西南林业大学学报, 2013, 33(3): 1-7. doi: 10.3969/j.issn.2095-1914.2013.03.001
    [24]

    CHEN H, LI D, XIAO K, et al. Soil microbial processes and resource limitation in karst and non‐karst forests[J]. Functional Ecology, 2018, 32(5): 1400-1409. doi: 10.1111/1365-2435.13069

    [25] 李荣坦. 低磷胁迫对番茄根系生长及根际土壤微生物多样性的影响[D]. 南宁: 广西大学, 2016.
    [26] 徐静, 张锡洲, 李廷轩, 等. 野生大麦对土壤磷吸收及其酸性磷酸酶活性的基因型差异[J]. 草业学报, 2015, 24(1): 88-98. doi: 10.11686/cyxb20150112
    [27] 郭轶敏, 文亦芾, 史亮涛, 等. 土壤磷水平对柱花草株高、根长及生物量的影响[J]. 草业与畜牧, 2014(1): 1-4.
    [28]

    LIU Y, MI G, CHEN F, et al. Rhizosphere effect and root growth of two maize ( Zea mays L. ) genotypes with contrasting P efficiency at low P availability[J]. Plant Science, 2004, 167(2): 217-223. doi: 10.1016/j.plantsci.2004.02.026

    [29] 王鑫, 郭平毅, 原向阳, 等. 2, 4−D丁酯对罂粟(Papaver somniferum L.)保护酶活性及脂质过氧化作用的影响[J]. 生态学报, 2008, 28(3): 1098-1103. doi: 10.3321/j.issn:1000-0933.2008.03.024
图(5)  /  表(2)
计量
  • 文章访问数:  690
  • HTML全文浏览量:  18
  • PDF下载量:  556
  • 被引次数: 0
出版历程
  • 收稿日期:  2021-04-22
  • 网络出版日期:  2023-05-17
  • 刊出日期:  2022-05-09

目录

/

返回文章
返回