Erdas非监督分类全过程.doc

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

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

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

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

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

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

7、指由图像文件整体（或其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 Classification

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

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

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

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

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

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

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

16、多项式几何校正模型 PolynomialOK .程序自动打开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 Map I

17、nformation，并关闭提示框。 .并自动打开采集控制点对话框。GCP的具体采集过程: 在图像几何校正过程中，采集控制点是一项非常重要和相当重要的工作，具体过程如下： .在GCP工具对话框中点select GCP 图标，进行GCP选择状态。 .在view#1中移动关联方框位置，寻找明显地物特征点，作为输入GCP。 .在GCP工具对话框中点击Great GCP图标，并在view#2中点击左键定点，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 Resample

20、对话框。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 图标，进入设置图象模式的状态，利用所提供的编辑工具，进行图象叠置组合调查。 .图象匹配设置，点击Ed

21、it /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 Mosaic

22、打开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对话框，并设置参数如下： setup Recod

24、e设置裁剪区域内新值为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图标面板工具条中，点击Interpreter|Utilities|Layer Stack Layer Selection and Stacking的对话框。第二步：在Layer Selection and Stacking对话框中，依此选择并加载（Add）单波段IMG图像： .输入单波段图像文件（Input File: *.img）:band1.imgAdd .输入单波段图像文件（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 ）、地形分析（Topogr

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

32、default.klb 的文件中，分别以33、55、77三组，每组又包括edge Detect/edge enhance/low pass/Highpass/Horizontal/vertical/summary等七种不同的处理方式。 具体执行过程如下： .ERDAS图标面板工具条，点击Interpreter/spatial Enhancement / convolution打开convolution对话框，并设置如下参数：2.3纹理分析 纹理分析通过在一定的窗口内进行二次变异分析或三次对称分析，使图像的纹理结构得 到增强，具体过程如下： ERDAS 图板面板工具条中，点击Interprete

33、r/spatial Enhancement/Texture打开Texture对话框，并设置参数如下：这一分析方法的关键是确定Window size 的定义操作函数Operator。3、辐射增强处理： 3.1辐射增强（Radiometric Enhancement） 辐射增强处理是对单个像元的灰度值进行变换达到图像增强的目的。 ERDAS IMAGINE供的辐射增强处理功能如表4.2所列。 表4.2 遥感图像辐射增强命令及其功能3.2查找表拉抻处理 是通过修改图像查找表，使输出图像值发生变化。通过定义，可实现线性拉伸、分段线性拉伸、非线性拉伸等处理。 .在ERDAD图标面标工具条中，点出Inte

34、rpreter/Radiometric Enhancement/LUT Stretch-打开LUT stretch对话框，并设置参数如下：其中关建是：点击View/custom Table 进入查找表编辑状态。根据需要修改查找表。 3.3直方图均衡化处理 该处理实质上是对图像进行非线性拉伸，重新分配图像像元值，使一定灰度范围内的像元的数量大致相等。 .在ERDAS图标面板工具条中，点击Interpreter/Radiometric Enhancement/Histogram Equalization-打开Histogram Equalization对话框,并设置参数如下：4.光谱增强处理 4.

35、1光谱增强处理是基于多波段数据对每个像元的灰度值进行变换，达到图像增强的目的。ERDSA IMAGINE 提供的光谱增强处理功能如下图所示 表4.3 遥感图像光谱增强命令及其功能4.2主成份变换 是一种常用的数据压缩方法，它可以将具有相关性的多波段数据压缩到完全独立的较少的几个波段上，便于图像分析与解译，具体过程如下： .在ERDAS 图标面标工具条中，点击Interpreter/spectral Enhancement/principal Comp-打开principal components对话框,并设置参数如下：上图为原始图像与主成份变换获得的主份量1、2、3的对比图原图像432(RGB

36、) 主成份432（RGB）合成图像4.3 缨穗变换 缨穗变换是针对植物学所关心的植被特征，对原始多波段图像数据进行空间旋转，获得具有物理意义的亮度、绿度、湿度等分量。 .在ERDAS图标面板工具条中，点击Interpret/spectral Enhancement/Tasseled Cap打开Tasseled cap对话框，并设置参数。 注，需要定义相关参数set coefficients,点击set coefficients按钮，打开对话框，确定区域内的类型。下图为原始图像与缨穗变换的亮度、绿度、湿度分量的比较图。4.4色彩变换 色彩变换是将区域图像从RGB的彩色空间转换到IHS作为定位参数

