高分四号数据计算表观反射率

高分四号数据计算表观反射率_第1张图片

前期准备:查找关于计算表观反射率的文献,总结出以上的几点:

(1)需要偏移和增益数据,在中国资源环境卫星获取。

需要注意的是,这个定标参数每年都会变化,需要根据数据日期选择,如下图所示。

高分四号数据计算表观反射率_第2张图片

(2)计算表观反射率的公式,那些参数需要根据自己使用的数据计算出,不可以直接套用我的计算出来的结果。

高分四号数据计算表观反射率_第3张图片

(3)所有的操作都是在ENVI上的band math的工具上进行的,因此,操作相对来说比较简单。

下面开始进行实践:

高分四号数据计算表观反射率_第4张图片

使用ENVI5.3软件完成,如上图所示:使用国产的小插件打开原始数据的PMS文件,打开之后如图所示。

高分四号数据计算表观反射率_第5张图片

查看数据的波段数量,这里需要注意的是,可以看出有5个波段,但是你用到什么波段,就对那个单波段处理,可以节省一大部分时间。

第一个波段是全色的,千万不要搞错。

在本次实验中,我只用到红色波段和近红外波段,因此,我只是用最后两个波段,B3是红包段,B4是近红外。

 

高分四号数据计算表观反射率_第6张图片

对红波段进行求辐射亮度值的操作,由表中可知,OFFSET的值都是0,因此,这个值就不加了。

所以,计算算法就是下图所示部分。

高分四号数据计算表观反射率_第7张图片

一定要选中红色波段输出,如下图所示。

 

高分四号数据计算表观反射率_第8张图片

高分四号数据计算表观反射率_第9张图片

命名是有所讲究的,至少让自己或者他人一下子就看懂就行。

高分四号数据计算表观反射率_第10张图片

输出后图像,上图所示。

同样的对近红外波段是同样的操作处理,但是近红外的Gains与红波段不是一样的,需要看清楚。

高分四号数据计算表观反射率_第11张图片

高分四号数据计算表观反射率_第12张图片

同样的,波段选择和命名要注意。

高分四号数据计算表观反射率_第13张图片

处理结果,如上图所示。

高分四号数据计算表观反射率_第14张图片

接下来计算表观反射率:输入计算公式(注:这个公式是自己根据影像数据算出来的,不是通用的),选择的波段是上个步骤的辐射亮度值。

高分四号数据计算表观反射率_第15张图片

波段选择和命名要注意,TOA是表观反射率的缩写!下图是输出结果。

高分四号数据计算表观反射率_第16张图片

近红外波段也是同样操作。

高分四号数据计算表观反射率_第17张图片

高分四号数据计算表观反射率_第18张图片

目前的结果都是单波段,可以使用Layerstacking 工具进行波段合成,使其变成两个波段的数据。

高分四号数据计算表观反射率_第19张图片

在工具栏里搜索即可。

高分四号数据计算表观反射率_第20张图片

选择两个波段,输出即可。本次实验只是针对两个波段进行,可以使用同样的操作对绿、蓝两个波段进行处理。

高分四号数据计算表观反射率_第21张图片

你可能感兴趣的:(GIS,遥感,GF4,表观反射率)