国土资源遥感, 2018, 30(4): 20-27 doi: 10.6046/gtzyyg.2018.04.04

基于SVM+SFS策略的多时相紧致极化SAR水稻精细分类

国贤玉1, 李坤,2, 王志勇1, 李宏宇3, 杨知4

1. 山东科技大学测绘科学与工程学院,青岛 266590

2. 中国科学院遥感与数字地球研究所,北京 100101

3. 中国地质大学(北京)地球科学与资源学院,北京 100083

4. 中国电力科学研究院输变电工程研究所,北京 100055

Fine classification of rice with multi-temporal compact polarimetric SAR based on SVM+SFS strategy

GUO Xianyu1, LI Kun,2, WANG Zhiyong1, LI Hongyu3, YANG Zhi4

1. College of Geomatics, Shandong University of Science and Technology, Qingdao 266590, China

2. Institute of Remote Sensing and Digital Earth, Chinese Academy of Sciences, Beijing 100101, China

3. School of Earth Science and Resources, China University of Geosciences(Beijing), Beijing 100083, China

4. Institute of Transmission and Transformation Engineering, China Electric Power Research Institute, Beijing 100055, China

通讯作者: 李坤(1984-),女,副研究员,主要从事雷达遥感、农业遥感等方面的研究。Email:likun@radi.ac.cn

责任编辑: 陈理

收稿日期: 2017-06-23   修回日期: 2017-09-3   网络出版日期: 2018-12-15

基金资助: 国家自然科学基金青年基金项目“基于紧致极化SAR水稻物候期反演方法研究”.  41401427
国家自然科学基金面上基金项目“复杂散射机制场景的SAR图像认知方法研究”.  61471358
国家自然科学基金重点基金项目“可控环境下多层介质目标微波特性全要素测量与散射机理建模”共同资助.  41431174

Received: 2017-06-23   Revised: 2017-09-3   Online: 2018-12-15

作者简介 About authors

国贤玉(1991-),男,硕士研究生,主要从事雷达遥感、农业遥感等方面的研究。Email:guoxianyu1@126.com。 。

摘要

种类和种植方式的不同会导致水稻长势、产量的差异。精细区分不同水稻品种与种植方式,能够为水稻长势监测与估产提供更精准的信息。紧致极化SAR(compact polarimetry synthetic aperture Radar,CP-SAR)作为新一代对地观测SAR系统的重要发展趋势之一,同时兼具相对丰富的极化信息和较大的幅宽,为大范围水稻精细制图提供了可能。本研究首先利用RADARSAT-2全极化SAR数据模拟CP-SAR数据,并提取了22个CP-SAR特征参数; 然后,针对CP-SAR多维特征信息,引入基于支持向量机和序列前进搜寻策略(support vector machine + sequential forward selection,SVM + SFS)的特征选择方法,构建基于决策树和SVM的水稻精细分类方法,得到了水稻精细分类的最优特征子集。实验结果表明,基于决策树的水稻精细分类方法可以获得较好的分类结果,总体精度达92.57%,Kappa系数达0.896,与全部特征参数进行分类的结果相比,总体精度高1.2%,Kappa系数大0.016。

关键词: CP-SAR ; SVM + SFS ; 水稻田 ; 决策树分类 ; 多时相

Abstract

Different types and planting methods can result in the discrepancy of rice growth and yield. It is of great importance to provide accurate growth vigor information for rice growth monitoring and estimation using fine distinction of different rice varieties and planting methods. As a new type of SAR sensor, compact polarimetry synthetic aperture Radar (CP-SAR) provides the possibility for the fine mapping of paddy land with abundant polarimetric information and large width. In this study, the authors firstly used RADARSAT-2 fully polarimetric SAR data to simulate CP-SAR data and extracted 22 types of feature parameters. In addition, on the basis of the multi-dimensional feature information CP - SAR data, the support vector machine and sequential forward selection (SVM + SFS) strategy were performed for feature selection, and the optimal feature subset of paddy land fine classification was obtained. Moreover, the decision tree and SVM method were used for paddy land fine classification based on feature subset. The results show that paddy land fine classification method based on decision tree can achieve better classification results. The overall classification precision and Kappa coefficient of optimal feature subset respectively are 92.57% and 0.896, which are higher than those of the set of all feature parameters by improving 1.2% in overall classification precision and 0.016 in Kappa coefficient.

Keywords: CP-SAR ; SVM+SFS ; paddy land ; decision tree classification ; multi-temporal

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

本文引用格式

国贤玉, 李坤, 王志勇, 李宏宇, 杨知. 基于SVM+SFS策略的多时相紧致极化SAR水稻精细分类. 国土资源遥感[J], 2018, 30(4): 20-27 doi:10.6046/gtzyyg.2018.04.04

GUO Xianyu, LI Kun, WANG Zhiyong, LI Hongyu, YANG Zhi. Fine classification of rice with multi-temporal compact polarimetric SAR based on SVM+SFS strategy. REMOTE SENSING FOR LAND & RESOURCES[J], 2018, 30(4): 20-27 doi:10.6046/gtzyyg.2018.04.04

0 引言

水稻是世界三大粮食作物之一,为我国一半以上的人口提供粮食来源。种类和种植方式的不同导致水稻长势、产量存在一定差异,传统的水稻制图(区分水稻和非水稻)已经难以满足高精度农业应用的需求。因此实现水稻精细制图,区分不同水稻品种与种植方式,为水稻长势监测提供更精准的信息,对于现代农业的发展具有重要意义。

