地表反射率计算.docx
《地表反射率计算.docx》由会员分享,可在线阅读,更多相关《地表反射率计算.docx(15页珍藏版)》请在冰豆网上搜索。
![地表反射率计算.docx](https://file1.bdocx.com/fileroot1/2023-2/22/f2006ad6-3bd4-4e93-af6d-dbc8c1d5cdd2/f2006ad6-3bd4-4e93-af6d-dbc8c1d5cdd21.gif)
地表反射率计算
黄石市地表反射率计算
一、数据预处理
1、打开:
用ENVI5.1将黄石市2000年遥感影像数据的3,、4、5波段打开
(1)用鼠标左键双击ENVI5.1图标,打开ENVI5.1程序;
(2)打开黄石市2000年遥感影像数据的3,、4、5波段。
ImageFile→选择黄石市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波段;
ImageFile→选图(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地表反射率