唔,最近在做作业的时候,一些实验内容涉及到了用ENVI处理遥感图像,然后自己手动操作软件一遍遍的输入各种参数神马的感觉挺无聊。。。。然后决定自己用python里面的opencv库写个脚本批处理图像反射率的计算试试~
核心步骤就是 遥感影像光谱辐射定标 →大气校正→计算反射率这三步了
①、遥感影像的光谱辐射定标
由遥感器的灵敏度特征引起的辐射畸变主要由其光学系统或光电转换系统的特征形成的,光电转换系统的灵敏性特征通常很重复,其校正一般是通过定期的地面测定值进行的。
遥感器光谱辐射定标时采用以下转换算式:
遥感器各波段偏移与增益值从论文找了找后,找到这么一张表~
那么这么个函数就能定标咯:
def computL(gain,Dn,bias):
return (gain*Dn+bias)
②、遥感影像的大气校正
任何一种依赖大气物理模型的大气校正方法都需要先进行遥感器的辐射校准。
公式是这个咯(Chavez P S,Jr. Image -Based Atmospheric Correction Revisited and Improved Photogrammetric Engineering and Remote Sensing, 1996,62,1025 -1036)
其中:Lhazel——大气层光谱辐射值;LI,min——遥感器每一波段最小光谱辐射值;LI,1%——反射率为1%的黑体辐射值。
关于LI,min和LI,1%的计算公式就省略了啊,感兴趣的同学可以自己去查