本发明专利技术公开了一种地表反照率产品的生成方法及系统,涉及卫星遥感领域。所述方法包括步骤:根据经过分类的POLDER?BRDF数据,得到POLDER反照率数据,转换POLDER反照率数据得到宽波段地表反照率数据;根据样本MODIS地表反射率数据和宽波段地表反照率数据建立第一转换模型,得到第一地表反照率数据;根据样本MODIS表观反射率数据和宽波段地表反照率数据建立第二转换模型,得到第二地表反照率数据;整合第一和第二地表反照率数据,及经网格划分的MODIS地表反照率数据,生成时空连续的地表反照率产品。本发明专利技术生产出了长时间序列并且时空连续、高时空分辨率的地表反照率产品,提高了地表反照率产品反演的精度。
【技术实现步骤摘要】
本专利技术涉及卫星遥感
,特别涉及一种地表反照率产品的生成方法及系统。
技术介绍
地表反照率反映了地球表面对太阳光辐射的反射能力,是地表辐射能量平衡以及地气相互作用中的驱动因子之一,其时空变化受到自然过程以及人类活动的影响,是全球环境变化的指示因子。目前一些卫星反照率产品已经业务化运行和发布,空间分辨率从250m_20km,时间分辨率从日到月,其中极轨卫星产品包括MODIS,MISR,CERES, POLDER, MERIS0地球同步卫星产品包括 Meteosat, and MSG。目前针对在轨运行的卫星遥感数据以及部分航空遥感数据进行了大量地表反照率反演研究,并进行了有效的反照率产品的地面验证,其中基于线性核驱动模型的地表 BRDF/反照率遥感反演模型算法是目前地表反照率遥感反演中应用最广泛的方法,已经在 M0DIS、MISR的地表反照率产品中得到了很好的应用。另外,POLDER/PARASOL系列传感器具有更好的多角度观测能力,空间分辨率稍低(6km),也发布了非常有特色的全球长时间序列反照率产品。其他如静止轨道气象卫星MSG、METE0SET,极轨卫星传感器AVHRR、VEGETATION 等传感器都有不同覆盖范围的反照率产品,我国气象卫星、环境减灾小卫星数据都可用于区域和全球范围的地表反照率反演,目前反照率产品生成算法正在研制过程中。已有的反照率产品的生产系统环境,均是采用个人计算机少量生产,并且缺少长时间序列的全球陆表宽波段发射率数据集。迄今为止,针对全球陆面变化研究与陆面模型研发,国际陆地遥感领域仍然缺乏长时间序列并且时空连续、高时空分辨率的全球陆表特征参量产品。国内的遥感产品生产均采用个人计算机生产小批量产品,无法满足长时间序列、高时空分辨率和高质量的遥感生产需求。
技术实现思路
(一)要解决的技术问题本专利技术要解决的技术问题是如何提供一种时空连续的地表反照率产品的生成方法及系统。( 二 )技术方案为解决上述技术问题,本专利技术提供一种地表反照率产品的生成方法,其包括步骤B 根据太阳天顶角、观测天顶角和相对方位角对MODIS地表反射率数据、MODIS表观反射率数据和MODIS地表反照率数据进行网格划分;C 对POLDER BRDF数据进行优选,按照地表类型对经过优选的POLDER BRDF数据进行分类;D 根据配套光谱数据模拟得到模拟POLDER BRDF数据和模拟MODIS地表反射率数据,建立所述模拟POLDER BRDF数据和所述模拟MODIS地表反射率数据之间的多元线性转换关系,根据所述多元线性转换关系和经过分类的POLDER BRDF数据,计算得到样本MODIS 地表反射率数据;E 根据经过分类的POLDER BRDF数据,得到POLDER反照率数据,转换所述POLDER 反照率数据得到宽波段地表反照率数据;F 根据所述样本MODIS地表反射率数据和所述宽波段地表反照率数据建立第一转换模型,根据所述第一转换模型和经过网格划分的MODIS地表反射率数据,得到第一地表反照率数据;G 根据所述POLDER BRDF数据得到地表反射特性参数,根据所述地表反射特性参数和大气状态参数计算得到样本MODIS表观反射率数据;H:根据所述样本MODIS表观反射率数据和所述宽波段地表反照率数据建立第二转换模型,根据所述第二转换模型和经过网格划分的MODIS表观反射率数据,得到第二地表反照率数据;I 整合所述第一地表反照率数据、所述第二地表反照率数据和经过网格划分的 MODIS地表反照率数据,生成时空连续的地表反照率产品。优选地,在所述步骤B之前还包括步骤A 输入POLDER BRDF数据、MODIS地表反射率数据、MODIS表观反射率数据、MODIS地表反照率数据和配套光谱数据。优选地,所述步骤I具体包括步骤Il 对所述第一地表反照率数据和所述第二地表反照率数据进行插值弥补;12:对所述第一地表反照率数据、所述第二地表反照率数据和经过网格划分的 MODIS地表反照率数据进行归一化处理;13 对归一化处理后的第一地表反照率数据、第二地表反照率数据和经过网格划分的MODIS地表反照率数据进行数据融合,得到融合数据;14 对所述融合数据进行平滑处理;15 对平滑处理后的融合数据进行反归一化处理,得到所述时空连续的地表反照率产品。优选地,所述反归一化处理的公式如下ξ = ξσ-^ι其中,ξ是所述平滑处理后的融合数据,ξ是所述时空连续的地表反照率产品,;; 是归一化处理后的第一地表反照率数据、第二地表反照率数据和经过网格划分的MODIS地表反照率数据的均值,G是归一化处理后的第一地表反照率数据、第二地表反照率数据和经过网格划分的MODIS地表反照率数据的消除随机误差后的标准差的平均值。优选地,所述步骤C具体包括步骤Cl 从所述POLDER BRDF数据中剔除由于云和气溶胶影响以及地表状态发生变化而不满足二向反射模型的数据;C2 对经过所述步骤Cl处理后的POLDER BRDF数据,按照植被、冰雪和裸地三种地表类型进行分类。 优选地,所述步骤C2具体包括步骤C21 判断POLDER BRDF数据的NDVI值是否大于0. 2,如果是,则判定所述POLDER BRDF数据对应的像元为植被,否则执行步骤C22 ;C22 判断POLDER BRDF数据是否其蓝光波段反射率大于0. 3或者红光波段反射率大于0. 3,如果是,则判定所述POLDER BRDF数据对应的像元为冰雪,否则,判定所述POLDER BRDF数据对应的像元为裸地。优选地,所述步骤E具体包括步骤El 对所述经过分类的POLDER BRDF数据进行半球积分,得到各波段的POLDER反照率数据;E2:根据窄波段向宽波段的转换公式,转换所述POLDER反照率数据得到所述宽波段地表反照率数据。本专利技术还提供一种地表反照率产品的生成系统,其包括网格划分模块用于根据太阳天顶角、观测天顶角和相对方位角对MODIS地表反射率数据、MODIS表观反射率数据和MODIS地表反照率数据进行网格划分;POLDER BRDF数据预处理模块用于对POLDER BRDF数据进行优选,按照地表类型对经过优选的POLDER BRDF数据进行分类;样本MODIS地表模块用于根据配套光谱数据模拟得到模拟POLDER BRDF数据和模拟MODIS地表反射率数据,建立所述模拟POLDER BRDF数据和所述模拟MODIS地表反射率数据之间的多元线性转换关系,根据所述多元线性转换关系和经过分类的P0LDERBRDF数据,计算得到样本MODIS地表反射率数据;POLDER转换模块用于根据经过分类的POLDER BRDF数据,得到POLDER反照率数据,转换所述POLDER反照率数据得到宽波段地表反照率数据;第一产品模块用于根据所述样本MODIS地表反射率数据和所述宽波段地表反照率数据建立第一转换模型,根据所述第一转换模型和经过网格划分的MODIS地表反射率数据,得到第一地表反照率数据;样本MODIS表观模块用于根据所述POLDER BRDF数据得到地表反射特性参数,根据所述地表反射特性参数和大气状态参数计算得到样本MODIS表观反射率数据;第二产品模块用于根据所述样本MODIS表观反射率数据本文档来自技高网...
【技术保护点】
【技术特征摘要】
【专利技术属性】
技术研发人员:赵祥,刘强,王立钊,张鑫,瞿瑛,梁顺林,
申请(专利权)人:北京师范大学,
类型:发明
国别省市:
还没有人留言评论。发表了对其他浏览者有用的留言会获得科技券。