空间降尺度降水数据检测方法、装置及电子设备与流程

本发明涉及气象信息技术领域，特别是涉及一种空间降尺度降水数据检测方法、装置及电子设备。

背景技术：

降水强烈的时空变化常使得常规地基台站的降水观测不能准确把握降水的空间分布和强度变化，而卫星微波遥感技术可以克服此局限，实现在全球范围内实现对降水的观测，而且相对于红外/可见光只能反映云厚、云高等信息而言，卫星微波能够穿透云体，利用云内降水粒子和云粒子与微波的相互作用对云、雨进行更为直接的探测，实现对降水数据的监测，然而，由于微波遥感技术发展历史较短，加之降水发生的高随机性与异质性，使得已有微波遥感降水数据分辨率较低，难以反映降水的空间细节特征。

技术实现要素：

基于此，本发明的目的在于，提供一种空间降尺度降水数据检测方法，其能够获取分辨率较高的降水数据。

本发明是通过如下方案实现的：

一种空间降尺度降水数据检测方法，包括如下步骤：

获取样本区域的白天地表温度和夜间地表温度，并根据所述白天地表温度和夜间地表温度获取样本区域的昼夜地表温度差；

获取样本区域的数字高程模型和样本区域的植被指数；

将第一分辨率的样本区域的白天地表温度、夜间地表温度、昼夜地表温度差、数字高程模型和植被指数重采样为与所述样本区域的卫星遥感降水数据相同的第二分辨率，其中，第二分辨率小于第一分辨率；

根据样本区域第二分辨率的白天地表温度、夜间地表温度、昼夜地表温度差、数字高程模型和植被指数，以及所述样本区域的卫星遥感降水数据，建立并训练最优随机森林算法模型，其中，所述样本区域的白天地表温度、夜间地表温度、昼夜地表温度差、数字高程模型和植被指数作为所述训练最优随机森林算法模型的输入样本，所述样本区域的卫星遥感降水数据作为所述训练最优随机森林算法模型的输出样本；

获取待检测区域的白天地表温度和夜间地表温度，并根据所述待检测区域的白天地表温度和夜间地表温度获取待检测区域的昼夜地表温度差；

获取待检测区域的数字高程模型和待检测区域的植被指数；

将所述第一分辨率的待检测区域白天地表温度、夜间地表温度、昼夜地表温度差、数字高程模型和植被指数输入最优随机森林算法模型，获取待检测区域第一分辨率的降水数据；

对所述第一分辨率的降水数据进行残差校正，获取待检测区域的空间降尺度降水数据。

本发明所述的空间降尺度降水数据检测方法，利用有卫星遥感降水数据的样本区域的白天地表温度、夜间地表温度、昼夜地表温度差、数字高程模型和植被指数以及样本区域的遥感降水数据，建立最优的随机森林算法模型，利用该模型，计算出待检测区域的高分辨率的降水数据，并对该降水数据进行残差校正，得到更为准确的空间降尺度的降水数据。

在一种实施例中，对所述第一分辨率的降水数据进行残差校正，获取待检测区域的空间降尺度降水数据，包括：

将待检测区域第一分辨率的降水数据重采样至第二分辨率；

计算待检测区域卫星遥感降水数据与重采样后的待检测区域第二分辨率的降水数据之间的差值，获取第一残差；

将所述第一残差空间内插至第一分辨率，获取第一分辨率的第二残差；

将待检测区域第一分辨率的降水数据与第二残差相加，获取待检测区域的空间降尺度降水数据。

在一种实施例中，所述样本区域与所述待检测区域的植被指数大于零。

在一种实施例中，所述第一分辨率为1km*1km，所述第二分辨率为25km*25km。

在一种实施例中，将第一分辨率的样本区域的白天地表温度、夜间地表温度、昼夜地表温度差、数字高程模型和植被指数重采样为与所述样本区域的卫星遥感降水数据相同的第二分辨率，包括：

计算第二分辨率像元范围内，所有第一分辨率像元的平均值。

进一步地，本发明还提供一种空间降尺度降水数据检测装置，包括：

第一数据采集模块，用于获取样本区域的白天地表温度和夜间地表温度，并根据所述白天地表温度和夜间地表温度获取样本区域的昼夜地表温度差，以及获取样本区域的数字高程模型和样本区域的植被指数；

