《地表反射率的计算》由会员分享,可在线阅读,更多相关《地表反射率的计算(12页珍藏版)》请在金锄头文库上搜索。
1、地表覆盖反射率的计算(6s 软件的应用)9 月 23 日首先在 envi 软件中打开已经处理好的真彩色影像(TM543 波段) ,我的影像因为没有居中,所以首先进行了裁剪,让影像满幅居中再操作。1、打开遥感影像,并裁减居中:先打开 7 个波段影像, ,选中 543 , 。合成, ,在弹出的对话框中点 ,按住 ctrl 再选中这 3 项,点 ok, 命名为 。打开矢量边界 , , ,选中 。建立掩膜 , , , , , 重命名为 ,点 ,形成掩膜文件。再应用掩膜 ,选 文件,点,选 , , , 重命名为 ,形成影像。所以接下来对背景进行裁剪, , ,选 , , ,选 , , , , ,重命名为
2、,形成影像。2、让#1 和 Scroll 中的红方框大致居中,在 #1 中任意位置双击弹出 “光标位置评估” ,或者右键找出也可。 (可是我不太清楚调出这个的目的?)3、寻找我们应用的黄石市遥感影像中头文件为 MTL.txt 的文件,以写字板的形式打开,方便查看遥感影像的具体信息。找到影像获取的时间即“DATE _ACQUIRED”,这个原始的影像获取时间才是我们需要的,不要被其他的信息误导。因为 6s 识别不了具体的时分秒,所以我们需要将具体时间换算成小时,即此处的 02:26:32 应转换成2.43 小时。4、打开 中的 ,我们在运行 6s 的一切操作,都是按着这个步骤来的,但是期间会出现
3、一些专业术语的特定要求,所以我们需要打开另外的文件 ,书名如下:打开到 35 页,IGEOM,从对应上我们找到 TM 影像,即 Landset 对应的数字为 7,接下来,我们运行 6s 软件。打开 中的 ,Geometrical conditions (几何条件)igeom 0-7:7(因为IGEOM,从对应上我们找到TM影像,即 Landset对应的数字为7);输好后只按一次enter键;4、接下来输入时间:月,日,小时,经度,纬度。5、Atmospheric model(大气模型)idatm 0-8: 2(中纬度夏季大气模式)6、Aerosol model(type)(气溶胶类型) iae
4、r0-12:37、 Aerosol model (concentration)(气溶胶浓度)the visibility:16 (输入光学厚度)8、The altitude of target ( 目标的高度)xps: -0.02 (目标高度0.02km)9、The sensor altitude ( 传感器的高度)xpp: -1000 (传感器在卫星高度)10、The spectral conditions ( 波段状况)iwave25-30:27 (25即TM1,26即TM2,27即TM3,28即TM4,29即TM5,30即TM7。但是tm6是热红外波段,不适用于我们的此时地表反射率的计算
5、。)11、 Ground reflectance (type)(地面反射率)inhomo: 0 (地面均匀)idirec: 0 (地面无方向影响)igroun: 0 (反射率不随波长变化)ro: 0 (地面波段反射率0)12、 Atmospheric correction mode (激活大气订正的方式)rapp: -0.2 (表观反射率)13、完成后将形成的文件 重命名为 ,波段4的操作依旧如此,只是在输入The spectral conditions ( 波段状况)时改为28即可,得到。14、打开 文件,根据这两个公式y=xa*(measured radiance)-xb; acr=y/(
6、1.+xc*y) ,输入的是影像所测的辐射亮度Radiance,得到的是反射率影像。将这两个公式合并为一个公式:y=xa*(measured radiance)-xb;acr=y/(1+xc*y)合并后:acr=xa*(measured radiance)-xb/1+xc*xa*(measured radiance)-xb15、打开envi的在波段计算器中输入(0.00336*b1-0.04143)/1+0.07179*(0.00336*b1-0.04143)并保存。 点“OK” 。选 3波段,形成一个新文件 ,选择 TM543 裁剪后的波段 3,保存,并以 的形式打开,得到我们需要的影像16、同样处理波段 4。打开 ,文件找到我们需要的参数选 4 波段。形成 TM4 公式:(0.00470*b1-0.01979)/1+0.04569*(0.00470*b1-0.01979)带入波段计算器:选中波段 4,重命名 ,并保存。打开: