地表反射率计算

黄石市地表反射率计算一、数据预处理1、打开：用ENVI5.1将黄石市2000年遥感影像数据的3,、4、5波段打开（1）用鼠标左键双击ENVI5.1图标，打开ENVI5.1程序；（2）打开黄石市2000年遥感影像数据的3,、4、5波段。 FileOpen Image File选择黄石市2000年遥感影像数据的3、4、5波段打开。2、 合成：对黄石市2000年遥感影像数据的3、4、5波段进行合成。 Basic ToolsLayer StackingImport File选择黄石市2000年遥感影像数据的3、4、5波段OkChoose命名（2000_band543_hecheng）打开Ok3、裁剪：用黄石市边界矢量数据裁剪合成后的2000黄石市遥感影像。（1）打开合成后的黄石市2000年遥感影像数据的3、4、5波段； FileOpen Image File选图（2000_band543_hecheng）打开Ok（2） 打开黄石市边界矢量数据；VectorOpen Vector File选图（黄石市边界范围.evf)打开 备注：建立掩膜时一定要将2000_band543_hecheng和黄石市矢量边界的影像打开。（3）以黄石市边界矢量数据建立掩膜；Basic ToolsMaskingBulid MaskDisplay #1Options Import EVFS选图（111）OkChoose命名（2000_band543_hecheng_yanmo)打开Apply（4）应用掩膜；Basic ToolsMaskingApply Mask2000_band543_hechengSelect Mask Bang2000_band543_hecheng_yanmoOkOkChoose命名（2000_band543_hecheng_clip)打开Ok （5）备注：若裁剪后Scroll窗口内黑色背景面积太大可以进行调整。Basic ToolsResize Data2000_band543_hecheng_clipSpatial SubsetROI/EVF选图（EVF:111）OkOkOkChoose命名（2000_band543_hecheng_clip _xiugai)打开Ok2、 地表反射率计算1、 利用6S模型进行大气辐射校正，从而得到计算地表反射率的公式和三个参数（1） 运行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页）iaer0-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 （ 波段状况）iwave25-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 radiance)-xb; acr=y/(1.+xc*y) ,输入的是影像所测的辐射亮度Radiance，得到的是反射率影像。将这两个公式合并为一个公式：y=xa*(measured radiance)-xb;acr=y/(1+xc*y)合并后：acr=xa*(measured radiance)-xb/1+xc*xa*(measured radiance)-xb2、 地表反射率计算打开ENVI软件Basic ToolBand Math输入(0.00336*b1-0.04143)/1+0.07179*(0.00336*b1-0.04143)Add to List（若添加成功说明公式正确）Ok选2000_band543_hecheng_clip_xiugai的第三波段Cho