2021遥感应用组三等奖:基于PolSAR和高光谱遥感的黄河口湿地协同分类研究
作品介绍
一、作品概述
受人类活动和气候变化的影响,滨海湿地正经历着显著的变化。滨海湿地生态系统的监测具有重要的现实意义。同时,随着卫星遥感技术的快速发展,可以利用高空间分辨率、高光谱分辨率和高时间分辨率的多源数据对滨海湿地进行高精度分类。
本作品利用珠海1号星座轨道高光谱卫星(OHS)和高分3号(GF-3)全极化合成孔径雷达(SAR)数据的光谱、指数、极化和纹理特征利用最大似然度(ML)、马氏距离(MD)和支持向量机(SVM)三种机器学习方法提取黄河三角洲7种滨海湿地。结果表明,两种数据协同分类的总体精度分别为97.3%、89.0%和97.2%,显著高于单一数据。
二、处理流程
本作品的主要技术流程如图1所示,主要包括5个模块,分别是:预处理、特征提取、数据融合、协同分类、分类后评估。
图1 本作品主要技术流程图
2.1 数据说明
本作品采用的数据是两种国产卫星遥感影像,GF-3全极化SAR影像和珠海1号OHS高光谱影像。其中GF-3和OHS影像的详细参数如下表1-3所示。
表1 GF-3全极化SAR影像参数
卫星名称 |
数据等级 |
波段 |
成像方式 |
极化方式 |
空间分辨率 (m) |
获取时间 |
GF-3 |
SLC |
C |
全极化条带Ⅰ |
HH+HV+VH+VV |
8 |
2019-03-13 |
表2 OHS高光谱影像参数
卫星名称 |
传感器 |
空间分辨率 (m) |
波段数 |
光谱范围 (nm) |
数据等级 |
获取时间 |
云量 |
OHS 2-D |
CMOSMSS |
10 |
32 |
400-1000 |
L1B |
2020-03-24 |
0% |
本研究的训练和验证样本主要来源于野外现场调查和谷歌地球、GF-2高分辨率光学遥感影像目视解译两种方法。其中GF-2影像的全色波段空间分辨率是0.8 m,多光谱波段空间分辨率为3.2 m,为了获取高分辨率的影像,本作品将GF-2全色和多光谱波段融合在一起。波段融合前,还需要将多波段和全色波段分别进行正射校正处理。正射校正后的影像借助ENVI 5.6中的“Gram-Schmidt Pan Sharpening”工具,多光谱数据作为低空间分辨率输入影像,全色数据作为高空间分辨率影像,传感器选择GF-2 PMS1,然后输出空间分辨率为0.8 m的多光谱影像数据。
图2 GF-2全色和多光谱融合
图3 GF-2全色和多光谱融合后自然真彩色影像(传感器:PMS1,R:Band 4,G:Band 3,B:Band 2,获取时间:2020-04-03)
2020年11月,实验小组在黄河三角洲国家级自然保护区进行实地调查,利用GPS RTK对不同湿地分布和生物类型进行定位、记录和拍照。以上两类数据构成样本数据和验证数据,用于建立分类器,验证准确性。表3列出了分类训练和验证样本的每个类样本点数。
表3 不同湿地类型分类训练和验证样本的数量
湿地类别 |
训练样本点数 |
验证样本点数 |
salt water |
59357 |
50710 |
farmland |
20783 |
12162 |
river |
8673 |
3047 |
shrub |
2225 |
2695 |
grass |
11401 |
6863 |
suaeda salsa |
269 |
124 |
tidal flat |
18552 |
22408 |
Total |
121260 |
98009 |
图4 分类样本和训练样本选择
2.2 影像预处理
两种国产数据分别使用ENVI 5.6进行影像预处理,其中GF-3数据的预处理过程包括辐射校正、极化滤波、极化矩阵转换;OHS数据的预处理过程包括辐射校正、大气校正。GF-3预处理过程主要通过ENVI 5.6中 SARScape插件完成。OHS预处理过程主要参考欧比特公司的官方数据处理文档《“珠海一号”OHS 高光谱数据预处理操作手册》。
预处理后得到OHS影像可借助软件中“Build 3D Cube”功能绘制光谱立方体。
图5 OHS光谱立方体
2.3 特征提取
2.3.1 GF-3特征提取
GF-3数据的特征集包括极化特征和纹理特征,其中极化特征是通过H-A-α极化分解得到,纹理特征是通过Freeman三分量极化分解和灰度共生矩阵计算得到纹理特征。
基于H-A-α极化分解得到的极化特征(alpha、anisotropy、entropy)如下图所示:
(a)alpha
(b)anisotropy
(c)entropy
图6 GF-3极化特征:(a) alpha;(b) anisotropy;(c) entropy
纹理特征基于ENVI 5.6 中灰度共生矩阵提取得到,操作示意图如下:
图7 GF-3基于灰度共生矩阵提取纹理特征
最后得到生成均值、方差、同质性、对比度、不相似性、熵、能量、相关性这8个纹理特征,结果如下:
图8 GF-3纹理特征假彩色影像(R:均值,G:方差,B:同质性)
2.3.2 OHS特征提取
OHS数据的特征集包括光谱特征和指数特征。其中光谱特征由PCA提取得到,并根据影像中心波长来计算指数特征NDVI和NDWI,如下图所示:
图9 OHS指数特征:(a) NDVI;(b)NDWI
图10 基于OHS数据绘制黄河口七种湿地地物类型的光谱曲线
2.4 特征融合
GF-3特征集经过重采样和地理编码处理之后波段叠加在一起;OHS特征集经过正射校正之后波段叠加在一起。其中波段叠加所采用的工具是ENVI 5.6中的“Build Layer Stack”工具,如下图所示:
图11 波段叠加
波段叠加的过程中可以对影像进行重采样、裁剪、重投影,把影像从地理坐标系投影到投影坐标系(WGS 1984 UTM Zone 50N)。
2.5 分类
融合的特征集基于训练样本集采用ENVI 5.6中的监督分类工具:最大似然、马氏距离、支持向量机分别进行分类。
(a)最大似然法分类
(b)马氏距离法分类
(c)支持向量机法分类
图12 监督分类:(a) 最大似然法分类;(b) 马氏距离法分类;(c) 支持向量机法分类
2.6 精度评估
多组分类结果基于ENVI 5.6中“Confusion Matrix Using Ground Truth ROIs”进行分类后精度评估,操作图如下所示:
图13 分类后精度评估
评估结果包括混淆矩阵、总体精度、Kappa系数、各地物的用户精度、制图精度。
三、研究结果
1、全极化GF-3数据分类效果要远远低于高光谱OHS数据和两种数据协同分类效果。
图14 分类结果图
2、三种分类器中支持向量机(SVM)分类效果最稳定,其次是最大似然法(ML),分类效果稳定性的是马氏距离法(MD)。
图15 各地物占比图
四、分析与讨论
1、精度评定
分类结果精度评定参数:
总体精度(OA)
Kappa系数
图16 分类结果精度图(1)
2、结论:
全极化GF-3数据分类结果精度最差;
高光谱OHS数据分类结果精度很高,高光谱影像数据在湿地地物分类上有很大的优势;
全极化GF-3数据和高光谱OHS数据协同分类的结果精度在两种单一数据源分类的精度上都有一定的提高。
图17 分类结果精度图(2)
GF-3得到的地物精度很低,大部分都在60%以下;
OHS得到的地物精度很高,基本都在90%以上,只有基于MD方法得到的碱蓬和草地的分类精度要低于80%;
协同分类的地物精度很高,基本上都在单一数据源上有一定的提升;
协同分类出现错分的情况要比单一的高光谱OHS数据少。
图18 混淆矩阵
2021遥感应用组三等奖:基于PolSAR和高光谱遥感的黄河口湿地协同分类研究相关推荐
- 【文本分类】基于预训练语言模型的BERT-CNN多层级专利分类研究
·阅读摘要: 本文利用BERT的预训练数据,结合CNN,应用于专利分类.由于专利的特性,本文还提出使用多层级方法来增强模型.(文章发表在<中文信息学报>,核心期刊) ·参考文献: ...
- java 知网 语义 相似度,基于知网语义相似度的中文文本分类研究 论文笔记
基于知网语义相似度的中文文本分类研究 1.传统的文本处理大部分是根据词频和逆向文档频率将文本表示成向量空间模型,实践证明这种模型确实简单高效并且得到了广泛应用,但这种模型表示缺乏对语义的理解,忽略了词 ...
- 基于熵,fft,psd特征,svm的分类研究(肌电信号分类研究)
肌电信号是一种生物电信号,肌电信号(EMG)是众多肌纤维中运动单元动作电位(MUAP)在时间和空间上的,他可以在一定程度上反映神经肌肉的活动:相对于针电极EMG,测量较为方便等优点.因而,SEMG在医 ...
- python遥感图像处理_基于Python的矿山遥感监测系统开发方法
目前,很多学者都是从宏观上讨论遥感和GIS一体化集成的可能性及集成的方法,但这些研究在GIS和RS方向只是对Python集成研究的思路或某一功能的介绍,并没有一个基于Python开发的集成GIS与RS ...
- Matlab高光谱遥感数据处理与混合像元分解实践技术
光谱和图像是人们观察世界的两种方式,高光谱遥感通过"图谱合一"的技术创新将两者结合起来,大大提高了人们对客观世界的认知能力,本来在宽波段遥感中不可探测的物质,在高光谱遥感中能被探测 ...
- Matlab高光谱遥感、数据处理与混合像元分解实践技术应用
光谱和图像是人们观察世界的两种方式,高光谱遥感通过"图谱合一"的技术创新将两者结合起来,大大提高了人们对客观世界的认知能力,本来在宽波段遥感中不可探测的物质,在高光谱遥感中能被探测 ...
- Matlab 高光谱遥感数据处理与混合像元分解
教程从基础.方法.实践三方面对高光谱遥感进行讲解.基础篇,站在学员的角度去理解"高光谱",用大家能听的懂的语言去讲述高光谱的基本概念和理论,帮助学员深入理解这项技术的底层科学机理. ...
- Matlab高光谱遥感、数据处理与混合像元分解及典型案例
站在学员的角度去理解"高光谱",用大家能听的懂的语言去讲述高光谱的基本概念和理论,帮助学员深入理解这项技术的底层科学机理.方法篇,将高光谱技术与MATLAB工具结合起来,采用MAT ...
- Matlab高光谱遥感数据处理
光谱和图像是人们观察世界的两种方式,高光谱遥感通过"图谱合一"的技术创新将两者结合起来,大大提高了人们对客观世界的认知能力,本来在宽波段遥感中不可探测的物质,在高光谱遥感中能被探测 ...
最新文章
- WCF服务一:WCF服务简介
- Jmeter接口测试进阶
- 《大规模分布式系统架构与设计实战》
- BZOJ 2844 | HYSBZ - 2844albus就是要第一个出场——线性基
- 华为手机连电脑_手机、电脑无网高速互传!华为神技逆天
- 产品狗,工作三年,转行AI应该怎样规划?
- stm32f10x单片机进阶--spi使用
- Linux 命令 之 【stat】 查看文件状态。 (包括修改时间)
- Echarts横向的柱状图
- 专栏全年主题合辑-代码中文命名相关实践 2018-11-10
- 如何利用网管软件管控SNMP协议的网络设备
- Photoshop CC 2019魔棒工具的抠图
- 邮件服务器搬家,邮件搬家操作方法
- kettle连接access问题总结
- 跳妹儿学编程之ScratchJr(三):什么是儿童编程语言?初识ScratchJr和Scratch
- 阿里云 root ssh远程登录 及 普通非root用户 ssh远程登录 Ubuntu1604
- 如何装计算机网络驱动,网卡驱动怎么安装,手把手教你电脑网卡驱动怎么安装...
- 微信生成海报 服务器性能,微信小程序生成清晰海报
- MEPR500+电子护照阅读器|识读模块嵌入式安装与应用说明
- 用Python实现原生爬取某牙直播平台数据