[发明专利]一种改善业务化核驱动二向性反射分布函数(BRDF)模型热点的方法有效
申请号: | 201310119889.0 | 申请日: | 2013-04-09 |
公开(公告)号: | CN103324827A | 公开(公告)日: | 2013-09-25 |
发明(设计)人: | 焦子锑;董亚冬 | 申请(专利权)人: | 北京师范大学 |
主分类号: | G06F19/00 | 分类号: | G06F19/00 |
代理公司: | 暂无信息 | 代理人: | 暂无信息 |
地址: | 100875 北京*** | 国省代码: | 北京;11 |
权利要求书: | 查看更多 | 说明书: | 查看更多 |
摘要: | 本发明涉及一种改善业务化核驱动二向性反射分布函数(BRDF)模型热点的方法。通过对现有MODIS业务化核驱动BRDF模型RTLSR的体散射核进行热点校正,提出一个考虑热点变化的新体散射核,该核与原几何光学核进行线性组合产生的RTCLSR模型对热点有更好的拟合能力,但不影响其它观测方向的拟合精度,我们用POLDER-3/BRDF数据库和机载CAR数据对RTCLSR模型的热点参数进行了标定和验证,RTCLSR新模型可用于研究热点效应随地表类型,NDVI和太阳天顶角的变化,同时,由于新模型很好保持了原模型的线性形式,没有大量增加模型运算的复杂性,因此,可能为MODIS BRDF/反照率业务化产品升级提供新算法基础和解决方案。本发明在空间信息技术领域,尤其在定量遥感方面具有重要的应用价值。 | ||
搜索关键词: | 一种 改善 业务 驱动 向性 反射 分布 函数 brdf 模型 热点 方法 | ||
【主权项】:
1.一种改善业务化核驱动二向性反射分布函数(BRDF)模型热点的方法,步骤包括: 步骤一:修正核驱动模型的体散射核 线性核驱动的二向性反射函数(BRDF)模型可以表示为: R(θv,θs,Δφ,λ)=fiso(λ)+fvol(λ)Kvol(θv,θs,Δφ)+fgeo(λ)Kgeo(θv,θs,Δφ)(1) 在此,R(θv,θs,Δφ,λ)是二向性反射分布函数,fiso(λ),fvol(λ)和fgeo(λ)是模型的三个参数,分别称为各向同性参数,体散射参数和几何光学参数;Kvol(θv,θs,Δφ)和Kgeo(θv,θs,Δφ)是核驱动模型的体散射和几何光学核函数,给定植被结构参数,两个核函数是观测和太阳位置的三角函数;目前,中分辨率成像光谱仪(MODIS)业务化的二向性反射和反照率产品采用了罗斯表层核(RossThick,KRT)和李氏稀疏互易核(LiSparseR,KLsR)线性复合的RTLSR模型,该模型的主要缺点是低估热点方向的反射率,这是由于罗斯表层核在热点方向上没有考虑光线和视线重合时的相关性问题,罗斯表层核表示如下:
该申请上述的罗斯表层核式(2)进行修正,通过增加了一个热点核,完善它在热点方向的拟合能力,构造了所谓的罗斯表层陈核(RossThickChen,KRTc),表示如下:
式(3)与式(2)相比:(i)增加了热点核
该热点核新增C1/C2两个参数分别用于调节热点高度和宽度,这两个参数需要提前标定;(ii)常数因子-π/4调整为-π/2,使罗斯表层陈核沿纵坐标下移π/4,该调整不会影响体散射核对多角度观测数据总体拟合能力,因为核的拟合能力只与核的形状有关,而与核的位置无关;同时,该调整可使模型各向同性参数(fiso)保留其物理意义,即,各向同性参数可作为太阳和观测在天顶方向时的反射率值。步骤二:标定热点核参数C1/C2对于热点核的两个新增参数C1/C2需要提前标定,由于大多数极轨卫星如,MODIS,MISR等很少获取热点方向的观测数据,所以我们用现有的POLDER-3/BRDF数据库进行模型参数标定,该数据库提供较为完整的热点观测,同时,这些观测按照已有地表生态类型(IGBP)进行分类整理,便于研究热点随地表类型的变化情况,修订后的模型我们称之为RTCLSR模型,可根据热点随IGBP类或NDVI变化所确定的C1/C2参数,用于MODIS现有BRDF和反照率业务化产品生产,标定时我们采用热点区域局部拟合误差最小原则,最小拟合误差公式如下:
其中,Robs(θv,θi,Δφ,λ)和Rmodel(θv,θi,Δφ,λ)分别为观测和模型的BRDF,n为观测数,因修订后的体散射核在热点区域以外的散射方向上和原体体散射核有很好吻合,所以我们取太阳与观测相对方位角ξ≤5°范围的观测数据进行最小拟合误差的计算,我们称之为局部最小拟合误差方法,通过求局部最小拟合误差,我们用迭代的方法求出最小的C1/C2值,在此考虑到C1/C2可能的变化范围和计算的复杂性,取C1=0.3-1.5,C2=0.5-1.5,分别求出整个POLDER-3/BRDF数据库对应的C1/C2值,数据库中每IGBP类(冰雪除外)的C1/C2值,以及不同NDVI分层的C1/C2值,这将有助于进一步分析热点随以上参数的变化; 步骤三:验证RTCLSR新模型 新模型验证分为两步:(i)分析C1/C2参数随地表类型(IGBP)、植被归一化指数(NDVI)和太阳天顶角(SZA)的变化情况;(ii)新模型(RTCLSR)、原模型(RTLSR)以及前人修订模型(RTJLSR)与观测数据的比较差异;基于前面求出的C1/C2值,我们分析了热点高度和宽度随地表类型,NDVI和太阳天顶角变化的规律,间接验证参数化热点的必要性;同时,我们用收集的机载飞行CAR数据和不同地表类型的POLDER观测数据,比较了新模型RTCLSR,原模型RTLSR和前人修订模型RTJLSR在热点方面的差异,对新模型进行间接和直接验证。
下载完整专利技术内容需要扣除积分,VIP会员可以免费下载。
该专利技术资料仅供研究查看技术是否侵权等信息,商用须获得专利权人授权。该专利全部权利属于北京师范大学,未经北京师范大学许可,擅自商用是侵权行为。如果您想购买此专利、获得商业授权和技术合作,请联系【客服】
本文链接:http://www.vipzhuanli.com/patent/201310119889.0/,转载请声明来源钻瓜专利网。
- 同类专利
- 专利分类
G06 计算;推算;计数
G06F 电数字数据处理
G06F19-00 专门适用于特定应用的数字计算或数据处理的设备或方法
G06F19-10 .生物信息学,即计算分子生物学中的遗传或蛋白质相关的数据处理方法或系统
G06F19-12 ..用于系统生物学的建模或仿真,例如:概率模型或动态模型,遗传基因管理网络,蛋白质交互作用网络或新陈代谢作用网络
G06F19-14 ..用于发展或进化的,例如:进化的保存区域决定或进化树结构
G06F19-16 ..用于分子结构的,例如:结构排序,结构或功能关系,蛋白质折叠,结构域拓扑,用结构数据的药靶,涉及二维或三维结构的
G06F19-18 ..用于功能性基因组学或蛋白质组学的,例如:基因型–表型关联,不均衡连接,种群遗传学,结合位置鉴定,变异发生,基因型或染色体组的注释,蛋白质相互作用或蛋白质核酸的相互作用
G06F 电数字数据处理
G06F19-00 专门适用于特定应用的数字计算或数据处理的设备或方法
G06F19-10 .生物信息学,即计算分子生物学中的遗传或蛋白质相关的数据处理方法或系统
G06F19-12 ..用于系统生物学的建模或仿真,例如:概率模型或动态模型,遗传基因管理网络,蛋白质交互作用网络或新陈代谢作用网络
G06F19-14 ..用于发展或进化的,例如:进化的保存区域决定或进化树结构
G06F19-16 ..用于分子结构的,例如:结构排序,结构或功能关系,蛋白质折叠,结构域拓扑,用结构数据的药靶,涉及二维或三维结构的
G06F19-18 ..用于功能性基因组学或蛋白质组学的,例如:基因型–表型关联,不均衡连接,种群遗传学,结合位置鉴定,变异发生,基因型或染色体组的注释,蛋白质相互作用或蛋白质核酸的相互作用