CN109448119B - 一种地理信息系统开发中关于数字高程模型（dem）应用的方法 - Google Patents

一种地理信息系统开发中关于数字高程模型（dem）应用的方法 Download PDF

Info

Publication number
CN109448119B
CN109448119B CN201811224002.3A CN201811224002A CN109448119B CN 109448119 B CN109448119 B CN 109448119B CN 201811224002 A CN201811224002 A CN 201811224002A CN 109448119 B CN109448119 B CN 109448119B
Authority
CN
China
Prior art keywords
dem
file
data
point
row
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
Application number
CN201811224002.3A
Other languages
English (en)
Other versions
CN109448119A (zh
Inventor
郭逸超
郭清
杨海霞
莫莉
刘治军
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Shenzhen Gongkan Geotechnical Group Co Ltd
Original Assignee
Shenzhen Gongkan Geotechnical Group Co Ltd
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Shenzhen Gongkan Geotechnical Group Co Ltd filed Critical Shenzhen Gongkan Geotechnical Group Co Ltd
Priority to CN201811224002.3A priority Critical patent/CN109448119B/zh
Publication of CN109448119A publication Critical patent/CN109448119A/zh
Application granted granted Critical
Publication of CN109448119B publication Critical patent/CN109448119B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

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)
  • Processing Or Creating Images (AREA)
  • Instructional Devices (AREA)

Abstract

本发明涉及数字高程模型(DEM)的技术领域，公开了一种地理信息系统开发中关于数字高程模型(DEM)应用的方法，包括以下步骤：(1)、将ASCII码的DEM数据集文件转换为二进制数据文件；(2)、在地理信息系统开发中，需要获知已知大地坐标(X，Y)的点位的高程值时，根据所述ASCII码的DEM数据集文件的头文件参数，计算该点位在所述DEM数据集文件中所在的位置(行×列)；(3)、根据计算出的该点位在所述DEM数据集文件中所在的位置(行×列)，计算该点位在所述二进制数据文件中的数据位；(4)、根据计算出的所述数据位直接从所述二进制数据文件中读取该数据位的高程值；直接从数据集文件中读取所需要的高程值，极大地提高了计算机的运行效率。

Description

