顾及实际天空光分布影响的地表反照率遥感模型的制作方法
【专利摘要】本发明公开了一种顾及实际天空光分布影响的地表反照率遥感模型，在构建过程中考虑了大气条件下天空光角度分布的影响，以及地表与大气之间的多次反射机制。模型用户需要通过所获得的空间传感器观测到的关于地面目标的多角度反射率数据，反演地表双向反射分布函数BRDF模型；利用太阳入射与传感器观测之间的几何角度关系计算BRDF模型核函数；通过实地测量或卫星观测反演得到的大气气溶胶的光学特性参数以及水汽、臭氧等成分的含量，实施地?气耦合辐射传输计算获得下行辐照度场中的漫总比参数、天空光角度分布函数、大气球面反照率等地表反照率遥感模型所需输入参数，同时计算直射与漫射反照率以及天空光角度分布影响下的漫射反照率。
【专利说明】
顾及实际天空光分布影响的地表反照率遥感模型
技术领域
[0001] 本发明涉及遥感科学地表反射特性提取建模领域，具体是一种顾及实际天空光分 布影响的地表反照率遥感模型。
【背景技术】
[0002] 空间光学载荷对地观测时，地表和大气的复合影响必须予以考虑。大气中关键成 分与粒子因其具有各向异性的散射相函数，使得大气层内散射辐射场具有非均一的空间分 布特性，即"天空光"的概念。一方面，地表反射改变了太阳入射辐射的重新分配，大多数的 自然地表均具有明显的方向性反射辐射特征，必须通过地表双向反射分布函数BRDF予以描 述。而在空间某一高度观测地面目标时，除了地表类型、材质、结构和粗糙度等因素的影响， "太阳-地表-观测"之间的相对几何位置关系是驱动BRDF模型计算结果的最主要影响因子。 "直射反照率"是通过获取的BRDF分别在观测空间进行角度积分，表征了绝对晴空条件下的 地表反射特性参量;而"漫射反照率"需要进一步的在入射空间进行角度积分，表征了阴天 条件下的地表反射特性参量。因此，掌握地表的方向性反射分布特征是顾及实际天空光分 布影响的地表反照率遥感模型构建的第一步。另一方面，天空光主要是大气分子和大气气 溶胶颗粒散射太阳辐射，以及地表和大气散射粒子之间的多次相互作用共同导致的结果， 它既包含了大气组分和状态的信息，也包含了地表的信息，同时还包含有地表与大气之间 的耦合信息。对这一物理量及其角度分布的研究将有助于更进一步的了解地-气系统相互 作用机制。
[0003] 地表反照率的计算大多针对太阳反射(0.4~2.5μπι)波段。目前国际上使用的地表 反照率计算式仍然依据基于天空光各向同性假设前提下得到的地表反照率模型，如美国 NASA MODIS BRDF研究小组建立的地表反照率(又称"蓝空"地表反照率)模型，该模型忽略 了现实大气条件下散射福射场的空间分布，并认为任何大气条件下天空光的分布可以用 "各向同性"予以近似。理论计算研究结果表明，这一模型只能用在太阳高度角大于40°时的 情况，而当太阳高度角低于40°时，大气浑浊度是影响模型精度的最主要因素。
【发明内容】
[0004] 本发明的目的是提供一种顾及实际天空光分布影响的地表反照率遥感模型，以解 决现有地表反射特性定量提取中的精度和应用局限性的问题。
[0005] 为了达到上述目的，本发明所采用的技术方案为：
[0006] 顾及实际天空光分布影响的地表反照率遥感模型，其特征在于：通过多角度卫星 观测地表目标的遥感数据，引入考虑遥感过程中存在的热点效应即后向反射峰值的Ross-Li-Maignan地表反射模型，简称RLM模型，并获得RLM模型的核系数；
[0007] 利用太阳入射角度、传感器观测角度计算地表朗伯散射核、Ross体散射核和Li-Sparse互易散射核函数的值，通过RLM模型计算任意入射角度条件下的地表双向反射比因 子BRF，再将得到的双向反射分布函数BRDF，BRDF = π-1 X BRF，分别通过视空间角度积分运算 获得直射反照率；
[0008]利用观测时刻获取的气溶胶和水汽等大气关键成分的光学特性参数，运行辐射传 输模型，获得太阳入射辐射场中的漫总比参数，并构建天空光角度分布函数；
[0009]通过将直射反照率在入射空间进行角度积分，并利用已获得的天空光角度分布函 数作为权重，在所有积分角度对直射反照率进行积分计算，从而获得实际天空光分布影响 下的漫射反照率；
[0010] 所构建的顾及实际天空光分布影响的地表反照率遥感模型如公式(1):
[0011]
[0012] 模型的公式（1)中，L表示辐射亮度，E表示辐射照度或辐射通量密度，符号丨/丨表示 福射传输的方向，分别对应向上/向下，下标ms表不多次散射，0 3和(J>s分别为太阳天顶角和 方位角，θν和φν则分别表示观测天顶角和方位角，f dt(A，T)表示下行辐射漫总比，adir(0s，A) 和<1-_ (：^(^，'〇分别表示直射和漫射反照率，6=(1-(1￡^_(：^€^ 111广1为地-气系统相互作用增 益因子，Gd和Gu分别表示地-气交界处向下和向上的增益因子，λ为传感器观测通道中心波 长，公式（1)中，a anidlf(A,T)用公式(2)计算：
[0013]
[0014] 公式（2)中，-为天空光角度分布 函数；
[0015] 模型的公式（1)在理论推导中运用了下面的关系式：
[0016]
[0017] 其中大气球面反照率的计算式如下：
[0018]
[0019] 公式⑷中，QS=(0S，(K)与Ων=(θ ν，φν)分别表示太阳与观测位置，P为地表高 度处的大气压强，Rrayl与Rrayl+aer。分别表不瑞利大气、瑞利分子和气溶胶混合大气的固有发 射率，为水汽柱含量，μ表示角度余弦值，至此便实现了高精度地表谱反照率的理论计 算和遥感模型的构建。
[0020] 所述的顾及实际天空光分布影响的地表反照率遥感模型，其特征在于:考虑了实 际天空光的分布特点。
[0021] 所述的顾及实际天空光分布影响的地表反照率遥感模型，其特征在于:考虑了实 际大气与地表的多次相互作用对地表水平处向下的辐射照度场和向上的反射辐射场的影 响机制。
[0022] 所述的顾及实际天空光分布影响的地表反照率遥感模型，其特征在于：可应用于 非均一分布的下垫面地表。
[0023] 所述的顾及实际天空光分布影响的地表反照率遥感模型，其特征在于:包含了太 阳入射角度和大气浑浊度这两个最关键的影响因素，所述模型对大角度入照条件和较浑浊 大气条件下地表反照率的高精度遥感具有技术优势。
[0024] 所述的顾及实际天空光分布影响的地表反照率遥感模型，其特征在于:有望解决 高炜度地区尤其是冬春季节地表反射特性精度遥感不足的问题。
[0025] 所述的顾及实际天空光分布影响的地表反照率遥感模型，其特征在于:为任何涉 及地表-大气耦合辐射传输的高精度计算和特征目标的辐射修正提供唯一下垫面地表反射 特性接口参数支持。
[0026] 所述的顾及实际天空光分布影响的地表反照率遥感模型，其特征在于:为提高卫 星遥感反照率产品与地基测量数据的一致性、减小数据交叉比对误差提供了桥梁衔接的作 用，有效支持空间高分辨卫星遥感数据的地面真实性检验和临近空间光电系统的实时标 校。
[0027] 本发明的原理是:依据大气散射辐射具有各向异性的本质特征，以及地表双向反 射BRDF理论，通过引入天空光角度分布函数来刻画下行漫射辐射在积分空间内对地表BRDF 的影响，通过漫总比来描述下行辐照度场中漫射成分所占比例，通过大气球面反照率、地表 直射反照率和漫射反照率来共同描述辐射在地-气系统中的多次相互作用机制。
[0028] 本发明研发了一种新型的顾及天空光实际分布影响的高精度地表反照率遥感模 型，以弥补现有模型在区域、季节适用性以及计算精度等方面的不足，更好地满足气候研究 对陆表参数遥感产品的应用需求，同时也为卫星陆表产品的地面真实性检验提供通用模型 支持。
[0029] 与现有技术相比本发明具有的积极效果在于：
[0030] 1)由于引入了天空光角度分布函数，因此对实际大气的区域适用性更加自由，既 可以应用于中低炜度地区，也可以应用于高炜度地区；
[0031] 2)由于构建的地表反照率遥感模型考虑了太阳位置和大气状况等主要参数，对于 高炜度尤其是冬春季节地表特性的定量提取，具有精度上的技术优势；
[0032] 3)由于使用了地-气耦合系统辐射传输理论，综合考虑了地表与大气之间的多次 反射机制，发现了太阳高度角在小于40°时，天空散射光对BRDF模型的体散射核和几何光学 散射核的影响规律。这一发现为临近空间对地观测目标的方向性校正和辐射校正提供了新 的视角，起到了积极的指导意义。
【附图说明】
[0033]图1是本发明的模块结构图。
[0034]图2是本发明的技术流程。
[0035] 图3是RLM BRDF模型积分核查找表图。
【具体实施方式】
[0036] 顾及实际天空光分布影响的地表反照率遥感模型，通过多角度卫星观测地表目标 的遥感数据，引入考虑遥感过程中存在的热点效应即后向反射峰值的R〇ss-Li-Maignan地 表反射模型，简称RLM模型，并获得RLM模型的核系数；
[0037] 利用太阳入射角度、传感器观测角度计算地表朗伯散射核、Ross体散射核和Li-Sparse互易散射核函数的值，通过RLM模型计算任意入射角度条件下的地表双向反射比因 子BRF，再将得到的双向反射分布函数BRDF，BRDF = jT1 X BRF，分别通过视空间角度积分运算 获得直射反照率；
[0038] 利用观测时刻获取的气溶胶和水汽等大气关键成分的光学特性参数，运行辐射传 输模型，获得太阳入射辐射场中的漫总比参数，并构建天空光角度分布函数；
[0039] 通过将直射反照率在入射空间进行角度积分，并利用已获得的天空光角度分布函 数作为权重，在所有积分角度对直射反照率进行积分计算，从而获得实际天空光分布影响 下的漫射反照率；
[0040] 所构建的顾及实际天空光分布影响的地表反照率遥感模型如公式(1):
[0041]
[0042] 模型的公式（1)中，L表示辐射亮度，E表示辐射照度或辐射通量密度，符号丨/丨表示 福射传输的方向，分别对应向上/向下，下标ms表不多次散射，0 3和(J>s分别为太阳天顶角和 方位角，θν和φν则分别表示观测天顶角和方位角，f dt(A，T)表示下行辐射漫总比，adir(0s，A) 和<1-_ (：^(^，'〇分别表示直射和漫射反照率，6=(1-(1￡^_(：^€^ 111广1为地-气系统相互作用增 益因子，Gd和Gu分别表示地-气交界处向下和向上的增益因子，λ为传感器观测通道中心波 长，公式(1)中，a ani_dlf(A，T)用公式(2)计算：
[0043] /(k
[0044] 公式（2)中，
％天空光角度分布 函数；
[0045] 模型的公式(1)在理论推导中运用了下面的关系式：
[0046]
[0047] 其中大气球面反照率的计算式如下：
[0048]
[0049] 公式(4)中，QS=(0S，(K)与Ων=(θ ν，φν)分别表示太阳与观测位置，P为地表高 度处的大气压强，Rrayl与Rrayl+aer。分别表不瑞利大气、瑞利分子和气溶胶混合大气的固有发 射率，CLg为水汽柱含量，μ表示角度余弦值，至此便实现了高精度地表谱反照率的理论计 算和遥感模型的构建。
[0050] 本发明顾及实际天空光分布影响的地表反照率遥感模型，考虑了实际天空光的分 布特点、实际大气与地表的多次相互作用对地表水平处向下的辐射照度场和向上的反射辐 射场的影响机制，可应用于非均一分布的下垫面地表。该模型，包含了太阳入射角度和大气 浑浊度这两个最关键的影响因素，模型对大角度入照条件和较浑浊大气条件下地表反照率 的高精度遥感具有技术优势。
[0051] 本发明有望解决高炜度地区尤其是冬春季节地表反射特性精度遥感不足的问题。 为任何涉及地表-大气耦合辐射传输的高精度计算和特征目标的辐射修正提供唯一下垫面 地表反射特性接口参数支持。为提高卫星遥感反照率产品与地基测量数据的一致性、减小 数据交叉比对误差提供了 "桥梁衔接"的作用，有效支持空间高分辨卫星遥感数据的地面真 实性检验和临近空间光电系统的实时标校。
[0052]参见图1，顾及实际天空光分布影响的地表反照率遥感模型，包括有6个大模块：
[0053]模块1是目标方向性反射率观测数据输入模块，包括近地面或空基光学平台探测 数据的提取与预处理。
[0054]模块2地表目标BRDF计算模块，该模块中涉及到几何角度计算、BRDF模型核函数计 算以及多角度BRDF模型反演方法。
[0055] 模块3是地表反射特性参数计算模块，包含有地表直射反照率、漫射反照率。
[0056] 模块4是大气光学特性参数测量与辐射传输计算，包括气溶胶光学特性参数测量、 水汽与臭氧等大气模式参数测量。
[0057] 模块5利用测量的大气光学数据作为输入参数，驱动大气辐射传输计算模型，得到 天空光角度分布函数、下行辐照度场漫总比和大气球面反照率等关键参数。
[0058] 模块6利用模块3和模块5分别得到的地表反射特性参数和大气辐射场参数，同步 输入至实际天空光分布影响下的漫射反照率计算模型。
[0059] 模块6是本发明的核心部分，亦即前述构建的顾及实际天空光分布影响的地表反 照率遥感模型，本发明在理论和原理上进行了详细推导。通过模块1~模块3的输出数据作 为顾及实际天空光分布影响的地表反照率遥感模型的地表参数输入，通过模块4~模块5的 输出数据作为顾及实际天空光分布影响的地表反照率遥感模型的大气参数输入，从而计算 瞬时地表反照率。
[0060] 本发明主要技术流程图如图2所示。其中，模块3中的中间数据的计算，可通过如图 3所示的查找表快速插值得到。由于前述反照率计算模型包含针对空间角度的四维积分计 算，为了提高顾及实际天空光分布影响的地表反照率遥感模型的计算效率，本发明预先生 成了RLM BRDF模型核函数积分查找表，计算效率提高约2000倍。
[0061]本发明主要实现了顾及实际天空光分布影响的地表反照率遥感模型，高精度的技 术优势主要体现在综合考虑了实际大气条件下天空光可能各向异性的分布特征，以及地表 与大气之间的多次反射机制两大方面的主要因素，高效率主要体现在预先建立的反照率积 分核查找表。该理论模型有助于推进国家空间高分辨卫星陆表参数产品的气候与地表定量 遥感应用水平。
【主权项】
1.顾及实际天空光分布影响的地表反照率遥感模型，其特征在于：通过多角度卫星观 测地表目标的遥感数据，引入考虑遥感过程中存在的热点效应即后向反射峰值的Ross-Li- Maignan地表反射模型，简称RLM模型，并获得RLM模型的核系数； 利用太阳入射角度、传感器观测角度计算地表朗伯散射核、Ross体散射核和Li-Sparse 互易散射核函数的值，通过RLM模型计算任意入射角度条件下的地表双向反射比因子BRF， 再将得到的双向反射分布函数BRDF，B畑F = 31-1X BRF，分别通过视空间角度积分运算获得直 射反照率； 利用观测时刻获取的气溶胶和水汽等大气关键成分的光学特性参数，运行福射传输模 型，获得太阳入射福射场中的漫总比参数，并构建天空光角度分布函数； 通过将直射反照率在入射空间进行角度积分，并利用已获得的天空光角度分布函数作 为权重，在所有积分角度对直射反照率进行积分计算，从而获得实际天空光分布影响下的 漫射反照率； 所构建的顾及实际天空光分布影响的地表反照率遥感模型如公式(1):模型的公式（1)中，L表示福射亮度，E表示福射照度或福射通量密度，符号?/U表示福射 传输的方向，分别对应向上/向下，下标ms表示多次散射，0S和Φβ分别为太阳天顶角和方位 角，目V和Φ V则分别表示观测天顶角和方位角，f dt (λ，τ )表示下行福射漫总比，Qdir (0S，λ)和 口 ani_dif(A，T)分别表示直射和漫射反照率，G=(l-<W_difaatmrl为地-气系统相互作用增益 因子，Gd和Gu分别表示地-气交界处向下和向上的增益因子，λ为传感器观测通道中屯、波长， 公式（1)中，aani_difa，T)用公式(2)计算：公式(2)中：为天空光角度分布函数； 模型的公式(1)在理论推导中运用了下面的关系式：公式(4)中，Qs=(0s，Φβ)与Ων=(θν，Φν)分别表示太阳与观测位置，P为地表高度处 的大气压强，Rrayl与Rrayl + aer。分别表示瑞利大气、瑞利分子和气溶胶混合大气的固有发射 率，为水汽柱含量，μ表示角度余弦值，至此便实现了高精度地表谱反照率的理论计算 和遥感模型的构建。2. 根据权利要求1所述的顾及实际天空光分布影响的地表反照率遥感模型，其特征在 于:考虑了实际天空光的分布特点。3. 根据权利要求1所述的顾及实际天空光分布影响的地表反照率遥感模型，其特征在 于:考虑了实际大气与地表的多次相互作用对地表水平处向下的福射照度场和向上的反射 福射场的影响机制。4. 根据权利要求1所述的顾及实际天空光分布影响的地表反照率遥感模型，其特征在 于:可应用于非均一分布的下垫面地表。5. 根据权利要求1所述的顾及实际天空光分布影响的地表反照率遥感模型，其特征在 于:包含了太阳入射角度和大气浑浊度运两个最关键的影响因素，所述模型对大角度入照 条件和较浑浊大气条件下地表反照率的高精度遥感具有技术优势。6. 根据权利要求1所述的顾及实际天空光分布影响的地表反照率遥感模型，其特征在 于:有望解决高缔度地区尤其是冬春季节地表反射特性精度遥感不足的问题。7. 根据权利要求1所述的顾及实际天空光分布影响的地表反照率遥感模型，其特征在 于:为任何设及地表-大气禪合福射传输的高精度计算和特征目标的福射修正提供唯一下 垫面地表反射特性接口参数支持。8. 根据权利要求1所述的顾及实际天空光分布影响的地表反照率遥感模型，其特征在 于:为提高卫星遥感反照率产品与地基测量数据的一致性、减小数据交叉比对误差提供了 桥梁衔接的作用，有效支持空间高分辨卫星遥感数据的地面真实性检验和临近空间光电系 统的实时标校。
【文档编号】G06F19/00GK105975777SQ201610297746
【公开日】2016年9月28日
【申请日】2016年5月4日
【发明人】崔生成, 朱文越, 李学彬, 张梓晗, 杨世植, 饶瑞中
【申请人】中国科学院合肥物质科学研究院