Erdas非监督分类全过程.docx

1、Erdas非监督分类全过程Erdas基础教程： 非监督分类来源:陕西师范大学旅游与环境学院1图像分类简介（Introduction to classification) 图像分类就是基于图像像元的数据文件值,将像元归并成有限几种类型、等级或数据集的过程。常规图像分类主要有两种方法：非监督分类与监督分类，专家分类方法是近年来发展起来的新兴遥感图像分类方法，下面介绍这三种分类方法。 非监督分类运用1SODATA（Iterative Self-Organizing Data Analysis Technique ）算法，完全按照像元的光谱特性进行统计分类，常常用于对分类区没有什么了解的情况。使用该方

2、法时.原始图像的所有波段都参于分类运算，分类结果往往是各类像元数大体等比例。由于人为干预较少，非监督分类过程的自动化程度较高。非监督分类一般要经过以下几个步骤：初始分类、专题判别、分类合并、色彩确定、分类后处理、色彩重定义、栅格矢量转换、统计分析。 监督分类比非监督分类更多地要求用户来控制，常用于对研究区域比较了解的情况。在监督分类过程中,首先选择可以识别或者借助其它信息可以断定其类型的像元建立模板，然后基于该模板使计算机系统自动识别具有相同特性的像元。对分类结果进行评价后再对模板进行修改，多次反复后建立一个比较准确的模板，并在此基础上最终进行分类.监督分类一般要经过以下几个步骤：建立模板(训

3、练样本）、评价模板、确定初步分类图、检验分类结果、分类后处理、分类特征统计、栅格矢量转换。 专家分类首先需要建立知识库，根据分类目标提出假设，井依据所拥有的数据资料定义支持假设的规则、条件和变量,然后应用知识库自动进行分类,ERDAS IMAG1NE图像处理系统率先推出专家分类器模块，包括知识工程师和知识分类器两部分，分别应用于不同的情况。 由于基本的非监督分类属于IMAGINE Essentia1s级产品功能、但在1MAGINE Professional级产品中有一定的功能扩展，而监督分类和专家分类只属于IMAGINE ProfeSsiona1级产品，所以，非监督分类命令分别出现在Data

4、Preparation菜单和classification菜单中，而监督分类和专家分类命令仅出现在Classification菜单中。 2 非监督分类（Unsupervised Classification） ERDAS IMAGINE使用ISODATA算法（基于最小光谱距离公式)来进行非监督分类。聚类过程始于任意聚类平均值或一个己有分类模板的平均值:聚类每重复一次,聚类的平均值就更新一次，新聚类的均值再用于下次聚类循环。 ISODATA实用程序不断重复，直到最大的循环次数已达到设定阈值或者两次聚类结果相比有达到要求百分比的像元类别已经不再发生变化。 21分类过程（classification

5、ProcedUre ） 第一步：调出非监督分类对话框 调出非监督分类对话框的方法有以下两种: 方法一:在ERDAS图标面板工具条中，点击Dataprep图标 Data Preparation unsupervised Classification Unsupervised Classification对话框如下：方法二： 在ERDAS图标面板工具条中点击Classifier图标C1assification Unsupervised Classification-unsupervised classification对话框如下：可以看到，两种方法调出的Unsupervised Classific

6、ation对话框是有一些区别的。 第二步：进行非监督分类 在Unsupervised classification对话框中：确定输出文件（Input Raster File）：lazhoucity.img（要被分类的图像） 确定输出文件（Output File）:lzisodat.img（即将产生的分类图像） 选择生成分类摸板文件: Output Signature Set（将产生一个模板文件) 确定分类摸板文件（Filename ): lzisodat。sig 对Clustering options选择Initialize from Statistics单选框 Initialize from

7、 Statistics指由图像文件整体(或其AOI区域）的统计值产主自由聚类，分出类别的多少由自己决定。Use Signature Means是基于选定的模板文件进行非监督分类，类别的数目由模板文件决定。确定初始分类数（Number of classes）: 18分出18个类别) 实际工作中一般将分类数取为最终分类数的2倍以上. .点击Initializing options按钮可以调出Fi1e Statistics Options对话框以设置ISODATA的一些统计参数， .点击Co1or Scheme Options按钮可以调出output color Scheme Options对话框以

