10. 遥感图像监督分类

10.1 内容介绍

监督分类,又称训练分类法,用被确认类别的样本像元去识别其他未知类别像元的过程。它就是在分类之前通过目视判读和野外调查,对遥感图像上某些样区中影像地物的类别属性有了先验知识,对每一种类别选取一定数量的训练样本,计算机计算每种训练样区的统计或其他信息,同时用这些种子类别对判决函数进行训练,使其符合于对各种子类别分类的要求, 随后用训练好的判决函数去对其他待分数据进行分类。使每个像元和训练样本作比较,按不同的规则将其划分到和其最相似的样本类,以此完成对整个图像的分类。

遥感影像的监督分类一般包括以下 6 个步骤,如下图10.1.1所示:

10.2 详细操作步骤

10.2.1 类别定义/特征判别

根据分类目的、影像数据自身的特征和分类区收集的信息确定分类系统;对影像进行特征判断,评价图像质量,决定是否需要进行影像增强等预处理。这个过程主要是一个目视查看的过程,为后面样本的选择打下基础。 启动ENVI5.3打开待分类数据。以 R:TM Band 5,G: TM Band 4,B:TM Band 3波段组合显示。 通过目视可分辨六类地物:林地、草地/灌木、耕地、裸地、沙地、其他六类。

10.2.2 样本选择

(1)在图层管理器 Layer Manager 中,can_tmr.img 图层上右键,选择“New Region Of Interest”,打开 Region of Interest (ROI) Tool 面板,下面学习利用选择样本。

【1】在 Region of Interest (ROI) Tool 面板上,设置以下参数:

图10.2.1 Region of Interest (ROI) Tool 面板上设置样本参数

【2】默认 ROIs 绘制类型为多边形,在影像上辨别林地区域并单击鼠标左键开始绘制多边形样本,一个多边形绘制结束后,双击鼠标左键或者点击鼠标右键,选择 Complete and Accept Polygon,完成一个多边形样本的选择;

【3】同样方法,在图像别的区域绘制其他样本,样本尽量均匀分布在整个图像上;

【4】这样就为林地选好了训练样本。

(2)在图像上右键选择 New ROI,或者在 Region of Interest (ROI) Tool 面板上,选择工具。重复“林地”样本选择的方法,分别为草地/灌木、耕地、裸地、沙地、其他 5 类选择样本;

(3)训练样本的选择

【1】计算样本的可分离性。在 Region of Interest (ROI) Tool 面板上,选择 Option>Compute ROI Separability,在 Choose ROIs 面板,将几类样本都打勾,点击 OK;

【2】表示各个样本类型之间的可分离性,用 Jeffries-Matusita, Transformed Divergence 参数表示,这两个参数的值在 0~2.0 之间,大于 1.9 说明样本之间可分离性好,属于合格样本;小于 1.8,需要编辑样本或者重新选择样本;小于 1,考虑将两类样本合成一类样本。如图10.2.2

图10.2.2 Merge ROIs 面板

(4)在图层管理器中,选择 Region of interest ,点击右键,save as,保存为.xml 格式的样本文件。

10.2.3 分类器选择

根据分类的复杂度、精度需求等确定哪一种分类器。目前 ENVI 的监督分类可分为基于传统统计分析学的,包括平行六面体、最小距离、马氏距离、最大似然,基于神经网络的,基于模式识别,包括支持向量机、模糊分类等,针对高光谱有波谱角(SAM),光谱信息散度,二进制编码。下面是几种分类器的简单描述。

(1)平行六面体(Parallelepiped)

根据训练样本的亮度值形成一个 n 维的平行六面体数据空间,其他像元的光谱值如果落在平行六面体任何一个训练样本所对应的区域,就被划分其对应的类别中。

(2)最小距离(Minimum Distance)

利用训练样本数据计算出每一类的均值向量和标准差向量,然后以均值向量作为该类在特征空间中的中心位置,计算输入图像中每个像元到各类中心的距离,到哪一类中心的距离最小,该像元就归入到哪一类。

(3)马氏距离(Mahalanobis Distance)

计算输入图像到各训练样本的协方差距离(一种有效的计算两个未知样本集的相似度的方法),最终技术协方差距离最小的,即为此类别。

(4)最大似然(Maximum Likelihood)

假设每一个波段的每一类统计都呈正态分布,计算给定像元属于某一训练样本的似然度,像元最终被归并到似然度最大的一类当中。

(5)神经网络(Neural Net)

指用计算机模拟人脑的结构,用许多小的处理单元模拟生物的神经元,用算法实现人脑的识别、记忆、思考过程。

(6)支持向量机(Support Vector Machine)

