武大遥感应用模型课程综合实习.docVIP

武汉大学教学实验报告 遥感信息工程学院 遥感科学与技术专业 2013年 12 月 8 日 实验名称 遥感应用模型实习 指导教师 XXX萌 姓名 XXX 班级 10011 学号 XXXX 成绩 预习部分 实验目的 实验基本原理 主要仪器设备（含必要的元器件、工具） 1 实验目的 《遥感应用模型实习》是一门综合设计型的实验教学课程，要求在学完《遥感应用模型》理论课的基础上，能够运用一些遥感模型的理论知识去解决实际具体的遥感应用问题。 通过运用遥感模型解决实际遥感应用问题的过程，锻炼发现问题、分析问题以及综合运用所掌握的遥感理论知识和技术解决实际问题的能力。加强对遥感应用模型》课程知识的理解和深化,训练实际编程和操作能力，培养具体解决某一遥感应用实际问题的能力。 大气溶胶光学厚度定量反演 1.2 气象因子与NDVI的相关性分析大气溶胶光学厚度定量反演’ Optique Atmospheri que）5S模式基础上发展而来的。6S模型的适用范围是太阳的短波0.25um-4.0um，适合计算无云大气的辐射传输，在其程序中大气模型可以自己定义，也可以按程序中的大气模型定义，其中程序自带的大气模型与 LOWTRAN中6种标准大气模型一致。 在大气气溶胶的反演中，就是利用6S模型来建立MODIS数据的第1波段和第3波段的查找表，其中包括了太阳天顶角、卫星天顶角、相对方位角、大气气溶胶光学厚度、地面反射率和表观反射率。太阳天顶角、卫星天顶角和相对方位角分别可以从角度数据的3、4、7波段的余弦图像中利用反余弦函数获得。表观反射率可以通过1、3波段图像，通过定标，乘以对应的定标系数获得。 由于第7波段是近红外波段，其通过定标后获得反射率即可以视作该波段的地表反射率。根据如下公式： （1） （2） 其中，表示第1波段的地面反射率，表示第3波段的地面反射率，表示第7波段的地面反射率。 根据式（1）、（2），可以很快的求出第1、3波段的地面反射率，因此结合第1、3波段的查找表，利用对应的五个参数来查找剩下的一个参数——光学厚度。 2.2 气象因子与NDVI的相关性分析大气溶胶光学厚度定量反演气象因子与NDVI的相关性分析 实验数据、表格及数据处理 实验操作过程（可用图表示） 结论 1 实验数据、表格及数据处理 1.1 大气溶胶光学厚度定量反演 大气溶胶光学厚度定量反演——把余弦值化为度数值。 1.2 气象因子与NDVI的相关性分析 气象因子与NDVI的相关性分析大气溶胶光学厚度定量反演 (一)数据预处理 格式转换：用ERDAS的Import模块中，选择Generic Binary，把格式为BIN的250m、500m和1000m的MODIS数据转换为IMG格式，利用在弹出的对话框中设置参数，把数据格式设为BSQ，数据设为Unsigned int 16 bit；对于角度数据，数据格式也设置为BSQ，把数据类型设为IEEE 32 Bit Float。在下面的Rows和 Cols中根据附给的Readme文件中的数据来设置。在Bands中输入波段数，250米的有两个波段，500米的有5个波段，1000米的有29个波段。然后点OK完成格式转换。 像元合并：用Interpreter模块下的 Utilities-Degrade来完成像元合并。像元合并是把250米和500米分辨率的数据转化为1000米分辨率的数据。在像元合并时，在X Scale 和Y Scale中设置把多少个像元合并为一个像元。对于250米，都设为4，对于500米都设为2，点击OK即可以完成像元的合并。 裁剪：是用DataPrep模块下Subset Image来完成裁剪，裁剪重要的是知道要裁剪的区域范围。由于格式转换后的IMG文件没有头文件，不能用地理位置数据来裁剪，裁剪是根据老师给的已经裁剪好的数据反算得到裁剪的范围。 图层分离：由于1、2波段为250m分辨率，3—7波段为500m分辨率，需要从250米和500米分辨率的数据提出1波段、3波段、7波段数据。用Interpreter模块下的 Utilities-Layer Stack来完成图层分离。选择250m分辨率图像的第一个图层，即第1波段，选择500m分辨率图像的第一和第五个图层，即第3波段和第7波段。角度数据，分别分离出3、4、7三个波段，即相对方位角、卫星方位角和太阳方位角的余弦值数据。 角度数据的处理：用Model模块下的Model Maker来把余弦值转化为度数。新建一