CN106910218B - 基于大区域栅格数字高程模型的地表表面面积计算方法 - Google Patents

基于大区域栅格数字高程模型的地表表面面积计算方法 Download PDF

Info

Publication number
CN106910218B
CN106910218B CN201710138021.3A CN201710138021A CN106910218B CN 106910218 B CN106910218 B CN 106910218B CN 201710138021 A CN201710138021 A CN 201710138021A CN 106910218 B CN106910218 B CN 106910218B
Authority
CN
China
Prior art keywords
pixel
dictionary
array
area
height
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
CN201710138021.3A
Other languages
English (en)
Other versions
CN106910218A (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.)
Beijing Institute of Surveying and Mapping
Original Assignee
Beijing Institute of Surveying and Mapping
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 Beijing Institute of Surveying and Mapping filed Critical Beijing Institute of Surveying and Mapping
Priority to CN201710138021.3A priority Critical patent/CN106910218B/zh
Publication of CN106910218A publication Critical patent/CN106910218A/zh
Application granted granted Critical
Publication of CN106910218B publication Critical patent/CN106910218B/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)
  • Image Generation (AREA)
  • Image Processing (AREA)
  • Image Analysis (AREA)

Abstract

本发明公开了一种基于大区域栅格数字高程模型的地表表面面积计算方法，它首先将原始栅格数字高程模型分割为若干个栅格数据，然后计算每个栅格数据的地表表面面积，最后累加所有栅格数据的表面面积得到原始栅格数字高程模型的表面面积。由于采用了大数据分块的策略，本发明能够高效出计算大区域栅格数字高程模型的地表表面面积，并且在像素级别尺度构建三角网，以及精确处理边缘像元，比现有计算方法计算精度更高，不仅提高了表面面积算法的精度，而且解决了计算大数据量数字高程模型表面面积计算机内存不足的问题。

Description

