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

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

Info

Publication number
CN109448119A
CN109448119A CN201811224002.3A CN201811224002A CN109448119A CN 109448119 A CN109448119 A CN 109448119A CN 201811224002 A CN201811224002 A CN 201811224002A CN 109448119 A CN109448119 A CN 109448119A
Authority
CN
China
Prior art keywords
dem
point
file
data
data collection
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.)
Granted
Application number
CN201811224002.3A
Other languages
English (en)
Other versions
CN109448119B (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

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 (4)

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

* 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
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技术平台

Patent Citations (8)

* 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
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 (7)

* Cited by examiner, † Cited by third party
Title
YANG HAIXIA ET AL.: "An Implementation about GIS Data Interoperability and Service Interoperability", 《2012 INTERNATIONAL CONFERENCE ON COMPUTER SCIENCE AND SERVICE SYSTEM》 *
刘晓艳等: "地理信息系统与虚拟现实之间的数据互操作研究", 《测绘通报》 *
刘燕等: "基于python的数字高程模型高程精度检测方法", 《青海国土经略》 *
徐晓庆等: "不同数字高程模型电子地图的数据格式转换", 《信息技术》 *
徐青著: "《地形三维可视化技术》", 30 August 2000, 北京：测绘出版社 *
杨海霞: "基于RS和GIS的建筑物空间分布格网化方法研究", 《地震》 *
蒋立辉等: "基于SRTM数字高程模型的WRF低空风场模拟研究", 《中国民航飞行学院学报》 *

Also Published As

Publication number Publication date
CN109448119B (zh) 2022-04-19

Similar Documents

Publication Publication Date Title
Vivoni et al. Generation of triangulated irregular networks based on hydrological similarity
Chang Introduction to geographic information systems
US20220067117A1 (en) System and method for linking data records for parcels
CN105893590B (zh) 一种用于数字地形分析建模知识案例化自动处理方法
CN115205481A (zh) 一种基于图神经网络的频谱地图构建方法及系统
CN105070185B (zh) 一种点要素群自动综合方法
CN109388843A (zh) 一种基于vtk的桁架天线的可视化系统及方法、终端
CN112066997A (zh) 高清航线地图的导出方法及系统
Donnay et al. Modelling geographical distributions in urban areas
Liu et al. Approximation theory applied to DEM vertical accuracy assessment
JP5403726B2 (ja) 浸水深調査システム及びプログラム
CN103280152A (zh) 一种基于经纬度坐标转化码的楼宇建筑编码方法
CN103345769B (zh) 全路径遥感影像专题制图法
CN109448119A (zh) 一种地理信息系统开发中关于数字高程模型（dem）应用的方法
CN111446968B (zh) 矢量空间数据多级压缩方法
CN110887495B (zh) 云平台实时路况在城市应急gis平台中应用的方法
CN112184900B (zh) 高程数据的确定方法、装置及存储介质
CN108615452A (zh) 一种基于人车轨迹点多分辨率处理的未知道路提取方法
CN117291000A (zh) 一种国土空间规划大数据分析辅助模型
Liu et al. The “M” in digital elevation models
JP2008309632A (ja) 浸水深現地調査補助システム
CN114443787B (zh) 大气污染特征识别方法和装置
Bednarczyk et al. Forest cover change or misinterpretation? On dependent and independent vectorisation approaches
KR100964317B1 (ko) 지알디 파일을 지오매니아 지시디 파일로 변환하는 시스템및 방법
CN106443732A (zh) 一种基于gps的路径图绘制方法和系统

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