37、的彩色空间，以便达到增强目的。 .在ERDAS图标面板工具条中，点击Interpreter/spectral Enhancement/RGB to IHS-打开RGB TO IHS对话框，并设置参数如下：上图为原始图像与色彩变换获得的亮度、色度、饱和度三个分量图的对比Erdas非督分类步骤一：分类模板定义来源：陕西师范大学旅游与环境学院定义分类模板（Define Signature Using signature Editor） ERDAS IMAGINE的监督分类是基于分类模板来进行的，而分类模板的生成、管理、评价、和编辑等功能是由分类模板编辑器来负责的。毫无疑问，分类模板生成器是进行监督分

38、类一个不可缺少的组件。 在分类模板生成器中，生成分类模板的基础是原图像和（或）其特征空间图像。因此，显示这两种图像的视窗也是进行监督分类的重要组件。 第一步：显示需要进行分类的图像 在视窗中显示ERDASHOMEexeciseljxtm.img （ Red4Grean5B1ue3、选择Fit to Frame，其它使用缺省设置）。 第二步：打开模板编辑器并调整显示字段 ERDAS图标面板工具子，点击C1assifier目标 Classification菜单 Signature Editor菜单项 Signature Editor对话框从上图中可以看到有很多字段，有些字段对分类的意义下大，我们希

39、望不显示这些这段，所以要进行如下调整： Signature Edit对话框菜单条： ViewColumnsview signature columns对话框点击最上一个字段的Co1unmn字段下拖拉直到最后一个段,此时，所有字段都被选择上，并用黄色（缺省色）标识出来。 按住shift键的同时分别点击Red、 Green、 B1ue三个字段Red、Green、 Blue三个字段将从选择集中被清除。 点击Apply按钮 点击Close按钮 从View Signature Co1umns对话框可以看到Red、Green、Blue三个字段将不再显示。 Erdas非督分类步骤二：获取分类模板信息（一）来

40、源：陕西师范大学旅游与环境学院可以分别应用AOI绘图工具、 AOI扩展工具、查询光标等三种方法，在原始图像或特征空间图像中获取分类模板信息。但在实际工作中也许只用一种方法就可以了。 （1）应用AOI绘图工具在原始图像获取分类模板信息 无论是在原图像还是在下面要讲的特征空间图像中，都是产主AOI区域来作为分类模板信息的来源 。首先练习如何用AOI绘图工具获取分类模板信息。 在显示有ljxtm.img图像的视窗： 点击图标（或者选择Raster菜单项选择Tools菜单） 打开Raster工具面板 点击Raster工具面板的图标 在视窗中选择红色区域，绘制一个多边形AOI 在Signature Ed

41、itr对话框，点击图标，将多边形AOI区域加载到Signature分类模板中 在Signature Editor中，改变刚才加入模板的Signature Name和Color。 重复上述操作过程以多选择几个红色区域AOI,并将其作为新的模板加入到Signature Editor当中，同时确定各类的名字及颜色。如果对同一个专题类型（如水体）采集了多个AOI并分别生成了模板，可以将这些模板合并，以便该分类模板具多区域的综合特性。具体做法是在Signature Editor对话框中，将该类的Signature 全部选定，然后点击合并图标，这时一个综合的新模板生成，原来的多个Signature 同时存

42、在（如果必要也可以删除）。Erdas非督分类步骤二：AOI扩展工具获取分类模板信息来源：陕西师范大学旅游与环境学院可以分别应用AOI绘图工具、 AOI扩展工具、查询光标等三种方法，在原始图像或特征空间图像中获取分类模板信息。但在实际工作中也许只用一种方法就可以了。 （2）应用AOI扩展工具在原始图像获取分类模板信息 扩展生成AOI的起点是一个种子像元。与该像元相邻的像元被按照各种约束条件来考察，如空间距离、光谱距离等。如果被接受，则与原种子一起成为新的种子像元组，并重新计算新的种子像元平均值（当然也可以设置为一直沿用原始种子的值。以后的相邻像元将以新的平均值来计算光谱距离。但空间距离一直是以最