第一重采样模块，用于将第一分辨率的样本区域的白天地表温度、夜间地表温度、昼夜地表温度差、数字高程模型和植被指数重采样为与所述样本区域的卫星遥感降水数据相同的第二分辨率；

随机森林训练模块，用于根据样本区域第二分辨率的白天地表温度、夜间地表温度、昼夜地表温度差、数字高程模型和植被指数，以及所述样本区域的卫星遥感降水数据，建立并训练最优随机森林算法模型，其中，所述样本区域的白天地表温度、夜间地表温度、昼夜地表温度差、数字高程模型和植被指数作为所述训练最优随机森林算法模型的输入样本，所述样本区域的卫星遥感降水数据作为所述训练最优随机森林算法模型的输出样本；

第二数据采集模块，用于获取待检测区域的白天地表温度和夜间地表温度，并根据所述待检测区域的白天地表温度和夜间地表温度获取待检测区域的昼夜地表温度差，以及获取待检测区域的数字高程模型和待检测区域的植被指数；

第一降水数据获取模块，用于将所述第一分辨率的待检测区域白天地表温度、夜间地表温度、昼夜地表温度差、数字高程模型和植被指数输入最优随机森林算法模型，获取待检测区域第一分辨率的降水数据；

残差校正模块，用于对所述第一分辨率的降水数据进行残差校正，获取待检测区域的空间降尺度降水数据。

本发明所述的空间降尺度降水数据检测装置，利用有卫星遥感降水数据的样本区域的白天地表温度、夜间地表温度、昼夜地表温度差、数字高程模型和植被指数以及样本区域的遥感降水数据，建立最优的随机森林算法模型，利用该模型，计算出待检测区域的高分辨率的降水数据，并对该降水数据进行残差校正，得到更为准确的空间降尺度的降水数据。

在一种实施例中，所述残差校正模块，包括：

第二重采样单元，用于将待检测区域第一分辨率的降水数据重采样至第二分辨率；

第一残差获取单元，用于计算待检测区域卫星遥感降水数据与重采样后的待检测区域第二分辨率的降水数据之间的差值，获取第一残差；

第二残差获取单元，用于将所述第一残差空间内插至第一分辨率，获取第一分辨率的第二残差；

降尺度降水数据获取单元，用于将待检测区域第一分辨率的降水数据与第二残差相加，获取待检测区域的空间降尺度降水数据。

在一种实施例中，所述第一重采样模块包括分辨率计算单元，用于计算第二分辨率像元范围内，所有第一分辨率像元的平均值。

进一步地，本发明还提供一种计算机可读介质，其上存储有计算机程序，该计算机程序被处理器执行时实现如上述任意一项空间降尺度降水数据检测方法。

进一步地，本发明还提供一种电子设备，包括存储器、处理器以及储存在所述存储器并可被所述处理器执行的计算机程序，所述处理器执行所述计算机程序时，实现如上述的任意一项空间降尺度降水数据检测方法。

为了更好地理解和实施，下面结合附图详细说明本发明。

附图说明

图1为一种实施例中空间降尺度降水数据检测方法流程图；

图2为一种实施例中残差校正流程示意图；

图3为另一种实施例中空间降尺度降水数据检测方法流程图；

图4为一种实施例中空间降尺度降水数据检测装置结构示意图；

图5为一种实施例中电子设备结构示意图。

具体实施方式

请参阅图1，在一种实施例中，空间降尺度降水数据检测方法包括如下步骤：

步骤s10：获取样本区域的白天地表温度和夜间地表温度，并根据所述白天地表温度和夜间地表温度获取样本区域的昼夜地表温度差。

所述样本区域为有卫星遥感降水数据的区域，所述样本区域的大小为所述卫星遥感降水数据的分辨率大小的整数倍。所述白天地表温度为地表温度在整个白天内的平均值，所述夜间地表温度为地表温度在整个夜间的平均值，在一种实施例中，所述白天地表温度和所述夜间地表温度通过卫星传感器modis(中分辨率成像光谱仪moderate-resolutionimagingspectroradiometer)获取，所述样本区域的昼夜温度差为所述白天地表温度和所述夜间地表温度之间的差值。

步骤s20：获取样本区域的数字高程模型和样本区域的植被指数。

