遥感数字图像处理-第五章 遥感数字图像的计算机分类(一)2016

1、1,第五章 遥感数字图像计算机分类,一、遥感数字图像计算机分类的一般原理 二、遥感分类的常用判别函数 三、遥感数字图像的分类方法 四、遥感数字图像的分类后处理 五、光谱特征分类中的辅助处理技术 六、遥感图像计算机分类新方法,2,一、遥感数字图像计算机分类的一般原理,遥感图像解译,3,一、遥感数字图像计算机分类的一般原理,遥感图像解译,遥感图像数据,专题图像 专题地图,专业应用,解译,方法,计算机分类,目视解译,光谱规律 地学规律 解译者的经验,解译标志（亮度、色调、位置、时间、纹理、结构等）,景物类型,属性识别、分类,提取信息、识别地物,4,图像分类的目的 将图像中每个像元根据其不同波段的光谱

2、亮度、空间结构特征或其他信息，按照某种规则或算法划分为不同的类别。 遥感图像分类 利用计算机技术来模拟人类的识别功能,对地球表面及其环境在遥感图像上的信息进行属性的自动判别和分类，达到提取所需地物信息的目的。,5,计算机分类 是通过模式识别理论，利用计算机将遥感图象自动分成若干地物类别的方法。 如土地覆盖/土地利用分类、森林类型分类、植被 类型分类、岩性分类、 数据信息 （遥感数据-地物信息）,主要数据类型： 1、原始光谱数据； 2、光谱变换后数据； 3、非遥感数据,6,理论依据 在理想条件下，图像中同类地物在相同的条件下应具有相同或相似的光谱信息特征和空间信息特征，从而表现出同类地物的某种内

3、在相似性，即同类地物像元的特征向量将集群在同一特征空间区域，而不同地区的光谱信息特征或空间信息特征应不同，因而将集群在不同的特征空间区域。,7,遥感图像分类的实质,图象分类过程的总目标是，将图象中所有像元自动地进行土地覆盖类型或土地覆盖专题的分类。,8,影响遥感图像分类精度的因素 大气状况的影响：吸收、散射。 下垫面的影响：下垫面的覆盖类型和起伏状态对分类具有一定的影响。 其他因素的影响：云朵覆盖；不同时相的光照条件不同，同一地物的电磁辐射能量不同；地物边界的多样性。,9,统计模式识别的概念和基本问题,模式（pattern）: 在多波段图象中，每个象元都具有一组对应取值，称为象元模式，即一个像

4、元对应多个值（地理信息的多维性）。 特征（feature）：在多波段图象中，每个波段都可看作一个变量，称为特征变量。特征变量构成特征空间。 波段：光谱波段 其它派生波段（纹理、上下文关系、波段比等） 辅助数据(ancillary data) (非遥感数据，如DEM、 土壤类型) 特征提取（feature extraction）：通过变换找出最能反映地物类别差异的特征变量用于分类的过程,10,模式识别 是人们沿着仿生学的道路，用计算机系统作为工具来模拟人类的感知和识别智能，它是人工智能的一个分支。,xn,分类器（或称判决器），可以根据一定的分类规则，把某一测量矢量X划入某一组预先规定的类别之中去

5、。,自然模式,接收器 (传感器),分类器 (判决器),x2,x1,结果,.,模式识别系统的模型,11,遥感图像目视判读与计算机分类的比较 相同点：目的一致； 不同点：目视判读直接利用人类的自然识别智能； 计算机分类利用计算机模拟人类的识别能力。,遥感图像的计算机分类 它是对地球表面及其环境在遥感图像上的信息进行属性的识别和分类，从而达到识别图像信息所相应的实际地物，提取所需地物信息。,12,遥感图像的模式特征及其分类方法 光谱特征 纹理特征 基于光谱特征的统计识别方法 基于纹理特征的统计识别方法 （主要） （次要、补充）,13,统计模式识别的基本含义 （1）把识别对象的特征的每一个观测量视为从

6、属于一定分布规律的随机变量； （2）在多维观测的情况下，把识别对象特征的各维观测值的总体视为一个随机矢量，每一个随机矢量在一个多维特征空间中都有一个特征点与之相对应； （3）所有特征点的全体在特征空间中将形成一系列的分布群体，每个分布群体中的特征点被认为具有相似特征，并可以划为同一类别； （4） 找到各个分布群体的边界线（面）或确定任意特征点落入每个分布群体中的条件概率，并以此依据来确定特征点（或相应的识别对象）的分类。,14,光谱特征向量与特征点集群,基于统计模式识别的分类方法 监督分类与非监督分类 分类处理与增强处理的异同 共同点：增强和提取遥感图像中目标物的信息 不同点： 增强处理增强视

