第8章遥感图像分类 76页PPT文档
合集下载
测绘学概论第8章 遥感科学与技术

遥
感
➢ 空间分辨率:通常指一个像素对应
图
地面的实际大小。一般遥感图像分
像
辨率指的是地面分辨率;
中
➢ 光谱分辨率:成象范围内波谱带数
的
目;
分
辨
➢ 时间分辨率:重复获取某地区图像
率
的周期;
➢ 温度分辨率(热红外):可探测的
温度变化幅度。
应用需求
资源
农业
环境
森林
国防
高光谱
多光谱
交通
城市
光谱分辨率
全色
100 m
公司 ESA NOAA
发射时间 1995 1994-1995
NASDA 1995 ISRO 1996-1997 中国 1997
位置 0o
75135oW 135oE 75oE 105oE
8.3 遥感信息获取
太阳同步极轨气象卫星系统
系统 公司
NOAA- NOAA 14
NOAA-K NOAA
FY-
中国
1A/1B
非图像方式(主动式和被动式)
雷达高度计 合成孔径雷达 微波辐射计 红外辐射计
…
8.2 遥感的电磁波谱
电磁波谱
8.2 遥感的电磁波谱
遥感技术使用的电磁波分类名称和波长范围
名称
波长范围
紫外线 可见光
近红外
100 A°~0.4μm 0.4~0.7μm 0.76~3.0μm
紫 0.38~0.43μm 蓝 0.43~0.47μm 青 0.47~0.50μm
radarsat
Landsat
SPOT
ers SEASAT
JERS
NOAA
航天遥感传感器搭载的主要平台是卫星。 上图是目前国外常用的遥感卫星。
第8章遥感图像自动识别分类

• 3、将男错判为女,或将女错判为男,其产 生消极后果没有区别,但现实中“弃真”造 成的消极后果一般要大于“存伪”错误。
• 本健康但错判断为癌症:存伪,虚惊一场;
第八章 遥感图像自动识别分类 §8.1 基础知识
三、模式识别系统
1、数据获取: 图像,波形,物理参量 2、预处理: 去噪、增强、退化复原 3、特征提取与选择 4、分类决策
4、其它方法:模拟退火法、遗传算法等。
12/13
第八章 遥感图像自动识别分类 §8.3 分类器设计
• 本课主要内容 • 分类器(判别准则)概念 • 最小错误率分类器(Bayes 准则) • 线性判别分类器 (Fisher准则) • 非线性判别分类器(最短距离法) • 本课重点内容 • Bayes判别准则 • 最短距离分类器
主分量变换的优良特性 1、变换后Y的协方差阵是对角阵,表明新特征矢量直接 彼此不相关 2、变换后,信息主要集中在前几个主成分中,根据统 计,对于landsat MSS四个波段的影像经KL变换后, PC1占90%的总信息量,PC2占7%的总信息量,PC3和 PC4共占3%的总信X息2 量。 Y2(第二分量)
第八章 遥感图像自动识别分类 §8.4 监督/非监督分类
• 上节主要内容 • 贝叶斯分类器 • 线性分类器(以Fisher准则为例) • 非线性分类器(以最短距离方法为例) • 本节主要内容 • 监督分类的思想 • 监督/非监督分类的区别 • 监督分类步骤 • K-均值聚类算法
第八章 遥感图像自动识别分类 §8.4 监督/非监督分类
参数估计:点估计/区间估计/最大似然估计 非参数估计:Parzen窗法、k-近邻估计法 Bayes分类器的理论作用大于实际作用 能否直接利用特征设计分类器呢?
第八章 遥感图像自动识别分类 §8.3 分类器设计
• 本健康但错判断为癌症:存伪,虚惊一场;
第八章 遥感图像自动识别分类 §8.1 基础知识
三、模式识别系统
1、数据获取: 图像,波形,物理参量 2、预处理: 去噪、增强、退化复原 3、特征提取与选择 4、分类决策
4、其它方法:模拟退火法、遗传算法等。
12/13
第八章 遥感图像自动识别分类 §8.3 分类器设计
• 本课主要内容 • 分类器(判别准则)概念 • 最小错误率分类器(Bayes 准则) • 线性判别分类器 (Fisher准则) • 非线性判别分类器(最短距离法) • 本课重点内容 • Bayes判别准则 • 最短距离分类器
主分量变换的优良特性 1、变换后Y的协方差阵是对角阵,表明新特征矢量直接 彼此不相关 2、变换后,信息主要集中在前几个主成分中,根据统 计,对于landsat MSS四个波段的影像经KL变换后, PC1占90%的总信息量,PC2占7%的总信息量,PC3和 PC4共占3%的总信X息2 量。 Y2(第二分量)
第八章 遥感图像自动识别分类 §8.4 监督/非监督分类
• 上节主要内容 • 贝叶斯分类器 • 线性分类器(以Fisher准则为例) • 非线性分类器(以最短距离方法为例) • 本节主要内容 • 监督分类的思想 • 监督/非监督分类的区别 • 监督分类步骤 • K-均值聚类算法
第八章 遥感图像自动识别分类 §8.4 监督/非监督分类
参数估计:点估计/区间估计/最大似然估计 非参数估计:Parzen窗法、k-近邻估计法 Bayes分类器的理论作用大于实际作用 能否直接利用特征设计分类器呢?
第八章 遥感图像自动识别分类 §8.3 分类器设计
《遥感图像分类》课件

