遥感图像解译实验报告.doc

实验一 非监督分类一、基本原理非监督分类的前提是假定遥感影像上同类地物在相同的条件下具有相同的光谱特征信息。非监督分类不需要先验类别知识的情况下，仅依靠影像本身的特征进行特征提取，根据统计特征及点群的分布情况来划分地物。ERDAS 的非监督分类是基于迭代自组织数据分析（ISODATA）算法实现的，算法的流程图如图1。结果输出确定类别中心像素聚类设定初始聚类中心开始图1 ISODATA算法流程图ISODATA算法实现步骤：（1）按照某个原则确定一些初始聚类中心。在实际操作时，要把初始聚类设定得大些，同时引入各种参数控制迭代的次数。（2）像素聚类与分析。计算像素与初始类别中心的距离，把像素分配到最近的类别中，从而获得每个初始类别的集群成员。（3）确定类别中心。在全部像素按照各类中心分类后，重新计算每一类新的均值，并作为下一次分类的中心，并执行上一步，循环迭代，直到达到迭代的次数。二、分类过程运用非监督分类对遥感图像进行分类的过程，主要分为以下两个步骤：1、动非监督分类模块，选择输入、输出影像（1）在ERDAS面板工具中选择Classifier-Classification-Unsupervised Classification，打开非监督分类对话框。（2）选择图像处理文件（Input Raster File）和输出文件（Output Cluster Layer Filename）设置被分类的图像和分类结果，并选择生成分类模板文件（Output Signature Set）产生一个模板文件。2、初始参数，执行非监督分类在非监督分类对话框中分别设置聚类选项（Clustering Options）和处理选项（Processing Options）：（1）选择Classifier-Unsupervised Classification命令，打开非监督分类对话框，设定输入输出数据，设置聚类选项（Clustering Options），确定初始聚类方法和分类数。（2）设置处理选项（Processing Options），确定循环次数和阈值。三、实验结果本次实验使用Initial from Statistics方法，分类数为6，其它参数取默认值。遥感图像原图如图2。图2 遥感图像原图非监督分类结果如图3生成的模板文件如图4图3 非监督分类结果图4 非监督分类模板文件对分类图像颜色进行调整后的分类结果如图5对应模板文件如图6图5 调整后的非监督分类结果图6 调整后的模板文件三、实验结果分析1、从图2、图3的对比，可以看出，经过非监督分类处理，水域、河流、道路、建筑区域、绿地这五大类别被分类出来的结果是很好的，但图像的亮度有些偏低。2、从图3、图4的对比，可以得出结论，经过颜色调整后的分类结果更清晰，细节更突出，颜色的表示更符合人们的习惯。3、非监督分类在没有经验的条件下进行，客观、实际、具体，但也面临着一些比较特殊的区域被分错误的问题。实验二 监督分类一、基本原理不同于非监督分类，监督分类基于先验知识，根据训练场提供的样本选择特征参数，建立判别函数，对待分类点进行分类。训练场地的选择是监督分类的关键，在监督分类中由于训练样本的不同，分类结果会出现极大的差异，因此应选取有代表性的样本，用于监督分类的训练样本应该是光谱特征比较均一的地区，在图像中根据均一的色调估计只有一类地物，且一类地物的训练样本可以选取一块以上。此外，用做样本的数目至少能满足建立分类用判别函数的要求，对于光谱特征变化比较大的地物，训练样本要足够得多，以反映其变化范围。一般情况下，要得到可靠的结果，每类至少选择10100个训练样本。监督分类的判别规则可以分为无变量和有变量两种。其中，无变量包括特征空间、平等六面体，有变量包括最大似然、最小距离和马氏距离。当然，非参数规则与参数规则可以同时使用，但要注意应用范围。另外，如果使用非参数型模板，还要确定叠加规则和未分类规则。二、实验过程1、定义分类模板定义分类模板操作包括分类模板的生成、管理、评价和编辑等，功能主要由分类模板编辑器（Signature Editor）完成，具体步骤包括：（1）需要分类的影像点击ERDAS-Viewer面板，打开需要进行分类的图像。（2）分类模板编辑器选择Classifier图标-Classification-Signature Editor命令，打开分类模板编辑器（Signature Editor）对话框。（3）属性字段在分类编辑器窗口中的分类属性表中有很多字段，可以对不需要的字段进行调整。选择分类编辑窗口的View-Column，打开View Signature Column对话框，选中需要显示的字段，单击Apply按钮，显示发生变化，单击Close按钮完成。（4）选取样本基于先验知识，需要对遥感图像选取训练样本，包括产生AOI、合并、命名，从而建立样本。考虑到同类地物颜色的差异，因此在采样过程中对每一地类的采样点不少于10个。（5）保存分类模板选择Signature Edit面板-File-Save命令，打开保存对话框，确定是保存所有模板或是只保存被选中的模板，并确定保存分类模板文件的目录和文件名，单击OK按钮即可。2、监督分类建立满意的分类模板后，就需要在一定的分类决策规则条件下，对像元进行聚类判断。在选择判别函数及相应的准则后，便可执行监督分类，具体步骤如下：选择Classifier图标-Supervised Classification命令，打开Supervised Classification对话框，参数设置如下：（1） 选择处理图像文件（Input Raster File）。（2） 确定输入分类模板（Input Signature）。（3） 定义输出分类文件（Classified File）。（4） 设置输出分类距离文件为Distance File。（5） 选择非参数规则（Non-Parametric Rule），一般选择Feature Space，即选择特征空间。（6） 选择叠加规则（Overlap Rule），一般为Parametric Rule。（7） 选择未分类规则（Unclassified Rule）为Parametric Rule。（8） 选择参数规则（Parametric Rule），一般选择Maximum Likelihood，即最大似然。（9） 还可以定义分类图的属性表项目，即单击Attribute Options按钮，进行选择。（10） 最后单击OK按钮，执行监督分类。三、实验结果根据上述实验过程，可得分类后结果如图8，其中分类模板如图9。图7 遥感影像原图图8 监督分类结果图9 监督分类模板四、实验结果分析1、由图8可以看出，河流、水域、绿色植物、道路、建筑物被分类出来，但在水域中出现了河流的错误分类，这是不完美的。2、道路一地物被分的很乱，这可能与邻近道路的像素与道路像素相近的原因引起的，这种问题可以通过增加分类各类的办法，将邻近道路的像素分成另一类的方法来解决。3、监督分类效果在一定程度上结合了人们所掌握的经验知识，但还是不能完美的将人们所掌握的知识与计算机机械的分类相结合，这是未来智能解译很面临的主要问题之一。实验三 分类后处理一、基本原理由于分类严格按照数学规则进行，分类后往往会产生一些只有几个像元甚至一两个像元的小图斑。这对分类图的分析、解译和制图，都是不利的，可通过几种分类后处理来解决。ERDAS中的分类后处理方法有：聚类统计、过滤分析、去除分析。1、 聚类统计聚类统计（Clump）是通过计算分类专题图像每个分类图斑的面积、记录相邻区域中最大图斑面积的分类值等操作，产生一个Clump类组输出图像，其中每个图斑都包含Clump类组属性。该图像是一个中间文件，用于进行下一步处理。2、 过滤分析Sieve功能是对经Clump处理后的Clump类组图像进行处理，按照定义的数值大小，删除Clump图像中较小的类组图斑，并给所有小图斑赋予新的属性值。显然，这里引出了一个新的问题，就是小图斑的归属问题。可以与原分类图对比确定其新属性，也可以通过空间建模方法、调用Delerows或Zonel工具进行处理。Sieve经常与Clump命令配合使用，对于无须考虑小图斑归属的应用问题，有很好的作用。3、 去除分析去除分析是用于删除原始分类图像中的小图斑或Clump聚类图像中的小Clump类组，与Sieve命令不同，将删除的小图斑合并到相邻的最大的分类当中。而且，如果输入图像是Clump聚类图像的话，经过Eliminate处理后，将小类图斑的属性值自动恢复为Clump处理前的原始分类编码。显然，Eliminate处理后的输出图像是对分类结果图像进行了制图综合。二、实验过程1、 聚类分析在ERDAS图标面板菜单条中，选择Main-Image Interpreter-GIS Analysis-Clump命令，打开Clump对话框，如图10，并设置下列参数。图10 Clump对话框（1） 选择处理图像文件（Input File）。（2） 定义输出文件（Output File）。（3） 选择文件坐标类型（Coordinate Type）为Map/File。（4） 处理范围确定（Subset Definition）：UL X/Y，LR X/Y。（5） 确定聚类统计邻域大小（Connected Neighbors）：统计分析将对每个像元四周的N个相邻像元进行。可以选择4个方向或者8个相邻的像元。（6） 点击OK按钮，关闭Clump对话框，执行聚类统计分析。2、 过滤分析在ERDAS图标面板菜单条中，选择Main-Image Interpreter-GIS Analysis-Sieve命令，打开Sieve对话框，如图11，确定下列参数：（1） 选择处理图像文件（Input File）。（2） 定义输出文件（Output File）。（3） 定义输出文件坐标类型（Coordinate Type）：Map/File。（4） 定义输出文件处理范围（Subset Definition）：UL X/Y，LR X/Y。（5） 选择处理的图层（Select Layer）。（6） 忽略输出统计零像元值Ignore Zero in Output Statistics。（7） 确定最小图斑大小（Minimum Size）。（8） 点击OK按钮，关闭Sieve对话框，执行过滤分析。图11 Sieve对话框3、 去除分析在ERDAS图标面板菜单条中，选择Main-Image Interpreter-GIS Analysis-Eliminate命令，打开Eliminate对话框，如图12，确定下列参数：（1） 选择处理图像文件（Input File）。（2） 定义输出文件（Output File）。（3） 定义输出文件坐标类型（Coordinate Type）：Map/File。（4） 定义输出文件处理范围（Subset Definition）：UL X/Y，LR X/Y。（5） 选择去除分析的图层（Select Layer）。（6） 忽略输出统计零像元值Ignore Zero in Output Statistics。（7） 确定最小图斑大小（Minimum Size）。（8） 点击OK按钮，关闭Elininate对话框，执行去除分析。图12 Eliminate对话框三、实验结果1、聚类统计处理结果如图13所示。图13 聚类统计处理结果2、过滤分析处理结果如图14所示。图14 过滤分析处理结果4、 去除分析处理结果如图15、图16所示。图15 去除分析处理结果图16 颜色调整后的去除分析处理结果四、实验结果分析1、经过聚类统计处理后，过滤分析过程完成较小的类组图斑的删除、归属过程，从监督分类后的结果图与过滤分析结果图对比可以看出，原来图像中较小的图斑消失，被划归到相邻大的图斑中去了，这