CN111274545B - 一种栅格尺度基于地形地貌的多模式产流计算方法 - Google Patents
一种栅格尺度基于地形地貌的多模式产流计算方法 Download PDFInfo
- Publication number
- CN111274545B CN111274545B CN202010116808.1A CN202010116808A CN111274545B CN 111274545 B CN111274545 B CN 111274545B CN 202010116808 A CN202010116808 A CN 202010116808A CN 111274545 B CN111274545 B CN 111274545B
- Authority
- CN
- China
- Prior art keywords
- grid
- runoff
- type
- underlying surface
- membership
- 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
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/18—Complex mathematical operations for evaluating statistical data, e.g. average values, frequency distributions, probability functions, regression analysis
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T17/00—Three dimensional [3D] modelling, e.g. data description of 3D objects
- G06T17/05—Geographic models
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Data Mining & Analysis (AREA)
- Theoretical Computer Science (AREA)
- Software Systems (AREA)
- Geometry (AREA)
- Mathematical Physics (AREA)
- Pure & Applied Mathematics (AREA)
- Mathematical Optimization (AREA)
- Mathematical Analysis (AREA)
- Computational Mathematics (AREA)
- Bioinformatics & Cheminformatics (AREA)
- Evolutionary Biology (AREA)
- Bioinformatics & Computational Biology (AREA)
- Life Sciences & Earth Sciences (AREA)
- Operations Research (AREA)
- Probability & Statistics with Applications (AREA)
- Remote Sensing (AREA)
- Algebra (AREA)
- Computer Graphics (AREA)
- Databases & Information Systems (AREA)
- General Engineering & Computer Science (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明公开了一种栅格尺度基于地形地貌的多模式产流计算方法,包括以下步骤:获取流域DEM数据并计算栅格的河面落差、坡度和地形湿度指数值;计算三个指标的累积分布函数,将非水体栅格分为洼地、坡地、高地超渗区和高地下渗区四种下垫面类型,得到四种类型的先验概率;利用最大似然分类法迭代计算各个栅格对每种下垫面类型的隶属度;将每种下垫面类型的产流作为一种产流模式,根据每个栅格的下垫面类型隶属度判别产流的模式和组成,将各产流模式对应的产流量加权叠加得到最终产流量。本发明通过对流域下垫面进行分类并基于栅格对各产流模式的隶属度加权计算得到最终产流量,减少了人为确定分类阈值带来的主观性,促进分布式产流模型的发展。
Description
技术领域
本发明属于水文预报领域,具体涉及一种栅格尺度基于地形地貌的产流模式最大似然判别与产流计算方法。
背景技术
产流是流域水文过程的重要组成部分,因此产流计算也是水文预报中的关键环节。降雨和下垫面因子的空间异质性导致了产流的空间不均匀性,而传统的集总式产流模型无法描述这种空间分布的复杂性。因此,分布式产流模型的构建是当前国内外研究的前沿和热点问题。依据地形、地质、土地利用、土壤类型等下垫面因子的水文和水力特征,分布式产流模型将流域划分为异质性空间计算单元,从而更加客观、合理地刻画产流分布。利用主导产流模式的识别划分流域水文单元是其中的方法之一,但人为确定分类阈值存在较强的主观性与经验性,且下垫面类型与产流模式的一对一构建无法刻画产流的复杂性。为此,需要设计一种新的技术方案给予解决。
发明内容
本发明的目的是提供一种栅格尺度基于地形地貌的产流模式最大似然判别与产流计算方法,利用数理统计原理计算流域栅格对各下垫面类型的隶属度,构建不同下垫面类型下的概念性产流模型,将不同模式下的产流量加权叠加得到各栅格的最终产流量,从而解决现有技术中存在的前述问题。
为了实现上述目的,本发明具体采用以下技术方案:
一种栅格尺度基于地形地貌的多模式产流计算方法,其特征在于,包括以下步骤:
步骤1,获取流域DEM数据并在此基础上计算栅格的河面落差(H)、坡度(Slope)和地形湿度指数(TWI)值;
步骤2,计算三个指标的累积分布函数,将非水体栅格分为洼地、坡地、高地超渗区和高地下渗区四类,得到四种类型的先验概率;
步骤3,利用最大似然分类法迭代得到各个栅格对每种下垫面类型的隶属度;
步骤4,在不同类型下垫面上建立不同的产流模式,将各类型对应的产流量根据下垫面类型隶属度加权叠加得到最终产流量。
所述步骤1中河面落差、坡度、地形湿度指数分别为:
H=d-d0 (1)
式中,H为河面落差,Slope为坡度,TWI为地形湿度指数,d为栅格所在高程,d0为沿着经过该栅格的流线方向所到达的第一个水体栅格的高程,和分别为栅格在水平和垂直方向的高程变化率,a为栅格单宽累积面积。
所述步骤2中计算三个指标的累积分布函数,将非水体栅格分为洼地、坡地、高地超渗区和高地下渗区四类,得到四种类型的先验概率,具体步骤包括:
步骤21,提取除水体外所有栅格的H、Slope和TWI值,分别从大到小排列得到各自的累积分布函数,从而得到栅格各指标的概率分布:
FH(x)=n1/Nt (4)
FSlope(x)=n2/Nt (5)
FTWI(x)=n3/Nt (6)
式中,FH(x),FSlope(x),FTWI(x)分别为该栅格H、Slope和TWI值的累积分布函数,n1,n2,n3分别为H、Slope、TWI值小于等于x的非水体栅格数,Nt为非水体栅格总数。
步骤22,计算非水体栅格对洼地、坡地、高地超渗区和高地下渗区四种类型的隶属度,将栅格类型定义为隶属度最高者对应的类型:
Mw=1-FH (7)
Mh=FH×FSlope (8)
Mpe=[FH×(1-FSlope)]×(1-FTWI) (9)
Mpp=[FH×(1-FSlope)]×FTWI (10)
式中,Mw、Mh、Mpe、Mpp分别为栅格对洼地、坡地、高地超渗区和高地下渗区的隶属度,FH、FSlope、FTWI分别为该栅格H、Slope和TWI值的概率分布;
步骤23,确定洼地、坡地、高地超渗区和高地下渗区四种类型的先验概率并计算每种类型下的栅格各指标的均值和方差(假定呈三维正态分布):
P(wi)=Ni/Nt (11)
式中,P(wi)为第i种类型下垫面的先验概率,wi为第i种下垫面类型,Ni为第i种类型的栅格数量。
所述步骤3中利用最大似然分类法迭代得到各个栅格对每种下垫面类型的隶属度,具体步骤包括:
步骤31,将步骤2中的分类结果作为初始场,用最大似然分类法对非水体栅格进行再分类;
步骤32,将前一次的分类结果作为后一次分类的初始场进行迭代;
步骤33,根据迭代终点,计算所有非水体栅格对四种类型的隶属度并将水体栅格定义为100%隶属于水体分类。
最大似然分类的判别函数及判别规则为:
如果di(X)=maxj=1,2,3,4dj(X),则X∈wi (13)
式中,di(X)为栅格对第i种类型的隶属度,X为栅格的三维指标向量,Di为第i种类型中栅格各指标的协方差矩阵,Ei为第i种类型的三维栅格指标均值向量。
最大似然分类法的迭代终点为:
|Pn(wi)-Pn-1(wi)|<0.001,i=1,2,3,4 (14)
式中,n为迭代次数,Pn(wi)为第n次迭代中第i中类型栅格所占流域的比例。
最终非水体栅格对四种分类的隶属度为:
式中,P(wi|X)为栅格对第i种类型的隶属度,P(X|wi)为第i种类型情况下该栅格出现的概率密度。
所述步骤4中在不同类型下垫面上建立不同的产流模式,将各类型对应的产流量根据隶属度加权叠加得到最终产流量,具体内容包括:水体、洼地、坡地、高地超渗区和高地下渗区五种下垫面类型对应的产流模式分别为:直接径流、饱和地面径流、快速壤中流、超渗地面径流和深层下渗(不产流)。各产流模式的计算公式为:
各非水体栅格最终产流量为:
Rs=P(w1|X)×Rsof+P(w3|X)×Reof (20)
Ri=P(w2|X)×Rssf (21)
Rg=P(w1|X)×Rg1+P(w2|X)×Rg2+P(w3|X)×Rg3 (22)
式中,Rsof、Rssf、Reof、Rg1、Rg2、Rg3分别为饱和地面径流、快速壤中流、超渗地面径流以及洼地、坡地、高地超渗区对应的地下径流,a1为Rsof与Rg1的分割比例,a2为Rssf与Rg2的分割比例,Pe为净雨量,S为蓄水容量,F为下渗能力,Rs、Ri、Rg分别为最终各栅格的地面径流、壤中流、地下径流产流量。
本发明的有益效果:本发明提供一种栅格尺度基于地形地貌的产流模式最大似然判别与产流计算方法,首先获取流域DEM数据并在此基础上计算栅格的H、Slope和TWI值;然后计算三个指标的累积分布函数,将非水体栅格分为洼地、坡地、高地超渗区和高地下渗区四种下垫面类型,得到四种类型的先验概率;之后利用最大似然分类法迭代得到各个栅格对每种类型的隶属度;最后在不同类型下垫面上建立不同的产流模式,将各类型对应的产流量根据隶属度加权叠加得到最终产流量。本发明通过数学方法对流域下垫面进行分类并基于栅格对各产流模式的隶属度加权计算得到最终产流量,减少了人为确定分类阈值带来的主观性,弥补了现阶段一种下垫面对应单一产流模式存在的缺陷,使得产流模型更具有普适性。
附图说明
图1是本发明提供的基于最大似然分类的栅格尺度多模式产流计算方法流程示意图;
图2是具体实施例中流域下垫面的H、Slope、TWI值分布图;
图3是具体实施例中流域下垫面初始分类图;
图4是具体实施例中流域下垫面栅格对每种类型的隶属度分布图。
具体实施方式
下面结合附图和具体实施例对本发明作进一步描述。
应当理解的是,此处所描述的具体实施方式仅仅用以解释本发明,并不用于限定本发明。
如图1所示,本发明实施例提供了一种栅格尺度基于地形地貌的产流模式最大似然判别与产流计算方法,包括以下步骤:
步骤1,获取流域DEM数据并在此基础上计算栅格的河面落差(H)、坡度(Slope)和地形湿度指数(TWI)值,包括:
H=d-d0 (1)
下载安徽省屯溪流域的DEM数据,利用Arcmap在DEM数据基础上计算H值、Slope值和TWI值分布图,如图2所示。
步骤2,计算三个指标的累积分布函数,将非水体栅格分为洼地、坡地、高地超渗区和高地下渗区四种下垫面类型。对于高地区,地形指数小的位置容易超渗,所以将其划分为高地超渗区,其余位置为高地下渗区。
根据划分的四种下垫面类型,得到四种类型的先验概率,具体步骤包括:
步骤21,提取除水体外所有栅格的H、Slope和TWI值,分别从大到小排列得到各自的累积分布函数,从而得到栅格各指标的概率分布,包括:
FH(x)=n1/Nt (4)
FSlope(x)=n2/Nt (5)
FTWI(x)=n3/Nt (6)
式中,FH(x),FSlope(x),FTWI(x)分别为该栅格H、Slope和TWI值的累积分布函数,n1,n2,n3分别为H、Slope、TWI值小于等于x的非水体栅格数,Nt为非水体栅格总数。
步骤22,计算非水体栅格对洼地、坡地、高地超渗区和高地下渗区四种类型的隶属度,将栅格类型定义为隶属度最高者对应的类型,如图3所示:
Mw=1-FH (7)
Mh=FH×FSlope (8)
Mpe=[FH×(1-FSlope)]×(1-FTWI) (9)
Mpp=[FH×(1-FSlope)]×FTWI (10)
式中,Mw、Mh、Mpe、Mpp分别为栅格对洼地、坡地、高地超渗区和高地下渗区的隶属度,FH、FSlope、FTWI分别为该栅格H、Slope和TWI值的概率分布;
步骤23,确定洼地、坡地、高地超渗区和高地下渗区四种类型的先验概率并计算每种类型下的栅格各指标的均值和方差(假定呈三维正态分布):
P(wi)=Ni/Nt (11)
式中,P(wi)为第i种类型下垫面的先验概率,wi为第i种下垫面类型,Ni为第i种类型的栅格数量。
步骤3,利用最大似然分类法迭代得到各个栅格对每种类型的隶属度,具体步骤包括:
步骤31,将步骤2中的分类结果作为初始场,用最大似然分类法对非水体栅格进行再分类;
步骤32,将前一次的分类结果作为后一次分类的初始场进行迭代;
步骤33,根据迭代终点,计算所有非水体栅格对四种下垫面类型的隶属度并将水体栅格定义为100%隶属于水体分类,如图4所示。
最大似然分类的判别函数及判别规则为:
i=1,2,3,4
如果di(X)=maxj=1,2,3,4dj(X),则X∈wi (13)
式中,di(X)为栅格对第i种类型的隶属度,X为栅格的三维指标向量,Di为第i种类型中栅格各指标的协方差矩阵,Ei为第i种类型的三维栅格指标均值向量。
最大似然分类法的迭代终点为:
|Pn(wi)-Pn-1(wi)|<0.001,i=1,2,3,4 (14)
式中,n为迭代次数,Pn(wi)为第n次迭代中第i中类型栅格所占流域的比例。
最终非水体栅格对四种分类的隶属度为:
式中,P(wi|X)为栅格对第i种类型的隶属度,P(X|wi)为第i种类型情况下该栅格出现的概率密度。
步骤4,在不同类型下垫面上建立不同的产流模式,将各类型对应的产流量根据隶属度加权叠加得到最终产流量。
所述步骤4中在不同类型下垫面上建立不同的产流模式,将各类型对应的产流量根据隶属度加权叠加得到最终产流量,具体内容包括:水体、洼地、坡地、高地超渗区和高地下渗区五种下垫面类型对应的产流模式分别为:直接径流、饱和地面径流、快速壤中流、超渗地面径流和深层下渗(不产流)。各产流模式的计算公式为:
各非水体栅格最终产流量为:
Rs=P(w1|X)×Rsof+P(w3|X)×Reof (20)
Ri=P(w2|X)×Rssf (21)
Rg=P(w1|X)×Rg1+P(w2|X)×Rg2+P(w3|X)×Rg3 (22)
式中,Rsof、Rssf、Reof、Rg1、Rg2、Rg3分别为饱和地面径流、快速壤中流、超渗地面径流以及洼地、坡地、高地超渗区对应的地下径流,a1为Rsof与Rg1的分割比例,a2为Rssf与Rg2的分割比例,Pe为净雨量,S为蓄水容量(与各点H值成正比),F为下渗能力,Rs、Ri、Rg分别为最终各栅格的地面径流、壤中流、地下径流产流量。
以上所述仅是本发明的优选实施方式,并非用以限定本发明。本领域技术人员在不脱离本发明构思的前提下做出若干改进和优化,都应当视为本发明保护的范围。
Claims (7)
1.一种栅格尺度基于地形地貌的多模式产流计算方法,其特征在于,包括以下步骤:
步骤1,获取流域DEM数据并在此基础上计算栅格的河面落差、坡度和地形湿度指数值;
步骤2,计算河面落差、坡度和地形湿度指数值三个指标的累积分布函数,将非水体栅格分为洼地、坡地、高地超渗区和高地下渗区四种下垫面类型,得到四种类型的先验概率;
步骤3,利用最大似然分类法迭代得到各个栅格对每种下垫面类型的隶属度;
步骤4,将每种下垫面类型的产流作为一种产流模式并建立相应计算方法,根据每个栅格下垫面类型的隶属度判别产流的模式和组成,将各产流模式对应的产流量根据下垫面类型隶属度加权叠加得到最终产流量;
所述步骤2包括:
步骤21,提取除水体外所有栅格的河面落差、坡度和地形湿度指数值,分别从大到小排列得到各自的累积分布函数,从而得到栅格各指标的概率分布:
FH(x)=n1/Nt (4)
FSlope(x)=n2/Nt (5)
FTWI(x)=n3/Nt (6)
式中,FH(x),FSlope(x),FTWI(x)分别为该栅格河面落差、坡度和地形湿度指数值的累积分布函数,n1,n2,n3分别为河面落差、坡度、地形湿度指数值小于等于x的非水体栅格数,Nt为非水体栅格总数;
步骤22,计算非水体栅格对洼地、坡地、高地超渗区和高地下渗区四种类型的隶属度,将栅格类型定义为隶属度最高者对应的类型:
Mw=1-FH (7)
Mh=FH×FSlope (8)
Mpe=[FH×(1-FSlope)]×(1-FTWI) (9)
Mpp=[FH×(1-FSlope)]×FTWI (10)
式中,Mw、Mh、Mpe、Mpp分别为栅格对洼地、坡地、高地超渗区和高地下渗区的隶属度,FH、FSlope、FTWI分别为栅格的河面落差、坡度和地形湿度指数值的概率分布;
步骤23,确定洼地、坡地、高地超渗区和高地下渗区四种类型的先验概率并计算每种类型下的栅格各指标的均值和方差:
P(wi)=Ni/Nt (11)
式中,P(wi)为第i种类型下垫面的先验概率,wi为第i种下垫面类型,Ni为第i种类型的栅格数量。
2.根据权利要求1所述的多模式产流计算方法,其特征在于,所述步骤3包括:
步骤31,将步骤2中的分类结果作为初始场,用最大似然分类法对非水体栅格进行再分类;
步骤32,将前一次的分类结果作为后一次分类的初始场进行迭代;
步骤33,根据迭代终点,计算最终非水体栅格对四种下垫面类型的隶属度并将水体栅格定义为100%隶属于水体分类。
4.根据权利要求3所述的多模式产流计算方法,其特征在于,最大似然分类法的迭代终点为:
|Pn(wi)-Pn-1(wi)|<0.001, (14)
i=1,2,3,4
式中,n为迭代次数,Pn(wi)为第n次迭代中第i种类型栅格所占流域的比例。
6.根据权利要求5所述的多模式产流计算方法,其特征在于,所述步骤4包括:水体、洼地、坡地、高地超渗区和高地下渗区五种下垫面类型对应的产流模式分别为:直接径流、饱和地面径流、快速壤中流、超渗地面径流和深层下渗。
7.根据权利要求6所述的多模式产流计算方法,其特征在于,各产流模式的计算公式为:
各非水体栅格最终产流量为:
Rs=P(w1|X)×Rsof+P(w3|X)×Reof (20)
Ri=P(w2|X)×Rssf (21)
Rg=P(w1|X)×Rg1+P(w2|X)×Rg2+P(w3|X)×Rg3 (22)
式中,Rsof、Rssf、Reof、Rg1、Rg2、Rg3分别为饱和地面径流、快速壤中流、超渗地面径流以及洼地、坡地、高地超渗区对应的地下径流,a1为Rsof与Rg1的分割比例,a2为Rssf与Rg2的分割比例,Pe为净雨量,S为蓄水容量,F为下渗能力,Rs、Ri、Rg分别为最终各栅格的地面径流、壤中流、地下径流产流量。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010116808.1A CN111274545B (zh) | 2020-02-25 | 2020-02-25 | 一种栅格尺度基于地形地貌的多模式产流计算方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010116808.1A CN111274545B (zh) | 2020-02-25 | 2020-02-25 | 一种栅格尺度基于地形地貌的多模式产流计算方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN111274545A CN111274545A (zh) | 2020-06-12 |
CN111274545B true CN111274545B (zh) | 2020-10-16 |
Family
ID=71002302
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202010116808.1A Active CN111274545B (zh) | 2020-02-25 | 2020-02-25 | 一种栅格尺度基于地形地貌的多模式产流计算方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN111274545B (zh) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112926418B (zh) * | 2021-02-08 | 2022-07-01 | 河海大学 | 一种利用高精度地形信息消除云层干扰的水体优化提取方法 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102508961A (zh) * | 2010-12-16 | 2012-06-20 | 南京大学 | 一种高分辨率的全分布式水文模型topx的设计方法 |
CN103236086A (zh) * | 2013-04-24 | 2013-08-07 | 武汉大学 | 一种顾及地表水文上下文的多尺度dem建模方法 |
CN106446394A (zh) * | 2016-09-19 | 2017-02-22 | 河海大学 | 通过地形指数提取流域自由水蓄水容量空间分布的方法 |
CN110570517A (zh) * | 2019-08-07 | 2019-12-13 | 河海大学 | 一种基于下垫面特征的重配置产流模拟方法 |
CN110598553A (zh) * | 2019-08-09 | 2019-12-20 | 中国科学院南京地理与湖泊研究所 | 基于遥感影像与地形数据的原真地貌采矿破坏区探测方法 |
-
2020
- 2020-02-25 CN CN202010116808.1A patent/CN111274545B/zh active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102508961A (zh) * | 2010-12-16 | 2012-06-20 | 南京大学 | 一种高分辨率的全分布式水文模型topx的设计方法 |
CN103236086A (zh) * | 2013-04-24 | 2013-08-07 | 武汉大学 | 一种顾及地表水文上下文的多尺度dem建模方法 |
CN106446394A (zh) * | 2016-09-19 | 2017-02-22 | 河海大学 | 通过地形指数提取流域自由水蓄水容量空间分布的方法 |
CN110570517A (zh) * | 2019-08-07 | 2019-12-13 | 河海大学 | 一种基于下垫面特征的重配置产流模拟方法 |
CN110598553A (zh) * | 2019-08-09 | 2019-12-20 | 中国科学院南京地理与湖泊研究所 | 基于遥感影像与地形数据的原真地貌采矿破坏区探测方法 |
Non-Patent Citations (3)
Title |
---|
DEM栅格单元异质性对地形湿度指数提取的影响分析;马建超 等;《地球信息科学学报》;20110415;第13卷(第2期);全文 * |
Hydrological regionalisation based on available hydrological information for runoff prediction at catchment scale;Li Qiaoling 等;《Proceedings of the International Association of Hydrological Sciences (IAHS)》;COPERNICUS GESELLSCHAFT MBH, BAHNHOFSALLE 1E, GOTTINGEN, 37081, GERMANY;20180615;第379卷;全文 * |
Remote sensing and GIS based watershed prioritization;Manish Pandey 等;《2017 IEEE International Geoscience and Remote Sensing Symposium (IGARSS)》;IEEE;20170728;全文 * |
Also Published As
Publication number | Publication date |
---|---|
CN111274545A (zh) | 2020-06-12 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106780089B (zh) | 基于神经网络元胞自动机模型的永久性基本农田划定方法 | |
CN110458048A (zh) | 顾及城镇格局特征的人口分布时空演变与认知 | |
CN103236086B (zh) | 一种顾及地表水文上下文的多尺度dem建模方法 | |
CN108536908B (zh) | 基于非点源氮磷流失风险对流域水环境安全评估的方法 | |
CN105118090B (zh) | 一种自适应复杂地形结构的点云滤波方法 | |
CN102184423B (zh) | 一种全自动的区域不透水面遥感信息精确提取方法 | |
CN106056591A (zh) | 一种融合光谱图像和激光雷达数据进行城市密度估计方法 | |
CN103150614B (zh) | 一种土地利用空间自动配置方法 | |
CN113887635B (zh) | 一种流域相似性分类方法及分类装置 | |
CN117332542B (zh) | 一种多尺度自适应选择的城市洪涝建模模拟方法 | |
CN106844979B (zh) | 一种城市虚拟公园交互式组合设计系统及方法 | |
CN116050163B (zh) | 一种基于气象站的生态系统碳水通量计算方法及系统 | |
CN116562583A (zh) | 一种多维度水资源供需预测方法及系统 | |
CN116205136A (zh) | 基于径流滞后信息的大尺度流域深度学习洪水预报方法 | |
CN102880754A (zh) | 基于遗传算法的土地利用分形维数作用尺度的识别方法 | |
CN111274545B (zh) | 一种栅格尺度基于地形地貌的多模式产流计算方法 | |
CN117892167B (zh) | 一种替代数值模拟的城市内涝代理模型构建方法 | |
CN106096129A (zh) | 一种基于山地汇水计算的山脚水面规模分析方法 | |
CN116977666A (zh) | 一种城市道路点云数据的绿度空间的评价方法及可读存储介质 | |
CN114611290B (zh) | 一种基于量变参数水文不确定性处理器的场次洪水水文模型实时预报方法 | |
Huang et al. | Mapping gully affected areas by using Sentinel 2 imagery and digital elevation model based on the Google Earth Engine | |
KR102408845B1 (ko) | 산림물지도를 활용한 수원함양기능평가방법 및 시스템 | |
CN112711917B (zh) | 一种基于人脸识别算法的实时洪水预报智能方法 | |
KR102408827B1 (ko) | 컴퓨터 프로그램을 이용한 산림물지도 제작방법 | |
CN112507635B (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 |