基于GF-6卫星的棉花种植面积提取研究——以山东省夏津县为例.pdf

基于G F 6卫星的棉花种植面积提取研究 以山东省夏津县为例收稿日期:2 0 2 2 1 1 1 0;修订日期:2 0 2 2 1 2 2 8;编辑:王敏基金项目:山东省农业农村厅2 0 2 1年度夏津县棉花绿色高质高效创建项目,鲁农计财字2 0 2 12 6号作者简介:马春莹(1 9 7 1),女,山东夏津人,高级工程师,主要从事自然资源方面的工作;Em a i l:2 9 8 1 3 1 0 9 7 3q q.c o m 马春莹(夏津县自然资源局,山东 德州 2 5 3 2 0 0)摘要:国产高分卫星6号(简称G F 6)具有高分辨率、宽覆盖、高质量和高效成像等特点,可为农业资源遥感监测提供可靠的数据支撑。本研究以G F 6卫星影像为主要数据源,综合利用遥感技术理论和方法开展对夏津县棉花种植面积的提取研究,利用同时相的哨兵2号(简称S e n t i n e l 2)影像解译结果进行对比分析。研究结果表明:基于G F 6卫星的夏津县棉花种植面积提取方案高效、准确,提取精度明显优于S e n t i n e l 2卫星,进一步证明了新增了红边波段的G F 6卫星能够显著增强作物的识别能力,其在大规模作物种植信息提取方面具有广阔的应用潜力。关键词:G F 6卫星;影像解译;棉花种植;山东省夏津县中图分类号:S 1 2 7 文献标识码:A d o i:1 0.1 2 1 2 8/j.i s s n.1 6 7 2 6 9 7 9.2 0 2 3.0 8.0 0 7引文格式:马春莹.基于G F 6卫星的棉花种植面积提取研究 以山东省夏津县为例J.山东国土资源,2 0 2 3,3 9(8):4 4 4 9.MAC h u n y i n g.S t u d yo nE x t r a c t i o no fC o t t o nP l a n t i n gA r e ai nX i a j i nC o u n t yB a s e do nG F 6S a t e l-l i t e S e t t i n gX i a j i nC o u n t y i nS h a n d o n gP r o v i n c ea sa nE x a m p l eJ.S h a n d o n gL a n da n dR e s o u r c e s,2 0 2 3,3 9(8):4 4 4 9.0 引言夏津是我国棉花种值大县,棉花种植面积和产量均稳居山东省前三位。近年来,随着乡村振兴战略的深入推进,夏津县高度重视棉花产业的发展,政府对棉花种植面积和空间分布的重视程度不断增强。如何准确、快速地获取大区域尺度下的棉花种植面积是掌握农作物基础信息的必要途径,同时对于落实棉花种植直补和各项惠农补贴、及时了解市场行情发挥着重要的导向作用1。早前,棉花种植面积主要通过人工实地调查的方式来获取,这种方法具有劳动强度大、速度慢、成本高等不足,且统计信息受人为因素影响大,在大尺度空间分布调查时很难得到精准的种植信息2。随着遥感技术的快速发展,卫星影像在现代农业管理中得到了广泛运用,早期卫星的空间分辨率较低,混合像元的存在限制了其在农作物面积提取方面的应用3。国产高分卫星6号(以下简称G F6)是国内第一颗搭载了能有效辨别作物类型的高空间分辨率遥感卫星,其在卫星传感器中首次配置了能有效反映农作物特有光谱特性的“红边”波段,该波段是反应植物生长状况的敏感性波段,能够显著提升作物精准识别的能力,为高精度、定量化的农业生产过程和农业资源环境要素监测提供了可靠的数据支撑。本文选取山东省夏津县为研究区,以G F6卫星影像为主要数据源,哨兵2号(以下简称S e n-t i n e l 2)卫星影像作为参照,采用面向对象的决策树型分类方法,并选用易康遥感影像分析软件(以下简称e C o g n i t o n)进行分类,实现了该县棉花种植面积的精准提取。1 数据源与方法1.1 研究区概况夏津县地处鲁西北平原、鲁冀两省交界处,北依44第3 9卷第8期 山 东 国 土 资 源 2 0 2 3年8月德州,南靠聊城,西临京杭大运河。夏津县属暖温带半湿润大陆性季风气候,冷热、干湿明显,四季变化明显,晚秋经常出现干旱,干旱对于喜温好光的棉花影响比其他作物小,这是夏津县种植棉花历史悠久的气候原因。近年来,夏津县强化要素集聚,科学划定棉花保护区,统筹高标准农田、小农水、土地整理等项目,重点提升棉田基建,助力绿色高质高效棉花生产。1.2 数据源本研究选取G F 6卫星影像作为遥感数据源,G F 6卫星配置2m全色/8m多光谱高分辨率相机、1 6m多光谱中分辨率宽幅相机,2m全色/8m多光谱相机观测幅宽9 0k m,1 6m多光谱相机观测幅宽8 0 0k m(表1)。为了检验G F 6卫星棉花种植面积的提取精度,同时选取同时相的S e n t i n e l 2影像数据进行提取,对比两者的解译精度。采用的其他数据如耕地分布图、向量行政区域规划图、农作物种植面积调查统计年鉴等资料均来自官方统计数据。利用E NV I 5.3软件对卫星影像数据进行预处理,主要包括正射校正、几何配准、影像融合、图像掩膜、辐射校正等步骤45。表1 G F 6卫星参数名称参数名称参数空间分辨率全色:优于2 m波段8多光谱:优于8 m量化值1 2 b i t光谱段范围P:0.4 50.9 0m覆盖宽度9 0k m蓝:0.4 50.5 2m信噪比全色:低端2 8 d B,高端4 7 d B绿:0.5 20.6 0m多光谱:低端2 0 d B,高端4 6 d B红:0.6 30.6 9m辐射定标精度绝对定标精度优于7%B 4:0.7 60.9 0m相对定标精度优于3%1.3 研究方法1.3.1 提取方案合理的作物识别方法可以提高作物识别的精度,有利于获取到准确的棉花空间分布信息与面积6。参照棉花的光谱特征和遥感影像颜色特征,采用构建决策树模型方式进行分析提取,得到基于监督分类方法下棉花的像元结果,按照棉花掩膜规则生成棉花掩膜后,将其应用在夏津县卫星影像上,由于其他作物基本被掩膜遮蔽,棉花光谱和掩膜区域内其他地物的光谱可分性较好,再使用e C o g n i-t o n软件进行分类,排除符合掩膜规则的农田以外地物,达到精确提取棉花种植面积的目的7,具体流程见图1。1.3.2 精度评价的方法将估测面积和统计数据进行对比,分析提取精度,再采用K a p p a系数对提取方法进行精度评价,K a p p a系数的计算基于混淆矩阵,K值的范围通常为11之间,一般大于08,K a p p a系数计算公式见式(1):K=P o P e1P e(1)其中:P o所有对角线元素之和除以整个矩阵元素之和得到的商值;P e每个真实样本乘以预测样本的和再除以样本个数的平方的值。图1 提取流程图2 棉花种植面积提取2.1 棉花光谱特征分析夏津县的棉花种植规律为4月份开始播种,此时的植被覆盖率较低,56月份,大多数的棉花还处于幼苗期,卫星影像呈现为淡绿色,此时地面植物的反射光谱特征基本相似,难以进行分类提取。8月初,棉花生长进入第一个花铃期,此阶段的棉花繁殖和发育最强烈,23周后棉花开始进入初开和盛花期,此时的植株长势达到最高峰,卫星图像呈现出鲜绿色。通过对棉花种植结构特征、遥感影像颜色特征、光谱特征的分析,可知棉花的最佳遥感识别期54第3 9卷第8期 技 术 方 法 2 0 2 3年8月为8月中下旬,此时地物间植被指数差异较大,棉花较容易的识别出来(图2)。图2 夏津县8月份几种地物的反射波谱曲线图2.2 植被指数分析参照夏津县棉花调查整理数据,利用其中有准确经纬度、种植作物信息的地面调查点,建立研究区作物样本图斑,分别建立棉花、玉米、水体、林地的训练样本和验证样本,统计遥感数据训练样本可分离系数的变化情况。本研究选取卫星影像中R、G、B和N I R共4个1 0m分辨率的波段图像进行融合,利用E NV 1 5.6和A r c G I S 1 0.7对融合后影像重新构建金字塔,然后提取相对应的每个采样点的光谱波段反射率值,结合卫星遥感指数已有算法,并利用E NV I 5.3平台的B a n dM a t h工具计算出不同土地利用类型的植被指数9。2.3 提取步骤本研 究 利 用e C o g n i t o n软 件 进 行 分 类 提 取,e C o g n i t o n是一款高分辨率遥感影像分析软件,主要通过面向对象的影像分割技术,来实现影像的分析与数据获取1 0,提取步骤如下:(1)棉花粗提。分类的基本过程中,先选定目标地物样本对象,由软件系统自动计算该类地物的属性值,然后建立特征空间结构1 1。以已确定样本的特点空间结构为核心,通过统计方法分析其余未分类型的摄影对象和该区位特性空隙之间的差距,那么影像对象离哪个区位的特性空隙最近,即会被分类在该类别中。8月棉花的平均反射率均高于树类和其他作物,而水域和灌木因为其独特的地理特点导致波谱特征几乎不会有太大明显变化,其他地类的波谱曲线的形态都有较大变动1 2,由此,参照不同作物的像元可以得到棉花粗提结果。(2)棉花掩膜。按照掩膜规则生成棉花掩膜后,将其应用在卫星影像上,由于其他作物基本被掩膜遮蔽,棉花光谱和掩膜区域内其他地物的光谱可分性较好。通过分割试验得到研究区域内的影像分割参数,调整各参数权值,包括分割尺寸、图形紧致度、平滑度和色彩参数权重等,得出在多种分割图层下解析不同地物最适宜的切割尺寸1 31 4。根据分割结果来调整样本点的疏密,使得遥感影像精度和分割尺度、样本集相匹配,最终使用的分割尺度参数为1 0 0、形状参数为0.3、精致度参数为0.7。使用支持向量机分类方法进行监督分类,排除符合掩膜规则的农田以外地物,达到精确识别棉花种植面积目的。(3)棉花提取。已经分离出的农作物图层中含有农田、苗木和少量的不透水地表(池塘等),首先在分割尺寸为1 0 0的影像对象层中进行粗提纯,再使用光谱特性“R a t i o r e d”和几何特性“L e n g t h/W i d t h”进行初步提纯;随后再将初步提纯后的作物类别传送到分割尺寸为1 0 0的影像对象层中再次进行提纯。提取过程如图3所示。3 结果与分析综合利用遥感数据源,采用面向对象的决策树型分类方法,并基于e C o g n i t o n的棉花种植信息分类提取方案,得出夏津县棉花识别的最佳解译结果(图4、图5),结合实地验证点、作物分类样本对解译精度及棉花提取精度进行验证1 5。3.1 遥感提取面积和官方统计数据对比根据分类结果,利用E NV I中的S t a t i s t i c(统计)工具,计算得出棉花所占的像元总数,然后计算其种植面积(种植面积=像元数像元面积)1 6,最后得到基于G F 6卫星和S e n t i n e l 2卫星提取的夏津县2 0 2 1年棉花种植面积分别为1 13 6 9.9 6h m2和1 11 0 9.9 6h m2。将估测结果与官方统计数据进行对比(统计结果见表2)。表2 提取面积与官方统计面积对比 单位:h m2面积G F 6卫星S e n t i n e l 2卫星官方统计面积1 1 8 5 6.6 61 1 8 5 6.6 6遥感提取面积1 1 3 6 9.9 61 1 1 0 9.9 6差值4 8 6.77 4 6.7从表2数据看出,基于G F 6卫星遥感提取的面积更加接近官方统计的面积。3.2 解译精度检验分别计算基于G F 6卫星和S e n t i n e l 2卫星64第3 9卷第8期 山 东 国 土 资 源 2 0 2 3年8月图3 棉花提取过程图4 基于G F 6卫星提取的夏津县棉花分布图的棉花种植信息提取的总体分类精度、K a p p a系数、生产精度和用户精度,统计数据见表3。表3 作物验证样本检验解译精度汇总表分类方法G F 6卫星S e n t i n e l 2卫星总体分类精度9 3.8 8%9 1.8 4%K a p p a系数0.9 2 2 30.8 9 0 3生产精度9 1.3 1%8 9.1 7%用户精度9 0.4 5%8 8.9 7%由表3可知,G F 6卫星在总体分类精度、生产精度、用户精度等方面的数值均高出S e n t i n e l 2卫星约2个百分点,证明了其对棉花形状、纹理、大小等特征均有更高的辨别度,红边波段的确能提高作物识别的能力1 6,提高作物的解译精度。3.3 实地对比分析为了评估分类的准确性,在夏津县棉花分布图中选取4 9个验证点与实地进行对比分析,由于验证点是随机选取,且未均匀分布于研究区域内,所以本次统计仅计算总体分类精度1 7。经实地调查统计,利用G F 6卫星提取的4 9个验证点中有4 6个分类正确,3个分类错误,总体精度达到9 3.8 8%;利用S e n t i n e l 2卫星提取的4 9个验证点中有4 5个分类正确,4个分类错误,总体精度达到9 1.8 4%。实地调查分类错误的原因,其中2个验证点周围种植有大豆和芦苇,提取结果将其误判为棉花;2个验证点套种有其他作物,提取结果产生漏分。G F 6卫星分类结果图斑具有更好的完整性,误判和漏分的情74第3 9卷第8期 技 术 方 法 2 0 2 3年8月图5 基于S e n t i n e l 2卫星提取的夏津县棉花分布图况要优于S e n t i n e l 2卫星。同时,选取某处验证样方对棉花种植面积提取的精度进行分析,检测方法为利用G P S R T K实地测量验证样方内的棉花种植面积,与G F 6卫星和S e n t i n e l 2卫星提取的验证样方内的棉花种植面积进行对比,统计结果见表4。从表4统计数据可以看出,G F 6卫星和S e n t i n e l 2卫星提取的棉花种植面积相对误差均小于1 0%,说明2种卫星提取的面积精度都比较高,G F 6卫星的全色分辨率达到2m,S e n t i n e l 2卫星的分辨率为1 0m,G F 6卫星对作物种类及边界匹配更加准确,其对验证样方中棉花地块的形状、纹理、大小等特征上的表现更加突出。表4 棉花种植面积提取相对误差统计表数据源提取结果/h m2实测结果/h m2相对误差/%G F 6卫星3 8.4 63 6.5 85.1S e n t i n e l 2卫星4 0.1 23 6.5 89.74 结语本研究以G F 6卫星影像为主要数据源,采用面向对象的决策树型分类方法,并基于e C o g n i t o n遥感影像分析软件成功提取出夏津县棉花种植面积,该方法在处理高维度、大体积遥感数据中展现出天然的优势,其对大区域尺度的历史遥感影像处理分析能够更高效、准确的获取区域种植空间分布信息,通过研究总结如下:(1)经检验,基于G F 6卫星的棉花面积提取结果在精度方面明显优于S e n t i n e l 2卫星,证明了G F 6卫星新增的红边波段能够显著增强作物的识别能力,红边波段特征与作物生长状况、植被特征的相关性更强,其8 5 0k m的宽幅数据可较好地用于中小区域尺度的地物提取,使得复杂地块作物面积监测结果更加客观合理,G F 6卫星能为遥感解译提供更加丰富的信息,其在作物面积提取方面具有广泛的应用潜力,极大地推动了卫星遥感技术的国产化应用。(2)随着对地观测技术的不断发展,高空间、时间以及光谱分辨率的卫星持续发射,其载荷类型不断丰富,空间信息网络通信能力不断提升,在轨数据处理能力不断加强,为农业资源调查等应用提供了更加丰富的数据源。考虑到单一遥感系统观测往往并不全面,后续通过多传感器进行数据融合研究能更有效发挥多平台互补观测的优势,即将包含同一目标或场景,在空间、时间、光谱上冗余或互补的多源遥感数据按照一定规则(或算法)进行运算处理,84第3 9卷第8期 山 东 国 土 资 源 2 0 2 3年8月获得比任何单一数据更精确、完整、有效的信息,实现更加精准、全面的作物分类提取。(3)为提高经济效益,当前的棉花田存在很多套种方式,呈现出“花花田”的现象,即在一块耕地上按照一定的行、株距和占地的宽窄比例种植几种作物,对棉花种植面积识别与提取造成较大的影响,在后续研究中,应当加大对多种作物光谱反射特征的研究,寻找作物之间光谱反射特征的差异,以进一步提高作物识别与提取的能力,为区域发展“精准农业”“智慧农业”建设提供重要的技术支撑。参考文献:1 潘力,夏浩铭,王瑞萌,等.基于G o o g l eE a r t hE n g i n e的淮河流域越冬作物种植面积制图J.农业工程学报,2 0 2 1,3 7(1 8):2 1 1 2 1 8.2 白雪武,武红旗,吕昱,等.基于L a n d s a t 8和高分一号影像的沙湾县作物种类识别研究J.山东农业科学,2 0 2 2(2):1 5 6 1 6 2.3 张悦琦,李荣平,穆西晗,等.基于多时相G F 6遥感影像的水稻种植面积提取J.农业工程学报,2 0 2 1(1 7):1 8 9 1 9 6.4 李亚妮,曹建君,杨树文,等.基于决策树的大尺度复杂地区夏收作物遥感提取与分析J.江苏农业学报,2 0 2 2(5):1 2 5 71 2 6 4.5 张华平,王增辉,谢颂诗,等.基于国产高分遥感数据的矿山环境监测系统研究J.山东国土资源,2 0 2 1,3 7(1 2):1 0 0 1 0 4.6 林娜,陈宏,李志鹏,等.基于G F 1号影像的南方棉花种植信息提取J.地理空间信息,2 0 2 1(3):6 0 6 3.7 赵文驰,宋伟东,陈敏.国产高分辨率遥感卫星融合方法比较J.测绘与空间地理信息,2 0 1 9(1 1):1 5 4 1 5 8.8 董文全,蒙继华.遥感数据时空融合研究进展及展望J.国土资源遥感,2 0 1 8(2):1 1 1.9 许青云,杨贵军,龙慧灵,等.基于MO D I SN D V I多年时序数据的农作物种植识别J.农业工程学报,2 0 1 4(1 1):1 3 4 1 4 4.1 0 张金盈,崔靓,徐凤玲,等.海量国产高分辨率卫星影像优化处理研究J.山东国土资源,2 0 2 0,3 6(9):6 5 6 9.1 1 李方杰,任建强,吴尚蓉,等.N D V I时序相似性对冬小麦种植面积总量控制的制图精度影响J.农业工程学报,2 0 2 1(9):1 2 7 1 3 9.1 2 玉苏甫买买提,吐尔逊艾山,买合皮热提吾拉木.新疆渭 库绿洲棉花种植面积遥感监测研究J.农业现代化研究,2 0 1 4(2):2 4 0 2 4 3.1 3 王霄煜,雷钧.基于高分一号卫星数据的新疆棉花种植面积信息提取研究J.新疆农业科技,2 0 2 1(1):2 3 2 6.1 4 刘焕军,闫岩,张新乐,等.面向农业区划的作物种植结构遥感提取J.中国农业资源与区划,2 0 1 7(8):4 3 5 4.1 5 李国庭,王德强,赵德良,等.基于多源异构数据的粮食作物种植面积提取J.山东国土资源,2 0 2 2,3 8(8):5 1 5 5.1 6 王克晓,周蕊,王茜,等.基于G F 1的山地主要农作物种植信息提取研究J.测绘与空间地理信息,2 0 2 0(6):3 3 3 6.1 7 何真,胡洁,蔡志文,等.协同多时相国产G F 1和G F 6卫星影像的艾草遥感识别J.农业工程学报,2 0 2 2,3 8(1):1 8 61 9 6.S t u d yo nE x t r a c t i o no fC o t t o nP l a n t i n gA r e a i nX i a j i nC o u n t yB a s e do nG F 6S a t e l l i t e S e t t i n gX i a j i nC o u n t y i nS h a n d o n gP r o v i n c ea sa nE x a m p l eMAC h u n y i n g(X i a j i nB u r e a uo fN a t u r a lR e s o u r c e s,S h a n d o n gD e z h o u2 5 3 2 0 0,C h i n a)A b s t r a c t:T h ed o m e s t i ch i g hr e s o l u t i o ns a t e l l i t e6(G F 6f o rs h o r t)h a s t h ec h a r a c t e r i s t i c so fh i g hr e s o l u-t i o n,w i d ec o v e r a g e,h i g hq u a l i t ya n de f f i c i e n ti m a g i n g.I tc a np r o v i d er e l i a b l ed a t as u p p o r tf o rr e m o t es e n s i n gm o n i t o r i n go f a g r i c u l t u r a l r e s o u r c e s.I nt h i ss t u d y,t a k i n gG F 6s a t e l l i t e i m a g ea st h em a i nd a t as o u r c e,b yu s i n gt h e t h e o r ya n dm e t h o do f r e m o t es e n s i n gt e c h n o l o g yc o m p r e h e n s i v e l y,t h ec o t t o np l a n t-i n ga r e a sh a v eb e e ne x t r a c t e d i nX i a j i nc o u n t y.A t t h es a m e t i m e,t h e i n t e r p r e t a t i o nr e s u l t so fS e n t i n e l 2(S e n t i n e l 2f o r s h o r t)i m a g ea r eu s e d f o r c o m p a r a t i v ea n a l y s i s.I t i s s h o w e d t h a t t h ee x t r a c t i o ns c h e m eo fc o t t o np l a n t i n ga r e a i nX i a j i nc o u n t yb a s e do nG F 6s a t e l l i t e i se f f i c i e n ta n da c c u r a t e,a n dt h ee x t r a c t i o na c c u r a c y i ss i g n i f i c a n t l yb e t t e r t h a nt h a to fS e n t i n e l 2s a t e l l i t e.I t i sp r o v e dt h a t t h en e w l ya d d e dG F 6s a t e l l i t ew i t hr e de d g eb a n dc a ns i g n i f i c a n t l ye n h a n c e t h e r e c o g n i t i o na b i l i t yo f c r o p s.I th a sb r o a da p p l i c a-t i o np o t e n t i a l i nl a r g e s c a l ec r o pp l a n t i n g i n f o r m a t i o ne x t r a c t i o n.K e yw o r d s:G F-6s a t e l l i t e;i m a g e i n t e r p r e t a t i o n;c o t t o np l a n t i n g;X i a j i nc o u n t y i nS h a n d o n gp r o v i n c e94第3 9卷第8期 技 术 方 法 2 0 2 3年8月