《遥感应用实习报告-康永正》由会员分享，可在线阅读，更多相关《遥感应用实习报告-康永正（21页珍藏版）》请在金锄头文库上搜索。
1、遥感应用实习报告姓名：康永正 班级：GIS 0701 学号：2007303204035实习目的：1.掌握ERDAS软件的基本使用技巧。2.会使用ERDAS来处理遥感图像以及根据特定的目的来分析遥感图像。3.掌握处理遥感影像的一般流程和步骤方法。4.会利用ERDAS进行一些简单的建模。实习安排：实习共分为两周；两部分：第一部分：遥感影像的处理与应用第二部分：遥感温度反演和高光谱遥感模型反演第一周(第一部分)：1.遥感影像的获取。2.遥感影像的校正。3.遥感影像的识别。4.土地利用变化分析。5.植被变化分析。6.立体土地利用图的制作。第二周（第二部分）：1.为期两天的野外实习。2.遥感温度反演和高
2、光谱遥感温度反演。实习一：遥感影像的处理与应用实验目的：利用ERDAS通过对原始遥感影像的获取，校正，识别得到进行最后的相关数据变量分析的成果图以及相关的数据。实习的过程中需要熟悉和了解目前主流的遥感图象处理软件ERDAS的主要功能模块，在此基础上，掌握视窗操作模块的功能和操作技能。在遥感图像的几何校正等后续实习中，通过实习操作，掌握遥感图像几何校正的基本方法和步骤，深刻理解遥感图像几何校正的意义。利用所得到的经过几何校正后的遥感影像进行监督分类和NDVI植被指数分析，并且能对所得到的结果进行相关数据处理与分析。实验步骤：(一) 遥感影像的获取。在ERDAS中使用预处理：Interpreter
3、UtilitiesLayer Stack将TIF格式的1992年的7个光谱波段，2001年的除了第六个的其他六个波段的影像数据组合多波段数据，并且需要忽略0值。组合后的图像*.tif格式转为*.img格式，得到1992和2001年的假彩色影像图。(二) 遥感影像的校正。1.在两个view窗口里分别打开之前获取的img格式的TM假彩色影像图，另个窗口里同时加入4块带有经纬度坐标的河流的shapefile文件。2.用ERDAS中的DATAPREP里的image geomatic correction功能模块。依次点选遥感影像和shape图，并选择polynomial在两幅图上选择对应的点进行对准校
4、正，获得校准模型以及每个点的选点误差。3.将校正输出后的遥感影像打开，叠加上河流的shape图，用utility里的功能进行校正质量检查。通过几何校正给遥感影像上的每个像元附上了相应的坐标。(三) 遥感影像的识别。在经过遥感影像的波段合成与几何校正等前期处理工作后，在此基础上进行遥感影像的识别监督分类。建立训练区矢量文件，定义分类模板为以下六类：(1)林地暗红色 、黑色 (2)旱地鲜红(3)水田黄绿、暗红色 (4)城镇用地蓝白色 (5)水域黑色 (6)荒地（河漫滩,仅2001年有，云）这里运用的是ERDAS里的classifier的signature editor模块；1.打开1992和200
5、1两副校正图，对应各类地物进行训练区选点，其中同类地物根据不同的颜色应该先区分归类。2.选择evaluate中的contingency查看六类地物的分类精度情况，在精度低的地方多取几个训练区。精度符合要求后，用classify-supervised将分类后的图像输出。3.利用Erdas中的interpretator-gisanalysis-recode，载入分类后的图像，将相同的地物归为一类，然后输出分类合并后的遥感影像，在raster-attributes里进行颜色的修改。得到两年的监督分类图后需要对监督分类结果进行分类精度评估（将分类图像中的特定象元与已知分类的参考象元进行比较）用clas
6、sifier-accuracy assessment打开分类后的图editinport user-defined points分别导入 check92.txt和check01.txt设置 comma，其他使用默认值，输出分类评价报告：reportaccuracy report报告的结果如下1992年分类评价报告： 2001年分类评价报告：(四) 土地利用变化分析由于两年的图像大小差异需要首先对监督分类后的1992和2001年的两幅图像使用interpretator-utilities-subsetimage进行遥感影像中相同的区域Subset(x:568121684371,y:34273593
7、288369)的裁剪。然后选择GIS Analysis-matrix用矩阵方法对年份不同的土地利用图做变化分析，得出一幅可以体现变化的成果图. 最终得到土地利用变化表如下: 2001年1992年林地旱地水田城镇水域荒地合计林地55331.77308.457701511259.66713.552260.26159888.6旱地11410310055120393738562.123091.714112.9494357.7水田20103811547427751658356.3666676340.14725391.4城镇29734.89979.1153435.921623.57451.464062.6
8、126287.4水域16078.74898.711772.39633.9648192191792492.66荒地383.042719.5310573.51631.161498.681842.318648.21合计416669.2240930.8634249.7141066.6153614.430535.21617066分析：根据土地利用变化的数据可以得出在2001年时，6类地物分别由1992年时的6类不同的地物所转变的面积大小的柱状图（上图）；再结合10年间6类地物向其他类型转变的变化比例图，可以看出1.水田和旱地的面积是比较大的，分别占了总面积的51.9%和19.8%。2.10年里林地，旱地
9、，城镇，荒地大部分都变成了水田，分别占各自总共的62.5%，61.07%和59.6%。3.林地，荒地还有一部分转变为城镇和旱地，自身一直保留下来的占少数。4.各种地物都有一部分转变成了城市，说明了城市在10年里的扩张。5.水域在10年里，主要是原来的水域形成的，保留原来的部分占72.24%，部分形成了城镇和水田。6.再根据表中的总计的值，计算10年内各地物的动态增减差值，可以看出林地，旱地是负增长，分别减少72.3%和55.26%；7.水田，城镇总面积是显著增加的，分别为33.58%和113%，水域略微有波动，幅度不大。从1992-2001十年中的土地利用变化可以看出，随着城镇的面积的不断扩张
10、，用于农业的水田面积整体上也有了显著的提高。(五) 植被变化分析首先将经过几何校正后的两年的遥感影像，同样取相同区域经行操作，用interpretator-utilities-subsetimage进行遥感影像中相同的区域Subset(x:568121684371,y:34273593288369)的裁剪。通过InterpreterSpectral enhancementindices得到经过归一化植被指数运算处理的影像，打开Erdas的核心功能modeler，分别载入影像，应用模型ndvi-cover.gmd计算植被覆盖度。在功能模块interpretator-utilities-opera
11、tors里进行两幅植被覆盖度影像的相减，得出植被覆盖度变化趋势图。92和01年植被覆盖度图如下：通过所得到的由92减去01年的植被覆盖度影像所得的植被覆盖度变化趋势图根据NDVI=（NIR-R）/（NIR+R）NIR近红外波段R 可见光波段-1=NDVI=1，负值表示地面覆盖为云、水、雪等，对可见光高反射；0表示有岩石或裸土等，NIR和R近似相等；正值，表示有植被覆盖，且随覆盖度增大而增大； 可以知道遥感影像的巨大多数区域有植被覆盖，部分是云，水与岩石；NDVI=0.2左右的对应的植被覆盖的面积最大，在0.5以后的面积就很小了。此地的植被覆盖度从92年到01年间，植被覆盖度在一定程度上降低了。
12、强植被覆盖度的地区面积不大，主要体现以水田为主的较稀疏的植被覆盖程度。(六) 立体土地利用图的制作选择dataprep-createsurface-surfacingtools将等高线shp图生成DEM的img图，其中要选择elev字段。在virtual gisvirtualgis viewer打开dem和土地利用图，叠加DEM与土地利用图，得到三维模型。实习二：遥感温度反演和高光谱遥感模型反演实验目的：通过使用在第一周里获得的1992年和2001年的经过几何校正以及NDVI指数处理后的遥感影像，使用Erdas软件，学会自主编写相应的模型，直接通过模型导出最后的反演结果图。懂得建立简单的模型，
13、结合所提供的数据资料，进行泥沙的高光谱反演。实验步骤：一、遥感温度反演遥感温度的反演过程主要包括数据采集TM/ETM影像， 归一化植被指数，计算比辐射率(Emissivity),计算温度反演结果与分析。由于前两步已经完成并有成果图，因此在建立温度反演模型时当做输入数据使用。我们第二组负责采集灌木样品，首先由老师给我们讲解了一下地物光谱观测的基本原理和基本步骤，以及我们学校唯一的一台地物光谱观测仪和它的相关的软件的使用方法。反射光谱分析方法：在物质成份分析以及物质鉴别中扮演者非常重要的角色。很多物质都有其特定的吸收峰或吸收带，我们只要测出该物质的反射，既可以知道该物质的吸收峰和吸收带，从而可以对
14、该物质进行一系列分析，如含何种成分、含量多少等。相对于细说光谱分析法，反射光谱对物质的定量分析要困难一些，它受到很多因素影响：表面的粗糙程度、表面大小和形状、背景光以及夹持装置等。根据温度反演的理论基础可知：1.地表比辐射率对地表温度反演精度影响很大，研究发现地表比辐射率与植被指数高度相关，建立关系模型：其中NDVI的取值范围 0.157,0.727。由于水体和岩石裸土的NDVI值小于等于0，因此当NDVI不在正值区间时，都取0。即2.进行温度反演时是运用Planck方程计算亮度温度。参考模型取3.反解Planck函数获取地表真实温度模型编写遥感温度反演模型：打开Erdas，使用里面的modeler模块，建立模型的框架，流程图如下所示：在n1_1992jh2输入的是1992年经过NDVI指数处理的影像图。1.地表比辐射率处理。 EITHER 0 IF ( $n1_1992jh2 0.727 ) OR (1.0094+(0.047 *LOG ( $n1_1992jh2 ) )OTHERWISE 。将结果保留进行下一步操作，同时在n9_1992jh中输入1992年经过几何校正的裁剪图。然后在下一 步的温度反演处理中使用其第六波段作为的值代入公式计算。得到 2. EITHER0IF($n5_memory=0)OR($n9_1992jh(6)-1.74)/(0.77