基于大区域栅格数字高程模型的地表表面面积计算方法
技术领域
本发明涉及一种基于大区域栅格数字高程模型的地表表面面积计算方法，属于地理信息空间分析技术领域。
背景技术
随着三维技术的飞速发展，三维地形数据大量出现在我们的工作和生活中，其中三维地形数字分析提供给我们丰富可靠的数据，其中地表表面面积是其中经常用到的一种数字分析。地表表面面积很难直接量测和计算，一般使用投影面积代替表面面积。在地形起伏较小的平原地区，地表真实面积和投影面积相差不大，但在地形起伏比较复杂的山区，两者差异明显。
目前有很多方法计算投影面积，运用投影面和坡度正割值计算表面面积、逐步分割多边形的方法计算表面面积、将一个像元分割为2个三角形统计DEM的面积等，这些方法存在一些缺点：计算精度不够或者没有考虑大数据量下的计算问题，并不能很好满足我们高效精确计算栅格数字高程模型的地表表面面积计的需求。
发明内容
本发明的目的是提供一种大区域栅格数字高程模型的地表表面面积计算方法，提高现有表面面积算法的精度，以及解决计算大数据量数字高程模型表面面积计算机内存不足的问题。
本发明解决其技术问题所采取的技术方案是：一种基于大区域栅格数字高程模型的地表表面面积计算方法，其特征是：它首先将原始栅格数字高程模型分割为若干个栅格数据，然后计算每个栅格数据的地表表面面积，最后得到原始栅格数字高程模型的表面面积。
进一步地，所述每个栅格数据的地表表面面积的计算具体过程为：
S2，查找每个栅格数据的边界像元并计算边界像元面积；
S3，对每个栅格数据构建三角网并计算三角网面积；
将步骤S2中计算得到的边界像元面积与步骤S3中计算得到的三角网面积进行相加，得到该栅格数据的地表表面面积。
进一步地，所述的原始栅格数字高程模型的表面面积通过累加所有栅格数据的表面面积得到。
进一步地，所述步骤S2的具体步骤如下：
S2.1，将每个栅格数据转化为高程数组Array[height][width]，with为栅格数据的水平方向像元个数，height为栅格数据的垂直方向像元个数，Array[x][y]表示栅格数据第x行第y列像元的高程值，如果像元[x,y]为空值，则取Array[x][y]为浮点数的最小值，记作noData；通过x,y和栅格数据左上角坐标计算出第x行第y列像元中心点的地理坐标；
记计算像元8邻域像元个数的函数为NeighbourCount，NeighbourCount(x,y)为计算位置[x-1,y-1]、[x-1,y-1]、[x-1,y-1]、[x-1,y-1]、[x-1,y-1]、[x-1,y-1]、[x-1,y-1]和[x-1,y-1]存在非空像元的个数，NeighbourCount(x,y)<8则代表像元(x,y)为边界像元；
声明边界像元面积字典Dictionary，键为边界点坐标，值为边界像元面积；Dictionary[x,y]＝S表示坐标为[x,y]的边界像元的面积为S；
获取一个像元的平面面积，记为S0；
设当前行i＝0,当前列j＝0：
S2.2，比较i和height的值：
如果i＝height，则退出计算；
如果i<height，则进入下一步；
S2.3，比较j和width的值：
如果j＝width，i＝i+1，j＝0，返回步骤S2.2；
如果j<width，则进入下一步；
S2.4，判断像元(i,j)是否为边界像元：
如果NeighbourCount(i,j)＝8，j，则j＝j+1，返回步骤S2.3；
如果NeighbourCount(i,j)<8，则进入下一步；
S2.5，在边界像元面积字典中加入像元[i,j]，并初始化面积为0，即Dictionary[i,j]＝0；
S2.6，j＝j+1,返回步骤S2.3。
进一步地，所述步骤S3的具体步骤如下：
S3.1，获取步骤S2.1高程数组Array[height][width]，像元空值noData和边界像元面积字典Dictionary；
构建三角网：提取每个像元的中心点作为顶点，每个顶点与相邻的8个像元中心点中的2个形成三角形；
计算边界像元面积：每个边界像元按照对角线和对边中位线分成8份，每一份面积为S0/8，计算在构建三角网之后剩余的份数，然后计算剩余份数的平面面积，即得到该边界像元的面积；
设当前行i＝0，当前列j＝0，栅格数据的地表表面面积为Sall；
S3.2，比较i和height的值：
如果i＝height，则进入步骤S3.8；
如果i<height，则进入下一步；
S3.3，比较j和width的值：
如果j＝width，i＝i+1，j＝0，返回步骤S3.2；
如果j<width，则进入下一步；
S3.4，比较Array[i][j]和noData的值：
如果Array[i][j]＝noData，j＝j+1，返回步骤S3.3；
如果Array[i][j]>noData，则进入下一步；
S3.5，判断第i-1行数据是否存在：
如果i-1<0，则进入步骤S3.6；
如果i-1>＝0，则进入下一步；
S3.5.1，判断像元Array[i-1][j+1]、Array[i][j+1]和Array[i-1][j]：
如果Array[i-1][j+1]和Array[i][j+1]存在且不等于noData，而且Array[i-1][j]等于noData，则以索引为[i,j]、[i,j+1]和[i-1,j+1]的像元的中心点构建三角网，每个顶点的高程即对应像元的高程值，计算三角形面积为Stra，Sall＝Sall+Stra，则进入下一步；
否则，进入步骤S3.5.3；
S3.5.2判断该三角形三个顶点是否为边缘像元，如果是，则计算其面积：
如果Dictionary中存在像元[i，j]，则Dictionary[i,j]＝Dictionary[i,j]+S0/8；
如果Dictionary中存在像元[i，j+1]，则Dictionary[i,j+1]＝Dictionary[i,j+1]+2*S0/8；
如果Dictionary中存在像元[i-1，j+1]，则Dictionary[i-1，j+1]＝Dictionary[i-1，j+1]+S0/8；
S3.5.3，判断像元Array[i-1][j]、Array[i-1][j+1]和Array[i][j+1]：
如果Array[i-1][j]和Array[i-1][j+1]存在且不等于noData，而且Array[i][j+1]不存在或者等于noData，则以索引为[i,j]、[i-1,j]和[i-1,j+1]的像元的中心点构建三角网，每个顶点的高程即对应像元的高程值，计算三角形面积为Stra，Sall＝Sall+Stra，则进入下一步；
否则，进入步骤S3.6
S3.5.4判断该三角形三个顶点是否为边缘像元，如果是，则计算其面积：
如果Dictionary中存在像元[i，j]，则Dictionary[i,j]＝Dictionary[i,j]+S0/8；
如果Dictionary中存在像元[i-1,j]，则Dictionary[i-1,j]＝Dictionary[i-1,j]+2*S0/8；
如果Dictionary中存在像元[i-1,j+1]，则Dictionary[i-1,j+1]＝Dictionary[i-1，j+1]+S0/8；
S3.6，判断第i+1行数据是否存在：
如果i+1>＝height，则退出计算；
如果i+1<height，则进入下一步；
S3.6.1，判断像元Array[i][j+1]和Array[i+1][j+1]：
如果Array[i][j+1]和Array[i+1][j+1]存在且不等于noData，则以索引为[i,j]、[i,j+1]和[i+1,j+1]的像元的中心点构建三角网，每个顶点的高程即对应像元的高程值，计算三角形面积为Stra，Sall＝Sall+Stra，则进入下一步；
否则，进入步骤S3.6.3；
S3.6.2，判断该三角形三个顶点是否为边缘像元，如果是，则计算其面积：
如果Dictionary中存在像元[i，j]，则Dictionary[i,j]＝Dictionary[i,j]+S0/8；
如果Dictionary中存在像元[i,j+1]，则Dictionary[i,j+1]＝Dictionary[i,j+1]+2*S0/8；
如果Dictionary中存在像元[i+1,j+1]，则Dictionary[i+1,j+1]＝Dictionary[i+1,j+1]+S0/8；
S3.6.3，判断像元Array[i+1][j]和Array[i+1][j+1]：
如果Array[i+1][j]和Array[i+1][j+1]存在且不等于noData，则以索引为[i,j]、[i+1,j]和[i+1,j+1]的像元的中心点构建三角网，每个顶点的高程即对应像元的高程值，计算三角形面积为Stra，Sall＝Sall+Stra，则进入下一步；
否则，进入步骤S3.7；
S3.6.4，判断该三角形三个顶点是否为边缘像元，如果是，则计算其面积：
如果Dictionary中存在像元[i，j]，则Dictionary[i,j]＝Dictionary[i,j]+S0/8；
如果Dictionary中存在像元[i+1,j]，则Dictionary[i+1,j]＝Dictionary[i+1,j]+2*S0/8；
如果Dictionary中存在像元[i+1,j+1]，则Dictionary[i+1,j+1]＝Dictionary[i+1,j+1]+S0/8；
S3.7，j＝j+1,返回步骤S3.3；
S3.8，对边界像元面积字典Dictionary中所有项求和得到Sdic，并加入到Sall，即得到整个栅格数据的地表表面面积Sall＝Sall+Sdic。
本发明提供的另一种基于大区域栅格数字高程模型的地表表面面积计算方法，其特征是：包括以下步骤：
S1，对原始栅格数字高程模型进行分块处理，生成若干个栅格数据；
S2，查找每个栅格数据的边界像元并计算边界像元面积；
S3，对每个栅格数据构建三角网并计算三角网面积；
S4，计算每个栅格数据的地表表面面积；
S5，累加所有栅格数据的表面面积即得到原始栅格数字高程模型的表面面积。
进一步地，在步骤S1中，对原始栅格数字高程模型进行格网化分割，生成若干个大小相同的栅格数据。
进一步地，在步骤S4中，每个栅格数据的地表表面面积为每块栅格数据的三角网面积和边界像元面积之和。
进一步地，所述步骤S2的具体步骤如下：
S2.1，将每个栅格数据转化为高程数组Array[height][width]，with为栅格数据的水平方向像元个数，height为栅格数据的垂直方向像元个数，Array[x][y]表示栅格数据第x行第y列像元的高程值，如果像元[x,y]为空值，则取Array[x][y]为浮点数的最小值，记作noData；通过x,y和栅格数据左上角坐标计算出第x行第y列像元中心点的地理坐标；
记计算像元8邻域像元个数的函数为NeighbourCount，NeighbourCount(x,y)为计算位置[x-1,y-1]、[x-1,y-1]、[x-1,y-1]、[x-1,y-1]、[x-1,y-1]、[x-1,y-1]、[x-1,y-1]和[x-1,y-1]存在非空像元的个数，NeighbourCount(x,y)<8则代表像元(x,y)为边界像元；
声明边界像元面积字典Dictionary，键为边界点坐标，值为边界像元面积；Dictionary[x,y]＝S表示坐标为[x,y]的边界像元的面积为S；
获取一个像元的平面面积，记为S0；
设当前行i＝0,当前列j＝0：
S2.2，比较i和height的值：
如果i＝height，则退出计算；
如果i<height，则进入下一步；
S2.3，比较j和width的值：
如果j＝width，i＝i+1，j＝0，返回步骤S2.2；
如果j<width，则进入下一步；
S2.4，判断像元(i,j)是否为边界像元：
如果NeighbourCount(i,j)＝8，j，则j＝j+1，返回步骤S2.3；
如果NeighbourCount(i,j)<8，则进入下一步；
S2.5，在边界像元面积字典中加入像元[i,j]，并初始化面积为0，即Dictionary[i,j]＝0；
S2.6，j＝j+1,返回步骤S2.3。
进一步地，所述步骤S3的具体步骤如下：
S3.1，获取步骤S2.1高程数组Array[height][width]，像元空值noData和边界像元面积字典Dictionary；
构建三角网：提取每个像元的中心点作为顶点，每个顶点与相邻的8个像元中心点中的2个形成三角形；
计算边界像元面积：每个边界像元按照对角线和对边中位线分成8份，每一份面积为S0/8，计算在构建三角网之后剩余的份数，然后计算剩余份数的平面面积，即得到该边界像元的面积；
设当前行i＝0，当前列j＝0，栅格数据的地表表面面积为Sall；
S3.2，比较i和height的值：
如果i＝height，则进入步骤S3.8；
如果i<height，则进入下一步；
S3.3，比较j和width的值：
如果j＝width，i＝i+1，j＝0，返回步骤S3.2；
如果j<width，则进入下一步；
S3.4，比较Array[i][j]和noData的值：
如果Array[i][j]＝noData，j＝j+1，返回步骤S3.3；
如果Array[i][j]>noData，则进入下一步；
S3.5，判断第i-1行数据是否存在：
如果i-1<0，则进入步骤S3.6；
如果i-1>＝0，则进入下一步；
S3.5.1，判断像元Array[i-1][j+1]、Array[i][j+1]和Array[i-1][j]：
如果Array[i-1][j+1]和Array[i][j+1]存在且不等于noData，而且Array[i-1][j]等于noData，则以索引为[i,j]、[i,j+1]和[i-1,j+1]的像元的中心点构建三角网，每个顶点的高程即对应像元的高程值，计算三角形面积为Stra，Sall＝Sall+Stra，则进入下一步；
否则，进入步骤S3.5.3；
S3.5.2判断该三角形三个顶点是否为边缘像元，如果是，则计算其面积：
如果Dictionary中存在像元[i，j]，则Dictionary[i,j]＝Dictionary[i,j]+S0/8；
如果Dictionary中存在像元[i，j+1]，则Dictionary[i,j+1]＝Dictionary[i,j+1]+2*S0/8；
如果Dictionary中存在像元[i-1，j+1]，则Dictionary[i-1，j+1]＝Dictionary[i-1，j+1]+S0/8；
S3.5.3，判断像元Array[i-1][j]、Array[i-1][j+1]和Array[i][j+1]：
如果Array[i-1][j]和Array[i-1][j+1]存在且不等于noData，而且Array[i][j+1]不存在或者等于noData，则以索引为[i,j]、[i-1,j]和[i-1,j+1]的像元的中心点构建三角网，每个顶点的高程即对应像元的高程值，计算三角形面积为Stra，Sall＝Sall+Stra，则进入下一步；
否则，进入步骤S3.6
S3.5.4判断该三角形三个顶点是否为边缘像元，如果是，则计算其面积：
如果Dictionary中存在像元[i，j]，则Dictionary[i,j]＝Dictionary[i,j]+S0/8；
如果Dictionary中存在像元[i-1,j]，则Dictionary[i-1,j]＝Dictionary[i-1,j]+2*S0/8；
如果Dictionary中存在像元[i-1,j+1]，则Dictionary[i-1,j+1]＝Dictionary[i-1，j+1]+S0/8；
S3.6，判断第i+1行数据是否存在：
如果i+1>＝height，则退出计算；
如果i+1<height，则进入下一步；
S3.6.1，判断像元Array[i][j+1]和Array[i+1][j+1]：
如果Array[i][j+1]和Array[i+1][j+1]存在且不等于noData，则以索引为[i,j]、[i,j+1]和[i+1,j+1]的像元的中心点构建三角网，每个顶点的高程即对应像元的高程值，计算三角形面积为Stra，Sall＝Sall+Stra，则进入下一步；
否则，进入步骤S3.6.3；
S3.6.2，判断该三角形三个顶点是否为边缘像元，如果是，则计算其面积：
如果Dictionary中存在像元[i，j]，则Dictionary[i,j]＝Dictionary[i,j]+S0/8；
如果Dictionary中存在像元[i,j+1]，则Dictionary[i,j+1]＝Dictionary[i,j+1]+2*S0/8；
如果Dictionary中存在像元[i+1,j+1]，则Dictionary[i+1,j+1]＝Dictionary[i+1,j+1]+S0/8；
S3.6.3，判断像元Array[i+1][j]和Array[i+1][j+1]：
如果Array[i+1][j]和Array[i+1][j+1]存在且不等于noData，则以索引为[i,j]、[i+1,j]和[i+1,j+1]的像元的中心点构建三角网，每个顶点的高程即对应像元的高程值，计算三角形面积为Stra，Sall＝Sall+Stra，则进入下一步；
否则，进入步骤S3.7；
S3.6.4，判断该三角形三个顶点是否为边缘像元，如果是，则计算其面积：
如果Dictionary中存在像元[i，j]，则Dictionary[i,j]＝Dictionary[i,j]+S0/8；
如果Dictionary中存在像元[i+1,j]，则Dictionary[i+1,j]＝Dictionary[i+1,j]+2*S0/8；
如果Dictionary中存在像元[i+1,j+1]，则Dictionary[i+1,j+1]＝Dictionary[i+1,j+1]+S0/8；
S3.7，j＝j+1,返回步骤S3.3；
S3.8，对边界像元面积字典Dictionary中所有项求和得到Sdic，并加入到Sall，即得到整个栅格数据的地表表面面积Sall＝Sall+Sdic。
本发明的有益效果是：本发明首先对栅格数字高程模型分块；然后通过查找每块数据边界像元以及对每块数据构建三角网来计算每块数据的地表表面面积；最后对所有分块数据面积求和，从而得到原始栅格数字高程模型的表面面积。由于采用了大数据分块的策略，本发明能够高效出计算大区域栅格数字高程模型的地表表面面积，并且在像素级别尺度构建三角网，以及精确处理边缘像元，比现有计算方法计算精度更高，不仅提高了表面面积算法的精度，而且解决了计算大数据量数字高程模型表面面积计算机内存不足的问题。
附图说明
图1是本发明的一种方法流程图；
图2是本发明的另一种方法流程图；
图3是本发明对原始栅格数字高程模型进行格网化分割示意图；
图4是本发明所述网格数据的示意图；
图5(a)为本发明所述网格数据中像元与其正左方和正下方像元构建三角形的示例图；
图5(b)为本发明所述网格数据中像元与其与正下方和右下方像元构建三角形示例图；
图5(c)为本发明所述网格数据中像元与其与正下方和左下方像元构建三角形示例图；
图5(d)为本发明所述网格数据中像元与其与正右方和正下方像元构建三角形的示例图；
图6为本发明的栅格数字高程模型与生成的三角网的二维示意图；
图7为本发明的栅格数字高程模型的三维示意图；
图8为本发明的生成三角网的三维示意图。
具体实施方式
为能清楚说明本方案的技术特点，下面通过具体实施方式，并结合其附图，对本发明进行详细阐述。下文的公开提供了许多不同的实施例或例子用来实现本发明的不同结构。为了简化本发明的公开，下文中对特定例子的部件和设置进行描述。此外，本发明可以在不同例子中重复参考数字和/或字母。这种重复是为了简化和清楚的目的，其本身不指示所讨论各种实施例和/或设置之间的关系。应当注意，在附图中所图示的部件不一定按比例绘制。本发明省略了对公知组件和处理技术及工艺的描述以避免不必要地限制本发明。
针对现有的评价公共设施均衡性方法的不足，为了实现城市公共设施的公平配置和分布均衡的目标，本发明合理地考虑了人与公共设施之间作用的跨区域性问题，既适用于区域之间也适用于区域内部，直观易懂，能够准确地评价一个城市的公共设施的服务辐射范围、服务能力和均衡性。
如图1、图3至图8所示，本发明的一种基于大区域栅格数字高程模型的地表表面面积计算方法，它首先将原始栅格数字高程模型分割为若干个栅格数据(栅格数据的大小可以根据计算机性能调整，一般取512*512)，然后计算每个栅格数据的地表表面面积，最后累加所有栅格数据的表面面积得到原始栅格数字高程模型的表面面积。
进一步地，所述每个栅格数据的地表表面面积的计算具体过程为：
S2，查找每个栅格数据的边界像元并计算边界像元面积；
S3，对每个栅格数据构建三角网并计算三角网面积；
将步骤S2中计算得到的边界像元面积与步骤S3中计算得到的三角网面积进行相加，得到该栅格数据的地表表面面积。
如图2至图8所示，本发明提供的另一种基于大区域栅格数字高程模型的地表表面面积计算方法，它包括以下步骤：
S1，对原始栅格数字高程模型进行分块处理，生成若干个栅格数据(栅格数据的大小可以根据计算机性能调整，一般取512*512)；
S2，查找每个栅格数据的边界像元并计算边界像元面积；
S3，对每个栅格数据构建三角网并计算三角网面积；
S4，计算每个栅格数据的地表表面面积；
S5，重复步骤S2至步骤S4直至遍历完所述生成的栅格数据为止，并累加所有栅格数据的表面面积即得到原始栅格数字高程模型的表面面积。
进一步地，在步骤S1中，对原始栅格数字高程模型进行格网化分割，生成若干个大小相同的栅格数据。
进一步地，在步骤S4中，每个栅格数据的地表表面面积为每块栅格数据的三角网面积和边界像元面积之和。
进一步地，在上述两种计算方法中，所述步骤S2的具体步骤如下：
S2.1，将每个栅格数据转化为高程数组Array[height][width]，with为栅格数据的水平方向像元个数，height为栅格数据的垂直方向像元个数，Array[x][y]表示栅格数据第x行第y列像元的高程值，如果像元[x,y]为空值，则取Array[x][y]为浮点数的最小值，记作noData；通过x,y和栅格数据左上角坐标计算出第x行第y列像元中心点的地理坐标；
记计算像元8邻域像元个数的函数为NeighbourCount，NeighbourCount(x,y)为计算位置[x-1,y-1]、[x-1,y-1]、[x-1,y-1]、[x-1,y-1]、[x-1,y-1]、[x-1,y-1]、[x-1,y-1]和[x-1,y-1]存在非空像元的个数，NeighbourCount(x,y)<8则代表像元(x,y)为边界像元；
声明边界像元面积字典Dictionary，键为边界点坐标，值为边界像元面积；Dictionary[x,y]＝S表示坐标为[x,y]的边界像元的面积为S；
获取一个像元的平面面积，记为S0；
设当前行i＝0,当前列j＝0：
S2.2，比较i和height的值：
如果i＝height，则退出计算；
如果i<height，则进入下一步；
S2.3，比较j和width的值：
如果j＝width，i＝i+1，j＝0，返回步骤S2.2；
如果j<width，则进入下一步；
S2.4，判断像元(i,j)是否为边界像元：
如果NeighbourCount(i,j)＝8，j，则j＝j+1，返回步骤S2.3；
如果NeighbourCount(i,j)<8，则进入下一步；
S2.5，在边界像元面积字典中加入像元[i,j]，并初始化面积为0，即Dictionary[i,j]＝0；
S2.6，j＝j+1,返回步骤S2.3。
进一步地，在上述两种计算方法中，所述步骤S3的具体步骤如下：
S3.1，获取步骤S2.1高程数组Array[height][width]，像元空值noData和边界像元面积字典Dictionary；
构建三角网：提取每个像元的中心点作为顶点，每个顶点与相邻的8个像元中心点中的2个形成三角形；
计算边界像元面积：每个边界像元按照对角线和对边中位线分成8份，每一份面积为S0/8，计算在构建三角网之后剩余的份数，然后计算剩余份数的平面面积，即得到该边界像元的面积；
设当前行i＝0，当前列j＝0，栅格数据的地表表面面积为Sall；
S3.2，比较i和height的值：
如果i＝height，则进入步骤S3.8；
如果i<height，则进入下一步；
S3.3，比较j和width的值：
如果j＝width，i＝i+1，j＝0，返回步骤S3.2；
如果j<width，则进入下一步；
S3.4，比较Array[i][j]和noData的值：
如果Array[i][j]＝noData，j＝j+1，返回步骤S3.3；
如果Array[i][j]>noData，则进入下一步；
S3.5，判断第i-1行数据是否存在：
如果i-1<0，则进入步骤S3.6；
如果i-1>＝0，则进入下一步；
S3.5.1，判断像元Array[i-1][j+1]、Array[i][j+1]和Array[i-1][j]：
如果Array[i-1][j+1]和Array[i][j+1]存在且不等于noData，而且Array[i-1][j]等于noData，则以索引为[i,j]、[i,j+1]和[i-1,j+1]的像元的中心点构建三角网，每个顶点的高程即对应像元的高程值，计算三角形面积为Stra，Sall＝Sall+Stra，则进入下一步；
否则，进入步骤S3.5.3；
S3.5.2判断该三角形三个顶点是否为边缘像元，如果是，则计算其面积：
如果Dictionary中存在像元[i，j]，则Dictionary[i,j]＝Dictionary[i,j]+S0/8；
如果Dictionary中存在像元[i，j+1]，则Dictionary[i,j+1]＝Dictionary[i,j+1]+2*S0/8；
如果Dictionary中存在像元[i-1，j+1]，则Dictionary[i-1，j+1]＝Dictionary[i-1，j+1]+S0/8；
S3.5.3，判断像元Array[i-1][j]、Array[i-1][j+1]和Array[i][j+1]：
如果Array[i-1][j]和Array[i-1][j+1]存在且不等于noData，而且Array[i][j+1]不存在或者等于noData，则以索引为[i,j]、[i-1,j]和[i-1,j+1]的像元的中心点构建三角网，每个顶点的高程即对应像元的高程值，计算三角形面积为Stra，Sall＝Sall+Stra，则进入下一步；
否则，进入步骤S3.6
S3.5.4判断该三角形三个顶点是否为边缘像元，如果是，则计算其面积：
如果Dictionary中存在像元[i，j]，则Dictionary[i,j]＝Dictionary[i,j]+S0/8；
如果Dictionary中存在像元[i-1,j]，则Dictionary[i-1,j]＝Dictionary[i-1,j]+2*S0/8；
如果Dictionary中存在像元[i-1,j+1]，则Dictionary[i-1,j+1]＝Dictionary[i-1，j+1]+S0/8；
S3.6，判断第i+1行数据是否存在：
如果i+1>＝height，则退出计算；
如果i+1<height，则进入下一步；
S3.6.1，判断像元Array[i][j+1]和Array[i+1][j+1]：
如果Array[i][j+1]和Array[i+1][j+1]存在且不等于noData，则以索引为[i,j]、[i,j+1]和[i+1,j+1]的像元的中心点构建三角网，每个顶点的高程即对应像元的高程值，计算三角形面积为Stra，Sall＝Sall+Stra，则进入下一步；
否则，进入步骤S3.6.3；
S3.6.2，判断该三角形三个顶点是否为边缘像元，如果是，则计算其面积：
如果Dictionary中存在像元[i，j]，则Dictionary[i,j]＝Dictionary[i,j]+S0/8；
如果Dictionary中存在像元[i,j+1]，则Dictionary[i,j+1]＝Dictionary[i,j+1]+2*S0/8；
如果Dictionary中存在像元[i+1,j+1]，则Dictionary[i+1,j+1]＝Dictionary[i+1,j+1]+S0/8；
S3.6.3，判断像元Array[i+1][j]和Array[i+1][j+1]：
如果Array[i+1][j]和Array[i+1][j+1]存在且不等于noData，则以索引为[i,j]、[i+1,j]和[i+1,j+1]的像元的中心点构建三角网，每个顶点的高程即对应像元的高程值，计算三角形面积为Stra，Sall＝Sall+Stra，则进入下一步；
否则，进入步骤S3.7；
S3.6.4，判断该三角形三个顶点是否为边缘像元，如果是，则计算其面积：
如果Dictionary中存在像元[i，j]，则Dictionary[i,j]＝Dictionary[i,j]+S0/8；
如果Dictionary中存在像元[i+1,j]，则Dictionary[i+1,j]＝Dictionary[i+1,j]+2*S0/8；
如果Dictionary中存在像元[i+1,j+1]，则Dictionary[i+1,j+1]＝Dictionary[i+1,j+1]+S0/8；
S3.7，j＝j+1,返回步骤S3.3；
S3.8，对边界像元面积字典Dictionary中所有项求和得到Sdic，并加入到Sall，即得到整个栅格数据的地表表面面积Sall＝Sall+Sdic。
此外，本发明的应用范围不局限于说明书中描述的特定实施例的工艺、机构、制造、物质组成、手段、方法及步骤。从本发明的公开内容，作为本领域的普通技术人员将容易地理解，对于目前已存在或者以后即将开发出的工艺、机构、制造、物质组成、手段、方法或步骤，其中它们执行与本发明描述的对应实施例大体相同的功能或者获得大体相同的结果，依照本发明可以对它们进行应用。因此，本发明所附权利要求旨在将这些工艺、机构、制造、物质组成、手段、方法或步骤包含在其保护范围内。