所述数字高程模型(digitalelevationmodel)，简称dem，是通过有限的地形高程数据实现对地面地形的数字化模拟，即地形表面形态的数字化表达，为一组有序数值阵列形式表示地面高程的一种实体地面模型，所述植被指数为根据植被的光谱特性，将卫星可见光和近红外波段进行组合，形成的植被覆盖指标，定性和定量评价植被覆盖及其生长活力。植被指数的值通常为-1至1，在冰雪覆盖、水体及沙漠区域，植被指数通常为小于零的常数，在本实施例中，所述样本区域为植被指数大于零的区域。

步骤s30：将第一分辨率的样本区域的白天地表温度、夜间地表温度、昼夜地表温度差、数字高程模型和植被指数重采样为与所述样本区域的卫星遥感降水数据相同的第二分辨率，其中，第二分辨率小于第一分辨率。

其中，在本实施例中，所述白天地表温度、夜间地表温度、昼夜地表温度差、数字高程模型和植被指数的第一分辨率为1km*1km，而卫星遥感降水数据的第二分辨率为25km*25km，因此，需要将第一分辨率的白天地表温度、夜间地表温度、昼夜地表温度差、数字高程模型和植被指数重采样为第二分辨率。

步骤s40：根据样本区域第二分辨率的白天地表温度、夜间地表温度、昼夜地表温度差、数字高程模型和植被指数，以及所述样本区域的卫星遥感降水数据，建立并训练最优随机森林算法模型。

所述样本区域的卫星遥感降水数据，为所述样本区域内，通过卫星遥感所监测到的降水数据，所述最优随机森林算法模型为经反复训练后，降水数据计算误差达到最小的随机森林算法模型。其中，所述样本区域的白天地表温度、夜间地表温度、昼夜地表温度差、数字高程模型和植被指数作为所述训练最优随机森林算法模型的输入样本，所述样本区域的卫星遥感降水数据作为所述训练最优随机森林算法模型的输出样本。

步骤s50：获取待检测区域的白天地表温度和夜间地表温度，并根据所述待检测区域的白天地表温度和夜间地表温度获取待检测区域的昼夜地表温度差。

在本实施例中，所述待监测区域的白天地表温度和所述夜间地表温度通过卫星传感器modis(中分辨率成像光谱仪moderate-resolutionimagingspectroradiometer)获取。

步骤s60：获取待检测区域的数字高程模型和待检测区域的植被指数。

在本实施例中，所述待检测区域为植被指数大于零的区域。

步骤s70：将所述第一分辨率的待检测区域白天地表温度、夜间地表温度、昼夜地表温度差、数字高程模型和植被指数输入最优随机森林算法模型，获取待检测区域第一分辨率的降水数据。

在本实施例中，从最优随机森林算法模型所获取的待检测区域的降水数据分辨率为1km*1km，但误差较大。

步骤s80：对所述第一分辨率的降水数据进行残差校正，获取待检测区域的空间降尺度降水数据。

所述残差为卫星遥感测定的降水数据与从最优随机森林算法模型所获取的待检测区域的降水数据之间的差值，对高分辨率的待检测区域的降水数据进行残差校正，能得到更为准确的高分辨率的待检测区域的降水数据。

本发明所述的空间降尺度降水数据检测方法，利用有卫星遥感降水数据的样本区域的白天地表温度、夜间地表温度、昼夜地表温度差、数字高程模型和植被指数以及样本区域的遥感降水数据，建立最优的随机森林算法模型，利用该模型，计算出待检测区域的高分辨率的降水数据，并对该降水数据进行残差校正，得到更为准确的空间降尺度的降水数据。

在一种实施例中，步骤s80中，对所述第一分辨率的降水数据进行残差校正，获取待检测区域的空间降尺度降水数据具体包括如下步骤：

步骤s81：将待检测区域第一分辨率的降水数据重采样至第二分辨率。

步骤s82：计算待检测区域卫星遥感降水数据与重采样后的待检测区域第二分辨率的降水数据之间的差值，获取第一残差。

步骤s83：将所述第一残差空间内插至第一分辨率，获取第一分辨率的第二残差。

步骤s84：将待检测区域第一分辨率的降水数据与第二残差相加，获取待检测区域的空间降尺度降水数据。

由于降水的高时空异质性，必然存在不能被最优随机森林算法模型有效指示的部分降水量，即残差。本发明通过将这部分残差进行降尺度并与1km分辨率降水模拟值相加的方式来最小化模型误差。由于待检测区域卫星遥感降水数据分辨率为25km*25km，因此可通过待检测区域卫星遥感降水数据值减去重采样后的待检测区域第二分辨率的降水数据计算得出25km*25km分辨率下的残差。由于残差产生的随机性，因此本实施例通过空间内插的方式对残差进行降尺度，本实施例中，使用薄板样条函数插值法(thin-platespline)对残差进行插值处理，来获取最优的降尺度结果。

