遥感技术基础07.ppt

遥感技术基础遥感技术基础学年学期：学年学期：2010-2011学年第一学期学年第一学期学学 院：测绘与地理科学学院院：测绘与地理科学学院班班 级：地理级：地理081-3班班教教 师：李师：李 玉玉2022/12/301课程内容课程内容一、分割（一、分割（SegmentationSegmentation）的一般原理的一般原理二、非监督分割二、非监督分割 （UnsupervisedUnsupervised）三、监督分割（三、监督分割（SupervisedSupervised）第七讲第七讲 遥感图像分割遥感图像分割2022/12/302一、一、分割（分割（SegmentationSegmentation）的一般原理）的一般原理n计算机判读：以计算机系统为支撑环境，利用模式识计算机判读：以计算机系统为支撑环境，利用模式识别技术与人工智能技术相结合，根据目标地物的影像别技术与人工智能技术相结合，根据目标地物的影像特征（颜色、形状、纹理、空间位置），结合专家知特征（颜色、形状、纹理、空间位置），结合专家知识库中目标地物的解译经验和成像规律等知识进行分识库中目标地物的解译经验和成像规律等知识进行分析和推理析和推理,实现对遥感图象的理解实现对遥感图象的理解,完成对遥感图象的完成对遥感图象的解译。解译。第七讲第七讲 遥感图像分割遥感图像分割2022/12/303n模式模式(Pattern)(Pattern)：即某种实物的标准形式或使人可以照着：即某种实物的标准形式或使人可以照着做的标准样式。举例：某地物的光谱特征曲线反映了该地做的标准样式。举例：某地物的光谱特征曲线反映了该地物的反射特性，所以，该光谱特征曲线就是该地物的一个物的反射特性，所以，该光谱特征曲线就是该地物的一个模式。模式。n模式识别模式识别 (Pattern Recognition)(Pattern Recognition)：即对需识别或分类的：即对需识别或分类的对象，进行一系列测量，例如，对未知类别的地物，测量对象，进行一系列测量，例如，对未知类别的地物，测量它在等波长处的反射特性。然后将这一系列测量所构成的它在等波长处的反射特性。然后将这一系列测量所构成的模式与已知类别的地物模式作比较，看它与那一个相同或模式与已知类别的地物模式作比较，看它与那一个相同或相似，即判认它是属于哪一个类别的地物。相似，即判认它是属于哪一个类别的地物。第七讲第七讲 遥感图像分割遥感图像分割2022/12/304n计算机判读与目视判读目的一样，但在判读和识别地计算机判读与目视判读目的一样，但在判读和识别地物的模式方面，有着明显的不同。物的模式方面，有着明显的不同。p目视判读：主要以图像的空间特征（地物的几何特征和光谱目视判读：主要以图像的空间特征（地物的几何特征和光谱特征的空间反映）为判读依据。特征的空间反映）为判读依据。p计算机判读：主要以图像像元的灰度（地物光谱特征的直接计算机判读：主要以图像像元的灰度（地物光谱特征的直接反映）为判读依据反映）为判读依据。第七讲第七讲 遥感图像分割遥感图像分割2022/12/3052022/12/306n图像分割图像分割:所谓图像分割是指将图像中具有特殊涵义的不同所谓图像分割是指将图像中具有特殊涵义的不同区域区分开来，这些区域是互相不交叉的，每一个区域都满区域区分开来，这些区域是互相不交叉的，每一个区域都满足特定区域的一致性。足特定区域的一致性。第七讲第七讲 遥感图像分割遥感图像分割2022/12/307n图像分割模型图像分割模型:对图像对图像 X=x1,xn 的分割就是将其划分为的分割就是将其划分为满足如下条件的子区域满足如下条件的子区域 S1,Sm,其中其中m m为子区域的个数。为子区域的个数。（1）Sj X,Sj ,j=1,m;（2）Sj Sk=,j k;（3）j=1,mSj=X。n图像分割的基本问题图像分割的基本问题p如何决定子区域个数；如何决定子区域个数；p特殊涵义的定义与建模；特殊涵义的定义与建模；p实现分割的方法。实现分割的方法。第七讲第七讲 遥感图像分割遥感图像分割2022/12/308n图像分割的理论基础图像分割的理论基础p同类地物具有相同（似）的光谱特征，不同地物的光谱特征具有明同类地物具有相同（似）的光谱特征，不同地物的光谱特征具有明显的差别。由于影响地物光谱特征的因素很多，所以影像的判读分显的差别。由于影响地物光谱特征的因素很多，所以影像的判读分类都是建立在统计分析的基础上的类都是建立在统计分析的基础上的。p同类地物的图像灰度概率在单波段（一维空间）符合正态分布规律。同类地物的图像灰度概率在单波段（一维空间）符合正态分布规律。p多维图像（即多波段）中的一个像元值（灰度）向量。在几何上相多维图像（即多波段）中的一个像元值（灰度）向量。在几何上相当于多维空间中的一个点，而同类地物的像元值，既不集中于一点，当于多维空间中的一个点，而同类地物的像元值，既不集中于一点，也绝非是杂乱无章的分布，而是相对地密集在一起，形成一个点群也绝非是杂乱无章的分布，而是相对地密集在一起，形成一个点群（一个点群就是地物的一种类别）。一般情况下，点群的边界不是（一个点群就是地物的一种类别）。一般情况下，点群的边界不是截然的，有少部分重叠和交错的情况截然的，有少部分重叠和交错的情况。第七讲第七讲 遥感图像分割遥感图像分割2022/12/309n图像分割的方法图像分割的方法p监督（监督（SupervisedSupervised）分割法（训练场地法）分割法（训练场地法/先学习后分割法）：先学习后分割法）：即先选择有代表性的实验区（训练区），用已知地面的各种地即先选择有代表性的实验区（训练区），用已知地面的各种地物光谱特征来训练计算机，取得识别判别规则，并以此做标准物光谱特征来训练计算机，取得识别判别规则，并以此做标准对未知地区的遥感数据进行自动分割识别。对未知地区的遥感数据进行自动分割识别。p非监督（非监督（UnsupervisedUnsupervised）分类法：即按照灰度值向量或波谱样）分类法：即按照灰度值向量或波谱样式在特征空间聚集的情况划分点群或类别。其类属是通过对各式在特征空间聚集的情况划分点群或类别。其类属是通过对各类光谱响应曲线进行分析以及与实地调查数据相比较后确定的。类光谱响应曲线进行分析以及与实地调查数据相比较后确定的。第七讲第七讲 遥感图像分割遥感图像分割2022/12/3010二、非监督分割二、非监督分割n以聚类分析的方法为例以聚类分析的方法为例n主要过程主要过程(1)(1)确定最初类别数和类别中心确定最初类别数和类别中心;(2)(2)计算每个像元多对应的特征矢量与各聚类中心的距离计算每个像元多对应的特征矢量与各聚类中心的距离;(3)(3)选与其中心距离最近的类别作为这一矢量选与其中心距离最近的类别作为这一矢量(像元像元)的所属类别的所属类别;(4)(4)计算新的类别均值向量计算新的类别均值向量;(5)(5)比较新的类别均值与原中心位置的变化比较新的类别均值与原中心位置的变化,形成新的聚类中心形成新的聚类中心;(6)(6)重复重复II-IV,II-IV,反复迭代反复迭代;(7)(7)如聚类中心不再变化如聚类中心不再变化,停止计算停止计算.第七讲第七讲 遥感图像分割遥感图像分割2022/12/3011n常用的非监督分割方法常用的非监督分割方法pK-MEANSK-MEANS方法方法pISODATAISODATA（Iterative Self-Organizing Data Analysis Technique，迭代自，迭代自组织数据分析技术组织数据分析技术）方法）方法nK-K-均值分割均值分割p在数据空间内平均计算各类的均值，使用最小距离法，将未在数据空间内平均计算各类的均值，使用最小距离法，将未知像元划分到某一类中，再调整均值，直到迭代次数满足，知像元划分到某一类中，再调整均值，直到迭代次数满足，或者距离阈值满足。或者距离阈值满足。第七讲第七讲 遥感图像分割遥感图像分割2022/12/3012pK-K-均值分割的步骤均值分割的步骤1.1.任意选择任意选择K K个聚类中心，一般选前个聚类中心，一般选前K K个样本；个样本；2.2.迭代，未知样本迭代，未知样本X X分到距离最近的类中；分到距离最近的类中；3.3.根据根据Step2Step2的结果，重新计算聚类中心；的结果，重新计算聚类中心；4.4.每一类的像元数变化少于选择的像元变化阈值或已经到了迭代的最多每一类的像元数变化少于选择的像元变化阈值或已经到了迭代的最多次数，算法结束。次数，算法结束。pK-K-均值分割的输入参数均值分割的输入参数 类别数目类别数目:要分成几类；要分成几类；变化阈值变化阈值Change threshold%(0-100):Change threshold%(0-100):如果每一类的像元数目变化小如果每一类的像元数目变化小于此数值，则迭代停止。于此数值，则迭代停止。第七讲第七讲 遥感图像分割遥感图像分割2022/12/3013nISODATA分割分割p在在K-K-均值算法的基础上，加入了试探性的步骤，能够吸取中间均值算法的基础上，加入了试探性的步骤，能够吸取中间结果的经验，在迭代的过程中可以进行类别的分离和合并，具结果的经验，在迭代的过程中可以进行类别的分离和合并，具有有“自组织自组织”性。性。第七讲第七讲 遥感图像分割遥感图像分割2022/12/3014pISODATA分割的步骤分割的步骤1.1.在整幅图像的特征空间上，随机地选择初始的集群组（即类别）在整幅图像的特征空间上，随机地选择初始的集群组（即类别）中心中心2.2.计算其它像元离这些中心的距离，按照最小距离法划入其对应的计算其它像元离这些中心的距离，按照最小距离法划入其对应的类别中；类别中；3.3.计算集群像元总数，如数目太少，删除该类别；计算集群像元总数，如数目太少，删除该类别；4.4.重新计算每个类别的均值，按照定义的参数合并或分开类别；重新计算每个类别的均值，按照定义的参数合并或分开类别；5.5.重复以上过程，直到有关参量的变动在容忍的范围内（每一类的重复以上过程，直到有关参量的变动在容忍的范围内（每一类的像元数变化少于选择的像元变化阈值或已经到了迭代的最多次数）像元数变化少于选择的像元变化阈值或已经到了迭代的最多次数）为止。为止。第七讲第七讲 遥感图像分割遥感图像分割2022/12/3015pISODATAISODATA分割的分割的的输入参数的输入参数 类别数目类别数目;变化阈值变化阈值;最长的时间或迭代次数最长的时间或迭代次数;每个集群组中最小的像元数量每个集群组中最小的像元数量;集群分裂标准：每类集群所容许的最大标准差集群分裂标准：每类集群所容许的最大标准差;集群合并标准：两集群均值间允许的最小距离集群合并标准：两集群均值间允许的最小距离.第七讲第七讲 遥感图像分割遥感图像分割2022/12/3016n人工神经网络法人工神经网络法p人工神经网络是基于生物神经系统分布式存储、并行处理、自人工神经网络是基于生物神经系统分布式存储、并行处理、自适应学习这些现象构造出具有一些低级智慧的人工神经网络系适应学习这些现象构造出具有一些低级智慧的人工神经网络系统。统。p目前代表性的神经网络模型有：目前代表性的神经网络模型有：BPBP（Back PropagationBack Propagation）神经网络）神经网络模糊自组织神经网络模糊自组织神经网络RBFRBF（Radial Basis FunctionRadial Basis Function）神经网络）神经网络KohonenKohonen自组织神经网络自组织神经网络第七讲第七讲 遥感图像分割遥感图像分割2022/12/3017第七讲第七讲 遥感图像分割遥感图像分割2022/12/3018n决策树分类法决策树分类法p决策树是一树状结构决策树是一树状结构,依据规则把遥感数据集依据规则把遥感数据集 一级级往下细分以定义一级级往下细分以定义 决策树的各个分支。决策树的各个分支。第七讲第七讲 遥感图像分割遥感图像分割2022/12/3019p基本思想：从基本思想：从“原级原级”（根结点）开始，利用表达式，每一个决策（根结点）开始，利用表达式，每一个决策将影像中的像元分成两类，使用另一表达式，每个新类又能被分成将影像中的像元分成两类，使用另一表达式，每个新类又能被分成另外的两个新类，如此不断地通过选择不同的特征用于进一步地有另外的两个新类，如此不断地通过选择不同的特征用于进一步地有效细分类，直到所要求的效细分类，直到所要求的“终极终极”（叶结点）类别分出为止（叶结点）类别分出为止.p在决策树分类中经常采用的特征在决策树分类中经常采用的特征光谱值；光谱值；通过光谱值算出来的指标（如通过光谱值算出来的指标（如NDVINDVI）；）；光谱值的算术运算值（如，和、差、比值等）；光谱值的算术运算值（如，和、差、比值等）；主成分主成分.p由于决策树分类法中的运算几乎都是由比较大小而组成的，所以与由于决策树分类法中的运算几乎都是由比较大小而组成的，所以与采用复杂计算公式的最大似然比分类法等相比，可以用很短的时间进采用复杂计算公式的最大似然比分类法等相比，可以用很短的时间进行分类处理行分类处理第七讲第七讲 遥感图像分割遥感图像分割2022/12/3020n专家系统方法专家系统方法p专家系统也是人工智能的一个分支，它采用人工智能语言将某专家系统也是人工智能的一个分支，它采用人工智能语言将某一领域（如地学知识）的专家分析方法或经验，对地物的多种一领域（如地学知识）的专家分析方法或经验，对地物的多种属性进行分析、判断，从而确定各地物的归属。属性进行分析、判断，从而确定各地物的归属。p利用这样的系统就可以把判读专家的经验性知识综合起来进行利用这样的系统就可以把判读专家的经验性知识综合起来进行分类。分类。n其它方法：除了以上常用的方法之外，还有多时相数据其它方法：除了以上常用的方法之外，还有多时相数据分类方法、模糊分类法、亚像元分类法，基于纹理的分分类方法、模糊分类法、亚像元分类法，基于纹理的分类法等等。类法等等。第七讲第七讲 遥感图像分割遥感图像分割2022/12/3021三、监督分割三、监督分割n常用的监督分类方法常用的监督分类方法p平行算法平行算法p最小距离法最小距离法p最大似然法最大似然法p马氏距离分类马氏距离分类p波谱角分类法波谱角分类法第七讲第七讲 遥感图像分割遥感图像分割2022/12/3022n训练区选择训练区选择p规则规则训练区要具有典型性和代表性；训练区要具有典型性和代表性；样本要保证一定数量；样本要保证一定数量；选取训练区使用的参考图件与分类数据在时间上和空间上最好保持选取训练区使用的参考图件与分类数据在时间上和空间上最好保持一致或相近。一致或相近。p方法方法实地调查实地调查屏幕选择：参考实地调查结果、研究区其它相关图件资料、辅助遥屏幕选择：参考实地调查结果、研究区其它相关图件资料、辅助遥感资料等信息，进行目视解译。感资料等信息，进行目视解译。第七讲第七讲 遥感图像分割遥感图像分割2022/12/3023n最小距离法（最小距离法（Minimum Distance）p利用训练样本中各类别在各波段的均值，根据各像元离训练样本利用训练样本中各类别在各波段的均值，根据各像元离训练样本平均值平均值 距离的距离的 大小来大小来 决定其决定其 类别。类别。第七讲第七讲 遥感图像分割遥感图像分割2022/12/3024第七讲第七讲 遥感图像分割遥感图像分割2022/12/3025n平行管道分类（平行六面体）平行管道分类（平行六面体）p规定每个特征的范围规定每个特征的范围第七讲第七讲 遥感图像分割遥感图像分割2022/12/3026n最大似然法（最大似然法（Maximum Likelihood）p计算每个像元属于各个类别的似然度，该像元分到似然度最大的计算每个像元属于各个类别的似然度，该像元分到似然度最大的类别中。类别中。p设设 P P(X X|i i)为某一类别为某一类别 i i 的特征矢量分布函数，它是把模式的特征矢量分布函数，它是把模式X X分分类到为最大的类别中的分类方法。类到为最大的类别中的分类方法。式中式中P P(i i)表示类别表示类别 i i 模式的先验概率，模式的先验概率，P P(X X|i i)表示条件概率表示条件概率密度函密度函 数，数，P P(i i|X X)为在观测模式为在观测模式 X X 时，该模式属于类别时，该模式属于类别 i i 的的似然度。似然度。这一方法称为最大似然法，它是理论上误差最小的分类这一方法称为最大似然法，它是理论上误差最小的分类方法。方法。第七讲第七讲 遥感图像分割遥感图像分割2022/12/3027p(a)(a)所示是一个用某值将特征空间分为两个区域的例子，因此而产所示是一个用某值将特征空间分为两个区域的例子，因此而产生生 的错误分类概率如图（的错误分类概率如图（b b）中的阴影部分的面积表示。）中的阴影部分的面积表示。T T值对应值对应的错分概的错分概 率最小率最小 。第七讲第七讲 遥感图像分割遥感图像分割(a)(b)2022/12/3028第七讲第七讲 遥感图像分割遥感图像分割2022/12/3029n监督分类与非监督分类监督分类与非监督分类p监督分类的优点监督分类的优点可充分利用对分类地区的先验知识，预先确定分类的类；可充分利用对分类地区的先验知识，预先确定分类的类；可控制训练样本的选择，并可通过反复检验训练样本，以提高分类可控制训练样本的选择，并可通过反复检验训练样本，以提高分类精度（避免分类中的严重错误）；精度（避免分类中的严重错误）；可避免非监督分类中对光谱集群组的重新归类。可避免非监督分类中对光谱集群组的重新归类。n监督分类的缺点监督分类的缺点人为主观因素较强；人为主观因素较强；训练样本的选取和评估需花费较多的人力、时间；训练样本的选取和评估需花费较多的人力、时间；只能识别训练样本中所定义的类别，对于因训练者不知或因数量太只能识别训练样本中所定义的类别，对于因训练者不知或因数量太少未被定义的类别，监督分类不能识别，从而影响分结果。少未被定义的类别，监督分类不能识别，从而影响分结果。第七讲第七讲 遥感图像分割遥感图像分割2022/12/3030p非监督分类的优点非监督分类的优点无需对分类区域有广泛地了解，仅需一定的知识来解释分类出的集群无需对分类区域有广泛地了解，仅需一定的知识来解释分类出的集群组组人为误差的机会减少，需输入的初始参数较少人为误差的机会减少，需输入的初始参数较少独特的、覆盖量小的类别均能够被识别。独特的、覆盖量小的类别均能够被识别。p非监督分类的缺点非监督分类的缺点后续处理比较麻烦后续处理比较麻烦;不一定分出想要的结果不一定分出想要的结果;不同图像的分类结果难以比较。不同图像的分类结果难以比较。第七讲第七讲 遥感图像分割遥感图像分割