地表反射率计算.docx

地表反射率计算.docx

《地表反射率计算.docx》由会员分享，可在线阅读，更多相关《地表反射率计算.docx（15页珍藏版）》请在冰点文库上搜索。

地表反射率计算

黄石市地表反射率计算

一、数据预处理

1、打开：

用ENVI5.1将黄石市2000年遥感影像数据的3,、4、5波段打开

（1）用鼠标左键双击ENVI5.1图标，打开ENVI5.1程序；

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

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

2、合成：

对黄石市2000年遥感影像数据的3、4、5波段进行合成。

BasicTools→LayerStacking→ImportFile→选择黄石市2000年遥感影像数据的3、4、5波段→Ok→Choose→命名（2000_band543_hecheng）→

打开→Ok

3、裁剪：

用黄石市边界矢量数据裁剪合成后的2000黄石市遥感影像。

（1）打开合成后的黄石市2000年遥感影像数据的3、4、5波段；

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

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

Vector→OpenVectorFile→选图（黄石市边界范围.evf）→打开

备注：

建立掩膜时一定要将2000_band543_hecheng和黄石市矢量边界的影像打开。

（3）以黄石市边界矢量数据建立掩膜；

BasicTools→Masking→BulidMask→Display#1→Options→ImportEVFS→选图（111）→Ok→Choose→命名（2000_band543_hecheng_yanmo）→打开→Apply

（4）应用掩膜；

BasicTools→Masking→ApplyMask→2000_band543_hecheng→

SelectMaskBang→2000_band543_hecheng_yanmo→Ok→Ok→Choose→命名（2000_band543_hecheng_clip）→打开→Ok

（5）备注：

若裁剪后Scroll窗口内黑色背景面积太大可以进行调整。

BasicTools→ResizeData→2000_band543_hecheng_clip→SpatialSubset→ROI/EVF→选图（EVF:

111）→Ok→Ok→Ok→Choose→命名（2000_band543_hecheng_clip_xiugai）→打开→Ok

2、地表反射率计算

1、利用6S模型进行大气辐射校正，从而得到计算地表反射率的公式和三个参数

（1）运行6S软件

（2）打开黄石市遥感图像头文件和6S操作手册查看要输入的数据

（3）输入数据

Geometricalconditions（几何条件）

igeom[0-7]:

7（在6S操作手册p1第34页对应上TM影像，即Landset对应的数字为7）

时间、经度、纬度（在黄石市头文件）：

时间均换算成以小时为单位

③Atmosphericmodel（大气模型）

idatm[0-8]:

2（中纬度夏季大气模式）（在6S操作手册p1第35页）

Aerosolmodel（type）（气溶胶类型）（在6S操作手册p1第36页）

iaer[0-12]:

3

Aerosolmodel（concentration）（气溶胶浓度）

thevisibility:

16（输入光学厚度）

Thealtitudeoftarget（目标的高度）

xps:

-0.02（目标高度＝0.02km）

Thesensoraltitude（传感器的高度）

xpp:

-1000（传感器在卫星高度）

Thespectralconditions（波段状况）

iwave[25-30]:

27（25即TM1,26即TM2，27即TM3，28即TM4，29即TM5,30即TM6。

但是tm6是热红外波段，不适用于我们的此时地表反射率的计算。

）

Groundreflectance（type）（地面反射率）

inhomo:

0（地面均匀）

idirec:

0（地面无方向影响）

igroun:

0（反射率不随波长变化）

ro:

0（地面波段反射率＝0）

Atmosphericcorrectionmode（激活大气订正的方式）

rapp:

-0.2（表观反射率）

（4）书写方程

①完成后将形成的文件

重命名为

。

波段4的操作同上，只是在输入Thespectralconditions（波段状况）时改为28即可，得到

。

②打开

文件，根据这两个公式y=xa*（measuredradiance）-xb;acr=y/（1.+xc*y）,输入的是影像所测的辐射亮度Radiance，得到的是反射率影像。

将这两个公式合并为一个公式：

y=xa*（measuredradiance）-xb;

acr=y/（1+xc*y）

合并后：

acr=[xa*（measuredradiance）-xb]/{1+xc*[xa*（measuredradiance）-xb]}

2、地表反射率计算

打开ENVI软件→BasicTool→BandMath→输入（0.00336*b1-0.04143）/[1+0.07179*（0.00336*b1-0.04143）]

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

3、相同步骤操作最后得到tm4地表反射率

波段4的数据

形成TM4公式：

（0.00470*b1-0.01979）/[1+0.04569*（0.00470*b1-0.01979）]

tm4地表反射率