[发明专利]一种分析蒸散发影响因子及建立多元回归模型的方法

权利要求书

1.一种分析蒸散发影响因子及建立多元回归模型的方法，包括如下步骤：

步骤1，确定研究流域，收集所述研究流域内的气象资料、遥感资料、流域下垫面资料、人口资料和植被覆盖率资料并对收集的资料进行预处理，所述气象资料、所述流域下垫面资料、所述人口资料和所述植被覆盖率资料均为逐日数据，遥感数据资料为瞬时数据；

步骤2，将步骤1预处理后的所述遥感资料按照研究流域进行裁剪，所述遥感资料用于所述研究流域地表温度的反演，基于SEBAL模型估算所述研究流域瞬时蒸散发量；将瞬时蒸散发量估算结果进行日尺度时间扩展，得到所述研究流域的日蒸散发量，并收集所述研究流域蒸散发的影响因子；

步骤3，利用所述遥感资料中的近红外和红外数据计算NDVI，基于NDVI与地表类型之间的关系，对所述研究流域进行分区；

步骤4，获取所述分区内所述蒸散发的影响因子，采用敏感性系数公式计算蒸散发对所述影响因子的敏感系数，采用相关系数公式计算蒸散发对所述影响因子的相关性，并对其相关性进行t检验；根据所述影响因子的所述敏感性系数、相关性系数和所述t检验结果筛选主要影响因子，借助R语言的lm函数对所述分区内蒸散发量及其所述主要影响因子建立多元回归模型，所述多元回归模型包括所述主要影响 因子的交互关系；依据AIC信息准则进行逐步筛选，结合所述AIC值和显著性检验结果建立最优多元回归模型。

2.根据权利要求1所述一种分析蒸散发影响因子及建立多元回归模型的方法，其特征在于：所述遥感资料的预处理主要为利用ENVI进行辐射定标、几何校正和大气校正，所述遥感资料为Landsat TM遥感资料；其他资料的预处理方法主要为剔除异常数据、统一数据尺度。

3.根据权利要求1所述一种分析蒸散发影响因子及建立多元回归模型的方法，其特征在于：所述影响因子包括植被覆盖率、气候干旱指数、二氧化碳含量、气温、风速风向、人口状况、有效灌溉面积、作物种值面积、土壤特性。

4.根据权利要求1所述一种分析蒸散发影响因子及建立多元回归模型的方法，其特征在于：利用所述遥感资料进行所述流域地表温度反演的方法为覃志豪单窗算法，其计算公式如下：

T_s＝[a6(1-C6)+[(b6-1)(1-C6-D6)+1]T6-D6T_a]/C6

其中：T_s为地表温度，K；T_a为大气平均温度，K；T6为亮度温度，K；a6和b6为常量，一般情况下地表温度在0～50℃范围内，取值a6＝-67.9542，b6＝0.45987；C6和D6是中间变量，分别为：

C6＝ε₆τ₆

D6＝(1-ε₆)[1+(1-ε₆)τ₆]

其中：ε₆为地表比辐射率，τ₆为大气透射率。

5.根据权利要求1所述一种分析蒸散发影响因子及建立多元回归模型的方法，其特征在于：所述SEBAL模型估算所述研究流域瞬时蒸散发量，是通过地表能量平衡相关分量，获得潜热通量，进而得到所述研究流域的瞬时蒸散发量；

所述地表能量平衡方程公式如下：

LE＝R_n-G+H

其中：R_n是净辐射通量，W/m²；G是土壤热通量，W/m²；H是感热通量，W/m²；λET为潜热通量，W/m²。

6.根据权利要求1所述一种分析蒸散发影响因子及建立多元回归模型的方法，其特征在于：所述瞬时蒸散发量估算结果进行日尺度时间扩展方法为恒定EF法，所述恒定EF法的公式如下：

其中：ET_d为日蒸散发量，mm/d；f为时间转换因子；λ_v为汽化潜热，J/kg；ρ_w为水的密度，kg/m³。