一种地理信息系统开发中关于数字高程模型（DEM）应用的 方法
技术领域
本发明涉及数字高程模型(DEM)的技术领域，尤其是一种地理信息系统开发中关于数字高程模型(DEM)应用的方法。
背景技术
数字高程模型(DEM)是指地面上按一定间距排列的规则格网上的高程数据集，它的数据标准包括头文件和数据集文件；它的主要应用是查询和分析地面高程、坡度、坡向、地面通视等三维分析应用，如土木工程的土方计算和无线电信号的覆盖分析等。数字高程模型一般是从野外实测地面特征点、特征线的高程和等高线，按一定的内插方法计算出规则格网点的高程值。从上世纪20 年代起，随着计算机技术的发展，使得DEM内插的数值计算采用计算机编程求解，有些原本属于数学领域的内插方法借助计算机技术广泛应用工程实践中。
目前DEM数据格式有多种，其中最常用的有USGS的*.dem格式和ESRI Arc/Info标准格网lattice转换的ASCII文件，这两种数据都是以ASCII码文件。
由于ASCII码文件是在计算机系统中是顺序文件，目前大多数地理信息系统(GIS)中调用DEM数据时都是将DEM数据集读入到计算机内存中，需要时从内存中提取，这样将占用大量的计算机内存资源。
发明内容
本发明的目的在于提供一种地理信息系统开发中关于数字高程模型(DEM) 应用的方法，旨在解决现有技术中地理信息系统(GIS)调用DEM数据占用计算机内存资源，计算机运行效率低的问题。
本发明是这样实现的，一种地理信息系统开发中关于数字高程模型(DEM) 应用的方法，包括以下步骤：
(1)、将ASCII码的DEM数据集文件转换为二进制数据文件(双精度或单精度)；
(2)、在地理信息系统开发中，需要获知已知大地坐标(X，Y)的点位的高程值时，根据所述ASCII码的DEM数据集文件的头文件参数，计算该点位在所述DEM数据集文件中所在的位置(行×列)；
(3)、根据计算出的该点位在所述DEM数据集文件中所在的位置(行×列)，计算该点位在所述二进制数据文件中的数据位；
(4)、根据计算出的所述数据位直接从所述二进制数据文件中读取该数据位的高程值。
进一步地，所述ASCII码的DEM数据集文件的前六行分别表示：DEM的总列数、DEM的总行数、DEM的左下角X坐标、DEM的左下角Y坐标、DEM的格网间距、DEM的无数据区数据(-9999)；所述ASCII码的DEM数据集文件的第七行开始则每一个数据代表一个标准格网的高程值，以左下角开始，逐列逐行存放。
进一步地，在所述步骤(2)中，计算该点位在所述DEM数据集文件中所在的位置(行×列)的计算公式为：
点位所在的列数：nhcell＝CInt[(X-xllcorner)/cellsize+0.5]+1；
点位所在的行数：nvcell＝CInt[(Y-yllcorner)/cellsize+0.5]+1；
式中，xllcorner为DEM的左下角X坐标，cellsize为DEM的格网间距，X 为该点位大地坐标(X，Y)中的X，yllcorner为DEM的左下角Y坐标，Y为该点位大地坐标(X，Y)中的Y。
进一步地，在所述步骤(3)中，计算该点位在所述二进制数据文件中的数据位的计算公式为：
点位在二进制数据文件中的数据位为：
ncell＝CLng[(nrows-nvcell)*ncols+nhcell+6]；
式中，nrows为DEM的总行数，nvcell为所述步骤(2)中所求得的点位所在的行数，ncols为DEM的总列数，nhcell为所述步骤(2)中所求得的点位所在的列数。
与现有技术相比，本发明提供的一种地理信息系统开发中关于数字高程模型(DEM)应用的方法，将DEM数据集以文件方式存放在任何地方，如计算机硬盘、服务器、网络空间等，在GIS系统中采用一定的算法，不需要将DEM数据集文件全部读入到计算机内存中，直接从数据集文件中读取所需要的高程值，极大地提高了计算机的运行效率。
附图说明
图1是本发明实施例提供的一种地理信息系统开发中关于数字高程模型 (DEM)应用的方法的流程示意图；
图2是本发明实施例提供的一种地理信息系统开发中关于数字高程模型 (DEM)应用的方法的具体应用示意图。
具体实施方式
为了使本发明的目的、技术方案及优点更加清楚明白，以下结合附图及实施例，对本发明进行进一步详细说明。应当理解，此处所描述的具体实施例仅仅用以解释本发明，并不用于限定本发明。
本实施例的附图中相同或相似的标号对应相同或相似的部件；在本发明的描述中，需要理解的是，若有术语“上”、“下”、“左”、“右”等指示的方位或位置关系为基于附图所示的方位或位置关系，仅是为了便于描述本发明和简化描述，而不是指示或暗示所指的装置或元件必须具有特定的方位、以特定的方位构造和操作，因此附图中描述位置关系的用语仅用于示例性说明，不能理解为对本专利的限制，对于本领域的普通技术人员而言，可以根据具体情况理解上述术语的具体含义。
以下结合具体实施例对本发明的实现进行详细的描述。
参照图1-2所示，为本发明提供较佳实施例。
本发明提供的一种地理信息系统开发中关于数字高程模型(DEM)应用的方法，用于解决现有技术中地理信息系统(GIS)调用DEM数据占用计算机内存资源，计算机运行效率低的问题。
一种地理信息系统开发中关于数字高程模型(DEM)应用的方法，包括以下步骤：
(1)、将ASCII码的DEM数据集文件转换为二进制数据文件(双精度或单精度)；
(2)、在地理信息系统开发中，需要获知已知大地坐标(X，Y)的点位的高程值时，根据ASCII码的DEM数据集文件的头文件参数，计算该点位在DEM 数据集文件中所在的位置(行×列)；
(3)、根据计算出的该点位在DEM数据集文件中所在的位置(行×列)，计算该点位在二进制数据文件中的数据位；
(4)、根据计算出的数据位直接从二进制数据文件中读取该数据位的高程值。
上述提供的一种地理信息系统开发中关于数字高程模型(DEM)应用的方法，将DEM数据集以文件方式存放在任何地方，如计算机硬盘、服务器、网络空间等，在GIS系统中采用一定的算法，不需要将DEM数据集文件全部读入到计算机内存中，直接从数据集文件中读取所需要的高程值，极大地提高了计算机的运行效率。
进一步地，ASCII码的DEM数据集文件的前六行分别表示：DEM的总列数、 DEM的总行数、DEM的左下角X坐标、DEM的左下角Y坐标、DEM的格网间距、DEM的无数据区数据(-9999)；ASCII码的DEM数据集文件的第七行开始则每一个数据代表一个标准格网的高程值，以左下角开始，逐列逐行存放。
ASCII码的DEM数据集文件的数据格式见下：
ncols 8
nrows 10
xllcorner 136342.5
yllcorner 25932.5
cellsize 5
NODATA_value-9999
data data data data data……
…………………………
进一步地，在步骤(2)中，计算该点位在DEM数据集文件中所在的位置(行×列)的计算公式为：
点位所在的列数：nhcell＝CInt[(X-xllcorner)/cellsize+0.5]+1(1)；
点位所在的行数：nvcell＝CInt[(Y-yllcorner)/cellsize+0.5]+1(2)；
式中，xllcorner为DEM的左下角X坐标，cellsize为DEM的格网间距，X 为该点位大地坐标(X，Y)中的X，yllcorner为DEM的左下角Y坐标，Y为该点位大地坐标(X，Y)中的Y。
进一步地，在步骤(3)中，计算该点位在二进制数据文件中的数据位的计算公式为：
点位在二进制数据文件中的数据位为：
ncell＝CLng[(nrows-nvcell)*ncols+nhcell+6] (3)；
式中，nrows为DEM的总行数，nvcell为步骤(2)中所求得的点位所在的行数，ncols为DEM的总列数，nhcell为步骤(2)中所求得的点位所在的列数。
以上所述仅为本发明的较佳实施例而已，并不用以限制本发明，凡在本发明的精神和原则之内所作的任何修改、等同替换和改进等，均应包含在本发明的保护范围之内。

