遥感导论实习报告_遥感导论实验报告
遥感导论实习报告由刀豆文库小编整理,希望给你工作、学习、生活带来方便,猜你可能喜欢“遥感导论实验报告”。
《遥感技术与应用》
实习报告
姓 名:侯会茹 学 号:04207113 系 别:07信息工程系 班 级:一班
专 业:地理信息系统
2010年 6 月17日 3S区域土地利用变化动态监测
一.实习目的及意义:
1.1了解土地利用动态监测的过程;
1.2学会获取遥感图像的方法,遥感图像的纠正,裁剪,合成,监督分类和非监督分类;
1.3掌握土地利用现状(1990年)及更新(2000年),变更后的影像内容处理;
1.4学会使用ERDAS进行图像的处理;
1.5能够运用ERDAS、ARCGIS、ARCVIEW进行综合成图,完成调查实习期间的用地类型的数量,质量,面积及分布状况。
二. 实习任务:
2.1对太原市土地利用类型进行分析,统计面积;
2.2生成太原市土地利用现状图,并完成土地利用变化分析;
三. 实习时间:
2010年5月—6月
四.实习地点:
中国地质大学长城学院图书馆 1111实验室
五.实习内容
5.1 数据的收集
打开FlashFXP文件,进行影像下载。5.2 数据的处理(1990年太原影像)5.2.1影像格式转换
点击Import,在对话框中进行设置,转化为.img格式
5.2.2图像融合(将六个波段的图像融合到一起)
在Interpreter中选择Utilities选项下的Layer Stack,进行六个波段的合成5.2.3提取合适的波段进行合成(1)对图像进行相关分析
点击Modeler下的Model Maker,在出现的对话框中进行如图所示添加
在上图中的图形上双击,从上到下依次设置如下格式,图一
图二
图三
设置好后运行,查看相关分析结果。在Tools工具下的Edit Text File,添加相关分析结果。
(3)选择出相关系数较小的三个波段,根据上图数据,选出3,4,5三个波段(4)对选择出的波段进行合成在Interpreter中选择Utilities选项下的Layer Stack
1990年 三波段合成5.2.4用行政区划图裁剪图像
(1)核对边界图的投影类型(必须与图像一致):然后进行核对(2)图像的叠加:在File中同时打开TM图像和矢量图
(3)裁剪:点击Interpreter中的Utilities菜单下的Vector to Raster,在对话框中进行添加,在Interper中选择Utilities下的Mask,对图像进行裁剪
1990年 裁剪后图像
5.2.5对图像进行投影变换(把UTM投影变换成等积投影,以保证准确的面积计算)
(1)打开图像,在视窗中点击,出现ImageInfo的对话框,点击edit下面的Add/Change Projection,出现Projection Chooser对话框,按下图设置:
设置完毕点击save保存,出现对话框,设置如下
(2)点击DataPrep下的Reproject Images,出现对话框,如下图设置:
图一
图二 1990年 投影后图像
5.2.6监督分类
(1)训练区样本的选择:点击Clifier下的Signature Editor,出现Signature Editor对话框。在View1中点击AOI下的Tools按钮,出现AOI框,点击在窗口的图像上画出一个样本,在Signature Editor对话框中点击按钮,将样本添加到对话框中。每个类型的利用类型画5个样本,如下图:(之后保存)
(2)监督分类:点击Claifier下的Supervised Claification,出现对话框并进设置,生成1990年 监督分类后图像
:5.2.7分类后评价
(1)点击Claifier下的Accuracy Aement,open打开监督分类后图像,出现对话框;原始图像和分类后图像的链接:点击在Accuracy Aement对话框中的在一起。
(2)设置随机点的色彩:在在Accuracy Aement对话框中,点击View下
按钮,然后在原始图像上点击一下,就把两个图像连接的Change Colors,出现对话框并且设置
(3)产生随机点:在Accuracy Aement对话框中,点击Edit下的Create/add Random Points,出现对话框,设置
(4)显示随机点:在Accuracy Aement对话框中,点击View下的Show All,原始图像上出现随机点,如图:
(5)显示随机点的类别:在Accuracy Aement对话框中,点击Edit下的Show Cla Values,上面的对话框改变为下图:
(增加了CLASS的值)
(6)输入随机点的实际类型:在Accuracy Aement对话框中的Reference一栏写入随机点的真实类型,如下:
(7)输出评价报告:在Accuracy Aement对话框中,在Report下的Options中的所有项打钩,之后选择Accuracy Report导出分类报告:
(8)分析评价报告,确认监督分类是否在合理的误差范围之内
(1990年为0.800、2000年的为0.850)
5.2.8分类后碎斑处理
(1)聚类统计:在Interpreter的GIS Analysis中,选择Clump选项,进行设置:
(2)在Interpreter的GIS Analysis中,选择Eliminate选项,进行设置:
(3)导出影像:运用Vector下的Raster to Vactor工具: 5.3 实习结果图的生成与分析
5.3.1在ARCGIS中打开碎斑处理后的图像,进行地图处理得到土地利用结 果图
5.3.2对两年的土地利用图进行叠加,生成土地利用变化图
(1)修改2000年土地利用图属性表中土地利用类型的列名称:在Arcmap中打开2000年的土地利用图属性表,新建属性列,其名为CODE,之后将原土地利用类型添加到此列中。
(2)在Arcmap中进行图像的叠加:选择Analysis中的Overlay工具下的Intersect,之后按图进行设置
(3)在ArcView中添加叠加后的图像
(4)生成土地利用变化属性表: 在ArcView中点击,在出现的属性表中添加新的一列,其名为“bianhuan”:在Table中点击Star Editing,之后在Edit中选择Add Field,设置。选中新添加的一列,点击之后点击new set。点击,给新的属性值起名为0808,点击OK,在出现的对话框中分别是土地利用的数值相等,对所有可能出现的组合进行如上操作后停止编辑。(5)生成土地利用变化图
显示图像,在图像标签处双击,出现的对话框中进行设置:
在VIEW中点击Layer生成土地利用变化图并进行必要的地图设置,之后在File中的Export下出图:
5.4转移矩阵的生成在ArcView中打开叠加后的图像,点击
打开属性表,之后点击,分别令Grid-Code为8、9、10、12、13(以Grid-Code=8为例)
点击New Set,选出Grid-Code=8的行,点击求和符号,进行如下设置:再点击,进行计算,先保存路径,field选择area,summarize by选择sum,然后点击add,即ok,然后把计算的结果用excel打开,进行转置的操作
用Excel打开生成的表,设计成最终的表格实习感想
在这段实习时间里,学到了很多在课本学习不到的东西,将所学知识运用到实践中去的能力,感谢我的指导老师,刘琦。此次实习学会了如何应用Erdas软件将平面地图转化成软件所需要的文件格式类型,基本掌握了对一个区域进行土地利用变化监测的步骤,熟悉了相关软件的使用,复习了ARCMAP软件的操作,由此掌握了专业软件的结合使用——ARCMAP,ARCVIW,ERDAS.了解各个软件功能的结合使用。