实验十监督分类实例与练习:某地区的遥感影像监督分类背景:现有某地区TM影像,需要进行土地利用类型划分,并对分类结果进行评价。
本例中使用监督分类的方法,实现土地利用类型的划分。
目的:通过练习,熟练掌握并理解监督分类的方法、评价方法的原理及实现过程,体会在具体应用中的适应性。
监督分类流程图:监督分类过程1 定义分类模板(1)在viewer窗口中选择打开smtm.img,在raster options选择fit to frame(2)单击classifer|classification|signature editor,打开分类模板编辑器(signature editor)(3)在viewer窗口中打开raster|tools,打开raster工具面板(4)选择AOI多边形绘制按钮,进入多边形AOI绘制状态,在图像上选择深蓝色区域,绘制一个多边形AOI,在signature editor窗口,单击按钮,将多边形AOI区域加载到signature editor分类模板属性表中(5)在图像上继续选择深蓝色区域,再绘制一个多边形AOI,在signature editor窗口,单击按钮,将多边形AOI区域加载到signature editor分类模板属性表中(6)用同样的方法加载9个深蓝色多边形AOI(7)在分类模板属性表中,依次单击这些AOI的class#字段下的分类编号(按住shift键),并单击,将所选中的模板合并成一个新的模板,则生成一个新的模板(8)单击其signature name属性进入编辑状态,输入water,单击color属性,选择深蓝色(9)在signature editor菜单条,单击edit|delete,删除合并前的模板(10)在图像上继续选择多个砖红色区域AOI(farmland),赭色区域AOI(forest),绿色区域AOI(grass),浅蓝色区域AOI(resident)(11)单击signatureedit|file|save命令,打开保存对话框,确定是保存所有模板(all),并保存分类模板文件的目录2评价分类模板(1)分类预警评价①选中water类别②在signature editor窗口,选择某类或者某几类模板,单击view|image alarm命令,打开signature alarm对话框③选中indicate overlap复选框,设置同时属于两个及以上的像元叠加预警显示,点击色框设置为黄色④点击edit parallelepiped limit|limit|set,设置计算方法(method):minimum/maximum,并选择使用的模板:current(当前模板)⑤设置完成后,单击ok按钮,返回limits对话框,单击close按钮,返回signature alarm对话框,单击Ok按钮执行分类预警评价,形成预警掩膜,单击close按钮,关闭signature alarm对话框⑥在viewer窗口,选择utility|flick命令,选择auto mode,speed设为600.⑦按照同样的方法分别对forest, farmland, grass类进行预警掩膜,查看准确性(2)可能性评价①选中signature editior属性表中的所有类别②单击evaluate|contingency, 打开contingency matrix对话框,nonparamatric选择feature space, overlap rules选择parametric rule,unclassified rule选择paramatric rule,paramatric rule 选择maximum likehood③ 设置完成后,点击ok按钮,便显示分类误差矩阵,若误差矩阵大于85%,结果令人满意。
(3)分类的分离性①在模板编辑器中选择water、forest②选择signature editor|evaluate|separability命令,打开signature separability对话框③组合数据层数(layer per combination)选择3,distance meature选择transformed divergence,output form选择ASCII,report type选择complete report④点击ok完成设置,计算其分离性⑤一次选择其他类的计算分类的分离性。
若TD值均大于1700,说明可以分开3 进行监督分类(1)选择处理图像文件为smtm(2)在input signature中选择super.img(3)在classified file中设置输出数据存储路径及其名称,这里为super.img (4)选中输出分类距离文件为distance file(5)在non-parametric rule中选择feature space(6)在overlap rule中选择parametric rule(7)在unclassified rule中选择parametric rule(8)在parametric rule中选择maximum likehood(9)单击ok,执行监督分类4 结果评价(1)分类叠加结束监督分类后,在viewer窗口打开smtm.img和super.img,在打开super.img时raster options中去除clear display选项。
在viewer窗口,选择utility|flick命令,选择auto mode,speed设为600。
检查分类结果的准确性。
(2)分类评估第一步:在视窗中打开原始图像在Viewer中打开分类前的原始图像,以便进行精度评估。
第二步:启动精度评估对话框ERDAS 图标面板菜单条:Main→Image Classification →Classification或ERDAS 图标面板工具条:点击Classifier 图标→Classification菜单→选择Accuracy Assessment菜单项→打开Accuracy Assessment对话框Accuracy Assessment 对话框中显示了一个精度评估矩阵(Accuracy Assessment Cellarray)。
精度评估矩阵中将包含分类图像若干像元的几个参数和对应的参考像元的分类值。
这个矩阵值可以使用户对分类图像中的特定像元与作为参考的已知分类的像元进行比较,参考像元的分类值是用户自己输入的。
矩阵数据存在分类图像文件中。
第三步:打开分类专题图像Accuracy Assessment 对话框菜单条:File→Open→打开Classified Image对话框→在Classified Image对话框中确定与视窗中对应的分类专题图像→OK(关闭Classified Image对话框)→返回Accuracy Assessment对话框第四步:将原始图像视窗与精度评估视窗相连接Accuracy Assessment对话框:→工具条:点击Select Viewer图标(或菜单条:选择View 菜单的Select Viewer)→将光标在显示有原始图像的视窗中点击一下→原始图像视窗与精度评估视窗相连接第五步:在精度评价对话框中设置随机点的色彩Accuracy Assessment对话框:→菜单条View →Change Colors菜单项→打开Change color面板→在Points with no Reference确定没有真实参考值的点的颜色→在Points with Reference确定有真实参考值的点的颜色→OK(执行参数设置)→返回Accuracy Assessment对话框第六步:产生随机点本步操作将在分类图像中产生一些随机的点,随机点产生之后,需要用户给出随机点的实际类别。
然后,随机点的实际类别与在分类图像的类别将进行比较。
Accuracy Assessment对话框:→Edit →Create/Add Random Points→打开Add Random Points对话框→在search Count中输入1024→在Number of Points中输入20→在Distribution Parameters选择Random单选框→OK(按照参数设置产主随机点)→返回Accuracy Assessment对话框可以看到在Accuracy Assessment对话框的数据表中出现了20个比较点,每个点都有点号、X\Y坐标值、Class、Reference等字段,其中点号、X/Y坐标值字段是有属性值的。
说明:在Add Random Point对话框中,search Count是指确定随机点过程中使用的最多分析像元数,当然这个数目一般都比Number of Point大很多,Number of Points设为20说明是产生20个随机点,如果是做一个正式的分类评价,必须产生20个以上的随机点。
选择Random意味着将产主绝对随机的点位,而不使用任何强制性规则。
Equalized Random是指每个类将具有同等数目的比较点。
stratified Random是指点数与类别涉及的像元数成比例,但选择该复选框后可以确定一个最小点数(选择Use Minimum Points),以保证小类别也有足够的分析点。
第七步:显随机点及其类别Accuracy Assessment对话框:→View. →Show All(所有随机点均以第五步年设置的颜色显示在视窗中)→Edit→Show Class Values(各点的类别号出现在数据表的class字段中)第八步:输入参考点的实际类别值Accuracy Assessment对话框:→在数据表的Reference字段输入各个随机点的实际类别值(只不要输入参考点的实际分类值,它在视窗中的色彩就变为第五步设置的Point With Reference颜色)第九步:设置分类评价报告输出环境及输出分类评价报告Accuracy Assessment对话框:→Report →Options→通过点击确定分类评价报告的参数→Report→Accuracy Report(产生分类精度报告)→Report→Cell Report(报告有关产生随机点的设置及窗口环境)→所有报告将显示在ERDAS文本编辑器窗口,可以保存为本文件→File→Save Table(保存分类精度评价数据表)→File→close (关闭Accuracy Assessment对话框)通过对分类的评价,如果对分类精度满意,保存结果。
如果不满意,可以进一步做有关的修改,如修改分类模板等,或应用其它功能进行调整。