8、决定输出的分类图像是彩色的还是黑白的。这两个设置项使用缺省值。 。定义最大循环次数(Maximum Iterations）: 24 最大循环次数（Maximum Iterations）是指ISODATA重新聚类的最多次数，这是为了避免程序运行时间太长或由于没有达到聚类标准而导致的死循环。一般在应用中将循环次数都取6次以上. 设置循环收敛阈值(Convergence Threshold）：0。95 收敛阈值（Convergence Threshold）是指两次分类结果相比保持不变的像元所占最大百分之此值的设立可以避免ISODATA无限循环下去。 点击OK按钮（关闭Unsupervised Cla

9、ssification对话框，执行非监督分类，获得一个初步的分类结果）2.2 分类评价（Evaluate Classification ) 获得一个初步的分类结果以后,可以应用分类叠加（Classification over1ay）方法来评价检查分类精度。其方法如下: 第一步：显示原图像与分类图像 在视窗中同时显示lanzhoucity.img和lz-isodat。img：两个图像的叠加顺序为lanzhoucity.img在下、lz-isodat.img在上,lanzhoucity显示方式用红（4)、绿（5)、蓝（3 )。 第二步：打开分类图像属性并调整字段显示顺序 在视窗工具条中:点击图标（

10、或者选择Raster菜单项-选择Tools菜单）打开Raster工具面板 点击RaSter工具面板的图标（或者在视窗菜单条：Rster-Attributes） 打开Raster Attribute Editor对话框(lz-isodat。img的属性表）属性表中的19个记录分别对应产生的18个类及Unclassified类,每个记录都有一系列的字段.如果想看到所有字段,需要用鼠标拖动浏览条,为了方便看到关心的重要字段，需要调整字段显示顺序。 Raster Attribute Editor对话框菜单条：EditColumn Properties column properties对话框在Colu

11、mns中选择要调整显示顺序的字段，通过Up、 Down、Top、Bottom等几个按钮调整其合适的位置,通过选择Display Width调整其显示宽度，通过Alignment调整其对齐方式。如果选择Editable复选框，则可以在Title中修改各个字段的名字及其它内容。在Column Properties对话框中调整字段顺序，最后使Histogram、opacity、 color、 class_names四个手段的显示顺序依次排在前面。 点击OK按钮(关闭Column properties对话框） 返回Raster Attribute Editor对话框（lzisodat。img的属性表）

12、 第三步：给各个类别赋相应的颜色(如果在分类时选择了彩色，这一步就可以省去） Raster Attribute Editor对话框（lzisodat.img的属性）： 点击一个类别的Row字段从而选择该类别 右键点击该类别的Color字段（颜色显示区） As Is菜单 选择一种颜色 重复以上步骤直到给所有类别赋予合适的颜色 第四步：不透明度设置 由于分类图像覆盖在原图像上面，为了对单个类别的判别精度进行分析，首先要把其它所有类别的不透明程度(Opacity）值设为0（即改为透明）,而要分析的类别的透明度设为1（即不透明）。 Raster Attribute Editor 对话框(lzisoda

13、t。img的属性表）： 右键点击Opacity字段的名字 Column Options菜单Formula菜单项 Formula对话框 在Formula对话框的Formula输入框中(用鼠标点击右上数字区）输入0 点击Apply按钮（应用设置） 返回Raster Attribute Editor 对话框(lz-isodat.img的属性表）： 点击一个类别的ROW字段从而选择该类别 点击该类别的Opacity字段从而进入输入状态 在该类别的Opacity 字段中输入1，并按回车键 此时,在视窗中只有要分析类别的颜色显示在原图像的上面，其它类别都是透明的.第五步：确定类别专题意义及其准确程度 视窗

14、菜单条：Utilityflickerviewer Flicker对话框Auto Mode 本小步是设置分类图像在原图像修背景上闪烁,观察它与背景图像之间的关系从而断定该类别的专题意义，并分析其分类准确与否。 第六步：标注类别的名称和相应颜色 Raster Attribute Editor对话框（lzisodat.img的属性表)： 点击刚才分析类别的ROW字段从而选择该类别 点击该类别的class Names字段从而进入输入状态 在该类别的Class Names字段中输入其专题意义(如居民区）,并按回车键 右键点击该类别的Color字段（颜色显示区) As Is菜单选择一种合适的颜色 重复以上

15、4、5、6三步直到对所有类别都进行了分析与处理.注意，在进行分类叠加分析时,一次可以选择一个类别,也可以选择多个类别同时进行。Erdas基础教程: 数据预处理来源：陕西师范大学旅游与环境学院在ERDAS中，数据预处理模块为Data preparation。图标面板工具条中，点击图标Data Preparation菜单1、 图象几何校正 第一步:显示图象文件 在视窗中打开需要校正的Landsat TM图象：lanzhoucity。img。 第二步：启动几何校正模块 在Viewer#1的菜单条中,选择RasterGeometric Correction。打开Set Geometric Model对

16、话框 。选择多项式几何校正模型 Polynomial-OK .程序自动打开Geo Correction Tools对话框和 Polynomial Model Properties对话框。先选择Close关闭Polynomial Model Properties对话框 .程序自动打开GCP Tool Reference Setup对话框。选择Keyboard Only .OK 。程序自动打开 Reference Map Information提示框。选择Map Units： Meters .添加地图投影参数，如下图：。选择OK 确定地图投影参数，并关闭上图。 。选择OK，确定 Reference

17、 Map Information，并关闭提示框。 。并自动打开采集控制点对话框.GCP的具体采集过程： 在图像几何校正过程中，采集控制点是一项非常重要和相当重要的工作，具体过程如下： 。在GCP工具对话框中点select GCP 图标,进行GCP选择状态. 。在view1中移动关联方框位置，寻找明显地物特征点，作为输入GCP. .在GCP工具对话框中点击Great GCP图标，并在view2中点击左键定点，GCP数据表将记录一个输入的GCP，包括编号、标识码、X、Y坐标。 。在GCP工具对话框中输入地图参数坐标X、Y。 。不断重复上述步骤,采集若干GCP，直到满足所选是的几何校正模型为止。采集

18、地面检查点 以上所采集的GCP为控制点，用于建立转换方模型及 多项式方程，地面检查点，则用于检验所建立的转换方程的精度和实用性,具体过程如下： 。在GCP TOOL菜单条中选择GCP类型：Edit/Set Point Type -check。 .在GCP TOOL 菜单条中确定GCP匹配参数:Edit/Point matching 打开GCP Matching 对话框，并确定参数。 .确定地面检查点，其操作与选择控制点完全一样。 .计算检查点误差：在GCP TOOL工具条中，点击Compute Error 图标，检查点的误差就会显示在GCP TOOL 的上方（如下图)，只有所有检查的误差小于一

19、个像元时,才能进行以下的步骤.。在Geo Correction Tools 对话框中选择Model Properties 图标打开选择或检查参数，然后选择close关闭。图像重采样 。在Geo Correction Tools 对话框中选择Image Resample 图标-打开Image Resample 对话框，并定义重采样参数。.输出图像文件名（output file)：rectify。img 。选择重采样方法(Resample Method):Nearest Neighbor 。定义输出图像在图与像元大小. 。设查输出统计中忽略零值. 。选择OK启动重采样进程，并关闭Image Res

20、ample 对话框.2.图象拼接处理 本练习将同一区域机邻的三幅TM图象进行拼接。 其过程如下： .启动图象拼接工具,在ERDAS图标面板工具条中,点击Dataprep/Data preparation/Mosaicc lmages-打开Mosaic Tool 视窗。加载Mosaic图像，在Mosaic Tool视窗菜单条中，Edit/Add images打开Add Images for Mosaic 对话框。依次加载窗拼接的图像。在Mosaic Tool 视窗工具条中,点击set Input Mode 图标，进入设置图象模式的状态，利用所提供的编辑工具，进行图象叠置组合调查。 。图象匹配设置

21、,点击Edit /Image Matching -打击Matching options 对话框，设置匹配方法:Overlap Areas。在Mosaic Tool视窗菜单条中，点击Edit/set Overlap Function打开set Overlap Function对话框设置以下参数： .设置相交关系(Intersection Method）：No Cutline Exists。 。设置重叠图像元灰度计算(select Function）：Average。 .Apply -close完成。 .运行Mosaic 工具 在Mosaic Tool视窗菜单条中，点击 Process/Run M

22、osaic 打开Run Mosaic对话框。设置下列参数：确定输出文件名:mosaic。img 确定输出图像区域:ALL OK进行图像拼接.3、图象分幅裁剪 在实际工作中,经常根据研究区的工作范围进行图像分幅裁剪，利用ERDAS 可实现两种图像分幅裁剪:规则分幅裁剪,不规则分幅裁剪. 1、 规则分幅裁剪 即裁剪的边界范围为一矩形，其具体方法如下: 在ERDAS图标面板工具条中，点击DataPrep/Data preparation/subset Image打开subset Image 对话框，并设置参数如下：说明：裁剪范围输入: 通过直接输入右上角、右下角的坐标值； 先在图像视窗中放置查询框,

23、然后在对话框中选择From Inquire Box 先在图像视窗中绘制AOL区域，然后在对话框中选择AIO功能，利用此方法也可实现不规则裁剪。2、不规则分幅裁剪 用AOI区域裁剪，与上述的的方法相同。用Arclnfo的多边形裁剪。 。将Arclnfo多边形转换成网格图象. .ERDAS图标面板工具条中，点击Vector/vector to Raster-打开vector to Raster对话框，并设置参数，并实现转换。 。通过掩膜运算实现图像不规则裁剪。 。ERDAS 图标面板工具条中,点击Interpreter/Utilities /Mask-打开Mask对话框，并设置参数如下： setu

24、p Recode设置裁剪区域内新值为1，区域外取0值。Erdas基础教程： 数据输入来源：陕西师范大学旅游与环境学院 1。 单波段二进制图像数据输入 在ERDAS图标面板工具条中，点击-打开输入输出对话框，如图所示。并做如下的选择： .选择数据输入操作：Import .选择数据输入类型(Type）为普通二进制：Generic Binary 。选择数据输入媒体(Media）为文件：File 。确定输入文件路径及文件名（Input File）：Band1.dat 。确定输出文件路径及文件名（Output File）band1.dat 。OK。打开Import Generic Binary Data

25、对话框，如图所示在Import Generic Binary Data对话框中定义下列参数: .数据格式（Data Format）：BSQ 。数据类型(Data Type）:Unsigned 8 Bit 。数据文件行数(Row）:5728 .数据文件列数(Cols）:6920 。文件波段数量(Bands）:1 .保存参数设置（Save Option）:。gen .OK退出Save Option File .OK执行输入操作。上图为进程状态条，OK完成数据输入。 重复上述过程，可依此将多波段数据全部输入，转换为 .IMG文件。 2、组合多波段数据 为了图像处理与分析，需要将上述转换的单波段IMG

26、文件组合为一个多波段图像文件. 第一步：在ERDAS图标面板工具条中,点击InterpreterUtilitiesLayer Stack Layer Selection and Stacking的对话框。第二步：在Layer Selection and Stacking对话框中,依此选择并加载（Add）单波段IMG图像: 。输入单波段图像文件（Input File： *。img）:band1.img-Add 。输入单波段图像文件(Input File: *.img）:band2。imgAdd .输入组合多波段图像文件(Output File:.img）：bandstack。img 。OK执行并

27、完成波段组合。Erdas基础教程:影像文件的打开和显示来源：陕西师范大学旅游与环境学院1、图像显示视窗（Viewer） 图像显示视窗（Viewer）是显示栅格图像、矢量图形、注记文件、AOI等数据层的主要窗口，每次启动 ERDAS IMAGING时，系统都会自动打开一个二维视窗（Viewer）如图2.1所示.在应用过程中可以随时打开新的视窗。 图2。1 二维视窗（Viewer）二维视窗(Viewer）主要由视窗菜单条、工具条、显示窗和状态条四部分组成。2、图像显示 第一步：启动程序 菜单上选择File | Open Raster LayerSelect Layer To Add对话框图2。2。

28、 或在工具条上选择Select Layer To Add对话框图2.2。 图2.2 Select Layer To Add对话框第二步：确定文件 图2。2中的File选项用于图像文件的确定，具体内容及实例如表所示.表2。1 图像文件的确定参数第三步：设置参数 在Select Layer To Add对话框中点击Raster Options, 就进入设置参数状态, 如图2.3所示。 图2.3 Select Layer To Add对话框各项参数具体内容及实例，如表2.2 所示。表2。2 图像文件显示参数第四步：打开图像 在Select Layer To Add对话框中，点击OK，打开所确定的图像

29、，视窗中显示该图像。图2.4视窗中显示该图像Erdas基础教程： 图象增强处理来源:陕西师范大学旅游与环境学院1图像解译功能简介（Introduction of Image Interpreter） 利用ERADS IMAGINE 进行图像增强主要采用ERADS IMAGINE的图像解译器（Image Interpreter)模块，该模块包含了50多个用于遥感图像处理的功能模块,这些功能模块在执行过程中都需要通过各种按键或对话框定义参数,多数功能都借助模型生成器(Model Maker）建立了图形模型算法，容易调用或编辑。 图像解译器(Image Interpreter或Interpreter

30、),可以通过两种途径启动： ERDAS图标面板菜单条： Main/Image Interpreter-Image Interpreter菜单 ERDAS图标面板工具条:点击Interpreter图标一Image Interpreter菜单从上图可以看出，ERDAS图像解译模块包含了8个方面的功能，依次是遥感图像的空间增强(Spatial Enhancement）、辐射增强(Radiometric Enhancement）、光谱增强(SpectralEnhancement）、高光谱工具（Hyper Spectral Tools）、傅立叶交换（Fourier Analysis ）、地形分析（Top

31、ographic Analysis）。地理信息系统分析（GlS Analysis）、以及其它实用功能（Utilities）.每一项功能菜单中又包含若干具体的遥感图像处理功能. 2、图象空间增强 21 空间增强(Spatial Enhancement） 空间增强技术是利用像元自身及其周围像元的灰度值进行运算，达到增强整个图像之目的。ERDAS IMAGINE 提供的空间增强处理功能如表4.1所列.表4。1 遥感图像空间增强命令及其功能2.2卷积增强处理 卷积增强是将整个图像按像元分块进行平均处理，用于改变图像的空间频率特征.卷积处 理的关键是卷积算子一系数矩阵的选择.ERDAS将常用的卷积积算子放在default。klb 的文件中,分别以33、55、77三组，每组又包括edg