地表反射率计算
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
算计射率石市地表反黄
一、数据预处理
1、打开：用 ENVI5.1 将黄石市 2000 年遥感影像数据的 3,、4、5 波段打开
（1）用鼠标左键双击 ENVI5.1 图标，打开 ENVI5.1 程序；
（2）打开黄石市 2000 年遥感影像数据的 3,、4、5 波段。
File→Open Image File→选择黄石市 2000 年遥感影像数据的 3、4、5
波段→打开。
波段进行合成。 4、5年遥感影像数据的 3、 2、合成：对黄石市 2000
感遥2000 年 File Basic Tools→LayerStacking→Import →选择黄石市
→2000_band543_hecheng→波段543影像数据的、、→Ok→Choose 命名（）
打开→Ok
黄石市遥感影像。、裁剪：用黄石市边界矢量数据裁剪合成后的20003
波段；5 3、4、遥感影像数据的（1）打开合成后的黄石市2000 年
）→打开→Ok2000_band543_hecheng File→OpenImage File→选图（）打开黄石市边界矢量数据；2（
→选图（黄石市边界范围.evf)→打开File Vector→OpenVector
备注：建立掩膜时一定要将 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 （
→2000_band543_hecheng_clip→SpatialBasic Tools→Resize Data名命→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)
度（在黄石市头文件）：时间均换算成以小时为单位②时间、经度、纬
大气模型）（?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
(
目标的高度） of target （⑥The altitude
0.02km)高度＝ (目标xps: -0.02
传感器的高度）（altitude ⑦The sensor
)传感器在卫星高度xpp: -1000
(
状况） conditions 波段（⑧The spectral
即即29TM5,30，即，即，即T25iwave[25-30]:27 （即
M1,26TM227TM328TM4
红外波段，不适用于我们的此时地表反射率的计算。）TM6是。但是tm6热
地面反射率）reflectance (type)（⑨Ground
)地面均匀inhomo: 0 (
)idirec: 0 (地面无方向影响
)(反射率不随波长变化igroun: 0
0)地面波段反射率＝ro: 0
(
（激活大气订正的方式）⑩Atmospheric correction mode
)（表观反射率rapp: -0.2
）书写方程4（
。重命名为成的文件形完成后将?
段状况）时（波 T波段4的操作同上，只是在输入he spectral conditions
。到即可，得改为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*(measured
radiance)-xb]/{1+xc*[xa*(measured、地表反射率计算2