Claims (2)

1.一种地理信息系统开发中关于数字高程模型(DEM)应用的方法，其特征在于，包括以下步骤：
(1)、将ASCII码的DEM数据集文件转换为二进制数据文件，所述二进制数据文件为双精度或单精度；
(2)、在地理信息系统开发中，需要获知已知大地坐标(X，Y)的点位的高程值时，根据所述ASCII码的DEM数据集文件的头文件参数，计算该点位在所述DEM数据集文件中所在的位置(行×列)；
(3)、根据计算出的该点位在所述DEM数据集文件中所在的位置(行×列)，计算该点位在所述二进制数据文件中的数据位；
(4)、根据计算出的所述数据位直接从所述二进制数据文件中读取该数据位的高程值；
在所述步骤(2)中，计算该点位在所述DEM数据集文件中所在的位置(行×列)的计算公式为：
点位所在的列数：nhcell＝CInt[(X-xllcorner)/cellsize+0.5]+1；
点位所在的行数：nvcell＝CInt[(Y-yllcorner)/cellsize+0.5]+1；
式中，xllcorner为DEM的左下角X坐标，cellsize为DEM的格网间距，X为该点位大地坐标(X，Y)中的X，yllcorner为DEM的左下角Y坐标，Y为该点位大地坐标(X，Y)中的Y；
在所述步骤(3)中，计算该点位在所述二进制数据文件中的数据位的计算公式为：
点位在二进制数据文件中的数据位为：
ncell＝CLng[(nrows-nvcell)*ncols+nhcell+6]；
式中，nrows为DEM的总行数，nvcell为所述步骤(2)中所求得的点位所在的行数，ncols为DEM的总列数，nhcell为所述步骤(2)中所求得的点位所在的列数。
2.如权利要求1所述的一种地理信息系统开发中关于数字高程模型(DEM)应用的方法，其特征在于，所述ASCII码的DEM数据集文件的前六行分别表示：DEM的总列数、DEM的总行数、DEM的左下角X坐标、DEM的左下角Y坐标、DEM的格网间距、DEM的无数据区数据，所述DEM的无数据区数据为-9999；所述ASCII码的DEM数据集文件的第七行开始则每一个数据代表一个标准格网的高程值，以左下角开始，逐列逐行存放。
CN201811224002.3A 2018-10-19 2018-10-19 一种地理信息系统开发中关于数字高程模型（dem）应用的方法 Active CN109448119B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201811224002.3A CN109448119B (zh) 2018-10-19 2018-10-19 一种地理信息系统开发中关于数字高程模型（dem）应用的方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201811224002.3A CN109448119B (zh) 2018-10-19 2018-10-19 一种地理信息系统开发中关于数字高程模型（dem）应用的方法

Publications (2)

Publication Number Publication Date
CN109448119A CN109448119A (zh) 2019-03-08
CN109448119B true CN109448119B (zh) 2022-04-19

Family

ID=65547511

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201811224002.3A Active CN109448119B (zh) 2018-10-19 2018-10-19 一种地理信息系统开发中关于数字高程模型（dem）应用的方法

Country Status (1)

Country Link
CN (1) CN109448119B (zh)

Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
KR100912936B1 (ko) * 2008-11-27 2009-08-20 삼아항업(주) 영상지도제작을 위한 수치표고모델 수정시스템 및 그 방법
CN103823981A (zh) * 2014-02-28 2014-05-28 武汉大学 一种数字高程模型辅助的卫星影像区域网平差方法
CN106650343A (zh) * 2016-10-19 2017-05-10 南京师范大学 一种dem置乱加密与还原方法
CN106652032A (zh) * 2016-09-30 2017-05-10 电子科技大学 一种基于Linux集群平台的DEM并行等高线生成方法
CN107024683A (zh) * 2017-06-14 2017-08-08 中国科学院遥感与数字地球研究所 基于dem的sar影像直接正射校正方法及系统
CN108230326A (zh) * 2018-02-08 2018-06-29 重庆市地理信息中心 基于gpu-cpu协同的卫星影像拉花变形快速检测方法
CN108287929A (zh) * 2018-03-06 2018-07-17 霍亮 基于WebGL的三维GIS技术平台

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20060279436A1 (en) * 2005-03-25 2006-12-14 Ambroziak Russell A Digital elevation model data structure for geospatial analysis

