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

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

（1）需要偏移和增益数据，在中国资源环境卫星获取。

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

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

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

下面开始进行实践：
使用ENVI5.3软件完成，如上图所示：使用国产的小插件打开原始数据的PMS文件，打开之后如图所示。
查看数据的波段数量，这里需要注意的是，可以看出有5个波段，但是你用到什么波段，就对那个单波段处理，可以节省一大部分时间。

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

在本次实验中，我只用到红色波段和近红外波段，因此，我只是用最后两个波段，B3是红包段，B4是近红外。
对红波段进行求辐射亮度值的操作，由表中可知，OFFSET的值都是0，因此，这个值就不加了。

所以，计算算法就是下图所示部
分。一定要选中红色波段输出，如下图所示

命名是有所讲究的，至少让自己或者他人一下子就看懂就行
。输出后图像，上图所示。

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

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

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

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

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

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

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

原文地址