43、早的种子像元来计算的。 应用AOI扩展工具在原始图像获取分类模板信息，首先必须设置种子像元特性，过程如下： 在显示有ljxtm.img图像的视窗中： AOI一Seed Properties菜单 Region Growing Properties对话框在Neighborhood选择按四个相邻像元扩展，表示被点击像元的上、下、左、右四个像元与被点击像元是相邻的。而表示其周围9个像元都与被点击像元相邻。这里选择。 在Geographic Constrains设置地理约束, Area确定每个AOI所包含的最多像元数（或者面积），而Distance确定AOI所包含像元距被点击像元的最大距离，这两个约束可

44、以只设置一个，也可以设置两个或者一个也不设。在此处只设置面积约束为300个像元。 在Spectral Euclidean Distance中设置波谱欧氏距离,本约束是指AOI可接受的像元值与种子像元平均值之间的最大光波欧氏距离（两个像元在各个波段数值之差的平方之和的二次根），大于该距离将不被接受。此处设置距离为：10 点击Options按钮，打开Region Grow Options面板以确定一些扩展设置 Region Grow Options面板上有三个复选框。在种子扩展的过程中可能会有些不符合条件的像元被符合条件的像元包围，选择Include Island Po1ygons使这些不符合条件

45、像元,将以岛的形式被删除出来，如果不选择则全部作为AOI的一部分。Update Region Mean是指每一次扩展后是否重新计算种子的平均值，如果选择该复选框则重新计算，如果不选择则一直以原始种子的值为平均值。Buffer Region Boundary复选框是指对AOI产生缓冲区，该设置在选择AOI编辑DEM数据时比较有用，可以避免高程的突然变化。这里选择Include Island Polygons和Update Region Mean。 至此完成了种子扩展特性的设置，下面将使用种子扩展工具产生一个AO1。 在显示有ljxtm.img图像的视窗中： 在视窗工具条中点击图标（或在视窗菜单条

46、： RasterTools） 打开Raster工具面板 点击Raster工具面板的图标点击视窗中的绿色区域绿色区域对应的是耕地，AOI自动扩展将生成一个针对耕地的AO1。如果扩展AOI不符合需要。可以修改Region Growing Properties直到满意为止，注意在Region Growing Properties对话框中修改设置之后，直接点击Redo按钮就可重新对刚才点击的像元生成新的扩展AO1。 在signature editor对话框，点击图标，将扩展AOI区域加载到signature分类模板中 在Signature Editor中，改变刚才加入模板的Signann1e的名字（N

47、ame ）和颜色（Color）。 重复上述操作步骤，选择多AOI区域，并将其作为新的模板加入到Signature Editor中，同时确定各类别的名字及颜色。 Erdas非督分类步骤二：查询光标扩展方法获取分类模板信息来源：陕西师范大学旅游与环境学院可以分别应用AOI绘图工具、 AOI扩展工具、查询光标等三种方法，在原始图像或特征空间图像中获取分类模板信息。但在实际工作中也许只用一种方法就可以了。 应用查询光标扩展方法获取分类模板信息 该方法与AOI扩展工具法大同小异，只不过AOI扩展工具法是在选择扩展工具后，用点击的方式在图像上确定种子像元，而本方法是要用查询光标（Inquire Curso

48、r）确定种子像元。种子扩展的设置与第AOI扩展工具法种方法完全相同。 在显示有ljxtm.img图像的视窗中点击： Utility一Inquire Cursor 在视窗中出现一个十字光标，十字交点可以准确定位一个像元的位置 将十字光标标交点移动到种子像元上 点击Region Growing Properties 对话框的Grow at Inquire按钮产生一个新的AOI 在Signature Editor对话框，点击图标，将AOI 区域加载到Signature分类模板中第四步：保存分类模板 以上分别用不同方法产生了分类模板，下面将该模板保存起来。 在Signature Editor对话框菜单

49、条： FileSave 打开Save Signature File As对话框 确定是保存所有的模板还是只保存被选中的模板 确定文件的目录和名字（Sjg文件） 点击OK按钮 Erdas非督分类步骤三：分类报警工具评价分类模板来源：陕西师范大学旅游与环境学院分类模板建立之后，就可以对其进行评价，删除、更名、与其它分类模板合并等操作。分类模扳的合并可使用户应用来自不同训练方法的分类模板进行综合复杂分类，这些模板训练方法包括监督、非监督、参数化和非参数化。 分类模板评价工具包括： Alarms：分类报警工具 Contingency matrix：可能性矩阵 Feature objects：特征对象

