[发明专利]一种大气校正方法以及系统在审
申请号: | 201610015091.5 | 申请日: | 2016-01-11 |
公开(公告)号: | CN105675016A | 公开(公告)日: | 2016-06-15 |
发明(设计)人: | 王中挺;厉青;王桥;周春艳;陈辉;张玉环;马鹏飞;张丽娟;毛慧琴;陈翠红 | 申请(专利权)人: | 环境保护部卫星环境应用中心 |
主分类号: | G01C25/00 | 分类号: | G01C25/00;G01C11/00 |
代理公司: | 北京路浩知识产权代理有限公司 11002 | 代理人: | 李相雨 |
地址: | 100029*** | 国省代码: | 北京;11 |
权利要求书: | 查看更多 | 说明书: | 查看更多 |
摘要: | |||
搜索关键词: | 一种 大气 校正 方法 以及 系统 | ||
1.一种大气校正方法,其特征在于,所述方法包括:
获取来自高分一号卫星的WFV相机的数据;
根据所述数据计算得到WFV相机各波段像元的表观反射率;
根据所述数据利用辐射传输模式和WFV相机的波段响应函数 进行相关计算,完成WFV相机遥感图像上不同区域大气分子散射差 异修正,得到修正后的大气参数,从而形成大气参数查找表;
根据所述表观反射率和所述大气参数查找表获得WFV相机各 波段像元的地表反射率,从而实现大气校正。
2.根据权利要求1所述的大气校正方法,其特征在于,所述根 据所述数据计算得到WFV相机各波段像元的表观反射率,具体包 括:
将所述数据的几何定位信息进行坐标转换,得到WFV相机各波 段像元的经纬度,并从所述数据的辅助数据中提取卫星过境时间;
根据所述经纬度和所述卫星过境时间通过下式计算得到各波段 像元的太阳天顶角:
cosθs=sin(lat)sinδ+cos(lat)cosδcost
根据所述卫星过境时间分季节选取辐射定标系数;
通过下式将各波段像元的亮度DN值转换为表观辐亮度:
L=DN/g+L0
根据所述太阳天顶角通过下式将所述表观辐亮度归一化为各波 段像元的表观反射率:
其中,lat表示各波段像元的地理纬度,δ表示太阳光与地球赤 道平面的夹角,t表示太阳的时角,θs表示所述太阳天顶角,g、L0分别表示辐射定标系数,DN表示所述亮度DN值,L表示所述表观 辐亮度,Eλ表示大气层顶太阳辐照度;ds表示日地距离修正因子,ρTOA表示所述表观反射率。
3.根据权利要求2所述的大气校正方法,其特征在于,所述方 法还包括:
结合高分一号卫星各WFV相机的联合观测方式,根据所述辅助 数据提供的图像中心点的观测天顶角,各WFV相机观测的倾角以及 各WFV相机的幅宽计算得到WFV相机各波段像元的观测天顶角。
4.根据权利要求3所述的大气校正方法,其特征在于,所述根 据所述数据利用辐射传输模式和WFV相机的波段响应函数进行相 关计算,完成WFV相机遥感图像上不同区域大气分子散射差异修 正,得到修正后的大气参数,从而形成大气参数查找表,具体包括:
根据各波段像元的经纬度从海拔数据库中查找各波段像元的海 拔值,并通过下式计算得到大气分子光学厚度:
τr(λ)=0.008569λ-4(1+0.0113λ-2+0.00013λ-4)e-0.125H
根据所述大气分子光学厚度、太阳天顶角、观测天顶角以及相 对方位角利用辐射传输模式和WFV相机的波段响应函数进行相关 计算,完成WFV相机遥感图像上不同区域大气分子散射差异修正, 得到修正后的大气参数,并存储在文件中,从而形成大气参数查找 表,
其中,λ表示光波波长,H表示所述海拔值,τr(λ)表示大气分 子光学厚度,e=2.718281828459,从所述辅助数据提取得到所述相 对方位角。
该专利技术资料仅供研究查看技术是否侵权等信息,商用须获得专利权人授权。该专利全部权利属于环境保护部卫星环境应用中心,未经环境保护部卫星环境应用中心许可,擅自商用是侵权行为。如果您想购买此专利、获得商业授权和技术合作,请联系【客服】
本文链接:http://www.vipzhuanli.com/pat/books/201610015091.5/1.html,转载请声明来源钻瓜专利网。
- 上一篇:损伤监测复合式压电薄膜传感器
- 下一篇:一种微机械陀螺仪零偏自动消除方法