资源与环境信息系统实验指导书汇.docx

资源环境信息系统实验指导书实验一 太湖叶绿素a反演数据预处理1、数据简介（1）遥感影像数据：2009年10月6日环境小卫星影像数据，共4个波段，详细信息如表1 表1 环境小卫星数据为HJ-1B-CCD1参数介绍波段编号波长（微米）分辨率（米）Band 10.43-0.5230Band 20.52-0.6030Band 30.63-0.6930Band 40.76-0.9030（2）监测点数据： 38组监测点数据及矢量图层（时间为2009年10月6日，水质参数为叶绿素a浓度）（3）基准影像：精校正的TM数据2、数据预处理2.1 辐射定标与波段融合（目的是消除传感器本身产生的误差）（1）打开ENVI，在右侧工具箱【Toolbox】中找到【Extensions】模块，双击打开【ENVI　HJ1A1B Tools】工具（2）传感器类型选择CCD，即【Sensor Type】选择CCD，【Input Path】输入遥感影像，点击【Search】，【Out Path】设置输出路径，勾选【Calibration】和【Layer Stacking】，即辐射定标和波段融合。（3）点击Apply，当进度条达到100%时，即完成辐射定标和波段融合，关闭窗口即可。（4）加载经过辐射定标和波段融合的.img格式影像。主菜单【File】【open】，点击波段前的勾选框，选择3个波段进行显示（每次最多显示3个波段），点击【Load Data】。2.2 区域裁剪（裁剪太湖区域）（1）在右侧工具箱【Toolbox】中找到【Raster Management】模块，双击打开【Resize Data】工具（2）在打开窗口中，选中输入影像，点击【spatiat】，在打开窗口点击【image】，通过移动缩放红色方框，确定裁剪区域，点击ok，继续点击ok，直到“Resize Data Parameters”对话框。（3）设置输出路径，在“Resize Data Parameters”对话框，点击【Choose】按钮，设置输出路径及文件名。（4）点击ok，格式转化完成将自动加载。2.3 转换格式（需将数据转换为bil）格式（1）在右侧工具箱 【Toolbox】中找到【Raster Management】模块，双击打开【Convert Interleave】工具（2）双击上述裁剪影像，在打开窗口，【Output Interleave】选择“BIL”。（3）设置输出路径，点击【Choose】按钮，设置输出路径及文件名，输入文件名时，必须同时输入扩展名. bil。（4）点击OK，转换完成，数据自动加载。2.4 大气校正（目的是消除大气和光照等因素对地物反射的影响。）（1）在右侧工具箱 【Toolbox】中找到【Radiometric Correction】模块下的【Atmospheric Correction】模块，双击打开【FLAASH Atmospheric Correction】工具 （2）加载.dil格式影像，点击【Input Radiance Image】按钮，双击.dil影像，在弹出的窗口，选择“Use Single scale factor for all bands”，“Single scale factor”可设为100。 （3）设置输出路径，点击【Output Reflectance File】，设置输出路径及文件名，输入文件名时，必须同时输入扩展名.img。（4）设置文件夹路径，点击【Output Directory for FLAASH Files】，设置与上一步相同的输出路径。 （5）记录影像中心位置经纬度，在主界面点击【Crosshair】工具，通过移动红色交叉线，将交叉点移动到影像中心位置，在小窗口可以查看经纬度，并将其记录下来。（6）输入影像中心经纬度，将刚刚记录的经纬度填写在大气校正窗口。（7）输入影像采集时间。在原数据中，打开.XML格式的头文件，找到影像采集时间，并输入到大气校正窗口。（８）添加波段相应函数，点击【Multispectral Settings】按钮，选择【Filter Function File】按钮，点击【Open】，【New File】，添加相应的波段响应函数。（通过查看原数据头文件，影像来源于HJ1B CCD1，因此添加对应的.sli格式的波段响应函数）（8）传感器高度和地面高程分别设为660km,0.004km（分别与传感器高度及研究区域有关）（9）点击0k，输出大气校正结果。2.5 几何校正（目的是使每一个采样点能最精准的定位到正确的像元） 提供了一幅经过了精校正的TM数据进行配准，可以选择在ARCGIS或ENVI、ERDAS中配准。3、采集监测点各多段像元值 首先将监测点矢量数据、及校正后的影像数据加载到ARC