支持向量机分类(Support Vector Machine 或 SVM)是一种建立在统计学习理论(Statistical Learning Theory 或 SLT)基础上的机器学习方法。SVM 可以自动寻找那些对分类有较大区分能力的支持向量,由此构造出分类器,可以将类与类之间的间隔最大化,因而有较好的推广性和较高的分类准确率。

(7)波谱角(Spectral Angle Mapper)

它是在 N 维空间将像元与参照波谱进行匹配,通过计算波谱间的相似度,之后对波谱之间相似度进行角度的对比,较小的角度表示更大的相似度。

10.2.4 影像分类

基于传统统计分析的分类方法参数设置比较简单,在 Toolbox/Classification/Supervised Classification 能找到相应的分类方法。这里选择支持向量机分类方法。在 toolbox 中选择 /Classification/Supervised Classification/Support Vector Machine Classification,选择待分类影像,点击 OK,按照默认设置参数输出分类结果。如图10.2.3和10.2.4

图10.2.3 支持向量机分类器参数设置

图10.2.4 支持向量机分类结果

10.2.5 分类后处理

包括更改类别颜色、分类后统计、小斑块处理、栅矢转换等

10.2.6 精度验证

对分类结果进行评价,确定分类的精度和可靠性。有两种方式用于精度验证:一是混淆矩阵,二是 ROC 曲线,比较常用的为混淆矩阵,ROC 曲线可以用图形的方式表达分类精度,比较抽象。

真实参考源可以使用两种方式:一是标准的分类图,二是选择的感兴趣区(验证样本区)。两种方式的选择都可以通过主菜单->Classification->Post Classification->Confusion Matrix 或者 ROC Curves 来选择。 真实的感兴趣区验证样本的选择可以是在高分辨率影像上选择,也可以是野外实地调查获取,原则是获取的类别参考源的真实性。由于没有更高分辨率的数据源,本例中就把原分类的TM 影像当作是高分辨率影像,在上面进行目视解译得到真实参考源。

  1. 在 Data Manager 中,分类样本上右键选择 Close,将分类样本从软件中移除
  2. 直接利用 ROI 工具,跟分类样本选择的方法一样,即重复第二步,在 TM 图上选择
  3. 类验证样本。如图10.2.5

图10.2.5 选择验证样本

(4)在 Toolbox 中,选择/Classification/Post Classification/Confusion Matrix Using Ground Truth ROIs,选择分类结果,软件会根据分类代码自动匹配,如不正确可以手动更改。点击 OK 后选择报表的表示方法(像素和百分比),点击 OK,就可以得到精度报表如图10.2.6和10.2.7。

图10.2.6 验证操作面板

图10.2.7 分类精度评价混淆矩阵

混淆矩阵中的评价指标:

等于被正确分类的像元总和除以总像元数。被正确分类的像元数目沿着混淆矩阵的对角线分布,总像元数等于所有真实参考源的像元总数,如本次精度分类精度表中的 Overall Accuracy = (1849/2346)  78.8150%。

它是通过把所有真实参考的像元总数(N)乘以混淆矩阵对角线(XKK)的和,再减去某一类中真实参考像元数与该类中被分类像元总数之积之后,再除以像元总数的平方减去某一类中真实参考像元总数与该类中被分类像元总数之积对所有类别求和的结果。Kappa 计算公式

指被分为用户感兴趣的类,而实际属于另一类的像元,它显示在混淆矩阵里面。本例中,林地有 419 个真实参考像元,其中正确分类 265,12 个是其他类别错分为林地(混淆矩阵中林地一行其他类的总和),那么其错分误差为 12/419=2.9%。

指本身属于地表真实分类,当没有被分类器分到相应类别中的像元数。如在本例中的耕地类,有真实参考像元 465 个,其中 462 个正确分类,其余 3 个被错分为其余类(混淆矩阵中耕地类中一列里其他类的总和),漏分误差为 3/465=0.6%

是指分类器将整个影像的像元正确分为 A 类的像元数(对角线值)与 A 类真实参考总数(混淆矩阵中 A 类列的总和)的比率。如本例中林地有 419 个真实参考像元,其中 265 个正确分类,因此林地的制图精度是 265/419=63.25%。

是指正确分到 A 类的像元总数(对角线值)与分类器将整个影像的像元分为 A 类的像元总数(混淆矩阵中 A 类行的总和)比率。如本例中林地有 265 个正确分类,总共划分为林地的有277,所以林地的用户精度是 265/277=95.67%。

