自然资源遥感, 2022, 34(3): 184-195 doi: 10.6046/zrzyyg.2021320

技术应用

基于多时相Sentinel-1水稻种植范围提取

查东平,1,2, 蔡海生,1, 张学玲1, 何庆港1

1.江西农业大学/江西省鄱阳湖流域农业资源与生态重点实验室,南昌 330045

2.江西旅游商贸职业学院,南昌 330100

Extraction of paddy fields using multi-temporal Sentinel-1 images

ZHA Dongping,1,2, CAI Haisheng,1, ZHANG Xueling1, HE Qinggang1

1. Key Laboratory of Poyang Lake Watershed Agricultural Resources and Ecology of Jiangxi Province, Jiangxi Agricultural University, Nanchang 330045, China

2. Jiangxi Tourism and Commerce Vocational College, Nanchang 330100, China

通讯作者: 蔡海生(1972-)男,教授,博士生导师,主要从事土地生态、土地经济研究。Email:chs@jxau.edu.cn。

责任编辑: 李瑜

收稿日期: 2021-09-30   修回日期: 2022-04-30  

基金资助: 国家自然科学基金项目“鄱阳湖流域土地集约利用与生态安全耦合发展及其综合响应机制研究”(31660140)
国家自然科学基金项目“亚热带山地草甸退化格局与土壤生态过程研究——以江西武功山为例”(31560150)
江西省“十三五”社科规划项“新时期共享型农业经营体系的政策和机制创新研究”(17YJ11)
江西省高校人文社科规划项目“新时期农村宅基地退出及其补偿机制研究——以江西省试点情况调研为例”(GL17113)
及江西省高校人文社科重点研究基地项目“基于乡村振兴战略下农村土地资源配置制度创新研究”(2018-32)

Received: 2021-09-30   Revised: 2022-04-30  

作者简介 About authors

查东平(1985-)男,博士研究生,主要从事农业资源与环境相关研究。Email: 345914421@qq.com。

摘要

通过遥感手段监测与提取水稻种植范围的是农业现代化管理的重要手段,南方地区春夏季节多云雨天气,难以获得有效的光学数据。为了准确提取多云雨地区的水稻种植范围,以江西省南昌县蒋巷镇为研究区,以Sentinel-1 多时相数据为数据源开展研究。通过早稻生长关键物候期Sentinel-1 SAR数据各地类后向散射系数,计算不同物候期条件下水稻田与其他地类的J-M距离,分析不同物候期SAR数据组合情况下水稻田与其他地类的J-M距离变化,获得提取水稻种植范围的最佳物候期影像。分别采用随机森林法、最大似然法、支持向量机和神经网络法进行分类,并对比验证分类精度。结果表明,孕穗期(6月14日)、三叶期(4月21日)、移栽期(5月3日)、二季晚稻移栽期(7月26日)组合SAR数据为水稻田提取最优时相组合。采用随机森林方法进行分类能够获得较高精度,研究区地物分类总体精度达到0.943,Kappa系数0.932。研究对于多云雨地区采用SAR数据开展水稻田制图,在时相选择和分类方法上有一定的借鉴意义和参考价值。

关键词： 遥感监测; 信息提取; SAR; Sentinel-1

Abstract

The monitoring and information extraction of paddy fields using remote sensing techniques is an important means for modern agricultural management. However, it is difficult to obtain effective optical monitoring data of south China due to the frequent cloudy and rainy weather in spring and summer in this area. To accurately extract information on paddy fields in areas subject to frequent cloudy and rainy weather, this study investigated the paddy fields in Jiangxiang Town in Nanchang County, Jiangxi Province, using multi-temporal Sentinel-1 SAR data as the data source. Specifically, this study calculated the J-M distance between paddy fields and other land types in different phenological periods, analyzed the changes in the distance based on the backscattering coefficients of various land types in key phenological periods, and then obtained the best phenological images for the information extraction of paddy fields. Afterward, this study conducted ground object classification using methods such as random forest, maximum likelihood, support vector machine, and neural network and then compared and verified the classification accuracy. The results are as follows. The combined SAR data of the different stages including booting stage (June 14), trefeil stage (April 21), transplantion period (May 3), and transplanting peried of second season late rice (July 26) is the optimal temporal combination for the information extraction of paddy fields. Higher classification accuracy of ground objects in the study area can be obtained using the random forest method, with overall classification accuracy of up to 0.943 and a Kappa coefficient of 0.932. This study conducted the mapping of paddy fields in areas with frequent cloudy and rainy weather using SAR data and will provide important references for the temporal selection and classification.

Keywords： remote sensing and monitoring; information extraction; SAR; Sentinel-1

PDF (7594KB) 元数据 多维度评价 相关文章 导出 EndNote| Ris| Bibtex  收藏本文

本文引用格式

查东平, 蔡海生, 张学玲, 何庆港. 基于多时相Sentinel-1水稻种植范围提取[J]. 自然资源遥感, 2022, 34(3): 184-195 doi:10.6046/zrzyyg.2021320

ZHA Dongping, CAI Haisheng, ZHANG Xueling, HE Qinggang. Extraction of paddy fields using multi-temporal Sentinel-1 images[J]. Remote Sensing for Land & Resources, 2022, 34(3): 184-195 doi:10.6046/zrzyyg.2021320

0 引言

