CN108010103B - 复杂河道地形快速精细生成方法 - Google Patents
复杂河道地形快速精细生成方法 Download PDFInfo
- Publication number
- CN108010103B CN108010103B CN201711195302.9A CN201711195302A CN108010103B CN 108010103 B CN108010103 B CN 108010103B CN 201711195302 A CN201711195302 A CN 201711195302A CN 108010103 B CN108010103 B CN 108010103B
- Authority
- CN
- China
- Prior art keywords
- section
- elevation
- interpolation
- point
- boundary
- 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
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T11/00—2D [Two Dimensional] image generation
- G06T11/20—Drawing from basic elements, e.g. lines or circles
- G06T11/206—Drawing of charts or graphs
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F16/00—Information retrieval; Database structures therefor; File system structures therefor
- G06F16/20—Information retrieval; Database structures therefor; File system structures therefor of structured data, e.g. relational data
- G06F16/29—Geographical information databases
Landscapes
- Engineering & Computer Science (AREA)
- Theoretical Computer Science (AREA)
- Databases & Information Systems (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Remote Sensing (AREA)
- Data Mining & Analysis (AREA)
- General Engineering & Computer Science (AREA)
- Instructional Devices (AREA)
- Processing Or Creating Images (AREA)
Abstract
本发明公开了一种河道地形快速精细生成方法。通过对复杂地形边界点的空间关系进行统一分块、归类、编号,将零散边界点拟合天然河道走向并生成新河道边界,解决河道大量边界点无序问题;通过有限河道监测断面地形,基于权重插值方法,实现河道地形快速生成;通过河道插值生成地形,利用待增补断面与河道已知采样高程断面的距离和比降变化,实现河道断面的精确增补。本发明根据河道零散边界点和地形高程采样数据无法满足河道水动力计算精度要求的特点,采用了拟合河道边界、河道地形插值和河道采样高程断面增补的方式,充分考虑到了天然河道走势和地势变化特点,有效减少了地形插值误差。
Description
技术领域
本发明属于地理信息***空间技术领域,涉及河道地形的拟合,特别涉及一种复杂河道地形快速精细生成方法。
背景技术
数字高程模型(DEM)是描述地表起伏形态特征的空间数据模型,是由地表规则网格点高程值构成的矩阵,形成栅格结构数据集。在水动力及水环境数值模拟中,计算网格上地形值能否反映真实地形非常重要,会直接影响计算精度。DEM是实际地形表面的再现,其对地形表达的可信程度,很大程度上取决于高程采样点的分布与密度。然而,高程采样点无法对地表面所有点进行观测,只能获取一定数量的高程采样断面数据,这些采样断面反映了地形的局部或者部分特征。尤其是对河道来讲,在许多实际情况下没有整个河道的连续地形数据,只有一定间距下的断面地形数据,数值模拟计算时需要通过插值方法由断面数据得到网格地形值,尤其河道转弯或河流交汇等起伏变化很大处的地形。如何快速生成河道地形,并满足高精度和高效率的计算需求,是当前环境管理与模拟分析所面临的重要问题,也是地形插值精度所面临的重要技术问题。
现有关于空间地形高程采样点加密的插值方法,主要集中在空间内插方法。空间内插是根据分布在内插点周围的样本点求出其值。现有内插方法包括最近邻法、反距离加权法、多项式法等,计算相对简单,但其忽略了样本点间存在的空间联系,插值结果往往受样本点的影响较大。克里金插值法以空间自相关为基础,利用原始样本点数据和半方差函数的结构性,对区域化变量的待内插点进行无偏最优估值。这类插值方法插值精度相关较高,但计算步骤较为繁琐,计算速度较慢。
空间内插方法给地形插值提供了一定的参考基础,高程内插应用上其通常考虑高程数值插值精度,忽略了插值河道边界以及监测断面不足的影响。目前,缺少一种在复杂河道快速精细的地形生成方法,尤其是在加密过程中,在确保高程数值精度和计算效率的情况下,如何依照河道走势和地势变化是获取高精度的插值地形的关键。本发明即针对该问题,实现了复杂河道地形的精细快速生成,能满足数值模拟计算精度和效率需求。
发明内容
本发明针对上述问题,面向地形高程采样点数据,设计了一种复杂地形快速精细生成方法。该方法涵盖了河道边界生成、区域高程采样点加密以及增补监测断面三大关键技术环节,能够有效解决高精度地形建模中高程采样点不足而导致地形表达失真的问题。
为实现上述发明目的,本发明采用的技术方案如下:
一种复杂地形快速精细生成方法,包括如下步骤:
第一步,河道边界趋势拟合:根据河道边界散点分布,按照走向趋势变化进行分块编码;对分块编码散点,进行趋势拟合,拟合河道边界的整体走向;整体走向判定无误,则执行第二步;若判定边界走向有误,则继续执行第一步;
第二步,对区域高程采样点进行加密:首先,根据已知监测断面高程采样点及河道边界要素,采用权重法,计算出加密点的空间位置信息;其次,依据地势变化,对加密点高程进行插值计算;最后,依据空间加密点的密集程度要求,进行进一步细化或粗化;
第三步,采样高程点断面增补:当已知监测断面不足时,易造成河道弯道处插值地形失真;读取弯道处待增补断面地形点的空间位置,根据待增补断面与河道上游已知监测断面的距离以及比降变化,得到待增补断面上插值点的高程,完成采样高程点断面增补,再按照第二步,完成整个河道插值地形加密处理。
上述第一步包括以下步骤:
步骤1,流域水系结构和边界的识别:将河道CAD图层导入GIS中,将线转换成线图层,或者基于高分辨率遥感影像数据,识别水系结构层中的流域边界范围层,导入GIS中转化为线图层;
步骤2,对线图层进行处理,包括河流边界的处理及坐标系转化:首先,对河流边界线图层进行编辑,对于少数复杂小枝杈存在的地方,可进行边界线平滑化处理,便于后续插值;其次,在GIS中,完成WGS 1984坐标系到WGS 1984UTM Zone 49N坐标系的坐标转化;最后,将处理后的线图层在GIS中进行线转点,分成河道左岸与右岸依次处理;
步骤3,河道边界散点的分块编码:按照河道散点坐标X和Y分布的走向,将其分为四种类型,即:X从小到大且Y从小到大,Y从小到大且X从大到小,X从大到小且Y从大到小,Y从大到小且X从小到大,四种情况分别赋予其属性值ND为①、②、③、④;通过勾画矩形框,将散点进行分块,河道分块后,要读取每块上面的最小的X1,Y1和最大的X2,Y2,即每个矩形块的左下角点和右上角点的两个点;按照从河道上游往下游依次分块,编号,构建河道边界散点关系表,该表中包含块的编号,块的范围,块的方向,块所属的左右岸边界线,以及左右岸边界线的分块块数;
步骤4,通过边界散点数据,判断边界散点所属的分块:基于每块的分块范围,搜索全局散点,判断边界散点的所属分块;基于每块内的散点坐标,以及块的方向属性,利用多项式曲线拟合,在误差最小情况下,拟合出当前分块的边界点趋势,依次类推,拟合出左右岸各分块的边界点趋势;
步骤5,分块的边界趋势拟合结束后,需将全部拟合河道边界点导入EXCEL中,按照拟合趋势连线,若无交叉或断线情况,则拟合无误,反之,需校核分块的拟合方向,重新进行拟合。
上述第二步包括以下步骤:
步骤1,以采样高程点为初始高程断面,基于生成河道拟合边界数据,通过权重插值方法生成加密点平面位置(X,Y);
步骤2,加密点高程插值计算:基于初始高程点断面及加密高程点平面位置,按照对应比例进行计算,计算包括:①以初始高程断面上采样点的起点距,计算每个采样点在整个断面上的相对位置比例;②通过上下游两个断面上采样点位置的权重和高程,按照X方向和Y方向设定的插值距离,将上下游两个断面重新插值,并获取高程;③在上下游两个断面之间,按照X和Y方向设定的插值距离,进行河段上空间插值,获取每点的插值坐标(X,Y);④根据插值坐标,计算每个插值断面上各插值点在断面上的相对位置比例;⑤搜索初始上下游高程断面采样点与插值断面插值点的权重相等的两点,根据两点间的距离以及河道的坡度变化大小,计算插值点位置的高程;
步骤3,完成细化或粗化计算:调整X和Y方向的插值距离,重复步骤1和步骤2,直到达到精确反应河道地形为止。
上述第三步包括以下步骤:
步骤1,河道区域弯道处采样高程断面的增补,根据弯道宽度在弯道横断面上选取两点,该两点所在断面即为待增补断面,读取待增补断面上插值点的空间位置,并运用GIS计算插值点的坐标;
步骤2,根据待增补断面的坐标,找到在已知插值完成断面中的权重位置,确定该权重处插值的原始河道地形断面上各插值点的起点距,从而得到待增补断面上各点的起点距;
步骤3,以待增补断面与上游已知采样高程断面的相对距离和比降变化,得到待增补断面上点的高程,其可通过下式求取:
Hi=hi+l*k
式中,Hi为待增补断面上起点距为i处的高程值;hi为上游采样高程断面上起点距为i处的高程值;l为待增补断面与上游采样高程断面的相对距离;k为待增补断面与上游采样高程断面之间的河道比降;
步骤4,完成待增补断面上的起点距和高程值,将待增补断面***到原始采样高程断面,合并成新的采样高程断面数据集,再按照第二步,完成整个区域干支流高程采样点的加密处理。
本发明与现有技术相比,具有以下优点和效果:
(1)本发明提出的河道边界拟合方法,兼顾了区域地形要素及线性地形要素,在边界点零散无序的情况下,将边界点重新分块定义并进行边界拟合,最大程度挖掘和保持地表地形特征,该方法为地形数据采集规范完善提供了技术支撑。
(2)本发明构建并实现的地形高程采样点加密方法,在顾及地表形态特征的基础上实现了稀疏样点向密集样点的转变,实现了高程采样方式的基本覆盖。该方法提高了地形建模数据精度,特别是复杂河道地形建模的数据精度。
(3)本发明构建了增补监测断面方法,在监测断面不足情况下,最大化利用已有监测断面信息和插值地形,结合地势变化趋势,最大限度精度插值支流地形。该方法极大地提高了地形插值精度,解决了采样高程断面缺乏问题。
附图说明
图1为河道地形插值技术框架;
图2为河道边界散点分块编码示意图;
图3为分块编码方向属性定义示意图;
图4为生成拟合河道边界示意图;
图5为地形插值加密过程示意图,(a)已知监测断面,(b)插补河道地形;
图6为河道地形插值2倍细化过程示意图,(a)原始插值地形,(b)2倍细化插值地形;
图7大宁河河道地形插值2倍细化,(a)大宁河河道地形插值细化前,(b)大宁河河道地形插值细化后;
图8大宁河增补采样高程断面,(a)大宁河DN14和DN15断面间插值地形,(b)大宁河DN14和DN15断面间增补断面后插值地形。
具体实施方式
下面结合附图对本发明作进一步详细描述:
(一)河道边界散点分块编码
①将大宁河河道的CAD地形导入GIS处理:需要对线图层进行处理,包括河流边界的处理及其与干流连接吻合度的处理。首先,在GIS中对河流边界线图层进行编辑,对于少数复杂小枝杈存在的地方,可进行边界线平滑化处理,便于后续插值。其次,在GIS左侧窗口中,“Layers”上右键,选择“Properties-Coordinate system-Predefined-Projectedcoordinatesystems-UTM-WGS 1984-Northern Hemisphere-WGS 1984UTM Zone 49N”,完成WGS1984坐标系到WGS 1984UTM Zone 49N坐标系的坐标转化。
②线转点:在GIS工具箱中“Data management Tools-Features-Feature topoint”,将线图层转成点图层,可分为左岸线和右岸线两个图层。在点图层属性表里新加两列(双精度),命名X和Y,重新计算坐标,完成后导出属性表为.dbf格式文件,即得到河流边界点坐标。
③分块编码:如图2所示,按照河道左右岸边界散点的坐标(X,Y),分为X从小到大且Y从小到大,Y从小到大且X从大到小,X从大到小且Y从大到小,Y从大到小且X从小到大四种形式,先对河道边界散点从上游往下游依次进行分块,并进行编号,先右岸再左岸,边界散点分块编码如图2所示。河道分块后,要读取每块上面的最小的X1,Y1和最大的X2,Y2。在GIS中就是每块的左下角点和右上角点的两个点。按照从大宁河从河道上游到下游依次读取每块的左下角点和右上角点两个点的坐标,同时判断每块的方向,分块编码方向属性定义如图3所示。
④分块关系表构建:基于上述三步的编码,构建如图3所示的河道边界散点关系表,如表1所示。表1中,分块方向:①、②、③、④分别表示X从小到大且Y从小到大,Y从小到大且X从大到小,X从大到小且Y从大到小,Y从大到小且X从小到大;左右岸边界线:1为左边界线,2为右边界线;
表1河道边界散点关系表
⑤拟合河道边界。在河道拟合完毕后,还有一点需要特别注意,将全部拟合河道边界点放入EXCEL中,进行点连成线,确保新生成的边界点是按照顺序的即无交叉或断线情况(图4),否则就会插值地形有很大影响。因此,在生成完河道边界点后,一定要放入EXCEL中校核一遍。
(二)河道地形插值
①以采样高程点为初始高程断面,基于生成河道拟合边界数据,通过权重插值方法生成加密点平面位置(X,Y)。如图5所示,在已知两个监测断面DNn与DNn+1之间进行河道地形插补(图5(a)),得到两个断面间的插值地形点高程(图5(b))。
1)以初始高程断面上采样点的起点距,计算每个采样点在整个断面上的相对位置比例;
2)通过上下游两个断面上采样点位置的权重和高程,按照X方向和Y方向设定的插值距离,将上下游两个断面重新插值,并获取高程;
3)在上下游两个断面之间,按照X和Y方向设定的插值距离,进行河段上空间插值,获取每点的插值坐标(X,Y);
4)根据插值坐标,计算每个插值断面上各插值点在断面上的相对位置比例;
5)搜索初始上下游高程断面采样点与插值断面插值点的权重相等的两点,根据两点间的距离以及河道的坡度变化大小,计算插值点位置的高程。
②完成细化或粗化计算:调整X和Y方向的插值距离,重复步骤1和步骤2,直到达到精确反应河道地形为止(图6)。如图7所示,在三峡利于大宁河局部河段上,通过将X和Y方向的插值距离有原来的40m、40m(图7(a))调整为20m、20m(图7(b)),更能反应出河道的真实地形变化。
(三)采样高程断面增补
①大宁河一共22个采样高程断面,一共21段,在第21段上,原始的采样高程断面为DN14和DN15,插值出的河道地形如图8(a)所示,需在弯道处进行增补采样高程断面,根据弯道宽度在弯道横断面上选取两点,该两点所在断面即为待增补断面DN14+1,读取待增补断面上上插值点的空间位置,并运用GIS计算其坐标(图8(b));
②根据待增补断面DN14+1的坐标,找到在已知DN14与DN15两个断面间完成的插值地形中沿河道方向的权重,从而确定该权重处的垂直河道插值断面上各插值点的起点距值,得到待增补断面DN14+1上各插值点的起点距;
③以待增补断面DN14+1与上游已知采样高程断面DN15的相对距离500m和比降变化0.0013,得到待增补断面DN14+1点的高程,其可通过下式求取:
Hi=hi+l*k (1)
式中,Hi为待增补断面上起点距为i处的高程值;hi为上游采样高程断面上起点距为i处的高程值;l为待增补断面与上游采样高程断面的相对距离;k为待增补断面与上游采样高程断面之间的河道比降。
④完成待增补断面上的起点距和高程值,将待增补断面***到原始采样高程断面,合并成新的采样高程断面数据集,一共23个断面,再按照第二步,完成整个区域干支流高程采样点的加密处理,如图8(b)所示。
Claims (1)
1.一种复杂地形快速精细生成方法,其特征在于,包括如下步骤:
第一步,河道边界趋势拟合:根据河道边界散点分布,按照走向趋势变化进行分块编码;对分块编码散点,进行趋势拟合,拟合河道边界的整体走向;整体走向判定无误,则执行第二步;若判定边界走向有误,则继续执行第一步;
第二步,对区域高程采样点进行加密:首先,根据已知监测断面高程采样点及河道边界要素,采用权重法,计算出加密点的空间位置信息;其次,依据地势变化,对加密点高程进行插值计算;最后,依据空间加密点的密集程度要求,进行进一步细化或粗化;
第三步,采样高程点断面增补:当已知监测断面不足时,易造成河道弯道处插值地形失真;读取弯道处待增补断面地形点的空间位置,根据待增补断面与河道上游已知监测断面的距离以及比降变化,得到待增补断面上插值点的高程,完成采样高程点断面增补,再按照第二步,完成整个河道插值地形加密处理;
所述的第一步包括以下步骤:
步骤1,流域水系结构和边界的识别:将河道CAD图层导入GIS中,将线转换成线图层,或者基于高分辨率遥感影像数据,识别水系结构层中的流域边界范围层,导入GIS中转化为线图层;
步骤2,对线图层进行处理,包括河流边界的处理及坐标系转化:首先,对河流边界线图层进行编辑,对于少数复杂小枝杈存在的地方,可进行边界线平滑化处理,便于后续插值;其次,在GIS中,完成WGS 1984坐标系到WGS 1984 UTM Zone 49N坐标系的坐标转化;最后,将处理后的线图层在GIS中进行线转点,分成河道左岸与右岸依次处理;
步骤3,河道边界散点的分块编码:按照河道散点坐标X和Y分布的走向,将其分为四种类型,即:X从小到大且Y从小到大,Y从小到大且X从大到小,X从大到小且Y从大到小,Y从大到小且X从小到大,四种情况分别赋予其属性值ND为①、②、③、④;通过勾画矩形框,将散点进行分块,河道分块后,要读取每块上面的最小的X1,Y1和最大的X2,Y2,即每个矩形块的左下角点和右上角点的两个点;按照从河道上游往下游依次分块,编号,构建河道边界散点关系表,该表中包含块的编号,块的范围,块的方向,块所属的左右岸边界线,以及左右岸边界线的分块块数;
步骤4,通过边界散点数据,判断边界散点所属的分块:基于每块的分块范围,搜索全局散点,判断边界散点的所属分块;基于每块内的散点坐标,以及块的方向属性,利用多项式曲线拟合,在误差最小情况下,拟合出当前分块的边界点趋势,依次类推,拟合出左右岸各分块的边界点趋势;
步骤5,分块的边界趋势拟合结束后,需将全部拟合河道边界点导入EXCEL中,按照拟合趋势连线,若无交叉或断线情况,则拟合无误,反之,需校核分块的拟合方向,重新进行拟合;
所述的第二步包括以下步骤:
步骤1,以采样高程点为初始高程断面,基于生成河道拟合边界数据,通过权重插值方法生成加密点平面位置(X,Y);
步骤2,加密点高程插值计算:基于初始高程点断面及加密高程点平面位置,按照对应比例进行计算,计算包括:①以初始高程断面上采样点的起点距,计算每个采样点在整个断面上的相对位置比例;②通过上下游两个断面上采样点位置的权重和高程,按照X方向和Y方向设定的插值距离,将上下游两个断面重新插值,并获取高程;③在上下游两个断面之间,按照X和Y方向设定的插值距离,进行河段上空间插值,获取每点的插值坐标(X,Y);④根据插值坐标,计算每个插值断面上各插值点在断面上的相对位置比例;⑤搜索初始上下游高程断面采样点与插值断面插值点的权重相等的两点,根据两点间的距离以及河道的坡度变化大小,计算插值点位置的高程;
步骤3,完成细化或粗化计算:调整X和Y方向的插值距离,重复步骤1和步骤2,直到达到精确反应河道地形为止;
所述的第三步包括以下步骤:
步骤1,河道区域弯道处采样高程断面的增补,根据弯道宽度在弯道横断面上选取两点,该两点所在断面即为待增补断面,读取待增补断面上插值点的空间位置,并运用GIS计算插值点的坐标;
步骤2,根据待增补断面的坐标,找到在已知插值完成断面中的权重位置,确定该权重处插值的原始河道地形断面上各插值点的起点距,从而得到待增补断面上各点的起点距;
步骤3,以待增补断面与上游已知采样高程断面的相对距离和比降变化,得到待增补断面上点的高程,其可通过下式求取:
Hi=hi+l*k
式中,Hi为待增补断面上起点距为i处的高程值;hi为上游采样高程断面上起点距为i处的高程值;l为待增补断面与上游采样高程断面的相对距离;k为待增补断面与上游采样高程断面之间的河道比降;
步骤4,完成待增补断面上的起点距和高程值,将待增补断面***到原始采样高程断面,合并成新的采样高程断面数据集,再按照第二步,完成整个区域干支流高程采样点的加密处理。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201711195302.9A CN108010103B (zh) | 2017-11-24 | 2017-11-24 | 复杂河道地形快速精细生成方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201711195302.9A CN108010103B (zh) | 2017-11-24 | 2017-11-24 | 复杂河道地形快速精细生成方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN108010103A CN108010103A (zh) | 2018-05-08 |
CN108010103B true CN108010103B (zh) | 2019-12-24 |
Family
ID=62053631
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201711195302.9A Active CN108010103B (zh) | 2017-11-24 | 2017-11-24 | 复杂河道地形快速精细生成方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN108010103B (zh) |
Families Citing this family (21)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108961403A (zh) * | 2018-06-22 | 2018-12-07 | 武汉大学 | 一种运用开源街道地图提取主干道的方法 |
CN108829997B (zh) * | 2018-06-25 | 2022-04-29 | 中国电建集团中南勘测设计研究院有限公司 | 一种批量提取河流大断面数据的方法及装置 |
CN108876872B (zh) * | 2018-06-25 | 2022-11-11 | 中国电建集团中南勘测设计研究院有限公司 | 一种基于断面要素的断面插值方法及装置 |
CN108986222B (zh) * | 2018-08-15 | 2022-10-14 | 四川大学 | 无汊河道数字地形生成方法 |
CN109031300B (zh) * | 2018-09-03 | 2020-06-05 | 中科卫星应用德清研究院 | 合成孔径雷达监测危岩体变形方法及*** |
CN110083673B (zh) * | 2019-04-30 | 2020-12-15 | 河海大学 | 一种基于高程散点的河道或堤防智能搜索方法 |
CN110120058B (zh) * | 2019-04-30 | 2022-08-12 | 河海大学 | 一种高程散点生成紧致外边界的方法 |
CN110232739A (zh) * | 2019-06-21 | 2019-09-13 | 国核电力规划设计研究院重庆有限公司 | 一种cad地形图数据转换为电力线路断面的方法 |
CN110309166B (zh) * | 2019-07-10 | 2021-10-01 | 中国电子科技集团公司第二十八研究所 | 一种可溯源的地理高程数据补齐方法 |
CN110751726B (zh) * | 2019-10-08 | 2023-04-18 | 江苏省水利科学研究院 | 一种河道工程质量检测方法 |
CN111145347B (zh) * | 2019-11-18 | 2023-04-25 | 长江水利委员会长江科学院 | 地形断面的数字高程数据的校正方法、装置及校正设备 |
CN110929415B (zh) * | 2019-12-06 | 2023-08-29 | 江苏师范大学 | 一种针对高原地区u型河谷水下未知地形精确模拟方法 |
CN111340950B (zh) * | 2020-02-18 | 2023-09-15 | 云南师范大学 | 数字高程模型重建方法及装置 |
CN111444561A (zh) * | 2020-03-04 | 2020-07-24 | 中电建生态环境集团有限公司 | 一种河道地形模型建立方法及装置 |
CN111681316B (zh) * | 2020-06-24 | 2023-05-30 | 河南省水利勘测设计研究有限公司 | 一种高精度河道地形插值方法 |
CN113656852B (zh) * | 2020-07-17 | 2023-07-04 | 长江水利委员会长江科学院 | 一种精细化河道地形快速生成方法 |
CN112182814A (zh) * | 2020-09-11 | 2021-01-05 | 河海大学 | 一种基于稀疏断面点据的河道水下地形建模方法 |
CN113139337B (zh) * | 2021-04-14 | 2023-09-22 | 扬州大学 | 一种用于湖泊地形模拟的分区插值处理方法与装置 |
CN113327323B (zh) * | 2021-06-09 | 2022-11-11 | 四川大学 | 基于散点数据的水体环境地形构建方法 |
CN115631218B (zh) * | 2022-12-21 | 2023-04-11 | 广东中浦科技有限公司 | 一种基于普通克里金插值的臭氧传输变化分析方法 |
CN115950679B (zh) * | 2023-03-11 | 2023-05-12 | 云南省生态环境科学研究院 | 一种河道底泥取样方法 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103425856A (zh) * | 2012-05-24 | 2013-12-04 | 泰瑞数创科技(北京)有限公司 | 一种依据河道测量断面数据自动生成河道地形的方法 |
CN103886187A (zh) * | 2014-03-06 | 2014-06-25 | 清华大学 | 一种基于数据同化的河道水沙实时预测方法 |
CN104598701A (zh) * | 2015-03-04 | 2015-05-06 | 广东工业大学 | 一种三维场景分段迭代的河网污染仿真推演方法 |
CN105631168A (zh) * | 2016-03-25 | 2016-06-01 | 中国水利水电科学研究院 | 一种实时高效的流域洪水演进可视化仿真方法 |
CN105844709A (zh) * | 2016-03-25 | 2016-08-10 | 中国水利水电科学研究院 | 复杂河道地形流域洪水演进虚拟仿真的淹没线追踪方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US7412362B2 (en) * | 2005-05-18 | 2008-08-12 | Microsoft Corporation | River modeling |
-
2017
- 2017-11-24 CN CN201711195302.9A patent/CN108010103B/zh active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103425856A (zh) * | 2012-05-24 | 2013-12-04 | 泰瑞数创科技(北京)有限公司 | 一种依据河道测量断面数据自动生成河道地形的方法 |
CN103886187A (zh) * | 2014-03-06 | 2014-06-25 | 清华大学 | 一种基于数据同化的河道水沙实时预测方法 |
CN104598701A (zh) * | 2015-03-04 | 2015-05-06 | 广东工业大学 | 一种三维场景分段迭代的河网污染仿真推演方法 |
CN105631168A (zh) * | 2016-03-25 | 2016-06-01 | 中国水利水电科学研究院 | 一种实时高效的流域洪水演进可视化仿真方法 |
CN105844709A (zh) * | 2016-03-25 | 2016-08-10 | 中国水利水电科学研究院 | 复杂河道地形流域洪水演进虚拟仿真的淹没线追踪方法 |
Non-Patent Citations (2)
Title |
---|
基于较高精度土壤库的三峡水库汇水区径流模拟;赵佳婧等;《水土保持研究》;20170630;第24卷(第3期);第45-52、58页 * |
太湖流域典型区域水生态环境***模拟平台研究;彭虹等;《中国水利水电科学研究院学报》;20110331;第9卷(第1期);第47-52页 * |
Also Published As
Publication number | Publication date |
---|---|
CN108010103A (zh) | 2018-05-08 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN108010103B (zh) | 复杂河道地形快速精细生成方法 | |
Mitasova et al. | Modelling spatially and temporally distributed phenomena: new methods and tools for GRASS GIS | |
CN103236086B (zh) | 一种顾及地表水文上下文的多尺度dem建模方法 | |
CN106504326B (zh) | 顾及形态精度的地形高程采样点加密方法 | |
CN108986222B (zh) | 无汊河道数字地形生成方法 | |
CN104574512A (zh) | 一种顾及地形语义信息的多尺度dem构建方法 | |
CN111681316B (zh) | 一种高精度河道地形插值方法 | |
CN102819568A (zh) | 基于地形采样点位置的海底地形数据建立方法 | |
CN111383336A (zh) | 三维地质模型构建方法 | |
CN104573234B (zh) | 一种基于ArcGIS的流域库容密度空间化方法 | |
Kenny et al. | Routing overland flow through sinks and flats in interpolated raster terrain surfaces | |
CN103927418A (zh) | 基于dem的城市道路渠网化排水通道制作方法 | |
CN114648617A (zh) | 一种基于数字高程模型dem的水系提取方法 | |
CN113235655A (zh) | 海底沉管隧道基槽分土质工程土方量计算方法 | |
CN110766792A (zh) | 基于ArcGIS水文分析工具网格演算次序编码方法 | |
KR102328326B1 (ko) | 수질예측 고도화를 위한 분포형 모델이 적용된 수리수문 해석 장치 | |
CN113987969A (zh) | 多高程尺度水流网的陆地水动态仿真模型 | |
CN108595572A (zh) | 一种适用于城市化区域高精度提取水系与流域单元方法 | |
Ye et al. | A parallel Python-based tool for meshing watershed rivers at continental scale | |
CN113327323B (zh) | 基于散点数据的水体环境地形构建方法 | |
CN110675501A (zh) | DEM架构下的TIN向Grid模型的高效转换方法 | |
Zhang et al. | Watershed characteristics extraction and subsequent terrain analysis based on digital elevation model in flat region | |
CN116050291A (zh) | 一种台风条件下风电场大尺度气象预测降尺度方法与*** | |
CN113656852B (zh) | 一种精细化河道地形快速生成方法 | |
CN112967354B (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 |