资源与的环境信息系统实验指导书.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 几何校正（目的是使每一个采样点能最精准的定位到正确的像元） 提供了一幅经