CN110288645B - 一种基于坡向约束的地形表面面积计算方法 - Google Patents
一种基于坡向约束的地形表面面积计算方法 Download PDFInfo
- Publication number
- CN110288645B CN110288645B CN201910354072.9A CN201910354072A CN110288645B CN 110288645 B CN110288645 B CN 110288645B CN 201910354072 A CN201910354072 A CN 201910354072A CN 110288645 B CN110288645 B CN 110288645B
- Authority
- CN
- China
- Prior art keywords
- data set
- point
- slope
- surface area
- pixel
- 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
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01C—MEASURING DISTANCES, LEVELS OR BEARINGS; SURVEYING; NAVIGATION; GYROSCOPIC INSTRUMENTS; PHOTOGRAMMETRY OR VIDEOGRAMMETRY
- G01C5/00—Measuring height; Measuring distances transverse to line of sight; Levelling between separated points; Surveyors' levels
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/60—Analysis of geometric attributes
- G06T7/62—Analysis of geometric attributes of area, perimeter, diameter or volume
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/30—Subject of image; Context of image processing
- G06T2207/30181—Earth observation
- G06T2207/30184—Infrastructure
Landscapes
- Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Geometry (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Theoretical Computer Science (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Processing Or Creating Images (AREA)
- Instructional Devices (AREA)
Abstract
本发明提供了一种基于坡向约束的地形表面面积计算方法,针对现有技术无法有效地直接存储地形信息数据,且忽视地形要素对于地形地貌的骨架控制作用等的问题,通过地形分析提取地形地貌要素,挖掘隐含在DEM数据中的地形地貌信息,通过地形因子辅助的作用提高了地形表面面积建模精度。
Description
技术领域
本发明属于地表面积计算技术领域,具体涉及一种基于坡向约束的地形表面面积计算方法。
背景技术
地形地貌是地球自然过程和人类活动共同作用的结果,是地表最复杂的现象之一。其形态千变万化但可以分解一系列规律的点、线、面要素。这些点、线、面要素决定了地貌形态的规模和类型,通过地形分析提取要素是地形地貌分析的基础。作用于地形地貌内外营力不同,所形成地形地貌的类型、尺度和规模会出现明显差异;经过长期的地理理论和实践研究,一般认为内营力是形成星体、巨型和大型地貌的主因,而外营力是形成中小型和微型地貌的主因。数字地形地貌分析是对DEM数据的深度信息的提取,由于DEM数据仅记录离散的高程点、线数据,无法有效地直接存储地形信息数据,因此有必要通过一系列提取方法挖掘隐含在DEM数据中的地形地貌信息。表面面积是众多地形因子中非常重要的一个,也是众多其他复杂地形因子的基础因子。随着遥感技术的提高,基于规则格网DEM格式存储的高精度高程数据得到广泛推广,而传统规则格网DEM数据的表面面积提取方法通过三角构网、曲面拟合等方法实现表面面积计算,忽视了地形要素对于地形地貌的骨架控制作用。
发明内容
本发明要解决的技术问题是:提供一种基于坡向约束的地形表面面积计算方法,用于通过地形因子辅助的作用提高地形表面面积的建模精度。
本发明为解决上述技术问题所采取的技术方案为:一种基于坡向约束的地形表面面积计算方法,包括以下步骤:
步骤S1:收集数字高程模型数据集A;
步骤S2:根据步骤S1得到的数据集A计算坡向数据集B;
步骤S3:根据数据集A创建表面面积数据集C;
步骤S4:以数据集A和步骤S2得到的数据集B作为输入数据,数据集C为输出参数,以像元矩阵作为滑动窗口,对数据集A、数据集B和数据集C进行遍历操作;
步骤S5:步骤S4创建的滑动窗口每次变换之后,以窗口内数据集A和数据集B作为输入参数,计算窗口中心像元表面面积,将结果计入数据集C;
步骤S6:累积数据集C中所有像元的数值之和作为所求的表面面积。
按上述方案,所述的步骤S2中,坡向计算采用基于数值积分的二次曲面拟合曲面法,设Aspect表示像元坡向,Slopesn表示像元在南北方向坡度,Slopewe在东西方向的坡度,则每个像元的坡度和坡向的计算公式为:
Aspect=Slopesn/Slopewe。
进一步的,所述的步骤S3中,具体步骤为:
步骤S31:分别设数据集A的行数为Rows、列数为Cols,设数据集A的像元尺寸为x-size和y-size,设数据集A的起始点的横坐标为x-tl、纵坐标为y-tl;
步骤S32:以Rows、Cols、x-size、y-size、x-tl、y-tl作为栅格参数,创建数据类型为浮点型的栅格数据集,记为C。
进一步的,所述的步骤S4中,具体步骤为:
步骤S41:设像元矩阵的行列数为3×3,设正整数Rc的取值范围为(1,Rows-1),遍历Rc的所有取值,每次递增1;
步骤S42:针对步骤S41中Rc的每一次取值,设正整数Cc的取值范围为(1,Cols-1),遍历Cc的所有取值,每次递增1;
步骤S43:针对步骤S42中Cc的每一次取值,分别读取数据集A、数据集B中行列号为(Rc-1,Cc-1)、(Rc-1,Cc)、(Rc-1,Cc+1)、(Rc,Cc-1)、(Rc,Cc)、(Rc,Cc+1)、(Rc+1,Cc-1)、(Rc+1,Cc)、(Rc+1,Cc+1)的像元值,作为当前窗口内数据集A和数据集B的读取值。
进一步的,所述的步骤S5中,具体步骤为:
步骤S51:将窗口内数据集A的像元中心点分别记为a、b、c、d、e、f、g、h、i;
步骤S52:将窗口内数据集B的中心像元值记为θ;
步骤S53:过e点以θ为方向角做直线l,l与a点和g点的连线相交于j点,与c点和i点的连线相交于k点;
步骤S54:分别连接b点和j点、h点和k点、b点和f点、d点和h点,将步骤S51得到的中心点围成的矩形区域划分为10个三角形区域;
步骤S55:计算步骤S54得到的各三角形面积之和再除以4,获得窗口中心点的表面积栅格像元值,将结果计入数据集C。
本发明的有益效果为:
1.本发明的一种基于坡向约束的地形表面面积计算方法通过地形因子辅助提高了地形表面面积的建模精度。
2.本发明能够挖掘隐含在DEM数据中的地形地貌信息。
附图说明
图1是本发明实施例的流程图。
图2是本发明实施例的坡向计算示意图。
图3是本发明实施例的坡向约束表面面积计算示意图。
图4是本发明实施例的研究区表面积结果图。
具体实施方式
下面结合附图和具体实施方式对本发明作进一步详细的说明。
以神农架林区作为研究,参见图1至图4,本发明的一种基于坡向约束的地形表面面积计算方法,包括以下步骤:
步骤S1:收集研究区数字高程模型数据集A。
步骤S2:根据步骤S1得到的数字高程模型数据集A计算研究区坡向数据集B;坡向计算采用基于数值积分的二次曲面拟合曲面法,设Aspect表示像元坡向,Slopesn表示像元在南北方向坡度,Slopewe在东西方向的坡度,则每个像元的坡度和坡向的计算公式为:
Aspect=Slopesn/Slopewe。
步骤S3:依据数据集A的坐标系统信息创建表面面积数据集C:
步骤S31:获取数据集A的行数为Rows=399、列数为Cols=474,获取数据集A的像元尺寸分别记为x-size=30m和y-size=30m,获取数据集A的起始点的横坐标为x-tl=110.3877°、纵坐标为y-tl=31.6819°;
步骤S32:以Rows、Cols、x-size、y-size、x-tl、y-tl作为栅格参数,创建数据类型为浮点型的栅格数据集,记为C。
步骤S4:以数据集A和步骤S2得到的数据集B作为输入数据,数据集C为输出参数,以像元矩阵作为滑动窗口,对数据集A、数据集B和数据集C进行遍历操作:
步骤S41:设像元矩阵的行列数为3×3,根据步骤S3设置的数据集A的行列数,设Rc为正整数,取值范围为(1,Rows-1=398),遍历Rc的所有取值,每次递增1;
步骤S42:针对步骤S41中Rc的每一次取值,设Cc为正整数,取值范围为(1,Cols-1=473),遍历Cc的所有取值,每次递增1;
步骤S43:针对步骤S42中Cc的每一次取值,分别读取数据A、B中行列号为(Rc-1,Cc-1)、(Rc-1,Cc)、(Rc-1,Cc+1)、(Rc,Cc-1)、(Rc,Cc)、(Rc,Cc+1)、(Rc+1,Cc-1)、(Rc+1,Cc)、(Rc+1,Cc+1)的像元值,作为当前窗口数据集A和数据集B的读取值。
步骤S5:步骤S4创建的滑动窗口每次变换之后,以窗口内数据集A和数据集B作为输入参数,计算窗口中心像元表面面积,将结果计入数据集C:
步骤S51:获取窗口内数据集A的像元中心点分别记为a、b、c、d、e、f、g、h、i;
步骤S52:获取窗口内数据B的中心像元值记为θ;
步骤S53:过e点以θ为方向角做直线l,l与a点和g点的连线相交于j点,与c点和i点的连线相交于k点;
步骤S54:分别连接b点和j点、h点和k点、b点和f点、d点和h点,将步骤S51得到的中心点围成的矩形区域划分为10个三角形区域;
步骤S55:计算步骤S54得到的各三角形面积之和再除以4,获得窗口中心点的表面积栅格像元值,将结果计入数据集C。
步骤6:累积数据集C中所有像元数值之和作为研究区表面面积。
综上所述,本发明的一种基于坡向约束的地形表面面积计算方法针对现有技术无法有效地直接存储地形信息数据,且忽视地形要素对于地形地貌的骨架控制作用等的问题,通过地形分析提取地形地貌要素,挖掘隐含在DEM数据中的地形地貌信息,通过地形因子辅助的作用提高了地形表面面积建模精度。
以上实施例仅用于说明本发明的设计思想和特点,其目的在于使本领域内的技术人员能够了解本发明的内容并据以实施,本发明的保护范围不限于上述实施例。所以,凡依据本发明所揭示的原理、设计思路所作的等同变化或修饰,均在本发明的保护范围之内。
Claims (4)
1.一种基于坡向约束的地形表面面积计算方法,其特征在于:包括以下步骤:
步骤S1:收集数字高程模型数据集A;
步骤S2:根据步骤S1得到的数据集A计算坡向数据集B;
所述的步骤S2中,坡向计算采用基于数值积分的二次曲面拟合曲面法,设Aspect表示像元坡向,Slopesn表示像元在南北方向坡度,Slopewe在东西方向的坡度,则每个像元的坡度和坡向的计算公式为:
Aspect=Slopesn/Slopewe;
步骤S3:根据数据集A创建表面面积数据集C;
步骤S4:以数据集A和步骤S2得到的数据集B作为输入数据,数据集C为输出参数,以像元矩阵作为滑动窗口,对数据集A、数据集B和数据集C进行遍历操作;
步骤S5:步骤S4创建的滑动窗口每次变换之后,以窗口内数据集A和数据集B作为输入参数,计算窗口中心像元表面面积,将结果计入数据集C;
步骤S6:累积数据集C中所有像元的数值之和作为所求的表面面积。
2.根据权利要求1所述的一种基于坡向约束的地形表面面积计算方法,其特征在于:所述的步骤S3中,具体步骤为:
步骤S31:分别设数据集A的行数为Rows、列数为Cols,设数据集A的像元尺寸为x-size和y-size,设数据集A的起始点的横坐标为x-tl、纵坐标为y-tl;
步骤S32:以Rows、Cols、x-size、y-size、x-tl、y-tl作为栅格参数,创建数据类型为浮点型的栅格数据集,记为C。
3.根据权利要求2所述的一种基于坡向约束的地形表面面积计算方法,其特征在于:所述的步骤S4中,具体步骤为:
步骤S41:设像元矩阵的行列数为3×3,设正整数Rc的取值范围为(1,Rows-1),遍历Rc的所有取值,每次递增1;
步骤S42:针对步骤S41中Rc的每一次取值,设正整数Cc的取值范围为(1,Cols-1),遍历Cc的所有取值,每次递增1;
步骤S43:针对步骤S42中Cc的每一次取值,分别读取数据集A、数据集B中行列号为(Rc-1,Cc-1)、(Rc-1,Cc)、(Rc-1,Cc+1)、(Rc,Cc-1)、(Rc,Cc)、(Rc,Cc+1)、(Rc+1,Cc-1)、(Rc+1,Cc)、(Rc+1,Cc+1)的像元值,作为当前窗口内数据集A和数据集B的读取值。
4.根据权利要求3所述的一种基于坡向约束的地形表面面积计算方法,其特征在于:所述的步骤S5中,具体步骤为:
步骤S51:将窗口内数据集A的像元中心点分别记为a、b、c、d、e、f、g、h、i;
步骤S52:将窗口内数据集B的中心像元值记为θ;
步骤S53:过e点以θ为方向角做直线l,l与a点和g点的连线相交于j点,与c点和i点的连线相交于k点;
步骤S54:分别连接b点和j点、h点和k点、b点和f点、d点和h点,将步骤S51得到的中心点围成的矩形区域划分为10个三角形区域;
步骤S55:计算步骤S54得到的各三角形面积之和再除以4,获得窗口中心点的表面积栅格像元值,将结果计入数据集C。
Priority Applications (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910354072.9A CN110288645B (zh) | 2019-04-29 | 2019-04-29 | 一种基于坡向约束的地形表面面积计算方法 |
AU2020100543A AU2020100543A4 (en) | 2019-04-29 | 2020-04-09 | Method for calculating surface area of terrain based on aspect constraints |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910354072.9A CN110288645B (zh) | 2019-04-29 | 2019-04-29 | 一种基于坡向约束的地形表面面积计算方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN110288645A CN110288645A (zh) | 2019-09-27 |
CN110288645B true CN110288645B (zh) | 2022-10-04 |
Family
ID=68001718
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910354072.9A Active CN110288645B (zh) | 2019-04-29 | 2019-04-29 | 一种基于坡向约束的地形表面面积计算方法 |
Country Status (2)
Country | Link |
---|---|
CN (1) | CN110288645B (zh) |
AU (1) | AU2020100543A4 (zh) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113379828B (zh) * | 2021-06-04 | 2023-02-10 | 西北农林科技大学 | 一种融合地表形态特征的坡长提取方法 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US7305983B1 (en) * | 2003-07-08 | 2007-12-11 | University Of Hawai'i | Assessment of solar energy potential on existing buildings in a region |
CN103592676A (zh) * | 2013-10-24 | 2014-02-19 | 中国石油天然气集团公司 | 一种基于地形因子的炮点偏移方法 |
CN107239664A (zh) * | 2017-06-08 | 2017-10-10 | 国家基础地理信息中心 | 一种地表覆盖表面积栅格化计算方法及装置 |
CN109492194A (zh) * | 2018-12-29 | 2019-03-19 | 南京泛在地理信息产业研究院有限公司 | 一种基于数学向量几何的dem二阶地形因子计算方法 |
-
2019
- 2019-04-29 CN CN201910354072.9A patent/CN110288645B/zh active Active
-
2020
- 2020-04-09 AU AU2020100543A patent/AU2020100543A4/en active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US7305983B1 (en) * | 2003-07-08 | 2007-12-11 | University Of Hawai'i | Assessment of solar energy potential on existing buildings in a region |
CN103592676A (zh) * | 2013-10-24 | 2014-02-19 | 中国石油天然气集团公司 | 一种基于地形因子的炮点偏移方法 |
CN107239664A (zh) * | 2017-06-08 | 2017-10-10 | 国家基础地理信息中心 | 一种地表覆盖表面积栅格化计算方法及装置 |
CN109492194A (zh) * | 2018-12-29 | 2019-03-19 | 南京泛在地理信息产业研究院有限公司 | 一种基于数学向量几何的dem二阶地形因子计算方法 |
Non-Patent Citations (2)
Title |
---|
Quantitative Analysis of Error in Extracting Slope Length from Grid-Based Digital Elevation Model;Bei Jin;《2010 International Conference on Multimedia Technology》;20101211;全文 * |
基于现代测量平差的InSAR三维形变估计理论与方法;胡俊;《中国博士学位论文全文数据库》;20150215;第6章 * |
Also Published As
Publication number | Publication date |
---|---|
AU2020100543A4 (en) | 2020-05-28 |
CN110288645A (zh) | 2019-09-27 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN108764255B (zh) | 一种冬小麦种植信息的提取方法 | |
CN110415265B (zh) | 基于无人机高精度dem数据的梯田自动提取方法 | |
CN110348324B (zh) | 一种基于遥感大数据的洪涝实时淹没分析方法及系统 | |
CN114387516B (zh) | 一种针对复杂地形环境下中小田块的单季稻sar识别方法 | |
Wang et al. | A Quantitative Study of Gully Erosion Based on Object‐Oriented Analysis Techniques: A Case Study in Beiyanzikou Catchment of Qixia, Shandong, China | |
Sun et al. | Hydrologic changes of Aral Sea: a reveal by the combination of radar altimeter data and optical images | |
CN111091079A (zh) | 基于tls的高寒脆弱区植被优势单株结构参数测定方法 | |
CN106846246B (zh) | 一种基于对象的遥感影像超分辨率制图方法 | |
CN110288645B (zh) | 一种基于坡向约束的地形表面面积计算方法 | |
Dinis et al. | Hierarchical object-based classification of dense urban areas by integrating high spatial resolution satellite images and lidar elevation data | |
CN117475314B (zh) | 一种地质灾害隐患立体识别方法、系统及介质 | |
Greco et al. | Temporal mapping of riparian landscape change on the Sacramento River, miles 196–218, California, USA | |
Faid et al. | Monitoring land-use change-associated land development using multitemporal Landsat data and geoinformatics in Kom Ombo area, South Egypt | |
Liu et al. | Vegetation mapping for regional ecological research and management: a case of the Loess Plateau in China | |
CN102663761A (zh) | 用于影像地图的线状矢量与遥感影像自动配准方法 | |
Zhou et al. | Modelling spatio-temporal pattern of landuse change using multi-temporal remotely sensed imagery | |
Na et al. | Mapping artificial terraces from image matching point cloud in Loess Plateau of China | |
CN117576485B (zh) | 融合正负地形和光照晕渲的沙丘提取方法 | |
CN117274844B (zh) | 利用无人机遥感影像的大田花生出苗情况快速提取方法 | |
CN112729145B (zh) | 一种确定滑坡位移与库区水位波动关系的方法 | |
CN117315471B (zh) | 一种基于遥感影像与机器学习的梯田识别方法 | |
Brock et al. | Characterizing surface water use for irrigation in the Mississippi Delta | |
Lidzhegu | A classification of large wetlands in Africa’s elevated drylands based on their formation, structure and hydrological functioning using Earth Observation and GIS (Unpublished PhD thesis) | |
Geyer et al. | A geo-archaeological approach to the study of hydro-agricultural systems in arid areas of Western Syria | |
Rademann et al. | CHARACTERIZATION AND ANALYSIS OF GULLY EROSION IN SOUTHERN BRAZIL WITH THE ASSISTANCE OF UNMANNED AERIAL VEHICLE |
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 |