及时准确地监测和获取水稻种植的空间分布范围对保障全国粮食安全具有重要意义[1-2]。传统的统计方法有抽样调查和逐级上报,这两种方法获取数据成本巨大,且存在很严重的信息滞后[3]。遥感技术为大尺度农作物种植范围监测和提取提供了有效的技术支撑[4-5]。目前多采用Landsat系列、MODIS或Sentinel-2等光学卫星数据对农作物信息进行提取,这些数据免费且易于获取,相关学者基于此类数据开展了大量研究并取得了较好的成效[6]。Zhu等[7]利用Landsat系列数据对我国南方水稻产区的水稻种植面积进行了监测,结果认为2000—2019年水稻总种植面积减少了171 474 km2; Cao等[8]采用Landsat数据识别早稻、晚稻、双季稻和中稻,总体精度在71.83%~85.81%之间; Li等[9]通过融合MODIS和Landsat数据,采用随机森林模型和决策树对农作物进行提取,精度分别达到90%和89.7%。利用光学卫星影像数据对水稻田的提取方法多要依靠多时相数据,单时相数据提取的精度往往比较低,但是南方地区水稻生长季节往往为多云雨天气,光学数据质量难以保障,影响了实际应用成效[10-11]。

合成孔径雷达(synthetic aperture Radar,SAR)不受天气影响,可以全天时全天候对地观测,有效解决了光学数据获取难的问题[12⇓-14]。基于SAR数据的水稻田分类虽然还没有大范围业务化应用,但是也开展了大量的研究,如Son等[15]利用Sentinel-1 SAR数据进行了水稻田提取,第一季和第二季水稻种植面积提取精度分别为84.2%和82.6%,用户精度分别为82.1%和85.3%; Graman等[16]利用Sentinel-1数据对Rabi地区水稻田进行提取,总体精度达到91.5%,Kappa系数为0.83。前人的研究证明了Sentinel-1数据VV(垂直)和VH(交叉)极化数据对提取水稻田的可行性。但是研究没有从水稻不同物候期数据去进行对比,筛选水稻田识别的最佳时相数据。目前成熟的,研究较多的分类方法包括最大似然法、支持向量机等,基于机器学习算法的分类方法应用在SAR图像水稻田提取的研究不多[17-18]。本研究拟通过分析不同物候期SAR数据水稻田与其他地类之间的可分离性,筛选出水稻田识别最佳时相组合,采用随机森林方法进行分类,以期得到较高精度的分类结果,为今后在水稻田识别和水稻生长监测中,因为云雨天气难以获得连续高质量的光学数据的情况下,采用SAR进行水稻监测提供方法上的参考。

1 研究区概况与数据源

1.1 研究区概况

研究区(图1)位于江西省南昌市南昌县北部、赣江下游蒋巷镇。地处E115°65',N28°32'。蒋巷镇四面环水,东临鄱阳湖、南北傍赣江、西接南昌市城郊,赣江南、中支流环抱而过。全镇国土总面积约为240 km2,南北宽约8 km,东西长约45 km。受亚热带季风气候影响,研究区内气候温暖湿润,光照充足,主要种植作物为水稻,是南昌县重要的水稻种植区。

图1

新窗口打开| 下载原图ZIP| 生成PPT

图1   研究区位置

Fig.1   The location of research area


1.2 SAR数据获取与处理

