[发明专利]遥感影像地形标准化方法有效

专利信息
申请号: 201510211271.6 申请日: 2015-04-29
公开(公告)号: CN104834814B 公开(公告)日: 2018-04-10
发明(设计)人: 张彦丽;李新;闻建光;柳钦火;阎广建;赵军 申请(专利权)人: 西北师范大学
主分类号: G06F19/00 分类号: G06F19/00
代理公司: 甘肃省知识产权事务中心62100 代理人: 周立新
地址: 730070 甘肃*** 国省代码: 甘肃;62
权利要求书: 查看更多 说明书: 查看更多
摘要:
搜索关键词: 遥感 影像 地形 标准化 方法
【权利要求书】:

1.一种遥感影像地形标准化方法,其特征在于,该方法具体为:

1)卫星传感器获得的大气顶辐射亮度LTOP(λ)由路径辐射Lp(λ)和来自地表反射的辐射L(λ),即:

LTOP(λ)= Lp(λ)+ L(λ)T(λ,θν) (1)

在山区,来自地表反射辐射为:

L(λ)=(1/π)(Edir(λ)+ Eaniso-dif(λ))×ρT(λ)(is,iν,Φs-ν)+ (1/π) (Eiso-dif(λ)+ Eref(λ))ρTsd(λ) (iν,Φν)(2)

式中:Edir(λ)、Eiso-dif(λ)、Eaniso-dif(λ)和Eref(λ)分别表示太阳直接辐射、各向同性散射辐射、各向异性散射辐射和来自周围地形的反射辐射;ρT(λ)和ρTsd(λ)分别表示坡地方向-方向反射率和半球-方向反射率;

对于相同地表覆盖类型而言,坡地与消除地形的平坦地表方向-方向反射率不同,其关系为:

ρT(λ)( is, ivs-v)=Ω(λ) ( is, ivs-v,θs,θvs-v) ×ρH(λ) (θs,θvs-v)(5)

式中,Ω(λ) ( is, ivs-v,θs,θvs-v)为坡地BRDF归一化方向反射因子;由于半球-方向反射率比较复杂,同时坡地各向同性散射辐射和来自周围地形反射辐射的能量较小,因此不区分坡地和平坦地表半球-方向反射率,即:

ρTsd(λ)(iv,Φv)=ρH(λ) (θs,θvs-v)(6)

因此,将公式(2)、公式(5)和公式(6)代入公式(1),地形标准化后的反射率为:

ρH(λ) (θs,θvs-v)=π(LTOP(λ)- Lp(λ))/T(λ,θv)〔(Edir(λ)+ Eaniso-dif(λ)) Ω(λ)+Eiso-dif(λ) + Eref(λ)〕(7)

当Ω(λ)=1时,地表假设为朗伯体反射特征;存在两种情况,一种是对遥感影像只进行大气校正,公式可变为:

ρH(λ) (θs,θvs-v)=π(LTOP(λ)- Lp(λ))/T(λ,θv)E0(λ)D0T(λ,θs) (8)

另外一种情况是对遥感影像同时进行大气校正与地形校正,即地形标准化,但只是去除地形对太阳照射角的影响,则地表反射率反演公式为:

ρH(λ) (θs,θvs-v)=π(LTOP(λ)- Lp(λ))/T(λ,θv)E(λ) (9)

2)大气透过率受气溶胶、水汽含量的影响较大,同时也受瑞利散射和臭氧吸收的影响,即:

T(λ,θs)= To3(λ,θs)Tw(λ,θs) Tr(λ,θs)Ta(λ,θs) (10)

式中:T(λ,θs)、To3(λ,θs)、Tw(λ,θs)、Tr(λ,θs)和Ta(λ,θs)分别表示总的大气波谱透过率、臭氧透过率、水汽透过率、瑞利散射透过率和气溶胶透过率, 由于臭氧吸收与瑞利散射随时空分布变化不大,臭氧光学厚度可采用多年平均的空间分布数据,瑞利散射可利用经验公式估算;然而在山区,气溶胶和水汽分布随时空变化较为强烈,水汽和气溶胶波谱透过率公式如下:

