一种基于低分辨率卫星遥感数据的环境变量动态筛选建模降尺度方法技术

本发明专利技术公开了一种基于低分辨率卫星遥感数据的环境变量动态筛选建模降尺度方法。本发明专利技术首先把1km的环境变量因子包括植被指数、数字高程模型、白天地表温、晚上地表温、地形湿度指数、坡度、坡向、坡长坡度8个数据进行聚合计算到25km，作为自变量，对应的25km分辨率的TRMM 3B43 v7降水数据作为因变量。M5的方法根据地理相似性将各个环境变量组成的数据集划分到不同的向量空间，然后在不同的向量空间里分别动态的筛选最有效的环境变量，并在相应的向量空间分别建立分区的多元回归模型；该模型最终应用到1km的环境变量上，最终得到1km分辨率的降水产品。通过分区与动态因子筛选的降尺度结果要明显优于基于常规回归模型的降尺度结果。

【技术实现步骤摘要】

本专利技术涉及一种基于低分辨率卫星遥感数据的环境变量动态筛选建模降尺度方法，具体涉及到一种TRMM 3B43 v7降水数据的分区环境因子动态筛选建模的降尺度方法。技术背景降水在水文学、气象学、生态学以及农业研究等领域担任了重要角色，特别是物质能量交换守恒的一个重要组成部分。地面观测站是一种应用广泛的降水测量手段，并且具有精度高和技术成熟的特点。但是地面观测站监测的降水量仅代表地表观测站及周边一定距离的降水状况，因此很难表述大面积降水分布特征，尤其是在地面观测站布网密度稀疏的高原地区。而卫星遥感技术能够提供较高时空分辨率的降水数据，覆盖空间范围更广，很好的克服了地面降水观测站和测雨雷达的局限，为全球降水监测提供了有力的数据支撑。近年来，随着气象卫星技术的发展，全球尺度高时空分辨率的测雨卫星产品应运而生，如美国热带降水测量卫星(Tropical Rainfall Measuring Mission)降水产品TRMM3B43 v7。TRMM降水卫星提供覆盖全球40°S～40°N以内的区域的降水数据。但是，TRMM卫星的原始分辨率较低(空间分辨率为0.25°，约25km)，在预测区域尺度降水方面具有一定的局限性和偏差，因此需要针对TRMM数据进行空间上的尺度转换，从而得到分辨率较高的降水测量值。但是目前尚没有一种方法能够比较精确地对复杂地区进行降水预测。
技术实现思路
本专利技术的目的在于解决现有技术中存在的问题，并提供一种基于低分辨率卫星遥感数据的环境变量动态筛选建模降尺度方法。本专利技术中的动态降尺度方法结合了M5方法与多元回归建模，根据环境变量因子与卫星降水数据组成的向量空间，采用M5方法进行分区和动态最优环境变量因子筛选，并建立分区的多元回归模型。本专利技术的具体技术方案如下：一种基于低分辨率卫星遥感数据的环境变量动态筛选建模降尺度方法，其特征在于，包括以下步骤：步骤1)数据获取：获取待测区域的TRMM 3B43 v7降水数据、MODIS卫星遥感影像数据以及ASTERGDEM卫星遥感影像数据，同时收集该待测区域内地面观测站点的日降水量观测值；其中MODIS卫星遥感影像数据包括MOD11A2数据产品和MOD13A2数据产品；步骤2)数据预处理：将步骤1)获取的TRMM 3B43 v7降水数据的时间分辨率处理为月；将ASTER GDEM卫星遥感影像数据进行聚合计算分别得到空间分辨率为1km和25km的DEM数据；从MOD11A2数据产品中提取白天地表温度和晚上地表温度参量，并通过聚合计算分别得到空间分辨率为1km和25km的白天地表温度数据以及空间分辨率为1km和25km的晚上地表温度数据；从MOD13A2数据产品中提取植被指数参量，经过异常值剔除处理后，通过聚合计算分别得到空间分辨率为1km和25km的植被指数数据；从ASTER GDEM卫星遥感影像数据中提取坡度、地形湿度指数、坡长坡度和坡向4个参量进行聚合计算分别得到1km和25km的坡度数据、地形湿度指数数据、坡长坡度数据和坡向数据；步骤3)进行分区及环境因子动态筛选回归建模：将步骤2)处理后的25kmTRMM 3B43 v7降水数据作为因变量，以空间分辨率为25km的白天地表温度数据、晚上地表温度数据、植被指数数据、DEM数据、坡度数据、坡向数据、坡长坡度数据和地形湿度指数数据作为自变量；将25km各个环境变量因子组成的数据集划分到不同的向量空间并通过计算不同环境变量因子的影响权重动态筛选出最有效的环境变量因子，然后对每个向量空间分别建立不同的多元回归模型。步骤4)降尺度预测：基于步骤3)确定1km尺度下相应的环境变量种类，然后通过25km环境变量因子组成的不同向量空间建立的多元回归模型计算得到空间分辨率为1km的降尺度降水产品；同时将空间分辨率为25km的降水回归残差值进行重采样得到空间分辨率为1km的降水回归残差值，并将其与空间分辨率为1km地面降水量预测值数据相加，得到空间分辨率为1km的TRMM 3B43 v7降水数据。所述的步骤1)中，TRMM 3B43 v7降水数据的空间分辨率为0.25°×0.25°，时间分辨率为3小时；所述的ASTER GDEM卫星遥感影像数据的空间分辨率为90m；所述的MODIS卫星遥感影像数据的空间分辨率为1km，时间分辨率为8天。所述的步骤2)中异常值剔除处理的具体步骤如下：将MOD13A2数据产品中提取的植被指数作为初始植被指数，首先删除初始植被指数中栅格值小于0的部分，再以10×10的窗口移动平滑植被指数，然后用初始植被指数减去平滑后的植被指数，再选择-0.1到0.1作为阈值范围对相减后的结果进行筛选，舍去超出阈值范围的栅格，同时剔除坡度为-1所对应的植被指数的范围，其余作为正常的植被指数点。进行分区及环境因子动态筛选回归建模采用M5决策树分类方法，具体为：根据地理相似性将各个环境变量组成的数据集划分到不同的向量空间，然后在不同的向量空间里分别动态筛选最有效的环境变量，并在相应的向量空间分别建立多元回归关系。所述的多元回归模型通用形式为：其中，N表示参数估算模型中自变量个数，具体视上述因子选择情况而定；an表示第n个环境变量的系数；a0表示模型参数的常数项系数；yn表示降雨量预测值；xn表示第n个环境变量；a0和an的计算公式如下： a n = Σ i = 1 k Σ n = 1 N ( x i n - x n ‾ ) ( y i - y ‾ ) Σ i = 1 k Σ n = 本文档来自技高网...