选取涵盖2020年早稻整个物候期的8期Sentinel-1地距影像产品数据(https://search.asf.alaska.edu/)。SAR过境时间分别对应水稻生长的出苗期(4月9日)、三叶期(4月21日)、移栽期(5月3日,直播稻为幼苗早期)、返青期(5月15日,直播稻为幼苗中期)、分蘖期(6月2日)、孕穗期(6月14日)、乳熟期(7月8日)、收割和二季晚稻移栽期(7月26日)。在ENVI SARscape 5.5软件中对SAR数据进行处理(图2),包括读取、多视处理、配准和辐射定标等操作,采用“Adaptive Non-Local SAR Fil-tering”方法进行滤波处理,减少相干斑噪声。

图2

新窗口打开| 下载原图ZIP| 生成PPT

图2   多时相Sentinel-1数据处理方法

Fig.2   Multi-temporal Sentinel-1 data processing method


1.3 外业数据采集

为获取研究区的解译样本及对结果进行验证分析,对研究区开展了现场踏勘调查,获得了研究区主要地物类型分布,此外于2020年7月17日采用了大疆DJI Mavic 2无人机对研究区进行了飞行录像(无人机飞行轨迹如图3所示)。为节约成本,无人机采用“抬头”拍摄,飞行过程中调整拍摄角度,从受局限的俯视过程过渡到开阔的视角,“抬头”拍摄镜头前方可观测距离超过1 km,视野范围影像截图如图4所示。

图3

新窗口打开| 下载原图ZIP| 生成PPT

图3   无人机飞行轨迹图

Fig.3   Flight trajectory map of UAV


图4

新窗口打开| 下载原图ZIP| 生成PPT

图4   无人机录像截图

Fig.4   Video screenshot of UAV


本文研究目的是提取水稻种植范围,结合土地利用/覆被类型特点,参考邓刚等[19]、张悦琦等[20]对水稻种植范围提取辅助分类,依据多次野外调查结果,确定分类系统包括: 水稻田、林地、水体、莲田、人工表面和湿地等6种地类。根据外业调查和无人机飞行录像数据,在同时相GF-1影像图上勾绘地物标志样本。每类地物人工勾绘采选200个样本点,选择手工勾绘确保在影像图上获得纯净像元,便于后期分析。

在ArcGIS中对样本点矢量数据进行编号,奇数序号样本点导出并命名为训练样本,偶数样本点导出并命名为验证样本,各地类训练样本点和验证样本均100个,并且大致均匀分布于研究区。训练样本点和验证样本点分布情况如图5所示。

图5

新窗口打开| 下载原图ZIP| 生成PPT

图5   样本点分布图

Fig.5   Sample points distribution map


2 研究方法

2.1 后向散射的一般特征

SAR使用微波作为发射电磁波,电磁波照射到地物上,一部分能量通过反射进入雷达接收天线,另外一部分能量则透过地物表面继续向下[21]。雷达照射表面如果是粗糙地物发生散射,能量朝各个方向,如果是光滑表面则发生镜面反射[22]。雷达接收天线接收的能量取决于雷达波长、入射角和地面粗糙程度[21]。一般来说长波比短波具有更好的穿透性,对于探测冠层下垫面具有更好的效果。入射角越小后向散射回的能量越大,大角度的回波信息中包含了地物粗糙程度信息[21,23]。

通常,水面和光滑的地表可认为光滑表面,电磁波照射后易发生镜面反射,雷达天线接收不到或者只能接收少部分回波信号,在SAR图像上呈现黑色[24]。随着入射角的变化,接收到的信号也会有所不同,当入射角达到90°时,反射信号直接进入雷达天线,图像上呈现明亮点[24]。不同极化数据,水平极化和垂直极化反射和散射后到接收器上的能量不同,两种方式的差别可以达到5 dB,同向极化(HH和VV)与交叉极化(HV和VH)的后向散射系数也不同,二者可相差10~15 dB[25]。除水体,大部分地物都可认为是程度不相同的粗糙表面[16,25]。根据瑞利散射原则,地物散射特性跟地物表面和波长有关,地物的散射十分复杂。通常,粗糙表面的散射是一条平缓的曲线,曲线斜率变化与入射角相关性较小。农作物在不同的物候期,表现出来的散射特性不同。SAR接收天线接收到的地物后向散射能量受地物特性影响,如植被的密度、植被的含水量、植被高度、下垫面类型等。不同时相、作物处于不同的物候期,会产生不同的雷达回波,是雷达监测及农作物制图的基础。

2.2 随机森林分类法

本文随机森林训练流程如图6所示。随机森林是一个包含多个决策树的分类器[26]。1995年贝尔实验室华裔美国人何天琴(Tin Kam Ho)提出了随机决策森林(random decision forest),2001年Leo Breiman和Adele Cutler发展并推出随机森林算法[27-28]。该算法通过构建多个决策树作为基本的分类器,每个决策树都是一个分类器,本质属于机器学习[26]。随机森林方法通过随机方式从原始训练集中选择K个训练集,对训练集构建分类决策树,集成多个分类投票的结果,将投票次数最多的分类作为最终地类并输出[26]。

图6

新窗口打开| 下载原图ZIP| 生成PPT

图6   随机森林训练流程图

Fig.6   Random forest training process


本研究随机森林分类通过对开源软件ENMAP-BOX v2.1.1源码封装后在ENVI 5.5软件中运行。

2.3 地物散射特征分析

2.3.1 水稻田散射特征分析

4月9日—7月26日水稻田的后向散射系数均值见表1,VV极化和VH极化后向散射系数变化趋势见图7。

表1   水稻田不同时相平均后向散射系数

Tab.1  Average backscattering coefficient of paddy field in different time-phase(dB)

极化4月9日4月21日5月3日5月15日6月2日6月14日7月8日7月26日
VV-10.34-12.23-10.82-9.63-10.65-11.96-9.64-11.47
VH-18.82-19.92-20.44-19.63-17.98-19.36-18.03-17.39

新窗口打开| 下载CSV


图7

新窗口打开| 下载原图ZIP| 生成PPT

图7   水稻田不同时相平均后向散射系数

Fig.7   Average backscattering coefficient of paddy field in different time-phase


水稻属于水生植物,大部分时间水稻根部被水体淹没,雷达波照射到水稻后大概可以分为3种散射过程:

1)水稻下垫面水体的镜面反射,水稻生长的早期,镜面反射是主要的散射过程,因此早期阶段水稻田的后向散射值较低。

2)水稻叶面的奇次散射,电磁波照射到水稻叶片后一部分电磁波发生奇次散射,直接进入雷达接收天线。

3)水稻叶片和茎秆的体散射,主要表现在雷达波到达水面后,水面反射到水稻茎秆或者叶片发生的多次散射,雷达波到达叶片、叶片散射到其他茎秆叶片或者水面发生的散射。

研究区水稻为直播稻,但水稻种植时间略有不同,4月初部分地区水稻已播种,部分地区水田尚未耕种长满野草,野草叶面积指数要比水稻种植区高,其后向散射系数也更高。4月21日,水稻处于三叶期,但水稻叶面积指数仍小于野草,后向散射系数要低于4月9日。随着水稻的生长,四月下旬三叶期水稻冠层本身散射和冠层散射逐渐增强。从极化情况来看,VV极化的穿透性要高于VH,其后向散射能量包含更多的植被冠层和植被下垫面的回波能量。4月下旬至5月中旬,随着水稻的生长VV极化数据呈现上升趋势,但是VH极化数据变化情况不大。根据地面数据采集情况,拔节、孕穗期之后水稻含水量逐渐下降,水稻的介电常数与含水量存在正相关性,结果受水稻的介电常数和降雨等情况的综合影响。7月26日已种植第二季晚稻,种植方式为插播,正处于返青期,此时仍受水面镜面反射主导。