紧致极化SAR(compact polarimetry synthetic aperture Radar,CP-SAR)降低了系统复杂度与能耗,缩小了传感器体积,已成为新一代对地观测SAR系统的重要发展趋势之一[1]。与全极化SAR相比,CP-SAR不仅能够保持丰富的极化信息,还能实现更大的幅宽与入射角范围。近年来,CP-SAR相关研究主要集中在3方面: ①CP-SAR系统接发模式研究[2,3]; ②CP-SAR模拟与伪极化(pseudo-quad-pol,PQ)SAR重建方法研究[4,5]; ③CP-SAR应用研究,如信息提取[6]、作物分类[7]、森林参数反演[8]、海冰和溢油[9,10]等。虽然目前基于CP-SAR的应用研究覆盖面很广,但还不够深入,以农业应用为例,大多数研究都集中在简单的作物制图上,对于种植方式和种类的区分研究很少。

目前SAR水稻制图方法主要依据有3类: ①后向散射特性的时相变化规律[11,12]; ②不同极化后向散射特性的差异[13]; ③全极化散射机理特点[14,15]。前2类方法都只利用后向散射强度信息,不包含雷达回波的相位信息。第3类方法精度高,普适性较强,对数据时相的要求也较低。虽然全极化SAR在水稻制图中具有较大优势,但全极化系统的脉冲重复频率是单双极化的2倍,相应的幅宽也小,限制了大范围水稻制图的应用。因此,在同时兼顾制图精度与面积的情况下,CP-SAR是最佳选择之一。2013年,Brisco等[16]基于CP-SAR开展水稻制图研究,对比分析了单双极化、CP-SAR与全极化SAR的制图效果,结果表明CP-SAR在水稻制图中的应用效果可与全极化相媲美,远优于单、双极化数据; 2015年,Uppala等[17]基于RISAT-1卫星CP-SAR数据利用监督分类进行水稻识别,得到了较高的制图精度。这些研究表明了CP-SAR在水稻制图中的应用潜力,但集中于区分水稻和非水稻,对于水稻种类以及种植方式的区分研究不足。

鉴于此,以江苏金湖地区为研究区,开展CP-SAR水稻精细制图方法研究。针对插秧籼稻/粳稻、撒播粳稻3类水稻田,考虑水稻植株分布特征、生理结构特点以及下垫面的影响,研究分析其CP-SAR响应特征以及时相变化规律,在此基础上,针对CP-SAR多维特征信息,引入基于支持向量机和序列前进搜寻(support vector machine and sequential forward selection,SVM + SFS)[18]策略的特征选择方法,构建基于决策树和SVM的水稻精细分类方法。

1 研究区概况与数据源

研究区位于江苏金湖(E118°41'34″~119°16'27″,N33°17'05″~33°56'39″),属于亚热带季风气候区,地势平坦,地块规则。该区水稻一年一熟(6—11月)。水稻种类为籼稻和粳稻,播种方式分为插秧和撒播,故水稻田可分为插秧籼稻田(TH)、撒播籼稻田、插秧粳稻田(TJ)和撒播粳稻田(DJ)4类。由于该区几乎没有撒播籼稻田,因此主要针对TH,TJ和DJ这3类(图1),开展精细制图方法研究。

图1

图1   3类水稻在不同物候期的特征

Fig.1   Three kinds of paddy land characteristics in different phonological stages


在研究区获取了9景RADARSAT-2精细全极化SAR数据,方位向和距离向空间分辨率分别为5.2 m和7.6 m。由于封行之前3类水稻田差异相对较大,因此,选择对应时段的SAR数据进行水稻精细分类方法研究,获取日期分别为2012年6月27日、7月11日和7月21日。首先基于3个时相的全极化SAR数据模拟CP-SAR数据。模拟数据为圆周极化发射线性极化接受模式(circular transimit and linear receive,CTLR),发射右旋圆(R)极化、接收水平(H)和垂直(V)极化[19],空间分辨率为30 m,噪声水平为-25 dB(图2)。获取SAR数据的同时,开展了地面实验,采集了水稻种类、种植方式和物候等信息,并利用高精度GPS获取了41块水稻样田的矢量数据,其中包括24块TH、6块TJ、11块DJ,还选择了8块水体和10块城镇建筑。

图2

图2   CP-SAR模拟数据在不同极化通道的假彩色合成影像
(CP-SAR RR(R),RV(G),RH(B)假彩色合成)

Fig.2   Color synthetic images of CP-SAR data in different polarization channels


2 研究方法

研究流程主要包括CP-SAR数据模拟与特征参数提取、数据预处理、基于SVM + SFS的CP-SAR特征参数优选以及基于优选特征利用决策树和SVM方法进行水稻田精细分类,具体技术流程如图3所示。

图3

图3   技术路线

Fig.3   Flow chart of technology


2.1 数据预处理

基于CP-SAR模拟数据,根据特征参数定义,提取22个CP特征参数(表1)。然后对特征参数进行辐射定标、几何纠正、研究区裁剪和斑点噪声滤波等预处理。通过比较选择Frost滤波方法,以7×7窗口进行降噪处理。在此基础上,基于地面样方,提取不同类型水稻田、水体和城镇建筑的CP-SAR特征参数。

表1   提取的22个CP特征参数

Tab.1  Twenty two CP characteristic parameters

