【精品】遥感数字影像处理：分类与特征提取82.ppt

遥感数字影像处理：分类与特征提取,钱乐祥广州大学地理科学学院地理信息科学系2005.2.22,分类与特征提取,概述分类原理与基本过程遥感图像分类方法遥感图像多种特征提取图像解译专家系统,概述,遥感图像计算机分类以遥感数字图像为研究对象，在计算机系统支持下，综合运用地学分析、遥感图像处理、地理信息系统、模式识别与人工智能技术，实现地学专题信息的智能化获取。其基本目标是将人工目视解译遥感图像发展为计算机支持下的遥感图像理解。,概述,利用计算机对遥感数字图像进行分类难度很大。遥感图像是从遥远的高空成像的，成像过程要受传感器、大气条件、太阳位置等多种因素的影响。影像中所提供的目标地物信息不仅不完全，而且或多或少地带有噪声，因此人们需要从不完全的信息中尽可能精确地提取出地表场景中感兴趣的目标物。遥感影像信息量丰富，与一般的图像相比，其包容的内容远比普通的图像多，因而内容非常“拥挤”。不同地物间信息的相互影响与干扰使得要提取出感兴趣的目标变得非常困难。遥感图像的地域性、季节性和不同成像方式更增加了计算机对遥感数字图像进行解译的难度。,概述,由于利用遥感图像可以客观、真实和快速地获取地球表层信息，这些现势性很强的遥感数据在自然资源调查与评价、环境监测、自然灾害评估与军事侦察上具有广泛应用前景。因此，利用计算机进行遥感图像智能化解译，快速获取地表不同专题信息，并利用这些专题信息迅速地更新地理数据库，这是实现遥感图像自动理解的基础研究之一，也是地理信息系统中数据采集自动化研究的一个方向，因此具有重要的理论意义和应用前景。,分类原理与过程,计算机遥感图像分类是统计模式识别技术在遥感领域中的具体应用。统计模式识别的关键是提取待识别模式的一组统计特征值，然后按照一定准则作出决策，从而对数字图像予以识别。遥感图像分类的主要依据是地物的光谱特征，即地物电磁波辐射的多波段测量值，这些测量值可以用作遥感图像分类的原始特征变量。然而，就某些特定地物的分类而言，多波段影像的原始亮度值并不能很好地表达类别特征，因此需要对数字图像进行运算处理（如比值处理、差值处理、主成分变换以及KT变换等），以寻找能有效描述地物类别特征的模式变量，然后利用这些特征变量对数字图像进行分类。分类是对图像上每个像素按照亮度接近程度给出对应类别，以达到大致区分遥感图像中多种地物的目的。,分类原理与过程,分类过程中采用的统计特征变量包括：全局统计特征变量和局部统计特征变量。全局统计特征变量是将整个数字图像作为研究对象，从整个图像中获取或进行变换处理后获取变量，前者如地物的光谱特征，后者如对TM的6个波段数据进行KT变换（缨帽变换）获得的亮度特征，利用这两个变量就可以对遥感图像进行植被分类。局部统计特征变量是将数字图像分割成不同识别单元，在各个单元内分别抽取的统计特征变量。例如，纹理是在某一图像的部分区域中，以近乎周期性或周期性的种类、方式重复其自身局部基本模式的单元，因此可以利用矩阵作为特征对纹理进行识别。,分类原理与过程,取n*n的窗口，有关矩阵为：关联矩阵 以偏离图像灰度为i的点一定位置（方向和距离）的点的灰度为j的概率P(i,j)，求出关联矩阵，从该矩阵中算出各种纹理的特征量（能量、熵、相关性等）旋转矩阵 以方向上灰度为i的点连续出现k个的频率P(i,k)为元素，求出旋转矩阵。从矩阵中算出各种纹理的特征量。,分类原理与过程,在很多情况下，利用少量特征就可以进行遥感图像的地学专题分类，因此需要从遥感图像n个特征中选取k个特征作为分类依据，我们把从n个特征中选取k个更有效特征的过程称为特征提取。特征提取要求所选择的特征相对于其他特征更便于有效地分类，使图像分类不必在高维特征空间里进行，其变量的选择需要根据经验和反复的实验来确定。,分类原理与过程,统计特征变量可以构成特征空间，多波段遥感图像特征变量可以构成高维特征空间。一般说来，高维特征空间数据量大，但这些信息中仅包含少量的样本分类信息。为了抽取这些最有效的信息，可以通过变换把高维特征空间所表达的信息内容集中在一到几个变量图像上。主成分变换可以把互相存在相关性的原始多波段遥感图像转换为相互独立的多波段新图像，而且使原始遥感图像的绝大部分信息集中在变换后的前几个组分构成的图像上，实现特征空间降维和压缩的目的。,分类原理与过程,遥感图像计算机分类的依据是遥感图像像素的相似度。相似度是两类模式之间的相似程度。在遥感图像分类过程中，常使用距离和相关系数来衡量相似度。距离:特征空间中象元数据和分类类别特征的相似程度。距离最小即相似程度最大。度量特征空间中的距离经常采用以下几种算法:绝对值距离欧氏距离,x为像元数据矢量,类别k的平均值矢量,分类原理与过程,马氏距离（Mahalanobis，既考虑离散度，也考虑各轴间的总体分布相关)混合距离（像元i到第g类类均值的距离）,为g类k变量的均值,g类的像元数,分类原理与过程,相关系数 是指像素间的关联程度。采用相关系数衡量相似度时，相关程度越大，相似度越大。两个像素之间的相关系数rij可以定义为：,像元i的第k个分量,均值,分类原理与过程,分类方法包括监督分类和非监督分类。监督分类方法。首先需要从研究区域选取有代表性的训练区作为样本。根据已知训练区提供的样本，通过选择特征参数（如像素亮度均值、方差等），建立判别函数，据此对样本像元进行分类，依据样本类别的特征来识别非样本像元的归属类别。非监督分类方法。是在没有先验类别（训练区）作为样本的条件下，即事先不知道类别特征，主要根据像元间相似度的大小进行归类合并（将相似度大的像元归为一类）的方法。,分类原理与过程,遥感数字图像计算机分类基本过程如下：首先明确遥感图像分类的目的及其需要解决的问题，在此基础上根据应用目的选取特定区域的遥感数字图像，图像选取时应考虑图像的空间分辨率、光谱分辨率、成像时间、图像质量等。根据研究区域，收集与分析地面参考信息与有关数据。为提高计算机分类的精度，需要对数字图像进行辐射校正和几何纠正（这部分工作也可能由提供数字图像的卫星地面站完成）。对图像分类方法进行比较研究，掌握各种分类方法的优缺点，然后根据分类要求和图像数据的特征，选择合适的图像分类方法和算法。根据应用目的及图像数据的特征制定分类系统，确定分类类别，也可通过监督分类方法，从训练数据中提取图像数据特征，在分类过程中确定分类类别。,分类原理与过程,找出代表这些类别的统计特征。为了测定总体特征，在监督分类中可选择具有代表性的训练场地进行采样，测定其特征。在无监督分类中，可用聚类等方法对特征相似的像素进行归类，测定其特征。对遥感图像中各像素进行分类。包括对每个像素进行分类和对预先分割均匀的区域进行分类。分类精度检查。在监督分类中把已知的训练数据及分类类别与分类结果进行比较，确认分类的精度及可靠性。在非监督分类中，采用随机抽样方法，分类效果的好坏需经实际检验或利用分类区域的调查材料、专题图进行核查。对判别分析的结果统计检验。,分类的方法,利用遥感图像进行分类（classification）是以区别图像中所含的多个目标物为目的的，对每个像元或比较匀质的像元组给出对应其特征的名称。在分类中注重的是各像元的灰度及纹理等特征。监督分类：最小距离分类法、多级切割分类法、特征曲线窗口法、最大似然比分类法非监督分类：多级集群法、动态聚类法等,最小距离分类法,最小距离分类法（minimum distance classifier）是用特征空间中的距离表示像元数据和分类类别特征的相似程度，在距离最小时（相似度最大）的类别上对像元数据进行分类的方法。包括：最小距离判别法最近邻域分类法,最小距离分类法,最小距离判别法 这种方法要求对遥感图像中每一个类别选一个具有代表意义的统计特征量（均值），首先计算待分像元与已知类别之间的距离，然后将其归属于距离最小的一类。最近邻域分类法 这种方法是上述方法在多波段遥感图像分类中的推广。在多波段遥感图像分类中，每一类别具有多个统计特征量。最近邻域分类法首先计算待分像元到每一类中每一个统计特征量间的距离，这样，该像元到每一类都有几个距离值，取其中最小的一个距离作为该像元到该类别的距离，最后比较该待分像元到所有类别间的距离，将其归属于距离最小的一类。,最小距离分类法,最小距离分类法原理简单，分类精度不高，但计算速度快，它可以在快速浏览分类概况中使用。,多级切割法,多级切割法（multi-level slice classifier）是根据设定在各轴上的值域分割多维特征空间的分类方法。这种方法要求通过选取训练区，详细了解分类类别（总体）的特征，并以较高的精度设定每个分类类别的光谱特征上限值和下限值，以便构成特征子空间。对于一个未知类别的像素来说，它的分类取决于它落入哪个类别特征子空间中。如落入某个特征子空间中，则属于该类，如落入所有特征子空间之外，则属于未知类型，因此多级切割分类法要求训练区样本的选择必须覆盖所有的类型，在分类过程中，需要利用待分类像素光谱特征值与各个类别特征子空间在每一维上的值域进行内外判断，检查其落入哪个类别特征子空间中，直到完成各像素的分类。,遥感图象处理：分类,用多级切割法分割三维特征空间,多级切割法,多级分割法分类便于直观理解如何分割特征空间，以及待分类像素如何与分类类别相对应。但它要求分割面总是与各特征轴正交，如果各类别在特征空间中呈现倾斜分布，就会产生分类误差。因此运用多级分割法分类前，需要先进行主成分分析，或采用其他方法对各轴进行相互独立的正交变换，然后进行多级分割。,特征曲线窗口法,特征曲线是地物光谱特征参数构成的曲线。由于地物光谱特征受到大气散射、天气状况等影响，即使同类地物，它们所呈现的特征曲线也不完全相同，而是在标准特征曲线附近摆动变化。因此以特征曲线为中心取一个条带，构造一个窗口，凡是落在此窗口范围内的地物即被认为是一类，反之，则不属于该类，这就是特征曲线法。特征曲线窗口法分类的依据是：相同的地物在相同的地域环境及成像条件下，其特征曲线是相同或相近的，而不同地物的特征曲线差别明显。,特征曲线窗口法,特征曲线选取的方法可以有多种，如地物吸收特征曲线，它将地物的标准吸收特征值连接成曲线，通过与其他像素吸收曲线比较，进行分类；也可以在图像训练区中选取样本，把样本地物的亮度值作为特征参数，连接该地物在每波段参数值即构成该类地物的特征曲线。特征曲线窗口法可以根据不同特征进行分类，如利用标准地物光谱曲线的位置、反射峰或谷的宽度和峰值的高度作为分类的识别点，给定误差容许范围，分别对每个像素进行分类；或者利用每一类地物的各个特征参数上、下限值构造一个窗口，判别某个待分像元是否落入该窗口，只要检查该像元各特征参数值是否落入到相应窗口之内。,特征曲线窗口法,特征曲线窗口法分类的效果取决于特征参数的选择和窗口大小。各特征参数窗口大小的选择可以不同，它要根据地物在各特征参数空间里的分布情况而定。,最大似然比分类法,最大似然比分类法（maximum likelihood classifier）求出像元数据对于各类别的似然度（likelihood），把该像元分到似然度最大的类别中去的方法。似然度是指，当观测到像元数据x时，它是从分类类别k中得到的（后验）概率。它假定训练区地物的光谱特征和自然界大部分随机现象一样，近似服从正态分布，利用训练区可求出均值、方差以及协方差等特征参数，从而可求出总体的先验概率密度函数。当总体分布不符合正态分布时，其分类可靠性将下降，这种情况下不宜采用最大似然比分类法。,最大似然比分类法,最大似然比分类法在多类别分类时，常采用统计学方法建立起一个判别函数集，然后根据这个判别函数集计算各待分像元的归属概率。x为待分像元，P（k)为类别k的先验概率，可以通过训练区来决定。由于上式中分母和类别无关，在类别间比较的时候可以忽略。,最大似然比分类法,最大似然比分类必须知道总体的概率密度函数P（x|k）。由于假定训练区地物的光谱特征和自然界大部分随机现象一样，近似服从正态分布（对一些非正态分布可以通过数学方法化为正态问题来处理），通过训练区，可求出其平均值及方差、协方差等特征参数，从而可求出总体的先验概率密度函数。此时，像素x归为类别k的归属概率Lk表示如下（这里省略了和类别无关的数据项）：,类别k的协方差矩阵,类别k的平均向量（n维）,最大似然比分类法,这种最大似然比分类法的特征是，在分类结果上具有概率统计的意义。但必须注意几点：（l）为了以较高精度测定平均值及方差、协方差，各个类别的训练数据至少也要为特征维数的2到3倍以上。（2）如果2个以上的波段相关性很强，那么方差协方差矩阵的逆矩阵就不存在，或非常不稳定。在训练数据几乎都取相同值的均质性数据组的情况下也是如此。此时，最好采用主成分分析法，把维数减到仅剩相互独立的波段。（3）当总体分布不符合正态分布时，不适于采用以正态分布的假设为基础的最大似然比分类法。其分类精度也将下降。,非监督分类,前提是假定遥感影像上同类物体在同样条件下具有相同的光谱信息特征。非监督分类方法不必对影像地物获取先验知识，仅依靠影像上不同类地物光谱信息（或纹理信息）进行特征提取，再统计特征的差别来达到分类的目的，最后对已分出的各个类别的实际属性进行确认。非监督分类主要采用聚类分析方法，聚类是把一组像素按照相似性归成若干类别，即“物以类聚”。它的目的是使得属于同一类别的像素之间的距离尽可能的小而不同类别上的像素间的距离尽可能的大。其常用方法有：分级集群法动态聚类法,分级集群法,当同类物体聚集分布在一定的空间位置上，它们在同样条件下应具有相同的光谱信息特征，其他类别的物体应聚集分布在不同的空间位置上。由于不同地物的辐射特性不同，反映在直方图上会出现很多峰值及其对应的一些灰度值，它们在图像上对应的像元分别倾向于聚集在各自不同灰度空间形成的很多点群，这些点群就叫做集群。分级集群法采用“距离”评价各样本（每个像元）在空间分布的相似程度，把它们的分布分割或者合并成不同的集群。每个集群的地理意义需要根据地面调查或者与已知类型的数据比较后方可确定。,分级集群法,分级集群法的分类过程如下：确定评价各样本相似程度所采用的指标，这里可以采用前面监督分类中介绍的几种距离。初定分类总数n。计算样本间的距离；根据距离最近的原则判定样本归并到不同类别。归并后的类别作为新类，与剩余的类别重新组合，然后再计算并改正其距离。在达到所要分类的最终类别数以前，重复样本间相似度的评价和归并，这样直到所有像素都归入到各类别中去。,确定采用的距离,确定分类总数n,找出距离最小的类别组,归并距离最小的类别,计算归并后新的个体间的距离,归并后的类别数,STOP,分级集群法,分级集群方法的特点是这种归并的过程是分级进行的，在迭代过程中没有调整类别总数的措施，如果一个像元被归入到某一类后，就排除了它再被归入到其他分支类别中的可能性，这样可能导致对一个像元的操作次序不同，会得到不同的分类结果，这是该方法的缺点。,动态聚类法,在初始状态给出图像粗糙的分类，然后基于一定原则在类别间重新组合样本，直到分类比较合理为止，这种聚类方法就是动态聚类。ISODATA(Iterative Orgnizing Data Analysize Technique迭代自组织数据分析技术）方法在动态聚类法中具有代表性。,ISODATA方法,按照某个原则选择一些初始类聚类中心。在实际操作中，要把初始聚类数设定得大一些，同时引入各种对迭代次数进行控制的参数，如控制迭代的总次数、每一类别最小像元数、类别的标准差、比较相邻两次迭代效果以及可以合并的最大类别对数等，在整个迭代过程中，不仅每个像元的归属类别在调整，而且类别总数也在变化。在用计算机编制分类程序时，初始聚类中心可按如下方式确定：设初始类别数为n，这样共有n个初始聚类中心，求出图像的均值M和方差，按下式可求出初始聚类中心:,k=1,2,n，为初始类中心编号，n为初始类总数。,ISODATA方法,计算像素与初始类别中心的距离，把该像素分配到最近的类别中。动态聚类法中类别间合并或分割所使用的判别标准是距离，待分像元在特征空间中的距离说明互相之间的相似程度，距离越小，相似性大，则它们可能会归入同一类。这里的距离可以采用前面介绍的几种距离。计算并改正重新组合的类别中心，如果重新组合的像素数目在最小允许值以下，则将该类别取消，并使总类别数减1。当类别数在一定的范围，类别中心间的距离在阈值以上，类别内的方差的最大值为阈值以下时，可以看做动态聚类的结束。当不满足动态聚类的结束条件时，就要通过类别的合并及分离，调整类别的数目和中心间的距离等，然后返回到上一步，重复进行组合的过程。,动态聚类法,动态聚类法中有类别的合并或分裂，这说明迭代过程中类别总数是可变的。如果两个类别的中心点距离近，说明相似程度高，两类就可以合并成一类；或者某类像元数太少，该类就要合并到最相近的类中去。类别的分裂也有两种情况：某一类像元数太多，就设法分成两类；如果类别总数太少，就将离散性最大的一类分成两个类别，可以先求出每个类别的均值和标准差，然后通过对每一个波段的标准偏差设定阈值来实现，标准差大于阈值，该类就要分裂。,分类后处理,无论是监督分类还是非监督分类，其结果都会产生一些面积很小的图斑。无论从专题制图的角度还是实际应用的角度，都有必要对这些小图斑进行剔除。处理方法：聚类统计(Clump)过滤分析(Sieve)去除分析(Eliminate)分类重编码（Recode),分类后处理,聚类统计：通过对分类专题图像计算每个分类图斑的面积、记录相邻区域中最大图斑面积的分类值等操作，产生一个Clump类组输出图像，其中每个图斑都包含Clump类组属性。这是一个中间结果，供下一步处理使用。,分类后处理,过滤分析：对经Clump处理后的Clump类组图像进行处理，按照定义的数值大小，删除Clump图像中较小的类组图斑，并给所有小图斑赋予新的属性值0。显然，这引出了一个小图斑归属问题。可以与原分类图对比确定新属性。,分类后处理,去除分析：用于删除原始分类图像中的小图斑或小Clump类组，与过滤不同，去除将删除的小图斑合并到相邻的最大分类中，而且如果输入图像是Clump聚类图像，经过去除处理后，将分类图斑的属性值自动恢复为Clump处理前的原始分类编码。即结果是简化的分类图像。,分类后处理,分类重编码：主要是针对非监督分类而言的，因在非监督分类过程中，用户一般要定义比最终需要多一定数量的分类数；在完全按照像元灰度值通过ISODATA聚类获得分类方案后，首先是将专题分类图像与原始图像对照，判断每个类别的专题属性，然后对相似或类似的分类通过图像重编码进行合并，并定义分类名称和颜色。分类重编码还可以用在其它方面，作用有所不同。,分类后处理,监督/非监督分类方法比较,根本区别点在于是否利用训练样区来获取先验的类别知识监督分类根据训练样区提供的样本选择特征参数，建立判别函数，对待分类点进行分类。因此，训练场地选择是监督分类的关键。对于不熟悉区域情况的人来说，选择足够数量的训练样区带来很大的工作量，操作者需要将相同比例尺的数字地形图叠在遥感图像上，根据地形图上的已知地物类型圈定分类用的训练样区。由于训练样区要求有代表性，训练样本的选择要考虑到地物光谱特征，样本数目要能满足分类的要求，有时这些还不易做到，这是监督分类不足之处。,监督/非监督分类方法比较,相比之下，非监督分类不需要更多的先验知识，它根据地物的光谱统计特性进行分类。因此，非监督分类方法简单，且分类具有一定的精度。严格说来，分类效果的好坏需要经过实际调查来检验。当光谱特征类能够和唯一的地物类型（通常指水体、不同植被类型、土地利用类型、土壤类型等）相对应时，非监督分类可取得较好分类效果。当两个地物类型对应的光谱特征类差异很小时，非监督分类效果不如监督分类效果好。,存在的问题,遥感图像计算机分类算法设计的主要依据是地物光谱数据。因此，存在着如下的问题：未充分利用遥感图像提供的多种信息遥感数字图像计算机分类的依据是像素具有的多光谱特征，并没有考虑相邻像素间的关系。例如，被湖泊包围的岛屿，通过分类仅能将陆地与水体区别，但不能将岛屿与临近的陆地（假定二者地面覆盖类型相同，具有同样的光谱特征）识别出来。这种方法的主要缺陷在于地物识别与分类中没有利用到地物空间关系等方面的信息。,存在的问题,未充分利用遥感图像提供的多种信息统计模式识别以像素作为识别的基本单元，未能利用图像中提供的形状和空间位置特征，其本质是地物光谱特征的分类。例如，根据水体的光谱特征，在分类过程中可以识别构成水体的像素，但计算机无法确定一定空间范围的水体究竟是湖泊还是河流。这个问题如果引入地物形状特征则可以识别。显然，遥感图像计算机分类未能充分利用遥感图像提供的多种信息。因此图像分类后，可以利用分类的结果，将这些目标对象进行重组，在区域分割或边界跟踪的基础上抽取遥感图像形态、纹理特征和空间关系等特征，然后利用这些特征对图像进行解译。,存在的问题,提高遥感图像分类精度受到限制 这里的分类精度是指分类结果的正确率。正确率包括地物属性被正确识别，以及它们在空间分布的面积被准确度量。遥感数字图像分类结果在没有经过专家检验和多次纠正的情况下，分类精度一般不超过90，其原因除了与选用的分类方法有关外，还存在着制约遥感图像分类精度的几个客观因素：大气状况的影响 不少人理想化地认为遥感图像只记录遥感观测区域内的地物电磁辐射能量，遥感图像的灰度大小及其变化只反映了地物的辐射光谱特征变化，这种观念是不正确的。,存在的问题,提高遥感图像分类精度受到限制大气状况的影响 地物辐射电磁波，必须经过大气层才能到达传感器，大气的吸收和散射会对目标地物的电磁波产生影响，其中大气吸收使得目标地物的电磁波辐射被衰减，到达传感器的能量减少，散射会引起电磁波行进方向的变化，非目标地物发射的电磁波也会因为散射而进入传感器，这样就导致遥感图像灰度级产生一个偏移量。对多时相图像进行分类处理时，由于不同时间大气成分以及湿度不同，散射影响也不同，因此遥感图像中的灰度值不完全反映目标地物辐射电磁波的特征。为了提高遥感图像分类的精度，必须在图像分类以前进行大气纠正。,存在的问题,提高遥感图像分类精度受到限制下垫面的影响 下垫面的覆盖类型和起伏状态对分类具有一定影响。下垫面的覆盖类型多种多样，受传感器空间分辨率限制，农田中的植被、土壤和水渠，石质山地稀疏的灌丛和裸露的岩石均可以形成混合像元，它们对遥感图像分类的精度影响很大。这种情况可以在分类前首先进行混合像元分解，把它们分解成子像元后再分类。分布在山区向阳面与背阳面的同一类地物，单位面积上接收太阳光能不同，地物电磁波辐射能量也不同，其灰度值也存在差异，容易造成分类错误。在地形起伏变化较大时，可以采用比值图像代替原图像进行分类，以消除地形起伏的影响。,存在的问题,提高遥感图像分类精度受到限制其他因素的影响图像中的云朵会遮盖目标地物的电磁波辐射，影响图像分类。对于图像中仅有少量云朵时，分类前可以采用去噪音方法进行清除。多时相图像分类时，不同景的图像由于成像时光照条件的差别，同一地物电磁波辐射量存在差别，这也会对分类产生影响。地物边界的多样性，使得判定类别的边界往往是很困难的事。例如，湖泊和陆地具有明确的界线，但森林和草地的界线则不明显，不少地物类型间还存在着过渡地带，要精确将其边界区别出来，并非是一件容易的事。因此，提高遥感图像分类精度，既需要对图像进行分类前处理，也需要选择合适的分类方法。,遥感图像多种特征提取,遥感图像解译，除了利用地物的光谱特征外，还需利用地物的形状特征和空间关系特征，因此需要提取图像的其他特征。对于高分辨率遥感图像，可以清楚地观察到丰富的结构信息，如城市是由许多街区组成的，每个街区又由多个矩形楼房构成，其中人造地物具有明显的形状和结构特征，如建筑物、厂房、农田田埂，因此可以设法去提取这类地物的形状特征及其空间关系特征，以作为结构模式识别的依据。结构模式识别也叫句法模式识别，它可以提取地物的形状特征和空间关系特征，在此基础上识别遥感影像上的目标地物。,遥感图像多种特征提取,地物边界跟踪法形状特征描述与提取地物空间关系特征描述与提取,地物边界跟踪法,地球表面的物体，从其分布特征来看，主要表现为三种形式，点状地物线状地物面状地物 随着遥感影像空间分辨率的变化，相同地物的分布特征是可变的。例如，一种地物在高分辨率图像上表现为面状地物，但在低分辨率图像上可能表现为点状地物。反之，点状地物在高分辨率图像上也可以表现为面状地物。从信息论的观点看，一个图像颜色（地物光谱）特征均一的地物单元，其边界信息最丰富，地物形状特征是通过边界信息表现出来的。,地物边界跟踪法,对地物边界跟踪，是获取地物形态特征的前提。针对不同地物分布特点，地物边界跟踪方法不同。第一种方法以图像像元作为跟踪的落脚点，跟踪点的连线作为地物的界线。这种跟踪方法适用于线状物体的跟踪。第二种方法认为地物的界线在相邻地物之间，因此边界跟踪的路径应该从两个相邻地物边界的像元中间穿过，这种方法适用于点状地物与面状地物的跟踪。,44|44,44 4,33 3,33|33,地物边界跟踪法,点状地物与面状地物的边界跟踪 点状地物是用一个像素或几个相邻像素表示的地物，因一个像素对应地表一定的面积，因此点状地物可视为缩小的面状地物。面状地物边界跟踪法同样适用于点状地物，只不过点状地物跟踪的方法更加简单。面状地物的边界跟踪法 面状地物边界跟踪是在遥感图像分类基础上进行的。遥感数字图像分类实际上是根据地物物理特性对像素的重新组合与归并。经过分类后的图像，每个地物单元内部是均一的，与其他地物的差异主要通过地物单元的边界来表现。据此可以认为，特征均一的地物单元（区域）的空间分别是由边界所确定的，边界由一个或多个孤段所构成，弧段是具有方向性的线段，线段由有限个有序的边界点所构成。,地物边界跟踪法,设G1为图像中第一种地物类型，G2为图像中第二种地物类型，像元A0，A1，A2，A3在二维图像空间分布中具有如下：若AG1或AG2 A=A0，A1，A2，A3则称 4个 像元在 2*2图像窗口的中间点为内部点，若存在AiG1,AjG2,ij,则称4个像元在2*2图像窗口的中间点为边界点，地物单元边界跟踪法是：起始点的确定：设图像数据具有m行和n列，分别在图像四周各增加一行或一列数据，增加的像元赋值为-1，构成（m2）*（n2）的数字图像，经过上述处理，图像周边的像元都成为边界点。我们规定，第一个地物单元的边界点在图像坐标原点，将坐标原点作为起始点，开始跟踪时，记下起始点的坐标，顺序扫描图像，查找边界点。,举例,地物边界跟踪法,下一个跟踪点的确定：查找下一个跟踪点，关键是确定下一个跟踪点的方向。通过对地物单元边界点在图像窗口各种可能出现状态的分析，发现在逆（或顺）时针方向查找时，当前区域的下一个边界点是唯一的，寻找当前区域的下一个边界点的方向依赖于当前点的位置和图像窗口内相邻像元的分布位置。利用已知的方向作指引，可以方便地查找到当前区域的下一个边界点。,举例,2 2 1 1 1 1 1 1 1*1 1 1 1 2 1 1 2 2 2,上一点位置,查找下一边界点的方向,当前点位置,地物边界跟踪法,利用方向指引，寻找到下一个边界点后，要求判别该点的坐标是否与起始点相同。不同时，继续寻找当前地物单元下一个边界点，相同时，记录下跟踪过程中边界点的各种信息，产生一个地物单元数据记录，然后利用同样跟踪方法，开始下一个地物单元边界的跟踪，一直到图像中所有地物单元跟踪完毕。上述跟踪方法可以拓广到多个地物类型的图像区域。,举例,地物边界跟踪法,线状地物信息检测与跟踪 线状地物检测是对线状地物像素的亮度与方向特性进行检测，找出可能的线状地物点的算法。线状地物检测的方法较多，简单的有阈值技术，复杂的有诸如假说检验程序等。线状地物检测的输出结果是线状地物种子点。线状地物信息检测首先对数字图像进行增强处理，然后对图像进行二值化处理，图像的二值化方法一般是通过设定某个阈值，并以该阈值为门限，把具有灰度级的图像变换为二值的黑白图像。通过二值化处理，将图像中的像元分为线状地物与背景，在此基础上进行噪声消除，从而获得较为正确的线状地物片断，作为下一步线状地物扩展与连接的起始点。数据经过以上处理后，得到一幅初始的线状地物图。,地物边界跟踪法,线状地物信息检测与跟踪 理想状态下，数字图像中代表线状地物的像素是连续的，这时可以采用逐行扫描的方法在二值图像上寻找每一行上代表线状地物的像素，根据像素的多寡计算出线状地物的宽度。在保证连通性的前提下对构成线段的像素进行剥皮细化，留下骨干像素，然后对下一行进行处理，寻找代表线状地物的像素，应用同样的方法，对线状地物进行跟踪。,地物边界跟踪法,线状地物信息检测与跟踪 在实际遥感图像计算机解译中，遥感图像中线状地物多以间断线状地物种子点形式出现，因此需要利用相关信息或全局约束来扩展线状地物种子点以形成线状地物片断。在此基础上，可以应用更多的全局信息及知识来进一步将线状地物片断扩展、连接成为更为完整的线状地物图。,形状特征描述与提取,地物形状特征的描述 通过边界跟踪可以获得一系列有序的边界点，这些边界点提供了地物单元形状特征的大量信息。这里主要介绍采用链码来记录和描述边界点的方法。链码是一组具有方向标志码的有序系列，A=a1，a2，a3，,an。它是由中心像素指向其8个邻点的方向来定义的。方向标志码可按逆（或顺）时针方向定义。,0,1,2,3,4,5,6,7,0,1,3,4,5,7,2,6,形状特征描述与提取,地物形状特征的描述 链码值在0,7，采用链码方法记录地物边界线，在地物单元跟踪时方便实用，它既可用来表示一条边界线相邻像素点之间的位置，又可以在边界跟踪过程中控制各像素邻点的检测顺序。为便于链码与地理坐标间的转换，规定地理坐标原点取在屏幕左上角，即X轴由左向右，y轴由上而下。对一个连通的像素序列来说，与某一像素连通的后续像素只可能是由07中的任一个方向码。因此，对图像边界来说，除像素序列的初始像素外，任一后续像素位置均可用07中的一个数值来唯一地确定。,形状特征描述与提取,地物形状特征的描述链码到像元的坐标换算方法为：x0,y0为边界起始点的坐标；xn,yn为当前链码所指向的边界点坐标；ajx，ajy分别为aj在x坐标和y坐标的分量。,形状特征描述与提取,地物形状特征的描述像元坐标到链码换算方法为：x1,y1分别为当前点在x坐标和y坐标分量；x0,y0分别为上一点在x坐标和y坐标分量；F(ai)为在不同状态时的对应函数值，即该点的链码，这里的链码是采用逆时针方向定义的。利用链码可以有效地记录地物边界信息，描述地物形态特征。,形状特征描述与提取,地物形态特征的提取地物单元周长P：设相邻像素间采用链码表示的长度为：n=Mod(2,i),i=1,2,7，为链码的方向。，j为地物边界像素点的个数。将所有链码段的长度相加，即为周长。运用计算周长的方法，也可以计算出线状地物的长度。地物面积S：可以利用边界点上的地理坐标求算面积，其面积提取方法如下：首先计算出面状地物每一行的面积Sk,像元地面分辨率,xi,xj为同一行边界点的x地理坐标，m为地物总列数。,形状特征描述与提取,地物形态特征的提取线状物体的曲率：面状地物的形状系数,线段长度,弦长,面状物体周长,面状地物面积,地物空间关系特征描述与提取,这里所讲的地物空间关系是指遥感数字图像中两个地物或多个地物之间在空间上的相互联系，这种联系是由地物的空间位置所决定的。不同地物之间的空间关系在二维空间，地物的空间关系主要表现为以下几种：方位关系：指两个地物之间方向与位置的相对关系。方位关系用来描述边界并不相互接触的两个物体。通常采用以一个物体为中心，描述另一个物体位于它的哪个方向上，距离它有多远。,地物空间关系特征描述与提取,不同地物之间的空间关系 方位关系的描述包括以下内容：距离关系：即一个物体到另一个物体的直线距离。由于空间分布的地物具有三种类型，因此，各种物体之间的距离关系定义也不相同。点状地物之间的距离则是两点间的距离，点状地物到线状地物的距离是该点到该线上某一点的最短垂直距离。点状地物到面状地物的最短距离为该点到面状地物边界的最短距离。线状地物到面状地物的最短距离是线上一点到面状地物边界点的最短距离，面状地物到面状地物的距离是两个面状地物边界点的最短距离。方向关系：即一个物体相对于另一个物体的方向。方向关系常用八个方向来描述，它们分别为：正北、东北、正东、东南、正南、西南、正西、西北。每个方向可以用方位角区间来定量表示。,地物空间关系特征描述与提取,不同地物之间的空间关系,不同地物类型之间的距离关系,地物空间关系特征描述与提取,不同地物之间的空间关系包含关系：一个物体位于另一个物体内部，并且边界不相邻。包含关系具有三种情况：点包含在面状地物内部，线状地物被包含在面状地物内部，一个小的面状物体被另一个大的面状物体所包围。相邻关系：指两个地物在边界上相邻。两个面状物体的相邻关系，存在着两种不同状况：外接邻域，内接邻域。点与面相邻是指点状地物位于面状地物的边界，线状地物与面状地物相邻是指线状地物上一点或多点位于面状地物边界。,地物空间关系特征描述与提取,不同地物之间的空间关系相交关系：两个地物在一点上交汇，它主要用来描述点状地物与线状地物，线状地物与线状地物的空间关系。它包含两种情况：点状地物位于线状地物的某一点，两条线状地物相交一点或相交多点。相贯关系：一个线状物体通过面状物体的内部，例如穿过林区的道路。,地物空间关系特征描述与提取,空间关系特征提取与描述方位关系的提取：计算距离DA、B为空间上两点。确定方位：设地物A、B分布在二维图像空间内，令x轴正向指向东，y轴正向指向北，分别对地物A、B作x与y方向正交投影，可以得到地物A、B在X轴上的坐标和在Y轴向的坐标，比较地物投影在X和Y方向的坐标大小，可以获得地物A相对于地物B的方位关系，或地物B相对于地物A的方位关系。,地物空间关系特征描述与提取,空间关系特征提取与描述包含关系特征提取与描述 提取点状地物与面状地物的包含关系，关键是判明点状地物是否为面状地物所包含，有两种方法可以判断点状地物是否在区域内：铅垂线法：设一个多边形F=（p1，p2，p3，pn）由有序的n个点P1，P2，P3，Pn联接构成其多边形边界。因为区域闭合，有p1=pn,设一个点状地物为p，由p作一条铅垂线，如果铅垂线与多边形相交的边数为偶数，则该点在多边形之外，否则在多边形之内。使用铅垂线法需要注意点状地物与多边形边界相交的异常情况。,地物空间关系特征描述与提取,空间关系特征提取与描述包含关系特征提取与描述射线法：若点状地物在多边形内部，由该点向任意方向作射线，必然与面状地物边界相交，据此可以判断点状地物在多边形内部，否则点状地物不被该面状地物包含。提取线状地物与面状地物的包含关系，关键是判明线状地物是否为面状地物所包含，假设线状地物为面状地物所包含，此时组成线段的各点也必然为面状地物所包含，这样可以将提取线状地物与面状地物的包含关系的运算简化，可以利用铅垂线法或射线法检测“线”上的一点是否在多边形内部。,地物空间关系特征描述与提取,空间关系特征提取与描述包含关系特征提取与描述 一个面状物体是否包含另一个面状物体，只需要判明这个面状物体的边界是否在另一个面状物体内部，此时把一个面状物体是否包含另一个面状物体的问题简化为线状地物与面状地物的包含关系，进一步简化为点状地物与面状地物的包含关系。这仍然可以采用铅垂线法或射线法来检测。,地物空间关系特征描述与提取,空间关系特征提取与描述相邻关系特征提取包括三种不同情况：点与面相邻，这可以通过检测“点”是否在多边形的边界上来确认。线状地物与面状地物相邻，首先需要了解线状地物与面状地物边界是否相交，如果存在相交，那么以相交点为裁剪点，将线状地物一分为二，分别检测这两个线段是否同时在面状地物的外部或在面状地物的内部。若同时在面状地物的外部或在面状地物的内部，说明线状地物与面状地物相邻，否则不是。,地物空间关系特征描述与提取,空间关系特征提取与描述相邻关系特征提取两个面状地物相邻问题，可以采取更为简单