2.3.2 其他地物散射特征分析

对莲田、林地、人工表面、水体和湿地等样本点在SAR时相序列后向散射数据中进行取值计算,各地物VV极化后向散射系数均值结果如表2和图8所示,地物散射系数箱型图见图9; 各地物VH极化后向散射系数均值如表3和图10所示,箱型图见图11。从VV极化各地物后向散射系数均值来看,各地物后向散射系数在-21~-5 dB之间,其中水体最低,人工表面最高。从图9可知,仅从VV极化数据,水体、莲田两种地类均具有较好的分离度,最佳区分的时相分别为5月15日、6月2日和6月14日。莲田与其他地类进行区分,最佳区分时相为4月9日和5月3日。值得一提的是,随着时间推移,莲田的VV后向散射系数逐渐提升,在5月3日之前,莲田的后向散射系数要低于水稻田、湿地等地类,6月14日之后,莲田后向散射系数要高于这两种地类,在6月26日甚至要高于林地的后向散射系数。林地的后向散射系数一直比较稳定,而且标准差较小,各时相均能与水稻田和莲田区分,与人工表面的最佳区分时相为4月21日,与莲田在4月9日、4月21日、5月3日和7月26日数据中均能较好区分,与水稻田在4月21日、6月14日、7月8日和7月26日均能较好区分; 与湿地在7月8日和7月26日具有较好的分离度。水稻田在水稻生长后期与其他地类能够较好区分,如6月14日、7月8日和7月26日。7月8日和7月26日湿地与人工表面、林地、莲田和水稻田具有较好的分离度,但是这2期数据湿地VV极化后向散射系数与水体极为相似,但是其他几个时相湿地可以与水体进行分离。

表2   VV极化各地物不同时相后向散射系数

Tab.2  Backscattering coefficients of VV polarized objects in different time-phase(dB)

地类4月9日4月21日5月3日5月15日6月2日6月14日7月8日7月26日
莲田-16.02-15.14-12.46-11.70-9.27-7.84-11.77-6.95
林地-9.58-8.88-9.70-8.25-7.96-8.43-8.21-8.53
人工表面-7.46-5.20-7.61-5.37-5.37-6.75-5.47-6.68
水体-19.42-18.26-18.18-20.20-17.47-19.11-20.35-18.95
水稻田-10.34-12.23-10.82-9.63-10.65-11.96-9.64-11.47
湿地-11.20-10.08-10.99-9.34-8.94-8.50-18.98-20.85

新窗口打开| 下载CSV


表3   VH极化各地物不同时相后向散射系数

Tab.3  Backscattering coefficients of VH polarized objects in different time-phase(dB)

地类4月9日4月21日5月3日5月15日6月2日6月14日7月8日7月26日
莲田-22.27-21.68-19.78-20.08-16.88-15.82-16.91-14.26
林地-15.59-15.01-16.43-14.92-14.02-14.72-14.01-14.98
人工表面-16.07-15.04-16.07-14.98-14.34-15.28-14.36-15.35
水体-25.92-25.32-24.34-25.02-23.91-24.28-24.03-25.52
水稻田-18.82-19.92-20.44-19.63-17.98-19.36-18.03-17.39
湿地-17.27-15.42-16.74-14.96-14.30-14.90-24.79-26.81

新窗口打开| 下载CSV


图8

新窗口打开| 下载原图ZIP| 生成PPT

图8   VV极化各地物不同时相后向散射系数图

Fig.8   Backscattering coefficient diagrams of VV polarized objects in different time-phase


图9

新窗口打开| 下载原图ZIP| 生成PPT

图9   不同时相数据地类VV散射系数箱型图

Fig.9   Box diagram of VV scattering coefficient for different time-phase data


图10

新窗口打开| 下载原图ZIP| 生成PPT

图10   VH极化各地物不同时相后向散射系数

Fig.10   Backscattering coefficient diagrams of VH polarized objects in different time-phase


图11

新窗口打开| 下载原图ZIP| 生成PPT

图11   不同时相数据地类VH散射系数箱型图

Fig.11   Box diagram of VH scattering coefficient for different time-phase data


从VH极化各地物后向散射系数均值来看,各地物后向散射系数在-27~-14 dB之间,其中水体最低,林地和人工表面数值较为接近,且相对其他地类后向散射系数值较高。

从图11中可以看出,水体与其他地类具有较好的分离度,除7月8日和7月26日与湿地较为相似,其他几个时相数据均能较好分离。在VH极化数据中人工表面、林地的后向散射系数比较接近,难以分离。4月9日—5月3日数据莲田与水稻田分离性较低,原因是该时间段,水稻和莲都处于生长初期,地物主要受到水体的镜面反射影响,二者区分度不大,但是6月2日之后,莲田与水稻田的后向散射系数值逐渐拉大, 7月8日和7月26日二者具有较好的区分度。从单时相数据来看,6月14日、7月8日和7月26日水稻田与其他地类具有较好的区分度。

2.4 水稻田提取可分离性分析

2.4.1 Jeffries-Matusita距离

Jeffries-Matusita距离(J-M) 通过计算地类样本之间的特征距离判断地类的可分离程度,是分类过程中样本可分离程度评价的常用指数[29]。采用J-M距离计算地类之间的差异是一种具有通用性的定量方式,它不需要假定不同地物后向散射系数呈正态分布[29-30]。J-M距离计算结果值在0~2之间,值越大可分离性越好,若大于1.8则样本可分离性非常好,在1.4~1.8之间样本选择合格,样本可分离性小于1.4地类之间可分离性较差,建议重新选择样本,可分离性小于1则建议将地类样本进行合并[29-30]。