7、觉效果提高图像的可解译性（定性） 图像分类地物类别的区分（定量信息）,15,光谱特征向量与特征点集群,图像波段总数,光谱特征向量 光谱特征变量：由于受外界各种因素的影响，使得同类地物的成像亮度值总是带有随机误差，导致图像亮度值（即光谱特征）的观测值为一个随机变量（x）。 光谱特征向量：同名地物点在不同波段图像中亮度的观测量构成一个多维的随机向量（X），称为光谱特征向量。即,16,光谱特征向量与特征点集群,光谱特征空间 为度量地物的光谱特征，建立的以各波段图像的亮度分布为子空间的多维空间。,17,光谱特征向量与特征点集群,特征点集群 每个地物点依其在各个波段所具有的光谱值可以在一个多维空间中找到

8、一个相应的特征点，但由于随机性所致，同类地物的各取样点在光谱特征空间中的特征点将不可能表现为一个点，而是表现为一个相对聚集的点集群，不同类地物点集群在特征空间中一般相互分离，这些点集群称为特征点集群。,18,特征点集群在空间中的分布情况,最理想情况 不同类别集群至少在一个特征子空间（某一波段图像）中的投影（即亮度范围）是完全可以相互区分开的可以用简单的图像密度分割实现。,19,典型情况 不同类别地物的集群，在任一个子空间都有重叠现象存在，但在总的特征空间中却是可以完全分开的。即单波段的图像不能实现图像的分类，只有利用多波段图像在多维空间中才能实现精确分类。,20,一般情况 不论是在总的特征空间

9、还是任一子空间，不同类别的集群之间总有重叠现象，这时重叠部分的特征点所相应的地物，在分类时会出现不同程度的分类误差，这是遥感图像分类中最常见的现象。,21,地物与光谱特征空间的关系,22,计算机分类的基本原理,基本原理 不同的地物具有不同的光谱特征，同类地物具有相同或相似的光谱特征将图像中的每个像元或区域划归为若干类别是的一种，即通过光谱特征分析选择特征参数将特征空间划分为不重叠的子空间将影像像元划分到各个子空间实现分类 图象分类：基于数字图象中反映的同类地物的光谱相似性和异类地物的光谱差异性。 依据是遥感图像像素的相似度。,距离 相关系数,23,计算机分类的基本原理,基本思想 同类地物的光谱

10、特征比较接近，在特征空间中的点聚集在该类的中心附近；多类目标聚集多个点族。,分类算法的核心 判别函数 判别准则,方法 判别函数： 找到一个函数使fAB(X)=0 判别准则： fAB(X)0，X为A类； fAB(X) n+1 ,甚至可以达到3n 。 B、代表性 训练样本应该反映该类地物的亮度特征；且需注意同类地物空间分布的不连续性（如，可以分别在块草地中选取训练样本）。 、亮度分布 从亮度特征角度而言，对于同一类地物具有不同亮度特征情况，通常都要选取（同物异谱）。,40,、纯度的把握 一般来说，在已知类别的中间部位纯度较高；如果过纯，则方差过小，导致许多此类地物的像元被拒分；如果不纯，则方差过大

11、，导致不是此类地物的像元被分到该类别中，导致出现“错分”的情况； 把握纯度的方法：读像元的亮度值，确定该类地物的亮度分布区间；同时，注意其它地物亮度值与其近似部分的像元亮度分布区间，依“错分损失最小”的原则，确定“分类界限”。,41,训练区的选择,42,选择训练样区-直接在图像上选取,43,选择训练样区在特征空间上选取,44,训练样本中的各个像素在特征空间中的分布情况,45,最小距离分类法 最小距离分类法是以特征空间中的距离作为像元分类的依据。原理简单，计算速度快，分类精度不高。,Step 2 for each unclassified pixel, calculate the distanc

12、e to average for each training area,特点：原理简单，容易理解，计算速度快。 缺点：因只考虑每一类样本的均值，而不用管类别内部的方差（每一类样本的分布），也不用考虑类别之间的协方差（类别和类别之间的相关关系），所以分类精度不高。 主要用途：一般不用它作为分类对精度有高要求的分类，但它可以在快速浏览分类概况中使用。,改进方法,多级切割分类法,通过设定在各轴上的一系列分割点，将多维特征空间划分成分别对应不同分类类别的互不重叠的特征子空间的分类方法。 对于一个未知类别的像素来说，它的分类取决于它落入哪个类别特征子空间中。 多级切割分类法要求训练区样本的选择必须覆盖所

13、有的类型。 运用多级分割法分类前，需要先进行主成分分析，或采用其他方法对各轴进行相互独立的正交变换。,用多级切割法分割三维特征空间,50,平行多面体分类,要求分割面总是与各特征轴正交，如果各类别在特征空间中呈倾斜分布，就会产生分类误差。可先采用各轴的正交变换后，再进行多级分割。,最大似然比分类法(Maximum Likelihood),通过求出每个像素对于各类别的归属概率，把该像素分到归属概率最大的类别中去的方法。 假定训练区地物的光谱特征和自然界大部分随机现象一样，近似服从正态分布。 利用训练区可求出均值、方差以及协方差等特征参数，从而可求出总体的先验概率密度函数。 然后，每一个像素归为某一

