月尺度农作物提取中GF-1 WFV纹理特征的应用及分析
Application and analyses of texture features based on GF-1 WFV images in monthly information extraction of crops
通讯作者: 赵红莉(1973-)女,教授级高级工程师,主要研究方向为水资源信息化。Email:zhaohl@iwhr.com。
责任编辑: 李瑜
收稿日期: 2020-10-21 修回日期: 2021-01-19
基金资助: |
|
Received: 2020-10-21 Revised: 2021-01-19
作者简介 About authors
王 镕(1993-)女,硕士,主要研究方向为农业与水资源遥感应用。Email:
农作物种植结构包含农作物种类、数量结构和空间分布特征等信息,是农业科学管理的基础。在不考虑农作物时间序列最佳窗口期的前提下,以石津灌区为研究区,基于高分一号(GF-1)WFV影像计算并分析纹理特征在农作物分类识别中的能力。并在纹理特征分类效果相对较差的时相内引入植被指数,从而弥补纹理在农作物表达上的缺陷。经过对比各组分类结果,可以发现: 在作物结构明显的4,8月份,单独纹理特征的分类精度可以达到80%以上,但是在5,6,7,9月等农作物最复杂的时间段内,分类精度仍低于80%。将植被指数与纹理特征组合后,这4个月份的分类结果有了很大改善,总体分类精度均大于80%,基本满足农业动态监测的需求; 与单独纹理相比,精度提高2.27%~9.75%, Kappa系数提高0.02~0.16; 利用夏玉米的验证样本进行验证,识别精度可以达到98%,识别效果相对完整,破碎程度达到最小化,与其他类别区分度也达到了最优。同时也证明了GF-1WFV纹理特征在农作物种植结构提取中的可用性,尤其是在作物结构相对明显的月份内,可以为影像的农作物提取提供一些有效的信息。
关键词:
The crop planting structure consists of information such as crop species, quantity structure, and spatial distribution characteristics, and it serves as the basis for agricultural scientific management. Taking the Shijin irrigation area, Hebei Province as the study area and on the premise of not considering the optimal window period of crop time series, this study calculates and analyzes the ability of texture features in crop classification and identification based on GF-1WFV images. Meanwhile, the vegetation index is introduced into the time phase in which the classification effects based on texture features are poor, in order to make up for the defects of texture in the expression of crops. According to the comparison of the classification results of various groups, the classification accuracy of individual texture features reached greater than 80% in April and August when the crop structure is obvious but was still less than 80% in May, June, July, and September when crops are the most complex. After combining the texture features with the vegetation index, the classification results of the crops in these four months were greatly improved. In detail, the overall classification accuracy was greater than 80%, which basically meets the need for agricultural dynamic monitoring. Meanwhile, the accuracy was improved by 2.27%~9.75 % and the Kappa coefficient was increased by 0.02~0.16 compared to the individual texture features. As verified using summer maize samples, the recognition accuracy reached up to 98%, the recognition effects were relatively complete, the fragmentation degree was the minimum, and the optimal discrimination from other crop categories was achieved. Meanwhile, it also proved that the texture features based on GF-1WFV images can be applied to the extraction of the crop planting structure, especially in the months when the crop structure is relatively obvious, and they can provide some effective information for the information extraction of crops base on images.
Keywords:
本文引用格式
王镕, 赵红莉, 蒋云钟, 何毅, 段浩.
WANG Rong, ZHAO Hongli, JIANG Yunzhong, HE Yi, DUAN Hao.
0 引言
农作物种植结构反映了人类农业生产在空间范围内利用农业生产资源的状况[1],是农作物种类、数量结构和空间分布特征等的信息组合,也是农业科学灌溉的基础。目前基于遥感的农作物种植结构提取多采用影像的光谱特征与纹理特征,光谱特征一般用于中低分辨率数据,如TM,SPOT,HJ等或NOAA,MODIS等[2]。中低分辨率的遥感数据获取来源较多,且影像覆盖范围广,但是一般适合大面积单一作物的识别[3]。针对复杂的地物类型,仅光谱特征难以准确的区分各类作物,最终影响其分类精度。如部分研究利用Landsat8[4], HJ-1A/B[5], MODIS-EVI[6]影像光谱信息构建小麦、水稻等作物的长时间序列变化特征,并提取相应地物的空间分布信息,精度可达到85%以上; 而纹理特征一般用于高分辨率的遥感影像识别过程,如QuickBird, SPOT[7], SAR, SuperView-1等。高分辨率影像拥有丰富的纹理、色调、形状和几何结构等特征信息,地物内部细节信息明显,边缘突出,具有较高的分辨精度与目标识别可靠性。如研究人员利用灰度共生矩阵[8]计算QuickBird,SPOT,SAR[9]等影像的纹理特征,实现植被类型的分布并实现长期的监测,精度均在90%以上。但是由于重访周期长、获取最佳时相困难、且价格高昂等缺点,难以满足农作物动态管理监测的需求。中高分辨率高分一号(GF-1)作为我国自主研发的卫星,因其重访周期短,分辨率高等特点,自发射以来,已开始应用于我国林业、农业、自然灾害、土地利用、定量化应用等众多领域。如黄健熙等[10]利用大豆等作物的最佳生长期与影像光谱信息的之间的变化规律实现了集中作物的分类,精度最高达到85%; 王利民等[11]从不同的分类方法入手实现冬小麦等大面积作物的提取,进而对比面向对象、决策树等分类方法的提取精度[12]。同时,也有部分研究通过多种卫星影像数据的对比,从而分析GF-1在分类中的优势。如欧阳玲[13]等等利用Landsat8,HJ-1A/B与GF-1数据,采用农作物遥感反演方法实现地物空间分布信息的提取,并基于影像的时效性与分类精度等多个方面分析GF-1影像的优势,最终推广国产卫星的应用范围[14]。
1 研究区概况及数据源
1.1 研究区概况
图1
图1
研究区位置示意图
Fig.1
The location of research area and the spectral characteristics of sample
1.2 试验设计与数据收集
1) 试验设计。利用光谱信息识别农作物,需要依赖作物的最佳窗口期,无法实现农作物的动态监测。为实现不依赖作物最佳窗口期的前提下动态获取农作物的变化特征,本文以2019年GF-1WFV多期影像为主,提出利用纹理特征以及光谱辅助纹理等两组实验提取农作物,通过对比农作物空间分布信息与精度评价结果,分析纹理特征在农作物种植结构识别中的能力。
2) 遥感数据。GF-1搭载了4台WFV相机,幅宽可达到800 km,重访周期为4 d,空间分辨率可达到16 m。根据对冬小麦[17]、夏玉米、棉花等作物的生长物候期分析以及本文对纹理特征在农作物识别效应的分析需求,选取2019年4月—2019年9月共计14景 GF-1WFV作物混杂程度不同的影像数据,实现了研究区内6期的全覆盖。对影像进行了辐射校正、大气校正、正射校正以及图像拼接裁剪等预处理。
图2
2 研究方法
2.1 影像纹理特征选取
GLCM通过两点之间的相关性表达了影像灰度在方向、间隔和变化幅度上的空间信息[25],但是在进行纹理分析时,还不能直接作为特征变量用于区别不同地物,需要通过在此基础上计算纹理特征量。常用的特征量一般有8种[26],通过不同角度刻画不同的空间纹理信息。具体如表1。表中,i, j为以影像左上角为坐标原点的各个像元的行列号; p(i, j)为像元的灰度值; μ是p(i, j)的均值; n为像元行列号差值的绝对值; L为影像的行或列的总数。其中考虑到窗口大小、距离d以及方向θ等对GLCM计算的影响,通过多组窗口的提取试验,并结合文献中关于纹理计算的研究,本文将移动窗口确定为5×5,同时默认θ=0,d=1对纹理进行计算。
表1 纹理特征量
Tab.1
统计量 | 计算公式 |
---|---|
Mean(平均值) | - |
Var(方差) | |
Con(对比度) | |
Hom(同质性) | |
Dis (非相似性) | |
Ent(信息熵) | - |
ASM(二阶矩) | |
Cor(相关性) |
2.2 分类方法选择
针对遥感影像的农作物种植结构识别,较为常用的方法有决策树、面向对象、神经网络、随机森林等[27]。如大量研究基于决策树实现研究区内的作物分布统计,但是该分类方法完全依赖人工操作,且需要操作人员具有较强的专业知识以及大量的经验积累; 还有部分研究基于面向对象实现土地利用六大类的提取,结果表明面向对象方法通过合适的分割尺度可以减少分类中的破碎斑块现象,但是无法实现自动化,后期需要大量的人为操作; 神经网络方法也常常用于影像的分类过程,该方法对专业性的要求不高,但是泛化性较差[28]; 而随机森林是一个较为完善的分类方法,由于参数较少、方法简单而被广泛使用。其本质上是对决策树的改进以达到分类精度提高的目的,是通过集成学习的思想将多棵CART决策树集成的机器学习方法,减少了大量的人为干预; 同时该方法也实现了样本与特征两者的随机性,从而避免了训练的过拟合现象,增强了方法的泛化能力[29]。利用随机森林方法可以有效的克服传统方法在分类中不足,具有较高的精确度。
3 基于纹理特征的农作物分类效果
在不考虑农作物时间序列最佳窗口期的影响,检验每期影像独立提取种植结构的精度,从而验证GF-1纹理特征在各期数据中的识别能力。首先利用纹理特征对2019年4—9月的6期影像进行分类,获取研究区内的农作物种植结构,分类结果统计如表2。
表2 纹理特征分类结果统计
Tab.2
日期 | 冬小麦 | 棉花 | 蔬菜 | 夏玉米 | 经济园林 |
---|---|---|---|---|---|
4月18 | 1 675.3 | 90.9 | 687.1 | ||
5月4 | 1 696.9 | 88.9 | 508 | 689.1 | |
6月30 | 90.4 | 562 | 755.8 | ||
7月21 | 105 | 521 | 1 785 | 658.1 | |
8月28 | 101 | 558 | 1 798 | 692.1 | |
9月24 | 129 | 648 | 1 830 | 742.5 |
各时相内的作物面积变化幅度均在2%~35%以内,如冬小麦在2期面积变化为12%,棉花6期相差3%~35%; 夏玉米3期面积相差3%~15%; 蔬菜5期面积相差3%~21%; 经济园林6期面积变化范围在2%~10%。与实际调查的作物面积统计数据相对保持一致,表明纹理可以为影像分类提供有效的信息,具有一定的应用价值。
利用纹理特征的分类结果具体如图3所示。可以发现,冬小麦、夏玉米种植较为集中,地块相对完整,主要集中在东部和南部,两者的空间分布基本一致; 北部和西部经济园林面积相对较大; 蔬菜等主要分布在城区周边; 棉花主要分布在东北部,基本符合实际的种植情况。六期数据中4月和8月2个月的分类结果较好,而其他月份随着作物类别增加与作物生长特征不明显等原因,大大增加了遥感分类的难度。如标注内5月份林地和冬小麦大面积混淆,7月份易混淆的是林地和蔬菜; 这也是由于两月份内几者的生长状态相似,纹理特征差异较小; 6月份影像有大片云层干扰,忽略云层影响,棉花和蔬菜出现大量的错分,这是由于棉花和蔬菜较为分散,且六月份均处于苗期,16 m的分辨率很难表达两者的差异; 而9月份的分类结果破碎斑块明显增加,夏玉米等作物错分现象明显。这4个月的结果也表明了单独纹理特征无法满足各个时相的分类需求。
图3
4 引入光谱特征后的农作物分类效果
根据上述的结论分析,可以得知基于纹理的分类结果在5月、6月、7月、9月等月份存在大量的错分误分现象。考虑到植被指数作物的分布及生长状态呈线性关系,是作物生长状态监测的最佳指标。在纹理特征相对较弱的时相内引入植被指数,最终确定适合不同月份的最佳分类特征。将光谱与纹理组合后的结果与训练样本作为随机森林分类器的输入,利用随机森林的回归分析并训练分类模型,获取组合后的分类结果面积统计如表3所示,各类作物的面积统计与年鉴统计调查结果高度一致。与单一的纹理特征相比,各类作物面积变化幅度相对较小,均保持在150 km2以内。根据实测样本对2种分类结果进行精度评价,结果如表4所示。综合6期分类结果数据,单独纹理特征分类结果总体精度最高可达到为92.95%,Kappa系数为0.91; 将纹理与光谱组合之后的总体分类精度最高可以达到95.22%,对应的Kappa系数为0.93。组合后6期的分类精度均在80%以上,分类结果较为理想; 而单独纹理特征在不同时期出现了低值; 如9月份的纹理特征信息较弱,大面积的林地减少,蔬菜等也出现了不同程度的破碎斑块。总体而言,两者组合后,分类特征增加,使得误分或错分的像元有所减少; 与单独纹理特征的分类精度相比,提高6%~10%。更加验证了多特征组合的分类结果要优于单独纹理的分类结果。
表3 组合后分类结果统计
Tab.3
日期 | 冬小麦 | 棉花 | 蔬菜 | 夏玉米 | 经济园林 |
---|---|---|---|---|---|
4月18 日 | 1 828.1 | 89.3 | 477.4 | ||
5月4日 | 1 853.7 | 83.4 | 603.4 | 532.4 | |
6月30日 | 94.9 | 668.5 | 649.1 | ||
7月21日 | 104 | 626.7 | 1 741 | 643.2 | |
8月28日 | 109 | 713.2 | 1 774 | 590.2 | |
9月24日 | 102 | 746.1 | 1 713 | 601.6 |
表4 分类精度评价表
Tab.4
日期 | 纹理特征 | 特征组合 | ||
---|---|---|---|---|
精度/% | Kappa | 精度/% | Kappa | |
4月18日 | 92.95 | 0.91 | 95.22 | 0.93 |
5月4日 | 81.29 | 0.74 | 85.24 | 0.79 |
6月30日 | 70.54 | 0.59 | 80.29 | 0.75 |
7月21日 | 78.69 | 0.71 | 81.24 | 0.73 |
8月28日 | 83.69 | 0.81 | 87.05 | 0.81 |
9月24日 | 82.91 | 0.79 | 86.39 | 0.80 |
在纹理的基础上引入植被指数生成研究区的主要农作物分布图(图4)。通过图中标注的区域与单独纹理分类结果中对应位置的对比,5月、6月、7月和9月的分类结果中各类别的错分和误分像元明显减少,而4月、8月份两组实验的结果相对一致,表明了光谱的引入大大提高了5月、6月、7月、9月份的分类精度; 尤其是9月份效果更加明显,破碎斑块明显减少。
图4
通过对比两组实验,发现单独纹理特征分类更适合4月、8月2个月份的分类。这是由于4月份作物结构相对简单,各类作物在影像上所呈现的纹理结构相对完整,分类过程中不会造成太大的干扰,更加有利于影像的识别; 而8月份是由于各类作物均处于最佳的生长状态,作物的形状、大小等差异也达到了最大化,在影像中纹理特征的区分度也有了很大的提升。因此,这2个月份内,单独纹理特征不仅可以达到农作物监测的精度要求,而且可以有效地降低分类的时间复杂度,极大程度地减少计算机的存储量。而针对其他时段,使用特征组合的方法要远远优于单独纹理特征的分类精度。这是由于5月、6月、7月、9月这4个月的作物大多处于萌芽阶段,形状及大小极为相近。作物区分难度增加,此时在纹理特征的基础上引入光谱特征,利用光谱的吸收和反射差异,增强不同作物的生长状态信息,从而提高分类特征的区分能力,精确的获取作物的种植结构,实现逐月的农作物动态监测。
5 结论与展望
本文基于GF-1WFV影像提取研究区主要种植农作物的空间分布信息。根据农作物由简单到复杂的生长期变化特点以及特征量的不同选择原则,选取2019年4月—2019年9月共6期影像进行2组实验,并利用实测样本数据进行验证,主要得到以下结论:
1)单独纹理特征可以为WFV影像分类提供许多有效的信息,如4月和8月份,分类精度可以达到80%以上,充分说明了GF-1WFV卫星影像数据的纹理特征在农作物识别过程中的可行性; 但是在5月、6月、7月、9月等农作物复杂的时间段内,分类精度仍低于80%。
2)光谱与纹理组合使得各时相内的分类精度均超过80%,与纹理的分类结果相比,精度提高2.27%~9.75%。
3)通过几组实验对比,在满足农作物监测的精度要求以及影像分类效率等的基础上,4月、8月采用单一纹理特征的识别方法,而5月、6月、7月、9月则采用纹理与光谱组合的提取方法。
GF-1卫星以4 d的周期实现宽幅并且高精度的数据,与SPOT,MODIS,Landsat等卫星相比,满足农业遥感逐月监测的需求,为农业灌溉面积监测以及农作物估产等应用提供了有力的支撑。但是本文在农作物信息提取过程中,也存在着一定的缺陷。比如缺少对特征的筛选。灰度共生矩阵计算的纹理特征量较多,部分特征量之间存在一定的相关性,从而出现数据冗余,增加了分类的时间复杂度,同时也影响了农作物提取的精度; 因此进一步的研究可以围绕特征的优选算法展开,通过动态的获取信息量最优的特征量并参与影像的分类,最终达到缩短分类时间、提高分类精度的目的。由于分类样本的不完整影响了分类的精度,由于种植时间、灌溉水量以及气候变化等存在差异,导致不同区域内农作物的长势存在差异,样本点无法完全覆盖作物的各个生长阶段,这在一定程度上影响了分类的精度。因此在进一步研究中,需要充分考虑各个时间段的农作物生长状态,尽可能的保证样本点的充分覆盖。
参考文献
Daily evapotranspiration estimation at the field scale:Using the modified SEBS model and HJ-1 data in a desert-oasis area,northwestern China
[J].DOI:10.3390/w10050640 URL [本文引用: 1]
基于多时相遥感影像的作物种植信息提取
[J].
Crop planting information extraction based on multi-temporal remote sensing images
[J].
利用多时相TM影像进行作物分类方法
[J].
Crop classification method using multi-temporal TM images
[J].
基于NDVI时间序列影像的天山北坡经济带农业种植结构提取
[J].
Extraction of agricultural planting structure based on NDVI time series images in the economic zone of the northern slope of the Tianshan Mountains
[J].
基于多时相HJ卫星的冬小麦面积提取
[J].
Winter wheat area extraction based on multi-temporal HJ satellites
[J].
基于典型物候特征的MODIS-EVI时间序列数据农作物种植面积提取方法——小区域冬小麦实验研究
[J].
Extraction method of crop planting area from MODIS-EVI time series data based on typical phenological characteristics:Experimental study on small area winter wheat
[J].
基于SPOT-VGT NDVI时间序列的农牧交错带植被物候监测
[J].
Vegetation phenology monitoring in the agro-pastoral ecotone based on SPOT-VGT NDVI time series
[J].
Texture analysis and classification of ERS SAR images for map updating of urban areas in The Netherlands
[J].
利用GLCM纹理分析的高分辨率SAR图像建筑区检测
[J].
High-resolution SAR image building area detection using GLCM texture analysis
[J].
基于GF-1 WFV数据的玉米与大豆种植面积提取方法
[J].
Extraction method of corn and soybean planting area based on GF-1 WFV data
[J].
基于GF-1卫星遥感数据识别京津冀冬小麦面积
[J].
Recognizing the area of winter wheat in Beijing-Tianjin-Hebei based on GF-1 remote sensing data
[J].DOI:10.3724/SP.J.1006.2018.00762 URL [本文引用: 1]
基于GF-1卫星数据的农作物种植面积遥感抽样调查方法
[J].
Remote sensing sampling survey method of crop planting area based on GF-1 satellite data
[J].
基于GF-1与Landsat8 OLI影像的作物种植结构与产量分析
[J].
Analysis crops planting structure and yield based on GF-1 and Landsat8 OLI images
[J].
基于高分1号杭州湾南岸滨海陆地土地覆盖信息提取方法研究
[J].
Research on the extraction method of land cover information based on the coastal land on the south coast of Hangzhou Bay of GF-1
[J].
基于GF-1影像的东江流域面向对象土地利用分类
[J].
Object-oriented land use classification of Dongjiang basin based on GF-1 image
[J].
基于光谱和纹理特征综合的农作物种植结构提取方法研究
[D].
Research on the extraction method of crop planting structure based on the integration of spectrum and texture features
[D].
冬小麦种植面积遥感提取方法研究
[J].DOI:10.6046/gtzyyg.2013.04.02 [本文引用: 1]
Research on remote sensing extraction method of winter wheat planting area
[J].DOI:10.6046/gtzyyg.2013.04.02 [本文引用: 1]
基于GF-1影像NDVI年度间相关分析的冬小麦面积变化监测
[J].
Monitoring of winter wheat area change based on inter-annual correlation analysis of GF-1 image NDVI
[J].
基于GF-1/WFVNDVI时间序列数据的作物分类
[J].
Crop classification based on GF-1/WFVNDVI time series data
[J].
Multiple Resolution Texture Analysis and Classification
[J].
基于分形和灰度共生矩阵纹理特征的种植型药用植物遥感分类
[J].
Remote sensing classification of planted medicinal plants based on fractal and gray-level symbiotic matrix texture features
[J].
基于小波纹理和随机森林的猕猴桃果园遥感提取
[J].
Remote sensing extraction of kiwifruit orchard based on wavelet texture and random forest
[J].
基于GF遥感数据纹理分析识别制种玉米
[J].
Recognition of seed production corn based on texture analysis of GF remote sensing data
[J].
一种农作物种植结构月尺度动态提取方法:中国, CN110321861A
[P].
A dynamic extraction method of crop planting structure on a monthly scale:China, CN110321861A
[P].
纹理特征优选的农作物种植结构月尺度动态提取方法:中国, CN110909652A
[P].
The monthly-scale dynamic extraction method of crop planting structure for texture feature optimization:China,CN110909652A
[P].
农作物遥感分类特征变量选择研究现状与展望
[J].
Current status and prospects of research on selection of feature variables for crop remote sensing classification
[J].
单变量特征选择的苏北地区主要农作物遥感识别
[J].
Remote sensing identification of main crops in northern Jiangsu based on univariate feature selection
[J].
基于多特征的高光谱遥感土地利用信息提取
[J].
Multi-feature-based hyperspectral remote sensing land use information extraction
[J].
面向对象的特色农作物种植遥感调查方法研究
[J].
DOI:10.12082/dqxxkx.2018.180116
[本文引用: 1]
宁夏自治区具有土地、光能、引黄灌溉等优势,为宁夏特色农作物(硒砂瓜、枸杞、大枣)的生长提供了先天条件。快速准确地获取特色农作物的种植信息不仅是宁夏特色农作物监测、估产和灾害评估的重要依据,同时也是分析特色农作物结构分布变化和评价区域特色农业生产影响的重要凭证。近年来,随着航天技术和卫星传感器的不断发展,越来越多的学者将遥感技术运用到农作物种植信息的提取研究中。但是传统的遥感调查模型都是基于中低分辨遥感数据建立的,对于新的高分数据没有完备的信息提取模型。此外,基于GF-1遥感影像对类似宁夏特色农作物(硒砂瓜、枸杞、大枣)的信息提取研究相对较少,决策条件和分类模型的选择也难以满足高分农业的需求。基于此,本文利用国产GF-1 PMS遥感影像,在分析3类特色农作物光谱特征和纹理特征的基础上,建立了面向对象的支持向量机(SVM)分类模型,总体分类精度达到94.94%,Kappa系数为0.9174。同时将分类结果与传统的SVM分类结果相比较,研究发现面向对象的SVM模型的精度更高,效果最好,纹理信息的引入使光谱特征差异较小的枸杞和大枣更容易区分,有效降低了模型错分和漏分误差,改善了模型分类结果。研究结果为实现宁夏特色农作物的快速自动化提取提供了有效途径,也为开展农作物承保和受灾定损评估体系建设提供技术支撑。
Research on object-oriented remote sensing survey method of characteristic crop cultivation
[J].
/
〈 |
|
〉 |