J-M距离j计算公式[29-30]为:

j=2(1-e-B)

式中B为某一特征维上的巴氏距离。

不同地类所选取的样本巴氏距离B计算公式[29-30]为:

B=18(m1-m2)22δ12+δ22+12lnδ12+δ122δ1δ2

式中:mi为地类特征平均值;δi2为某地类特征的方差[29-30]。

2.4.2 时相优选

为达到理想的分类精度,同时尽可能少地选用数据,降低数据获取和处理成本[31],研究对时相数据进行优选分析,确定水稻田识别的最佳时相。统计了不同时相水稻田与其他地物的可分离性,结果如表4所示。从单时相上来看,除了4月9日和7月26日的湿地,5月15日和6月14日的水体,7月26日的莲田具有可分离性外,其他时相水稻田与其他地物之间的J-M距离小于1.4,可分离性不佳,证明双极化数据单时相在地类分类过程中存在一定的局限性。

表4   单时相数据水稻田与各地物J-M距离

Tab.4  Single-phase data J-M value between paddy filed and other land types

地类水稻田
4月9日4月21日5月3日5月15日6月2日6月14日7月8日7月26日
莲田1.0690.5160.5270.2780.0780.8991.0011.455
林地0.8171.2361.2641.2801.2511.3901.1731.219
人工表面0.4550.9150.7150.8790.8571.0930.9950.930
水体1.2671.1061.3181.5951.0691.6041.2071.673
湿地1.5750.9431.1131.1741.1030.5361.1871.998

新窗口打开| 下载CSV


4月9号水稻刚刚出苗,易与杂草相混淆,7月26号处于第二季水稻返青时期,第二季水稻与第一季水稻种植面积有所差异。为了探测第一季水稻的种植范围,根据表4,从4月9日—7月8日第一季水稻生长过程中的SAR影像数据可分离程度最高的时相为6月14日。通过在6月14日垂直极化数据的基础上加上其他极化条件,每次加入的数据以J-M距离作为标准。表5给出了不同时相条件组合条件下水稻与其他地物的J-M距离。结果可以看出,增加时相数据均能增加地类的可分离性。这表明采用多时相数据可有效提高地类的分离度。

表5   多时相数据水稻田与各地物J-M距离

Tab.5  Multi-temporal data J-M value between paddy field and other land types

时相数据组合类型莲田林地人工表面水体湿地
(VV+VH)6.14+(VV/VH)6.14水稻田1.3051.6461.5231.5991.541
(VV+VH)6.14+(VV/VH)6.14+(VV+VH)4.09水稻田1.6441.8411.6721.9021.743
(VV+VH)6.14+(VV/VH)6.14+(VV+VH)4.21水稻田1.5011.9141.7861.7471.852
(VV+VH)6.14+(VV/VH)6.14+(VV+VH)5.03水稻田1.6351.9521.7431.8271.919
(VV+VH)6.14+(VV/VH)6.14+(VV+VH)5.15水稻田1.6061.9141.7131.8171.891
(VV+VH)6.14+(VV/VH)6.14+(VV+VH)6.02水稻田1.4931.8401.6591.7651.767
(VV+VH)6.14+(VV/VH)6.14+(VV+VH)7.08水稻田1.5241.7981.6651.7781.966
(VV+VH)6.14+(VV/VH)6.14+(VV+VH)7.26水稻田1.6751.8581.6821.9541.998

新窗口打开| 下载CSV


从表5中可以看出,加入7月26日影像数据水稻田与莲田可分离性最大达到1.675; 加入5月3日数据林地与水稻田可分离性最大,达到1.952; 加入4月21日数据人工表面与水稻田可分离性性最大,达到1.786; 加入7月26日数据水稻田与水体和湿地可分离性最大,分别达到1.954和1.998。

将4月21日、5月3日和7月26日数据与6月14日数据进行叠加生成一个包含不同时相的数据文件,计算水稻田与其他地类之间的J-M距离,结果如表6所示。

表6   水稻田与各地物J-M距离

Tab.6  J-M value between paddy field and other land types

类型莲田林地人工表面水体湿地
J-M距离1.8321.9791.8901.9751.999

新窗口打开| 下载CSV


采用以上数据进行组合的SAR数据水稻田与其他地类的可分离度均超过1.8,具有较好的分离度。研究区为实现各种地类有效区分,可选择6月14日、4月21日、5月3日和7月26日这4个时相数据。

2.4.3 水稻田识别精度评价

采用总体精度和F测度来分析和评价地物提取的精度,其中F测度包括准确率、召回率和F指数3个指标,准确率指在水稻田识别结果中水稻田被正确识别的比例,召回率指样本中水稻田被正确识别的比例,F指数指准确率P和召回率R的协调平均值[31],总体精度A是所有验证样本中被正确识别的比例,计算公式为:

P=TPTP+FP
R=TPTP+FN
F=2×P·RP+R
A=TP+TNTP+FP+TN+FN

式中: TP为水稻田样本被正确识别的个数; FP为非水稻田样本被识别为水稻田的个数; FN为水稻田样本被识别为非水稻田的个数; TN为非水稻田样本被正确识别的个数[31]。

3 结果分析与验证

3.1 时相优选前的水稻田识别

时相优选前,采用水稻生长所有物候期影像数据进行分类。分类方法包括随机森林法、最大似然法、支持向量机和神经网络法4种。