50、Feature Space to image masking：特征空间到图像掩模 Histograms：直方图方法 Signature separability: 分类的分离性 Statistics:分类统计分析 当然，不同的评价方法各有不同的应用范围。例如不能用Separability工具对非参数化（由特征空间产生）分类模板进行评价，而且分类模板中至少应具有5个以上的类别。 或多个类别进行。如果没有在Signature Editor中选择类别，那么当前活动类别(Signature Editor中“”符号旁边的类别)就被用于进行报警。具体使用过程如下： 在Signature Editor对话框

51、：View Image Alarm 打开Signature Alarm对话框 选中Indicate Overlap 点击Edit Parallelepiped Limits按钮 Limits对话框 点击SET按钮 打开Set Parallelepiped Limits对话框 设置计算方法（Method）:Minimum/Maximum 选择使用的模板（Signature）:Current OK（关闭set Parallelepiped Limits对话框） 返回Limits对话框 Close（关闭Limits对话框） 返回Signature Alarm对话框 OK（执行报警评价，形成报警掩膜）

52、 Close（关闭signature Alarm对话框） 根据Signature Editor中指定的颜色，选定类别的像元显示在原始图像视窗中，并覆盖在原图像之上， 形成一个报警掩膜。第二步，利用Flicker功能查看报警掩膜 第三步：删除分类报警掩膜 视窗菜单条： ViewArrange Layers菜单 打开Arrange Layers对话框 右键点击Alarm Mask图层 弹出Layer Options菜单选择Delete Layer Alarm Mask图层被删除 App1y （应用图层删除操作） 提示verify Save on C1ose N0 Close （关闭Arrange

53、Layers对话框）Erdas非督分类步骤三：特征空间到图像掩模法评价分类模板来源：陕西师范大学旅游与环境学院分类模板建立之后，就可以对其进行评价，删除、更名、与其它分类模板合并等操作。分类模扳的合并可使用户应用来自不同训练方法的分类模板进行综合复杂分类，这些模板训练方法包括监督、非监督、参数化和非参数化。 分类模板评价工具包括： Alarms：分类报警工具 Contingency matrix：可能性矩阵 Feature objects：特征对象 Feature Space to image masking：特征空间到图像掩模 Histograms：直方图方法 Signature separ

54、ability: 分类的分离性 Statistics:分类统计分析 由特征空间模板产生图像掩膜 只有产生于特征空间Signature才可使用本工具，使用时可以基于一个或者多个特征空间模板。如果没有选择集，则当前处于活动状态（位于“”符号旁边）的模板将被使用。如果特征空间模板被定义为一个掩膜，则图像文件会对该掩膜下的像元作标记，这些像元在视窗中也将被显示表达出来（Highlighted）。因此可以直观地知道哪些像元将被分在特征空间模板所确定的类型之中。必须注意，在本工具使用过程中视窗中的图像必须与特征空间图像相对应。 下面是本工具的使用过程： signature Editor对话框。 选择要分析

55、的特征空间模板 FeatureMaskingFeature Space to Image 打开FS to Image Masking对话框不选择Indicate overlay复选框（选择Indicate Over1ay复选框意味着“属于不只一个特征空间模板的像元”将用该复选框后边的颜色显示） Apply （应用参数设置，产主分类掩膜） Close （关闭FS to Image Masking对话框） 模板对象图示 模板对象图示工具可以显示各个类别模板（无论是参数型辽是非参数型）的统计图，以便比较不同的类别，统计图以椭圆形式显示在特征空间图像中：每个椭园都是基于类别的平均值及其标准差。可以同时产主一个类别或多个类别的图形显示。如果没有在模板编辑器中选择类别。那么当前处于活动状态（位于“”符号旁边）的类别就被应用，模板对象图示工具还可以同时显示两个波段类别均值、平行六面体和标识（Label）。由于是在特征空间图像中绘画椭圆，所以特征空间图像必须处于打开状态。 Signature Editor对话框： 菜单条： Feature Objects 打开Signature objects对话框 确定特征空间图像视窗（Viewer）: 2 确定绘制分类统计椭圆：选择Plot Ellipses选框 确定统计标准差（Std.dev.）4 OK（执行模板对象图示，绘制分