地表反射率计算教案分析.doc

黄石市地表反射率计算 一、数据预处理 1、打开：用ENVI5.1将黄石市2000年遥感影像数据的3,、4、5波段打开 （1）用鼠标左键双击ENVI5.1图标，打开ENVI5.1程序； （2）打开黄石市2000年遥感影像数据的3,、4、5波段。 File→Open Image File→选择黄石市2000年遥感影像数据的3、4、5波段→打开。 合成：对黄石市2000年遥感影像数据的3、4、5波段进行合成。 Basic Tools→Layer Stacking→Import File→选择黄石市2000年遥感影像数据的3、4、5波段→Ok→Choose→命名（2000_band543_hecheng）→ 打开→Ok 3、裁剪：用黄石市边界矢量数据裁剪合成后的2000黄石市遥感影像。 （1）打开合成后的黄石市2000年遥感影像数据的3、4、5波段； File→Open Image File→选图（2000_band543_hecheng）→打开→Ok 打开黄石市边界矢量数据； Vector→Open Vector File→选图（黄石市边界范围.evf)→打开 备注：建立掩膜时一定要将2000_band543_hecheng和黄石市矢量边界的影像打开。 （3）以黄石市边界矢量数据建立掩膜； Basic Tools→Masking→Bulid Mask→Display #1→Options →Import EVFS→选图（111）→Ok→Choose→命名（2000_band543_hecheng_yanmo)→打开→Apply （4）应用掩膜； Basic Tools→Masking→Apply Mask→2000_band543_hecheng→ Select Mask Bang→2000_band543_hecheng_yanmo→Ok→Ok→Choose→命名（2000_band543_hecheng_clip)→打开→Ok （5）备注：若裁剪后Scroll窗口内黑色背景面积太大可以进行调整。 Basic Tools→Resize Data→2000_band543_hecheng_clip→Spatial Subset→ROI/EVF→选图（EVF:111）→Ok→Ok→Ok→Choose→命名（2000_band543_hecheng_clip _xiugai)→打开→Ok 地表反射率计算 利用6S模型进行大气辐射校正，从而得到计算地表反射率的公式和三个参数 运行6S软件 （2）打开黄石市遥感图像头文件和6S操作手册查看要输入的数据 （3）输入数据 ①Geometrical conditions (几何条件) igeom [0-7]:7(在6S操作手册p1第34页对应上TM影像，即Landset对应的数字为7) ②时间、经度、纬度（在黄石市头文件）：时间均换算成以小时为单位 (Atmospheric model（大气模型） idatm [0-8]: 2（中纬度夏季大气模式）（在6S操作手册p1第35页） ④Aerosol model(type)（气溶胶类型） （在6S操作手册p1第36页） iaer[0-12]:3 ⑤ Aerosol model (concentration)（气溶胶浓度） the visibility:16 (输入光学厚度) ⑥The altitude of target （ 目标的高度） xps: -0.02 (目标高度＝0.02km) ⑦The sensor altitude （ 传感器的高度） xpp: -1000 (传感器在卫星高度) ⑧The spectral conditions （ 波段状况） iwave[25-30]:27 （25即TM1,26即TM2，27即TM3，28即TM4，29即TM5,30即TM6。但是tm6是热红外波段，不适用于我们的此时地表反射率的计算。） ⑨Ground reflectance (type)（地面反射率） inhomo: 0 (地面均匀) idirec: 0 (地面无方向影响) igroun: 0 (反射率不随波长变化) ro: 0 (地面波段反射率＝0) ⑩Atmospheric correction mode （激活大气订正的方式） rapp: -0.2 （表观反射率) （4）书写方程 (完成后将形成的文件重命名为。 波段4的操作同上，只是在输入The spectral conditions （ 波段状况）时改为28即可，得到。 (打开文件，根据这两个公式y=xa*(measured radia