Processing math: 0%
    • EI
    • CSA
    • CABI
    • 卓越期刊
    • CA
    • Scopus
    • CSCD
    • 核心期刊

河北省耕地多功能评价与关键权衡协同分析

武子豪, 郝晋珉, 陈航, 谭永忠

武子豪,郝晋珉,陈航,等. 河北省耕地多功能评价与关键权衡协同分析[J]. 农业工程学报,2024,40(14):199-209. DOI: 10.11975/j.issn.1002-6819.202403052
引用本文: 武子豪,郝晋珉,陈航,等. 河北省耕地多功能评价与关键权衡协同分析[J]. 农业工程学报,2024,40(14):199-209. DOI: 10.11975/j.issn.1002-6819.202403052
WU Zihao, HAO Jinmin, CHEN Hang, et al. Multifunctional evaluation and key trade-offs and synergy relationships of cultivated land in Hebei province of China[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2024, 40(14): 199-209. DOI: 10.11975/j.issn.1002-6819.202403052
Citation: WU Zihao, HAO Jinmin, CHEN Hang, et al. Multifunctional evaluation and key trade-offs and synergy relationships of cultivated land in Hebei province of China[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2024, 40(14): 199-209. DOI: 10.11975/j.issn.1002-6819.202403052

河北省耕地多功能评价与关键权衡协同分析

基金项目: 国家自然科学基金项目(42071269);浙江省哲学社会科学规划项目(21WZQH12YB);中央高校基本科研业务费专项资金资助(S20240008)
详细信息
    作者简介:

    武子豪,博士生,研究方向为国土空间利用与优化。Email:12322005@zju.edu.cn

    通讯作者:

    谭永忠,教授,主要研究方向为耕地资源利用与保护。Email:tanyongzhong@zju.edu.cn

  • 中图分类号: F323.211

Multifunctional evaluation and key trade-offs and synergy relationships of cultivated land in Hebei province of China

  • 摘要:

    耕地多功能评价有助于显化耕地多功能价值,通过探究功能之间权衡协同的关键区域与作用类型,提供更合理的耕地利用与优化措施。该研究以河北省为研究区域,构建功能评价指标体系分类评价县域耕地生产、生态和生活功能,基于Spearman秩相关和空间自相关衡量功能相关关系,并判断相关关系的区域差异,由此划分耕地功能分区与核心区。结果表明:1)2010―2020年河北省耕地生产功能维持稳定,生态功能有所下降,生活功能逐步提升,三类功能之间均为协同关系;2)功能协同的关键区域多为功能高值区,部分功能之间的相关关系会随着区域功能强度的不同发生变化;3)粮食生产、景观美学与粮食安全保障子功能的相互作用对于三类功能协同关系贡献较大;4)基于耕地功能强度划分6类耕地功能优势区,进一步基于空间自相关分析识别出3类功能核心区。分析耕地功能之间权衡协同的区域差异,确定发挥重要贡献的子功能作用类型,以此划分耕地功能分区,以期提供更精确的耕地利用优化方案。

    Abstract:

    In 2021, China achieved a comprehensive victory in the battle against poverty, and fully advancing the work of rural revitalization became the new focus of efforts. As an important component of rural areas, cultivated land is the material foundation and important support for rural revitalization.The multifunctional evaluation of cultivated land can help to reveal the multifunctional value of cultivated land, providing more reasonable cultivated land use and optimization measures by exploring the key areas and types of trade-offs and synergy between functions. It is of great significance for advancing the rural revitalization strategy and achieving the modernization of agriculture and rural areas. Taking Hebei Province as the research area, this paper constructs a functional evaluation index system to classify and evaluate cultivated land's production, ecological and living functions and subfunctions. Based on Spearman rank correlation and spatial autocorrelation, we examined the trade-offs and synergy among different cultivated land functions, used the natural discontinuity classification method to divide cultivated land functions into secondary, tertiary and fourth-level classifications, compared the Spearman rank correlation coefficients between the functions in each group of data, determined regional differences in functional correlation, identified key functional trade-off and synergy regions as well as sub-functional role types, and finally cultivated land functional zones and core areas were divided. The results show that: (1) From 2010 to 2020, the production function of cultivated land in Hebei Province remained stable, the ecological function decreased, and the living function gradually improved. The three functions are all in a synergy relationship and apart from a weakening in the synergy between living and ecological functions from 2015 to 2020, the degree of synergy among other types of functions has gradually increased. (2) The key areas for functional synergy are mostly high-value areas. In the zoning based on the level of production function, the relationship between production function and ecological function shifted from low-value synergy to high-value trade-off. And in the zoning based on the level of living function, the synergy between the living function and the other two types of functions decreased as the living function improved. (3) There are more complex trade-offs and synergy among the subfunctions of cultivated land in Hebei Province. The synergy between landscape aesthetic subfunction and the three types of production subfunctions is the key to the synergy between production function and ecological function, and the synergy between landscape aesthetic subfunction and food security subfunction contributes the most to the synergistic relationship between ecological function and living function, and the synergy between grain production subfunction, large-scale production subfunction and food security subfunction is the key synergy between production function and living function. (4) Based on the values of cultivated land function, 6 types of cultivated land function advantage areas are divided, and 3 types of core functional areas are identified by combining spatial autocorrelation analysis. Our research analyzes the regional differences in trade-offs and synergy of cultivated land functions, identifies the important types of subfunctions that make significant contributions to functional synergy and divides cultivated land functional zones accordingly, to provide a more accurate optimization scheme for cultivated land use.

  • 中国是全球最大的桃生产国,2021年中国桃子种植面积和产量分别占全球的54.83%和64.08%,桃产业总产值近千亿元[1]。然而,桃树缩叶病是桃树的主要病害之一,其主要发生在嫩叶,严重时会危害嫩枝、花朵和果实,造成桃树减产,甚至影响桃树次年产量。桃树缩叶病的快速和准确的识别,对防治措施具有重要意义。目前桃树缩叶病一般是果农依靠经验识别,这种方式不仅识别准确率和效率低,还会受到果农主观因素影响。因此,桃树缩叶病发病情况的智能检测是桃产业急需实现的目标。

    近年来随着计算机硬件技术的发展,基于深度学习的目标检测方法被广泛应用于农业领域的识别任务中[2-4]。王磊磊等[5]在YOLOv5的模型中添加注意力模块,并改进损失函数构建了OMM-YOLO模型用于平菇目标检测,改进模型对各个成熟度的平菇检测精度均有提高。WANG等[6]提出了一种基于两阶段移动视觉的级联害虫检测方法,该方法将作物分类模型作为预训练模型训练害虫检测,通过预先对作物种类进行区分再进行害虫识别解决物种间图像数据不平衡的问题。LI等[7]提出了一种基于改进YOLOv5的蔬菜病害检测方法,该方法采用多尺度特征融合策略提高了特征提取能力,减少了由于复杂背景造成的漏检和误检。兰玉彬等[8]在YOLOv5s模型中加入CA(coordinate attention,CA)注意力减少小目标特征丢失的情况,并引入GhostNet网络中的Ghost模块轻量化网络,实现了对自然场景下生姜叶片病虫害识别准确率的提升。SU等[9]在YOLOv5模型中加入SE(squeeze-and-excitation,SE)解决了传统目标检测模型无法有效筛选出芸豆褐斑病的关键特征的问题。XIAO等[10]利用增强的YOLOv7和边缘计算对荔枝病害进行实时轻量级检测,相比原始YOLOv7模型检测速度更快,模型参数更少。HOU等[11]在Faster-RCNN网络中引入特征金字塔结构和ISResNet(inception Squeeze-and-Excitation ResNet)构建了一种深度学习模型,该模型在检测苹果叶片病害方面具有较高的准确性和通用性。李颀[12]采用基于回归的轻量型检测算法RFBNet(reception field block network,RFBNet)检测桃树的病虫害,其方法识别速度快,对桃树缩叶病的识别精度达到了76.58%。YANG等[13]提出了一种基于YOLOv8s深度学习算法并结合LW-Swin Transformer模块的LS-YOLOv8s的草莓成熟度检测和分级模型,增加了Swin Transformer模块,利用多头自注意力机制提高模型的泛化能力。ZHANG等[14]提出了一种基于深度学习的苹果病害检测算法,设计了一种Bole卷积模块(Bole convolution module,BCM)减少苹果叶病图像特征提取过程中的冗余特征信息,提出了一种交叉注意力模块(cross-attention module,CAM)降低图像背景信息对模型的干扰。上述研究在农作物病害检测等方面取得了较好的效果,但仍存在如下问题:

    多数研究仅在单一叶片上进行,对实际场景中出现的遮挡、光线等影响因素考虑不足;多数研究仅对病害进行分类,并未区分病害发展情况;针对桃树缩叶病识别的研究较少。

    针对以上问题,本文以果园中桃树为目标,采集并搭建桃树缩叶病目标检测数据集。以YOLOv5su为基础模型,通过引入可变形自注意力机制、可分离大核卷积注意力和轻量化权重自适应下采样模块提高模型对特征的拟合能力,提出了一种基于改进YOLOv5su的桃树缩叶病检测模型DLL-YOLOv5su,为实现桃树缩叶病发生的实时监测提供一种有效方法。

    本文的桃树缩叶病图像采集地点为中国重庆市农业科学院桃树种植基地,以“紫金红”油桃树为主要研究对象。采集时间为2023年4—6月,采集设备与方式分别为数码相机Canon EOS M50 m2拍摄照片和手机SAMSUNG SM-G9910拍摄视频,设备距离目标20~80 cm,拍摄不同时段和视角下的桃树缩叶病图像,包含从桃树缩叶病发生早期到后期的各阶段图像。拍摄得到291张分辨率为6 000×4 000像素的桃树图像和5分38秒视频,视频帧率为60帧/s,分辨率为3 840×2 160像素。数码相机采集到的图像分辨率高,且包含的内容多,故对291张高分辨率图像进行裁剪,筛选出包含缩叶病目标的图像1 358张。对视频进行抽帧和筛选,剔除过于模糊和不包含目标的背景图片得到162张桃树缩叶病图像。合并两个图像集得到1 520张桃树缩叶病图像。

    采用轻量级图像注释工具“LabelImg”将采集到的桃树缩叶病图像按照YOLO格式标注,并以8:2:1的比例随机划分为训练集1 105张、验证集277张和测试集138张。由于数据集均于同一场景中拍摄,为避免过拟合,提高模型鲁棒性和泛化能力,采用添加高斯噪声、改变亮度、Cutout、旋转、翻转、平移等图像处理方法随机组合对训练集进行数据增广,扩充后的训练集包含6 630张图片。

    图1所示,在桃树缩叶病数据集中,每一张图片均包含病害目标,且病害分布随机,包含多种光照强度与方向。其中早期缩叶病标签数一共为19 419,中期缩叶病标签数为11 172,后期缩叶病标签数为10 069。

    图  1  桃树缩叶病不同发生时期示例
    Figure  1.  Examples of peach leaf curl disease in different periods

    YOLO是目前流行的目标检测和图像分割模型之一,由华盛顿大学的Joseph Redmon和Ali Farhadi开发,因其高速度和高精度的优点而迅速流行起来。YOLOv5主要由输入(Input)、骨干(Backbone)、颈部(Neck)和头部(Head)组成。

    YOLOv5u[15]在YOLOv5模型的基础架构上,集成了解耦头部结构,将分类任务与回归任务解耦。同时采用了无需预定义锚框匹配真实框的Anchor-free策略避免预定义锚框尺寸不合理的问题。YOLOv5u根据网络深度和宽度由小到大分为YOLOv5nu、YOLOv5su、YOLOv5mu、YOLOv5lu、YOLOv5xu五种不同结构。

    本研究的试验图像皆取自室外自然环境中的桃树缩叶病发病区域。为了减轻光照和复杂背景环境等因素的不利影响,提高目标检测模型的鲁棒性,对YOLOv5su算法进行了改进,构建了DLL-YOLOv5su。具体改进措施如下:

    1)在骨干网络最后一层C3模块的Bottleneck结构中加入Version Transformer [16]的可变形自注意(DA)模块[17]组成新的C3-DA模块。

    2)在快速空间金字塔池化模块(SPPF)中加入了Version Attention Network [18]中的可分离大核卷积注意力模块[19]构成SPPF-LSKA模块。

    3)在感受野注意力卷积(receptive-field attention convolution,RFAConv)[20]的基础上提出了一种轻量化权重自适应下采样模块(lightweight adaptive weighted downsampling,LAWD)替换部分初始网络结构中的CBS模块。

    改进后的DLL-YOLOv5su模型结构如图2所示。

    图  2  DLL-YOLOv5su模型结构图
    注:Upsample为上采样模块,Concat为从通道维度进行特征图拼接,Conv2 d为卷积层,BN表示批量归一化层,SiLU为激活函数,Bbox.Loss为定位损失,Cls.Loss为分类损失。
    Figure  2.  DLL-YOLOv5su model structure diagram
    Note: Upsample is the upsampling module, Concat is the feature map concatenation from the channel dimension, Conv2 d is the convolutional layer, BN represents the batch normalization layer, SiLU is the activation function, Bbox.Loss is the regression loss, and Cls.Loss is the classification loss.

    自注意力最早在自然语言处理模型Transformer[21]中被提出。视觉转换器(version transformer,ViT)将自注意力机制用于视觉领域,构造了多头注意力模块(multi-head self-attention)。自注意力有三个重要元素:查询(query)、键(key)和值(value)。其输出式如下:

    \text{Attention(}\boldsymbol{q}{,}{}\boldsymbol{k}\text{,}{}\boldsymbol{v}\text{)=Softmax(}\boldsymbol{q}{\boldsymbol{k}}^{\text{T}}/\sqrt{{{d}}_{\boldsymbol{k}}}\text{)}\boldsymbol{v} (1)

    式中qkv分别表示查询向量、键向量和值向量, {d}_{k} k的维度,Softmax表示归一化。

    自注意力可以捕获全局信息,并建立特征通道和目标之间的长期依赖关系。然而简单地扩大感受野,会导致内存和计算资源的浪费,并且特征会受到背景区域影响[17]。为了有效地建模和捕获全局语义信息,在主干网络最后一层C3模块中加入可变形自注意力模块强化缩叶病目标信息,抑制背景信息。可变形自注意力使用多组可变形的采样点确定特征图的重要区域,并以此重要区域为依据对桃树缩叶病不同特征之间的关系进行建模。这种灵活的方案使变形自注意模块能够专注于相关区域并捕获更多信息特征。可变形自注意力结构如图3所示。

    图  3  变形注意力模块与偏移网络结构图
    注:q表示查询向量, \tilde{\boldsymbol{k}} \tilde{\boldsymbol{v}} 分别表示变形的键和值向量, {\boldsymbol{W}}_{\boldsymbol{q}}{{、}\boldsymbol{W}}_{\boldsymbol{k}}\mathrm{、}{\boldsymbol{W}}_{\boldsymbol{v}}\mathrm{和}{\boldsymbol{W}}_{\boldsymbol{o}} 分别表示查询向量、键向量、值向量和输出向量的投影矩阵,z表示输出特征图,HWC分别表示图像高度、宽度和通道数,r为下采样缩减比例,GroupConv为组卷积,LayerNorm为层归一化,GELU为激活函数。
    Figure  3.  Deformable attention module and offset network structure diagram
    Note: q is the query vector, \tilde{\boldsymbol{k}} and \tilde{\boldsymbol{v}} represent the deformed key and value vectors, \tilde{\boldsymbol{x}} represents the feature map sampling result, {\boldsymbol{W}}_{\boldsymbol{q}}{{,}\boldsymbol{W}}_{\boldsymbol{k}}{,}{\boldsymbol{W}}_{\boldsymbol{v}}\;{{\mathrm{and}}}\;{\boldsymbol{W}}_{\boldsymbol{o}} denote the projection matrices of the query vector, key vector, value vector and output vector respectively, z is the output feature map, H, W and C represent the image height, width, and number of channels, r is the downsampling reduction ratio, GroupConv is a group convolution, LayerNorm is a layer normalization, and GELU is an activation function.

    图3a所示,变形注意力模块首先对输入特征图 \boldsymbol{x} 生成一个均匀的二维点网格 {p} ,网格中的点作为采样点。将特征图线性投影到查询向量,然后馈送到偏移子网络 \theta_{\text{offset}} 以生成偏移量 \Delta {p} ,使用此偏移量与采样点相加得到变形点。在变形点处利用双线性插值函数进行采样得到采样特征 \tilde{\boldsymbol{x}} 。其式为:

    \boldsymbol{q}=\boldsymbol{x}{{W}}_{\boldsymbol{q}}{,}\;\;\tilde{\boldsymbol{k}}=\tilde{\boldsymbol{x}}{{W}}_{\boldsymbol{k}}{,}\;\;\tilde{\boldsymbol{v}}=\tilde{\boldsymbol{x}}{{W}}_{\boldsymbol{v}} (2)
    \Delta {p}\text=\theta_{\text{offset}}{(}\boldsymbol{q}\text{)} (3)
    \tilde{\boldsymbol{x}}=\phi\text{(}\boldsymbol{x}\text{;}{p}\text{+}\Delta {p}\text{)} (4)
    \phi {(}\boldsymbol{z}{;(}{{p}}_{{x}}{,}{{p}}_{{y}}{))}={\sum }_{{(}{{r}}_{{x}}{,}{{r}}_{{y}}{)}}{g}{(}{{p}}_{{x}}{,}{{r}}_{{x}}{)}{g}{(}{{p}}_{{y}}{,}{{r}}_{{y}}{)}\boldsymbol{z}{[}{{r}}_{{y}}{,}{{r}}_{{x}}{,:]} (5)

    式中 {\boldsymbol{W}}_{\boldsymbol{q}}{{、}\boldsymbol{W}}_{\boldsymbol{k}}{、}{\boldsymbol{W}}_{\boldsymbol{v}} 分别为查询向量、键向量和值向量的投影矩阵, {g}{(}{a}{,}{}{b}{)}{=}{{\mathrm{max}}(0,1-}{|}{a}{-}{b}{|}{)} ( {{r}}_{{x}}{,}{{r}}_{{y}} )索引\boldsymbol{z}\in {{R}}^{{H\times W\times C}} 上的所有位置。由于 {g}{(}{a}{,}{}{b}{)} 仅在最接近 {(}{{p}}_{{x}}{,}{{p}}_{{y}}{)} 的4个积分点上非零,可以简化式(5)为

    \begin{split} & {\boldsymbol{z}}^{{(}{m}{)}}= {\sigma}\left(\frac{{\boldsymbol{q}}^{{(}{m}{)}}{\tilde{\boldsymbol{k}}}^{{(}{m}{)T}}}{\sqrt{{d}}}+\phi {(}\widehat{{B}}{;}{R})\right){\tilde{\boldsymbol{v}}}^{{(}{m}{)}}\\ & {m=}{1,…,}{M}{.} \end{split} (6)
    \boldsymbol{z}\text{=Concat(}{\boldsymbol{z}}^{\text{(1)}}\text{)}{\boldsymbol{W}}_{\boldsymbol{o}} (7)

    式中 {\boldsymbol{z}}^{{(}{m}{)}} 表示第m个头部输出向量, \sigma (·)表示Softmax函数,M表示多头注意力头部数量,d=C/M表示每个头部的维度, {\boldsymbol{q}}^{{(}{m}{)}}{{、}\,\tilde{\boldsymbol{k}}}^{{(}{m}{)}}{、}\,{\tilde{\boldsymbol{v}}}^{{(}{m}{)}}\in {{R}}^{{N}{\times}{d}} 分别表示查询、键和值向量的对应头的元素, {\boldsymbol{W}}_{\boldsymbol{o}} 表示输出向量投影矩阵。 \phi {(}\widehat{{B}}{;}{R}{)}\in {{R}}^{{HW}{\times}{HW}} 对应于Swin Transformer[22]中的位置嵌入。将每个头部的特征连接在一起并通过 {\boldsymbol{W}}_{\boldsymbol{o}} 投影得到输出z。将变形注意力模块加入到C3模块中组成C3-Deformable attention模块和Bottleneck-DA模块结构,如图4所示。

    图  4  C3-DA模块结构图
    Figure  4.  C3-DA module structure diagram

    大卷积核注意力(large kernel attention, LKA)模块[18]在VAN中发挥了卓越的性能,使得VAN在一系列基于视觉的任务中的表现超过了ViT和CNN等模型。LKA吸收了卷积和自注意力的优点,包括局部结构信息、长程依赖和适应性,同时避免了通道维度上忽略适应性等缺点。大核卷积模块包括3个部分:深度卷积(depth-wise convolution, DWC)、深度扩张卷积(depth-wise dilation convolution, DWDC)和通道卷积(1×1 convolution)。LKA模块可以表示为

    {\mathrm{LKA}}(\boldsymbol{x})={{C}}^{\text{1×1}} ({\mathrm{DWD}}\_{\mathrm{C}}({\mathrm{DW}}\_{\mathrm{C}}( \boldsymbol{x} ))) (8)

    式中 \boldsymbol{x} 是输入特征图, {\mathrm{LKA}}(\boldsymbol{x}) 是注意力输出特征图, {{C}}^{{1\times 1}} 表示1×1卷积,DW_C(·)表示深度卷积,DWD_C(·)表示深度扩张卷积。

    LSKA将LKA中的深度卷积层的二维卷积核分解为级联的水平和垂直一维卷积核。与标准的LKA相比,LSKA可以提供相当的性能,并且计算复杂度和内存占用更低。为了提高SPPF的特征信息获取能力,在SPPF层中引入可分离大卷积核注意力,构成SPPF-LSKA层,其模型如图5所示。

    图  5  SPPF-LSKA 模块结构图
    注: MaxPool2d为最大池化层,DW-Conv为深度卷积,DW-D-Conv为深度扩张卷积,k表示卷积核,s表示步长,p表示像素填充,d2表示卷积核元素间隔为1。
    Figure  5.  SPPF-LSKA module structure diagram
    Note: MaxPool2 d is the maximum pooled layer, DW-Conv is the deep-wise convolution, DW-D-Conv is the deep-wise dilation convolution, k represents the convolution kernel, s represents the step size, p represents the pixel filling, and d2 represents the convolution kernel element interval of 1.

    标准卷积运算提取特征图是通过卷积核与相同尺寸的感受域相乘然后求和得到。每个感受域相同空间位置处的特征共享相同的卷积核参数。因此,标准卷积操作没有考虑不同位置所包含的信息差异,这在一定程度上限制了卷积神经网络的性能。ZHANG等[23]观察到组卷积可以减少模型的参数和计算开销,但是组内信息之间交互不足会影响网络性能;HUANG等[24]通过重用特征进行融合来改进特征信息,以解决网络梯度消失的问题;DAI等[25]提出了可变形卷积,通过学习偏移来改变卷积核的采样位置,进一步提高了卷积神经网络的性能。ZHNAG等[20]提出了感受野注意卷积(receptive-field attention convolution, RFAConv),该方法动态地确定每个特征在感受野中的权重。

    本文为减少标准卷积的卷积核参数共享导致的空间信息损失与轻量化模型,在RFAConv的基础上提出轻量化权重自适应下采样(lightweight adaptive weight downsampling, LAWD)替换模型中的部分卷积模块。其结构如图6所示。

    图  6  轻量级权重自适应下采样结构图
    注: AvgPool为平均池化层,Rearrange为张量维度调整,K1KC表示第1~C的卷积核,Sum表示求元素和。
    Figure  6.  Lightweight adaptive weight downsampling structure diagram
    Note: Avgpool is the average pooling layer, Rearrange is the tensor dimension adjustment, K1-KC represents the 1st to C convolution kernel, and Sum is the element sum.

    LAWD上半部分为权重路线,首先采用的平均池化聚合每个感受域信息,再通过1×1卷积交互不同通道间的空间信息,提高网络性能[26-28]。对得到的权重矩阵进行维度转换,这样在既不损失权重信息的前提下可以有效减少网络的参数量。

    下半部分为特征提取路线,首先通过组卷积,对输入特征图进行切片操作,将输入特征图的WH维度信息保存到通道中,再进行卷积操作得到的二倍下采样特征图可以很好地保留原始图像信息。相较于普通卷积,采用分组大小为16的分组卷积可以将参数量减少到普通卷积的1/16。再对特征图进行维度变换,并与权重相乘。最后对第四个维度进行求和得到最终下采样的特征图。LAWD可以通过廉价的操作为不同的感受域生成权重矩阵,避免了卷积操作的卷积核参数共享问题,并且减少了模型参数量。LAWD计算过程可以表示为

    \begin{split} {\mathrm{LAWD}}(\boldsymbol{x})=\,& {\mathrm{Sum}}({\mathrm{Softmax}}( {{C}}^{{1 \times 1}} ({\mathrm{Avgpool}}( \boldsymbol{x} )))\\ \,& \times {\mathrm{SiLU}}({\mathrm{BN}}( {{g}}^{\text{3×3}} ( \boldsymbol{x} )))) \end{split} (9)

    式中x表示输入特征图, {{g}}^{\text{3×3}} 表示3×3组卷积。

    改进后的DLL-YOLOv5su桃树缩叶病检测整体网络结构中,LAWD模块替换了骨干网络后三层CBS模块和颈部网络的两层下采样卷积。

    本文研究中针对桃树缩叶病数据集的训练和测试的试验平台环境为:操作系统为64位Windows11,处理器型号为12 th Gen Inter(R) Core(TM)i5-12400F 2.5 GHz,显卡型号为Nvidia GeForce RTX 4060Ti(8G),内存为32GB(3200 MHz),深度学习框架为PyTorch 2.0, CUDA版本为11.7,编程平台为PyCharm,编程语言为Python3.11。

    本文所有试验均在相同环境下进行,选择YOLOv5su作为桃树缩叶病检测的原始模型,为避免过高的学习率导致训练精度震荡,模型收敛困难,超参数文件选择low版本。训练采用的批量大小设置为8,初始学习率为0.01,动量设置为0.937,训练图片尺寸为640×640像素,训练轮数设置为300轮。

    本文主要以准确率(precision, P)、召回率(recall, R)、平均检测精度(average precision, AP),平均检测精度均值(mean average precision,mAP)、每秒帧数(frames per second,FPS)和权重大小(MB)作为评价指标。

    YOLOv5su是在YOLOv5s基础上提出的,为了选择合适的基线模型,本文对部分版本进行了对比试验,试验结果如表1所示。

    表  1  基线模型对比结果
    Table  1.  Baseline model comparison results
    模型
    Models
    准确率
    Precision P/%
    召回率
    Recall R/%
    平均精度均值
    Mean average precision/%
    权重大小
    Model size/MB
    帧速率
    Frame per second
    FPS/(帧∙ {\mathrm{s}}^{-1} )
    平均精度
    Average precision AP/%
    mAP50 mAP50~95 早期
    Early
    中期
    Mid
    后期
    End
    YOLOv5nu 75.3 68.9 73.8 44.3 7.6 105.7 75.1 76.2 70.1
    YOLOv5su 76.5 70.7 76.1 46.8 18.5 119.7 78.3 75.3 74.6
    YOLOv5 mu 77.1 72.4 78.6 48.0 50.5 78.0 79.9 80.2 75.6
    YOLOv5s 74.4 69.6 73.0 44.8 14.5 83.0 78.0 70.1 70.7
    注: mAP50为IoU阈值为0.5时的平均精度均值,mAP50~95为IoU阈值为0.5到0.95,步长为0.05时mAP的平均值。
    Note: mAP50 is the average precision when the intersection over union threshold is 0.5, mAP50~95 is the average value of mAP for IoU thresholds of 0.5 to 0.95 in steps of 0.05.
    下载: 导出CSV 
    | 显示表格

    表1结果可以看出,YOLOv5su模型检测效果优于YOLOv5nu,模型权重大小为18.5 MB,远小于YOLOv5 mu的50.5 MB。对比YOLOv5s,YOLOv5su在准确率、召回率和平均检测精度方面均优,模型增加了4 MB。因此,引入无锚框的分体式头部结构有效地提升了模型性能。综合以上因素,本文选择YOLOv5su作为基线模型。

    为了验证本文改进和提出的方法对于桃树缩叶病检测是否有效,以及评估各个方法所达到的效果,设计了消融试验。试验结果如表2所示,模型改进训练过程参数如图7所示。

    表  2  消融试验结果
    Table  2.  Ablation test results
    模型
    Models
    准确率
    Precision P/%
    召回率
    Recall R/%
    平均精度均值
    Mean average precision/%
    权重大小
    Model size/MB
    帧速率
    Frame per second
    FPS/(帧∙ {\mathrm{s}}^{-1} )
    平均精度
    Average precision AP/%
    mAP50 mAP50~95 早期
    Early
    中期
    Mid
    后期
    End
    YOLOv5su 76.5 70.7 76.1 46.8 18.5 119.7 78.3 75.3 74.6
    YOLOv5su+C3-DA 76.7 74.2 78.4 48.3 19.1 96.1 79.9 80.4 75.1
    YOLOv5su+LSKA 79.0 71.0 78.4 49.0 20.7 114.8 80.4 79.2 75.7
    YOLOv5su+LAWD 76.9 69.7 77.8 47.8 15.0 96.5 79.1 78.5 75.9
    DLL-YOLOv5su 80.7 73.1 80.4 50.4 17.6 83.0 80.6 81.5 79.1
    下载: 导出CSV 
    | 显示表格
    图  7  改进模型性能对比曲线
    Figure  7.  Performance comparison curve of the improved model

    由消融试验结果可以看出各个改进方法的效果:改进C3-DA模块之后,模型相较于初始模型召回率提高了3.5个百分点,平均精度均值mAP50提高了2.3个百分点,但是模型权重增大了0.6 MB;引入LSKA模块之后,模型在准确率和平均精度mAP50方面分别提高了2.5和2.3个百分点,模型权重增加了2.2 MB;使用LAWD替换部分CBS模块之后的模型召回率降低了1个百分点,但准确率和平均精度均值分别上升了0.4和1.7个百分点,模型大小减少了18.9个百分点。

    图7可以看出,改进后的DLL-YOLOv5su模型收敛速度稍慢,但在训练200轮左右时效果已经超越初始版本的YOLOv5su。最终本研究提出的DLL-YOLOv5su桃树缩叶病检测模型准确率、召回率和平均精度均值分别达到了80.7%、73.1%和80.4%。相比原YOLOv5su模型分别提升了4.2、2.4和4.3个百分点。检测速度为83.0帧/s,可以满足桃树缩叶病实时检测的要求。

    为验证本文提出模型的有效性,将其他主流目标检测模型与改进后的模型进行对比试验。所有模型均在相同试验环境下,对相同数据集进行试验,试验结果如表3所示。

    表  3  对比试验结果
    Table  3.  Contrast test results
    模型
    Models
    准确率
    Precision P/%
    召回率
    Recall R/%
    mAP50/% 权重大小
    Model size/MB
    平均精度
    Average precision AP/%
    早期
    Early
    中期
    Mid
    后期
    End
    YOLOv5su 76.5 70.7 76.1 18.5 78.3 75.3 74.6
    Faster R-CNN - - 51.9 110.9 31.5 62.2 62.1
    YOLOv3-tiny 73.9 67.2 68.6 17.5 69.3 66.2 70.4
    YOLOv7 80.1 74.6 78.3 74.8 82.4 78.5 74.0
    YOLOv8s 80.0 70.6 76.3 22.5 78.0 77.9 73.1
    DLL-YOLOv5su 80.7 73.1 80.4 17.6 80.6 81.5 79.1
    下载: 导出CSV 
    | 显示表格

    表3可以看出,YOLOv5su原始模型对比YOLOv7和YOLOv8,在桃树缩叶病识别平均精度均值接近的情况下,权重最小,证明了本文选择其作为基线模型的可行性。在对原始模型进行针对桃树缩叶病特征识别进行改进之后,改进的DLL-YOLOv5su模型平均精度均值mAP50相较于Faster R-CNN、YOLOv3-tiny、YOLOv7和YOLOv8分别高出28.5、11.8、2.1和4.1个百分点。召回率为73.1%,仅低于YOLOv7模型1.5个百分点,而准确率80.7%高于YOLOv7,因为对与目标检测任务,召回率与准确率难以同时兼顾。同时改进的DLL-YOLOv5su模型权重为17.6MB仅为YOLOv7的23.5%。

    为进一步验证DLL-YOLOv5su对桃树缩叶病检测的效果,随机选择测试集图像对模型测试,测试结果如图8所示。从图中可以看出,对于桃树缩叶病早期的小目标检测任务,改进的DLL-YOLOv5su模型检测精度更高。Faster R-CNN和YOLOv7均出现了不同程度的误检;对于模糊目标,仅有改进的模型检测到了两个正确的目标,其余模型均仅检出一个目标或出现误检;对于遮挡的目标识别,YOLOv5su模型将桃树枝干误识别为缩叶病目标,而改进的DLL-YOLOv5su模型识别准确率高,且误检率低。

    图  8  不同环境下各模型对桃树缩叶病识别效果
    注: 绿色方框处为正确检出,红色圆圈处为漏检,蓝色圆圈处为误检,黄色圆圈处为重复检出。
    Figure  8.  Recognition effect of each model on peach leaf curl disease in different environments
    Note: The green boxes are correctly checked out, the red circle is missed, the blue circle is false, and the yellow circle is double-checked.

    试验结果表明,本文提出的DLL-YOLOv5su模型相较于主流目标检测模型,对不同发生时期的桃树缩叶病具有更好的识别效果;相较于文献[12]中的桃树病害检测模型,对桃树缩叶病的检测精度高出了3.82个百分点,且能够反映桃树缩叶病的发生情况。DLL-YOLOv5su模型能够为桃树缩叶病的针对性措施提供技术支持。

    本文通过对自然环境下采集的桃树缩叶病图像进行增强处理,搭建了桃树缩叶病数据集,提出了一种基于改进的YOLOv5su网络的目标检测模型DLL-YOLOv5su以提高桃树缩叶病在自然环境中的检测准确率。研究结果如下:

    1)本文以YOLOv5su为基础网络,在其骨干网络的中引入可变形自注意力机制和可分离大内核注意模块,并在感受野卷积的基础上提出了权重自适应下采样模块替换网络中部分卷积模块,实现了在提高识别准确率和平均精度的前提下,轻量化模型。改进的DLL-YOLOv5su模型的准确率、召回率和平均精度均值mAP50分别达到了80.7%、73.1%和80.4%,相比初始模型分别提高了4.2、2.4和4.3个百分点。改进模型在桃树缩叶病3个发生时期的检测效果均高于初始模型,说明了改进方法的有效性。

    2)通过对比试验与Faster R-CNN、YOLOv3-tiny、YOLOv7、YOLOv8s等主流目标检测模型进行比较,改进的DLL-YOLOv5su模型的准确率最高,平均精度均值mAP50分别高出28.5、11.8、2.1、4.1个百分点。与桃树病害检测模型相比,DLL-YOLOv5su精度也更高,进一步展示了DLL-YOLOv5su桃树缩叶病检测模型的优势。DLL-YOLOv5su模型能够识别自然环境下桃树缩叶病发生的3个阶段,对早期目标的识别效果提升可以更好地达到预警的效果,对于模糊和遮挡目标也具备一定推理能力。

  • 图  1   河北省耕地多功能空间分布图

    Figure  1.   Spatial distribution of multifunctional cultivated land in Hebei Province

    图  2   基于耕地功能四级分区的Spearman秩相关系数走向

    注:A,生产-生态功能;B,生产-生活功能;C,生态-生活功能。

    Figure  2.   Spearman correlation coefficient trend based on the four-level zoning of cultivated land function

    Note: A, Production-Ecological function; B, Production-Living function; C, Ecological-Living function.

    图  3   河北省耕地多功能双变量局部空间自相关集聚图

    Figure  3.   Bivariate Local Moran’s I gradation map of multifunctional cultivated land in Hebei Province

    图  4   河北省耕地多功能分区

    Figure  4.   Multifunctional zoning of cultivated land in Hebei Province

    图  5   河北省耕地功能核心区

    Figure  5.   Cultivated land function core area in Hebei Province

    表  1   河北省耕地多功能评价指标体系

    Table  1   Multifunctional evaluation index system of cultivated land in Hebei Province

    功能
    Function
    子功能
    Subfunction
    指标
    Index
    单位
    Unit
    性质
    Attribute
    指标解释
    Index description
    生产功能
    Production function
    粮食生产功能 粮食生产水平 t·hm−2 + 粮食产量/粮食作物播种面积,表征耕地粮食生产能力
    蔬菜生产功能 蔬菜生产水平 t·hm−2 + 蔬菜产量/蔬菜播种面积,表征耕地蔬菜生产能力
    规模保障功能 土地垦殖率 % + 耕地面积/土地总面积,表征区域耕地相对规模大小
    生态功能
    Ecological function
    生态多样性功能 农田生态系统多样性指数 % + -\sum {b}_{i}\ln{b}_{i} {b}_{i} 为各品种作物播种面积与作物总播种面积之比
    环境承载功能 化肥施用强度 t·hm−2 - 化肥施用量(折纯量)/耕地面积,反映耕地利用过程造成的化学污染
    景观美学功能 耕地景观聚集度 - + 通过fragstats软件进行计算,反映耕地斑块的连片聚集程度
    生活功能
    Living function
    就业保障功能 劳动力承载量 % + 农林牧渔业从业人员数量/乡村从业人员数量,反映耕地就业保障能力
    经济保障功能 农业经济占比 % + 农业增加值/地区生产总值,反映农业经济保障能力
    粮食安全保障功能 人均自产粮食保障率 % + 粮食产量/(常住人口数 \times 400 kg),表征耕地粮食安全保障能力
    注:农田生态系统多样性指数综合考虑河北省农作物种植类别,选取包括稻谷、小麦、玉米、大豆4种主要粮食作物和油料、棉花、糖料、蔬菜4种主要经济作物的农作物播种面积。
    Note: The farmland ecosystem diversity index comprehensively considers the crop categories in Hebei Province, and selects the sown area of four major food crops including rice, wheat, corn and soybean and four main cash crops of oil seeds, cotton, sugar and vegetables.
    下载: 导出CSV

    表  2   河北省耕地多功能评价指标权重

    Table  2   The weight of the multifunctional evaluation index of cultivated land in Hebei Province

    功能
    Function
    权重
    Weight
    指标
    Index
    指标权重Index weight
    熵权法
    Entropy
    weight
    method
    层次
    分析法
    Analytic
    hierarchy
    process
    平均
    权重
    Average
    weight
    生产功能
    Production function
    0.4600 粮食生产水平 0.1652 0.5396 0.3524
    蔬菜生产水平 0.1667 0.2970 0.2319
    耕地垦殖率 0.6681 0.1634 0.4158
    生态功能
    Ecological function
    0.3189 农田生态系统多样性指数 0.4213 0.4000 0.4107
    化肥施用强度 0.1257 0.2000 0.1629
    耕地景观聚集度 0.4530 0.4000 0.4265
    生活功能
    Living function
    0.2211 劳动力承载量 0.1170 0.3750 0.2460
    农业经济占比 0.5170 0.2500 0.3835
    人均自产粮食保障率 0.3660 0.3750 0.3705
    下载: 导出CSV

    表  3   河北省耕地各项功能间的Spearman秩相关系数

    Table  3   Spearman correlation coefficients among cultivated land functions in Hebei Province

    年份
    Year
    生产-生态
    Production-Ecological
    生产-生活
    Production-Living
    生态-生活
    Ecological-Living
    20100.611**0.368**0.274**
    20150.727**0.397**0.384**
    20200.755**0.469**0.336**
    注:**表示在0.01级别(双尾)的水平上显著相关。
    Note: ** indicate significant correlation at the level of 0.01 (double tails).
    下载: 导出CSV

    表  4   河北省耕地子功能之间的Spearman秩相关系数

    Table  4   Spearman correlation coefficients among cultivated land subfunctions in Hebei Province

    项目
    Item
    粮食生产
    Grain production
    蔬菜生产
    Vegetable production
    规模保障
    Large-scale production
    就业保障
    Job security
    经济保障
    Financial security
    粮食安全保障
    Grain security
    环境承载
    Environment-
    bearing
    景观美学
    Landscape aesthetics
    生态多样性
    Ecological diversity
    粮食生产 1 0.343** 0.560** -0.309** 0.168* 0.607** -0.566** 0.516** 0.238**
    蔬菜生产 0.113 1 0.104 -0.306** -0.002 0.181* -0.341** 0.131 0.176*
    规模保障 0.590** 0.099 1 -0.297** 0.425** 0.611** -0.327** 0.808** 0.342**
    就业保障 -0.174* 0.049 -0.088 1 0.312** -0.291** 0.351** -0.267** -0.303**
    经济保障 0.129 0.353** 0.367** 0.455** 1 0.316** -0.072 0.368** 0.088
    粮食安全保障 0.450** 0.111 0.793** 0.195* 0.594** 1 -0.298** 0.624** 0.190*
    环境承载 -0.498** -0.155* -0.243** 0.203** -0.053 -0.061 1 -0.318** -0.203*
    景观美学 0.497** -0.017 0.856** -0.112 0.266** 0.679** -0.299** 1 0.358**
    生态多样性 0.308** 0.215** 0.355** -0.141 0.247** 0.163* -0.388** 0.298** 1
    注:**、*分别表示在0.01、0.05(双尾)的水平上显著相关,白色背景区域和灰色背景区分别为2010、2020年耕地子功能相关系数值。
    Note: ** and * indicate significant correlation at the level of 0.01 and 0.05 (double tails), the white background area and the gray background area are the value of the sub-functional correlation coefficient of cultivated land in 2010, 2020, respectively.
    下载: 导出CSV

    表  5   河北省耕地多功能分区标准

    Table  5   Zoning standard of multifunctional cultivated land

    耕地功能分区
    Functional zoning of
    cultivated land
    耕地功能
    Cultivated land functions
    耕地功能核心区
    Cultivated land
    function core area
    耕地功能空间相关状况
    Spatial correlation of cultivated land functions
    生产功能
    Production function
    生态功能
    Ecological function
    生活功能
    Living function
    生产-生态
    Production-ecological
    生产-生活
    Production-living
    生态-生活
    Living-ecological
    综合功能优势区 生产-生态协同核心区 HH - -
    综合功能优势区 生产-生活协同核心区 - HH -
    综合功能优势区 生态-生活协同核心区 - - HH
    生产-生态功能优势区 生产-生态协同核心区 HH - -
    生产-生活功能优势区 生产-生活协同核心区 - HH -
    生态-生活功能优势区 生态-生活协同核心区 - - HH
    生产功能优势区 非核心区 - - -
    生态功能优势区 非核心区 - - -
    生活功能优势区 非核心区 - - -
    功能低值区 非核心区 - - -
    注:在耕地功能核心区划分中,需要满足耕地功能强度和空间相关要求,同时存在两类及以上核心区则为综合功能协同核心区。
    Note: In the division of cultivated land functional core areas, it is necessary to meet the requirements for the intensity of cultivated land functions and spatial relevance. At the same time, if there are two or more types of core areas, it is considered a comprehensive functional synergy core area.
    下载: 导出CSV
  • [1] 范业婷,金晓斌,项晓敏,等. 苏南地区耕地多功能评价与空间特征分析[J]. 资源科学,2018,40(5):980-992.

    FAN Yeting, JIN Xiaobin, XIANG Xiaomin, et al. Evaluation and spatial characteristics of arable land multifunction in southern Jiangsu[J]. Resources Science, 2018, 40(5): 980-992. (in Chinese with English abstract)

    [2] 施园园. 大都市区耕地多功能评价[D]. 北京:中国地质大学(北京),2016. SHI Yuanyuan. Assessment of Arable Land Multifunction in Metropolitan Areas: Case Study in Beijing[D]. Beijing: China University of Geosciences, 2016. (in Chinese with English abstract
    [3]

    GROOT R D. Function-analysis and valuation as a tool to assess land use conflicts in planning for sustainable, multi-functional landscapes[J]. Landscape and Urban Planning, 2005, 75(3): 175-186.

    [4]

    GROOT R S D, WILSON M A, BOUMANS R MJ. A typology for the classification, description and valuation of ecosystem functions, goods and services[J]. Ecological Economics, 2002, 41(3): 393-408. doi: 10.1016/S0921-8009(02)00089-7

    [5]

    ASSESSMENT M E. Ecosystems and Human Well-being[M]. Washington DC: Island Press, 2005.

    [6] 刘超,许月卿,孙丕苓,等. 土地利用多功能性研究进展与展望[J]. 地理科学进展,2016,35(9):1087-1099. doi: 10.18306/dlkxjz.2016.09.004

    LIU Chao, XU Yueqing, SUN Piling, et al. Progress and prospects of multi-functionality of land use research[J]. Progress in Geography, 2016, 35(9): 1087-1099. (in Chinese with English abstract) doi: 10.18306/dlkxjz.2016.09.004

    [7] 宋小青,欧阳竹. 耕地多功能内涵及其对耕地保护的启示[J]. 地理科学进展,2012,31(7):859-868.

    SONG Xiaoqing, OU Yangzhu. Connotation of multifunctional cultivated land and its implications for cultivated land protection[J]. Progress in Geography, 2012, 31(7): 859-868. (in Chinese with English abstract)

    [8]

    JIN X B, ZHANG Z H, WU X W, et al. Co-ordination of land exploitation, exploitable farmland reserves and national planning in China[J]. Land Use Policy, 2016, 57: 682-693. doi: 10.1016/j.landusepol.2016.06.036

    [9] 李怡,方斌,李裕瑞,等. 城镇化进程中耕地多功能权衡/协同关系演变及其驱动机制[J]. 农业工程学报,2022,38(8):244-524.

    LI Yi, FANG Bin, LI Yurui, et al. Trade-off and synergy evolution of farmland functions and its dynamic mechanism in the process of urbanization[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2022, 38(8): 244-254. (in Chinese with English abstract)

    [10] 钱凤魁,迟艳茹,徐欢,等. 2006-2020年沈阳市耕地多功能权衡协同关系演变分析[J]. 中国土地科学,2022,36(10):31-41.

    QIAN Fengkui, CHI Yanru, XU Huan, et al. Study on evolution of trade-off and synergy relationship of multifunctional cultivated land from 2006 to 2020: A case of Shenyang City[J]. China Land Science, 2022, 36(10): 31-41. (in Chinese with English abstract)

    [11] 王成,彭清,唐宁,等. 2005~2015年耕地多功能时空演变及其协同与权衡研究——以重庆市沙坪坝区为例[J]. 地理科学,2018,38(4):590-599.

    WANG Cheng, PENG Qing, TANG Ning, et al. Spatio-temporal evolution and the synergy and trade-off relationship of cultivated land multi-function in 2005-2015: A case of Shapingba District, Chongqing City[J]. Scientia Geographica Sinica, 2018, 38(4): 590-599. (in Chinese with English abstract)

    [12] 卫新东,林良国,罗平平,等. 耕地多功能耦合协调发展时空格局与驱动力分析[J]. 农业工程学报,2022,38(4):260-269.

    WEI Xindong, LIN Liangguo, LUO Pingping, et al. Spatiotemporal pattern and driving force analysis of multi-functional coupling coordinated development of cultivated land[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2022, 38(4): 260-269. (in Chinese with English abstract)

    [13]

    FISHER B, TURNER R K, MORLING P. Defining and classifying ecosystem services for decision making[J]. Ecological Economics, 2009, 68(3): 643-653. doi: 10.1016/j.ecolecon.2008.09.014

    [14]

    FU B J, SU C H, WEI Y P, et al. Double counting in ecosystem services valuation: Causes and countermeasures[J]. Ecological Research, 2011, 26(1): 1-14. doi: 10.1007/s11284-010-0766-3

    [15] 陈丽. 黄淮海平原耕地多功能效应及保护补偿研究[D]. 北京:中国农业大学,2016. CHEN Li. Study on Multifunction Effect of Cultivated Land and Compensation for Cultivated Land Protection in Huang-Huai-Hai Plain[D]. Beijing: China Agricultural University, 2016. (in Chinese with English abstract
    [16] 高星,宋昭颖,李晨曦,等. 城乡梯度下的耕地多功能价值空间分异特征[J]. 农业工程学报,2021,37(16):251-259.

    GAO Xing, SONG Zhaoying, LI Chenxi, et al. Spatial differentiation characteristics of cultivated land multifunctional value under urban-rural gradient[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2021, 37(16): 251-259. (in Chinese with English abstract)

    [17] 熊昌盛,张永蕾,王雅娟,等. 中国耕地多功能评价及分区管控[J]. 中国土地科学,2021,35(10):104-114.

    XIONG Changsheng, ZHANG Yonglei, WANG Yajuan, et al. Multi-function evaluation and zoning control of cultivated land in China[J]. China Land Science, 2021, 35(10): 104-114. (in Chinese with English abstract)

    [18] 祝汉收,翟俊,侯鹏,等. 生态系统服务权衡与协同视角下的重点生态功能区保护特征[J]. 地理学报,2022,77(5):1275-1288.

    ZHU Hanshou, ZHAI Jun, HOU Peng, et al. The protection characteristics of key ecological functional zones from the perspective of ecosystem service trade-off and synergy[J]. Acta Geographica Sinica, 2022, 77(5): 1275-1288. (in Chinese with English abstract)

    [19] 方莹,王静,孔雪松,等. 耕地利用多功能权衡关系测度与分区优化——以河南省为例[J]. 中国土地科学,2018,32(11):57-64.

    FANG Ying, WANG Jing, KONG Xuesong, et al. Trade-off relation measurement and zoning optimization of multi-functionality of cultivated land use: A case study of Henan Province[J]. China Land Science, 2018, 32(11): 57-64. (in Chinese with English abstract)

    [20] 牛丽楠,邵全琴,宁佳,等. 西部地区生态状况变化及生态系统服务权衡与协同[J]. 地理学报,2022,77(1):182-195.

    NIU Linan, SHAO Quanqin, NING Jia, et al. Ecological changes and the tradeoff and synergy of ecosystem services in western China[J]. Acta Geographica Sinica, 2022, 77(1): 182-195. (in Chinese with English abstract)

    [21] 张春悦,白永平,杨雪荻,等. 多情景模拟下宁夏平原生态系统服务簇识别研究[J]. 地理研究,2022,41(12):3364-3382.

    ZHANG Chunyue, BAI Yongping, YANG Xuedi, et al. Identification of ecosystem service bundles in Ningxia Plain under multi-scenario simulation[J]. Geographical Research, 2022, 41(12): 3364-3382. (in Chinese with English abstract)

    [22] 方露露,许德华,王伦澈,等. 长江、黄河流域生态系统服务变化及权衡协同关系研究[J]. 地理研究,2021,40(3):821-838.

    FANG Lulu, XU Dehua, WANG Lunche, et al. The study of ecosystem services and the comparison of trade-off and synergy in Yangtze River Basin and Yellow River Basin[J]. Geographical Research, 2021, 40(3): 821-838. (in Chinese with English abstract)

    [23] 刘玉,郜允兵,潘瑜春,等. 基于多源数据的乡村功能空间特征及其权衡协同关系度量[J]. 地理研究,2021,40(7):2036-2050.

    LIU Yu, GAO Yunbing, PAN Yuchun, et al. Spatial differentiation characteristics and trade-off/synergy relationships of rural multi-functions based on multi-source data[J]. Geographical Research. 2021, 40(7): 2036-2050. (in Chinese with English abstract)

    [24] 董鹏宇,赵华甫. 耕地多功能权衡与协同关系研究——以上海市青浦区为例[J]. 长江流域资源与环境,2019,28(2):368-375.

    DONG Pengyu, ZHAO Huafu. Study on trade-off and synergy relationship of cultivated land multifunction: A case of Qingpu District, Shanghai[J]. Resources and Environment in the Yangtze Basin, 2019, 28(2): 368-375. (in Chinese with English abstract)

    [25] 聂名萱,黄思华,濮励杰,等. 快速城镇化地区生态系统服务的时空动态及权衡与协同分析——以苏锡常地区为例[J]. 长江流域资源与环境,2021,30(5):1088-1099.

    NIE Mingxuan, HUANG Sihua, PU Lijie, et al. Spatial and temporal dynamics and trade-offs and synergies analysis of ecosystem services in rapidly urbanizing areas: A case study of the Su-Xi-Chang region[J]. Resources and Environment in the Yangtze Basin, 2021, 30(5): 1088-1099. (in Chinese with English abstract)

    [26] 齐麟,张月,许东,等. 东北森林屏障带生态系统服务权衡与协同关系[J]. 生态学杂志,2021,40(11):3401-3411.

    QI Lin, ZHANG Yue, XU Dong, et al. Trade-offs and synergies of ecosystem services in forest barrier belt of Northeast China[J]. Chinese Journal of Ecology, 2021, 40(11): 3401-3411. (in Chinese with English abstract)

    [27]

    YANG J, HUANG X. The 30 m annual land cover dataset and its dynamics in China from 1990 to 2019[J]. Earth System Science Data, 2021, 13(8): 3907-3925. doi: 10.5194/essd-13-3907-2021

    [28] 安悦,谭雪兰,李印齐,等. 洞庭湖地区耕地功能时空演变特征及影响因素研究[J]. 地理科学,2022,42(7):1272-1282.

    AN Yue, TAN Xuelan, LI Yinqi, et al. Spatio-temporal evolution characteristics and influencing factors of cultivated land function in the Dongting Lake aera[J]. Scientia Geographica Sinica, 2022, 42(7): 1272-1282. (in Chinese with English abstract)

    [29] 施园园,赵华甫,郧文聚,等. 北京市耕地多功能空间分异及其社会经济协调模式解释[J]. 资源科学,2015,37(2):247-257.

    SHI Yuanyuan, ZHAO Huafu, YUN Wenju, et al. Analysis on spatial differentiation of arable land multifunction and socio-economic coordination model in Beijing[J]. Resources Science, 2015, 37(2): 247-257. (in Chinese with English abstract)

    [30] 杨雪,谈明洪. 近年来北京市耕地多功能演变及其关联性[J]. 自然资源学报,2014,29(5):733-743.

    YANG Xue, TAN Minghong. Changes and relationships of arable land functions in Beijing in recent years[J]. Journal of Natural Resources, 2014, 29(5): 733-743. (in Chinese with English abstract)

    [31] 朱庆莹,胡伟艳,赵志尚. 耕地多功能权衡与协同时空格局的动态分析——以湖北省为例[J]. 经济地理,2018,38(7):143-153.

    ZHU Qingying, HU Weiyan, ZHAO Zhishang. Dynamic analysis on spatial-temporal pattern of trade-offs and synergies of multifunctional cultivated land—Evidence from Hubei Province[J]. Economic Geography, 2018, 38(7): 143-153. (in Chinese with English abstract)

    [32] 路昌,耿昊,刘蕊,等. 山东省耕地多功能耦合协调度时空格局分析[J]. 农业机械学报,2021,52(3):265-274.

    LU Chang, GENG Hao, LIU Rui, et al. Spatial-temporal pattern analysis of cultivated land multi-functional coupling coordination degree in Shandong Province[J]. Transactions of the Chinese Society for Agricultural Machinery, 2021, 52(3): 265-274. (in Chinese with English abstract)

    [33] 张一达,刘学录,任君,等. 基于耕地多功能权衡与协同分析的耕地利用转型研究——以北京市为例[J]. 中国农业资源与区划,2020,41(6):25-33.

    ZHANG Yida, LIU Xuelu, REN Jun, et al. Study on cultivated land use transformation based on multifunctional trade-offs and collaborative analysis of arable land: taking Beijing as an example[J]. Chinese Journal of Agricultural Resources and Regional Planning, 2020, 41(6): 25-33. (in Chinese with English abstract)

    [34]

    ANSELIN L, SYABRI I, KHO Y. GeoDa : An introduction to spatial data analysis[J]. Geographical Analysis, 2006, 38(1): 5-22.

    [35] 文高辉,王夏玙,胡贤辉,等. 耕地规模、细碎化对化肥面源污染的影响[J]. 中国生态农业学报(中英文),2023,31(7):1081-1093.

    WEN Gaohui, WANG Xiayu, HU Xianhui, et al. Impact of cultivated land scale and fragmentation on chemical fertilizer non-point source pollution[J]. Chinese Journal of Eco-Agriculture, 2023, 31(7): 1081-1093. (in Chinese with English abstract)

  • 期刊类型引用(3)

    1. 杨麒. 北方严寒地区农田旱改水工程施工技术探析. 建筑施工. 2024(04): 502-506 . 百度学术
    2. 贾新乐,石舟,李锐,张国海,耿端阳,兰玉彬,王伯龙. 丘陵山区小型农机底盘自动调平系统设计与试验. 农业机械学报. 2024(S1): 108-115 . 百度学术
    3. 孔德航,张学东,崔巍,吴海华,孙星,王志伟,王春雷,宁义超. 基于模糊PID的顶夹式取苗装置苗盘定位控制方法. 农业机械学报. 2024(S1): 207-216+229 . 百度学术

    其他类型引用(2)

图(5)  /  表(5)
计量
  • 文章访问数:  179
  • HTML全文浏览量:  23
  • PDF下载量:  91
  • 被引次数: 5
出版历程
  • 收稿日期:  2024-03-06
  • 修回日期:  2024-04-27
  • 刊出日期:  2024-07-29

目录

/

返回文章
返回