式中,k(λ)、mw分别表示水汽有效吸收系数、大气质量及大气可降水厚度;αβ是angstrom波长指数与浑浊度;对于大气可降水厚度,MODIS水汽产品PW直接作为水汽透过率输入参数;AOD产品中的波长指数α也作为气溶胶透过率的输入参数;β表示为:

Β=τ0.55(0.55)a (14)

在山区,总的太阳波谱辐射由太阳直接辐射、太阳散射和周围地形反射辐射三部分组成,均为波长的函数;考虑到太阳方向对散射辐射的影响,将其又分解为各向同性散射辐射和各向异性散射辐射,即:

E(λ)= Edir(λ)+Eiso-dif(λ)+Eaniso-dif(λ)+Eref(λ)(15)

对于Landsat TM每个波段的太阳直接辐射与散射辐射同时受大气与地形作用的影响;

周围地形反射辐射估算方法如下:

式中,Fij为地形结构因子,表示来自周围可见像元的反射辐射能量与可见像元自身总的太阳波谱辐射; Li(λ)是周围可见像元的辐射亮度,由Landsat TM影像测量而得;

3)坡地BRDF 归一化方向反射因子是遥感影像地形标准化的重要参数,是坡地与消除地形的平坦地表方向-方向反射率之比,而这两种方向反射率均遵循地表BRDF规律;半经验线性核驱动模型由一定物理意义核的线性组合来拟合地表二向反射特性,其公式如下:

ρH(λ) (θs,θvs-v)=fiso(λ)+ fvol(λ)Kvol(θs,θvs-v)+ fgeo(λ)Kgeo(θs,θvs-v)(17)

式中,Kvol(θs,θvs-v)和Kgeo(θs,θvs-v)表示核函数,是太阳照射角与传感器观测角度的函数, fiso(λ)、fvol(λ) 和fgeo(λ)表示BRDF模型各向同性、体散射与几何光学核函数的权重系数;

以坡元实际太阳照射角代替太阳天顶角,以坡元传感器实际观测角代替传感器天顶角,得到坡元在传感器过境时刻的反射率:

ΡT(λ) (is, iv,Φs-v)=fiso(λ)+ fvol(λ)Kvol(is, iv,Φs-v)+ fgeo(λ)Kgeo(is, iv,Φs-v)(18)

其中,Φs-v=∣Φs -Φv∣,

式中:Φs-v为太阳与传感器相对方位角,Φs表示由于地形影响坡元上实际来自太阳的照射角,Φv表示由于地形影响坡元上实际来自传感器方向的观测角。

下载完整专利技术内容需要扣除积分,VIP会员可以免费下载。

该专利技术资料仅供研究查看技术是否侵权等信息,商用须获得专利权人授权。该专利全部权利属于西北师范大学,未经西北师范大学许可,擅自商用是侵权行为。如果您想购买此专利、获得商业授权和技术合作,请联系【客服

本文链接:http://www.vipzhuanli.com/pat/books/201510211271.6/1.html,转载请声明来源钻瓜专利网。

×

专利文献下载

说明:

1、专利原文基于中国国家知识产权局专利说明书;

2、支持发明专利 、实用新型专利、外观设计专利(升级中);

3、专利数据每周两次同步更新,支持Adobe PDF格式;

4、内容包括专利技术的结构示意图流程工艺图技术构造图

5、已全新升级为极速版,下载速度显著提升!欢迎使用!

请您登陆后,进行下载,点击【登陆】 【注册】

关于我们 寻求报道 投稿须知 广告合作 版权声明 网站地图 友情链接 企业标识 联系我们

钻瓜专利网在线咨询

周一至周五 9:00-18:00

tel code back_top