特征变换
将原始特征进行变换,生成新的特征,以更好地 反映地物类别之间的差异。
分类器设计
监督分类
利用已知样本的训练集设计分类器,对未知样本进行分类。
非监督分类
对未知样本进行聚类分析,将相似的样本归为同一类。
混合分类
结合监督分类和非监督分类的优势,提高分类精度和稳定性。
分类结果评价
精度评价
通过比较分类结果与实际地物类别, 计算分类精度、混淆矩阵等指标。
THANKS
感谢观看
分类器。
多源遥感数据融合问题
多源遥感数据融合可以提高分类精度和可靠性,但同时也带 来了数据匹配、融合算法选择等问题。
解决多源遥感数据融合问题的策略包括使用先进的融合算法 ,如基于深度学习的融合方法,以及优化数据匹配方法。
遥感图像分类技术的发展趋势
01
遥感图像分类技术正朝着高精度、高效率和自动化的方向发展 。
可靠性评价
评估分类结果的稳定性、可靠性以及 抗干扰能力。
应用价值评价
根据分类结果在实际应用中的价值, 如土地利用、资源调查、环境监测等
,对分类方法进行综合评价。
04
CATALOGUE
遥感图像分类的挑战与展望
数据质量问题
遥感图像常常受到噪声、失真和 模糊等影响,导致数据质量下降
。
数据质量问题还表现在不同传感 器获取的图像之间的差异,以及 不同时间获取的图像之间的变化
遥感图像分类的应用
遥感图像分类在多个领域有广泛应用,如环境保护、城市规划、资源调查、军事 侦察等。
通过遥感图像分类,可以快速获取大范围的地物信息,为相关领域的决策提供科 学依据。
02
CATALOGUE
遥感图像分类的方法
将原始特征进行变换,生成新的特征,以更好地 反映地物类别之间的差异。
分类器设计
监督分类
利用已知样本的训练集设计分类器,对未知样本进行分类。
非监督分类
对未知样本进行聚类分析,将相似的样本归为同一类。
混合分类
结合监督分类和非监督分类的优势,提高分类精度和稳定性。
分类结果评价
精度评价
通过比较分类结果与实际地物类别, 计算分类精度、混淆矩阵等指标。
THANKS
感谢观看
分类器。
多源遥感数据融合问题
多源遥感数据融合可以提高分类精度和可靠性,但同时也带 来了数据匹配、融合算法选择等问题。
解决多源遥感数据融合问题的策略包括使用先进的融合算法 ,如基于深度学习的融合方法,以及优化数据匹配方法。
遥感图像分类技术的发展趋势
01
遥感图像分类技术正朝着高精度、高效率和自动化的方向发展 。
可靠性评价
评估分类结果的稳定性、可靠性以及 抗干扰能力。
应用价值评价
根据分类结果在实际应用中的价值, 如土地利用、资源调查、环境监测等
,对分类方法进行综合评价。
04
CATALOGUE
遥感图像分类的挑战与展望
数据质量问题
遥感图像常常受到噪声、失真和 模糊等影响,导致数据质量下降
。
数据质量问题还表现在不同传感 器获取的图像之间的差异,以及 不同时间获取的图像之间的变化
遥感图像分类的应用
遥感图像分类在多个领域有广泛应用,如环境保护、城市规划、资源调查、军事 侦察等。
通过遥感图像分类,可以快速获取大范围的地物信息,为相关领域的决策提供科 学依据。
02
CATALOGUE
遥感图像分类的方法
遥感图像分类ppt课件

