地表反射率计算.docx

地表反射率计算.docx

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

地表反射率计算

算计射率石市地表反黄

一、数据预处理

1、打开：

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

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

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

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

波段→打开。

波段进行合成。

4、5年遥感影像数据的3、2、合成：

对黄石市2000

感遥2000年FileBasicTools→LayerStacking→Import→选择黄石市

→2000_band543_hecheng→波段543影像数据的、、→Ok→Choose命名（）

打开→Ok

黄石市遥感影像。

、裁剪：

用黄石市边界矢量数据裁剪合成后的20003

波段；53、4、遥感影像数据的

（1）打开合成后的黄石市2000年

）→打开→Ok2000_band543_hechengFile→OpenImageFile→选图（）打开黄石市边界矢量数据；2（

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

备注：

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

打开。

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

BasicTools→Masking→BulidMask→Display#1→Options→Import

EVFS→选图（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（

→2000_band543_hecheng_clip→SpatialBasicTools→ResizeData名命→Ok→Ok→Choose→→OkSubset→ROI/EVF→选图（EVF:

111）

开→Ok打_xiugai）→（2000_band543_hecheng_clip

算计地表反射率2、

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

件6S）运行（1软

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

）输入数据3（

）（几何条件conditions①Geometrical

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

（气溶胶浓度）Aerosol⑤model（concentration））输入光学厚度thevisibility:

16（

目标的高度）oftarget（⑥Thealtitude

0.02km）高度＝（目标xps:

-0.02

传感器的高度）（altitude⑦Thesensor

）传感器在卫星高度xpp:

-1000（

状况）conditions波段（⑧Thespectral

即即29TM5,30，即，即，即T25iwave[25-30]:

27（即M1,26TM227TM328TM4

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

）TM6是。

但是tm6热

地面反射率）reflectance（type）（⑨Ground

）地面均匀inhomo:

0（

）idirec:

0（地面无方向影响

）（反射率不随波长变化igroun:

0

0）地面波段反射率＝ro:

0（

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

）（表观反射率rapp:

-0.2

）书写方程4（

。

重命名为成的文件形完成后将?

段状况）时（波T波段4的操作同上，只是在输入hespectralconditions

。

到即可，得改为28

radiance）-公式y=xa*（measured个文件，根据这两打开?

，得到的是xb;acr=y/（1.+xc*y）输入的是影像所测的辐射亮度,Radiance

反射率影像。

为一个公式：

将这两个公式合并

radiance）-xb;y=xa*（measured

acr=y/（1+xc*y）

合并后：

radiance）-xb]}acr=[xa*（measuredradiance）-xb]/{1+xc*[xa*（measured、地表反射率计算2

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

（若添加成功说明公式正确）→Ok→选toList→Add

地表tm3_波段→Choose→命名为的2000_band543_hecheng_clip_xiugai第三

→Ok反射率→打开

地表反射率tm4、相同步骤操作最后得到3的数据4波段

公式：

TM4形成

（0.00470*b1-0.01979）/[1+0.04569*（0.00470*b1-0.01979）]地表反射率tm4