地表反射率的计算.docx

《地表反射率的计算.docx》由会员上传分享，免费在线阅读，更多相关内容在学术论文-天天文库。

1、地表覆盖反射率的计算(6s软件的应用)9月23日首先在envi软件中打开已经处理好的真彩色影像（TM543波段），我的影像因为没有居中，所以首先进行了裁剪，让影像满幅居中再操作。1、打开遥感影像，并裁减居中：先打开7个波段影像，，选中543，。合成，，，在弹出的对话框中点，按住ctrl再选中这3项，点ok，命名为。打开矢量边界，，，选中。建立掩膜，，，，，，重命名为，点，形成掩膜文件。再应用掩膜，选文件，点，选，，，重命名为，形成影像。所以接下来对背景进行裁剪，，，选，，，选，，，，，重命名为，形成影像。2、让#1和Scroll中的红方框大致居中，在#1中

2、任意位置双击弹出“光标位置评估”，或者右键找出也可。（可是我不太清楚调出这个的目的？）3、寻找我们应用的黄石市遥感影像中头文件为MTL.txt的文件，以写字板的形式打开，方便查看遥感影像的具体信息。找到影像获取的时间即“DATE_ACQUIRED”，这个原始的影像获取时间才是我们需要的，不要被其他的信息误导。因为6s识别不了具体的时分秒，所以我们需要将具体时间换算成小时，即此处的02:26:32应转换成2.43小时。4、打开中的,我们在运行6s的一切操作，都是按着这个步骤来的，但是期间会出现一些专业术语的特定要求，所以我们需要打开另外的文件，书名如下：打开

3、到35页，IGEOM，从对应上我们找到TM影像，即Landset对应的数字为7，接下来，我们运行6s软件。打开中的，Geometricalconditions(几何条件)igeom[0-7]:7(因为IGEOM，从对应上我们找到TM影像，即Landset对应的数字为7);输好后只按一次enter键；4、接下来输入时间：月，日，小时，经度，纬度。5、Atmosphericmodel（大气模型）idatm[0-8]:2（中纬度夏季大气模式）6、Aerosolmodel(type)（气溶胶类型）iaer[0-12]:37、Aerosolmodel(concent

4、ration)（气溶胶浓度）thevisibility:16(输入光学厚度)8、Thealtitudeoftarget（目标的高度）xps:-0.02(目标高度＝0.02km)9、Thesensoraltitude（传感器的高度）xpp:-1000(传感器在卫星高度)10、Thespectralconditions（波段状况）iwave[25-30]:27（25即TM1,26即TM2，27即TM3，28即TM4，29即TM5,30即TM7。但是tm6是热红外波段，不适用于我们的此时地表反射率的计算。）11、Groundreflectance(type)（地

5、面反射率）inhomo:0(地面均匀)idirec:0(地面无方向影响)igroun:0(反射率不随波长变化)ro:0(地面波段反射率＝0)12、Atmosphericcorrectionmode（激活大气订正的方式）rapp:-0.2（表观反射率)13、完成后将形成的文件重命名为，波段4的操作依旧如此，只是在输入Thespectralconditions（波段状况）时改为28即可，得到。14、打开文件，根据这两个公式y=xa*(measuredradiance)-xb;acr=y/(1.+xc*y),输入的是影像所测的辐射亮度Radiance，得到的是反

6、射率影像。将这两个公式合并为一个公式：y=xa*(measuredradiance)-xb;acr=y/(1+xc*y)合并后：acr=[xa*(measuredradiance)-xb]/{1+xc*[xa*(measuredradiance)-xb]}15、打开envi的在波段计算器中输入(0.00336*b1-0.04143)/[1+0.07179*(0.00336*b1-0.04143)]并保存。点“OK”。选3波段，形成一个新文件，选择TM543裁剪后的波段3，保存，，并以的形式打开，得到我们需要的影像16、同样处理波段4。打开，文件找到我们需要

7、的参数选4波段。形成TM4公式：(0.00470*b1-0.01979)/[1+0.04569*(0.00470*b1-0.01979)]带入波段计算器：选中波段4，重命名，并保存。打开：