符号参数名称参数物理意义公式
g=[g0,g1,g2,g3]TStokes矢量表征散射回波强度和极化状态文献[19]公式2
σ0RH
σ0RV
σ0RL
σ0RR
RH,RV,RL,RR极化后向散射系数表征目标在RH,RV,RL,RR极化的回波强度σ0RH=[1,1,0,0]·g
σ0RV=[1,-1,0,0]·g
σ0RL=[1,0,0,1]·g
σ0RR=[1,0,0,-1]·g
δRH和RV相位角RH和RV极化散射回波之间的相对相位角δ=arctan(g3/g2)
m-δ_db
m-δ_vol
m-δ_s
m-δ分解二次散射、体散射和面散射表征地物二次散射、体散射和面散射文献[20]公式7
m-χ_db
m-χ_vol
m-χ_s
m-χ分解二次散射体散射和面散射表征地物二次散射、体散射和面散射文献[19]公式5
μC圆极化比Stokes向量派生参数μC=(g0-g3)/(g0+g3)
μ一致性系数与目标散射机理有关μ=2Im<sRHs*RV><sRHs*RV>+<sRHs*RV>,
式中:sRHsRV分别为极化散射矩阵元素; Im
表示取复数的虚部; <>表示内积运算; *表示
复共轭
m极化度表征目标散射回波去相关程度m=g12+g22+g32g0,0&lt;m&lt;1
ρRH-RV相关系数表征散射矩阵元素sRHsRV的相关性ρ=&lt;|sRHs*RV|&gt;&lt;sRHs*RV&gt;+&lt;|sRHs*RV|&gt;
α平均散射角与目标散射机理有关α=12arctan[g12+g22/(±g3)]
Hi香农熵表征相干矩阵的强度Hi=3lgπetr[T]3,
式中:T为3×3相干矩阵;
tr[T]表示相干矩阵T的迹
符号参数名称参数物理意义公式
Hp香农极化度表征Bakarat(pT)的极化度,
pT=1-27det[T]tr[T]3
式中det[T]表示T的行列式值
Hp=lg27det[T]tr[T]3

新窗口打开| 下载CSV


2.2 SVM + SFS策略特征选择方法

为了充分挖掘CP-SAR多维特征信息,同时保证分类方法的简洁性,引入基于SVM + SFS的特征选择方法,对22个CP-SAR参数进行优选。把每一特征参数看作由一个向量和一个标记组成,即 Di=(xi,yi), x=[x1,,xi,,xn]为训练数据向量,n为训练数据个数, yi为分类标记( yi取-1或1)。定义函数和超平面分别为

g(xi)=<w,x>+b,i[1,n]

<w,x>+b=0

式中: w为系数向量,其维度为n; b为常数变量。若使分类数据被超平面分成2类,超平面必须满足 yi(<w,x>)1。SVM思想是使所求最优超平面能够具有最大的分类间隔,分类间隔 δi表示为

δi=1wg(xi),i[1,n]

式中: w为向量 w的范数; g(xi)g(xi)的绝对值。这等同于求二次规划问题,即

minΦ(w)=12<w,w>

yi(<w,x>+b)1,i[1,n]

引入Lagrange算子 α*,令 α*0,满足式(6)有唯一解,即

αi*[yi(<w,x>+b*)-1]=0,i[1,n]

式中 b*最优化的常数变量。当样本点到超平面距离为最短距离,则 yi(<w,x>+b)=1α*0,否则 yi(<w,x>+b)>1α*=0α*=0的样本称为支持向量(support vector,SV),样本的总个数称为SV个数(number of SV,NSV)。在SVM分类算法中,可分性的优劣就是由NSV判断,NSV越小,可分性越好。