采用随机森林方法对影像进行分类,需要设置决策树个数n和节点分裂过程中抽取的特征个数m。特征变量m抽取时,默认选择总特征个数的算术平方根作为不同方案的特征个数; 理论上,决策树n的个数越多分类精度越高,但时间成本也越高。有研究发现,决策树n≥100时,所有方案的袋外(out of bag,OOB)误差逐渐收敛并趋于稳定[26,32]。因此选取n=100作为生成决策树的数量。最大似然法核函数选择径向基核函数(radial basis function),惩罚参数设置100,其他参数默认; 最大似然法概率阈值设置单值(single value),其他参数默认; 神经网络法训练阈值贡献0.9,训练率0.2,其他参数默认。各种方法执行分类后在ENVI软件中通过主要分析(majority)处理后导出,结果如图12所示。

图12

新窗口打开| 下载原图ZIP| 生成PPT

图12   时相优选前分类结果

Fig.12   Classification results before time phase optimization


随机森林、最大似然法和支持向量机对水稻田识别的总体精度分别为0.937,0.927和0.933,但最大似然法解译的人工表面比实际情况偏大,验证样本中3个水稻田样本被错分为人工表面,2个样本被错分为莲田,支持向量机2个样本被错分为人工表面,1个样本被错分为林地和2个样本被错分为莲田。神经网络法分类法被错分为莲田的面积最大,其水稻田验证样本中1个样本被错分为人工表面,6个样本被错分为莲田。

精度验证结果如表7所示。随机森林的准确率、召回率、F指数和总体精度均为最高,而神经网络分类结果各项指标都要比其他3种方法差,从水稻田识别的初步结果来看,在时相优选前,采用水稻生长各物候期的影像进行分类,随机森林分类方法所得到的结果总体最高,达到0.937。

表7   时相优选前结果精度比较

Tab.7  Accuracy comparison of classification results before time phase optimization

分类器准确率(P)召回率(R)F指数(F)总体精度(A)
随机森林0.9140.9600.9370.937
最大似然法0.9130.9500.9310.927
支持向量机0.8960.9500.9220.933
神经网络0.7420.9300.8250.770

新窗口打开| 下载CSV


3.2 时相优选结果分析

针对时相优选结果,同时使用随机森林、最大似然法、支持向量机和神经网络4种分类器进行分类并对比结果,解译结果如图13所示。解译结果总体上与图12保持一致。不过由于时相减少,优选后的解译结果有所变化。4种不同分类器对水稻田识别的结果见表8,随机森林、最大似然法和神经网络法解译结果的总体精度比时相优选前的结果均有所提高,说明时相优选后保留的影像数据在水稻田识别冗余信息减少。时相优选后,准确率最高的是神经网络法,达到0.936,随机森林法略低,为0.913; 从召回率上来看最大似然法、支持向量机和神经网络3种分类器在时相优选后均有提升,解译结果的漏检率减少,也说明时相优选后对水稻田识别有用的特征增多。从总体精度上来看,随机森林方法进行水稻田识别效果最佳,达到0.943,Kappa系数0.932。

图13

新窗口打开| 下载原图ZIP| 生成PPT

图13   时相优选后分类结果

Fig.13   Classification results after time phase optimization


表8   时相优选后分类结果精度比较

Tab.8  Accuracy comparison of classification results after time phase optimization

分类器准确率(P)召回率(R)F指数(F)总体精度(A)
随机森林0.9130.9500.9310.943
最大似然法0.9150.9700.9420.930
支持向量机0.8620.9700.9130.923
神经网络0.9360.9700.9530.907

新窗口打开| 下载CSV


4 结论与讨论

本文采用多时相Sentinel -1数据,多种分类方法提取南昌县蒋巷镇水稻种植范围,得到如下结论:

1)从J-M距离上来看,单时相Sentinel-1数据用于地物可分离度不高,水稻田容易与其他地类相混淆。Sentinel-1的VV和VH双极化单时相数据在不同的物候期,水稻田与其他地类J-M距离无法全部达到1.4以上。

2)通过增加时相数据可以增加地物的可分离性,选择最优时相数据不但可以降低成本、提升效率,还能降低数据冗余、提升数据解译精度。早稻物候期中采用三叶期、移栽期、孕穗期和二季晚稻移栽期SAR数据组合水稻田与莲田、林地、人工表面和水体的J-M距离分别可达到1.832,1.979,1.890和1.975。解译结果上来看,时相优选后,解译精度普遍得到了提升。

3)随机森林法、最大似然法、支持向量机和神经网络法在水稻田识别的潜力不同,其中无论是时相优选前还是时相优选后,随机森林分类方法的总体精度最高,时相优选后采用随机森林方法解译总体精度达到0.943,Kappa系数为0.932。

本文对比研究了随机森林、最大似然法、支持向量机和神经网络方法对水稻田进行提取,结果表明,随机森林方法提取异质性较强的水稻田具有较高的精度。但是在实际使用过程中不可避免地存在一些局限,如Sentinel-1 数据处理后的空间分辨率为20 m,在小尺度上可能难以满足应用需求,因此后续如果需要提取更加精细的地类分布范围,可以考虑将SAR数据与高空间分辨率数据进行融合再对数据进行提取。

参考文献

Lu Z, Wang D Q, Deng Z D, et al.

Application of red edge band in remote sensing extraction of surface water body:A case study based on GF-6 WFV data in arid area

