CN106570936B - 一种基于栅格dem数据的等距离权重内插加密方法 - Google Patents
一种基于栅格dem数据的等距离权重内插加密方法 Download PDFInfo
- Publication number
- CN106570936B CN106570936B CN201610998262.0A CN201610998262A CN106570936B CN 106570936 B CN106570936 B CN 106570936B CN 201610998262 A CN201610998262 A CN 201610998262A CN 106570936 B CN106570936 B CN 106570936B
- Authority
- CN
- China
- Prior art keywords
- point
- data
- grid
- interpolation
- dem data
- 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
- G06T17/00—Three dimensional [3D] modelling, e.g. data description of 3D objects
- G06T17/05—Geographic models
Landscapes
- Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Geometry (AREA)
- Software Systems (AREA)
- Remote Sensing (AREA)
- Computer Graphics (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Image Processing (AREA)
- Image Analysis (AREA)
Abstract
一种基于栅格DEM数据的等距离权重内插加密方法,涉及地理信息技术应用领域。本发明首先根据实际应用,选择初始DEM数据,经过计算后,将该栅格DEM数据转换为各方格点的高程信息,根据地表高程连续这一特征,采用内插算法,在点矩阵中的x、y方向上选择最近邻点进行等距离内插,同时依据内插点与初始点矩阵中的最近两点间的距离确定该两已知点高程对内插点的影响权重,以此得到点密度更高的点矩阵,最后按照栅格数据特征,由栅格中四个点的高程信息生成该栅格的高程信息,从而生成精度更高的栅格DEM数据。本发明实现了简单、快速、高效的内插DEM数据的目的,扩展了现有的DEM数据集,同时丰富了DEM数据的应用范围。
Description
技术领域
本发明涉及地理信息技术应用领域,尤其是一种基于栅格DEM数据内插重采样、误差分析、种子蔓延扩散算法的DEM数据加密方法。
背景技术
地理科学和计算机科学的快速发展和有机结合,产生了基于地理信息与计算机技术的地理信息系统,它是在计算机软、硬件系统的支持下,对整个或部分地球表层(包括大气层)空间中有关地理分布数据进行采集、存储、管理、运算、分析和描述的技术系统。近30年来,地理信息系统取得了惊人的发展,其广泛应用于资源调查、环境评估、灾害预测、国土管理、水利水电等领域,其中DEM(Digital Elevation Model数字高程模型)是地理信息系统中所特有的数据模型,生成DEM的方法有很多种:(1)直接地面测量,例如用GPS、全站仪、野外测量等;(2)根据航空或航天影像通过摄影测量途径获取;(3)从现有地形图上采集。DEM数据的分布式是DEM刻化地形精确程度的重要指标,同样也是决定其使用范围的一个主要影响因素,DEM数据的分辨率是指DEM最小单元格的长度,因为DEM是离散的数据,所有(x,y)坐标都是代表一个一个的小方格,每个小方格标识出其高程,小方格长度值越小,栅格DEM的网格越密,刻化的地形程度就越精确,同时数据量也呈几合级数增长。
发明内容
本发明目的是旨在基于现有DEM数据精度的基础上,提供一种简单、快速、高效生成更高网格密度的栅格DEM数据算法,从而丰富现有DEM数据的数据集及扩展现其应用范围。随着计算机技术的发展,尤其是大数据技术的普及,使得现有的具体行业应用中使用高栅格密度的DEM数据成为可能。因此,快速高效地基于现有DEM数据生成更高栅格密度的DEM数据具有实际应用意义。
为达到上述目的,本发明采用以下技术方案:首先根据实际应用,选择初始DEM数据,经过计算后,将该栅格DEM数据转换为各方格点的高程信息(点矩阵),根据地表高程连续这一特征,采用内插算法,在点矩阵中的x、y方向上选择最近邻点进行等距离内插,同时依据内插点与初始点矩阵中的最近两点间的距离确定该两已知点高程对内插点的影响权重,以此得到点密度更高的点矩阵,最后按照栅格数据特征,由栅格中四个点的信息生成栅格信息,并以此来代替为该处高程值,从而生成精度更高的栅格DEM数据。
具体包括如下步骤:
(1)以实际DEM数据应用为导向,选择初始精度的栅格DEM数据,先将栅格信息转换为点信息,具体计算算法为将与DEM数据中栅格点相邻的四个栅格的值求平均值,以该均值代替为该点的高程值,逐行逐列依次对所有栅格进行计算,并将最后整个栅格数据所得到的点信息保存为点矩阵。
(2)对获得的点矩阵在x、y方向上进行最近邻点之间的等距离插值,通过内插点与其最近已知栅格点的距离为权重决定已知点高程对内插点高程的影响,以已知栅格点的高程信息为基础,获得插值点的高程信息,生成插值后的点密度更高的点矩阵。
(3)根据插值后所获得的密度更高的点矩阵,逆操作步骤(1)中的计算过程,将栅格中的四个端点信息转换栅格信息,以此获得连续的栅格数据,从而得到栅格密度更高的DEM数据。
(4)将初始栅格数据(父栅格数据)的高程与内插后所得到的各个栅格数据(子栅格数据)分别进行误差检验,以此来检验所得数据的一致性与均一性,通过开始时所设定的误差函数范围,对超限区域的栅格基于种子蔓延算法重新进行修正计算,以此获得拟合度更符合实际地表类型的栅格DEM数据模型。
本发明采用上述技术方案,具有以下优点:
(1)本发明以现有DEM数据为基础,利用点面数据的转换以及内插等算法来获得更高栅格密度的数据,较其他精细DEM生成方法的成本更低,算法简单高效,在实际应用中具有很高的可操作性。
(2)本发明所使用的内插算法核心思想是将原始DEM数据中的基础栅格内插出更细小的栅格数据,内插点的高程值由与其最近的两个点的高程内插所得,因此能较好的模拟出地表连续的特征。
(3)本发明中设计的算法具有很好的自适应性和可移植性,如将a*b的栅格DEM数据内插成ka*kb的栅格数据,k的值可在算法中设置成变量(在实际的应用中,一般k的最佳值为2和3),以此来自动生成原始数据k倍栅格密度的产品。同时算法中对每个原始栅格和其所生成栅格数据的高程信息逐一进行误差检测,对超限的栅格点,采用种子蔓延算法确定该栅格的最优值,进行局部平滑,可获得拟合度更好的DEM数据。
附图说明
图1是原始栅格数据生成原始点矩阵的示意图。
图2是原始点矩阵内插出高密度点矩阵的示意图。
图3是高密度点矩阵生成高密度栅格数据的示意图。
图4是插值后的栅格高程提取示意图。
图5是原始栅格数据与生成后的栅格数据误差检验示意图。
具体实施方式
下面结合附图对本发明的技术方案进行详细说明:
具体包含以下步骤:
S-1本实施例以小流域洪水淹没模拟所需的DEM数据为应用导向,选择精度为30m的栅格DEM数据,以a*b栅格中每个栅格的高程信息为基准,根据栅格数据中各栅格点与其相邻栅格的关系,以与点相邻的栅格平均值为该栅格点的高程值,得到由栅格转换而成的(a+1)*(b+1)点矩阵。具体实施过程如图1所示,图中截取为3*4的DEM数据,其中H1,H2,···,H12为DEM数据中12个栅格的高程信息,由于边缘缺少相应的参考信息,因此A点高程为H1,B点高程为(H1+H2)/2,先逐行求出第一行C、D、E的值,再依次求出第二行F、G、H、I、J的值,最后输出该栅格数据的(4*5)点矩阵信息。
S-2原始栅格DEM数据中的栅格为正方形,因此在x、y轴上的单个栅格边上的内插数量相同。根据a*b栅格数据转换所得到的(a+1)*(b+1)点矩阵,进行k值逐行逐列内插时,可得到(ka+1)*(kb+1)的加密点矩阵,其中k为内插系数,本发明算法中的关于内插系数k的初始边界条件为k=2,k=k+1(实际应用中,k的最佳值为2、3,k值越大,冗余信息越多,计算越慢)。在本实施例中,在进行第一次内插时(k=2),依据原始4*5点矩阵的点信息得到7*9点矩阵,内插过程以及内插后得到的点矩阵信息如图2所示。
S-3根据所得到的(ka+1)*(kb+1)点矩阵重新生成栅格数为ka*kb的DEM数据,此时栅格DEM数据的高程值由每个栅格的四个点的值来决定。
S-3-1本实施例中,根据7*9点矩阵,得到只有点信息的栅格图像,其中内插点的命名规则以原始点顺时针方向开始命名,如栅格ABGF中,AB边上的内插点为a1,FC上的点为f1,a1和f1上的内插点为a2,AF的内插点为a3,栅格中各点的信息对应内插后点矩阵中各点的信息,具体如图3所示。
S-3-2将S-3-1中所生成的DEM栅格点信息转换为栅格信息,每个栅格的高程为其四个点的平均值,内插后的子栅格中的H11由点A、点a1、点a2、点a3的平均值来决定,其中栅格命名规则为以原始栅格名为基准,对内插后的栅格进行顺时针命名,如H1栅格内插后生成H11、H12、H13、H14,按照算法依次逐行逐列得出如图4中所列栅格H11、H12、H13、H14、···、H121、H122、H123、H124的值。
S-4在实际应用中可根据本发明算法中k=2,k=k+1条件,基于30m的DEM数据提取15m、10m、5m的DEM数据。在本实施例中,初始k=2时,在完成S-2、S-3操作后,即已基于原始的3*4栅格数据生成6*8的DEM栅格数据,单个栅格的分辨率由30m提升为15m。
S-4-1本实施例中的k=2为内插第一步,由原始30m分辨率DEM数据内插出四个15m分辨率DEM数据,如图5所示,取其中位于同一个地理位置的栅格做误差分析,对插值前的H1与插值后的H11、H12、H13、H14进行比较。最后逐行逐列对全局栅格进行误差分析。
S-4-2根据S-4-1中的全局误差统计,对超出误差限定(由具体应用需求确定,本实验中进行了2倍插值,因此设置15m为超限值)的子栅格,先判断其为坡谷点还是坡峰点,若为坡谷点或坡峰点,则以该点为中心,搜寻其两个栅格距离以内的24个栅格,选择与待校正点处于同一线上的两个坡度最大的相邻栅格,以该两点的高程差为该待校正点的地形坡度,拟合出待校正点高程。如果为非坡谷点或非坡峰点,则基于种子蔓延算法,从其周边8个相邻栅格中选择与待校正栅格处在同一坡度线上的两个栅格,通过该坡度线上两个栅格点拟合出待校正点的值,达到平滑全局栅格DEM数据的目的。
S-4-3将平滑后的数据再进行误差检验,通过则进入下一步,否则在S-4-2中选择次级策略(次级坡度大小)进行校正。
S-5将上述结果与具体应用相结合,如果符合实际应用要求,则进入下一步。否则继续本发明算法中的k=3内插,重复上述步骤,在原始栅格边上进行三等分,内插新生成的两个点的高程以端点值为基础,以待定点与已知点的距离为权重决定因子确定已知点对内插点高程的影响程度,从而得出内插点的高程值。
S-6快速生成更高栅格密度的DEM数据内插方法可以简单、快速、高效地生成格网密度更高的栅格DEM数据,克服现有的高栅格密度DEM数据生产及应用成本过高的缺点,同时根据k值的不同,可生成不同细密度的栅格DEM数据,通过设置精度检查环节,采用种子扩散蔓延算法原理,进一步修正加密后的栅格数据,使得模拟的结果更加符合实际地表情况,极大地扩展了现有的栅格DEM数据集,在一定程度上满足实际具体应用中如洪水模拟、水文模型分析中对高栅格密度数据的应用需求,在具体的分析中具有较好的应用价值。
Claims (5)
1.一种基于栅格DEM数据的等距离权重内插加密方法,其特征在于:首先选择初始DEM数据,经过计算后,将该栅格DEM数据转换为各方格点的高程信息,根据地表高程连续这一特征,采用内插算法,在点矩阵中的x、y方向上选择最近邻点进行等距离内插,同时依据内插点与初始点矩阵中的最近两点间的距离确定初始点矩阵中最近两点的高程对内插点的影响权重,以此得到点密度更高的点矩阵,最后按照栅格数据特征,由栅格中四个点的高程信息得到该栅格的高程信息,从而生成分辨率更高的栅格DEM数据。
2.根据权利要求1所述的基于栅格DEM数据的等距离权重内插加密方法,其特征在于上述将该栅格DEM数据转换为各方格点的高程信息的计算方法为:将与DEM数据中栅格点相邻的四个栅格的值求平均值,以该均值代替为对应方格点的高程值,逐行逐列依次对所有栅格进行计算,并将最后整个栅格数据所得到的点高程信息保存为点矩阵。
3.根据权利要求1所述的基于栅格DEM数据的等距离权重内插加密方法,其特征在于上述内插算法为:对获得的点矩阵在x、y方向上进行最近邻点间的等距离插值,通过以内插点与其最近已知栅格点的距离为权重决定已知点高程对内插点高程的影响,以已知栅格点的高程信息为基础,获得插值点的高程信息,生成插值后的点密度更高的点矩阵。
4.根据权利要求1或2或3所述的基于栅格DEM数据的等距离权重内插加密方法,其特征在于根据插值后所获得的密度更高的点矩阵,逆操作将与DEM数据中栅格点相邻的四个栅格的值求平均高程值,以该平均高程值代替该点的高程值,逐行逐列依次对所有栅格进行计算,并将最后整个栅格数据所得到的点信息保存为点矩阵;将栅格中的四个端点信息转换单个栅格信息,以此获得连续的栅格数据,从而得到栅格密度更高的DEM数据。
5.根据权利要求4所述的基于栅格DEM数据的等距离权重内插加密方法,其特征在于将初始栅格数据的高程与内插后所得到的各个栅格数据分别进行误差检验,以此来检验所得数据的一致性与均一性,通过开始时所设定的误差函数范围,对超限区域的栅格基于种子蔓延算法重新进行修正计算,以此获得拟合度更符合实际地表类型的栅格DEM数据模型。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610998262.0A CN106570936B (zh) | 2016-11-14 | 2016-11-14 | 一种基于栅格dem数据的等距离权重内插加密方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610998262.0A CN106570936B (zh) | 2016-11-14 | 2016-11-14 | 一种基于栅格dem数据的等距离权重内插加密方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN106570936A CN106570936A (zh) | 2017-04-19 |
CN106570936B true CN106570936B (zh) | 2019-07-02 |
Family
ID=58542759
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201610998262.0A Active CN106570936B (zh) | 2016-11-14 | 2016-11-14 | 一种基于栅格dem数据的等距离权重内插加密方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106570936B (zh) |
Families Citing this family (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107860375A (zh) * | 2017-11-03 | 2018-03-30 | 广州地理研究所 | 一种基于三维激光扫描技术的滑坡灾害体积快速提取方法 |
CN108280880B (zh) * | 2018-01-24 | 2018-11-16 | 长春工程学院 | 一种利用遥感影像提高山体的数字高程数据分辨率的方法 |
CN109345445B (zh) * | 2018-09-03 | 2022-11-25 | 中国电子科技集团公司第二十八研究所 | 一种地理数据区域化迁移方法及系统 |
CN113077380B (zh) * | 2021-04-09 | 2023-12-08 | 余凡 | 一种栅格数据快速插值方法 |
CN117011406A (zh) * | 2023-06-30 | 2023-11-07 | 国家基础地理信息中心 | 地图制作方法、装置、电子设备及计算机可读存储介质 |
CN117456123B (zh) * | 2023-12-26 | 2024-03-22 | 北京云庐科技有限公司 | 基于点数据生成数字高程模型的方法 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102324106A (zh) * | 2011-06-02 | 2012-01-18 | 武汉大学 | 一种顾及地表光谱信息的sfs三维重建加密稀疏dem方法 |
CN102521882A (zh) * | 2011-12-05 | 2012-06-27 | 西北工业大学 | 基于离散高程和自适应混合加权得到海床地形数据的方法 |
CN103399990A (zh) * | 2013-07-18 | 2013-11-20 | 北京工业大学 | 城市排水仿真模拟系统的道路精细离散网格构建方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20140257700A1 (en) * | 2013-03-08 | 2014-09-11 | The Government Of The United States Of America, As Represented By The Secretary Of The Navy | System and method for estimating uncertainty for geophysical gridding routines lacking inherent uncertainty estimation |
-
2016
- 2016-11-14 CN CN201610998262.0A patent/CN106570936B/zh active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102324106A (zh) * | 2011-06-02 | 2012-01-18 | 武汉大学 | 一种顾及地表光谱信息的sfs三维重建加密稀疏dem方法 |
CN102521882A (zh) * | 2011-12-05 | 2012-06-27 | 西北工业大学 | 基于离散高程和自适应混合加权得到海床地形数据的方法 |
CN103399990A (zh) * | 2013-07-18 | 2013-11-20 | 北京工业大学 | 城市排水仿真模拟系统的道路精细离散网格构建方法 |
Non-Patent Citations (1)
Title |
---|
不同采样与插值方法的DEM质量分析;刘沁萍 等;《测绘科学技术学报》;20091031;第26卷(第5期);第360-362页 |
Also Published As
Publication number | Publication date |
---|---|
CN106570936A (zh) | 2017-04-19 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106570936B (zh) | 一种基于栅格dem数据的等距离权重内插加密方法 | |
CN108763825B (zh) | 一种模拟复杂地形的风场的数值模拟方法 | |
US20120089381A1 (en) | Automated method and system for predicting high resolution tidal heights and currents in coastal zones | |
CN102819568A (zh) | 基于地形采样点位置的海底地形数据建立方法 | |
CN107784165B (zh) | 基于光伏电站的地表温度场多尺度资料同化方法 | |
CN110045432A (zh) | 基于3d-glq的球坐标系下重力场正演方法及三维反演方法 | |
CN108664705B (zh) | 一种基于OpenFOAM的模拟复杂地形地表粗糙度的方法 | |
CN112069449B (zh) | 基于初值集合的天气预报方法及装置 | |
Herzog et al. | Least-cost kernel density estimation and interpolation-based density analysis applied to survey data | |
Rabah et al. | The use of minimum curvature surface technique in geoid computation processing of Egypt | |
Ligas et al. | Geostatistical prediction of a local geometric geoid-kriging and cokriging with the use of EGM2008 geopotential model | |
CN116842877A (zh) | 一种基于多源数据综合利用的小尺度三维风场重构算法 | |
Fazilova et al. | Comparative study of interpolation methods in development of local geoid | |
CN106546958B (zh) | 一种优化的雷达数据同化方法 | |
CN117517213A (zh) | 一种矿区土壤碳含量的采样方法 | |
CN110287620B (zh) | 适用于地表观测面的球坐标系密度界面正演方法及系统 | |
Wijemannage et al. | Comparison of spatial interpolation methods for rainfall data over Sri Lanka | |
Du et al. | Real-time tropospheric delay map retrieval using sparse GNSS stations | |
Zhao et al. | Rectification methods comparison for the ASTER GDEM V2 data using the ICESat/GLA14 data in the Lvliang mountains, China | |
CN115841512A (zh) | 一种复杂地表真实面积的计算方法 | |
Zhang et al. | Spatial interpolation of air temperature with ANUSPLIN in Three Gorges Reservoir Area | |
CN107220314A (zh) | 一种构建数字高程模型的方法 | |
Dawod et al. | Accuracy assessment of the PGM17 global geopotential model: a case study of Egypt and Northeast Africa | |
CN118068453B (zh) | 基于多源预报产品的大释用短期降水预报方法和系统 | |
Wang et al. | A new method of generating grid DEM from contour lines |
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 |