【技术保护点】
一种基于低分辨率卫星遥感数据的环境变量动态筛选建模降尺度方法，其特征在于，包括以下步骤：步骤1)数据获取：获取待测区域的TRMM 3B43 v7降水数据、MODIS卫星遥感影像数据以及ASTERGDEM卫星遥感影像数据，同时收集该待测区域内地面观测站点的日降水量观测值；其中MODIS卫星遥感影像数据包括MOD11A2数据产品和MOD13A2数据产品；步骤2)数据预处理：将步骤1)获取的TRMM 3B43 v7降水数据的时间分辨率处理为月；将ASTER GDEM卫星遥感影像数据进行聚合计算分别得到空间分辨率为1km和25km的DEM数据；从MOD11A2数据产品中提取白天地表温度和晚上地表温度参量，并通过聚合计算分别得到空间分辨率为1km和25km的白天地表温度数据以及空间分辨率为1km和25km的晚上地表温度数据；从MOD13A2数据产品中提取植被指数参量，经过异常值剔除处理后，通过聚合计算分别得到空间分辨率为1km和25km的植被指数数据；从ASTER GDEM卫星遥感影像数据中提取坡度、地形湿度指数、坡长坡度和坡向4个参量进行聚合计算分别得到1km和25km的坡度数据、地形湿度指数数据、坡长坡度数据和坡向数据；步骤3)进行分区及环境因子动态筛选回归建模：将步骤2)处理后的25kmTRMM 3B43 v7降水数据作为因变量，以空间分辨率为25km的白天地表温度数据、晚上地表温度数据、植被指数数据、DEM数据、坡度数据、坡向数据、坡长坡度数据和地形湿度指数数据作为自变量；将25km各个环境变量因子组成的数据集划分到不同的向量空间并通过计算不同环境变量因子的影响权重动态筛选出最有效的环境变量因子，然后对每个向量空间分别建立不同的多元回归模型；步骤4)降尺度预测：基于步骤3)确定1km尺度下相应的环境变量种类，然后通过25km环境变量因子组成的不同向量空间建立的多元回归模型计算得到空间分辨率为1km的降尺度降水产品；同时将空间分辨率为25km的降水回归残差值进行重采样得到空间分辨率为1km的降水回归残差值，并将其与空间分辨率为1km地面降水量预测值数据相加，得到空间分辨率为1km的TRMM 3B43 v7降水数据。...