– 假设遥感图像有K个波段,则(i,j)位置的像素在
每个波段上的灰度值可以构成表示为X=(x1,
T
5
8.1 概述
• 8.1.2 分类方法
– 根据是否需要分类人员事先提供已知类别及其 训练样本,对分类器进行训练和监督,可将遥 感图像分类方法划分为监督分类和非监督分类。
– 事先己经知道类别的部分信息(即类别的先验知 识),对未知类别的样本进行分类的方法称之为 监督分类(Supervised Classification)。事先没 有类别的先验知识,对未知类别的样本进行分 类的方法称之为非监督分类(Unsupervised Classification)
14
8.4 非监督分类
• 非监督分类,是指人们事先对分类过程不
施加任何的先验知识,仅凭据遥感影像地 物的光谱特征的分布规律,随其自然地进 行盲目的分类。其分类的结果,只是对不 同类别进行了区分,并不确定类别的属性, 其属性是通过事后对各类的光谱响应曲线 进行分析,以及与实地调查相比较后确定 的。
• 非监督分类的理论依据:遥感图像上的同
4
8.1 概述
• 8.1.1 基本原理
– 同类地物在相同的条件下(光照、地形等)应该 具有相同或相似的光谱信息和空间信息特征。 不同类的地物之间具有差异根据这种差异,将 图像中的所有像素按其性质分为若干个类别 (Class)的过程,称为图像的分类。
– 遥感图像分类以每个像素的光谱数据为基础进 行。
9
8.2 相似性度量
3.马氏(Mahalanobis)距离
马氏距离是一种加权的欧氏距离,它通 过协方差矩阵来考虑变量的相关性。这 是由于在实际中,各点群的形状是大小 和方向各不相同的椭球体,如图所示, 尽管K点距MA的距离DA比距MB的距离 DB小,即DA<DB ,但由于B点群比A点 群离散得多,因而把K点划入B类更合 理。加权可以这样理解,计算的距离与 各点群的方差有关。方差愈大,计算的 距离就愈短。如果各个点群具有相同的 方差,则马氏距离是欧氏距离的平方。
遥感图像分类分析PPT课件

“Max stdev from Mean”文本框中输入用于限定相对于均值的标 准差的大小。 ➢要 为 每 一 类 别 设 置 不 同 的 阈 值 :
➢A. 在类别列表中,点击想要设置不同阈值的类别。 ➢B. 点击“Multiple Values”来选择它。 ➢C. 点击“Assign Multiple Values”按钮。 ➢D. 在出现的对话框中,点击一个类别选中它,然后在对话框底部的文本
➢选择Classification > Super vised > Maximum Likelihood ➢设 定 似 然 度 的 阈 值 , 范 围 0 - 1 ➢数 据 比 例 系 数 : 这 个 比 例 系 数 是 一 个 比 值 系 数 , 用 于 将 整 型 反 射 率 或 辐 射 率
多辅助方法,如上面的可以显示不同的假彩色合成窗口,也可以进行主成分分析后进行 假彩色合成,由于去除了波段间的相关性,不同地物区分的更加明显;还可以借助 Google Earth辅助解译。
第8页/共53页
ENVI/IDL
6.2 监督分类
➢第二种方法,在散点图上进行选择 ➢(1)在主图像上,选择tools > 2D scatter plots,将1波段作为X,4波段作为Y,原理
第2页/共53页
ENVI/IDL
6.1 分类类型
❖1、监督分类
❖监 督 分 类 : 又 称 训 练 分 类 法 , 用 被 确 认 类 别 的 样 本 像 元 去 识 别 其他未知类别像元的过程。它就是在分类之前通过目视判读和野 外调查,对遥感图像上某些样区中影像地物的类别属性有了先验 知识,对每一种类别选取一定数量的训练样本,计算机计算每种 训练样区的统计或其他信息,同时用这些种子类别对判决函数进 行训练,使其符合于对各种子类别分类的要求,随后用训练好的 判决函数去对其他待分数据进行分类。使每个像元和训练样本作
➢A. 在类别列表中,点击想要设置不同阈值的类别。 ➢B. 点击“Multiple Values”来选择它。 ➢C. 点击“Assign Multiple Values”按钮。 ➢D. 在出现的对话框中,点击一个类别选中它,然后在对话框底部的文本
➢选择Classification > Super vised > Maximum Likelihood ➢设 定 似 然 度 的 阈 值 , 范 围 0 - 1 ➢数 据 比 例 系 数 : 这 个 比 例 系 数 是 一 个 比 值 系 数 , 用 于 将 整 型 反 射 率 或 辐 射 率
多辅助方法,如上面的可以显示不同的假彩色合成窗口,也可以进行主成分分析后进行 假彩色合成,由于去除了波段间的相关性,不同地物区分的更加明显;还可以借助 Google Earth辅助解译。
第8页/共53页
ENVI/IDL
6.2 监督分类
➢第二种方法,在散点图上进行选择 ➢(1)在主图像上,选择tools > 2D scatter plots,将1波段作为X,4波段作为Y,原理
第2页/共53页
ENVI/IDL
6.1 分类类型
❖1、监督分类
❖监 督 分 类 : 又 称 训 练 分 类 法 , 用 被 确 认 类 别 的 样 本 像 元 去 识 别 其他未知类别像元的过程。它就是在分类之前通过目视判读和野 外调查,对遥感图像上某些样区中影像地物的类别属性有了先验 知识,对每一种类别选取一定数量的训练样本,计算机计算每种 训练样区的统计或其他信息,同时用这些种子类别对判决函数进 行训练,使其符合于对各种子类别分类的要求,随后用训练好的 判决函数去对其他待分数据进行分类。使每个像元和训练样本作
遥感图像类型与特性PPT课件

3.投影性质 动态多中心投影
(385个投影中心)
185 km
480 m
185 km
动态多中心投影的影像亦存在像点位移。像 点位移的大小与卫星平台运行高度、地表起伏高 差以及扫描角有关。
由于卫星平台运行高度较高,总扫描角较小 (11.56°),所以当地表相对高差较小且成图精 度要求不高时,可将图像近似看作是垂直投影 (正射投影)。当成图精度要求较高时(如 TM 1 : 5万成图),应根据DEM进行几何精校正。
像元
像元数目:
TM 1-5、7 6166×6166
TM 6
1542×1542
像元大小(对应地面面积):
TM 1-5、7 30m×30m
TM 6
120m×120m
像元是在扫描成像过程中通过采样而形成的 采样点,是扫描影像中最小可分辨面积,也是进 行计算机处理时的最基本单元。
② 量化
图像函数数值离散化
Landsat/TM
(专题制图仪)
TM1 0.45~0.52 μm (30m)
Landsat/MSS
(多光谱扫描仪)
————
TM2 0.52~0.60 μm (30m) MSS4 0.5~0.6 μm (79m)
TM3 0.63~0.69 μm (30m) MSS5 0.6~0.7 μm (79m)
按遥感器成像方式和工作波段划分常规摄影像片非常规摄影像片光机扫描图像固体自扫描图像天线扫描图像黑白全色像片可见光天然彩色像片可见光黑白红外像片近红外彩色红外像片部分可见光近红外紫外像片紫外多波段像片紫外近红外全景像片可见光近红外红外扫描图像中热远红外多波段扫描图像紫外远红外超多波段扫描图像可见光远红外固体自扫描图像可见光近红外成像雷达图像微波航空摄影像片航天摄影像片热红外图像landsattm图像成像波谱仪图像spothrv图像sar图像1
遥感分类PPT课件

第10页/共76页
概念及原理 特点和原则 分类顺序 分类统计量
10
一、遥感影像分类
概念及原理 特点和原则 分类顺序 分类统计量
混合像元的存在常使得同类地物的特征向 量也不尽相同。
30m R
10m R
11
第11页/共76页
一、遥感影像分类
概念及原理 特点和原则 分类顺序 分类统计量
影像分类的关键问题之一是选择适当的分 类规则(或叫分类器),通过分类器把影像数据划 分为尽可能符合实际情况的不同类别。
迭代次数进行控制的参数,如控制迭代的总次数、每一类别最小像 元数、类别的标准差、比较相邻两次迭代效果以及可以合并的最大 类别数等,在整个迭代过程中,不仅每个像元的归属类别在调整, 而且类别总数也在变化。在用计算机编制分类程序时,初始聚类中 心可按如下方式确定:设初始类别数为n,这样共有n个初始聚类 中心,求出图像的均值M和方差σ,按下式可求出初始聚类中心:
52
第52页/共76页
三、监督分类
(1)
最小距离分类法 最大似然比分类法
53
第53页/共76页
三、监督分类
(2)
最小距离分类法 最大似然比分类法
第54页/共76页
11
54
三、监督分类
(2)
最小距离分类法 最大似然比分类法
第55页/共76页
11
55
三、监督分类
(3)
最小距离分类法 最大似然比分类法
解决由于地形因素造成的“同 物异谱”现象的一种有效的方 法是采用“同类多组法”来选 取训练样本。
68
第68页/共76页
四、遥感影像分类处 理中的几个问题
训练场地和样本选择 地形因素影响 混合像元问题 特征变量的选择问题 空间信息的应用问题 后期处理问题
概念及原理 特点和原则 分类顺序 分类统计量
10
一、遥感影像分类
概念及原理 特点和原则 分类顺序 分类统计量
混合像元的存在常使得同类地物的特征向 量也不尽相同。
30m R
10m R
11
第11页/共76页
一、遥感影像分类
概念及原理 特点和原则 分类顺序 分类统计量
影像分类的关键问题之一是选择适当的分 类规则(或叫分类器),通过分类器把影像数据划 分为尽可能符合实际情况的不同类别。
迭代次数进行控制的参数,如控制迭代的总次数、每一类别最小像 元数、类别的标准差、比较相邻两次迭代效果以及可以合并的最大 类别数等,在整个迭代过程中,不仅每个像元的归属类别在调整, 而且类别总数也在变化。在用计算机编制分类程序时,初始聚类中 心可按如下方式确定:设初始类别数为n,这样共有n个初始聚类 中心,求出图像的均值M和方差σ,按下式可求出初始聚类中心:
52
第52页/共76页
三、监督分类
(1)
最小距离分类法 最大似然比分类法
53
第53页/共76页
三、监督分类
(2)
最小距离分类法 最大似然比分类法
第54页/共76页
11
54
三、监督分类
(2)
最小距离分类法 最大似然比分类法
第55页/共76页
11
55
三、监督分类
(3)
最小距离分类法 最大似然比分类法
解决由于地形因素造成的“同 物异谱”现象的一种有效的方 法是采用“同类多组法”来选 取训练样本。
68
第68页/共76页
四、遥感影像分类处 理中的几个问题
训练场地和样本选择 地形因素影响 混合像元问题 特征变量的选择问题 空间信息的应用问题 后期处理问题
- 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。
12
8.3 工作流程
图像的预处理 选择分类方法
1.确定工作范围 4.辐射校正
2.多源图像的几何配准 5.几何精纠正
3.噪声处理 6.多图像融合
特征选择和提取
选择合适的分类 参数进行分类
特征选择是从众多特征中挑选出可以参加分类运算的若干个特征。
特征提取是在特征选择以后,利用特征提取算法(如主成分分析算法) 从原始特征中求出最能反映地物类别性质的一组新特征。
16
8.4 非监督分类
主要过程如下: (1)确定初始类别参数,即确定最初类别数和类别中心 (集群中心) (2)计算每一个像元所对应的特征矢量与各集群中心的 距离。 (3)选与中心距离最短的类别作为这一矢量的所属类别。 (4)计算新的类别均值向量。 (5)比较新的类别均值与原中心位置上的变化,若位置 发生了改变,则以新的类别均值作为聚类中心,再从第2 步开始重复,进行反复迭代操作。 (6)如果聚类中心不再变化,计算停止。
式中:m为已形成的初始类别数。
4)在此基础上,再对所有各剩余点的最小距离进行相互比较,
取其中最大者,并选择与该最大的最小距离相应的抽样点(例
如
X 11
)作为一个新的初始类别中心点(例如
)。
Z3 X11
重复以上步骤,直到初始类别的个数达到需要的个数为止。
25
8.4 非监督分类
8.4.1 初始类别参数的选定 4.局部直方图峰值定心法
22
8.4 非监督分类
8.4.1 初始类别参数的选定
3.最大最小距离选心法
该法的选心原则是使各初始类别之间,尽可能地保持远距
离,为做到这一点,首先在整幅图像中按一定方式(如等间隔
方式)获取一个抽样点像素集合
{
X
}
{X}{X1X2X3 Xn}
式中n为抽样个数
23
8.4 非监督分类
Z q im ii[2 (q 1 )/Q ( 1 ) 1 ]
例如, Q=4时
Z 1i
mi
;
i
Z 3i
m
+
i
1 3
;
i
Z 2i
mi
1 3
i;
Z 4i
m
+
i
;
i
用作图法不难证明,该法所选定的初始类别中心通过整体集群
中心的一条直线均匀散布的,并基本上包括在整体集群的范围 之内
运用距离判别函数时,要求各个类别点群的中心位置己知。 对于光谱特征空间中的任一点k,计算它到各类中心点的距离 d。若di<dj,则k像素属于i类而不属于j类,此处,判别准则 为di<dj
8
8.2 相似性度量
1.绝对距离
dik为当前像素i到类k的距离, P为波段数, xij为像素i在j波段的像素值, Mjk为类k在波段j的均值。
特征选取的原则:尽可能选取数量适当有助于增加同类地物的类聚性 (即类内个体间离散性小)、能扩大地物类别之间差异性的特征变量。
分类后处理
成果输出
13
8.4 非监督分类 根据事先指定的某一准则,让计算机自动进行判
别归类,无需人为干预,分类后需确定地面类别。在非 监督分类中,先确定光谱可分的类别(光谱类)然后定 义它们的信息类。
17
8.4 非监督分类
8.4.1 初始类别参数的选定 8.4.2 常用非监督分类算法
18
8.4 非监督分类
8.4.1 初始类别参数的选定 初始类别参数是指基准类别集群中心(数学期望Mi),以 及集群分布的协方差矩阵Σi 。因为无论采用何种判别函数, 都要预先确定其初始类别的参量。以下为几种常用的方法。
7
8.2 相似性度量
距离是常用的相似性度量的概念。分类是确定像素距离哪个 点群中心较近,或落入哪个点群范围可能性大的问题。像素 与点群的距离越近,那么,属于该点群的可能性越高。按照 一定的准则,当距离小于一定值时,像素被划分给最近的点 群。每个点群为一个类。
同一类别中点间的距离一般来说比不同类别点间距离要小。 也可以认为,一个点属于某一类,一定与这个类别中心的距 离比与其他类别中心的距离小。因此,在点群(团)中心已知的 情况下,以每个点与点群中心的距离作为判定的准则,就可 以完成分类工作。
遥感数字图像处理
林金堂 闽江学院地理科学系
1
第8章 遥感图像分类
2
第8章 遥感图像分类
8.1 概述 8.2 相似性度量 8.3 工作流程 8.4 非监督分类 8.5 监督分类 8.6 其他分类方法 8.7 分类后处理 8.8 分类精度分析
3
8.1 概述
8.1.1 基本原理 8.1.2 分类方法
5
8.1 概述
8.1.2 分类方法 根据是否需要分类人员事先提供已知类别及其训练样本, 对分类器进行训练和监督,可将遥感图像分类方法划分为 监督分类和非监督分类。 事先己经知道类别的部分信息(即类别的先验知识),对未 知类别的样本进行分类的方法称之为监督分类(Supervised Classification)。事先没有类别的先验知识,对未知类别的 样本进行分类的方法称之为非监督分类(Unsupervised Classification) 根据分类使用的统计数学方法可以分为随机统计方法和模 糊数学方法分类。前者以随机数学理论为基础,包括K-均 值分类、最大相似性分类等。后者以模糊数学理论为基础, 主要是模糊分类。
N
( xij
j 1
mi
)
2
1 N
mi
N
xij
j 1
式中,i为波段号;j为像素点号;x ij 为像素i在第j波段的亮度 值;n为波段数;N为像素总数。
21
8.4 非监督分类
8.4.1 初始类别参数的选定 2.总体直方图均匀选心法
现假设我们需要Q个初始类别,每个初始类别集群中心可 按下式确定:
则
Z2 X 7
24
8.4 非监督分类
8.4.1 初始类别参数的选定 3.最大最小距离选心法
3)对于剩余的每个抽样点,计算它到已有各初始类别中心的距
离 Dij(i1,2,,m) m为已有初始类别数,并取其中的最小距离
作为该点的代表距离 D j
D j mD i1jn ,D 2 { j, ,D m }j
• 1.像素光谱特征的比较法 • 2.总体直方图均匀选心法 • 3.最大最小距离选心法 • 4.局部直方图峰值定心法
19
8.4 非监督分类
8.4.1 初始类别参数的选定 1.像素光谱特征的比较法
• 首先,在遥感图像中定义一个抽样集,它可以是整幅图像的 所有像素,但通常是按一定间隔抽样的像素。
在进行聚类分析时,首先要确定基准类别的参量。在非监督 分类的情况下,并无基准类别的先验知识可以利用,因而, 只能先假定初始的参量,并通过预分类处理来形成集群,再 由集群的统计参数来调整预制的参量,接着再聚类,再调整。 如此不断地迭代,直到有关参数达到允许的范围为止。
所以说,非监督算法的核心问题是初始类别参数的选定,以 及它的迭代调整问题。
2.欧氏(Euclidean)距离
该距离是平面上两点之间的halanobis)距离
马氏距离是一种加权的欧氏距离,它通 过协方差矩阵来考虑变量的相关性。这 是由于在实际中,各点群的形状是大小 和方向各不相同的椭球体,如图所示, 尽管K点距MA的距离DA比距MB的距离 DB小,即DA<DB ,但由于B点群比A点 群离散得多,因而把K点划入B类更合 理。加权可以这样理解,计算的距离与 各点群的方差有关。方差愈大,计算的 距离就愈短。如果各个点群具有相同的 方差,则马氏距离是欧氏距离的平方。
• 然后,选定抽样集中任一像素作为第一个类别; • 给定一个光谱相似性比较阈值; • 依次把抽样集中每个像素的光谱特征与已建起的初始类别比
较,若该像素与其中一个初始类别相似,则作为该类集群中 的一个成员,若不与任何一个相似,则以该像素建立一个新 的初始类别。 • 当以上聚类过程完成后,每个已经建立起来的初始类别都包 含了一定成员,依据此可统计其集群中心(数学期望Mi)和 协方差矩阵Σi 。
4
8.1 概述
8.1.1 基本原理 同类地物在相同的条件下(光照、地形等)应该具有相同或 相似的光谱信息和空间信息特征。不同类的地物之间具有 差异根据这种差异,将图像中的所有像素按其性质分为若 干个类别(Class)的过程,称为图像的分类。 遥感图像分类以每个像素的光谱数据为基础进行。 假设遥感图像有K个波段,则(i,j)位置的像素在每个波段上 的灰度值可以构成表示为X=(x1,x2,…,xk)T,包含X的K 维空间称为特征空间,这样K个波段的多光谱图像便可以 用K维特征空间中的一系列点来表示。在遥感图像分类问 题中,常把图像中的某一类地物称为模式,而把属于该类 中的像素称为样本, X=(x1,x2,…,xk)T 可以称为样本的 观测值。
10
8.2 相似性度量
3.马氏(Mahalanobis)距离
11
8.2 相似性度量
4.相似系数 又称为余弦距离,其数学表达式为:
α为图像像素光谱与参照光谱之间的夹角(光谱角) X为图像像素光谱曲线向量 Y为参照光谱曲线向量。 当cosα的值接近1的时候有最好的估计光谱值和类别分类结果。
光谱角值以很小的α弧度角来表示,它代表了光谱曲线之间的 相似性,其变化范围是[0,π/2]
由于整幅遥感图像总体直方图的分布是由各类别直方图叠 加而形成的。同时,每个类别的集群中心一般位于本类别直方 图的峰值位置,而在该位置处的总体直方图往往会出现局部峰 值。因而,该法便以搜索总体直方图局部峰值为原理来选定初 始类别的中心。其基本过程包括:图像数据抽样集的获取,建 立总体直方图以及搜索直方图局部峰值3个步骤。
8.3 工作流程
图像的预处理 选择分类方法
1.确定工作范围 4.辐射校正
2.多源图像的几何配准 5.几何精纠正
3.噪声处理 6.多图像融合
特征选择和提取
选择合适的分类 参数进行分类
特征选择是从众多特征中挑选出可以参加分类运算的若干个特征。
特征提取是在特征选择以后,利用特征提取算法(如主成分分析算法) 从原始特征中求出最能反映地物类别性质的一组新特征。
16
8.4 非监督分类
主要过程如下: (1)确定初始类别参数,即确定最初类别数和类别中心 (集群中心) (2)计算每一个像元所对应的特征矢量与各集群中心的 距离。 (3)选与中心距离最短的类别作为这一矢量的所属类别。 (4)计算新的类别均值向量。 (5)比较新的类别均值与原中心位置上的变化,若位置 发生了改变,则以新的类别均值作为聚类中心,再从第2 步开始重复,进行反复迭代操作。 (6)如果聚类中心不再变化,计算停止。
式中:m为已形成的初始类别数。
4)在此基础上,再对所有各剩余点的最小距离进行相互比较,
取其中最大者,并选择与该最大的最小距离相应的抽样点(例
如
X 11
)作为一个新的初始类别中心点(例如
)。
Z3 X11
重复以上步骤,直到初始类别的个数达到需要的个数为止。
25
8.4 非监督分类
8.4.1 初始类别参数的选定 4.局部直方图峰值定心法
22
8.4 非监督分类
8.4.1 初始类别参数的选定
3.最大最小距离选心法
该法的选心原则是使各初始类别之间,尽可能地保持远距
离,为做到这一点,首先在整幅图像中按一定方式(如等间隔
方式)获取一个抽样点像素集合
{
X
}
{X}{X1X2X3 Xn}
式中n为抽样个数
23
8.4 非监督分类
Z q im ii[2 (q 1 )/Q ( 1 ) 1 ]
例如, Q=4时
Z 1i
mi
;
i
Z 3i
m
+
i
1 3
;
i
Z 2i
mi
1 3
i;
Z 4i
m
+
i
;
i
用作图法不难证明,该法所选定的初始类别中心通过整体集群
中心的一条直线均匀散布的,并基本上包括在整体集群的范围 之内
运用距离判别函数时,要求各个类别点群的中心位置己知。 对于光谱特征空间中的任一点k,计算它到各类中心点的距离 d。若di<dj,则k像素属于i类而不属于j类,此处,判别准则 为di<dj
8
8.2 相似性度量
1.绝对距离
dik为当前像素i到类k的距离, P为波段数, xij为像素i在j波段的像素值, Mjk为类k在波段j的均值。
特征选取的原则:尽可能选取数量适当有助于增加同类地物的类聚性 (即类内个体间离散性小)、能扩大地物类别之间差异性的特征变量。
分类后处理
成果输出
13
8.4 非监督分类 根据事先指定的某一准则,让计算机自动进行判
别归类,无需人为干预,分类后需确定地面类别。在非 监督分类中,先确定光谱可分的类别(光谱类)然后定 义它们的信息类。
17
8.4 非监督分类
8.4.1 初始类别参数的选定 8.4.2 常用非监督分类算法
18
8.4 非监督分类
8.4.1 初始类别参数的选定 初始类别参数是指基准类别集群中心(数学期望Mi),以 及集群分布的协方差矩阵Σi 。因为无论采用何种判别函数, 都要预先确定其初始类别的参量。以下为几种常用的方法。
7
8.2 相似性度量
距离是常用的相似性度量的概念。分类是确定像素距离哪个 点群中心较近,或落入哪个点群范围可能性大的问题。像素 与点群的距离越近,那么,属于该点群的可能性越高。按照 一定的准则,当距离小于一定值时,像素被划分给最近的点 群。每个点群为一个类。
同一类别中点间的距离一般来说比不同类别点间距离要小。 也可以认为,一个点属于某一类,一定与这个类别中心的距 离比与其他类别中心的距离小。因此,在点群(团)中心已知的 情况下,以每个点与点群中心的距离作为判定的准则,就可 以完成分类工作。
遥感数字图像处理
林金堂 闽江学院地理科学系
1
第8章 遥感图像分类
2
第8章 遥感图像分类
8.1 概述 8.2 相似性度量 8.3 工作流程 8.4 非监督分类 8.5 监督分类 8.6 其他分类方法 8.7 分类后处理 8.8 分类精度分析
3
8.1 概述
8.1.1 基本原理 8.1.2 分类方法
5
8.1 概述
8.1.2 分类方法 根据是否需要分类人员事先提供已知类别及其训练样本, 对分类器进行训练和监督,可将遥感图像分类方法划分为 监督分类和非监督分类。 事先己经知道类别的部分信息(即类别的先验知识),对未 知类别的样本进行分类的方法称之为监督分类(Supervised Classification)。事先没有类别的先验知识,对未知类别的 样本进行分类的方法称之为非监督分类(Unsupervised Classification) 根据分类使用的统计数学方法可以分为随机统计方法和模 糊数学方法分类。前者以随机数学理论为基础,包括K-均 值分类、最大相似性分类等。后者以模糊数学理论为基础, 主要是模糊分类。
N
( xij
j 1
mi
)
2
1 N
mi
N
xij
j 1
式中,i为波段号;j为像素点号;x ij 为像素i在第j波段的亮度 值;n为波段数;N为像素总数。
21
8.4 非监督分类
8.4.1 初始类别参数的选定 2.总体直方图均匀选心法
现假设我们需要Q个初始类别,每个初始类别集群中心可 按下式确定:
则
Z2 X 7
24
8.4 非监督分类
8.4.1 初始类别参数的选定 3.最大最小距离选心法
3)对于剩余的每个抽样点,计算它到已有各初始类别中心的距
离 Dij(i1,2,,m) m为已有初始类别数,并取其中的最小距离
作为该点的代表距离 D j
D j mD i1jn ,D 2 { j, ,D m }j
• 1.像素光谱特征的比较法 • 2.总体直方图均匀选心法 • 3.最大最小距离选心法 • 4.局部直方图峰值定心法
19
8.4 非监督分类
8.4.1 初始类别参数的选定 1.像素光谱特征的比较法
• 首先,在遥感图像中定义一个抽样集,它可以是整幅图像的 所有像素,但通常是按一定间隔抽样的像素。
在进行聚类分析时,首先要确定基准类别的参量。在非监督 分类的情况下,并无基准类别的先验知识可以利用,因而, 只能先假定初始的参量,并通过预分类处理来形成集群,再 由集群的统计参数来调整预制的参量,接着再聚类,再调整。 如此不断地迭代,直到有关参数达到允许的范围为止。
所以说,非监督算法的核心问题是初始类别参数的选定,以 及它的迭代调整问题。
2.欧氏(Euclidean)距离
该距离是平面上两点之间的halanobis)距离
马氏距离是一种加权的欧氏距离,它通 过协方差矩阵来考虑变量的相关性。这 是由于在实际中,各点群的形状是大小 和方向各不相同的椭球体,如图所示, 尽管K点距MA的距离DA比距MB的距离 DB小,即DA<DB ,但由于B点群比A点 群离散得多,因而把K点划入B类更合 理。加权可以这样理解,计算的距离与 各点群的方差有关。方差愈大,计算的 距离就愈短。如果各个点群具有相同的 方差,则马氏距离是欧氏距离的平方。
• 然后,选定抽样集中任一像素作为第一个类别; • 给定一个光谱相似性比较阈值; • 依次把抽样集中每个像素的光谱特征与已建起的初始类别比
较,若该像素与其中一个初始类别相似,则作为该类集群中 的一个成员,若不与任何一个相似,则以该像素建立一个新 的初始类别。 • 当以上聚类过程完成后,每个已经建立起来的初始类别都包 含了一定成员,依据此可统计其集群中心(数学期望Mi)和 协方差矩阵Σi 。
4
8.1 概述
8.1.1 基本原理 同类地物在相同的条件下(光照、地形等)应该具有相同或 相似的光谱信息和空间信息特征。不同类的地物之间具有 差异根据这种差异,将图像中的所有像素按其性质分为若 干个类别(Class)的过程,称为图像的分类。 遥感图像分类以每个像素的光谱数据为基础进行。 假设遥感图像有K个波段,则(i,j)位置的像素在每个波段上 的灰度值可以构成表示为X=(x1,x2,…,xk)T,包含X的K 维空间称为特征空间,这样K个波段的多光谱图像便可以 用K维特征空间中的一系列点来表示。在遥感图像分类问 题中,常把图像中的某一类地物称为模式,而把属于该类 中的像素称为样本, X=(x1,x2,…,xk)T 可以称为样本的 观测值。
10
8.2 相似性度量
3.马氏(Mahalanobis)距离
11
8.2 相似性度量
4.相似系数 又称为余弦距离,其数学表达式为:
α为图像像素光谱与参照光谱之间的夹角(光谱角) X为图像像素光谱曲线向量 Y为参照光谱曲线向量。 当cosα的值接近1的时候有最好的估计光谱值和类别分类结果。
光谱角值以很小的α弧度角来表示,它代表了光谱曲线之间的 相似性,其变化范围是[0,π/2]
由于整幅遥感图像总体直方图的分布是由各类别直方图叠 加而形成的。同时,每个类别的集群中心一般位于本类别直方 图的峰值位置,而在该位置处的总体直方图往往会出现局部峰 值。因而,该法便以搜索总体直方图局部峰值为原理来选定初 始类别的中心。其基本过程包括:图像数据抽样集的获取,建 立总体直方图以及搜索直方图局部峰值3个步骤。