以区分3种水稻田为目的,利用SVM+SFS特征选择方法步骤如下: ①原始特征集 F={fj,(j=1,2,,22)},依次计算特征参数 fj的NSV; ②每一特征参数均得到3个NSV,对该参数的3个NSV求平均值,得 N-jSV=NSV(fj); ③将22个特征参数的NSV按升幂排列, F'为排序特征集,即 F'={f'j,(j=1,2,,22)}; ④设结果特征集 F″,f'1为第一个结果特征即 f1; ⑤ f'j(j=2,3,,22)作为待选特征参数; ⑥计算NSV( F″f'j)和NSV( F″),若NSV( F″f'j)<NSV( F″),则将 f'j选为结果特征参数,否则舍弃 f'j; ⑦遍历所有待选特征参数得到结果特征集 F″

除了3类水稻田的最优特征,利用上述方法还选出了区分水稻与非水稻的最优特征。

2.3 3类水稻田CP-SAR响应规律

面向水稻田精细分类,利用SVM + SFS方法,优选出的CP特征参数如表2所示。图4给出了3类水稻田在优选参数上的差异,且将优选特征参数分为2类: ①强度特征参数; ②非强度特征参数。

表2   利用SVM + SFS方法优选的CP-SAR特征参数

Tab.2  Optimal CP-SAR characteristic parameters by SVM + SFS

影像日期CP-SAR特征参数
2012年6月27日μ,δ,α,m-δ_db,g0,g3,σ0RL,σ0RV,σ0RH,σ0RR,m-δ_vol,m-χ_vol
2012年7月11日m-χ_s,σ0RL,σ0RV,σ0RH,m-δ_vol,Hi,g1,m-χ_vol
2012年7月21日δ,σ0RL,σ0RV,σ0RH,σ0RR,g0,m-χ_s,μC

新窗口打开| 下载CSV


图4

图4   3种水稻强度和非强度极化特征参数

Fig.4   Intensity and non-intensity polarization characteristic parameters about three kinds of paddy land


相对于TH和TJ,DJ水稻植株密度更大,因此其后向散射和体散射都比较大; 而TH和TJ的下垫面为水面,引起镜面反射使其后向散射和体散射较小,这导致DJ与TH,TJ的后向散射和体散射差异较大。由于 σ0RHσ0RV主要来自体散射的去极化作用,因此DJ的 σ0RHσ0RV大于TH和TJ(如图4(a)所示),差值约为0.8 dB; 而 g0,g1与后向散射密切相关,2参数对于区分DJ与TH,TJ有较大贡献。TH下垫面为水面,且籼稻幼苗植株更高且粗壮,下垫面与植株垂直结构更容易形成二面角,因此TH的二次散射更强; TJ植株高度较小,DJ下垫面为土壤,因此TJ和DJ的二次散射相对较弱。由于二次散射在RR上的响应较强,因此TH的 σ0RR大于TJ和DJ,差值约为1.3 dB; m-χ_dbm-δ_db表征二次散射的强度,因此TH的 m-χ_dbm-δ_db强度值大于TJ和DJ,差值约为3.8 dB。所以 σ0RR, m-χ_dbm-δ_db对于区分TH与TJ,DJ有较大贡献。DJ下垫面为土壤,其面散射贡献最大; TJ植株相对弱小,下垫面粗糙面散射贡献较大,TH植株相对高而粗壮,面散射最弱,由于 σ0RL, g3, m-χ_s,m-δ_s与面散射密切相关,因此3类水稻田对应的这4个参数差异较大(如图4(a)所示)。以 g3为例,其差值约为3 dB,对于区分3类水稻田贡献较大。TJ植株密度相对较小,而且粳稻植株相对弱小,因此其体散射相对于DJ和TH较小。 m-δ_volm-χ_vol表征地物体散射,故TJ的这2个参数小于DJ和TH,其差值约为0.7 dB,对于区分TJ与DJ,TH贡献较大。

Hi表征散射机制的复杂程度,由于体散射更为复杂,因此体散射贡献越大 Hi越大。通过前面3类水稻田的散射机理分析,DJ的体散射贡献最大,TH次之,TJ最小,由图4(b)可以看出,DJ的 Hi大于TH且远大于TJ,因此, Hi对于区分DJ和TJ贡献较大。 μα都与散射机理密切相关, μ从大到小分别表示面散射、体散射和二次散射; 而 α反之。因此DJ的 μ值大于TJ和TH,而DJ的 α值小于TJ和TH。 μC也与目标的散射机理密切相关,其值与面散射的贡献成反比,DJ对应的 μC值小于TJ和TH,对于区分TH与DJ贡献较大。

2.4 基于CP-SAR优选特征的水稻精细分类

基于SVM + SFS方法优选CP-SAR特征,分别采用决策树和SVM方法进行水稻精细分类。另外,将3类水稻田、城镇建筑和水体样方分为训练和验证样本2部分,TH、水体和城镇建筑的训练和验证样本各占一半,二者之间没有重叠。由于TJ和DJ的样本数较少,训练和验证样本之间约有30%的重叠。

2.4.1 决策树分类

首先利用CP-SAR优选特征区分水稻与非水稻,再进行3类水稻田的区分,最终实现精细分类,决策树分类如图5所示。图中变量的数字后缀代表影像获取日期。

图5

图5   3类水稻田的分类决策树

Fig.5   Classification decision tree of three kinds of paddy land


研究区非水稻区域主要包括城镇建筑和水体等,水稻与水体、建筑的二次散射贡献差异很大(图6),而RR极化对二次散射敏感,因此首先根据 σ0RR,区分水稻和非水稻。水体的 m-δ_db约为-35 dB,小于其他非水稻区域,因此利用 m-δ_db区分水体; 最后再利用 m-δ_dbσ0RR将城镇建筑与其他非水稻区分开。针对3类水稻田,先利用6月27日(幼苗期) m-χ_db_0627, μC_0627m-δ_vol_0627区分不同种植方式,即DJ与TH,TJ。因为撒播田下垫面为土壤,且植株矮小,二次散射比插秧田弱,而由于植株密度较大,其体散射较弱大于TJ,小于TH。另外,针对TH与TJ可分性较弱,且二者种植方式相同,田块结构相似,只能依靠水稻植株形态差异进行区分。7月11日,TH在RR上的响应较强,21日由于冠层密度增加,衰减增大,TH发生二次散射的能量减少,在RR上的响应减弱; 而TJ刚好与之相反,因此利用二者在2个时相上的差异来实现区分。

图6

图6   不同地物极化特征参数

Fig.6   Polarimetric parameters in different surface features


2.4.2 SVM分类

基于CP-SAR优选特征,利用SVM进行分类。选择径向基核函数(radial basis function,RBF),其Gamma值为输入图像波段的倒数,惩罚参数为100; 分级处理等级为0,以原图像空间分辨率进行分类处理; 分类概率阈值为0。

3 结果与分析

通过设计4组对比实验进行结果分析: ①利用6月27日12个CP优选参数进行SVM分类,并与全部22个参数SVM分类结果进行比较; ②考虑时相信息,利用3个时相28个优选参数进行SVM分类,并与全部66个参数分类结果进行比较; ③利用决策树方法进行TH,DJ与非水稻区的区分; ④利用决策树进行3类水稻田与非水稻区的区分。最后利用验证样本对分类结果进行精度评价(表3)。可以看出,水体和城镇建筑分类效果较好,生产者精度和用户精度均在90%以上,不同方法、不同时相组合对应的分类结果差异不大。

表3   2种分类方法的分类精度比较

Tab.3  Comparison of classification accuracy between the two methods

方法分类参数THTJDJ水体城镇建筑总体精
度/%
Kappa系数
生产者精
度/%
用户精
度/%
生产者精
度/%
用户精
度/%
生产者精
度/%
用户精
度/%
生产者精
度/%
用户精
度/%
生产者精
度/%
用户精
度/%
SVMT1-12-399.4241.901.1858.2364.5678.3899.8299.8792.0899.7583.720.778
T1-22-389.9049.780.5710078.7466.3210010093.7299.8485.300.798
T3-28-310059.6324.1695.2091.4485.6610010095.1310091.390.880
T3-66-398.8561.4046.4873.5873.9889.1810010095.6910091.380.880
决策树T3-28-294.7387.8995.4492.4310010095.8999.6297.440.962
T3-28-396.2569.5745.7468.0686.3988.5310010095.3099.0392.570.896

新窗口打开| 下载CSV


对于水稻来说,单一时相12个优选参数SVM分类,TH平均精度约为70.6%; TJ的生产者和用户精度都很低。多时相28个优选参数SVM分类,3类水稻的精度都有所提高,因此多时相对区分水稻种类、种植方式具有一定贡献; 但是TJ生产者精度只有24.16%,即大部分TJ被错分成TH,说明水稻种类的区分能力仍然不高。多时相全部66个参数SVM分类,总体精度为91.38%,Kappa系数为0.880,与多时相28个优选参数SVM分类结果相近,可见对于水稻精细分类,基于SVM + SFS策略优选的28个特征参数能够与全部66个参数达到同样效果,避免了数据赘余、提高了运算效率。

利用多时相28个优选参数,进行决策树方法,区分TH和DJ,总体精度为97.44%,Kappa系数为0.962; 区分3类决策树分类总体精度达到92.57%,Kappa系数达到0.896。与优选参数SVM分类对比总体精度提高1%~9%,Kappa系数提高0.01~0.12。TJ生产者精度为45.74%,比SVM分类精度提高了0~40%。总体来看,SVM + SFS策略优选参数决策树分类要优于SVM分类,并且分类速度更快。

从分类精度来看,TJ生产者精度较低,区分效果不好主要有以下几方面原因: ①由于播种方式的不同,幼苗期插秧稻田种植稀疏,植株成行成垄,而撒播稻田植株较稠密,在雷达响应上表现差异性大,因此DJ容易与TH,TJ区分。幼苗期TH和TJ这2类水稻具有水稻共性,并且幼苗期水稻植株小,导致在雷达响应上表现差异性小; 分蘖期和拔节期2类水稻植株表现出差异性,但随着植株生长,植株间的缝隙减小,这种差异性又淹没在水稻群体中,导致在雷达响应上差异性小,使得TJ区分效果不好; ②研究区TJ种植面积少,在研究区地面获取的样方也少,影响TJ生产者精度; ③本研究使用CP-SAR模拟数据,空间分辨率为30 m,噪声水平为-25 dB,空间分辨率和噪声水平与真实SAR数据(以RADARSAT-2全极化为例,空间分辨率8 m,噪声水平约为-32 dB)存在一定的差异。

采用SVM和决策树分类方法,3个时相28个参数分类结果如图7所示。

图7

图7   3个时相28个参数分类结果

Fig.7   Classification results of 28 parameters at three periods


图7(a)中可看出,城镇建筑和水体被明显分出,这与城镇建筑、水体与水稻的散射特性差异性大有关。除水体和城镇建筑外,水稻田分为3类,TH多分布在东南区,TJ和DJ多分布在西北部,以TH分布最为广泛。这与研究区实际水稻种植分布现状基本相符; 在图7(b)中,利用决策树分类比SVM分类效果更细,将村庄道路也区分出来,从整体来看,依然是TH分布在金湖地区东南部,TJ和DJ分布在西北部,城镇多分布在研究区南部。

4 结论

利用CP-SAR模拟数据提取多维特征信息,引入基于SVM + SFS的特征选择方法,构建了基于决策树和SVM的水稻精细分类方法,为水稻长势监测与估产提供了更精准的信息。具体结论如下:

1)利用多时相CP-SAR模拟数据,分析了不同种植方式、不同品种的3类水稻田的CP-SAR响应特征、散射机理及其时相变化规律。