【技术特征摘要】
1.一种基于低分辨率卫星遥感数据的环境变量动态筛选建模降尺度方法，其特征在于，包括以下步骤：步骤1)数据获取：获取待测区域的TRMM 3B43 v7降水数据、MODIS卫星遥感影像数据以及ASTERGDEM卫星遥感影像数据，同时收集该待测区域内地面观测站点的日降水量观测值；其中MODIS卫星遥感影像数据包括MOD11A2数据产品和MOD13A2数据产品；步骤2)数据预处理：将步骤1)获取的TRMM 3B43 v7降水数据的时间分辨率处理为月；将ASTER GDEM卫星遥感影像数据进行聚合计算分别得到空间分辨率为1km和25km的DEM数据；从MOD11A2数据产品中提取白天地表温度和晚上地表温度参量，并通过聚合计算分别得到空间分辨率为1km和25km的白天地表温度数据以及空间分辨率为1km和25km的晚上地表温度数据；从MOD13A2数据产品中提取植被指数参量，经过异常值剔除处理后，通过聚合计算分别得到空间分辨率为1km和25km的植被指数数据；从ASTER GDEM卫星遥感影像数据中提取坡度、地形湿度指数、坡长坡度和坡向4个参量进行聚合计算分别得到1km和25km的坡度数据、地形湿度指数数据、坡长坡度数据和坡向数据；步骤3)进行分区及环境因子动态筛选回归建模：将步骤2)处理后的25kmTRMM 3B43 v7降水数据作为因变量，以空间分辨率为25km的白天地表温度数据、晚上地表温度数据、植被指数数据、DEM数据、坡度数据、坡向数据、坡长坡度数据和地形湿度指数数据作为自变量；将25km各个环境变量因子组成的数据集划分到不同的向量空间并通过计算不同环境变量因子的影响权重动态筛选出最有效的环境变量因子，然后对每个向量空间分别建立不同的多元回归模型；步骤4)降尺度预测：基于步骤3)确定1km尺度下相应的环境变量种类，然后通过25km环境变量因子组成的不同向量空间建立的多元回归模型计算得到空间分辨率为1km的降尺度降水产品；同时将空间分辨率为25km的降水回归残差值进行重采样得到空间分辨率为1km的降水回归残差值，并将其与空间分辨率为1km地面降水量预测值数据相加，得到空间分辨率为1km的TRMM 3B43 v7降水数据。2.如权利要求1所述的一种基于低分辨率卫星遥感数据的环境变量动态筛选建模降尺度方法，其特征在于，所述的步骤1)中，TRMM 3B43 v7降水数据的空间分辨率为0.25°×0.25°，时间分辨率为3小时；所述的ASTER GDEM卫星遥感影像数据的空间分辨率为90m；所述的MODIS卫星遥感影像数据的空间分辨率为1km，时间分辨率为8天。3.如权利要求1所述的一种基于低分辨率卫星遥感数据的环境变量动态筛选建模降尺度方法，其特征在于，所述的步骤2)中异常值剔除处理的具体步骤如下：将MOD13A2数据产品中提取的植被指数作为初始植被指数，首先删除初始植被指数中栅格值小于0的部分，再以10×10的窗口移动平滑植被指数，然后用初始植被指数减去平滑后的植被指数，再选择-0.1到0.1作为阈值范围对相减后的结果进行筛选，舍去超出阈值范围的栅格，同时剔除坡度为-1所对应的植被指数的范围，其余作为正常的植被指数点。4.如权利要求1所述的一种基于低分辨率卫星遥感数据的环境变量动态筛选建模降尺度方法，其特征在于，进行分区及环境因子动态筛选回归建模采用M5决策树分类方法，具体为：根据地理相似性将各个环境变量组成的数据集划分到不同的向量空间，然后在不同的向量空间里分别动态筛选最有效的环境变量，并在相应的向量空间分别建立多元回归关系。5.如权利要求1所述的一种基于低分辨率卫星遥感数据的环境变量动态筛选建模降尺度方法，其特征在于，所述的多元回归模型形式为：其中，N表示参数估算模型中自变量个数；an表示第n个环境变量的系数；a0表示模型参数的常数项系数；yn表示降雨量预测值；xn表示第n个环境变量；a0和an的计算公式如下： a n = Σ i = 1 k Σ n = 1 N ( x ...

【专利技术属性】
技术研发人员：史舟，马自强，吕志强，刘用，
申请(专利权)人：浙江大学，
类型：发明
国别省市：浙江;33