Claims (4)

1.基于大区域栅格数字高程模型的地表表面面积计算方法，其特征是：它首先将原始栅格数字高程模型分割为若干个栅格数据，然后计算每个栅格数据的地表表面面积，最后得到原始栅格数字高程模型的表面面积；
所述每个栅格数据的地表表面面积的计算具体过程为：
S2，查找每个栅格数据的边界像元并计算边界像元面积；
S3，对每个栅格数据构建三角网并计算三角网面积；
将步骤S2中计算得到的边界像元面积与步骤S3中计算得到的三角网面积进行相加，得到该栅格数据的地表表面面积；
所述步骤S2的具体步骤如下：
S2.1，将每个栅格数据转化为高程数组Array[height][width]，width为栅格数据的水平方向像元个数，height为栅格数据的垂直方向像元个数，Array[x][y]表示栅格数据第x行第y列像元的高程值，如果像元[x,y]为空值，则取Array[x][y]为浮点数的最小值，记作noData；通过x,y和栅格数据左上角坐标计算出第x行第y列像元中心点的地理坐标；
记计算像元8邻域像元个数的函数为NeighbourCount，NeighbourCount(x,y)为计算位置[x-1,y-1]、[x,y-1]、[x+1,y-1]、[x-1,y]、[x+1,y]、[x-1,y+1]、[x,y+1]和[x+1,y+1]存在非空像元的个数，NeighbourCount(x,y)<8则代表像元(x,y)为边界像元；
声明边界像元面积字典Dictionary，键为边界点坐标，值为边界像元面积；Dictionary[x,y]＝S表示坐标为[x,y]的边界像元的面积为S；
获取一个像元的平面面积，记为S0；
设当前行i＝0,当前列j＝0：
S2.2，比较i和height的值：
如果i＝height，则退出计算；
如果i<height，则进入下一步；
S2.3，比较j和width的值：
如果j＝width，i＝i+1，j＝0，返回步骤S2.2；
如果j<width，则进入下一步；
S2.4，判断像元(i,j)是否为边界像元：
如果NeighbourCount(i,j)＝8，则j＝j+1，返回步骤S2.3；
如果NeighbourCount(i,j)<8，则进入下一步；
S2.5，在边界像元面积字典中加入像元[i,j]，并初始化面积为0，即Dictionary[i,j]＝0；
S2.6，j＝j+1,返回步骤S2.3；
所述步骤S3的具体步骤如下：
S3.1，获取步骤S2.1高程数组Array[height][width]，像元空值noData和边界像元面积字典Dictionary；
构建三角网：提取每个像元的中心点作为顶点，每个顶点与相邻的8个像元中心点中的2个形成三角形；
计算边界像元面积：每个边界像元按照对角线和对边中位线分成8份，每一份面积为S0/8，计算在构建三角网之后剩余的份数，然后计算剩余份数的平面面积，即得到该边界像元的面积；
设当前行i＝0，当前列j＝0，栅格数据的地表表面面积为Sall；
S3.2，比较i和height的值：
如果i＝height，则进入步骤S3.8；
如果i<height，则进入下一步；
S3.3，比较j和width的值：
如果j＝width，i＝i+1，j＝0，返回步骤S3.2；
如果j<width，则进入下一步；
S3.4，比较Array[i][j]和noData的值：
如果Array[i][j]＝noData，j＝j+1，返回步骤S3.3；
如果Array[i][j]>noData，则进入下一步；
S3.5，判断第i-1行数据是否存在：
如果i-1<0，则进入步骤S3.6；
如果i-1>＝0，则进入下一步；
S3.5.1，判断像元Array[i-1][j+1]、Array[i][j+1]和Array[i-1][j]：
如果Array[i-1][j+1]和Array[i][j+1]存在且不等于noData，而且Array[i-1][j]等于noData，则以索引为[i,j]、[i,j+1]和[i-1,j+1]的像元的中心点构建三角网，每个顶点的高程即对应像元的高程值，计算三角形面积为Stra，Sall＝Sall+Stra，则进入下一步；
否则，进入步骤S3.5.3；
S3.5.2判断该三角形三个顶点是否为边缘像元，如果是，则计算其面积：
如果Dictionary中存在像元[i，j]，则Dictionary[i,j]＝Dictionary[i,j]+S0/8；
如果Dictionary中存在像元[i，j+1]，则Dictionary[i,j+1]＝Dictionary[i,j+1]+2*S0/8；
如果Dictionary中存在像元[i-1，j+1]，则Dictionary[i-1，j+1]＝Dictionary[i-1，j+1]+S0/8；
S3.5.3，判断像元Array[i-1][j]、Array[i-1][j+1]和Array[i][j+1]：
如果Array[i-1][j]和Array[i-1][j+1]存在且不等于noData，而且Array[i][j+1]不存在或者等于noData，则以索引为[i,j]、[i-1,j]和[i-1,j+1]的像元的中心点构建三角网，每个顶点的高程即对应像元的高程值，计算三角形面积为Stra，Sall＝Sall+Stra，则进入下一步；
否则，进入步骤S3.6
S3.5.4判断该三角形三个顶点是否为边缘像元，如果是，则计算其面积：
如果Dictionary中存在像元[i，j]，则Dictionary[i,j]＝Dictionary[i,j]+S0/8；
如果Dictionary中存在像元[i-1,j]，则Dictionary[i-1,j]＝Dictionary[i-1,j]+2*S0/8；
如果Dictionary中存在像元[i-1,j+1]，则Dictionary[i-1,j+1]＝Dictionary[i-1，j+1]+S0/8；
S3.6，判断第i+1行数据是否存在：
如果i+1>＝height，则退出计算；
如果i+1<height，则进入下一步；
S3.6.1，判断像元Array[i][j+1]和Array[i+1][j+1]：
如果Array[i][j+1]和Array[i+1][j+1]存在且不等于noData，则以索引为[i,j]、[i,j+1]和[i+1,j+1]的像元的中心点构建三角网，每个顶点的高程即对应像元的高程值，计算三角形面积为Stra，Sall＝Sall+Stra，则进入下一步；
否则，进入步骤S3.6.3；
S3.6.2，判断该三角形三个顶点是否为边缘像元，如果是，则计算其面积：
如果Dictionary中存在像元[i，j]，则Dictionary[i,j]＝Dictionary[i,j]+S0/8；
如果Dictionary中存在像元[i,j+1]，则Dictionary[i,j+1]＝Dictionary[i,j+1]+2*S0/8；
如果Dictionary中存在像元[i+1,j+1]，则Dictionary[i+1,j+1]＝Dictionary[i+1,j+1]+S0/8；
S3.6.3，判断像元Array[i+1][j]和Array[i+1][j+1]：
如果Array[i+1][j]和Array[i+1][j+1]存在且不等于noData，则以索引为[i,j]、[i+1,j]和[i+1,j+1]的像元的中心点构建三角网，每个顶点的高程即对应像元的高程值，计算三角形面积为Stra，Sall＝Sall+Stra，则进入下一步；
否则，进入步骤S3.7；
S3.6.4，判断该三角形三个顶点是否为边缘像元，如果是，则计算其面积：
如果Dictionary中存在像元[i，j]，则Dictionary[i,j]＝Dictionary[i,j]+S0/8；
如果Dictionary中存在像元[i+1,j]，则Dictionary[i+1,j]＝Dictionary[i+1,j]+2*S0/8；
如果Dictionary中存在像元[i+1,j+1]，则Dictionary[i+1,j+1]＝Dictionary[i+1,j+1]+S0/8；
S3.7，j＝j+1,返回步骤S3.3；
S3.8，对边界像元面积字典Dictionary中所有项求和得到Sdic，并加入到Sall，即得到整个栅格数据的地表表面面积Sall＝Sall+Sdic。
2.根据权利要求1所述的基于大区域栅格数字高程模型的地表表面面积计算方法，其特征是，所述的原始栅格数字高程模型的表面面积通过累加所有栅格数据的表面面积得到。
3.基于大区域栅格数字高程模型的地表表面面积计算方法，其特征是：包括以下步骤：
S1，对原始栅格数字高程模型进行分块处理，生成若干个栅格数据；
S2，查找每个栅格数据的边界像元并计算边界像元面积；
S3，对每个栅格数据构建三角网并计算三角网面积；
S4，计算每个栅格数据的地表表面面积；
S5，累加所有栅格数据的表面面积即得到原始栅格数字高程模型的表面面积；
在步骤S1中，对原始栅格数字高程模型进行格网化分割，生成若干个大小相同的栅格数据；
所述步骤S2的具体步骤如下：
S2.1，将每个栅格数据转化为高程数组Array[height][width]，width为栅格数据的水平方向像元个数，height为栅格数据的垂直方向像元个数，Array[x][y]表示栅格数据第x行第y列像元的高程值，如果像元[x,y]为空值，则取Array[x][y]为浮点数的最小值，记作noData；通过x,y和栅格数据左上角坐标计算出第x行第y列像元中心点的地理坐标；
记计算像元8邻域像元个数的函数为NeighbourCount，NeighbourCount(x,y)为计算位置[x-1,y-1]、[x,y-1]、[x+1,y-1]、[x-1,y]、[x+1,y]、[x-1,y+1]、[x,y+1]和[x+1,y+1]存在非空像元的个数，NeighbourCount(x,y)<8则代表像元(x,y)为边界像元；
声明边界像元面积字典Dictionary，键为边界点坐标，值为边界像元面积；Dictionary[x,y]＝S表示坐标为[x,y]的边界像元的面积为S；
获取一个像元的平面面积，记为S0；
设当前行i＝0,当前列j＝0：
S2.2，比较i和height的值：
如果i＝height，则退出计算；
如果i<height，则进入下一步；
S2.3，比较j和width的值：
如果j＝width，i＝i+1，j＝0，返回步骤S2.2；
如果j<width，则进入下一步；
S2.4，判断像元(i,j)是否为边界像元：
如果NeighbourCount(i,j)＝8，则j＝j+1，返回步骤S2.3；
如果NeighbourCount(i,j)<8，则进入下一步；
S2.5，在边界像元面积字典中加入像元[i,j]，并初始化面积为0，即Dictionary[i,j]＝0；
S2.6，j＝j+1,返回步骤S2.3；
所述步骤S3的具体步骤如下：
S3.1，获取步骤S2.1高程数组Array[height][width]，像元空值noData和边界像元面积字典Dictionary；
构建三角网：提取每个像元的中心点作为顶点，每个顶点与相邻的8个像元中心点中的2个形成三角形；
计算边界像元面积：每个边界像元按照对角线和对边中位线分成8份，每一份面积为S0/8，计算在构建三角网之后剩余的份数，然后计算剩余份数的平面面积，即得到该边界像元的面积；
设当前行i＝0，当前列j＝0，栅格数据的地表表面面积为Sall；
S3.2，比较i和height的值：
如果i＝height，则进入步骤S3.8；
如果i<height，则进入下一步；
S3.3，比较j和width的值：
如果j＝width，i＝i+1，j＝0，返回步骤S3.2；
如果j<width，则进入下一步；
S3.4，比较Array[i][j]和noData的值：
如果Array[i][j]＝noData，j＝j+1，返回步骤S3.3；
如果Array[i][j]>noData，则进入下一步；
S3.5，判断第i-1行数据是否存在：
如果i-1<0，则进入步骤S3.6；
如果i-1>＝0，则进入下一步；
S3.5.1，判断像元Array[i-1][j+1]、Array[i][j+1]和Array[i-1][j]：
如果Array[i-1][j+1]和Array[i][j+1]存在且不等于noData，而且Array[i-1][j]等于noData，则以索引为[i,j]、[i,j+1]和[i-1,j+1]的像元的中心点构建三角网，每个顶点的高程即对应像元的高程值，计算三角形面积为Stra，Sall＝Sall+Stra，则进入下一步；
否则，进入步骤S3.5.3；
S3.5.2判断该三角形三个顶点是否为边缘像元，如果是，则计算其面积：
如果Dictionary中存在像元[i，j]，则Dictionary[i,j]＝Dictionary[i,j]+S0/8；
如果Dictionary中存在像元[i，j+1]，则Dictionary[i,j+1]＝Dictionary[i,j+1]+2*S0/8；
如果Dictionary中存在像元[i-1，j+1]，则Dictionary[i-1，j+1]＝Dictionary[i-1，j+1]+S0/8；
S3.5.3，判断像元Array[i-1][j]、Array[i-1][j+1]和Array[i][j+1]：
如果Array[i-1][j]和Array[i-1][j+1]存在且不等于noData，而且Array[i][j+1]不存在或者等于noData，则以索引为[i,j]、[i-1,j]和[i-1,j+1]的像元的中心点构建三角网，每个顶点的高程即对应像元的高程值，计算三角形面积为Stra，Sall＝Sall+Stra，则进入下一步；
否则，进入步骤S3.6
S3.5.4判断该三角形三个顶点是否为边缘像元，如果是，则计算其面积：
如果Dictionary中存在像元[i，j]，则Dictionary[i,j]＝Dictionary[i,j]+S0/8；
如果Dictionary中存在像元[i-1,j]，则Dictionary[i-1,j]＝Dictionary[i-1,j]+2*S0/8；
如果Dictionary中存在像元[i-1,j+1]，则Dictionary[i-1,j+1]＝Dictionary[i-1，j+1]+S0/8；
S3.6，判断第i+1行数据是否存在：
如果i+1>＝height，则退出计算；
如果i+1<height，则进入下一步；
S3.6.1，判断像元Array[i][j+1]和Array[i+1][j+1]：
如果Array[i][j+1]和Array[i+1][j+1]存在且不等于noData，则以索引为[i,j]、[i,j+1]和[i+1,j+1]的像元的中心点构建三角网，每个顶点的高程即对应像元的高程值，计算三角形面积为Stra，Sall＝Sall+Stra，则进入下一步；
否则，进入步骤S3.6.3；
S3.6.2，判断该三角形三个顶点是否为边缘像元，如果是，则计算其面积：
如果Dictionary中存在像元[i，j]，则Dictionary[i,j]＝Dictionary[i,j]+S0/8；
如果Dictionary中存在像元[i,j+1]，则Dictionary[i,j+1]＝Dictionary[i,j+1]+2*S0/8；
如果Dictionary中存在像元[i+1,j+1]，则Dictionary[i+1,j+1]＝Dictionary[i+1,j+1]+S0/8；
S3.6.3，判断像元Array[i+1][j]和Array[i+1][j+1]：
如果Array[i+1][j]和Array[i+1][j+1]存在且不等于noData，则以索引为[i,j]、[i+1,j]和[i+1,j+1]的像元的中心点构建三角网，每个顶点的高程即对应像元的高程值，计算三角形面积为Stra，Sall＝Sall+Stra，则进入下一步；
否则，进入步骤S3.7；
S3.6.4，判断该三角形三个顶点是否为边缘像元，如果是，则计算其面积：
如果Dictionary中存在像元[i，j]，则Dictionary[i,j]＝Dictionary[i,j]+S0/8；
如果Dictionary中存在像元[i+1,j]，则Dictionary[i+1,j]＝Dictionary[i+1,j]+2*S0/8；
如果Dictionary中存在像元[i+1,j+1]，则Dictionary[i+1,j+1]＝Dictionary[i+1,j+1]+S0/8；
S3.7，j＝j+1,返回步骤S3.3；
S3.8，对边界像元面积字典Dictionary中所有项求和得到Sdic，并加入到Sall，即得到整个栅格数据的地表表面面积Sall＝Sall+Sdic。
4.根据权利要求3所述的基于大区域栅格数字高程模型的地表表面面积计算方法，其特征是，在步骤S4中，每个栅格数据的地表表面面积为每块栅格数据的三角网面积和边界像元面积之和。
CN201710138021.3A 2017-03-09 2017-03-09 基于大区域栅格数字高程模型的地表表面面积计算方法 Active CN106910218B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201710138021.3A CN106910218B (zh) 2017-03-09 2017-03-09 基于大区域栅格数字高程模型的地表表面面积计算方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201710138021.3A CN106910218B (zh) 2017-03-09 2017-03-09 基于大区域栅格数字高程模型的地表表面面积计算方法

