土地变化信息提取.doc

《土地变化信息提取.doc》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。

1、遥感实验报告课程名称：遥感概论年级：实验日期：姓名：学号：班级：实验名称：土地变化信息提取实验序号：10成员人数：实验目的：1、了解常用提取土地变化信息的方法，掌握分类后比较方法提取土地利用变化信息。2、通过实习内容掌握土地利用变化信息的基本流程（变化信息发现----变化信息提取----变化信息类型表示）。实验内容：1、已学过的监督分类方法对两时相影像进行土地利用类型分类，提取土地利用变化信息2、完成一份实验报告实验方法和步骤：（1）观察图像，建立土地覆被分类体系，利用ERDASAOI画训练样区；（

2、2）选择适当方法进行监督分类，及后处理，分别获land1990land2006两个时相的土地覆被分类信息。(3)对两个时相的土地覆被分类结果进行精度检验，如果精度不合格，重新进行分类，检验，直到精度满足为止。(4)打开ERDAS中的图像解译（Interpreter）-GIS分析工具-matrix，，提取变化信息，导入到ARCGIS进行专题制图。（5）计算土地覆被变化矩阵。实验的过程和结果：（一）分别对两张图像进行监督分类1.第一步：定义分类模板第一步：打开分类的图像，南宁市1990年9月16日TM

3、，目视判断该遥感图像中南宁市土地利用类型，确定土地利用分类体系为：耕地、灌草地、林地、水域、建设用地、裸地。第二步：打开模板编辑器并调整显示字段点击主菜单上的classifier打开classification对话框,选择signatureeditor。如图1-1：图1-1第三步：获取分类模板信息，点击AOI，利用AOI-tools中的多边形工具绘制某一地类的样区。将画好的耕地AOI添加到模板。signatureeditor-edit-add.重复步骤第三步，在图中采集多个耕地样本。选择所有耕地样本

4、模板，按merge按纽合并这组分类模板。合并后将模板取名为耕地。利用同样的方法，依次做好其灌草地、林地、水域、建设用地、裸地土地覆盖类型模板。如图1-2：1990年图像的分类信息模板：图1-22006年分类信息模板：图1-3第四步：保存分类模板。2.评价分类模板第一步：点signatureeditor-Evaluate-contingency，利用可能性矩阵方法评价分类模板精度。达到90以上即为精度满足要求，否则重新选择训练样区，再次进行精度评价，直到精度满足。1990年的模板评价结果：Classi

5、fiedDataforestandwaterforesturban--------------------------------------------------forestand99.830.001.141.66water0.00100.000.000.00forest0.000.0090.900.00urban0.000.000.0098.34bareland0.170.000.000.00formland0.000.008.960.00ColumnTotal604638703602Refe

6、renceData--------------ClassifiedDatabarelandformlandRowTotal----------------------------------------forestand0.000.00621water0.000.00638forest0.004.88662urban0.000.00592bareland100.000.00612formland0.0095.12648ColumnTotal6116153773-----EndofErrorMatri

7、x-----2006年影像的模板评价结果：3.执行监督分类打开Classifier-surperviseclassification,选择分类决策规则，进行监督分类。1990年的结果如图1-4：图1-42006年的结果如图1-5：图1-5（一）分类后处理1.去噪：去噪声，interpreter-gisanalysis，然后在弹出的窗口中选择neighborhood,1990年结果如图2-1：图2-12006年的结果如图2-2：图2-22.分类重编码分类重编码主要是针对之前非监督分类结果不理想的一种修

8、正，可以将同种用地类型但属于不同类别的类合并在一起，具体操作如下：点击ERDAS图标面板工具条下的interpreter图标，选择GISAnalysis/下的Recode命令，打开Recode对话框：1990年的结果如图2-3：图2-32006年的结果如图2-4：图2-4（一）分类精度评价1.ERDAS图标面板工具条：点击Classifier图标→Classification菜单→选择AccuracyAssessment菜单项→打开AccuracyAssessment