一种基于光照强度的遥感反射率校正方法与流程

1.本发明涉及一种遥感反射率校正方法，具体涉及一种基于光照强度的遥感反射率校正方法，属于遥感数据处理计算领域。


背景技术：

2.地面光谱数据的采集作为遥感定量反演的重要部分，直接影响遥感反演结果的准确性。光照强度的变化对地面光谱数据采集质量有明显影响，在进行固定测点、固定时间、长时间地面光谱数据采集过程中光照强度变化严重影响光谱数据的准确性。目前针对该问题，主要的措施是在光照条件充足的情况下，开展数据测量，保证数据质量。但是这些方法无法有效地处理光照条件不足情况下的光谱数据，也无法有效地保障数据的长时间观测。


技术实现要素：

3.本发明的目的在于克服上述不足而提供一种基于光照强度的遥感反射率校正方法。该方法主要校正光照强度变化对固定站点的光谱实测数据的影响，消除光照强度变化对测量结果的干扰，提高获取数据的有效率。
4.本发明采取的技术方案为：
5.一种基于光照强度的遥感反射率校正方法，包括步骤如下：
6.(1)观测数据准备
7.准备长时间、固定站点观测的光照强度数据以及遥感反射率数据，并记录观测时间；根据测量时间对数据进行整理，每年相同日期相同时间点的数据视为同一组；
8.(2)确定基准数据
9.光照强度主要受所在地的经纬度、观测时间影响，可以根据所在地的经纬度、测量时间计算出理论上该观测点该时间的光照强度，但是受其它因素干扰，很难达到理论值；在历史数据库中选择实测光照强度最接近理论光照强度的数据，作为基准光照强度数据(t_l
x_j_q
)，最接近理论光照强度下测得的遥感反射率值作为基准遥感反射率(t_r
rs_j
)；
10.(3)将光照强度进行归一化处理
11.由于光照强度在不同情况下变化幅度非常大，不利于构建光照强度与原始遥感反射率关系构建，因此首先对光照强度数据进行归一化处理，处理方式如式(一)、(二)所示：
[0012][0013][0014]
式中t_l
x
为t时间归一化计算之后的光照强度，t_l
x_0
为t时间准备数据中观测的光照强度，t_l
x_j_q
为基准光照强度，t_l
x_j
为基准归一化光照强度，t_l
x_min
为t时间准备数据中观测的光照强度最小值，t_l
x_max
为t时间准备数据中观测的光照强度最大值；
[0015]
(4)建立原始遥感反射率与归一化光照强度的相关关系
[0016]
对相同测量时间的归一化光照强度(t_l
x
)与原始遥感反射率(t_r
rs
)进行主成分分析，分析结果表明(t_l
x
)2、与原始遥感反射率(t_r
rs
)的相关性最高，故选择(t_l
x
)2、与t_r
rs
进行多元相关性分析，利用相关性分析工具，经过分析其相关关系如式(三)所示：
[0017][0018]
式中t_r
rs
为t时间原始遥感反射率，β、γ为校正参数，t_l
x
为t时间归一化光照强度，b为常数项，其中相同观测时间的校正参数β、γ相同；
[0019]
(5)计算基准遥感反射率与基准归一化光照强度的相关关系
[0020]
通过式(三)，计算基准遥感反射率(t_r
rs_j
)与基准归一化光照强度(t_l
x_j
)的相关关系，结果如式(四)所示：
[0021][0022]
式中：β
′
、γ
′
为校正参数，b
′
为常数项；
[0023]
(6)建立校正模型
[0024]
通过式(三)与式(四)计算t_r
rs_j
与t_r
rs
的比值，得到的结果如式(五)所示：
[0025][0026]
通过式(六)的变化，令t_k为校正参数，
[0027][0028]
故校正模型如式(七)所示：
[0029][0030]
式中t_r
rs_h
为遥感反射率校正后的结果，t_r
rs
为实测遥感反射率；
[0031]
(7)数据校正
[0032]
根据待校正数据中的测量时间，选定式子中对应的时间t，获取传入数据光照强度对应的归一化光照强度值和基准归一化光照强度，如果传入的光照强度大于历史数据的最大值或小于历史数据的最小值时，视为异常数据，不进行校正；若数据正常时，利用归一化光照强度值和基准归一化光照强度由式(六)计算校正参数，再利用校正参数和传入数据中的遥感反射率数据由式(七)计算得到校正后的遥感反射率数据。
[0033]
本发明通过主成分分析，选择归一化光照强度中与遥感反射率相关关系最强的三组主成分，构建遥感反射率与归一化光照强度之间的关系，实现弱光照强度下获取的遥感反射率校正为正常光照强度下的遥感反射率。本发明能消除光照强度对观测数据的干扰，能够将光照条件不满足使用要求的数据通过校正可以正常使用；实现校正后的数据保障了点位观测数据的连续性；确保了建模的精度和反演结果的准确性。
附图说明
[0034]
图1为本发明方法的流程图。
具体实施方式
[0035]
为了更好地阐述本发明技术方案，下面结合具体实施例进一步说明。
[0036]
实施例1
[0037]
一种基于光照强度的遥感反射率校正方法，包括步骤如下：
[0038]
1.准备长时间、固定站点观测的光照强度数据以及遥感反射率数据，并记录观测时间；根据测量时间对数据进行整理。
[0039]
2.1确定基准数据
[0040]
本次需要校正的遥感数据是x月x日t点的某地(xx
°
xx
′
n、xxx
°
xx
′
e)获取的实测光照强度t_l
x_0
下的遥感反射率数据t_r
rs
。
[0041]
首先根据经纬度和时间信息确定该地该时间理论的光照强度为l_0，通过历史数据库中发现同期数据中最接近的光照强度为此时测得的遥感反射率为t_r
rs_j
，由于受其它因素干扰，t_l
x_j_q
略小于l_0，通过长期的数据观测，t_l
x_j_q
可以作为该地、该时间的基准光照强度，t_l
x_j_q
对应下测得的遥感反射率t_r
rs_j
作为基准遥感反射率。
[0042]
2.2将光照强度进行归一化处理
[0043]
由于光照强度在不同情况下变化幅度非常大，不利于构建光照强度与原始遥感反射率关系构建，因此首先对光照强度数据进行归一化处理，综合历史数据中选择t时刻的实测光照强度最大值和最小值，分别作为t_l
x_max
、t_l
x_min
，通过判断t_l
x_min
＜t_l
x_0
＜t_l
x_max
数据在正常范围内，可以进行后续计算。
[0044]
对t_l
x_0
数据进行归一化处理，归一化后的光照强度基准归一化后光照强度
[0045]
2.3建立原始遥感反射率与归一化光照强度相关关系
[0046]
对相同测量时间的归一化光照强度(t_l
x
)与原始遥感反射率(t_r
rs
)进行主成分分析，分析结果表明(t_l
x
)2、与原始遥感反射率(t_r
rs
)的相关性最高，相关关系式为：
[0047][0048]
β、γ为校正参数，b为常数项。
[0049]
2.4计算基准遥感反射率与基准归一化光照强度相关关系
[0050]
类似步骤2.3，计算基准遥感反射率(t_r
rs_j
)与基准归一化光照强度(t_l
x_j
)的相关关系，相关关系式为：
[0051][0052]
β
′
、γ
′
为校正参数，b
′
为常数项。
[0053]
2.5建立校正模型
[0054]
通过步骤2.3与2.4结果计算t_r
rs_j
与t_r
rs
的比值，得到的结果：
[0055][0056]
校正参数t_k的值为：
[0057][0058]
故校正后的数据为：
[0059]
本发明的以上描述用于使本领域的任何普通技术人员能够实现或使用本发明。本发明并不限于本文所述的实例和设计，而是与本文公开的原理和新颖性特性的最广范围相一致。