在一种实施例中，将第一分辨率的白天地表温度、夜间地表温度、昼夜地表温度差、数字高程模型和植被指数重采样为第二分辨率，是通过计算第二分辨率像元范围内，所有第一分辨率像元的平均值来实现的，即计算25km*25km范围内，每个1km*1km的平均值。

在一种实施例中，将待检测区域第一分辨率的降水数据重采样至第二分辨率，是通过计算第二分辨率像元范围内，所有第一分辨率像元的平均值来实现的，即计算25km*25km范围内，每个1km*1km的平均值。

请参阅图3，在一个具体的实施例中，本发明空间降尺度降水数据检测方法包括如下步骤：

步骤s301：获取样本区域的白天地表温度和夜间地表温度，计算昼夜地表温度差。

步骤s302：获取样本区域的数字高程模型和样本区域的植被指数。

步骤s303：将所述样本区域的第一分辨率的白天地表温度、夜间地表温度、昼夜地表温度差、数字高程模型和植被指数进行重采样。

步骤s304：根据样本区域的白天地表温度、夜间地表温度、昼夜地表温度差、数字高程模型和植被指数建立原始样本集s。

步骤s305：通过bootstrap方法在原始样本集s中抽取k个训练样本集。

步骤s306：对k个训练集进行学习，以此生成k个决策树模型。在决策树生成过程中，共有4个输入变量，从4个变量中随机抽取n个变量，各个内部节点均是利用这n个特征变量上最优的分裂方式来分裂，且n值在随机森林模型的形成过程中为恒定常数。

步骤s307：将k个决策树的结果进行组合，经反复训练，形成最优随机森林算法模型。

步骤s308：获取待检测区域的白天地表温度和夜间地表温度，计算待检测区域的昼夜地表温度差。

步骤s209：获取待检测区域的数字高程模型和待检测区域的植被指数。

步骤s210：将所述待检测区域的第一分辨率的白天地表温度、夜间地表温度、昼夜地表温度差、数字高程模型和植被指数输入最优随机森林算法模型，获取待检测区域第一分辨率的降水数据。

步骤s211：将待检测区域第一分辨率的降水数据重采样至第二分辨率。

步骤s212：计算待检测区域卫星遥感降水数据与重采样后的待检测区域第二分辨率的降水数据之间的差值，获取第一残差。

步骤s213：将所述第一残差空间内插至第一分辨率，获取第一分辨率的第二残差。

步骤s214：将待检测区域第一分辨率的降水数据与第二残差相加，获取待检测区域的空间降尺度降水数据。

本发明所述的空间降尺度降水数据检测方法，利用有卫星遥感降水数据的样本区域的白天地表温度、夜间地表温度、昼夜地表温度差、数字高程模型和植被指数以及样本区域的遥感降水数据，建立最优的随机森林算法模型，利用该模型，计算出待检测区域的高分辨率的降水数据，并对该降水数据进行重采样，与卫星遥感降水数据进行残差校正，得到高精度的空间降尺度的降水数据。

请参阅图4，在一个实施例中，本发明空间降尺度降水数据检测装置40包括：

第一数据采集模块41，用于获取样本区域的白天地表温度和夜间地表温度，并根据所述白天地表温度和夜间地表温度获取样本区域的昼夜地表温度差，以及获取样本区域的数字高程模型和样本区域的植被指数；

第一重采样模块42，用于将第一分辨率的样本区域的白天地表温度、夜间地表温度、昼夜地表温度差、数字高程模型和植被指数重采样为与所述样本区域的卫星遥感降水数据相同的第二分辨率；

随机森林训练模块43，用于根据样本区域第二分辨率的白天地表温度、夜间地表温度、昼夜地表温度差、数字高程模型和植被指数，以及所述样本区域的卫星遥感降水数据，建立并训练最优随机森林算法模型；

第二数据采集模块44，用于获取待检测区域的白天地表温度和夜间地表温度，并根据所述待检测区域的白天地表温度和夜间地表温度获取待检测区域的昼夜地表温度差，以及获取待检测区域的数字高程模型和待检测区域的植被指数；