[J]. Hydrology Research, 2021, 52(6):1526-1541.

DOI:10.2166/nh.2021.050      URL     [本文引用: 1]

吴炳方, 张淼, 曾红伟, 等.

全球农情遥感速报系统20年

[J]. 遥感学报, 2019, 23(6):1053-1063.

[本文引用: 1]

Wu B F, Zhang M, Zeng H W, et al.

Twenty years of crop watch:Progress and prospect

[J]. Journal of Remote Sensing, 2019, 23(6):1053-1063.

[本文引用: 1]

杨知. 基于极化SAR的水稻物候期监测与参数反演研究[D]. 北京: 中国科学院大学, 2017.

[本文引用: 1]

Yang Z. Rice phenology estimation and parameter retrieval based on polarimetric synthetic aperture radar[D]. Beijing: The University of Chinese Academy of Sciences, 2017.

[本文引用: 1]

Perros N, Kalivas D, Giovos R.

Spatial analysis of agronomic data and UAV imagery for rice yield estimation

[J]. Agriculture, 2021, 11(9):809.

DOI:10.3390/agriculture11090809      URL     [本文引用: 1]

陈仲新, 任建强, 唐华俊, 等.

农业遥感研究应用进展与展望

[J]. 遥感学报, 2016, 20(5):748-767.

[本文引用: 1]

Chen Z X, Ren J Q, Tang H J, et al.

Progress and perspectives on agricultural remote sensing research and applications in China

[J]. Journal of Remote Sensing, 2016, 20(5):748-767.

[本文引用: 1]

陈怀亮, 李颖, 张红卫.

农作物长势遥感监测业务化应用与研究进展

[J]. 气象与环境科学, 2015, 38(1):95-102.

[本文引用: 1]

Chen H L, Li Y, Zhang H W.

Operational application and research review of crop growth monitoring with remote sensing

[J]. Meteorological and Environmental Sciences, 2015, 38(1):95-102.

[本文引用: 1]

Zhu L H, Liu X N, Wu L, et al.

Detection of paddy rice cropping systems in southern China with time series Landsat images and phenology-based algorithms

[J]. GIScience & Remote Sensing, 2021, 58(5):1-23.

[本文引用: 1]

Cao J J, Cai X L, Tan J W, et al.

Mapping paddy rice using Landsat time series data in the Ganfu Plain irrigation system,Southern China,from 1988-2017

[J]. International Journal of Remote Sensing, 2021, 42(4):1556-1576.

DOI:10.1080/01431161.2020.1841321      URL     [本文引用: 1]

Li R Y, Xu M Q, Chen Z Y, et al.

Phenology-based classification of crop species and rotation types using fused MODIS and Landsat data:The comparison of a random-forest-based model and a decision-rule-based model

[J]. Soil and Tillage Research, 2021, 206:104838.

DOI:10.1016/j.still.2020.104838      URL     [本文引用: 1]

Chandna P, Mondal S.

Analyzing multi-year rice-fallow dynamics in Odisha using multi-temporal Landsat-8 OLI and Sentinel-1 Data

[J]. GIScience & Remote Sensing, 2020, 57(4):431-449.

[本文引用: 1]

杨辉山, 谢萍.

华南地区季度性全覆盖卫星影像获取方案优化

[J]. 测绘与空间地理信息, 2019, 42(5):159-162.

[本文引用: 1]

Yang H S, Xie P.

Optimized solution of the high frequency and full coverage of satellite images acquisition in southern China

[J]. Geomatics & Spatial Information Technology, 2019, 42(5):159-162.

[本文引用: 1]

Gasnier N, Denis L, Fjrtoft R, et al.

Narrow river extraction from SAR images using exogenous information

[J]. IEEE Journal of Selected Topics in Applied Earth Observations and Remote Sensing, 2021, 14:5720-5734.

DOI:10.1109/JSTARS.2021.3083413      URL     [本文引用: 1]

Chang L, Chen Y T, Wang J H, et al.

Rice-field mapping with sentinel-1A SAR time-series data

[J]. Remote Sensing, 2021, 13(1):103.

DOI:10.3390/rs13010103      URL     [本文引用: 1]

Wu X X, Washaya P, Liu L, et al.

Rice yield estimation based on spaceborne SAR:A review from 1988 to 2018

[J]. IEEE Access, 2020(8):157462-157469.

[本文引用: 1]

Son N T, Chen C F, Chen C R, et al.

A phenological object-based approach for rice crop classification using time-series Sentinel-1 synthetic aperture Radar (SAR) data in Taiwan

[J]. International Journal of Remote Sensing, 2021, 42(7):2722-2739.

DOI:10.1080/01431161.2020.1862440      URL     [本文引用: 1]

Graman M, Kaliaperumal R, Pazhanivelan S, et al.

Rice area estimation using parameterized classification of sentinel 1A SAR data

[J]. ISPRS - International Archives of the Photogrammetry,Remote Sensing and Spatial Information Sciences, 2019, XLII-3/W6:141-147.

[本文引用: 2]

Zhu H L, Wang W Y, Leung R.

SAR target classification based on Radar image luminance analysis by deep learning

[J]. IEEE Sensors Letters, 2020, 4(3):1-4.

[本文引用: 1]

Bianchi F, Espeseth M, Borch N.

Large-scale detection and categorization of oil spills from SAR images with deep learning

[J]. Remote Sensing, 2020, 12(14):2260.

DOI:10.3390/rs12142260      URL     [本文引用: 1]

邓刚, 唐志光, 李朝奎, 等.

