地表反射率的计算

地表覆盖反射率的计算(6s 软件的应用) 9 月 23 日首先在 envi 软件中打开已经处理好的真彩色影像（TM543 波段） ，我的影像因为没有居中，所以首 先进行了裁剪，让影像满幅居中再操作。 1、打开遥感影像，并裁减居中：先打开 7 个波段影像，选中 543，。合成，在弹出的对话框中点，按住 ctrl 再选中这 3 项，点 ok，命名为。打开矢量边界，选中。建立掩膜，重命名为，点，形成掩膜文件。再应用掩膜，选文件，点，选，重命名为，形成影像。所以接下来对背景进行裁剪， 选，选，重命名为，形成影像。2、让#1 和 Scroll 中的红方框大致居中，在#1 中任意位置双击弹出“光标位置评估” ，或者右键找出也可。（可是我不太清楚调出这个的目的？） 3、寻找我们应用的黄石市遥感影像中头文件为 MTL.txt 的文件，以写字板的形式打开，方便查看遥感影 像的具体信息。找到影像获取的时间即“DATE _ACQUIRED” ，这个原始的影像获取时间才是我们需要的，不要被其他 的信息误导。 因为 6s 识别不了具体的时分秒，所以我们需要将具体时间换算成小时，即此处的 02:26:32 应转换成 2.43 小时。4、打开中的,我们在运行 6s 的一切操作，都是按着这个步骤来的，但是期间会出现一些专业术语的特定要求，所以我们需要打开另外的文件，书名如下：打开到 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)（气溶胶类型） iaer0-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是热红外 波段，不适用于我们的此时地表反射率的计算。）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/(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，重命名，并保存。打开：