2)针对CP-SAR多维特征参数,引入基于SVM + SFS的特征选择方法,建立了面向水稻田精细分类的CP-SAR最优特征集,并结合物理意义分析了这些特征在不同水稻田区分中的优势。

3)基于优选的CP-SAR特征参数,建立了不同种植方式、不同品种的3类水稻田的精细分类方法,TH与DJ的分类精度较好,平均精度分别达到88%和82%。TJ的分类结果相对较差,平均精度达到60%。

4)当利用3个时相CP-SAR数据水稻精细分类时,基于SVM + SFS优选特征的分类结果优于全部特征的分类结果。

但是TJ分类精度不高,应继续分析TJ与TH,DJ的差异,充分利用CP-SAR数据,提高TJ分类精度将是我们下一步的工作重点。

参考文献

郭华东, 李新武 .

新一代SAR对地观测技术特点与应用拓展

[J]. 科学通报, 2011,56(15):1155-1168.

[本文引用: 1]

Guo H D, Li X W .

Technical characteristics and potential application of the new generation SAR for earth observation

[J]. Chinese Science Bulletin, 2011,56(15):1155-1168.

[本文引用: 1]

Souyris J C, Imbo P, Fjortoft R , et al.

Compact polarimetry based on symmetry properties of geophysical media:The /spl pi//4 mode

[J]. IEEE Transactions on Geoscience and Remote Sensing, 2005,43(3):634-646.

DOI:10.1109/TGRS.2004.842486      URL     [本文引用: 1]