Publications (2)

Publication Number Publication Date
CN106910218A CN106910218A (zh) 2017-06-30
CN106910218B true CN106910218B (zh) 2019-07-09

Family

ID=59187808

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201710138021.3A Active CN106910218B (zh) 2017-03-09 2017-03-09 基于大区域栅格数字高程模型的地表表面面积计算方法

Country Status (1)

Country Link
CN (1) CN106910218B (zh)

Families Citing this family (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110378951A (zh) * 2018-09-10 2019-10-25 天津京东深拓机器人科技有限公司 确定堆垛物品支撑面积的方法和装置
CN110298778B (zh) * 2019-06-10 2022-12-27 东南大学 一种基于栅格数据和面积比修正的交通小区人口统计方法
CN112686945A (zh) * 2021-01-11 2021-04-20 阿坝州自然资源与科技信息研究所 一种基于Web的三维地形表面积测量方法
CN116645482B (zh) * 2023-06-01 2024-03-22 中国铁路设计集团有限公司 一种大范围地面连续高程提取方法

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102842104A (zh) * 2012-07-16 2012-12-26 长江水利委员会长江科学院 面向海量dem数据的高精度河道洪水淹没区生成方法
CN102902844A (zh) * 2012-09-03 2013-01-30 南京师范大学 基于大数据量dem数据的子流域划分方法
CN103256914A (zh) * 2012-12-27 2013-08-21 北京地拓科技发展有限公司 一种基于dem计算淤地坝淹没面积的方法及***
CN106123812A (zh) * 2016-08-14 2016-11-16 覃泽林 基于遥感影像获取起伏地表甘蔗种植面积的方法及装置

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102842104A (zh) * 2012-07-16 2012-12-26 长江水利委员会长江科学院 面向海量dem数据的高精度河道洪水淹没区生成方法
CN102902844A (zh) * 2012-09-03 2013-01-30 南京师范大学 基于大数据量dem数据的子流域划分方法
CN103256914A (zh) * 2012-12-27 2013-08-21 北京地拓科技发展有限公司 一种基于dem计算淤地坝淹没面积的方法及***
CN106123812A (zh) * 2016-08-14 2016-11-16 覃泽林 基于遥感影像获取起伏地表甘蔗种植面积的方法及装置

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
基于地理坐标的地理国情统计单元表面面积精确计算;谢成磊 等;《遥感信息》;20140831;第48页左栏第2段，第49页左栏第3段-右栏第2段及图3

Also Published As

Publication number Publication date
CN106910218A (zh) 2017-06-30

Similar Documents

Publication Publication Date Title
CN106910218B (zh) 基于大区域栅格数字高程模型的地表表面面积计算方法
Arnold A new approach for dealing with depressions in digital elevation models when calculating flow accumulation values
US7983474B2 (en) Geospatial modeling system and related method using multiple sources of geographic information
CN104952107A (zh) 基于车载LiDAR点云数据的桥梁三维重建方法
Ermolaev et al. Automated construction of the boundaries of basin geosystems for the Volga Federal District
CN109241846A (zh) 遥感影像的时空变化估测方法、装置与存储介质
da Paz et al. Errors in river lengths derived from raster digital elevation models
CN104463970B (zh) 一种基于遥感影像确定城市三维重心的方法及其应用
CN106408581B (zh) 一种快速的三维点云直线提取方法
CN103559506B (zh) 一种基于矢量边界的亚像元制图方法
CN102122395A (zh) 一种保持地形特征的自适应尺度dem建模方法
CN108919295A (zh) 机载LiDAR点云道路信息提取方法及装置
Wang et al. The isotropic organization of DEM structure and extraction of valley lines using hexagonal grid
US7778808B2 (en) Geospatial modeling system providing data thinning of geospatial data points and related methods
Gumus et al. Comparison of spatial interpolation methods and multi-layer neural networks for different point distributions on a digital elevation model
Brenner Scalable estimation of precision maps in a mapreduce framework
CN110207670A (zh) 一种使用二维图像获取人工林森林冠幅参数的方法
JP4619504B2 (ja) ３次元デジタル地図作成装置
Adreani et al. Digital twin framework for smart city solutions
Sun et al. Building outline extraction from aerial imagery and digital surface model with a frame field learning framework
Zhao et al. A new algorithm for the automatic extraction of valley floor width
de Smith Distance transforms as a new tool in spatial analysis, urban planning, and GIS
Adreani et al. A photorealistic 3D city modeling framework for smart city digital twin
Zhu et al. Uncertainty of slope length derived from digital elevation models of the Loess Plateau, China
CN114969944B (zh) 一种高精度道路dem构建方法

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