《遥感技术与应用遥感专题信息提取精品文稿.ppt》由会员分享，可在线阅读，更多相关《遥感技术与应用遥感专题信息提取精品文稿.ppt（47页珍藏版）》请在taowenge.com淘文阁网|工程机械CAD图纸|机械工程制图|CAD装配图下载|SolidWorks_CaTia_CAD_UG_PROE_设计图分享下载上搜索。
1、遥感技术与应用遥感专题信息提取第1页，本讲稿共47页 二、图像解译的背景知识1专业知识 指需要熟悉所解译的学科及相关学科的知识。包括对地物成因联系、空间分布规律、时相变化以及地物与其他环境要素间的联系等知识。比如遥感地质探矿，首先需具备地层、构造、蚀变带等与找矿直接相关的地质知识和经验，还需要了解植物分布、土壤等相关知识。2地理区域知识 指区域特点、人文自然景观等。每个区域有其独特的区域特征，即地域性，它影响到图像的图型结构等，它能直接帮助识别地物和现象3遥感系统知识 解译者必须了解每一图像是怎样生成的；不同遥感器是如何描述景观特征的，它使用的何种电磁波段，空间分辨率是多少等等，怎样从图像中获
2、得有用的信息等遥感图像解译与日常目视观察有3点不同：A 遥感图像通常为“顶视-鸟瞰”，而不同于日常生活中的透视B 遥感图像常用可见光以外的电磁波段，而大多数我们熟悉的特征是在可见光谱段，可以表现得十分不同C 遥感图像常用一种不熟悉或变化的比例和分辨率描述地球表面第2页，本讲稿共47页因此，对于初学者需要多对照地形图、实地或熟悉地物观测，以增强立体感和和景深印象，纠正视觉误差，积累图像判读经验遥感图像的解译过程是地面目标成像过程的逆过程，即从地面实况的模拟影像中提取遥感信息、反演地面原型的过程三、遥感图像解译的两种方式1目视解译（Visual Interpretation）借助人眼的观察和人的智
3、能，结合一些量测工具（测微尺、放大镜、立体镜等）来识别图像中的目标。特点：解译者的学识和经验在判读中起主要作用，精度高，但难以对海量空间信息的定量化分析2 计算机图像理解（Computer Imagery Understanding）以计算机软硬件为支撑，利用模式识别技术（Pattern Recognition）和人工智能（Artificial Intelligence）技术，建立图像解译专家系统（Experts System)从而让计算机模拟人工解译过程，“读取遥感图像上的特征”，进而确定图像上的目标。特点：速度快，处理方式灵活多样，整个处理过程通常是以人机交互方式进行，对计算机技术和算法要
4、求较高，识别的精度通常不及目视解译第3页，本讲稿共47页第二节 遥感图像目视解译标志一、解译标志的概念解译标志（Interpretation Mark）：遥感图像上能帮助人们识别不同目标的那些影像特征；直接解译标志：能在遥感影像上直接看到可供判读的影像特征称为直接判读标志，如形状、大小、阴影、纹理、色调等；间接解译标志：指运用某些直接解译标志，根据地物的相关属性等地学知识，间接推断出的影像标志，如根据道路与河流相交处的特殊影像特征，可以判断渡口；根据植被、地貌与土壤的关系，来识别土壤类型和分布（如温带针叶林区多为灰化土）二、常用的解译标志1色调（Tone）指地物反射、辐射电磁波的性质和能量强弱
5、在影像上的表现说明两点：a 解译者必须了解该解译图像中影像色调的支配因素 如可见光和近红外图像均反映地物“反射”波谱特征的差异，涉及地物的物质组成、水分含量等；而热红外图像则反映地物“发射”波谱特征的差异，是地物温度差异的记录；雷达图像反映地物后向散射能量的差异，涉及到地物介电常数、表面粗糙程度等物理性质 b 影像色调受到多种因素的影响 除了受目标本身的波谱特征因时因地、因环境变化而变化外，还受到成像高度、成像时间（光照角度、强度）、遥感器第4页，本讲稿共47页观察角度、成像后处理等多种因素的影响，因此使用色调标志时要特别小心影响影像色调变化的因素：（1）地物本身的颜色（2）地物表面的结构 表
6、面结构不一致，其光滑程度有差异，反光能力就不同，反映 在影像上的色调就不相同 a 光滑表面（镜面反射 Specular Reflection）若传感器照射方向与太阳入射 方向一致，此时影像呈暗色调；若传感器照射方向逆着镜面反射的方向，则 影像呈亮色调 b 无光泽表面（粗糙表面，漫反射Diffuse Reflection）各方向上光线被均匀 散射开来，因此无论哪个方向得到的影像色调均一 c 起伏不平的表面 阳坡（Adret-Sunny Slope）影像色调亮 阴坡（Ubac-Shady Slope）影像色调暗（3）地物本身的反光能力（4）湿度的大小 通常湿度越大则影像色调越暗 利用此特征可以监测
7、土壤湿度（5）摄影季节的不同 春季植物色调浅，夏季色调深第5页，本讲稿共47页2阴影（Shadow）指因倾斜照射，地物自身遮挡光源而造成影像上的暗色调阴影分为本影和投落阴影，前者是地物背光面在影像上呈现的暗色调构像，有助于获得目标的立体感；后者是地物背光方向地物投射到地面的暗斑在影像上的构像，它反映地物的侧面形态，可以根据侧影的长度及照射角度反算地物的高度 但是，因为阴影的遮盖会掩盖一些信息，给解译工作带来不便3大小（Size）指地物长、宽、高、面积、体积等在影像上的记录 地物影像的大小不仅能求得地物的数量特征，而且对判读地物的性质也有帮助。例如，单轨铁路和双轨铁路从形状上往往不易区分，但通过
8、量算其宽度则很容易区分地物影像的大小取决于影像的比例尺大小。因此，解译前要获得影像的比例尺4形状（Shape）是目标物的外形、轮廓 a 遥感图像记录的是目标物的顶面形状，不同于地面看的侧面形状 b 自然形成物往往形状复杂，不规则，轮廓边界呈现较圆滑自然的曲线 如河曲 冲击扇、火山锥等；人为活动的产物一般具有规则的几何形状，拐角明显呈 一定角度，如田块、果园、城镇建筑物等 c 影像的比例尺和分辨率对形状影响较大，小比例尺时，往往呈现的是综合体的形状第6页，本讲稿共47页5纹理（Texture）是指具有不同色调和形状的细微影像按某种规则排列组合的集合，这种细微影像称为纹理基元。它往往是指图像上地物
9、表面质感（平滑、粗糙、细腻等）如“花格子衬衫”（强调单个事物的重复规则排列，如果园果树的株行距、集约化程度高的田块排列，贵州山区梯田的排列）6图案（Pattern），即图型结构，是指多个人工目标重复排列的空间形式，它反映地物目标的空间分布特征。如建筑物、道路、绿化带形成的住宅区的图案7位置（Site）即地理位置，它反映地物所处的地点与环境，地物与周边的关系 如菜地多分布于居民点周围与河流两侧；机场多在大城市郊区的平坦地；堤坝与道路（形色难分），堤在河流两侧并与之平行，道路与居民点相连；有的植物生长于高地，有些只能生长在湿地8组合（Association）指某些目标的特殊表现和空间组合关系，它不
10、同于那种严格按图型结构显示的空间排列，而指物体间一定的位置关系和排列方式，即空间配置和布局如砖场由高烟囱、取土坑、堆砖场等组合而成；军事目标可能有雷达站、军车、军营及周围的配套的军事设施等第7页，本讲稿共47页第三节 遥感图像的目视解译一、目视解译的原则1从应用的目的出发，总体观察，全面分析图像特征。坚持先易后难，由粗入细，由 整体到局部的原则2应充分利用各种解译标志，包括直接标志和间接标志，相互补充，彼此验证3随着多光谱和多平台遥感技术的发展，应尽可能创造条件开展多波段、多时相、多类型遥感图像的对比分析二、目视解译的基本过程第8页，本讲稿共47页1准备阶段 根据解译的目的和经费的支撑选择适当
11、空间分辨率、时相、光谱分辨率的遥感数据并进行波段组合进行彩色增强；尽量收集到质量好、现势性较强的地形图及有关专业图件和文字资料；编写解译标志表及说明草稿（这是一编检索表的过程，工作量大）；野外选择典型地区进行粗查，为建立合理的解译标志表打好基础。2室内图像解译阶段 按解译的基本原则进行，对无法解译或把握不大的区域应记录下来，有待野外验证3野外验证与补充解译阶段4成果整理与制图总结阶段三、遥感解译过程的复杂性（Complexity）1 遥感图像显示的是某一特定地理环境的综合体，这种综合表现在两个方面：a 地理要素的综合 它反映的是地质、地貌、水文、土壤、植被、社会生态等 多种自然、人为要素的综合
12、，这些因子是密切相关、交织在一起，往往难以区分 b 遥感信息本身的综合 遥感图像可以是不同时间、空间、光谱、辐射分辨率的综合遥感数据对应的地理环境又是一个复杂的、多要素的、多层次的、具有动态结构和 第9页，本讲稿共47页明显地域差异的开放巨系统，它在时间和空间上是不断变化的，因而遥感信息中的诸要素是相互关联、复杂变化的因此，遥感解译过程是复杂的，主要体现在：1 地物波谱特征是复杂的，它受多种因素控制，本身也因时因地而变化2 自然界存在着大量的“同物异谱”和“同谱异物”现象3 地物的时空属性和地学规律是错综复杂的，各要素、各类别之间的关系是多种类型的地带性规律：由于太阳辐射随纬度分布的规律性，造
13、成沿纬度的水平地带性现象，由于温度、湿度等随地形高度分布的规律性造成沿高度的垂直地带性现象，植物从播种到成熟的季节性变化规律等；有的具有随机性、不确定性，如自然灾害；有的具有模糊性存在过渡渐变关系，如气候带、草场类型的变化等为了提高解译过程的正确性、可靠性，必须补充必要的辅助数据（Ancillary）和先验知识（Priori Knowledge）-指地学、生物学、物理学、数学等专业知识，在GIS的支持下，发展一系列相关的、多层次的、综合的应用分析方法，进行遥感与地学综合分析第10页，本讲稿共47页第四节 模式识别（Pattern Recognition）概述遥感图像的计算机分类是模式识别技术在
14、遥感领域的具体应用，是遥感图像应用处理的重要内容和关键技术之一。首先提及关于模式识别的基本概念一、基本概念1类别（Category）是与概念（Concepts）或原型（Prototype）相关的“自然”状态或 对象种类。假定有c个类别，也称为解释空间（Interpretation Space）。解释空间有概念驱动（Concept-Driven）的属性，如类别的联合、交叉等2 模式（Pattern）是对象的“物理”描述。通常是信号、图像或简单的数值表。也常把模式称为对象、案例或样本。3特征（Feature）是从模式得到的对分类有用的度量、属性或基元（Primitive）所以特征选取是门艺术。二、
15、模式识别的方法1数据聚类（Clustering）是用某种相似性度量的方法将数据组织成有意义的和有用的各组数据。数据聚类不需要利用已知类的信息，它是一种非监督学习的方法，解决方案是数据驱动的第11页，本讲稿共47页2统计分类（Statistical Classification）最为经典的模式识别方法，在选取训练区的基础上，估计出每类型的统计特征参数（均值向量和协方差矩阵），从而确定每类型的概率分布密度函数，进而建立判别函数，实现对待分像元进行分类。这种分类方法是建立在多元变量的多元正态分布的假设基础上的。这是一种监督学习方法。3神经元网络（Artificial Neural Network）是
16、以数学算法结合计算机程序来模拟人的物理神经元之间的联系方式和思维过程，通过反馈控制的方式来调整神经元之间的联系权重，进而使神经元的实际输出无限的逼近期望输出，这一过程通过迭代方式来训练神经元网络。这是一种非线性的数学方法，并且是与“模型无关”，也就是它无须对数据的分布作出假设。存在的不足（1）关于神经网络的数学解释很复杂（2）对于一些重要的决策，设计者常常没有理论依据，因此必须进行大量的试验，并从试验的失败中得到启发从而不断改进神经网络（3）从神经网络得不到任何的语义的信息第12页，本讲稿共47页 4 结构模式识别 其基本思想是：任何一个具有意义的物理模式都具有良好的结构性，它是由许多已知类别
17、和性质的子模式或源模式（也称模式基元），按照一定的规则构成的组合件，所以一旦知道生成规则和子模式、源模式，便可由“模式描述”语句的文法推理出一个用字符串来表示的模式集合。若所研究的模式被包含在此集合内，就认为该模式被识别出来了。此法至今在遥感图像识别分类中尚未得到实际应用，有待进一步开发 5 模糊模式识别 人在进行识别时可以根据一些模糊的印象就能做到比较准确的判断。然而，计算机是建立在0、1二值逻辑基础上的，它对事物的分析的结论是“非假”即“真”。这种二值逻辑不适宜处理模糊事物。根据人辨识事物的思维逻辑，吸取人脑的识别特点，模糊集合论（Fuzzy Set Theory）把数学从二值逻辑转向连续
18、逻辑这就更接近人类大脑的识别活动了。由此产生了一种相当独特的识别方法：模糊识别第13页，本讲稿共47页第五节 遥感图像计算机分类概述一、遥感图像计算机分类的概念就是将遥感图像中每个像元根据其在不同波段的光谱亮度、空间结构特征或其他信息按照某种规则或算法划分为不同的土地覆盖类型专题的技术过程。二、遥感图像计算机分类（模式识别）的种类1光谱模式识别（Spectral PR）是指根据像元到像元的光谱信息来自动划分土地覆盖类型的分类过程的总称。它只考虑了地物在光谱辐射方面的差异。2空间（结构）模式识别（Spatial or Structural PR）是指图像像元根据和它周围像元的空间关系来进行图像分
19、类。空间关系主要是指：图像的结构、像元的近似度特征的大小、形状、方向性、分形维、扁率、重复度和上下关系等空间特征。这种方法比光谱模式识别过程更复杂，而且计算算法上更精深。3时间模式识别（Temporal PR）是指在特征鉴别中将时间作为一种辅助信息，利用多时相的遥感数据动态的进行土地覆盖类型的分类说明：遥感图像分类没有唯一“正确”的方法，因此可以将它们结合在一起运用，形成混合模式（Hybrid Pattern）第14页，本讲稿共47页三、遥感图像计算机分类的方式按照在分类过程中分类者的参与程度，将分类分为两种方式：1监督分类（Supervised Classification）2非监督分类（U
20、nsupervised Classification）四、遥感图像计算机分类的原理这样，图像中的任何一个像元就可以用一n维的光谱值向量来描述：式中：xji代表第j个波段上第i号 像元的亮度值 i=1,2,N，N代表每个 波段图像中像元的个数这样，任何一像元点的光谱向量对应着n维欧氏空间的一个点第15页，本讲稿共47页若某两个像元属于同一地物类型，则它们在欧氏空间中就应该集结在一起；相反，若两个像元分别属于不同的地物类型，则它们在欧氏空间中就应该彼此分离。分类图像（类码标记图像）第16页，本讲稿共47页假设分割A、B两个点集的曲面（图中是两维变量因此是直线）表达式为：则方程 称为A、B两类的判别
21、边界（Decision Boundary）在 已经确定以后，特征空间中任意一点是属于A类？还是属于B类？根据几何学知识可知：当 时 当 时 上述的确定样本点所属类别的规则称为判别准则（Decision Criteria），称为判别函数（Decision Function）遥感图像分类算法的核心就是确定判别函数 和相应的判别准则，为了保证所确定的判别函数能够将各类地面目标在特征空间中的点分割开来，通常是在一定的准则（如Bayes分类器中误分概率最小准则等）下求判别函数和相应的判别准则。第17页，本讲稿共47页五、特征选择与特征变换在遥感图像分类之前，为了压缩参与分类多光谱图像数据的维度和找到那些
22、对识别感兴趣的地物类型最为有效的“特征图像”，通常需要进行特征选择（Feature Selection）和特征提取（Feature Extraction）。在分类时我们总是希望用最少的影像数据取得最好的分类结果，这样就需要从多光谱图像中选择最佳的特征影像组合进行分类，这就称为特征选择。将原始多光谱图像通过一定的数学变换（如主成分、缨帽变换、植被指数变换等）生成一组进行了信息融合的新的特征图像，这种技术手段就称为特征变换。特征选择时，除了凭借研究者的定性判断进行特征选择外，通常可用一些定性的指标计算来帮助进行特征选择，这些指标主要有：1距离测度（Distance Measure）距离是最基本的类
23、别可分性的测度，如果所选择的特征图像能使感兴趣的目标类别的类内距离最小，而与其它类别的类间距离最大，则用这组特征图像设计的分类器分类效果一定很好第18页，本讲稿共47页如果将所有的特征图像都计算这一参数值，通过比较找到参数值最大的那个特征图像，即是最佳的分类特征图像2变换离散度（Transform Divergence）衡量遥感数据各地物类型可分性好坏的数量指标常采用变换离散度和JM距离。离散度（Divergence）是具有分布重叠度的两个类型之间可分性的数量测度。若遥感图像中各类型的光谱矢量遵从或近似遵从多元正态分布，则图像中每两个类型之间的离散度定义为：式中：表示求矩阵的迹，和 分别是第i
24、类型和第j类型的协方差矩阵 和 分别为第i类型和第j类型的光谱均值向量实际应用时，常采用变换离散度公式定义光谱类型间的可分性，变换离散度定义为：第19页，本讲稿共47页当需要研究M个类型之间的可分性时，通常是考察M个类别间的平均离散度，它定义为：式中：M是选定地物的类型数，是第i类型的的先验概率。3最优指数公式OIF（Optimum Index Formula）在前面假彩色合成时已介绍第20页，本讲稿共47页Classification Accuracy AssessmentAnother area that is continuing to receive increased attenti
25、on by remote sensing specialists is that of classification accuracy assessment.Historically,the ability to produce digital land cover classifications far exceeded the ability to meaningfully quantify their accuracy.In fact,this problem sometimes precluded the application of automated land cover clas
26、sification techniques even when their cost compared favorably with more traditional means of data collection.The lesson to be learned here is embodied in the expression“A classification is not complete until its accuracy is assessed”.第21页，本讲稿共47页Classification error matrixOne of the most common mean
27、s of expressing classification accuracy is the preparation of a classification error matrix(sometimes called a confusion matrix or a contingency table).Error matrices compare,on a category-by-category basis,the relationship between known reference data(ground truth)and the corresponding results of a
28、n automated classification.Such matrices are square,with the number of rows and columns equal to the number of categories whose classification accuracy is being assessed.Table 7.3 is an error matrix that an image analyst has prepared to determine how well a classification has categorized a represent
29、ative subset of pixels used in the training process of a supervised classification.This matrix stems from classifying the sampled training set pixels and listing the known cover types used for training(columns)versus the pixels actually classified into each land cover category by the classifier(rows
30、)第22页，本讲稿共47页第23页，本讲稿共47页Several characteristics about classification performance are expressed by an error matrix.For example,one can study the various classification errors of omission(inclusion).Note in the above table that the training set pixels that are classified into the proper land cover ca
31、tegories are located along the major diagonal of the error matrix(running from upper left to low right).All non-diagonal elements of the matrix represent errors of omission or commission.Omission errors correspond to non-diagonal column elements(e.g.,16 pixels that should have been classified as“san
32、d”were omitted from that category).Commission errors are represented by non-diagonal row elements(e.g.,38“urban”pixels plus 79“hay”pixels were improperly included in the“corn”category).Several other descriptive measures can be obtained from the error matrix.For example,the overall accuracy is comput
33、ed by dividing the total number of correctly classified pixels(i.e.,the sum of the elements along the major diagonal)by the total number of reference pixels.Likewise,the accuracies of individual categories can be calculated by dividing the number of pixels in the corresponding row or column.第24页，本讲稿
34、共47页What are often termed producers accuracies result from dividing the number of correctly classified pixels in each category(on the major diagonal)by the number of training set pixels used for that category(the column total).This figure indicates how well training set pixels of the given cover typ
35、e are classified.Users accuracies are computed by dividing the number of correctly classified pixels in each category by the total number of pixels that were classified in that category(the row total).This figure is a measure of commission error and indicates the probability that a pixel classified
36、into a given category actually represents that category on the ground.Note that the error matrix indicates an overall accuracy of 84%.However,producers accuracies range from just 51%to 100%and users accuracies vary from 72%to 99%.Furthermore,this error matrix is based on training data.It should be r
37、emembered that such procedures only indicate how well the statistics extracted from these areas can be used to categorize the same areas.第25页，本讲稿共47页If the results are good,it means nothing more than that the training areas are homogeneous,the training classes are spectrally separable,and the classi
38、fication strategy being employed works well in the training areas.This aids in the training set refinement process,but it indicates little about how the classifier performs elsewhere in a scene.One should expect training area accuracies to be overly optimistic,especially if they are derived from lim
39、ited data sets.A further point to be made about interpreting classification accuracies is the fact that ever a completely random assignment of pixels to classes will produce percentage correct values in the error matrix.In fact,such a random assignment statistic is a measure of the difference betwee
40、n the actual agreement between reference data and an automated classifier and the chance agreement between the reference data and a random classifier.Conceptually,K=(observed accuracy chance agreement)/(1-chance agreement)第26页，本讲稿共47页This statistic serves as an indicator of the extent to which the p
41、ercentage correct values of an error matrix are due to“true”agreement versus“chance”agreement.As true agreement(observed)approaches 1 and chance agreement between 0 and 1.For example,a k value of 0.67 can be thought of as an indication that an observed classification is 67 percentage better than one
42、 resulting from chance.A k value of 0 suggests that a given classification is no better than a random assignment of pixels.In cases where chance agreement is large enough,K can take on negative values-an indication of very poor classification performance.(Because the possible range of negative value
43、s depends on the specific matrix,the magnitude of negative values should not be interpreted as an indication of relative classification performance)第27页，本讲稿共47页六、分类结果的质量评价在遥感图像的分类问题中常以混淆矩阵（Confusion Matrix）来表示分类结果的精度混淆矩阵定义如下：(Usually,this matrix was derived from the validation dataset(ground truth)t
44、he user specified.)式中：mij表示本应属于第i类结果被分入 第j类中去的像元总数 n表示总的类别数所以，该矩阵的主对角线上的元素值越大则说明分类的效果越好第28页，本讲稿共47页若通过分类获取了上述混淆矩阵，可以计算基本的精度估计量：1 总体分类精度 正确分类的总像元数（主对角线上元素的和）与参加分类评价的 总像元数（矩阵中所有元素之和）之比，表示为2用户精度 每一类别被正确分类的像元数量除以被分作该类的总像元数 （列元素之和），表示为3 生产者精度 每一类别被正确分类的像元数量除以该类用作分类训练的总像元 （行元素之和），表示为 第29页，本讲稿共47页一、无监分类的概念
45、在对研究区不了解的情况下，只是依据每一类型地物所具有的相似性（Similarity，类似度），把反映各类型地物特征值的分布按相似分割和概率统计理论将各像元归并成不同的空间集群，然后结合地面实地调查来确定各集群的地物类型从而达到识别分类目的的分类方法。其数学基础是：多元统计的聚类分析（Clustering Analysis）第六节 遥感图像计算机无监分类4Kappa分析第30页，本讲稿共47页二、非监督分类的特点1分类者不需要对研究区有广泛的了解和熟悉，但分类结束后仍需要一定的知识或 地面调查来确定各集群的地物覆盖类型2人为误差的机会减少，因为分类者参与的程度低，只需要事先设定几个参数即可3非监
46、督分类产生的光谱集群组并不一定对应于分析者想要的类别4分析者较难对产生的类别进行控制三、非监督分类的算法1K-Mean算法 其基本思想是通过迭代，逐次移动各类别的中心，直到后一次获得的各类的中心与前一次聚类获得的中心不再发生移动为止。具体步骤：假设图像上的目标要分成m类，m事先由分类者指定（1）适当地选取m个类的初始中心 ，初始中心的选择对聚类的 结果有一定的影响，初始中心的选择一般有以下方法：第31页，本讲稿共47页 1）根据问题的性质，根据经验确定类别数m，从数据中找出从直观上看来比较合适 的m个类的初始中心 2）将全部数据随机地分成m个类别，计算每个类别的重心，将这些重心作为m个类的 初
47、始聚类中心（2）在第K次迭代中，对任一样本X按如下方法把它调整到m个类别中的某一类别中去。对于所有的 ，如果：（3）由第（2）步得到 类新的中心第32页，本讲稿共47页（4）对于所有的 则迭代结束，否则转到第二步继续进行迭代YN开始选择m个聚类中心z1,z2,zm所有像元分到m个聚类中心计算新的聚类中心聚类中心不变结束K均值算法框图说明：这种方法的聚类结果受所选聚类中心的数目及其初始位置以及模式分布的几何性质和读入次序等因素影响，并且在迭代过程中没有调整类别数的措施第33页，本讲稿共47页2ISODATA算法 ISODATA(Iterative Self-Organizing Data Ana
48、lysis Techniques Algorithm，迭代自组织数据分析技术算法），它与K-Mean算法有两点不同：第一，它不是每调整一个样本的类别就重新计算一次各类样本的均值，而是在每次 把所有样本都调整完后才重新计算一次各类样本均值，前者被成为“逐个样本 修正法”，后者被成为“成批样本修正法”第二，ISODATA算法可以自动进行类别的“合并”与“分裂”，从而得到类别数比较合理 的聚类结果，所以类别数是动态的A 分类前预定义参数（1）打算分类的最大类别数Cmax和最小类别数Cmin（2）迭代过程中，最大类别保持不变的像元数量百分比，当达到这个百分比时，ISODATA算法停止（3）每个类别中最
49、小的像元数量、最大的标准方差（4）最小的集群均值间距离，如果两个类别间距离小于这个值，则两类合并（5）集群分散值，通常为0第34页，本讲稿共47页B ISODATA算法的基本步骤（1）初始随机选取Cmax个中心（2）计算其它像元到这些中心的距离，按照最小距离规则将各像元划分到其对应的 类别中去（3）重新计算每个类别的均值，按照前面定义的参数合并或分裂相关类别（4）重复（2）、（3），直到达到最大不变像元百分比或最大迭代次数为止第七节 遥感图像的计算机有监分类一、有监分类的概念利用对研究区已有类别的先验知识从遥感图像上选取若干有代表性的训练区作为样本据此估计出各类别的统计特征参数（主要是指均值向
50、量和协方差矩阵）进而建立判别函数，然后利用判别函数实现对待分类像元进行分类的方法。有监分类的数学基础：多元统计的判别分析（Discriminant Analysis）第35页，本讲稿共47页训练区（Training Area）：是分类者事先已知地物的类别然后从遥感图像上相应位置“圈定”的一定数量像元所构成的样区，其目的是用来对计算机进行“训练”，让计算机学习到描述不同类别地物特征（主要是统计特征，均值向量和协方差矩阵）为建立基于统计学的判别函数打下基础。训练区选取的好坏决定分类的成功与否，因此选训练区是“艰辛的、科学的、艺术的”分类者选训练区的要求：A 代表性 通常在每种类型地物所对应图像的中