Patent Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
KR100912936B1 (ko) * 2008-11-27 2009-08-20 삼아항업(주) 영상지도제작을 위한 수치표고모델 수정시스템 및 그 방법
CN103823981A (zh) * 2014-02-28 2014-05-28 武汉大学 一种数字高程模型辅助的卫星影像区域网平差方法
CN106652032A (zh) * 2016-09-30 2017-05-10 电子科技大学 一种基于Linux集群平台的DEM并行等高线生成方法
CN106650343A (zh) * 2016-10-19 2017-05-10 南京师范大学 一种dem置乱加密与还原方法
CN107024683A (zh) * 2017-06-14 2017-08-08 中国科学院遥感与数字地球研究所 基于dem的sar影像直接正射校正方法及系统
CN108230326A (zh) * 2018-02-08 2018-06-29 重庆市地理信息中心 基于gpu-cpu协同的卫星影像拉花变形快速检测方法
CN108287929A (zh) * 2018-03-06 2018-07-17 霍亮 基于WebGL的三维GIS技术平台

Non-Patent Citations (6)

* Cited by examiner, † Cited by third party
Title
An Implementation about GIS Data Interoperability and Service Interoperability;Yang Haixia et al.;《2012 International Conference on Computer Science and Service System》;20121231;第1603-1606页 *
不同数字高程模型电子地图的数据格式转换;徐晓庆等;《信息技术》;20150525(第5期);第48-52,56页 *
地理信息系统与虚拟现实之间的数据互操作研究;刘晓艳等;《测绘通报》;20030225(第2期);第17-20页 *
基于python的数字高程模型高程精度检测方法;刘燕等;《青海国土经略》;20171226(第6期);第67-69页 *
基于RS和GIS的建筑物空间分布格网化方法研究;杨海霞;《地震》;20150715;第35卷(第3期);第136-146页 *
基于SRTM数字高程模型的WRF低空风场模拟研究;蒋立辉等;《中国民航飞行学院学报》;20150315;第27卷(第2期);第13-17,22页 *

Also Published As

Publication number Publication date
CN109448119A (zh) 2019-03-08

Similar Documents

Publication Publication Date Title
CN110473251B (zh) 基于网格空间索引的自定义范围空间数据面积统计方法
CN111260784B (zh) 一种城市三维空间网格压缩编码方法、装置及终端设备
CN106777302B (zh) 空间地理坐标的转换方法及装置
CN110633262A (zh) 基于Spark的地图相交区域面积计算方法及系统
CN106802958B (zh) Cad数据到gis数据的转换方法及系统
CN103839291B (zh) 构建三维地形的方法及装置
CN111446968B (zh) 矢量空间数据多级压缩方法
CN109448119B (zh) 一种地理信息系统开发中关于数字高程模型（dem）应用的方法
CN110989886B (zh) 基于空间地图的三维空间网格选取方法及装置
CN108052629B (zh) 一种基于高精度dem数据的快速海陆判断方法
CN111737375B (zh) 一种基于Geohash分区的区域内容空间对象索引方法
CN110887495B (zh) 云平台实时路况在城市应急gis平台中应用的方法
CN112883130A (zh) 基于空间数据库输出地图高程数据的方法及设备、介质
CN112184900B (zh) 高程数据的确定方法、装置及存储介质
CN113806601B (zh) 一种周边兴趣点检索方法及存储介质
CN116450958A (zh) 检测点选址方法、装置、电子设备及可读存储介质
CN116542060A (zh) 一种基于法定图则规划的排水防涝规划地形模型构建方法
CN103049583B (zh) 土地监测数据的转换方法和转换接口
CN107992555B (zh) 一种栅格数据的存储及读取方法
CN111506576A (zh) 一种基于区域四叉树的地块编码方法及装置
CN113656852A (zh) 一种精细化河道地形快速生成方法
KR100964317B1 (ko) 지알디 파일을 지오매니아 지시디 파일로 변환하는 시스템및 방법
CN115239894A (zh) 一种适用于大尺度地理坐标系栅格数据的ls因子提取方法
CN117893636B (zh) 线状基础地理实体与专题地理实体边界融合方法及系统
CN117807154B (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