基于MODIS时序数据的湖南省水稻种植面积提取及时空变化分析

[J]. 国土资源遥感, 2020, 32(2):177-185.doi: 10.6046/gtzyyg.2020.02.23.

[本文引用: 1]

Deng G, Tang Z G, Li C K, et al.

Extraction and analysis of spatiotemporal variation of rice planting area in Hunan Province based on MODIS time-series data

[J]. Remote Sensing for Land and Resources, 2020, 32(2):177-185.doi: 10.6046/gtzyyg.2020.02.23.

[本文引用: 1]

张悦琦, 李荣平, 穆西晗, 等.

基于多时相GF-6遥感影像的水稻种植面积提取

[J]. 农业工程学报, 2021, 37(17):189-196.

[本文引用: 1]

Li Y Q, Li R P, Mu X H, et al.

Extraction of paddy rice planting areas based on multi-temporal GF-6 remote sensing image

[J]. Transactions of the Chinese Society of Agricultural Engineering, 2021, 37(17):189-196.

[本文引用: 1]

Woodhouse I. Introduction to microwave remote sensing[M]. USA: CRC Press, 2017:157-158.

[本文引用: 3]

Xu Z.

Wavelength-resolution SAR speckle model

[J]. IEEE Geoscience and Remote Sensing Letters, 2022, 19:1-5.

[本文引用: 1]

Hu C B, Laurent L F, Kuang G Y.

Ship discrimination using polarimetric SAR data and coherent time-frequency analysis

[J]. Remote Sensing, 2013, 5(12):6899-6920.

DOI:10.3390/rs5126899      URL     [本文引用: 1]

Luo Y J, Zhao Y S, Li X W, et al.

Research and application of multi-angle polarization characteristics of water body mirror reflection

[J]. Science in China(Series D:Earth Sciences), 2007, 50(6):946-952.

[本文引用: 2]

Jang J W, Kim T Y, Lim S B.

RF Interference analysis and verification in the synthetic aperture Radar satellite system

[J]. Aerospace Engineering and Technology, 2009, 8(1):187-196.

[本文引用: 2]

吴琳琳, 李晓燕, 毛德华, 等.

基于遥感和多源地理数据的城市土地利用分类

[J]. 自然资源遥感, 2022, 34(1):127-134.doi: 10.6046/zrzyyg.2021061.

[本文引用: 4]

Wu L L, Li X Y, Mao D H, et al.

Urban land use classification based on remote sensing and multi-source geographic data

[J]. Remote Sensing for Natural Resources, 2022, 34(1):127-134.doi: 10.6046/zrzyyg.2021061.

[本文引用: 4]

桂预风, 李振平, 万爽, 等.

收敛性随机森林模型及其遥感应用

[J]. 数学的实践与认识, 2015, 45(18):207-212.

[本文引用: 1]

Gui Y F, Li Z P, Wan S, et al.

Convergent random forests model and its application in remote sensing

[J]. Mathematics in Practice and Theory, 2015, 45(18):207-212.

[本文引用: 1]

吕岚.

基于随机森林的快速兴趣点检测

[J]. 自动化技术与应用, 2016, 35(9):95-100.

[本文引用: 1]

Lyu L.

A fast interest point detector based on random forest

[J]. Techniques of Automation and Applications, 2016, 35(9):95-100.

[本文引用: 1]

韩冰冰, 陈圣波, 曾庆鸿, 等.

基于J-M距离的多时相Sentinel-1农作物分类

[J]. 科学技术与工程, 2020, 20(17):6977-6982.

[本文引用: 6]

Han B B, Chen S B, Zeng Q H, et al.

Time-series classification of Sentinel-1 data based on J-M distance

[J] Science Technology and Engineering, 2020, 20(17):6977-6982.

[本文引用: 6]

刘警鉴, 李洪忠, 华璀, 等.

基于Sentinel-1A数据的临高县早稻面积提取

[J]. 国土资源遥感, 2020, 32(1):191-199.doi: 10.6046/gtzyyg.2020.01.26.

[本文引用: 5]

Liu J J, Li H Z, Hua C, et al.

Extraction of early paddy rice area in Lingao county based on sentinel-1A data

[J]. Remote Sensing for Land and Resources, 2020, 32(1):191-199.doi: 10.6046/gtzyyg.2020.01.26.

[本文引用: 5]

马腾耀, 肖鹏峰, 张学良, 等.

基于特征优选的GF-3全极化数据积雪识别

[J]. 遥感技术与应用, 2020, 35(6):1292-1302.

[本文引用: 3]

Ma T Y, Xiao P F, Zhang X L, et al.

Recognition of snow cover based on features selection in GF-3 fully polarimetric data

[J]. Remote Sensing Technology and Application, 2020, 35(6):1292-1302.

[本文引用: 3]

张颖, 高倩倩.

基于随机森林分类算法的巢湖水质评价

[J]. 环境工程学报, 2016, 10(2):992-998.

[本文引用: 1]

Zhang Y, Gao Q Q.

Water quality evaluation of Chaohu Lake based on random forest method

[J]. Chinese Journal of Environmental Engineering, 2016, 10(2):992-998.

[本文引用: 1]

/

〈 〉

京ICP备05055290号-2
版权所有 © 2015 《自然资源遥感》编辑部
地址:北京学院路31号中国国土资源航空物探遥感中心 邮编：100083
电话:010-62060291/62060292 E-mail:zrzyyg@163.com
本系统由北京玛格泰克科技发展有限公司设计开发