地表覆盖反射率的计算(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)(气溶胶类型)
iaer[0-12]:3
7、 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 (波段状况)
iwave[25-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)-xb]}
15、打开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,重命名,并保存。
打开:。