We assess the performance of synthetic aperture radar (SAR) compact polarimetry architectures based on mixed basis measurements, where the transmitter polarization is either circular or orientated at 45/spl deg/(/spl pi//4), and the receivers are at horizontal and vertical polarizations with respect to the radar line of sight. An original algorithm is proposed to reconstruct the full polarimetric (FP) information from this architecture. The performance assessment is twofold: it first concerns the level of information preserved in comparison with FP, both for point target analysis and crop fields classification, using L-band SIRC/XSAR images acquired over Landes forest and Jet Propulsion Laboratory AIRSAR images acquired over Flevoland. Then, it addresses the space implementation complexity, in terms of processed swath, downloading features, power budget, calibration, and ionospheric effects. The polarization uniqueness in transmission of this mixed basis mode, hereafter referred to as the /spl pi//4 mode, maintains the standard lower pulse repetition frequency operation and hence maximizes the coverage of the sensor. Because of the mismatch between transmitter and receiver basis, the power budget is deteriorated by a factor of 3 dB, but it can partly be compensated.

Raney R K .

Dual-polarized SAR and stokes parameter

[J]. IEEE Geoscience and Remote Sensing Letters, 2006,3(3):317-319.

DOI:10.1109/LGRS.2006.871746      URL     [本文引用: 1]

For a given transmission polarization, the four-element Stokes vector captures all of the information inherent to the dual-polarized backscattered signals. Stokes parameters are linear combinations of the like-polarized power, the cross-polarized power, and the cross product between the complex image amplitudes in the two receive channels. Most so-called dual-polarized synthetic aperture radar (SAR) systems provide only the like- and cross-product images, ignoring the cross product. In so doing, potentially valuable information is thrown away. For most dual-polarized SARs, data sufficient to generate the Stokes parameters can be generated with small marginal cost. A practical and sufficient output data format would be a mapping of the (complex) cross product between the two channels, as well as their respective images

Raney R K .

Hybrid-polarity SAR architecture

[J]. IEEE Transactions on Geoscience and Remote Sensing, 2007,45(11):3397-3404.

DOI:10.1109/TGRS.2007.895883      URL     [本文引用: 1]

Collins M J, Denbina M, Atteia G .

On the reconstruction of quad-pol SAR data from compact polarimetry data for ocean target detection

[J].IEEE Transactions on Geoscience and Remote Sensing 2013,(51) 1:591-600.

DOI:10.1109/TGRS.2012.2199760      URL     [本文引用: 1]

Circular-transmit/linear-receive compact polarimetry synthetic aperture radar systems combine coherent dual polarization with wide-swath imaging. The polarization information in these data may be represented as a Stokes vector, or one can reconstruct several quadpolarized covariance elements. Two reconstruction algorithms have been published in the literature: one by Souyris et al. and a refined algorithm by Nord et al. We investigated the application of these two algorithms for reconstructing ocean clutter for the purpose of detecting targets. We tested the assumptions underlying these algorithms and found that they were not valid for ocean scenes. We present a simple empirical reconstruction model whose reconstruction and target detection performance is superior to the two published algorithms.

张红, 谢镭, 王超 , .

简缩极化SAR数据信息提取与应用

[J]. 中国图象图形学报, 2013,18(9):1065-1073.

[本文引用: 1]

Zhang H, Xie L, Wang C , et al.

Information extraction application of compact polarimetric SAR data

[J]. Journal of Image and Graphics, 2013,18(9):1065-1073.

[本文引用: 1]

Ainsworth T L, Kelly J P, Lee J S .

Classification comparisons between dual-pol,compact polarimetric and quad-pol SAR imagery

[J]. ISPRS Journal of Photogrammetry and Remote Sensing, 2009,64(5):464-471.

DOI:10.1016/j.isprsjprs.2008.12.008      URL     [本文引用: 1]

We present a study of the polarimetric information content of dual-pol imaging modes and dual-pol imaging extended by polarimetric scattering models. We compare Wishart classifications both among the partial polarimetric datasets and against the full quad-pol dataset. Our emphasis is the inter-comparisons between the classification results based on dual-pol modes, compact polarimetric modes and scattering model extensions of the compact polarimetric modes. We primarily consider novel dual-pol modes, e.g. transmitting a circular polarization and receiving horizontal and vertical polarizations, and the pseudo-quad-pol data derived from polarimetric scattering models based on dual-pol data. We show that the overall classification accuracy of the pseudo-quad-pol data is essential the same as the classification accuracy obtained directly employing the underlying dual-pol imagery.

Garestier F ,

Dubois-Fernandez P C,Papathanassiou K.Pine forest height inversion using single-pass X-band PolInSAR data

[J]. IEEE Transactions on Geoscience and Remote Sensing, 2008,46(1):59-68.

DOI:10.1109/TGRS.2007.907602      URL     [本文引用: 1]

A sparse pine forest is investigated at X-band on a single-pass polarimetric synthetic aperture radar interferometry (PolInSAR) data set using HH and HV channels. These first preliminary results show that the associated phase centers present a significant vertical separation (about 6 m) allowed by penetration through gaps in the canopy. Forest parameter inversion using the random volume over ground (RVoG) model is evaluated and adapted at this frequency. The forest height can be retrieved accurately by supposing a high mean extinction coefficient (around 1.6 (dB/m). The penetration depth is estimated to be around 4 m, based on the forest height ground measurements. Finally, a time-frequency analysis using a sublook decomposition is performed to increase the vertical separation of the polarimetric phase centers. As a consequence, RVoG-inversion performance is improved, and a penetration depth that is in better accordance with a previous work (of the order of 2 m) is found. This paper has shown that the height inversion of a pine forest was possible using PolInSAR X-band data and that the performance was more dependent on the forest density than at lower frequencies.

Zhang X, Dierking W, Zhang J , et al.

Retrieval of the thickness of undeformed sea ice from simulated C-band compact polarimetric SAR images

[J]. The Cryosphere, 2016,10(5):1529-1545.

DOI:10.5194/tc-10-1529-2016      URL     [本文引用: 1]

In this paper we introduce a parameter for the retrieval of the thickness of undeformed first-year sea ice that is specifically adapted to compact polarimetric??(CP) synthetic aperture radar??(SAR) images. The parameter is denoted as the "CP ratio". In model simulations we investigated the sensitivity of the CP ratio to the dielectric constant, ice thickness, ice surface roughness, and radar incidence angle. From the results of the simulations we deduced optimal sea ice conditions and radar incidence angles for the ice thickness retrieval. C-band SAR data acquired over the Labrador Sea in circular transmit and linear receive??(CTLR) mode were generated from RADARSAT-2 quad-polarization images. In comparison with results from helicopter-borne measurements, we tested different empirical equations for the retrieval of ice thickness. An exponential fit between the CP ratio and ice thickness provides the most reliable results. Based on a validation using other compact polarimetric SAR images from the same region, we found a root mean square (rms) error of 8???cm and a maximum correlation coefficient of??0.94 for the retrieval procedure when applying it to level ice between 0.1??and 0.8???m thick.

Truongloi M L, Freeman A, Dubois-Fernandez P C , et al.

Estimation of soil moisture and Faraday rotation from bare surfaces using compact polarimetry

[J]. IEEE Transactions on Geoscience and Remote Sensing, 2009,47(11):3608-3615.

DOI:10.1109/TGRS.2009.2031428      URL     [本文引用: 1]

The potential of compact polarimetry (CP) mode at longer wavelengths in a space environment for surface parameter estimation is investigated. CP consists of transmitting a single polarization while receiving two polarizations. At longer wavelengths, one of the main challenges associated with CP from space is Faraday rotation (FR) estimation and correction. In this paper, an estimation procedure...

Choudhury I, Chakraborty M .

SAR signature investigation of rice crop using RADARSAT data

[J]. International Journal of Remote Sensing, 2006,27(3):519-534.

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

The current study investigates the potential of multi‐temporal RADARSAT ScanSAR Narrow Beam B (SCNB) data to monitor rice crop growth and condition where special emphasis was given to the signature analysis of the crop. The study area is located in the Baleshwar and Bhadrak districts of Orissa. The temporal variations of radar backscatter of all land‐cover classes were analysed as a function of time. The analysis of the Synthetic Aperture Radar (SAR) backscatter coefficient (σ0) of rice crop shows significant temporal behaviour and a large dynamic range during its growth period, which is due to the interaction of microwave radiation with the crop canopy, increasing from the transplanting stage to the reproductive stage. This temporal variation of SAR backscatter clearly differentiates rice fields from other land‐cover classes. Separability studies among different class pairs carried out using t‐test and Bhattacharya distance show that all the rice classes are separable from each other except the early rice, which was mixed with far sides of hills and shadow. Knowledge‐based decision rule classifier based on the temporal evaluation of SAR backscatter was attempted to classify rice and non‐rice areas, and achieved more than 98% accuracy in the case of rice class. The results are promising and confirm the possibility of operational use of RADARSAT data for rice crop growth monitoring.

Bouvet A, Toan T L .

Use of ENVISAT/ASAR wide swath data for timely rice fields mapping in the Mekong River Delta

[J]. Remote Sensing of Environment, 2011,115(4):1090-1101.

DOI:10.1016/j.rse.2010.12.014      URL     [本文引用: 1]

78We develop a large-scale rice fields mapping method using wide-swath SAR data. 78Multi-track datasets are used to increase the observation frequency. 78The temporal change of backscatter in each track is an effective classifier. 78Precise planted area estimations are hampered by the coarse resolution. 78Operational application is possible in integrated systems using ancillary data.

Bouvet A, Toan T L, Lam-Dao N .

Monitoring of the rice cropping system in the Mekong Delta using ENVISAT/ASAR dual polarization data

[J]. IEEE Transactions on Geoscience and Remote Sensing, 2009,47(2):517-526.

DOI:10.1109/TGRS.2008.2007963      URL     [本文引用: 1]

The rice cropping system in Asia is undergoing major changes to cope with increasing demography and changing climate, making rice monitoring a critical issue. Past studies have demonstrated the use of C-band synthetic aperture radar (SAR) data to map rice areas. The methods were based on the temporal change of intensity backscattering coefficient of vertically or horizontally co-polarized data (VV or HH). In this paper, we assess the use of the HH/VV polarization ratio derived from Advanced SAR (ASAR) data from ENVISAT data for the production of rice paddy maps. The approach is based on past knowledge on the polarization behavior of rice canopy, i.e., VV backscattering is much lower than HH during a large part of the rice season, due to the attenuation of the wave by the vertical structure of the plants. The methodology is developed for the Mekong Delta, Vietnam, where a complex cropping pattern is found (one to three crops of rice per year). The approach includes a statistical analysis of the HH/VV distributions of rice and non-rice classes at different dates. The analysis results confirm that HH/VV can be used as classifier and point out the need for relevant speckle filtering prior to classification. A classification method is developed and applied to single- and multidate data sets. The methods are tested at one district of the province of An Giang and extended to the whole province. Comparisons of the mapping results to geographic-information-system land-use data and official agricultural statistics show very good agreement. The method will be further applied to the entire Mekong Delta.

Li K, Zhang F, Shao Y , et al.

Polarization signature analysis of paddy rice in southern China

[J]. Canadian Journal of Remote Sensing, 2011,37(1):122-135.

DOI:10.5589/m11-018      URL     [本文引用: 1]

Rice is one of the three largest food grains in the world. Radar remote sensing has been proven to be an effective tool for rice monitoring. With the emergence of spaceborne polarimetric synthetic aperture radar (SAR) satellites, research on polarization backscatter behaviors and identification methods for rice are of great significance and attract the attention of remote sensing communities. The Zhazuo area located in Guizhou province of southwest China was selected as the test site for this study. The RADARSAT-2 polarimetric data acquired on four different dates were used to analyze the polarization backscatter behaviors and temporal variation of rice. Identification methods for rice based on polarization combinations and polarization decompositions were developed. The results indicated that HH with HH/HV was the optimal polarization combination for rice identification, with an accuracy of up to 86.65%. Based on the Pauli decomposition, rice can be discriminated effectively, with an accuracy of about 87.00%. Furthermore, using combinations of different polarization decompositions the identification results were greatly improved. The combination of the Pauli decomposition and the parameter 脦虏 derived from eigenvector-eigenvalue-based decomposition was best for rice identification, with an accuracy of up to 93.51%.

Li K, Brisco B, Yun S , et al.

Polarimetric decomposition with RADARSAT-2 for rice mapping and monitoring

[J]. Canadian Journal of Remote Sensing, 2012,38(2):169-179.

DOI:10.5589/m12-024      URL     [本文引用: 1]

Synthetic aperture radar for paddy rice crop mapping is operational in many countries in southeast Asia due to the unique temporal backscatter signature when the rice grows above the water surface. Although the rapid increase of up to 10 dB or more allows for the identification and subsequent mapping of the rice crop extent the saturation of the backscatter signal during the later growth stages makes yield determination difficult. The development of polarimetric decomposition techniques may improve the capability to estimate rice production as the phase centre of the backscatter is expected to change as the crop matures and the water distribution in the plant parts changes. These changes can potentially be related to crop yield using growth models. This capability for monitoring growth stage and enhancing yield estimation is explored with a multitemporal RADARSAT-2 polarimetric data set of a rice growing region in China. The Freeman urden, the Cloude ottier, and the Touzi decompositions were used for classification and then evaluated as a function of the growth stage, which was monitored in the rice paddies during the image acquisition period. The Touzi decomposition produced the best classification accuracy when six or more parameters were used in the classification, but all polarimetric decompositions produced suitable accuracies for rice mapping. This indicated only one date of imagery is needed for rice identification when polarimetric data are available. The results also demonstrated some sensitivity to the growth stage in the polarimetric decompositions, especially with the Freeman urden and the Touzi parameters. A controlled experiment is planned to further develop the methodology for using polarimetric decompositions for rice monitoring and yield estimation.

Brisco B, Li K, Tedford B , et al.

Compact polarimetry assessment for rice and wetland mapping

[J]. International Journal of Remote Sensing, 2013,34(6):1949-1964.

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

Polarimetric RADARSAT-2 data of rice and wetlands are used to simulate compact polarimetry (CP) mode data from the upcoming RADARSAT Constellation Mission (RCM). The simulated CP data are then used to evaluate the information content for rice and wetland mapping using supervised classification, and the results are compared for linear and circular polarization combinations and polarimetric decompositions from the fully polarimetric data and the simulated CP data. The results are consistent for both rice and wetlands and show that the classification accuracy increases as one goes up the polarization hierarchy. The circular polarizations produced the best classification results for the polarization combinations. This result requires further research to verify. Although the CP data did not perform as well as the fully polarimetric data, the results were better than for dual polarization, and this mode may offer the best option for rice and wetland mapping applications because of swath coverage. Note that both the compact simulations and the fully polarimetric data produced operationally suitable classification accuracies. Additional research is underway to evaluate the monitoring capability of this new CP mode. This article describes the approach used for the analyses and the classification results for both rice and wetlands.

Uppala D, Kothapalli R V, Poloju S , et al.

Rice crop discrimination using single date RISAT1 Hybrid (RH,RV) polarimetric data

[J]. Photogrammetric Engineering and Remote Sensing, 2015,81(7):557-563.

DOI:10.14358/PERS.81.7.557      URL     [本文引用: 1]

Rice is the most important food grain crop in India and contributes to more than 40 percent of the country’s food grain production. Spaceborne remote sensing offers economically viable and accurate production and area statistics. The utility of optical remote sensing in mapping rice cropped area is limited by persistent cloud cover during monsoon season. Temporal availability of SAR data has facilitated an operational procedure to monitor the rice crop. The current study discriminates rice crop, using single date hybrid polarimetric data available from RISAT-1 SAR. This was subjected to Raney m-δ, m-χ decompositions, and supervised classification was performed. The accuracy was estimated using the field points. The results were compared with rice map generated using optical sensor Resourcesat-2 LISS-IV and statistical data. The spatial agreement between the estimate from RISAT-1 and LISS-IV data was found to be 85 percent. The class kappa value was 0.94 and 0.92 for LISS-IV and RISAT-1, respectively.

吴永辉, 计科峰, 郁文贤 .

SVM全极化SAR图像分类中的特征选择

[J]. 信号处理, 2007,23(6):877-881.

DOI:10.3969/j.issn.1003-0530.2007.06.019      URL     [本文引用: 1]

本文将SVM用于全极化SAR图像分类,并提出一种新的应用于 SVM分类的特征选择算法.该算法以支持向量个数作为特征评估准则,利用顺序前进法加入特征.基于NASA/JPL实验室AIRSAR系统的L波段荷兰 Flevoland全极化数据的与RELIEF-F算法的对比实验表明,在特征个数更少(或相当)的情况下,本文特征选择算法能在更广泛的SVM参数取值 范围内获得更高的分类精度.

Wu Y H, Ji K F, Yu W X .

A new feature selection algorithm for SVM-based fully polarimetric SAR image classification

[J]. Signal Processing, 2007,23(6):877-881.

[本文引用: 1]

Raney R K, Cahill J T S, Patterson G W , et al.

The m-chi decomposition of hybrid dual-polarimetric Radar data with application to lunar craters

[J]. Journal of Geophysical Research Atmospheres, 2012,117(E12):5093-5096.

DOI:10.1029/2011JE003986      URL     [本文引用: 3]

[1] We introduce a new technique derived from the classical Stokes parameters for analysis of polarimetric radar astronomical data. This decomposition is based on m (the degree of polarization) and chi (the Poincar ellipticity parameter). Analysis of the crater Byrgius A demonstrates how m-chi can more easily differentiate materials within ejecta deposits and their relative thicknesses. We use Goldschmidt crater to demonstrate how m-chi can differentiate coherent deposits of water ice. Goldschmidt crater floor is found to be consistent with single bounce Bragg scattering suggesting the absence of water ice and further corroborating adsorbed H to mineral grains or an H2O frost as plausible explanations for a H2O/OH detection by near-infrared instruments.

Charbonneau F J, Brisco B, Raney R K , et al.

Compact polarimetry overview and applications assessment

[J]. Canadian Journal of Remote Sensing, 2010,36(s2):298-315.

DOI:10.5589/m10-062      URL     [本文引用: 1]

A synthetic aperture radar (SAR) with hybrid-polarity (CL-pol) architecture transmits circular polarization and receives two orthogonal, mutually coherent linear polarizations, which is one manifestation of compact polarimetry. The resulting radar is relatively simple to implement and has unique self-calibration features and low susceptibility to noise. It also enables maintenance of a larger swath coverage than fully polarimetric SAR systems. A research team composed of various departments of the Government of Canada evaluated this compact polarimetry mode configuration for application to soil moisture estimation, crop identification, ship detection, and sea-ice classification. This paper presents an overview of compact polarimetry, the approach developed for evaluation, and preliminary results for applications important to the Government of Canada. The implications of the results are also discussed with respect to future SAR missions such as the Canadian RADARSAT Constellation Mission, the American DESDynI, and India鈥檚 RISAT.

/

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