ENVI学习总结(十)——遥感图像监督分类相关推荐

  1. RSE2020/云检测:基于弱监督深度学习的高分辨率遥感图像精确云检测

    RSE2020/云检测:Accurate cloud detection in high-resolution remote sensing imagery by weakly supervised ...

  2. 基于深度学习的高分辨率遥感图像目标检测技术目前的研究现状

    参考   基于深度学习的高分辨率遥感图像目标检测技术目前的研究现状 - 云+社区 - 腾讯云 目录 一.概述 二.通用目标检测方法 1.类不平衡问题 2.复杂背景 3.目标的尺度变化 4.特殊视角 5 ...

  3. 目标检测YOLO实战应用案例100讲-基于深度学习的光学遥感图像目标检测及价值评估

    目录 基于深度学习的光学遥感图像目标检测方法研究 传统的目标检测方法

  4. 论文总结(一)—基于深度学习的普通遥感图像质量改进

    论文名称: A universal remote sensing image quality improvement method with deep learning 链接:https://ieee ...

  5. ENVI:如何进行遥感图像的分类?(决策树模型)

    目录 01 决策树模型的简单介绍 02 实操 2.1 加载遥感影像以及对应的DEM影像 2.2 新建决策树 2.2.1 找到 新建决策树 工具并点击打开 2.2.2 决策树的编辑 2.2.3 NDVI ...

  6. OpenCV图像处理学习二十,图像直方图均衡化原理与实现

    一.图像直方图的概念 图像直方图,是指对整个图像在灰度范围内的像素值(0~255)统计出现频率次数,据此生成的直方图,称为图像直方图.直方图反映了图像灰度的分布情况,是图像的统计学特征.图像的灰度直方 ...

  7. win10+Python3.7.3+OpenCV3.4.1入门学习(十二 图像轮廓)————12.4 Hu矩

    文章目录 12.4 Hu矩 12.4.1 Hu矩函数 12.4.2 形状匹配 Python版本是Python3.7.3,OpenCV版本OpenCV3.4.1,开发环境为PyCharm 12.4 Hu ...

  8. win10+Python3.7.3+OpenCV3.4.1入门学习(十二 图像轮廓)————12.7 利用形状场景算法比较轮廓

    Python版本是Python3.7.3,OpenCV版本OpenCV3.4.1,开发环境为PyCharm 文章目录 12.7 利用形状场景算法比较轮廓 12.7.1 计算形状场景距离 12.7.2 ...

  9. 2022TGRS/云检测:用于遥感图像云检测的无监督域不变特征学习Unsupervised Domain-Invariant Feature Learning for Cloud Detection

    2022TGRS/云检测:Unsupervised Domain-Invariant Feature Learning for Cloud Detection of Remote Sensing Im ...

  10. 遥感图像预处理与土地利用动态监测

    实验报告目录 1.实验目的 对池州市 2013.2017 年 Landsat 8 OLI 的遥感图像进行土地利用分类,对分类结果进行修正和精度评价,并对各类土地利用类型的面积进行统计.基于计算机分类结 ...

最新文章

  1. Homography matrix(单应性矩阵)在广告投放中的实践
  2. Android IJKPlayer缓冲区设置以及播放一段时间出错解决方案
  3. TypeScript学习笔记3:运算符
  4. boost::geometry::util::is_implemented2用法的测试程序
  5. 【SpringBoot零基础案例09】【IEDA 2021.1】SpringBoot将核心配置文件中的自定义配置映射到一个对象
  6. hbase过滤器查询
  7. 深度学习神经网络:改善与优化的方法
  8. 软件包管理器(树链剖分)
  9. python 引用文件中的类 报错_Python学习笔记7 头文件的添加规则(转载)
  10. udt编写高性能服务器,基于UDT协议的Oracle数据库远程备份的设计和实现
  11. Debug Current Instruction Pointer 删除debug蓝色小箭头
  12. PHP Ctype函数(转)
  13. 博客开园了~~~~~~
  14. oracle打开dmp文件乱码,oracle中导入dmp字符乱码分析和解决方案
  15. linux的yum是什么意思,yum是什么?(linux命令)
  16. Python进阶实战 1.0 :爬取豆瓣TOP250电影排行榜
  17. 华硕路由器(RT-AC86U)架设简易git服务器
  18. 【数字IC验证快速入门】11、Verilog TestBench(VTB)入门
  19. wp/wordpress文章页面添加阅读量/点击量,后台并显示阅读量
  20. 【数据压缩作业9】H.264码流分析仪

热门文章

  1. 如何设计一个秒杀系统?
  2. 公文流转 java_javaweb 公文流转系统制作
  3. 100个超好用的Excel快捷键
  4. win10 -- 注册机认为是病毒,死活不让下载 --那就关了它的自带杀毒软件
  5. 麒麟 linux下安装显卡驱动,中标麒麟7安装英伟达显卡驱动过程简述
  6. 泊松分布、二项分布与正态分布
  7. svn —— 版本回退
  8. P2P网贷谋求“去担保”:引入保险 收益下滑
  9. 宇视NVR如何宇视云解绑?
  10. libtorrent-bittorrent