CN115659853B - 一种非线性混合效应变系数降尺度方法和系统 - Google Patents
一种非线性混合效应变系数降尺度方法和系统 Download PDFInfo
- Publication number
- CN115659853B CN115659853B CN202211688080.5A CN202211688080A CN115659853B CN 115659853 B CN115659853 B CN 115659853B CN 202211688080 A CN202211688080 A CN 202211688080A CN 115659853 B CN115659853 B CN 115659853B
- Authority
- CN
- China
- Prior art keywords
- precipitation
- model
- coefficient
- downscaling
- resolution
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
- 238000000034 method Methods 0.000 title claims abstract description 72
- 238000001556 precipitation Methods 0.000 claims abstract description 207
- 230000000694 effects Effects 0.000 claims abstract description 29
- 238000007637 random forest analysis Methods 0.000 claims abstract description 29
- 230000007613 environmental effect Effects 0.000 claims abstract description 14
- 230000004044 response Effects 0.000 claims abstract description 13
- 239000013598 vector Substances 0.000 claims description 53
- 239000011159 matrix material Substances 0.000 claims description 45
- 230000003993 interaction Effects 0.000 claims description 12
- 238000000354 decomposition reaction Methods 0.000 claims description 8
- 230000008859 change Effects 0.000 claims description 7
- 238000001914 filtration Methods 0.000 claims description 7
- 230000009466 transformation Effects 0.000 claims description 6
- 238000010276 construction Methods 0.000 claims description 4
- 238000002790 cross-validation Methods 0.000 claims description 4
- 230000001131 transforming effect Effects 0.000 claims description 4
- 238000007476 Maximum Likelihood Methods 0.000 claims description 3
- 230000002301 combined effect Effects 0.000 claims 1
- 238000002156 mixing Methods 0.000 claims 1
- 230000009467 reduction Effects 0.000 description 12
- 230000008569 process Effects 0.000 description 11
- 238000004088 simulation Methods 0.000 description 6
- 238000003066 decision tree Methods 0.000 description 5
- 238000012216 screening Methods 0.000 description 5
- 238000010586 diagram Methods 0.000 description 4
- 238000011156 evaluation Methods 0.000 description 4
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 4
- 230000001276 controlling effect Effects 0.000 description 3
- 238000011160 research Methods 0.000 description 3
- 238000004458 analytical method Methods 0.000 description 2
- 230000008901 benefit Effects 0.000 description 2
- 238000004422 calculation algorithm Methods 0.000 description 2
- 238000004364 calculation method Methods 0.000 description 2
- 230000002596 correlated effect Effects 0.000 description 2
- 230000007812 deficiency Effects 0.000 description 2
- 230000004907 flux Effects 0.000 description 2
- 238000007726 management method Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000005855 radiation Effects 0.000 description 2
- 238000012952 Resampling Methods 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 239000008358 core component Substances 0.000 description 1
- 230000008878 coupling Effects 0.000 description 1
- 238000010168 coupling process Methods 0.000 description 1
- 238000005859 coupling reaction Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000002474 experimental method Methods 0.000 description 1
- 239000000284 extract Substances 0.000 description 1
- 230000006870 function Effects 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 230000010354 integration Effects 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 230000007246 mechanism Effects 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 239000002245 particle Substances 0.000 description 1
- 230000002093 peripheral effect Effects 0.000 description 1
- 239000005436 troposphere Substances 0.000 description 1
Images
Classifications
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
- Y02A90/00—Technologies having an indirect contribution to adaptation to climate change
- Y02A90/10—Information and communication technologies [ICT] supporting adaptation to climate change, e.g. for weather forecasting or climate simulation
Landscapes
- Complex Calculations (AREA)
Abstract
本申请涉及应用电子设备进行识别的方法或装置技术领域,提供一种非线性混合效应变系数降尺度方法和系统。该方法包括:获取降水影响因素数据集和第一分辨率降水遥感数据;基于非线性混合效应变系数模型和随机森林模型,根据第一分辨率降水遥感数据,构建降水空间降尺度模型;其中,随机森林模型用于从降水影响因素数据集中确定降水关键主控要素;降水关键主控要素用于组成降水空间降尺度模型的解释变量;非线性混合效应变系数模型是在空间变系数回归模型引入随机效应,并将空间变系数回归模型中的局部线性特征扩展为随环境要素的非线性响应特征得到的;基于降水空间降尺度模型对第一分辨率降水遥感数据进行降尺度,得到第二分辨率降水遥感数据。
Description
技术领域
本申请涉及应用电子设备进行识别的方法或装置技术领域,特别涉及一种非线性混合效应变系数降尺度方法和系统。
背景技术
降水是描述某地区乃至全球气候系统变化的关键指标,是全球水分和能量循环的核心组成部分。作为流域生态水文过程的关键输入参数,降水分布的准确估算直接影响着流域内水文水资源分析、水资源规划与管理及生态环境治理等。高精度高分辨率的空间化降水信息对认识流域生态系统与水文过程及其相互作用的途径和机理,建立流域生态水文过程模拟模型和水资源决策支持系统,提高流域水-生态-经济耦合系统演变的综合分析与预测能力等方面具有重要的科学意义和应用价值。
传统的降水获取方式主要基于站点数据采用一定的插值方法进行计算。然而,基于站点的插值方法往往受限于站点的分布及密度限制,使得插值方法具有很大的不确定性。随着遥感卫星事业的快速发展及现代获取空间数据技术手段的进步,当前出现了大量的卫星遥感数据产品,这些数据为站点系数、地形复杂地区提供了连续的降水分布信息,能一定程度上弥补缺测资料地区地面观测信息的不足,已被广泛应用于不同研究领域。然而受到传感器性能、云层特性、反演算法等所限,其定量误差比较突出。此外,遥感数据往往空间分辨率较低,不能满足精细尺度局地模拟的需求。
将低分辨率信息通过一定技术手段得到高分辨率数据信息并同时提高其精度的过程被称之为降尺度(Downscaling)。降尺度方法大体上可以分为动力降尺度和统计降尺度两大类,相对于动力降尺度方法,统计降尺度因其具有计算量小节省时间、方法众多形式灵活、更易于操作等优点而被广泛应用于卫星遥感降水产品的降尺度研究中。在以往统计降尺度方法中,地理加权回归模型是最具代表性应用最为广泛的降尺度方法之一。地理加权回归作为局部空间变系数回归方法能够解析地理环境要素的空间异质性及空间变化的依赖关系,但其受共线性影响显著,并且该方法要求回归系统具有一致的平滑度。而实际问题中,降水与环境变量的关系并不一定是空间非平稳,这使得回归系数的变异程度有所差异。
因此,需要提供一种针对上述现有技术不足的全新技术方案。
发明内容
本申请的目的在于提供一种非线性混合效应变系数降尺度方法和系统,以解决或缓解上述现有技术中存在的问题。
为了实现上述目的,本申请提供如下技术方案:
本申请提供了一种非线性混合效应变系数降尺度方法,包括:
获取降水影响因素数据集和第一分辨率降水遥感数据;
基于非线性混合效应变系数模型和随机森林模型,根据所述第一分辨率降水遥感数据,构建降水空间降尺度模型;
其中,所述随机森林模型用于从所述降水影响因素数据集中确定降水关键主控要素;所述降水关键主控要素用于组成所述降水空间降尺度模型的解释变量;所述非线性混合效应变系数模型是在空间变系数回归模型引入随机效应,并将所述空间变系数回归模型中的局部线性特征扩展为随环境要素的非线性响应特征得到的;
基于所述降水空间降尺度模型对所述第一分辨率降水遥感数据进行降尺度,得到第二分辨率降水遥感数据;其中,所述第二分辨率降水遥感数据的空间分辨率高于所述第一分辨率降水遥感数据的空间分辨率。
优选地,所述非线性混合效应变系数模型的表达式如下:
式中,y为待求变量;x p 为第p个解释变量;α p 为第p个回归系数;P为解释变量的个数;x i 为第i个解释变量;x j 为第j个解释变量;为哈达玛积;为空间变异系数构成的向量;ε为残差;N表示正态分布;为方差;I为单位矩阵;为系数向量;为个特征向量构成的矩阵;为常数项;为变化项;控制第个回归系数的空间平滑度,控制回归系数的变异程度。
优选地,所述非线性混合效应变系数模型是在空间变系数回归模型引入随机效应,并将所述空间变系数回归模型中的局部线性特征扩展为随环境要素的非线性响应特征得到的;
所述空间变系数回归模型的表达式如下:
式中:y为待求变量;x p 为第p个解释变量;P为解释变量的个数;为空间变异系数构成的向量;为哈达玛积;为空间变异系数构成的向量;ε为残差;N表示正态分布;为方差;I为单位矩阵;为系数向量;为个特征向量构成的矩阵;为常数项;为变化项。
优选地,所述空间变系数回归模型为基于特征向量分解的空间滤波模型构建的;
所述基于特征向量分解的空间滤波模型的表达式如下:
优选地,所述降水影响因素数据集至少包括:地理地形要素、天气系统要素、第一分辨率降水遥感数据,以及由所述地理地形要素、所述天气系统要素、所述第一分辨率降水遥感数据形成的数据集中两两组合得到的交互项影响因素。
优选地,所述基于所述降水空间降尺度模型对所述第一分辨率降水遥感数据进行降尺度,得到第二分辨率降水遥感数据,具体为:
将所述降水空间降尺度模型进行变换处理,得到降水空间降尺度模型变换表达式:
式中:y为待求变量,为解释变量构成的矩阵,D为特征向量构成的矩阵,x p 为第p个解释变量,P为解释变量的个数,α为回归系数,ε为残差,为哈达玛积,0 LP 为值为0的向量,为哈达玛积,I LP 为阶单位矩阵,为对角矩阵,其第个对角线元素为;为方差;
基于残差最大似然方法,对降水空间降尺度模型变换表达式中的参数进行求解,以得到第二分辨率降水遥感数据。
优选地,还包括:
基于气象站点观测数据,对所述第二分辨率降水遥感数据进行交叉验证,以检验所述第二分辨率降水遥感数据的精度。
本申请实施例还提供一种非线性混合效应变系数降尺度系统,包括:
获取单元,配置为获取降水影响因素数据集和第一分辨率降水遥感数据;
构建单元,配置为基于非线性混合效应变系数模型和随机森林模型,根据所述第一分辨率降水遥感数据,构建降水空间降尺度模型;
其中,所述随机森林模型用于从所述降水影响因素数据集中确定降水关键主控要素;所述降水关键主控要素用于组成所述降水空间降尺度模型的解释变量;所述非线性混合效应变系数模型是在空间变系数回归模型引入随机效应,并将所述空间变系数回归模型中的局部线性特征扩展为随环境要素的非线性响应特征得到的;
降尺度单元,配置为基于所述降水空间降尺度模型对所述第一分辨率降水遥感数据进行降尺度,得到第二分辨率降水遥感数据;其中,所述第二分辨率降水遥感数据的空间分辨率高于所述第一分辨率降水遥感数据的空间分辨率。
有益效果:
本申请提供的技术方案中,通过在空间变系数回归模型引入随机效应,并将空间变系数回归模型中的局部线性特征扩展为随环境要素的非线性响应特征,得到非线性混合效应变系数模型,使得所得到的非线性混合效应变系数模型中各回归系数具有不同程度的空间变异特征及光滑性,能够很好地刻画出由于区域差异导致的变量空间异质性,同时降低了局部线性特征对非线性混合效应变系数模型的影响,能够很好地刻画随环境要素的非线性响应特征,进而提高了降尺度的精度;在此基础上,结合随机森林模型筛选降水关键主控要素,构建基于非线性混合效应变系数模型和随机森林模型的降水空间降尺度模型;通过该降水空间降尺度模型对低分辨率的降水遥感数据进行降尺度处理,得到高分辨率的降水遥感数据。
附图说明
构成本申请的一部分的说明书附图用来提供对本申请的进一步理解,本申请的示意性实施例及其说明用于解释本申请,并不构成对本申请的不当限定。其中:
图1为根据本申请的一些实施例提供的非线性混合效应变系数降尺度方法的逻辑示意图;
图2为根据本申请的一些实施例提供的非线性混合效应变系数降尺度方法的流程示意图;
图3为根据本申请的一些实施例提供的基于随机森林方法筛选主控要素的逻辑示意图;
图4为根据本申请的一些实施例提供的非线性混合效应变系数降尺度系统的结构示意图。
具体实施方式
下面将参考附图并结合实施例来详细说明本申请。
在以下描述中,所涉及的术语“第一/第二/第三”仅仅是区别类似的对象,不代表对对象的特定排序。
除另有定义,本文所使用的所有的技术和科学术语与属于本公开的技术领域的技术人员通常理解的含义相同。本文中所使用的术语只是为了描述本公开实施例的目的,不是旨在限制本公开。
为了便于理解本申请的技术方案,下面对相关技术进行简要说明。
如背景技术所述,降尺度能够提高气候要素的空间分辨率。动力降尺度法是建立在区域气候模式(Regional Climate Model,RCM)的基础之上的,使用全球模式(GlobalClimate Model,GCM)为区域气候模式提供初始和边界条件,通过高分辨率区域气候模式的数值积分获得高分辨率降尺度结果。其中,区域气候模式采用数学物理方程描述气候系统内部的各种动力和热力学过程。虽然动力降尺度方法能够较好地保持气候要素各个变量之间的依赖关系,但是其需要大量的计算资源,随着模式分辨率的提高,计算量呈指数形式增长。统计降尺度方法是通过转换函数,建立大尺度(低分辨率/粗分辨率)降水因子与区域尺度气候预报变量之间的统计函数关系。其计算量小,方法众多形式灵活、更易于操作,能够节约大量的时间。
为了克服背景技术所述的传统空间变系数回归模型(比如地理加权回归模型,Geographically Weighted Regression,GWR)受共线性影响显著以及不能反映区域差异导致的变量空间异质性的缺陷,本申请实施例提供一种非线性混合效应变系数降尺度方法和系统,通过在空间变系数回归模型引入随机效应,并将空间变系数回归模型中的局部线性特征扩展为随环境要素的非线性响应特征,构建非线性混合效应变系数模型(NMSVC),从而使得NMSVC模型中各回归系数具有不同程度的空间变异特征及光滑性,能够精准刻画区域差异造成的地理环境异质性,同时能很好的刻画随环境要素的非线性响应特征,从而有效地提高建模精度。在此基础上,结合随机森林(RandomForest,RF)模型筛选关键主控要素,构建基于RF和NMSVC的降水空间降尺度模型,并基于该降水空间降尺度模型,以低分辨率降水遥感数据作为输入得到高分辨率降水遥感数据,极大提高了降水数据降尺度结果的精度。
示例性方法
本申请实施例提供一种非线性混合效应变系数降尺度方法,如图1至图3所示,该方法包括:
步骤S101、获取降水影响因素数据集和第一分辨率降水遥感数据。
基于前述说明可知,降尺度是将低分辨率信息通过一定技术手段得到高分辨率数据信息并同时提高其精度的过程,本申请实施例所提供的方法,能够应用于各种低分辨率(也称为粗分辨率)地理环境要素的降尺度模拟,为地学环境要素模拟提供新的方案参考。为了便于理解,以将低分辨率降水数据降尺度为高精度降水数据为例,对本申请提供的方法进行详细说明。
为了便于描述,将低分辨率的降水遥感数据称为第一分辨率降水遥感数据,将精度提高后得到的高分辨率降水遥感数据称为第二分辨率降水遥感数据。
考虑到时间尺度以及各个遥感卫星的表现性能差异,本申请实施例选择GSMaP逐日降水产品作为低分辨率降水遥感数据(即第一分辨率降水遥感数据),对其进行降尺度,以得到高分辨率降水遥感数据。
需要说明的是,GSMaP逐日降水产品是全球降水观测计划(Global PrecipitationMeasurement mission,GPM)提供的数据产品。其中,GPM携带多种卫星传感器,能够为遥感水文科学社区提供高时空分辨率(1-hour,10-km)和更高精度的全球卫星降水观测。GSMaP逐日降水产品是利用多传感器多卫星多算法结合卫星网络和雨量计反演得到的降水数据产品。
在对粗分辨率降水数据进行降尺度之前,需要先确定模型的解释变量的备选集合,即降水影响因素,为此,一些实施例中,降水影响因素数据集至少包括:地理地形要素、天气系统要素、第一分辨率降水遥感数据,以及由地理地形要素、天气系统要素、第一分辨率降水遥感数据形成的数据集中两两组合得到的交互项影响因素。
具体地,以逐日降水的降尺度模拟为例进行说明。传统的降尺度方法中,降水影响因素(即解释变量)的选取主要集中在地理地形等要素,例如经度、纬度、海拔、坡度等。而经度、纬度、海拔、坡度等地理地形特征往往与长时间尺度(年尺度甚至更长时间)气候特征有关,很难精确刻画天气尺度现象,不足以反应局地尺度降水与陆表、大气之间的交互特征。特别针对日尺度降水,由于受大气环流特征影响显著,逐日降水的降尺度模拟需要在以往降尺度研究基础上考虑其他的天气层面影响因素。为此,本申请实施例中,针对逐日降水影响因素,在传统降水影响因素(解释变量)基础上增加了天气系统层面因素、粗分辨率降水数据(即第一分辨率降水遥感数据),以及各要素交互项,以构成本申请实施例的解释变量的备选集合。
图3示出了本申请实施例中降水影响因素的筛选过程,参见图3,对于降水而言,传统影响因素包括地理、地形要素,本申请实施例新增天气系统要素、粗分辨率降水数据以及各要素交互项,作为降水影响因素。应当理解,每一种降水影响因素均包括一种或多种降水影响要素,比如降水影响因素为地理、地形,其又具体包括经纬度、海拔、坡度等;天气系统要素包括:云量、云光学厚度、云粒子有效半径、云顶温度、云顶气压、云水路径、对流层中下层500hPa、600 hPa、700 hPa、800 hPa、850 hPa、900 hPa、950 hPa及1000 hPa的位势高度、空气温度、潜热通量、感热通量、短波辐射、长波辐射、相对湿度、最大相对湿度、最小相对湿度、比湿、海平面气压、风速、高程、坡度、经度、纬度、到海岸线的距离等;粗分辨率降水数据包括叶面积指数NDVI、粗分辨率网格的降水值及其周边网格的降水值等。将上述多种降水影响要素组成初始降水影响因素数据集。此外,本申请还引入上述多种降水影响要素相互之间的交互项,用于反应局地尺度降水与陆表、大气之间的交互特征,以进一步提高降尺度的精度。
其中,多种降水影响要素相互之间的交互项通过地理地形要素、天气系统要素、粗分辨率降水数据组成的初始降水影响因素数据集中的要素两两组合得到。在此基础上,将降水影响因素数据集中各个降水影响要素作为备选解释变量,用于后续步骤中构建降水空间降尺度模型。
步骤S102、基于非线性混合效应变系数模型和随机森林模型,根据第一分辨率降水遥感数据,构建降水空间降尺度模型。
其中,随机森林模型用于从降水影响因素数据集中确定降水关键主控要素;降水关键主控要素用于组成降水空间降尺度模型的解释变量;非线性混合效应变系数模型是在空间变系数回归模型引入随机效应,并将空间变系数回归模型中的局部线性特征扩展为随环境要素的非线性响应特征得到的。
下面对本申请的非线性混合效应变系数模型进行详细说明。
本申请实施例中,非线性混合效应变系数模型是在空间变系数回归模型的基础上构建的,一些实施例中,空间变系数回归模型为基于特征向量分解的空间滤波模型构建的。
具体地,特征向量空间滤波是将空间自相关的变量通过移除空间模式转化为空间无关的变量的方法,使用地理连接矩阵特征向量的线性组合,可以将原变量分为两个合成部分,一个空间相关部分和一个非空间相关部分。这样,这个变量可以独立于观察相邻区域的影响进行研究。基于特征向量分解的空间滤波模型的表达式如下:
进一步地,上述方法可发展为一种空间变系数回归模型:
对式(2)进行变换处理,可以得到最终的空间变系数回归模型,表达式如下:
式中:y为待求变量;x p 为第p个解释变量;P为解释变量的个数;为空间变异系数构成的向量;为哈达玛积;为空间变异系数构成的向量;ε为残差;N表示正态分布;为方差;I为单位矩阵;为系数向量;为个特征向量构成的矩阵;为常数项;为变化项。
然而,上述公式(3)所表达的空间变系数回归模型中,特征向量分解模型为固定效应模型,不能很好的刻画由于区域差异所导致的变量空间异质性,为此申请实施例将随机效应引入上述公式(3)所表达的空间变系数回归模型,构建混合效应变系数模型,以精准刻画区域差异造成的地理环境异质性。同时,考虑到当前空间变系数回归模型的局部线性特征,并不能很好的刻画随环境要素的非线性响应特征,进一步的将所构建的混合效应变系数模型扩展为非线性混合效应变系数模型,最终构建了基于特征向量分解的非线性混合效应变系数模型(NMSVC),以有效提高建模精度:
本申请实施例的非线性混合效应变系数模型用公式(4)表示,公式(4)如下:
式中,y为待求变量;x p 为第p个解释变量;α p 为第p个回归系数;P为解释变量的个数;x i 为第i个解释变量;x j 为第j个解释变量;为哈达玛积;为空间变异系数构成的向量;ε为残差;N表示正态分布;为方差;I为单位矩阵;为系数向量;为个特征向量构成的矩阵;为常数项;为变化项;控制第个回归系数的空间平滑度,用于控制回归系数的变异程度。
下面对从解释变量的备选集合中筛选降水关键主控要素的过程进行详细说明。
本申请实施例中,随机森林模型用于从地理地形、天气系统、粗分辨率降水遥感数据等包含的降水影响要素以及各要素两辆之间的交互项组成的降水影响因素数据集中确定降水关键主控要素;降水关键主控要素用于组成降水空间降尺度模型的解释变量。
其中,随机森林是一种利用自助法重抽样方法从原始样本中抽取多个样本,组合很多个相对独立的决策树,建立决策树“森林”,进行“投票”决定最终预测结果的统计学习方法。本申请实施例中,降水影响因素由多个数据集组成,每个数据集对应一个解释变量(即降水影响要素),为了减少特征维度,采用随机森林对降水影响因素数据集进行筛选,以确定对模型影响较大的解释变量,以此来减少降水空间降尺度模型的计算量。
具体的,基于随机森林模型对降水影响因素数据集进行筛选可以包括以下步骤:将降水影响因素数据集输入随机森林模型,并判断降水影响因素数据集中每个降水影响要素在随机森林中每棵决策树上的贡献度,接着基于每个降水影响要素对每棵决策树的贡献度,计算贡献度均值,随后对降水影响要素的贡献度均值进行对比,选取贡献度均值降序排序后得到的排序列表中位置处于靠前的降水影响要素作为降水关键主控要素。
示例性地,可以采用随机森林模型对降水影响因素数据集中每个降水影响要素进行重要性评估,并量化每个降水影响要素对各个决策树分类性能的贡献度,然后对贡献度进行排序,筛选出降序排序后前10的降水影响要素作为降水关键主控要素,并由该10个降水关键主控要素组成降水空间降尺度模型的解释变量。
其中,贡献度可以采用基尼指数或袋外(OutOfBag,OOB)数据错误率作为评估指标,或者其他的评估指标,本申请对贡献度评估指标的选择不作限定。
步骤S103、基于降水空间降尺度模型对第一分辨率降水遥感数据进行降尺度,得到第二分辨率降水遥感数据。
其中,第二分辨率降水遥感数据的空间分辨率高于第一分辨率降水遥感数据的空间分辨率。
具体地,在进行降尺度时,基于降水空间降尺度模型对第一分辨率降水遥感数据进行降尺度得到第二分辨率降水遥感数据可以包括模型变换的步骤和模型求解的步骤。详细步骤如下:
首先对公式(4)所表达的降水空间降尺度模型进行改写,得到公式(5),公式(5)如下:
然后,将降水空间降尺度模型进行进一步变换处理,得到降水空间降尺度模型变换表达式:
式中:y为待求变量,为解释变量构成的矩阵,D为特征向量构成的矩阵,x p 为第p个解释变量,P为解释变量的个数,α为回归系数,ε为残差,0 LP 为值为0的向量,为哈达玛积,I LP 为阶单位矩阵,为对角矩阵,其第个对角线元素为;为方差。
最后,基于残差最大似然方法,对变换后表达式(6)中的参数进行求解,以得到第二分辨率降水遥感数据。具体可以由公式(7)求解得到,公式(7)如下:
一些实施例中,该方法还包括:基于气象站点观测数据,对第二分辨率降水遥感数据进行交叉验证,以检验第二分辨率降水遥感数据的精度。
具体地,基于上述降水空间降尺度模型,对GSMaP降水产品进行降尺度,并基于交叉验证方法将降尺度结果与气象站点观测数据进行比较,同时与经典的基于地理加权回归的降尺度方法进行对比。实验表明,本申请提供的方法能够大幅度提高降尺度精度,从而为地理环境要素降尺度以及地学要素的建模提供了一种更有效的、全新的局部回归方法。
综上所述,本申请实施例中,通过在空间变系数回归模型引入随机效应,并将空间变系数回归模型中的局部线性特征扩展为随环境要素的非线性响应特征,得到非线性混合效应变系数模型,使得所得到的非线性混合效应变系数模型中各回归系数具有不同程度的空间变异特征及光滑性,能够很好地刻画出由于区域差异导致的变量空间异质性,同时降低了局部线性特征对非线性混合效应变系数模型的影响,能够很好地刻画随环境要素的非线性响应特征,进而提高了降尺度的精度;在此基础上,结合随机森林模型筛选降水关键主控要素,构建基于非线性混合效应变系数模型和随机森林模型的降水空间降尺度模型;通过该降水空间降尺度模型对低分辨率的降水遥感数据进行降尺度处理,得到高分辨率的降水遥感数据,该方法是一种更有效的、全新的局部回归方法,能够应用于地理环境要素降尺度以及地学要素的建模,并能提高降尺度和地学要素的建模精度。
示例性系统
本申请实施例提供一种非线性混合效应变系数降尺度系统,如图4所示,该系统包括:获取单元401、构建单元402和降尺度单元403。其中:
获取单元401,配置为获取降水影响因素数据集和第一分辨率降水遥感数据。
构建单元402,配置为基于非线性混合效应变系数模型和随机森林模型,根据第一分辨率降水遥感数据,构建降水空间降尺度模型;
其中,随机森林模型用于从降水影响因素数据集中确定降水关键主控要素;降水关键主控要素用于组成降水空间降尺度模型的解释变量;非线性混合效应变系数模型是在空间变系数回归模型引入随机效应,并将空间变系数回归模型中的局部线性特征扩展为随环境要素的非线性响应特征得到的。
降尺度单元403,配置为基于降水空间降尺度模型对第一分辨率降水遥感数据进行降尺度,得到第二分辨率降水遥感数据;其中,第二分辨率降水遥感数据的空间分辨率高于第一分辨率降水遥感数据的空间分辨率。
一些实施例中,非线性混合效应变系数模型的表达式如下:
式中,y为待求变量;x p 为第p个解释变量;α p 为第p个回归系数;P为解释变量的个数;x i 为第i个解释变量;x j 为第j个解释变量;为哈达玛积;为空间变异系数构成的向量;ε为残差;N表示正态分布;为方差;I为单位矩阵;为系数向量;为个特征向量构成的矩阵;为常数项;为变化项;控制第个回归系数的空间平滑度,控制回归系数的变异程度。
一些实施例中,空间变系数回归模型的表达式如下:
式中:y为待求变量;x p 为第p个解释变量;P为解释变量的个数;为空间变异系数构成的向量;为哈达玛积;为空间变异系数构成的向量;ε为残差;N表示正态分布;为方差;I为单位矩阵;为系数向量;为个特征向量构成的矩阵;为常数项;为变化项。
本申请实施例提供的非线性混合效应变系数降尺度系统能够实现上述任一实施例提供的非线性混合效应变系数降尺度方法的流程、步骤,并达到相同的技术效果,在此不做一一赘述。
以上所述仅为本申请的优选实施例,并不用于限制本申请,对于本领域的技术人员来说,本申请可以有各种更改和变化。凡在本申请的精神和原则之内,所作的任何修改、等同替换、改进等,均应包含在本申请的保护范围之内。
Claims (8)
1.一种非线性混合效应变系数降尺度方法,其特征在于,包括:
获取降水影响因素数据集和第一分辨率降水遥感数据;
基于非线性混合效应变系数模型和随机森林模型,根据所述第一分辨率降水遥感数据,构建降水空间降尺度模型;
其中,所述随机森林模型用于从所述降水影响因素数据集中确定降水关键主控要素;所述降水关键主控要素用于组成所述降水空间降尺度模型的解释变量;所述非线性混合效应变系数模型是在空间变系数回归模型引入随机效应,并将所述空间变系数回归模型中的局部线性特征扩展为随环境要素的非线性响应特征得到的;
所述非线性混合效应变系数模型的表达式如下:
式中:y为待求变量,x p 为第p个解释变量,α p 为第p个回归系数,P为解释变量的个数,x i 为第i个解释变量,x j 为第j个解释变量,为哈达玛积,为非线性混合效应变系数模型的空间变异系数构成的向量,ε为残差,N表示正态分布,为方差,I为单位矩阵,为系数向量,为个特征向量构成的矩阵,为常数项,为变化项,控制第个回归系数的空间平滑度,控制回归系数的变异程度;
基于所述降水空间降尺度模型对所述第一分辨率降水遥感数据进行降尺度,得到第二分辨率降水遥感数据;其中,所述第二分辨率降水遥感数据的空间分辨率高于所述第一分辨率降水遥感数据的空间分辨率。
4.根据权利要求1所述的非线性混合效应变系数降尺度方法,其特征在于,所述降水影响因素数据集至少包括:地理地形要素、天气系统要素、第一分辨率降水遥感数据,以及由所述地理地形要素、所述天气系统要素、所述第一分辨率降水遥感数据形成的数据集中两两组合得到的交互项影响因素。
5.根据权利要求1所述的非线性混合效应变系数降尺度方法,其特征在于,所述基于所述降水空间降尺度模型对所述第一分辨率降水遥感数据进行降尺度,得到第二分辨率降水遥感数据,具体为:
将所述降水空间降尺度模型进行变换处理,得到降水空间降尺度模型变换表达式:
式中:y为待求变量,为解释变量构成的矩阵,D为特征向量构成的矩阵,x p 为第p个解释变量,P为解释变量的个数,α为回归系数,ε为残差,0 LP 为值为0的向量,为哈达玛积,I LP 为阶单位矩阵,为对角矩阵,其第个对角线元素为;为方差;
基于残差最大似然方法,对降水空间降尺度模型变换表达式中的参数进行求解,以得到第二分辨率降水遥感数据。
6.根据权利要求1所述的非线性混合效应变系数降尺度方法,其特征在于,还包括:
基于气象站点观测数据,对所述第二分辨率降水遥感数据进行交叉验证,以检验所述第二分辨率降水遥感数据的精度。
7.一种非线性混合效应变系数降尺度系统,其特征在于,包括:
获取单元,配置为获取降水影响因素数据集和第一分辨率降水遥感数据;
构建单元,配置为基于非线性混合效应变系数模型和随机森林模型,根据所述第一分辨率降水遥感数据,构建降水空间降尺度模型;
其中,所述随机森林模型用于从所述降水影响因素数据集中确定降水关键主控要素;所述降水关键主控要素用于组成所述降水空间降尺度模型的解释变量;所述非线性混合效应变系数模型是在空间变系数回归模型引入随机效应,并将所述空间变系数回归模型中的局部线性特征扩展为随环境要素的非线性响应特征得到的;
所述非线性混合效应变系数模型的表达式如下:
式中,y为待求变量;x p 为第p个解释变量;α p 为第p个回归系数;P为解释变量的个数;x i 为第i个解释变量;x j 为第j个解释变量;为哈达玛积;为非线性混合效应变系数模型的空间变异系数构成的向量;ε为残差;N表示正态分布;为方差;I为单位矩阵;为系数向量;为个特征向量构成的矩阵;为常数项;为变化项;控制第个回归系数的空间平滑度,控制回归系数的变异程度;
降尺度单元,配置为基于所述降水空间降尺度模型对所述第一分辨率降水遥感数据进行降尺度,得到第二分辨率降水遥感数据;其中,所述第二分辨率降水遥感数据的空间分辨率高于所述第一分辨率降水遥感数据的空间分辨率。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202211688080.5A CN115659853B (zh) | 2022-12-28 | 2022-12-28 | 一种非线性混合效应变系数降尺度方法和系统 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202211688080.5A CN115659853B (zh) | 2022-12-28 | 2022-12-28 | 一种非线性混合效应变系数降尺度方法和系统 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN115659853A CN115659853A (zh) | 2023-01-31 |
CN115659853B true CN115659853B (zh) | 2023-03-10 |
Family
ID=85022863
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202211688080.5A Active CN115659853B (zh) | 2022-12-28 | 2022-12-28 | 一种非线性混合效应变系数降尺度方法和系统 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN115659853B (zh) |
Families Citing this family (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN115952743B (zh) * | 2023-03-14 | 2023-05-05 | 中国科学院地理科学与资源研究所 | 耦合随机森林和hasm的多源降水数据协同降尺度方法和系统 |
CN116070792B (zh) * | 2023-03-28 | 2023-07-25 | 中国科学院地理科学与资源研究所 | 一种多源降水数据的融合方法、装置、存储介质和设备 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106776481A (zh) * | 2016-11-29 | 2017-05-31 | 河海大学 | 一种作用于卫星降水数据的降尺度校正方法 |
CN109325540A (zh) * | 2018-09-30 | 2019-02-12 | 河南大学 | 一种针对遥感每天降水量数据的空间降尺度方法 |
CN112699959A (zh) * | 2021-01-11 | 2021-04-23 | 中国科学院地理科学与资源研究所 | 基于能量泛函模型的多源多尺度降水数据融合方法和装置 |
CN114461983A (zh) * | 2022-02-09 | 2022-05-10 | 中国科学院、水利部成都山地灾害与环境研究所 | 一种基于水量平衡原理的卫星降水产品空间降尺度方法 |
Family Cites Families (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110738252B (zh) * | 2019-10-14 | 2020-08-14 | 广州地理研究所 | 空间自相关的机器学习卫星降水数据降尺度方法、系统 |
CN114324410A (zh) * | 2021-12-31 | 2022-04-12 | 黄陵县农产品质量安全检验检测站 | 一种多地形的微波遥感土壤湿度降尺度方法 |
-
2022
- 2022-12-28 CN CN202211688080.5A patent/CN115659853B/zh active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106776481A (zh) * | 2016-11-29 | 2017-05-31 | 河海大学 | 一种作用于卫星降水数据的降尺度校正方法 |
CN109325540A (zh) * | 2018-09-30 | 2019-02-12 | 河南大学 | 一种针对遥感每天降水量数据的空间降尺度方法 |
CN112699959A (zh) * | 2021-01-11 | 2021-04-23 | 中国科学院地理科学与资源研究所 | 基于能量泛函模型的多源多尺度降水数据融合方法和装置 |
CN114461983A (zh) * | 2022-02-09 | 2022-05-10 | 中国科学院、水利部成都山地灾害与环境研究所 | 一种基于水量平衡原理的卫星降水产品空间降尺度方法 |
Non-Patent Citations (4)
Title |
---|
Kai Tang等.Spatial Downscaling of Land Surface Temperature over Heterogeneous Regions Using Random Forest Regression Considering Spatial Features.2021,全文. * |
Na Zhao.An efficient downscaling scheme for high-resolution precipitation estimates over a high mountainous watershed.2021,第13卷(第13期),全文. * |
徐彬仁 ; 魏瑗瑗 ; .基于随机森林算法对青藏高原TRMM降水数据进行空间统计降尺度研究.2018,(03),全文. * |
王晨亮 等.高精度曲面建模的中国气候降尺度模型.2012,第14卷(第5期),全文. * |
Also Published As
Publication number | Publication date |
---|---|
CN115659853A (zh) | 2023-01-31 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Feloni et al. | Flood vulnerability assessment using a GIS‐based multi‐criteria approach—The case of Attica region | |
CN115659853B (zh) | 一种非线性混合效应变系数降尺度方法和系统 | |
Agarwal et al. | Delineation of groundwater potential zone: An AHP/ANP approach | |
Kalteh et al. | Review of the self-organizing map (SOM) approach in water resources: Analysis, modelling and application | |
Cherkassky et al. | Computational intelligence in earth sciences and environmental applications: Issues and challenges | |
Langella et al. | High-resolution space–time rainfall analysis using integrated ANN inference systems | |
Solaiman et al. | Development of probability based intensity-duration-frequency curves under climate change | |
CN111665575B (zh) | 一种基于统计动力的中长期降雨分级耦合预报方法及系统 | |
CN114254561A (zh) | 一种内涝预测方法、系统及存储介质 | |
Muneer et al. | Modeling of runoff in the arid regions using remote sensing and geographic information system (GIS) | |
Kanishka et al. | Streamflow estimation in ungauged basins using watershed classification and regionalization techniques | |
CN115795399B (zh) | 一种多源遥感降水数据自适应融合方法和系统 | |
Zhu et al. | Internal and external coupling of Gaussian mixture model and deep recurrent network for probabilistic drought forecasting | |
Mouatadid et al. | A machine learning approach to non-uniform spatial downscaling of climate variables | |
Alhamshry et al. | Prediction of summer rainfall over the source region of the Blue Nile by using teleconnections based on sea surface temperatures | |
Chu et al. | Middle-and long-term streamflow forecasting and uncertainty analysis using Lasso-DBN-Bootstrap model | |
Jaberalansar et al. | Spatial downscaling of climate variables using three statistical methods in Central Iran | |
Devi et al. | Predictive modelling of land use land cover dynamics for a tropical coastal urban city in Kerala, India | |
Salman Aal-shamkhi et al. | Extraction and modeling of urban sprawl development in Karbala City using VHR satellite imagery | |
Min et al. | Spatially downscaling IMERG at daily scale using machine learning approaches over Zhejiang, southeastern China | |
Karamouz et al. | Soil moisture data using citizen science technology cross-validated by satellite data | |
CN114372631A (zh) | 一种基于小样本学习和lstm的缺资料地区径流预测方法 | |
Yuan et al. | A novel method based on deep learning model for national-scale landslide hazard assessment | |
Zhang et al. | Merging multisatellite precipitation products using stacking method and the censored-shifted gamma ensemble model output statistics in china's Beimiaoji basin | |
CN117830381A (zh) | 湖泊面积变化空间分析模型、构建与面积变化预测方法及装置 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |