前期准备：查找关于计算表观反射率的文献，总结出以上的几点：
（1）需要偏移和增益数据，在中国资源环境卫星获取。
需要注意的是，这个定标参数每年都会变化，需要根据数据日期选择，如下图所示
2）计算表观反射率的公式，那些参数需要根据自己使用的数据计算出，不可以直接套用我的计算出来的结果。 3）所有的操作都是在ENVI上的band math的工具上进行的，因此，操作相对来说比较简单。
下面开始进行实践： 使用ENVI5.3软件完成，如上图所示：使用国产的小插件打开原始数据的PMS文件，打开之后如图所示。 查看数据的波段数量，这里需要注意的是，可以看出有5个波段，但是你用到什么波段，就对那个单波段处理，可以节省一大部分时间。
第一个波段是全色的，千万不要搞错。
在本次实验中，我只用到红色波段和近红外波段，因此，我只是用最后两个波段，B3是红包段，B4是近红外。 对红波段进行求辐射亮度值的操作，由表中可知，OFFSET的值都是0，因此，这个值就不加了。
所以，计算算法就是下图所示部 分。一定要选中红色波段输出，如下图所示 命名是有所讲究的，至少让自己或者他人一下子就看懂就行 。输出后图像，上图所示。
同样的对近红外波段是同样的操作处理 ，但是近红外的Gains 与红波段不是一样的，需要看清楚。
同样的，波段选择和命名要注意。
接下来计算表观反射率：输入计算公式（注：这个公式是自己根据影像数据算出来的，不是通用的），选择的波段是上个步骤的辐射亮度值。
波段选择和命名要注意，TOA是表观反射率的缩写！下图是输出结果。 近红外波段也是同样操作。
目前的结果都是单波段，可以使用Layerstacking 工具进行波段合成，使其变成两个波段的数据。
选择两个波段，输出即可。本次实验只是针对两个波段进行，可以使用同样的操作对绿、蓝两个波段进行处理。
原文地址
