Effects of rice-amaranth intercropping on rice yield，rhizosphere soil available nutrient and As accumulation
-
摘要:目的
探究水稻与苋菜间作对水稻产量、根区速效养分及水稻各部位吸收累积As的影响。
方法通过小区试验研究水稻单作(对照)、水稻与圆叶白苋菜间作、水稻与籽粒苋间作、水稻与红苋菜间作条件下水稻产量、根际土壤速效养分与有机质含量、As累积量的变化。
结果4个处理水稻收获期地上部、地下部、谷壳和糙米中As含量均表现为：水稻与圆叶白苋菜间作>水稻与籽粒苋间作>水稻与红苋菜间作>水稻单作。与对照相比，水稻与圆叶白苋菜间作显著提高了水稻地上部、地下部、谷壳和糙米As含量，增幅分别为88.50%、65.38%、68.25%和45.45%；对土壤As的提取量为5.44 mg·m−2，显著高于对照，为对照的1.67倍；水稻单株生物量显著高于对照、与红苋菜间作和与籽粒苋间作，地上部生物量比对照提高63.79%。与水稻单作相比，间作处理水稻单位面积产量比对照显著下降，间作中单位面积产量最高的为与圆叶白苋菜间作，为对照的59.15%；间作处理水稻根际土壤速效养分、有机质含量均有不同程度的增加。
结论水稻与苋菜间作，水稻单位面积产量降低，地上部、地下部、谷壳和糙米中As含量以及根区土壤速效养分与有机质含量显著提高。水稻与圆叶白苋菜间作能从As污染土壤中提取最多的As，生产出的苋菜符合国家安全标准，为最优方案。
Abstract:ObjectiveTo explore the effects of rice-amaranth intercropping on rice yield, rhizosphere soil available nutrient and As accumulation in various parts of rice.
MethodThrough field plot experiment, four treatments including rice monocropping (control), rice-white amaranth intercropping, rice-grain amaranth intercropping and rice-red amaranth intercropping were setted to study the changes of rice yield and As accumulation, rhizosphere soil available nutrient and organic matter contents.
ResultAs contents of shoot, root, husk and brown rice during harvest time in four treatments all behaved as rice-white amaranth intercropping > rice-grain amaranth intercropping > rice-red amaranth intercropping > rice monocropping. Compared with the control, rice-white amaranth intercropping significantly increased As contents of shoot, root, husk and brown rice by 88.50%, 65.38%, 68.25%, 45.45% respectively. As extract amount from soil in rice-white amaranth intercropping treatment reached 5.44 mg·m −2, which was 1.67 times significantly higher than control. The biomass per rice plant in rice-white amaranth intercropping treatment was significantly higher than control, rice-red amaranth intercropping and rice-grain amaranth intercropping, and 63.79% higher than control. Compared with control, intercropping treatments significantly reduced rice yield per unit area. Rice-white amaranth intercropping had the highest yield per unit area, which was 59.15% of control. The available nutrient and organic matter contents in rice rhizosphere soil increased in different degrees compared with monoculture.
ConclusionRice-amaranth intercropping reduces rice yield per unit area, and significantly increases As contents in shoot, root, husk and brown rice, and available nutrient and organic matter contents in rhizosphere soil. Rice-white amaranth intercropping is the optimal scheme, as it extract the most As in As contaminated soil, and the grown amaranth meets national safety standards.
-
Keywords:
- rice /
- amaranth /
- intercropping /
- As accumulation /
- As contaminated soil
-
芒果Mangifera indica作为著名的热带亚热带水果，以其色香味俱佳而享有“热带果王”的美誉，具有极高的营养价值和经济价值，在热带及南亚热带地区有广泛栽培[1]。我国是芒果原生产地之一，种植历史悠久，在海南、广东、广西、福建、云南、四川等省份建立的规模化种植的芒果园，成为重要的地区特色农业经济产业。
由于长期施用过量的化学肥料和农药，土壤环境不断地发生恶化，果类病毒产生新的抗药性，芒果出现新的感染状况并产生诸多病症，如海绵组织病、心腐病、软鼻子病、空心病等生理性病害，对芒果的品质产生了严重的影响。在部分芒果产区，由此类病害造成的芒果损失率甚至达到了30%以上。患有此类病害的芒果，表面与没有患病的芒果相比并无差异，凭肉眼无法观测出来，经验丰富的果农也无法准确对此类芒果做到挑拣和分类，这也就对芒果果实的品质检测和分级技术提出了新的应用方向与挑战[2-5]。
现代农业无损检测技术正朝着智能化方向发展，结合计算机图形图像分析与深度学习的技术正成为水果品质检测的新方向和研究热点[5-7]。研究采用可见/近红外光谱、核磁共振和X射线等检测技术，无损地获得果类果实的内部组织特征图像，然后结合图像处理和模式识别等相关技术，进行组织提取和病状类型分析识别，从而实现果实品质检测和分级[8-9]。
国内外众多学者在农产品果类检测方向开展了大量的研究。2018年余心杰[10]基于可见光和近红外光(Visible-near infrared, VIS-NIR)高光谱反射成像技术，研究开发了一种由堆栈自编码器(Stacked auto-encoder，SAE)和全连接神经网络组成的深度学习方法，用于预测收获的库尔勒香梨的硬度和可溶性固形物含量(Soluble solids content, SSC)，证明了深度学习方法与高光谱成像技术相结合，可用于库尔勒香梨的快速无损检测硬度和SSC。2019年Nasiri等[11]采用深卷积神经网络模型(Convolutional neural network, CNN)VGG-16在枣果图像数据集上进行训练和测试，来鉴别健康枣果和缺陷枣果并用来预测健康枣的成熟期，VGG-16模型的分类准确率达到了96.98%，证明了深卷积神经网络模型优于传统的基于图像特征的枣果图像分类方法。2019年Wan等[12]提出了一种改进的快速R-CNN深度学习网络模型，用于机器人视觉系统，实现了苹果、芒果和桔子等多类水果检测，达到了91%以上的检测正确率，该深度学习模型提高了图像处理速度，也更适合于实际场景。2020年Osako等[13]通过对预训练的VGG16网络进行微调，构建了荔枝品种识别模型，该模型对荔枝水果图像进行品种识别，实现了98.33%的准确性，研究表明深度学习可用于区分荔枝品种。
迁移学习是一种机器学习的方法[14-17]，是解决训练数据不足这一基本问题而发展起来的重要方法，其基本原理是通过已学习的相关任务中知识转移来改进新任务的学习[17-20]。近年来深度学习技术在图像识别分析领域取得了长足的进步，但是必须要有海量的训练数据作为支撑，来理解数据的潜在模式，才能够实现网络模型的训练[20-22]。深度迁移学习将在原域中预先训练好的深度学习网络结构重新利用，将其用于目标领域的深度神经网络的一部分，因为原网络已经学习了图像丰富特征，有很好的泛化性，通过网络微调可以实现相应的目标领域学习，大大加快了学习速度[23-25]。
本文基于芒果果实计算机层析成像(Computed tomography，CT)设备所得的芒果CT序列图像，提出了传统图像处理技术和深度迁移技术相结合的芒果病状识别判断方法。工作内容包括：1)芒果果实CT序列图像内部坏损区域的提取方法研究；2)利用AlexNet和GoogLeNet[26-28]预训练网络开展芒果病状的识别研究。
1. 研究方法
1.1 芒果CT序列图像的获取
基于平板探测器的计算机层析成像系统(图1)采集芒果CT序列图像。该系统主要由X射线源分系统、平板探测器分系统、机械扫描分系统、图像处理分系统以及安全及辅助分系统组成[29]。该系统采用的平板探测器为VARIAN PaxScan2520V，其像元尺寸为0.127 mm，A/D转换器位数为14位。
![图 1 基于平板探测器的计算机层析成像(CT)系统]()
图 1 基于平板探测器的计算机层析成像(CT)系统Figure 1. Computer tomography (CT) system based on flat panel detector选择外表无破损、内部组织坏损的芒果类型。实验室环境要求常温，湿度小于82%，磁场干扰较小，所有芒果样品图像采集的过程中实验室条件、采集方法和参数设置必须一致，这样才能减少误差。设备参数：平板探测器的积分时间为200 ms，X光机的电压为45~50 kV，电流为1 mA，CT检测时根据检测对象大小设置热区参数，扫描方式为标准连续扫描。
1.2 芒果CT序列图像预处理
本文采集海绵组织病和空心病芒果16位DICOM格式CT序列图像(512×512)，海绵组织病序列图像650张，空心病序列图像379张，有效图像分别为624张和365张。海绵组织病和空心病芒果CT序列图像如图2所示。
![图 2 海绵组织病芒果和空心病芒果]()
图 2 海绵组织病芒果和空心病芒果Figure 2. Mango with spongy tissue disease and mango with hollow disease芒果CT图像的像素点的灰度值集中分布在0~35以内，芒果外区域集中在0~10。图3中的a和b分别是对海绵组织病芒果和空心病的CT序列图像作灰度值统计，验证了芒果有效像素值分布在灰度值较低范围，截取256灰度级的原始图像灰度。
![图 3 海绵组织病芒果和空心病芒果CT图像灰度值分布]()
图 3 海绵组织病芒果和空心病芒果CT图像灰度值分布Figure 3. Gray value distribution in CT images of mango with spongy tissue disease and mango with hollow disease灰度变换是实现图像对比度增强的重要手段，最终的目的是为了增强图像清晰度，让图像细节更加突出。在样本的处理中发现未坏损区域的灰度值接近于恒定在14~35，根据灰度的分布特点，结合实际的图像效果，本文采用了分段函数法来实现图像对比度增强。
灰度变换原理是把原图像的灰度函数
$f(x,y)$ ，经过1个变换函数T变换成1个新的图像函数$g(x,y)$ ，即$$ g(x,y) = T[f(x,y)], $$ (1) 通过变换，达到对比度增强的效果，在变换时对每个像素
$(x,y)$ 都经过了同样的处理，因此该方法又叫做图像点处理。分段函数法又称为分段变换，为了突出感兴趣的目标或灰度区间，抑制不感兴趣的目标或灰度区间，常采用分段线性变换法。芒果图像未坏损区域的灰度值接近于恒定在14~35，采用分段函数将未坏损区域的灰度值拉伸变大，坏损区域的灰度值保持不变，公式如下：
$$g(x,y) = \left\{ {\begin{array}{*{20}{l}}{f(x,y) \times 7},&{ 14 \leqslant f(x,y) \leqslant 35}\\{f(x,y) \times 1 },&{f(x,y) < 14}\end{array}} \right.{\text{。}}$$ (2) ![图 4 海绵组织病芒果图像采用分段函数法变换]()
图 4 海绵组织病芒果图像采用分段函数法变换Figure 4. Piecewise function to transform image of mango with spongy tissue disease![图 5 空心病芒果图像采用分段函数法变换]()
图 5 空心病芒果图像采用分段函数法变换Figure 5. Piecewise function to transform image of mango with hollow disease1.3 图像降噪
芒果CT序列图像在获取的过程中，因为受到温湿度、磁场干扰、信号传输过程中的损耗、振动噪声等试验条件的限制，降低了图像的质量，造成最终的成像效果失真，这样必然影响后续的图像分析与研究，采用图像滤波处理的方法来进行降噪处理。根据获取的芒果CT序列图像的成像特点，研究采用空间域滤波处理的方式。空间域滤波包括线性滤波和非线性滤波。最终采用了中值滤波和双边滤波相结合的方法，对图像实现降噪处理。中值滤波有去除脉冲噪声和椒盐噪声，保留图像的边缘细节的作用。其算法原理是基于滤波器所在图像区域中像素灰度值
${Z_k}$ 的排序，由排序结果决定的值R代替中心像素的值，公式如下：$$ R = {\rm{mid}}\left\{ {{Z_k}|k = 1,2,3, \cdots } \right\}{\text{。}} $$ (3) 双边滤波是高斯滤波的改进，高斯滤波中只有关于空间距离的权重，而双边滤波在此基础上添加了关于灰度距离的权重，可以消除图像噪声，实现图像平滑，同时保留边缘信息。实现方法：1个与空间距离相关的高斯函数与1个灰度距离相关的高斯函数相乘，得到最终的卷积模板，双边滤波需要每个中心点邻域的灰度信息来确定其系数，计算量增长速度为核大小的平方。空间距离指的是当前点与中心点的欧式距离。空间域高斯函数其数学形式为：
$$d(i,j,k,l) = {\rm{exp}}\Bigg( - \frac{{\left( {{{(i - k)}^2} + {{(j - l)}^2}} \right)}}{{2\sigma _{\rm{d}}^2}}\Bigg),$$ (4) 式中，
$(\left. {i,j} \right)$ 代表输出点，$(\left. {k,l} \right)$ 代表(多个)输入点；${\sigma _{\rm{d}}}$ 为空间域标准差。灰度距离：当前点灰度值与中心点灰度值之差的绝对值。值域高斯函数其数学形式为：$$r(i,j,k,l) = {\rm{exp}}\Bigg( - \frac{{\left. {\left\| {f{{(i,j)}^2} - f{{(k,l)}^2}} \right.} \right\|}}{{2\sigma _{\rm{r}}^2}}\Bigg),$$ (5) 式中，
$(\left. {i,j} \right)$ 代表输出点；$(\left. {k,l} \right)$ 代表(多个)输入点；${\sigma _{\rm{r}}}$ 为值域标准差。将上述2个公式相乘就可以得到依赖数据的双边滤波权重函数：$$\omega (i,j,k,l) = \exp \left( { - \frac{{{{(i - k)}^2} + {{(j - l)}^2}}}{{2\sigma _{\rm{d}}^2}} - \frac{{\parallel f{{(i,j)}^2} - f{{(k,l)}^2} \parallel }}{{2\sigma _{\rm{r}}^2}}} \right),$$ (6) 将得到的权重与邻域像素灰度值相乘并叠加，然后再除以权重和，最终获得该像素点的灰度值：
$$g(i,j) = \frac{{\displaystyle\sum\limits_{k,l} f (k,l)\omega (i,j,k,l)}}{{\displaystyle\sum\limits_{k,l} \omega (i,j,k,l)}}{\text{。}}$$ (7) 由于X射线在穿透物体时，光子会发生散射、折射和衍射等现象，在芒果CT图像区域会形成大量的椒盐噪声(如图6所示)。要求去除椒盐噪声，尽量保留坏损区域的边缘信息和实现图像平滑，确定先双边滤波，再进行中值滤波的滤波处理方案。选定中值滤波窗口尺寸为4×4，双边滤波参数：滤波半径=5，全局方差=8，局部方差=0.2，达到了较好的图像处理效果。
从整体和细节放大角度进行了对比分析，以海绵组织病芒果的第300层CT序列图像为例，进行了滤波效果展示。通过图7、图8分析，中值滤波结合双边滤波的图像滤波处理方案是可行的，能够去除椒盐噪声，实现图像平滑，坏损区域边缘信息得到了保留。文中采用的均值滤波器和高斯滤波器模板如图9所示。
![图 7 芒果CT图像多方案滤波效果对比]()
图 7 芒果CT图像多方案滤波效果对比Figure 7. Effect comparison of several filtering schemes in mango CT images![图 8 芒果CT图像多方案滤波效果细节放大对比]()
图 8 芒果CT图像多方案滤波效果细节放大对比Figure 8. Effect comparison of several filtering schemes with detail amplification in mango CT images1.4 图像分割
图像分割是为了提取芒果CT图像内部组织坏损区域，对芒果果实内部的各种病害进行特征分析。不同病状芒果以及同一芒果不同层级序列图像呈现出不同的差异特点，使用图像二值化、形态学填充和差影法相结合的图像分割算法。
因为芒果CT序列图像上的不同部分具有不同亮度，局部自适应阈值法可以在同一幅图像上的不同区域采用的是不同的阈值，亮度较高的图像区域的二值化阈值通常会较高，而亮度较低的图像区域的二值化阈值则会相应地变小，从而在亮度不同的情况下得到更好的结果。本文选择了局部自适应阈值法，进行图像的二值化处理。局部自适应阈值法步骤：
1)均值平滑图像(本文采用4×4模板)。
2)求平滑图像梯度图：
$$\begin{array}{*{20}{c}} {{G_x}(x,y)}{ = P(x + 1,y) - P(x - 1,y)} ,\\ {{G_y}(x,y)}{ = P(x,y + 1) - P(x,y - 1)} ,\\ {F(x,y)}{ = \sqrt {{G_x}{{(x,y)}^2} + {G_y}{{(x,y)}^2}} } , \end{array}$$ (8) ${G_x}(x,y)$ ：x方向梯度，${G_y}(x,y)$ ：y方向梯度，$P(x,y)$ 像素点灰度值，$F(x,y)$ ：图像梯度。3)平滑图像作Laplacian运算：
$$\begin{array}{*{20}{c}} {R(x,y) = P(x,y + 1) + P(x,y - 1) + P(x - 1,y) + } \\ {P(x + 1,y) - 4P(x,y)} , \end{array}$$ (9) 找到局部最大阈值的点，这些点的原始灰度值作为候选局部阈值，
$R(x,y)$ 表示Laplacian算子，查找经Laplacian运算后图像中灰度值为零的点，将这些点对应的原始图像灰度值作为候选局部阈值。4)新建与原图像同样大小的零矩阵图像，用原始图像候选点中的灰度值代替零矩阵图像中对应像素的灰度值。
5)多次迭代求阈值面：
$$\mathop {\lim }\limits_{n \to \infty } {P_n}(x,y) = {P_{n - 1}}(x,y) + \frac{{\beta \cdot {R_n}(x,y)}}{4},$$ (10) 当
$\;\beta = 0$ ，残差消失，$1 < \beta < 2$ ，收敛更快，本文中$\;\beta = 1$ ；因为实际中n值不可能取无限大，本文程序中取n值不大于1 500；当${P_n}(x,y)$ 与${P_{n - 1}}(x,y)$ 相减所得残差小于${10^{{\rm{ - }}3}}$ ，停止迭代。6)根据阈值面作二值分割：
$$g(x,y) = \left\{ {\begin{array}{*{20}{l}} {0,}&{f(x,y) < {P_n}(x,y)} \\ {255,}&{f(x,y) \geqslant {P_n}(x,y)} \end{array}} \right.{\text{。}}$$ (11) 本文对海绵组织病和空心病芒果所有有效层CT图像做了二值化处理。以图10中海绵组织病芒果二值化图像为例，局部自适应阈值法可以准确地将坏损区域和正常果肉区域分开，各相近层提取的坏损区域之间表现出明显的连续性。
![图 10 海绵组织病芒果二值化序列图像]()
图 10 海绵组织病芒果二值化序列图像a、b、c、d、e、f、g、h分别是芒果的第300、304、308、312、500、504、508、512层二值化图像Figure 10. Binary sequence images of mangoes with spongy tissue diseasea, b, c, d, e, f, g, h are binary images for the 300th, 304th, 308th, 312th, 500th, 504th, 508th and 512th layers of mango, respectively研究的病状芒果果实样品都是外部完好、内部出现坏损的类型，果实的边缘部分都是封闭的，所以本文中采用了基于边界封闭的种子填充形态学方法进行芒果果实图像内部区域填充[30]。种子填充算法定义：
$${X_k} = \left( {{X_{k - 1}} \oplus B} \right) \cdot {A^c},\quad k = 1,2,3, \cdots , $$ (12) 点
$ {X}_{0} $ 像素灰度值为1，$B$ 表示结构元，$ \oplus $ 表示形态学膨胀运算，如果${X_k} = {X_{k - 1}}$ ，则算法在迭代的第k步结束，用1填完所有孔洞。${X_k}$ 和A的并集包含被填充的集合和它的边界。利用Ac的交集将结果限制在感兴趣区域内，实现条件膨胀。本文选择种子填充算法，孔洞填充后效果如图11b、图12b所示。差影法多应用于同一背景的图片中搜索差异。图像f1被填充，其背景、位置在填充过程中没有发生改变，整个区域的值都为255，得到填充图像f2；未填充的图像f1保留着内部坏损信息(坏损区域值为0)，与f2做差影，得到图像f3。差影法公式：![图 11 空心病芒果图像差影法分割效果(第300层)]()
图 11 空心病芒果图像差影法分割效果(第300层)Figure 11. Image segmentation effect of mango with hollow disease using background subtraction (The 300th layer)![图 12 海绵组织病芒果图像差影法分割效果(第300层)]()
图 12 海绵组织病芒果图像差影法分割效果(第300层)Figure 12. Image segmentation effect of mango with spongy tissue disease using background subtraction (The 300th layer)$${\rm{f}}2 - {\rm{f}}1 = {\rm{f}}3{\text{。}}$$ (13) 空心病芒果序列图像的第300层和海绵组织病芒果序列图像的第300层的图像分割效果，如图11、图12所示。
2. 分类算法
2.1 分类模型
深度迁移学习模型，针对未做图像处理和已做图像处理的芒果图像数据开展训练和测试，分析分类识别效果。基于AlexNet和GoogLeNet深度学习网络模型，芒果序列图像数据集作为训练数据开展迁移学习，调整超参数，完成训练过程的网络微调，实现分类任务。在不同模型下，对比了未做图像处理和做了图像处理的芒果图像测试集在模型的分类结果。
AlexNet和GoogLeNet已基于超过一百万个图像进行训练，将图像分为1 000个对象类别，学习了丰富的特征表示。AlexNet是在LeNet的基础上加深了网络的结构，学习更丰富更高维的图像特征，其特点包括：更深的网络结构；使用层叠的卷积层和池化层来提取图像的特征；使用Dropout抑制过拟合；使用数据增强抑制过拟合；使用Relu替换之前的Sigmoid的作为激活函数。GoogLeNet相比与AlexNet网络，增加了网络层数和Inception模块增强卷积模块功能，在增加网络深度和宽度的同时减少了参数，其参数数量约为AlexNet的1/12，提升了计算性能。
2.2 模型算法流程
修改预先训练好的网络并使用迁移学习进行训练以执行新识别任务，微调预先训练好的网络比构造和训练新网络更快、更容易，使用更少数量的训练图像迅速迁移学习到新任务。迁移学习的算法流程，如图13所示。
最后3层替换为全连接层、softmax层和分类输出层，将层迁移到新分类任务。根据芒果数据指定新的全连接层的选项，将全连接层分类值设置为分类数目。对于迁移学习，要保留预训练网络的较浅层中的特征(迁移层权重)，通过增大全连接层的学习率因子，加快新的最终层中的学习速度，这种学习率设置组会加快新层中的学习速度；执行迁移学习时，所需的训练轮数相对较少，一轮训练是对整个训练数据集的一个完整训练周期。
网络评估的目的是为了提高模型在应用过程的准确率，可以继续迭代，重新研究训练方案、检查数据以及重新配置模型参数。评估分类模型指标有精准率(Precision)、召回率(Recall)、F1分值(F1-score)、准确率(Accuracy)和宏平均(Macro-average)。
$$ {\rm{Precision}} = {\rm{TP}}/({\rm{TP}} + {\rm{FP}}), $$ (13) $$ {\rm{Recall}} = {\rm{TP}}/({\rm{TP}} + {\rm{FN}}), $$ (14) $${{F}}1 {\text{-}} {\rm{score }} = \frac{{ {\rm{ }}2{\rm{Recall}} \cdot {\rm{Precision}}}}{{{\rm{Recall}} + {\rm{Precision}}}}, $$ (15) $$ {\rm{Accuracy}} = \frac{{{\rm{TP + TN}}}}{{{\rm{TP + TN + FP + FN}}}}, $$ (16) $$ {\rm{Macro {\text{-}} average}} = \frac{{\displaystyle\sum {{{F}}1 {\text{-}} {\rm{score}}} }}{n}, $$ (17) 式中，TP表示样本属于正类，且被正确地归为正类；FP表示样本不属于正类，但被错误地归为正类；FN表示样本属于正类，但被错误地归为其他类；TN表示样本属于其他类，且被正确地归为其他类，n表示类的个数。Precision衡量正样本的分类准确率，Recall表示分类正确的正样本占总的正样本的比例，F1-score是Precision和Recall的调和平均，Accuracy衡量所有样本被分类准确的比例。Macro-average则是综合了性能指标的Precision、Recall和F1-score的信息，常和Accuracy一起用来评价模型整体的分类性能。
3. 结果与分析
3.1 芒果图像数据集
芒果数据集分为训练集、验证集和测试集，打上标签。芒果图像数据集分为两大类，即未经图像处理数据集和经图像处理后的数据集。其中未经图像处理数据集包括3类：无病症芒果CT图片(2 168张)，海绵组织病芒果CT图片(1 924张)，空心病芒果CT图像(864张)。经图像处理后的数据集包括5类：无病症芒果图片(2 168张)、海绵组织病芒果图片(1 924张)、空心病芒果图片(864张)、海绵组织病坏损区域图片(1 920张)、空心病坏损区域图片(864张)。通过水平翻转、垂直翻转、顺时针旋转90°得到了更多的图像数据，图14和图15分别展示了未经图像处理数据集和经图像处理后的数据集其变换后的形式。
3.2 加载未经图像处理数据集
未经图像处理数据集分为训练集、验证集和测试集，随机取各类图片数目的60%作为训练集，20%作为验证集，20%作为测试集。训练集是作为模型拟合的数据样本，验证集用于调整模型的超参数和对分类模型的能力进行初步评估，测试集是用于评估分类模型的泛化能力。
深度迁移学习模型训练阶段监视模型的训练准确率、训练损失、验证准确率和验证损失，在训练过程中，通过调整超参数来优化分类模型，调整的超参数：迭代轮数(Epoch)、最小批(Mini batch)和学习速率(Learning rate)。迭代1轮是指训练算法完全通过整个训练集，训练迭代轮数越多，网络训练时间越长，准确率会随着每次迭代而提高。最小批是指在CPU上处理的训练数据集的子集，最小批越大，训练速度越快。学习速率是控制训练精度的主要参数，学习速率越低，训练结果准确率越高，但网络训练时间会更长。基于未处理的芒果图像数据集进行训练时，迭代轮数设为3，最小批值设置为64，不断调整学习速率来优化分类模型。基于AlexNet和GoogLeNet网络通过设置学习率不同值，开展训练，现分别取值为0.000 8、0.000 5和0.000 2，进行分类模型性能的评估。表1中展示了各个模型针对未处理芒果图像的分类性能参数。所有模型的迭代次数均为138。
表 1 深度迁移学习模型性能参数(未处理图像)Table 1. Performance parameters of deep transfer learning model (Unprocessed image)模型
Model学习速率
Learning rate图片类型
Image type精准率/%
Precision召回率/%
RecallF1分值/%
F1-score模型准确率/%
Accuracy of modelAlexNet 0.0008 无病症芒果 Healthy mango 87.58 92.63 90.03 84.88 海绵组织病芒果
Mango with spongy tissue disease94.90 77.40 85.26 空心病芒果 Mango with hollow disease 64.84 82.08 72.45 0.0005 无病症芒果 Healthy mango 96.02 94.47 95.24 90.22 海绵组织病芒果
Mango with spongy tissue disease93.71 85.19 89.25 空心病芒果 Mango with hollow disease 73.02 90.75 80.93 0.0002 无病症芒果 Healthy mango 97.22 96.77 97.00 94.15 海绵组织病芒果
Mango with spongy tissue disease96.70 91.43 93.99 空心病芒果 Mango with hollow disease 82.63 93.64 87.80 GoogLeNet 0.0008 无病症芒果 Healthy mango 95.55 94.01 94.77 90.73 海绵组织病芒果
Mango with spongy tissue disease95.40 86.23 90.59 空心病芒果 Mango with hollow disease 82.05 92.49 82.05 0.0005 无病症芒果 Healthy mango 98.12 93.61 97.21 95.36 海绵组织病芒果
Mango with spongy tissue disease97.32 94.29 95.78 空心病芒果 Mango with hollow disease 85.49 95..38 90.16 0.0002 无病症芒果 Healthy mango 99.77 99.08 99.42 98.79 海绵组织病芒果
Mango with spongy tissue disease98.96 98.96 98.96 空心病芒果 Mango with hollow disease 96.02 97.69 96.85 文中加载未经图像处理数据集，通过训练深度迁移学习模型，构建了多种分类算法模型。AlexNet模型在学习率0.0008、0.0005和0.0002下，Accuracy分别为84.88%、90.22%和94.15%，Macro-average分别为83.18%、88.47%和92.93%。GoogLeNet模型在学习率0.0008、0.0005和0.0002下，Accuracy分别为90.73%、95.36%和98.79%，Macro-average分别为89.14%、94.38%和98.41%。当学习率为0.0005和0.0002，深度迁移学习模型分类效果较优；学习率值变小，深度迁移学习模型的分类性能有较大提升；在同一学习率下，GoogLeNet网络的分类效果明显优于AlexNet网络。
3.3 加载图像处理数据集
与未处理图像的分类比例相同，随机取各类图片数目的60%作为训练集，20%作为验证集，20%作为测试集。图像处理数据集与未经处理数据集相比，增加了海绵组织病坏损区域图片和空心病坏损区域图片，为保证对比效果，其他各类图片数量保持不变。对处理的数据集进行模型训练时，将迭代轮数值固定为3，将最小批值设置为64，不断调整学习速率值来优化分类模型。基于AlexNet和GoogLeNet网络训练模型的学习率分别取值为0.0008、0.0005和0.0002，进行性能的评估。表2中展示了各个模型针对未处理芒果图像的分类性能参数。所有的模型的迭代次数均为216。
表 2 深度迁移学习模型性能参数(处理图像)Table 2. Deep transfer learning model performance parameters (Processed image)模型
Model学习速率
Learning rate图片类型
Image type精准率/%
Precision召回率/%
RecallF1分值/%
F1-score模型准确率/%
Accuracy of modelAlexNet 0.0008 无病症芒果 Healthy mango 100 96.48 98.21 95.39 海绵组织病坏损区域
Damaged area from spongy tissue disease96.16 97.92 97.03 空心病坏损区域
Damaged area from hollow disease86.46 95.95 90.96 海绵组织病芒果 Mango with spongy tissue disease 93.67 96.10 94.87 空心病芒果 Mango with hollow disease 96.71 84.97 90.46 0.0005 无病症芒果 Healthy mango 97.42 97.42 98.57 97.40 海绵组织病坏损区域
Damaged area from spongy tissue disease98.44 98.44 98.44 空心病坏损区域
Damaged area from hollow disease100 100 95.58 海绵组织病芒果 Mango with spongy tissue disease 98.18 98.18 96.92 空心病芒果 Mango with hollow disease 90.75 90.75 95.15 0.0002 无病症芒果 Healthy mango 100 100 100 100 海绵组织病坏损区域
Damaged area from spongy tissue disease100 100 100 空心病坏损区域
Damaged area from hollow disease100 100 100 海绵组织病芒果 Mango with spongy tissue disease 100 100 100 空心病芒果 Mango with hollow disease 100 100 100 GoogLeNet 0.0008 无病症芒果 Healthy mango 97.93 99.77 98.84 98.05 海绵组织病坏损区域
Damaged area from spongy tissue disease99.21 98.70 98.96 空心病坏损区域
Damaged area from hollow disease93.51 100 96.63 海绵组织病芒果 Mango with spongy tissue disease 98.42 96.88 97.64 空心病芒果 Mango with hollow disease 100 93.06 96.41 0.0005 无病症芒果 Healthy mango 100 100 100 100 海绵组织病坏损区域
Damaged area from spongy tissue disease100 100 100 空心病坏损区域
Damaged area from hollow disease100 100 100 海绵组织病芒果 Mango with spongy tissue disease 100 100 100 空心病芒果 Mango with hollow disease 100 100 100 0.0002 无病症芒果 Healthy mango 100 100 100 100 海绵组织病坏损区域
Damaged area from spongy tissue disease100 100 100 空心病坏损区域
Damaged area from hollow disease100 100 100 海绵组织病芒果 Mango with spongy tissue disease 100 100 100 空心病芒果 Mango with hollow disease 100 100 100 3.4 模型分类结果分析
加载图像处理数据集，训练深度迁移学习模型。深度迁移学习模型在图像处理数据集下的分类模型分类指标较未处理的数据集下有较大的提升。
AlexNet模型在学习率0.0008、0.0005和0.0002下，Accuracy分别为95.39%、97.40%和100%，Macro-average分别为94.95%、96.93%和100%。GoogLeNet模型在学习率0.0008、0.0005和0.0002下，Accuracy分别为98.05%、100%和100%，Macro-average分别为97.70%、100%和100%。基于同一数据集并且超参数一致时，GoogLeNet网络的分类效果明显优于AlexNet网络。GoogLeNet网络的训练时间较长。在学习率为0.000 5和0.000 2时，基于图像处理数据集得到的GoogLeNet分类模型，达到了100%的分类精度。相较于芒果图像识别算法的计数识别准确率82.3%[24]，基于改进Faster R-CNN的芒果识别平均识别精度均值为88.94%[28]，有极大提升。
综合考虑模型的性能参数和训练时间，选用基于GoogLeNet网络深度迁移学习模型在学习率为0.0002、迭代轮数为3、最小批值为64时，开展训练，将所得模型作为最终的分类模型。
4. 结论
1)针对芒果CT序列图像坏损区域获取的问题，本文采用分段函数法、中值滤波结合双边滤波的图像处理方案实现图像增强，通过局部阈值自适应法开展图像二值化处理，然后通过形态学种子填充法进行填充，最后利用差影法提取了芒果内部组织的坏损区域，能够准确地分割出坏损区域，分割效果较好。
2)本文训练AlexNet和GoogLeNet分类模型，通过比较Accuracy和Macro-average分类性能参数，GoogLeNet网络对芒果图像分类识别的效果最优。当学习速率为0.0002，基于未处理数据集和处理过的数据集进行分类，GoogLeNet网络分类模型性能极佳，可以作为本文芒果数据的分类模型。深度迁移学习模型在水果图像分类领域大有可为。
3)本文分割提取后的芒果坏损区域图像，特征明显，易于区分，在类型识别过程中可达到极高准确度，体现了结合图像处理技术开展深度学习分类识别，具有较高的实际意义。
4)算法的局限性。文中的芒果数据并不能代表所有的海绵组织病和空心病芒果的坏损情况，主要是由于芒果有自己的生长周期，并且坏损情况受到温湿度、土壤水平和种植手段等多方面的影响，在未来的研究中，可以通过增加数据量来探究某类型病状芒果内部坏损的演变过程。
-
表 1 单作与间作处理水稻生物量、产量及单株产量1)
Table 1 Biomass, yield per unit area and yield per plant of rice in monoculture and intercropping treatments
处理
Treatment生物量/g Biomass 总有效穗数
Total effective panicles每穗粒数
Grains per panicle千粒质量/g
1 000-grain weight产量/
(g·m−2)
Yield单株产量/g
Yield per plant地上部
Shoot地下部
Root水稻单作(对照)
Rice monocropping (control)11.63±1.24b 0.87±0.14a 72.00±2.08a 117.30±2.49a 21.93±0.70a 185.34±10.23a 11.58±0.64ab 水稻与圆叶白苋菜间作
Rice-white amaranth intercropping18.99±0.62a 0.74±0.06a 43.33±1.86b 117.17±3.62a 21.59±0.34a 109.65±6.50b 13.71±0.81a 水稻与籽粒苋间作
Rice-grain amaranth intercropping8.36±0.88c 0.61±0.10a 35.00±0.58c 116.43±1.51a 21.01±0.59a 85.67±3.72b 10.71±0.47b 水稻与红苋菜间作
Rice-red amaranth intercropping10.76±0.43bc 0.93±0.19a 34.33±2.40c 111.61±1.72a 22.52±0.53a 85.78±6.57b 10.72±0.82b 1) 同列数据后的不同小写字母表示差异显著(P<0.05，Duncan’ s 法)
1) Different lowercase letters in the same column indicate significant differences(P<0.05, Duncan’ s method)
下载: 导出CSV
表 2 单作与间作处理水稻各部位As含量1)
Table 2 As contents in various parts of rice in monoculture and intercropping treatments
w(As)/(mg·kg−1) 处理
Treatment地上部
Shoot地下部
Root谷壳
Husk糙米
Brown rice水稻单作(对照)
Rice monocropping (control)11.74±1.16b 103.63±11.66b 0.63±0.10b 0.44±0.03c 水稻与圆叶白苋菜间作
Rice-white amaranth intercropping22.13±1.97a 171.38±12.95a 1.06±0.21a 0.64±0.02a 水稻与籽粒苋间作
Rice-grain amaranth intercropping20.55±1.76a 167.43±13.21a 0.86±0.04a 0.52±0.01b 水稻与红苋菜间作
Rice-red amaranth intercropping18.81±1.03a 133.33±12.60ab 0.80±0.04ab 0.52±0.01b 1)同列数据后的不同小写字母表示差异显著(P<0.05，Duncan’ s 法)
1) Different lowercase letters in the same column indicate significant differences(P<0.05, Duncan’ s method)
下载: 导出CSV
表 3 不同品种苋菜单株生物量及As含量1)
Table 3 Biomass per plant and As contents in various parts of different amaranth varieties
品种
Variety生物量/g Biomass w(As)/(mg·kg−1) 地上部
Shoot地下部
Root地上部
Shoot地下部
Root圆叶白苋菜 White amaranth 26.57±2.10a 5.20±1.35a 1.54±0.07b 16.14±0.81b 籽粒苋 Grain amaranth 14.56±0.23b 0.44±0.07b 3.18±0.13a 13.76±0.81b 红苋菜 Red amaranth 24.81±2.46a 1.48±0.02b 1.22±0.18b 25.42±0.37a 1)同列数据后的不同小写字母表示差异显著(P<0.05，Duncan’ s 法)
1) Different lowercase letters in the same column indicate significant differences(P<0.05, Duncan’ s method)
下载: 导出CSV
表 4 单作与间作处理植株不同部位As累积量及土壤中As含量1)
Table 4 As accumulations in various parts of plants and As contents in soil in monoculture and intercropping treatments
处理
TreatmentAs累积量/(mg·m−2) As accumulations w(As)土壤/
(mg·kg−1)
As concent in soil水稻地上部
Rice shoot水稻地下部
Rice root苋菜地上部
Amaranth shoot苋菜地下部
Amaranth root总和
Sum水稻单作(对照)
Rice monocropping (control)2.21±0.29b 1.04±0.09a 3.25±0.38b 63.53±0.96a 水稻与圆叶白苋菜间作
Rice-white amaranth intercropping3.36±0.24a 1.05±0.16a 0.33±0.03a 0.69±0.16a 5.43±0.34a 59.81±1.49a 水稻与籽粒苋间作
Rice-grain amaranth intercropping1.39±0.17b 0.84±0.04a 0.37±0.01a 0.05±0.01b 2.65±0.20b 65.18±1.68a 水稻与红苋菜间作
Rice-red amaranth intercropping1.63±0.12b 0.76±0.07a 0.24±0.01b 0.30±0.00b 2.93±0.16b 63.64±2.03a 1)同列数据后的不同小写字母表示差异显著(P<0.05，Duncan’ s 法)
1) Different lowercase letters in the same column indicate significant differences(P<0.05, Duncan’ s method)
下载: 导出CSV
表 5 单作与间作水稻根际土壤速效养分与有机质含量1)
Table 5 Available nutrient and organic matter contents in rice rhizosphere soil in monoculture and intercropping treatments
处理
Treatmentw(铵态氮)/(mg·kg−1)
Ammonia nitrogen contentw(速效磷)/(mg·kg−1)
Available phosphorous contentw(速效钾)/(mg·kg−1)
Available potassium contentw(有机质)/(g·kg−1)
Organic matter content水稻单作(对照)
Rice monocropping (control)25.86±0.89b 126.81±3.64b 91.67±5.37b 36.09±0.12c 水稻与圆叶白苋菜间作
Rice-white amaranth intercropping35.73±0.45a 139.35±2.12a 106.60±1.27a 45.92±1.84a 水稻与籽粒苋间作
Rice-grain amaranth intercropping27.37±0.56b 128.44±3.00ab 99.17±3.87ab 40.69±0.36b 水稻与红苋菜间作
Rice-red amaranth intercropping27.90±0.89b 130.83±3.93ab 98.40±0.75ab 41.26±1.85b 1)同列数据后的不同小写字母表示差异显著(P<0.05，Duncan’ s 法)
1) Different lowercase letters in the same column indicate significant differences(P<0.05, Duncan’ s method)
下载: 导出CSV
-
[1] 环境保护部国土资源部. 全国土壤污染状况调查公报[EB/OL]. [2014-04-17] http://www.mep.gov.cn/gkml/hbb/qt/201404/t20140417_270670.Htm. [2] 刘小诗, 李莲芳, 曾希柏, 等. 典型农业土壤重金属的累积特征与源解析[J]. 核农学报, 2014, 28(7): 1288-1297. doi: 10.11869/j.issn.100-8551.2014.07.1288 [3] 钟松雄, 何宏飞, 陈志良, 等. 水淹条件下水稻土中砷的生物化学行为研究进展[J]. 土壤学报, 2018, 55(1): 1-17. doi: 10.11766/trxb201704250028 [4] 赵宏亮, 王麒, 孙羽, 等. 秸秆还田下灌溉方式对水稻产量及水分利用率的影响[J]. 核农学报, 2018, 32(5): 959-969. doi: 10.11869/j.issn.100-8551.2018.05.0959 [5] LIU C P, LUO C L, GAO Y, et al. Arsenic contamination and potential health risk implications at an abandoned tungsten mine, southern China[J]. Environ Pollut, 2010, 158(3): 820-826. doi: 10.1016/j.envpol.2009.09.029
[6] 黄亚涛. 我国稻米中无机砷的污染分布研究及风险评估[D]. 北京: 中国农业科学院, 2014: 38-44. [7] 胡拥军, 王海娟, 王宏镔, 等. 砷胁迫下不同砷富集能力植物内源生长素与抗氧化酶的关系[J]. 生态学报, 2015, 35(10): 3214-3224. [8] 沈孝辉, 李仁英, 徐向华, 等. 土壤−水稻系统砷迁移累积的影响因素及调控措施[J]. 土壤通报, 2014, 45(5): 1273-1280. [9] 吴川, 莫竞瑜, 薛生国, 等. 不同渗氧能力水稻品种对砷的耐性和积累[J]. 生态学报, 2014, 34(4): 807-813. [10] TALUKDER A S M H M, MEISNER C A, SARKAR M A R, et al. Effect of water management, arsenic and phosphorus levels on rice in a high-arsenic soil-water system: II: Arsenic uptake[J]. Ecotoxicol Environ Saf, 2012, 80: 145-151. doi: 10.1016/j.ecoenv.2012.02.020
[11] ZHAO S T, LI H C, SHI Y H. Speciation analysis of antimony and arsenic in soil and remediation of antimony and arsenic in contaminated soils[J]. Adv Mater Res, 2015, 1088: 578-582. doi: 10.4028/www.scientific.net/AMR.1088.578
[12] CHEN J, WANG S, ZHANG S, et al. Arsenic pollution and its treatment in Yangzonghai lake in China: In situ remediation[J]. Ecotoxicol Environ Saf, 2015, 122: 178-185. doi: 10.1016/j.ecoenv.2015.07.032
[13] KLABER N S, BARKER A V. Accumulation of phosphorus and arsenic in two perennial grasses for soil remediation[J]. Commun Soil Sci Plant Anal, 2014, 45(6): 810-818. doi: 10.1080/00103624.2013.857681
[14] 申红玲, 何振艳, 麻密. 蜈蚣草砷超富集机制及其在砷污染修复中的应用[J]. 植物生理学报, 2014, 50(5): 591-598. [15] 宋书巧, 周永章, 周兴, 等. 土壤砷污染特点与植物修复探讨[J]. 热带地理, 2004, 24(1): 6-9. doi: 10.3969/j.issn.1001-5221.2004.01.002 [16] AMPIAH-BONNEY R J, TYSON J F, LANZA G R. Phytoextraction of arsenic from soil by Leersia oryzoides[J]. Int J Phytoremediation, 2007, 9(1/2/3): 31-40.
[17] 张少斌, 梁开明, 张殷, 等. 水稻与水合欢间作对作物群体产量、氮素吸收及土壤氮素的影响[J]. 生态环境学报, 2016, 25(11): 1856-1864. [18] 梁开明, 傅玲, 章家恩, 等. 水稻/再力花单、间作系统中作物对镉胁迫的叶绿素荧光响应及镉积累特征[J]. 华南农业大学学报, 2014, 35(4): 35-41. doi: 10.7671/j.issn.1001-411X.2014.04.008 [19] 宁川川, 杨荣双, 蔡茂霞, 等. 水稻−雍菜间作系统中种间关系和水稻的硅、氮营养状况[J]. 应用生态学报, 2017, 28(2): 474-484. [20] 吴华杰, 李隆, 张福锁. 水稻/小麦间作中种间相互作用对镉吸收的影响[J]. 中国农业科技导报, 2003, 5(5): 43-46. doi: 10.3969/j.issn.1008-0864.2003.05.010 [21] 张骞, 曾希柏, 苏世鸣, 等. 不同品种苋菜对砷的吸收能力及植株磷砷关系研究[J]. 农业环境科学学报, 2016, 35(10): 1888-1894. doi: 10.11654/jaes.2016-0687 [22] 中华人民共和国农业部. 酸性土壤铵态氮、有效磷、速效钾的测定联合浸提−比色法: NY/T 1849—2010[S]. 北京: 中国农业出版社, 2010. [23] 中华人民共和国国家卫生和计划生育委员会, 国家食品药品监督管理总局. 食品中污染物限量: GB 2762—2017[S]. 北京: 中国标准出版社, 2017. [24] GOVE B, HUTCHINSON J J, YOUNG S D, et al. Uptake of metals by plants sharing a rhizosphere with the hyperaccumulator Thlaspi caerulescens[J]. Int J Phytoremediation, 2002, 4(4): 267-281. doi: 10.1080/15226510208500087
[25] WHITING S N, DE SOUZA M P, TERRY N. Rhizosphere bacteria mobilize Zn for hyperaccumulation by Thlaspi caerulescens[J]. Environ Sci Technol, 2001, 35(15): 3144-3150. doi: 10.1021/es001938v
[26] 赵冰, 沈丽波, 程苗苗, 等. 麦季间作伴矿景天对不同土壤小麦−水稻生长及锌镉吸收性的影响[J]. 应用生态学报, 2011, 22(10): 2725-2731. [27] 李凝玉, 李志安, 丁永祯, 等. 不同作物与玉米间作对玉米吸收积累镉的影响[J]. 应用生态学报, 2008, 19(6): 1369-1373. [28] 张静静, 黄河, 吴海霞, 等. 不同土壤类型中水稻间作红蛋对镉累积的影响[J]. 西南农业学报, 2016, 29(12): 2871-2876. [29] 刘维涛, 倪均成, 周启星, 等. 重金属富集植物生物质的处置技术研究进展[J]. 农业环境科学学报, 2014, 33(1): 15-27. doi: 10.11654/jaes.2014.01.002 [30] 章家恩, 高爱霞, 徐华勤, 等. 玉米/花生间作对土壤微生物和土壤养分状况的影响[J]. 应用生态学报, 2009, 20(7): 1597-1602. [31] 曾露苹, 周飞, 陈玥如, 等. 木薯和花生间作对Cd吸收及根区速效养分的影响[J]. 环境科学研究, 2018, 31(2): 303-309. [32] IKERRA S T, MAGHEMBE J A, SMITHSON P C, et al. Soil nitrogen dynamics and relationships with maize yields in a gliricidia-maize intercrop in Malawi[J]. Plant Soil, 1999, 211(2): 155-164. doi: 10.1023/A:1004636501488
[33] LI B, LI Y Y, WU H M, et al. Root exudates drive interspecific facilitation by enhancing nodulation and N2 fixation[J]. P Natl Acad Sci USA, 2016, 113(23): 6496-6501. doi: 10.1073/pnas.1523580113
[34] NORTON G J, ADOMAKO E E, DEACON C M, et al. Effect of organic matter amendment, arsenic amendment and water management regime on rice grain arsenic species[J]. Environ Pollut, 2013, 177: 38-47. doi: 10.1016/j.envpol.2013.01.049
[35] ROWLAND H A L, BOOTHMAN C, PANCOST R, et al. The role of indigenous microorganisms in the biodegradation of naturally occurring petroleum, the reduction of iron, and the mobilization of arsenite from west bengal aquifer sediments[J]. J Environ Qual, 2009, 38(4): 1598-1607. doi: 10.2134/jeq2008.0223
-
期刊类型引用(10)
1. 蔡若妍，赵珠宇，谭力，薛锦春，张晓. 硫化铜矿酸性环境调节技术研究. 环境监测管理与技术. 2024(05): 1-6+12 . 
百度学术
2. 薛迎斌，宋佳，李枭艺，李小豪，陈经烨，伍萍珍，朱胜男，刘颖. 大豆GmMADS4基因克隆、亚细胞定位及功能分析. 华南农业大学学报. 2023(03): 420-429 . 本站查看
3. 何春梅，王飞，黄毅斌，刘彩玲，曹卫东，李清华，张辉. 紫云英翻压条件下化肥减施对黄泥田水稻产量与土壤肥力的影响. 植物营养与肥料学报. 2023(06): 1070-1081 . 百度学术
4. 黄家新，薛锦春，王伟伟. 城门山铜矿尾矿库酸性废水处理试验. 现代矿业. 2023(08): 199-203 . 百度学术
5. 杜梦柯，连文婷，张晓，李欣欣. 氮处理对大豆根瘤固氮能力及GmLbs基因表达的影响. 植物学报. 2021(04): 391-403 . 百度学术
6. 董艳鑫，蒋景龙，李丽，阳妮，焦成瑾，陈德经. 健康和患根腐病西洋参根水提液化感作用差异性分析. 北方园艺. 2021(12): 9-16 . 百度学术
7. 王静，刘晓静，郝凤，汪雪. 不同氮效率紫花苜蓿各生育期氮利用特征研究. 草地学报. 2021(11): 2461-2469 . 百度学术
8. 蔡秋亮，钟宁，方显瑞，陈玲，张传芳. 硫酸铵处理下茶园土壤氮素转化和交换态金属离子含量的动态变化. 江西农业学报. 2021(11): 83-90 . 百度学术
9. 朱亚琼，简大为，郑伟，王朴，黎松松，郝帅，娜尔克孜，刘岳含，艾丽菲热. 不同种植模式下豆科绿肥对土壤改良效果的影响. 草业科学. 2020(05): 889-900 . 百度学术
10. 蔡媛，刘浩，王勇庆，谢景，黄建华，劳嘉，贺炜，张水寒. 基于16S rRNA基因高通量测序方法分析多花黄精内生细菌群落结构及多样性. 湖南中医药大学学报. 2020(07): 846-852 . 百度学术
其他类型引用(14)
