文章信息
- 李鹏, 黎达辉, 李振洪, 王厚杰
- LI Peng, LI Dahui, LI Zhenhong, WANG Houjie
- 黄河三角洲地区GF-3雷达数据与Sentinel-2多光谱数据湿地协同分类研究
- Wetland Classification Through Integration of GF-3 SAR and Sentinel-2B Multispectral Data over the Yellow River Delta
- 武汉大学学报·信息科学版, 2019, 44(11): 1641-1649
- Geomatics and Information Science of Wuhan University, 2019, 44(11): 1641-1649
- http://dx.doi.org/10.13203/j.whugis20180258
-
文章历史
收稿日期: 2018-11-04
2. 青岛海洋科学与技术试点国家实验室海洋地质过程与环境功能实验室, 山东 青岛, 266061;
3. 华东师范大学地理科学学院, 上海, 200241;
4. 英国纽卡斯尔大学工程学院, 英国 纽卡斯尔, NE1 7RU
2. Laboratory of Marine Geology, Pilot National Laboratory for Marine Science and Technology(Qingdao), Qingdao 266061, China;
3. School of Geographic Sciences, East China Normal University, Shanghai 200241, China;
4. School of Engineering, Newcastle University, Newcastle Upon Tyne NE1 7RU, UK
湿地系统是全球自然生态环境的重要组成部分,与海洋、森林等并称为全球三大生态系统,其在调节气候、补充地下水、蓄洪防旱、降解污染、土壤侵蚀控制等方面发挥着重要作用[1]。黄河三角洲地处陆海相互作用的敏感地带,其河口堆积来自陆域径流的大量泥沙,形成了许多新生湿地,同时又受到海水波流动力侵蚀影响,湿地资源变化较为频繁[2]。
目前,黄河三角洲湿地面临断流、污染、海水入侵等威胁,生态破损严重,湿地面积严重萎缩,生态功能逐渐退化[3]。因此,对黄河口三角洲湿地进行分类、识别,了解黄河三角洲湿地的动态变化及其影响因素,对该地区湿地资源的合理利用、开发保护具有极为重要的意义。
合成孔径雷达(synthetic aperture radar, SAR)数据具有全天候、全天时、高空间分辨率等特点,重访观测周期通常可达一两周至一两个月左右,空间上可以连续获取,其侧视成像、多极化能力可以反映地物纹理特征、地表植被散射机制,上述特点有利于开展湿地分类、识别与监测。然而,SAR作为一种相干测量系统,其信号受到相干斑噪声的干扰,进而会影响其地物识别精度。
由于光学数据能反映地物的光谱特征信息,在一定程度上可以弥补极化SAR数据在地物识别方面的缺陷。与单一数据源相比,将光学与雷达数据结合用于地物分类,能减少或消除目标地物可能存在的多义性、不完全性、不确定性和差异性[4]。
已有研究表明,结合光学与微波雷达影像各自的优点,可增加数据所含的有用信息,以增强对地物的识别能力。例如,文献[5]融合RADARSAT-1 SAR和Landsat专题制图仪(thematic mapper, TM)数据,与原始TM影像分类结果相比,分类精度有较大的提高;文献[6]结合环境卫星多光谱影像与Envisat ASAR单极化影像,协同分类精度比单独利用多光谱数据提高约5%;文献[7]利用高分三号(GF-3)极化SAR和Sentinel-2A光学数据对洞庭湖盆地一年一熟、一年两熟水稻等作物进行分类,其总体精度达到了85.3%,且高于两者单独分类的总体精度。
本文利用黄河三角洲地区Sentinel-2B光学和GF-3雷达影像,结合该湿地地区7类典型地物,在辐射校正、配准等预处理基础上,通过主成分变换、波段指数计算、极化分解等方式提取不同类型的特征协同集成,开展最大似然法(maximum likelihood, ML)、基于分类回归树(classification and regression tree, CART)算法的决策树(decision tree, DT)、支持向量机(support vector machine, SVM)方法的有监督分类,并对分类结果进行分析与评价。
1 研究区域与数据 1.1 研究区域黄河三角洲位于山东省东营市黄河入海口,东邻莱州湾,北靠渤海湾。黄河三角洲湿地主要指1976年5月黄河尾闾改道至清水沟流路入海至今形成的新生湿地,介于118°42′E~119°20′E、37°34′N~37°59′N之间,面积达1 144.8 km2。黄河三角洲位于北半球暖温带欧亚大陆东部,受海陆热力性质差异影响形成温带季风气候,四季分明,夏季高温多雨,冬季寒冷干燥[8]。
本文研究区域范围如图 1中红色框所示,图 1中的底图为Landsat-8陆地成像仪(operational land imager, OLI)5/6/4波段合成影像。
1.2 研究数据目前,黄河三角洲地区的湿地分类研究主要基于光学遥感卫星和国外SAR卫星(如RADARSAT-2、ALOS-2等)[9-10]。
Sentinel-2B卫星由欧洲空间局于2017年3月发射,与2015年6月发射的Sentinel-2A组成星座,具有多光谱成像(13个波段)、宽幅(290 km)、高分辨率(10~60 m)等特性,侧重于监测全球陆地和沿海地区的地物覆盖变化。
国产GF-3 SAR卫星具备1~500 m空间分辨率、10~650 km成像幅宽、12种成像模式的观测能力,是我国实施海洋开发、陆地环境资源监测和防灾减灾的重要支撑[11-13]。
如表 1、表 2所示,本文选取2017-10获取的Sentinel-2B多光谱影像与GF-3全极化条带模式SAR影像作为实验数据。
卫星名称 | 数据等级 | 波段 | 成像方式 | 极化方式 | 空间分辨率(方位向)/m | 获取时间 | 入射角/(°) |
GF-3 | SLC | C | 全极化条带I (QPSI) | HH/HV/VH/VV | 8 | 2017-10-12 | 22.27~25.00 |
Sentinel-2B L1C级数据是经过几何精校正的正射影像,但是没有进行辐射定标和大气校正,可以通过欧洲空间局专门用于实现Sentinel-2 L1C级到L2A级产品数据转换的Sen2Cor插件获得经过辐射定标和大气校正的大气底层反射率数据。考虑到GF-3分辨率为8 m,本文仅提取Sentinel-2B分辨率为10 m的4个波段(波段2/3/4/8,分别对应蓝光、绿光、红光、近红外波段)用于分类特征的提取。
本文采用的Sentinel-2B预处理及地物分类详细流程如图 2所示,图 2中的PCA1、PCA2为主成分变换(principal component analysis, PCA)后的首要两个主成分分量,INDVI为归一化植被指数,IBI为土壤亮度指数,INDWI2为第二类归一化水体指数。
2.2 GF-3预处理为了获得较好的分类特征,同样需要对GF-3全极化(HH、HV、VH、VV)单视复数影像(single look complex, SLC)数据进行预处理。本文首先进行辐射定标获得后向散射系数,然后使用改进型Lee滤波进行极化滤波处理,以减弱相干斑等噪声的影响,以便获得准确的强度、纹理、极化散射等特征信息。本文采用的GF-3预处理及地物分类详细流程如图 3所示,为了将Sentinel-2B预处理结果与GF-3预处理结果进行对比,将两者预处理结果放在了一起,结果如图 4所示。
3 特征提取与分类 3.1 分类算法本文采用目前常用的3种经典监督分类算法进行分类[14]。
1)基于统计模型的最大似然法。该方法建立在贝叶斯判决规则上,通过似然函数判断类别。由于遥感信息统计分布具有高度的复杂性和随机性,使用该方法常会产生较大误差。
2)基于知识的决策树分类。该方法通过对训练样本的特征属性和目标变量建立适当的树状分类结构进行分类,具有灵活、直观、效率高等特点。常用的决策树算法有ID3、C4.5、CART等,本文使用CART决策树分类算法。
3)基于统计学习理论的支持向量机法[15]。该方法利用非线性变换将输入空间的线性不可分问题转化到高维空间,根据间隔最大化准则构造最优分类超平面[16]。由于在处理非线性数据和高维特征空间识别方面具有独特优势,该方法已广泛用于遥感影像分类和多源遥感影像信息融合等领域。
3.2 Sentinel-2B特征提取与分类如图 2所示,Sentinel-2B数据经过预处理后的4个波段之间信息并不独立。为了消除波段间的相关性,减少信息冗余,通过PCA提取两个独立的主成分分量PCA1、PCA2(贡献率大于99%)。结合研究区植被、裸地、水体多的特点,利用提取得到的4个波段计算3种指数特征,分别是INDVI、IBI和INDWI2。其中,INDVI可以反映植被的密度和健康程度;IBI对土壤的亮度敏感,与土壤的湿度及盐度相关,可在一定程度上反映裸土的分布;INDWI2可以检测陆地表面水的范围,即使在湿地环境下,仍有一定的识别水体能力。
3种指数特征的计算方式如下:
${{I}_{\text{NDVI}}}=\frac{{{\rho }_{\text{NIR}}}-{{\rho }_{\text{RED}}}}{{{\rho }_{\text{NIR}}}+{{\rho }_{\text{RED}}}}=\frac{{{B}_{8}}-{{B}_{4}}}{{{B}_{8}}+{{B}_{4}}}$ | (1) |
${{I}_{\text{BI}}}=\sqrt{\frac{{{\rho }_{\text{GREEN}}}+{{\rho }_{\text{RED}}}}{2}}=\sqrt{\frac{{{B}_{3}}+{{B}_{4}}}{2}}$ | (2) |
$ {{I}_{\text{NDWI}2}}=\frac{{{\rho }_{\text{GREEN}}}-{{\rho }_{\text{NIR}}}}{{{\rho }_{\text{GREEN}}}+{{\rho }_{\text{NIR}}}}=\frac{{{B}_{3}}-{{B}_{8}}}{{{B}_{3}}+{{B}_{8}}} $ | (3) |
式中,ρNIR、ρRED、ρGREEN表示近红外、红、绿波段的反射率;B3、B4、B8表示3、4、8波段的波长
结合文献[1-3, 17]及实地调查资料,黄河三角洲湿地的典型地物主要包括:
1)河流;
2)水塘,包括人工养殖塘、自然池塘、水库等;
3)居民地;
4)草地;
5)裸地,包括滩涂、盐碱地等空地;
6)林地;
7)庄稼地。
对应的训练样本分布如图 5所示。
根据图 2将Sentienl-2B得到的5个分类特征(PCA1、PCA2、NDVI、BI、NDWI2)和样本数据输入到ML、DT、SVM等3种监督分类器,分类结果如图 6(a)、6(b)、6(c)所示。
3.3 GF-3特征提取与分类GF-3极化散射特征主要利用Freeman三分量分解法和H-A-α分解得到。
Freeman三分量分解法是以物理实际为基础,对极化协方差矩阵C3建立3种基本散射机制模型,分别为适度粗糙表面散射、随机方向偶极子组成的体散射、由正交平面构成的二次散射[18-19]。假设表面散射PS、二次散射PD和体散射分量PV互不相关,由每种散射机制成分的权重值(fS、fD、fV)可以得到目标极化总功率PSPAN,即上述每种独立散射机制成分统计量的叠加之和[20]:
$ \begin{array}{*{35}{l}} {{P}_{\text{SPAN}}}={{\left| {{S}_{\text{HH}}} \right|}^{2}}+2{{\left| {{S}_{\text{HV}}} \right|}^{2}}+{{\left| {{S}_{\text{VV}}} \right|}^{2}}= \\ \ \ \ \ \ \ \ {{P}_{S}}+{{P}_{D}}+{{P}_{V}} \\ \end{array} $ | (4) |
其中,
总功率PSPAN是由HH、HV、VH、VV 4通道强度叠加得到,强度特征较为强烈,叠加后部分正负噪声干扰得到了消除或减弱。考虑到不同地物纹理特征存在差异,结合纹理特征可以提升地物分类精度,因此本文利用灰度共生纹理矩阵作用于极化总功率交换器端口分析器(switched port analyzer,SPAN)图像来提取纹理特征,共得到8个纹理特征[21]。由于这些纹理特征并非相互独立,为了减少数据冗余量,通过7类地物感兴趣区域(region of interest, ROI)样本均值来提取纹理特征。另外,纹理特征单位或数量级存在差异,需要将其归一化才能区分差异。
图 7给出了GF-3 ROI样本纹理特征均值曲线。由ROI样本区域均值的统计分析可知,均值、均匀性、非相似性3个特征对7类样本的区分较为明显,故选择这3个纹理特征参与分类。
H-A-α分解利用散射矩阵S转换得到相干矩阵〈[T3]〉,其中〈[T3] 〉是半正定埃尔米特矩阵。根据矩阵分解理论,可将其扩展为3个独立目标之和,即用3个散射矩阵来表示3个独立目标[22],即:
$ \left\langle \left[ {{T}_{3}} \right] \right\rangle =\sum\limits_{j=1}^{3}{{{\lambda }_{j}}}{{\mathrm{u}}_{j}}\mathrm{u}_{j}^{\text{*T}}=\sum\limits_{j=1}^{3}{{{\lambda }_{j}}}{{T}_{j}} $ | (5) |
式中,〈[·]〉表示多视处理;上标*表示共轭;上标T表示转置;Tj、λi和uj分别表示第j个独立目标、特征值、特征向量。
H-A-α分解的3个二阶参数是〈[T3]〉的特征值和特征矢量函数,其定义如下。
1)极化散射熵H:
$ H=-\sum\limits_{k=1}^{3} P_{k} \log _{3}\left(P_{k}\right) $ | (6) |
其中,
2)平均α角:
$ \alpha ={{P}_{1}}{{\alpha }_{1}}+{{P}_{2}}{{\alpha }_{2}}+{{P}_{3}}{{\alpha }_{3}}$ | (7) |
式中,α1、α2、α3的大小指示主要的散射机制
3)极化各向异性度A:
$ A=\frac{{{\lambda }_{2}}-{{\lambda }_{3}}}{{{\lambda }_{2}}+{{\lambda }_{3}}}$ | (8) |
根据图 3所示流程图,本文将GF-3 SAR影像得到的强度特征SPAN、3个Freeman三分量分解极化散射特征、3个H-A-α分解极化散射特征、3个纹理特征作为敏感特征与训练样本输入ML、DT、SVM等3种监督分类器,分类结果见图 6。
3.4 Sentinel-2B和GF-3协同分类Sentinel-2B和GF-3的精确配准是实现两者协同分类的关键,本文通过野外实地调查及目视选择控制点实现Sentinel-2B和GF-3数据准确配准。将配准后的Sentinel-2B和GF-3影像特征(共15个)进行叠加,与训练样本一起输入到ML、DT、SVM等3种监督分类器,实现两者的协同分类,分类结果如图 6(g)、6(h)、6(i)所示。
3.5 结果分析与讨论本文采用生产者精度、用户精度、总体精度、Kappa系数作为定量评估不同数据与分类算法的性能指标[21, 23],结合图 6和表 3可知:
类别 | 最大似然法(ML) | 决策树(DT) | 支持向量机(SVM) | |||||||||
Sentinel-2B | GF-3 | 两者协同 | Sentinel-2B | GF-3 | 两者协同 | Sentinel-2B | GF-3 | 两者协同 | ||||
生产者精度/% | 河流 | 88.22 | 50.41 | 91.19 | 91.46 | 85.55 | 95.95 | 88.35 | 87.82 | 95.75 | ||
水塘 | 87.20 | 56.01 | 91.44 | 88.17 | 65.12 | 95.65 | 88.11 | 57.12 | 95.48 | |||
居民地 | 80.31 | 81.98 | 93.09 | 36.62 | 48.88 | 68.14 | 34.83 | 57.07 | 79.26 | |||
草地 | 55.72 | 84.89 | 79.53 | 84.77 | 94.43 | 95.74 | 83.84 | 94.47 | 94.05 | |||
裸地 | 95.04 | 92.95 | 96.13 | 98.97 | 92.84 | 99.08 | 99.25 | 92.12 | 98.92 | |||
林地 | 75.95 | 87.95 | 89.98 | 84.86 | 92.60 | 94.39 | 82.40 | 94.14 | 95.92 | |||
庄稼地 | 92.21 | 28.84 | 94.00 | 41.69 | 8.30 | 49.99 | 67.92 | 6.15 | 63.19 | |||
用户精度/% | 河流 | 98.22 | 66.11 | 98.17 | 97.67 | 75.20 | 98.57 | 98.5 | 72.00 | 97.62 | ||
水塘 | 88.12 | 69.82 | 92.46 | 91.99 | 79.46 | 97.69 | 90.67 | 83.02 | 94.86 | |||
居民地 | 23.02 | 26.35 | 31.82 | 83.93 | 82.80 | 90.29 | 70.43 | 77.53 | 87.22 | |||
草地 | 89.45 | 84.72 | 94.74 | 84.20 | 83.43 | 89.00 | 85.83 | 84.84 | 93.61 | |||
裸地 | 94.43 | 78.42 | 95.86 | 92.41 | 94.81 | 97.36 | 91.18 | 94.59 | 97.89 | |||
林地 | 84.79 | 95.22 | 94.03 | 81.73 | 94.7 | 93.52 | 83.62 | 94.20 | 94.29 | |||
庄稼地 | 16.32 | 19.77 | 32.64 | 68.27 | 56.15 | 84.14 | 45.02 | 48.26 | 78.26 | |||
总体精度/% | 81.91 | 76.77 | 90.38 | 89.45 | 86.48 | 95.44 | 88.71 | 86.07 | 95.65 | |||
Kappa系数 | 0.777 0 | 0.706 2 | 0.880 1 | 0.866 3 | 0.828 7 | 0.942 3 | 0.857 3 | 0.823 4 | 0.945 1 |
1)由3种数据的分类结果来看,极大似然法的精度最低,决策树和支持向量机的精度相对较高,且两者差别微弱,但仅凭精度难以区分后两个算法的好坏,因为取不同参数时,不同分类算法的分类精度也会有变化。
2)从生产者精度和用户精度综合来看,居民地和庄稼地的分类精度较低,这主要是因为它们在图中范围分布数量较小(影像像素少),且往往是局部分布,容易受到其邻近地物拓扑关联(邻接、包含等)分布的影响。同样地,尽管两者分类精度较低,但它们占比少,对于总体精度、Kappa系数影响小,只要其他地类的分类精度足够高,仍然可以得到较高的总体精度和Kappa系数。
3)其他5类地物的分类精度相对较高,特别是裸地在所有分类方案下的精度普遍达到90%以上,这可能与该地物在该地区面积较大且连片分布有关。水塘和河流、草地和林地之间存在部分混淆,主要由于其特征较为相似且存在邻接边界。
4)如图 6(d)、6(e)、6(f)所示,GF-3大部分海域被分为裸土,这是由于裸地训练样本包含了南部潮间带地区,该地区裸地定期受涨潮影响而被海水覆盖,而训练样本没有考虑海域部分,并且两者均有较强的表面散射特征,因此海域可能被错分为裸土区域。
5)总体上,GF-3的总体分类精度和Kappa系数略低于Sentinel-2B分类,而其居民地分类精度要优于Sentinel-2B,这可能是由于居民地在GF-3影像二次散射较强烈,有利于将其区分出来。
6)两者协同分类后各类地物的分类精度比单独分类得到一定程度的提升,尤其是Sentinel-2B和GF-3分类精度都较低的居民地和庄稼地。协同分类的总体精度和Kappa系数均达到较高数值,主要是由于光学数据与SAR数据的协同分类得到的是“1+1 > 2”的信息融合与互补,更加有利于地物识别。
4 结语本文利用Sentinel-2B多光谱数据和GF-3全极化SAR数据分别提取了黄河三角洲地区7类湿地地物的光谱、指数、极化散射、纹理等特征,分别利用最大似然法、决策树法、支持向量机3种算法实现了有监督协同分类,总体精度可达90%以上,均高于两者单独分类的总体精度,表明了基于两种不同种类的卫星遥感数据源进行湿地协同分类的可靠性和优势。
本文后续研究将考虑加入不同波段的多光谱数据和SAR数据(如L/X波段等)开展更细致的地物分类,以便准确反映该地区的湿地资源状况。结合黄河三角洲地区长时间序列多源遥感数据,有望实现该地区地物覆盖的动态分布监测。
致谢: 感谢国家卫星海洋应用中心提供高分三号SAR数据以及欧洲空间局提供Sentinel-2B多光谱数据。[1] |
Han Mei, Zhang Xiaohui, Liu Liyun. Research Progress on Wetland of the Yellow River Delta[J]. Ecology and Enviroment, 2006, 15(4): 872-875. (韩美, 张晓惠, 刘丽云. 黄河三角洲湿地研究进展[J]. 生态环境学报, 2006, 15(4): 872-875. DOI:10.3969/j.issn.1674-5906.2006.04.041 ) |
[2] |
Xiao Duning, Han Mukang, Li Xiaowen, et al. Sea Level Rising Around Bohai Sea and Deltaic Wetlands Protection[J]. Quaternary Sciences, 2003, 23(3): 237-246. (肖笃宁, 韩慕康, 李晓文, 等. 环渤海海平面上升与三角洲湿地保护[J]. 第四纪研究, 2003, 23(3): 237-246. DOI:10.3321/j.issn:1001-7410.2003.03.001 ) |
[3] |
Xiu Changjun, Wang Xiaohui. Environmental Impact and Management for the Yellow River Delta Swamp with Shengli Oilfield Development[J]. China environmental Management, 2003(3): 59-60. (修长军, 王晓慧. 胜利油田开发对黄河三角洲湿地的环境影响及环境管理[J]. 中国环境管理, 2003(3): 59-60. ) |
[4] |
Sun Jiabing. Principles and Applications of Remote Sensing[M]. 3rd ed. Wuhan: Wuhan University Press, 2009. (孙家抦. 遥感原理与应用[M]. 3版. 武汉: 武汉大学出版社, 2009. )
|
[5] |
Zhang Hailong, Jiang Jianjun, Wu Hong'an, et al. The BP Neural Network Classification Based on the Fusion of SAR and TM Images[J]. Acta Geodaetica et Cartographica Sinica, 2006, 35(3): 229-233. (张海龙, 蒋建军, 吴宏安, 等. SAR与TM影像融合及在BP神经网络分类中的应用[J]. 测绘学报, 2006, 35(3): 229-233. DOI:10.3321/j.issn:1001-1595.2006.03.007 ) |
[6] |
Jia Kun, Li Qiangzi, Tian Yichen, et al. Accuracy Improvement of Spectral Classification of Crop Using Microwave Backscatter Data[J]. Spectroscopy and Spectral Analysis, 2011, 31(2): 483-487. (贾坤, 李强子, 田亦陈, 等. 微波后向散射数据改进农作物光谱分类精度研究[J]. 光谱学与光谱分析, 2011, 31(2): 483-487. DOI:10.3964/j.issn.1000-0593(2011)02-0483-05 ) |
[7] |
Gao H, Wang C C, Wang G Y, et al. A Crop Classification Method Integrating GF-3 PolSAR and Sentinel-2A Optical Data in the Dongting Lake Basin[J]. Sensors, 2018, 18(9): 1-19. DOI:10.1109/JSEN.2018.2815438 |
[8] |
Zhang Lei, Gong Zhaoning, Wang Qiwei, et al. Wetland Mapping of Yellow River Delta Wetlands Based on Multi-feature Optimization of Sentinel-2 Images[J]. Journal of Remote Sensing, 2019, 23(2): 313-326. (张磊, 宫兆宁, 王启为, 等. Sentinel-2影像多特征优选的黄河三角洲湿地信息提取[J]. 遥感学报, 2019, 23(2): 313-326. ) |
[9] |
Buono A, Nunziata F, Migliaccio M, et al. Classification of the Yellow River Delta Area Using Fully Polarimetric SAR Measurements[J]. International Journal of Remote Sensing, 2017, 38(23): 6714-6734. DOI:10.1080/01431161.2017.1363437 |
[10] |
Gou S P, Li X F, Yang X F. Coastal Zone Classification with Fully Polarimetric SAR Imagery[J]. IEEE Geoscience and Remote Sensing Letters, 2016, 13(11): 1616-1620. DOI:10.1109/LGRS.2016.2597965 |
[11] |
Chang Y L, Li P X, Yang J, et al. Polarimetric Calibration and Quality Assessment of the GF-3 Satellite Images[J]. Sensors, 2018, 18(2): 1-12. |
[12] |
Zhang Qingjun. System Design and Key Technologies of the GF-3 Satellite[J]. Acta Geodaetica et Cartographica Sinica, 2017, 46(3): 269-277. (张庆君. 高分三号卫星总体设计与关键技术[J]. 测绘学报, 2017, 46(3): 269-277. ) |
[13] |
Yin Junjun, Yang Jian, Zhang Qingjun. Assessment of GF-3 Polarimetric SAR Data for Physical Scattering Mechanism Analysis and Terrain Classification[J]. Sensors, 2017, 17(12): 1-11. DOI:10.1109/JSEN.2017.2693641 |
[14] |
Guo Jian, Zhang Jixian, Zhang Yonghong, et al. Study of the Comparison of Land Cover Classification for Multitemporal MODIS Images[J]. Acta Geodaetica et Cartographica Sinica, 2009, 38(12): 88-92. (郭健, 张继贤, 张永红, 等. 多时相MODIS影像土地覆盖分类比较研究[J]. 测绘学报, 2009, 38(12): 88-92. ) |
[15] |
Vapnik V N. An Overview of Statistical Learning Theory[J]. IEEE Transactions on Neural Networks, 1999, 10(5): 988-999. DOI:10.1109/72.788640 |
[16] |
Wu Yonghui, Ji Kefeng, Yu Wenxian. Comparison of Classification Porfermance of Full-, Dual- and Single-Polarization SAR Images Using SVM[J]. Journal of Remote Sensing, 2008, 12(1): 46-53. (吴永辉, 计科峰, 郁文贤. 利用SVM的全极化、双极化与单极化SAR图像分类性能的比较[J]. 遥感学报, 2008, 12(1): 46-53. ) |
[17] |
Huang Jin. Classification of the Yellow River Estuary Wetland Based on Multiband and Multipolarization SAR Data[D]. Qingdao: China University of Petroleum(East China), 2011 (黄瑾.基于多波段多极化SAR数据的黄河口湿地分类研究[D].青岛: 中国石油大学(华东), 2011) http://d.wanfangdata.com.cn/Thesis/Y1875773
|
[18] |
Zhao Quanhua, Guo Shibo, Li Xiaoli, et al. Polarimetric SAR Sea Ice Classification Based on Target Decompositional Features[J]. Acta Geodaetica et Cartographica Sinica, 2018, 47(12): 1609-1620. (赵泉华, 郭世波, 李晓丽, 等. 利用目标分解特征的全极化SAR海冰分类[J]. 测绘学报, 2018, 47(12): 1609-1620. DOI:10.11947/j.AGCS.2018.20170551 ) |
[19] |
Freeman A, Durden S L. A Three-component Scattering Model for Polarimetric SAR Data[J]. IEEE Transactions on Geoscience Remote Sensing, 1998, 36(3): 963-973. DOI:10.1109/36.673687 |
[20] |
Liu Xiuguo, Huang Xiaodong, Chen Qihao, et al. Inconsistency Elimination in Freeman Decomposition for Polarimetric SAR Data on Normalization Method[J]. Geomatics and Information Science of Wuhan University, 2013, 38(3): 257-261. (刘修国, 黄晓东, 陈启浩, 等. 归一化法消除全极化SAR影像Freeman分解中的不一致[J]. 武汉大学学报·信息科学版, 2013, 38(3): 257-261. ) |
[21] |
Li Fenling, Chang Qingrui, Liu Jiaqi, et al. SVM Classification with Multi-texture Data of ZY-1 02C HR Image[J]. Geomatics and Information Science of Wuhan University, 2016, 41(3): 455-461, 486. (李粉玲, 常庆瑞, 刘佳岐, 等. 基于多纹理和支持向量机的ZY-1 02C星HR数据分类[J]. 武汉大学学报∙信息科学版, 2016, 41(4): 455-461, 486. ) |
[22] |
Hou Fei, Hu Zhaoling. Land Cover Classification Based on Fully Polarmetric SAR Data Supported by Topography Information[J]. Bulletin of Surveying and Mapping, 2016(3): 40-43. (侯飞, 胡召玲. 地形信息辅助下的全极化SAR数据土地覆盖分类[J]. 测绘通报, 2016(3): 40-43. ) |
[23] |
Wang Chunyan, Liu Jiaxin, Xu Aigong, et al. A New Method of Fuzzy Supervised Classification of High Resolution Remote Sensing Image[J]. Geomatics and Information Science of Wuhan University, 2018, 43(6): 922-929. (王春艳, 刘佳新, 徐爱功, 等. 一种新的高分辨率遥感影像模糊监督分类方法[J]. 武汉大学学报·信息科学版, 2018, 43(6): 922-929. ) |