14、类别的归属概率可求出。,54,使用最大似然比分类法的注意事项 各类别的训练数据至少要为特征维数的2-3倍以上，这样才能统计出具有较高精度的均值、方差与协方差； 如果2个以上的波段相关性很强，那么方差协方差矩阵的逆矩阵就不存在，或非常不稳定。在训练数据几乎都取相同值的均质性数据组的情况下也是如此。此时，最好采用主成分分析法，把维数减到仅剩相互独立的波段。 当总体分布不符合正态分布时，不适于采用以正态分布的假设为基础的最大似然比分类法。其分类精度也将下降。,55,分类结果,56,非监督分类法不必使用训练样区数据作为分类的基础。通过检验大量的未知像元并根据多光谱空间中的自然集群将它们分为若干类。 基

15、本前提：同种地物类型的值在测量空间上应该相互接近，而不同的地物类型应该具有比较明显的可分性，即：假定遥感影像上同类物体在同样条件下具有相同的光谱信息特征。 主要方法：采用聚类分析方法，聚类是把一组像素按照相似性归成若干类别，即“物以类聚”。,非监督分类(unsupervised classification),57,主要目的是使得同一类别的像素之间的距离尽可能的小而不同类别上的像素间的距离尽可能的大。 分类标准的确定不需要人的参与，由计算机根据自动按照某一标准（例如距离最短）自动进行。需要确定要分几种类别，或者类似的输入条件。分类后的结果，还需要再给出具体的涵义，即先分类后确定属性。 区别：监

16、督分类法确定有用信息类别，然后再检验它们的光谱可分性；而非监督分类法先确定光谱的可分性，然后再确定它们的有用信息。,58,分类算法的思路 先选择若干个模式点作为聚类的中心，每一个中心代表一个类别，按照某种相似性度量方法将各模式归于各聚类中心所代表的类别，形成初始分类。然后由聚类准则判断初始分类是否合理，如果不合理就修改分类，如此反复迭代运算，直到合理为止。,59,根据事先指定的某一准则，而进行计算机自动判别归类，无须人为干预，分类后需确定地面类别。 在非监督分类中，先确定光谱可分的类别，然后定义它们的信息类。,60,非监督分类的核心：初始类别参数的选定（基准类别集群中心（数学期望），集群分布的

17、协方差矩）；迭代调整。 非监督分类的主要过程,61,非监督分类基本过程,62,初始类别参数的选定方法 像素光谱特征的比较法,63,总体直方图均匀定心法在总体直方图的基础上选定类别中心,64,65,最大最小距离选心法（效果最好） 选心原则：各初始类别之间尽可能保持远距离,66,先抽样，设有n个样本，计算步骤：,67,非监督分类的主要算法 ISODATA法（迭代自组织数据分析技术） 基本思想：在初始类别参数选定的基础上，通过逐步趋近求出实际的基准类别参数。,68,算法流程：,69,70,71,72,总结：以初始类别为“种子”进行自动迭代聚类过程，它可以自动地进行类别的“合并”与“分裂”，其各个参数

18、也在不断地聚类调整中逐渐确定，并最终构成所需要的判别函数。 因此，基准类别参数的确定过程，也正是利用光谱特征本身的统计性质对判别函数的不断调整和“训练”的过程。,73,K-mean算法 基本思想：通过迭代，移动各个基准类别的中心，得到最好的聚类结果。 聚类准则（收敛条件）：使每一聚类中，多模式点到该类别的中心的距离的平方和最小。即：对于图像中互不相交的任意一类，计算该类中的像素值与该类均值差的平方和。将图像中所有类的差的平方和相加，并使相加后的值达到最小！,74,75,K-mean算法的计算步骤：,76,K-mean算法与ISODATA算法的区别： （1）ISODATA算法它不是每调整一个样本

19、的类别就重新计算一次各样本的均值，而是在把所有的样本都调整完之后才重新计算。 （2） K-mean算法是逐个样本修正法， ISODATA算法是成批样本修正法； （3） ISODATA算法不仅可以通过调整样本所属类别完成样本的聚类分析，而且可以自动地进行类别“合并”与“分裂”，从而得到类数比较合理的聚类结果。,77,4-3-2假彩色合成图像,非监督分类(聚类)结果: 15个光谱类,非监督分类,78,非监督分类 1. 4-3-2假彩色合成图象(香港九龙); 2. 聚类结果(10类),79,3. 聚类结果合并(5类); 4. 最终结果 (类别颜色改变),监督分类与非监督分类的比较,根本区别：是否利用训练场地来获取先验的类别知识。 监督分类的关键是选择训练场地。 监督分类法优