[发明专利]地表反照率反演方法及系统

说明书

技术领域

本发明涉及遥感技术领域，特别涉及一种地表反照率反演方法及系统。

背景技术

在众多地表参数中，地表双向反射分布函数(BidirectionalReflectance Distribution Function，BRDF)/反照率的准确遥感反演从狭义上讲常常是准确估算其他陆地表面参数如植被参数、土地利用和土地覆盖的先决条件。因此，作为影响地球气候系统的关键变量，地表反照率是数值气候模型和地表能量平衡方程中的一个重要参数。准确计算下垫面反照率可以揭示局地和区域气候变化的内在机制，可以更好的对生态系统模型进行参数化，提高中、长期气候预报的水平。虽然地表BRDF/反照率遥感反演具有非常重要的意义，但是在BRDF/反照率反演过程中，面临主要的问题表现在：

1)传感器获得的多角度的地表二向反射率数据利用效率低，如何协同利用多源遥感数据有待加强

BRDF反演需要多个角度的信息，且反照率作为所有短波反射能量与入射能量之比，包括地表半球空间所有方向的反射能量。多数情况下由于方向反射的原因，不同方向地表反射率不同，不同地表组成，在各方向的反射率也有差异，故精确进行地表BRDF/反照率遥感反演需要所有半球方向的多角度的地表二向反射率数据。而这对于遥感技术是较难实现的。目前大多数传感器只在一个方向进行观测。即使是对于大视场传感器，在假定地表短期变化不大的情况下，利用同一地表目标的多时相的重复观测可以获取多个角度的信息。但是所需周期较长，如中分辨率成像光谱仪(Moderate-Resolution ImagingSpectroradiometer，MODIS)需要半个月的时间生成BRDF/反照率的产品。而且由于时间跨度大，很难保证地表的不变性。虽然如此，即使将所有角度都用上也无法满足直接进行反照率提取的需求。多源数据融合有待加强，充分利用多光谱信息，将多角度遥感信息与多光谱遥感信息相结合。

2)波谱信息利用形式单一，与多角度数据结合有待加强

地表反照率还是全波段光谱反照率之和，但是目前传感器观测地表的波段数据有限，传感器无法覆盖全波段，因此用有限的波段去计算全波段的反照率具有比较大的误差。在全波段反照率获取时，最早的扫描辐射仪AVHRR。只能提供两个可用波段，MSS四个波段，目前专题绘图仪(Thematic Mapper，TM)可以提高到6个波段，即使是MODIS也只能提供20个波段。目前的BRDF/反照率(ALBEDO)模型主要是针对传感器设置波段，反演时也是分别实施，而没有充分利用多光谱信息的优势，而不同传感器之间波段设置区域与响应的差距也给在多源遥感数据协同反演时带来困难，使得来自不同平台传感器的数据不能有效地结合，如何有效的综合利用多角度多光谱的遥感数据还有待于加强，因为两者结合可为我们提供更多的关于植被本身的信息，也可以在反演中降低对多角度的地表二向反射率数据量的要求。

3)地表的二向反射获取

目前虽有比较多的模型来描述地表的二向反射特性。但由于观测方向的有限性，目前很难精确进行地表二向反射的获取。目前物理模型有些过于复杂，参数反演较为困难。作为半经验的核驱动模型也需要多角度的地表二向反射率数据，并且进行迭代运算获取模型的参数，运算量大而且不太稳定。

总而言之，地表BRDF/反照率遥感反演需要多角度信息的支持，目前只能通过多角度/大视场传感器如MODIS、POLDER、MISR等实现。单一传感器利用多角度信息进行的BRDF/反照率遥感反演往往具有信息量不足，反演精度不能满足要求，特别是在多角度的地表二向反射率数据角度分布不合理条件下，反演存在极大的风险性，有时会生成极端错误的结果；同时由于传感器波段设置的局限，单一的传感器往往不能覆盖全波段，用极有限的波段求取全波段的地表反照率也具有比较大的误差。在获取宽波段反照率时需要将反演得到的窄波段反照率转化到宽波段上，目前针对这个问题多数研究采用了给定经验权重的算法，这一过程无疑又会带来额外的误差，而且权重系数也随不同地表类型及不同卫星而变化，不同学者给出的结果也存在较大的差异，地域性很强，难以推广，同时这种差异性也制约着综合利用多源遥感数据协同反演，不能更好的综合利用有限的数据资源。

发明内容

(一)要解决的技术问题

本发明要解决的技术问题是如何解决观测信息量不足，传统算法中窄波段向宽波段转化误差很大，以及遥感数据利用率低的问题。

(二)技术方案

为解决上述技术问题，本发明提供了一种地表反照率反演方法，包括以下步骤：

S1：通过观测卫星的遥感数据上的每个像元获取多角度的地表二向反射率数据；