第一降水数据获取模块45，用于将所述第一分辨率的待检测区域白天地表温度、夜间地表温度、昼夜地表温度差、数字高程模型和植被指数输入最优随机森林算法模型，获取待检测区域第一分辨率的降水数据；

残差校正模块46，用于对所述第一分辨率的降水数据进行残差校正，获取待检测区域的空间降尺度降水数据。

在一种实施例中，残差校正模块46包括：

第二重采样单元461，用于将待检测区域第一分辨率的降水数据重采样至第二分辨率；

第一残差获取单元462，用于计算待检测区域卫星遥感降水数据与重采样后的待检测区域第二分辨率的降水数据之间的差值，获取第一残差；

第二残差获取单元463，用于将所述第一残差空间内插至第一分辨率，获取第一分辨率的第二残差；

降尺度降水数据获取单元464，用于将待检测区域第一分辨率的降水数据与第二残差相加，获取待检测区域的空间降尺度降水数据。

在一种实施例中，第一重采样模块42包括第一分辨率计算单元421，用于计算第二分辨率像元范围内，所有第一分辨率像元的待检测区域白天地表温度、夜间地表温度、昼夜地表温度差、数字高程模型和植被指数的平均值。

在一种实施例中，第二重采样单元461包括第二分辨率计算单元4611，用于计算第二分辨率像元范围内，所有第一分辨率像元的第一残差的平均值。

本发明还提供一种计算机可读介质，其上存储有计算机程序，该计算机程序被处理器执行时实现上述任意一项实施例中的基于随机森林算法的降水数据估算方法。

请参阅图5，在一种实施例中，本发明的电子设备50包括存储器51和处理器52，以及储存在所述存储器51并可被所述处理器52执行的计算机程序，所述处理器52执行所述计算机程序时，实现如上述任意一项实施例中的基于随机森林算法的降水数据估算方法。

在本实施例中，控制器52可以是一个或多个应用专用集成电路(asic)、数字信号处理器(dsp)、数字信号处理设备(dspd)、可编程逻辑器件(pld)、现场可编程门阵列(fpga)、控制器、微控制器、微处理器或其他电子元件。存储介质51可采用在一个或多个其中包含有程序代码的存储介质(包括但不限于磁盘存储器、cd-rom、光学存储器等)上实施的计算机程序产品的形式。计算机可读储存介质包括永久性和非永久性、可移动和非可移动媒体，可以由任何方法或技术来实现信息存储。信息可以是计算机可读指令、数据结构、程序的模块或其他数据。计算机的存储介质的例子包括但不限于：相变内存(pram)、静态随机存取存储器(sram)、动态随机存取存储器(dram)、其他类型的随机存取存储器(ram)、只读存储器(rom)、电可擦除可编程只读存储器(eeprom)、快闪记忆体或其他内存技术、只读光盘只读存储器(cd-rom)、数字多功能光盘(dvd)或其他光学存储、磁盒式磁带，磁带磁磁盘存储或其他磁性存储设备或任何其他非传输介质，可用于存储可以被计算设备访问的信息。

所属领域的技术人员可以清楚地了解到，为描述的方便和简洁，上述描述的系统，装置和单元的具体工作过程，可以参考前述方法实施例中的对应过程，在此不再赘述。

在本申请所提供的几个实施例中，应该理解到，所揭露的系统，装置和方法，可以通过其它的方式实现。例如，以上所描述的装置实施例仅仅是示意性的，例如，所述单元的划分，仅仅为一种逻辑功能划分，实际实现时可以有另外的划分方式，例如多个单元或组件可以结合或者可以集成到另一个系统，或一些特征可以忽略，或不执行。另一点，所显示或讨论的相互之间的耦合或直接耦合或通信连接可以是通过一些接口，装置或单元的间接耦合或通信连接，可以是电性，机械或其它的形式。

所述作为分离部件说明的单元可以是或者也可以不是物理上分开的，作为单元显示的部件可以是或者也可以不是物理单元，即可以位于一个地方，或者也可以分布到多个网络单元上。可以根据实际的需要选择其中的部分或者全部单元来实现本实施例方案的目的。

以上所述实施例仅表达了本发明的几种实施方式，其描述较为具体和详细，但并不能因此而理解为对发明专利范围的限制。应当指出的是，对于本领域的普通技术人员来说，在不脱离本发明构思的前提下，还可以做出若干变形和改进，这些都属于本发明的保护范围。