地表反射率计算

黄石市地表反射率计算

一、数据预处理

1、打开：用ENVI5.1将黄石市2000年遥感影像数据的3,、4、5波段打开（1）用鼠标左键双击ENVI5.1图标，打开ENVI5.1程序；

（2）打开黄石市2000年遥感影像数据的3,、4、5波段。

Image File→选择黄石市2000年遥感影像数据的3、4、5波段→打开。

2、合成：对黄石市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波段；

Image File→选图（2000_band543_hecheng）→打开→Ok

（2）打开黄石市边界矢量数据；

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

二、地表反射率计算

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页）

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 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)-xb]} 2、地表反射率计算

打开ENVI软件→Basic Tool→Band Math→输入

(0.00336*b1-0.04143)/[1+0.07179*(0.00336*b1-0.04143)]

→Add to List（若添加成功说明公式正确）→Ok→选2000_band543_hecheng_clip_xiugai的第三波段→Choose→命名为tm3_地表反射率→打开→Ok