CN107564100A - 一种等高线生成平滑实体的方法 - Google Patents

一种等高线生成平滑实体的方法 Download PDF

Info

Publication number
CN107564100A
CN107564100A CN201710798962.XA CN201710798962A CN107564100A CN 107564100 A CN107564100 A CN 107564100A CN 201710798962 A CN201710798962 A CN 201710798962A CN 107564100 A CN107564100 A CN 107564100A
Authority
CN
China
Prior art keywords
contour
point
triangle
smooth
delaunay triangulation
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
CN201710798962.XA
Other languages
English (en)
Other versions
CN107564100B (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.)
Guodian United Power Technology Co Ltd
Original Assignee
Guodian United Power Technology 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 Guodian United Power Technology Co Ltd filed Critical Guodian United Power Technology Co Ltd
Priority to CN201710798962.XA priority Critical patent/CN107564100B/zh
Publication of CN107564100A publication Critical patent/CN107564100A/zh
Application granted granted Critical
Publication of CN107564100B publication Critical patent/CN107564100B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Landscapes

  • Processing Or Creating Images (AREA)

Abstract

本发明公开了一种等高线生成平滑实体的方法，包括步骤：(1)检验并纠错待使用的等高线；(2)以Delaunay三角剖分算法将所述等高线转换成Delaunay三角网；(3)查询所述Delaunay三角网中的钝角三角形，对所述钝角三角形采用Grid格点加密，使其转换成Grid矩阵网格，则生成平滑实体。其中步骤(3)的具体步骤为：遍历所有三角形，建立索引地理编号；当三角形为钝角三角形且钝角大于120°时，根据索引地理编号对所述钝角三角形和周边进行Grid格点按照行列矩阵加密；删除所述钝角三角形，即得平滑实体。本发明生成的实体表面为局部光滑的曲面，且能提供更为明显的地表特征，大大提高了该平滑实体高性能集群计算的资源利用率和结果收敛率。

Description

一种等高线生成平滑实体的方法
技术领域
本发明涉及等高线技术领域，特别是涉及一种等高线生成平滑实体的方法。
背景技术
等高线作为重要的空间地理信息和赖以进行地形分析的核心数据，已在测绘、资源与环境，灾害防治、国防等科研及经济领域发挥着巨大的作用。等高线转换成数字高程模型，其数字高程模型质量的优劣是核心问题，然而转换后的数字高程模型能否逼真的表现地表是问题的关键所在。
以Delaunay方式构成的不规则三角网，是常用来表现地形表面的主要方法，其具有整场信息冗余量小，地表网格逼近性强等优势。但该方法存在表面模型单一、表面不光滑、一阶导数不连续的问题。
本发明就是针对其上述所存在的问题，创设一种新的等高线生成平滑实体的方法，使其生成的地表实体更平滑、地标特征更明显。
发明内容
本发明要解决的技术问题是提供一种等高线生成平滑实体的方法，使其生成的地表实体更平滑、地标特征更明显，从而克服现有的生成地表实体方法的不足。
为解决上述技术问题，本发明提供一种等高线生成平滑实体的方法，所述方法包括如下步骤：
(1)检验并纠错待使用的等高线；
(2)以Delaunay三角剖分算法将所述等高线转换成Delaunay三角网；
(3)查询所述Delaunay三角网中的钝角三角形，对所述钝角三角形
采用Grid格点加密，使所述钝角三角形转换成Grid矩阵网格，则生
成平滑实体。
作为本发明的一种改进，所述步骤(1)中，遍历所有待使用的等高线的高程值，去除所述等高线中具有错误信息的等高线。
进一步改进，所述等高线的错误信息包括空值、等高线交叉信息。
进一步改进，所述步骤(2)中，按照空心圆法则生成所述Delaunay三角网的具体方法步骤如下：
A1.以任意等高线拐点为起点，查找距离此点最近的第二点，连线所述起点与第二点；
A2.继续寻找距离所述第二点最近的第三点，连线所述第二点与第三点以及起点与第三点，生成三角形；
A3.以所述步骤A2中生成的三角形任一点作为第二点，继续按照所述步骤A2寻找第三点，依次外推，生成Delaunay三角网。
进一步改进，所述步骤(3)中，对所述钝角三角形采用Grid格点加密的具体步骤如下：
B1.遍历所有Delaunay三角网中的三角形，建立索引地理编号；
B2.当三角形为钝角三角形且钝角大于120°时，根据索引地理编号对所述钝角三角形和周边进行Grid格点按照行列矩阵加密；
B3.删除所述钝角三角形；
B4.所述行列矩阵***按照空心圆法则构建三角网，内部生成Grid网格。
进一步改进，所述步骤(3)中，生成的平滑实体按照STL实体网格格式表示，所述STL实体网格按照三维坐标自动生成x格点、y格点和z格点，则形成集群仿真模拟。
采用这样的设计后，本发明至少具有以下优点：
本发明通过采用Grid格点加密的方法消除Delaunay三角网中过于细长的钝角三角形，以行列矩阵来替代细长的钝角三角形，使生成的实体表面为局部光滑的曲面，且能提供更为明显的地表特征，克服了现有技术中采用不光滑的Delaunay三角网可能产生大量的运算耗费和难于收敛的情况，大大提高了其高性能集群计算的资源利用率和结果收敛率。
附图说明
上述仅是本发明技术方案的概述，为了能够更清楚了解本发明的技术手段，以下结合附图与具体实施方式对本发明作进一步的详细说明。
图1是本发明等高线生成平滑实体的方法的流程示意图。
具体实施方式
参照附图1所示，本实施例等高线生成平滑实体的方法，包括如下步骤：
(1)检验并纠错待使用的等高线
遍历所有等高线高程值，要求去除等高线错误信息，如空值、等高线交叉等，保证每条等高线具有明确的实际物理意义。
(2)以Delaunay三角剖分算法将所述等高线转换成Delaunay三角网；
具体的，按照空心圆法则生成所述Delaunay三角网的具体步骤如下：
A1.以任意等高线拐点为起点，查找距离此点最近的第二点，连线所述起点与第二点；
A2.继续寻找距离所述第二点最近的第三点，连线所述第二点与第三点以及起点与第三点，生成三角形；
A3.以所述步骤A2中生成的三角形任一点作为第二点，继续按照所述步骤A2寻找第三点，依次外推，生成Delaunay三角网。
(3)查询所述Delaunay三角网中的钝角三角形，对所述钝角三角形采用Grid格点加密，使所述钝角三角形转换成Grid矩阵网格，则生成平滑实体。
其中，对所述钝角三角形采用Grid格点加密的具体步骤如下：
B1.遍历所有Delaunay三角网中的三角形，建立索引地理编号；
B2.当三角形为钝角三角形且钝角大于120°时，根据索引地理编号对所述钝角三角形和周边进行Grid格点按照行列矩阵加密；
B3.删除所述钝角三角形；
B4.所述行列矩阵***按照空心圆法则构建三角网，内部生成Grid网格。
还有步骤(3)中，生成的平滑实体按照STL实体网格格式表示，所述STL实体网格按照三维坐标自动生成x格点(水平坐标x方向)、y格点(水平坐标y方向)和z格点(垂向坐标z方向)，则形成集群仿真模拟。
本发明为一种将实体局部采用Grid格点行列对Delaunay三角网自动加密的方法，用于将等高线转化为光滑实体，为解决Delaunay三角网生成的实体表面不平滑，计算效率不高的问题提供了有利条件，使生成的实体表面光滑，且提供的地表特征更为明显，大大提高了其高性能集群计算的资源利用率和结果收敛率。
以上所述，仅是本发明的较佳实施例而已，并非对本发明作任何形式上的限制，本领域技术人员利用上述揭示的技术内容做出些许简单修改、等同变化或修饰，均落在本发明的保护范围内。

Claims (6)

1.一种等高线生成平滑实体的方法，其特征在于，所述方法包括如下步骤：
(1)检验并纠错待使用的等高线；
(2)以Delaunay三角剖分算法将所述等高线转换成Delaunay三角网；
(3)查询所述Delaunay三角网中的钝角三角形，对所述钝角三角形采用Grid格点加密，使所述钝角三角形转换成Grid矩阵网格，则生成平滑实体。
2.根据权利要求1所述的一种等高线生成平滑实体的方法，其特征在于，所述步骤(1)中，遍历所有待使用的等高线的高程值，去除所述等高线中具有错误信息的等高线。
3.根据权利要求2所述的一种等高线生成平滑实体的方法，其特征在于，所述等高线的错误信息包括空值、等高线交叉信息。
4.根据权利要求1所述的一种等高线生成平滑实体的方法，其特征在于，所述步骤(2)中，按照空心圆法则生成所述Delaunay三角网的具体方法步骤如下：
A1.以任意等高线拐点为起点，查找距离此点最近的第二点，连线所述起点与第二点；
A2.继续寻找距离所述第二点最近的第三点，连线所述第二点与第三点以及起点与第三点，生成三角形；
A3.以所述步骤A2中生成的三角形任一点作为第二点，继续按照所述步骤A2寻找第三点，依次外推，生成Delaunay三角网。
5.根据权利要求1所述的一种等高线生成平滑实体的方法，其特征在于，所述步骤(3)中，对所述钝角三角形采用Grid格点加密的具体步骤如下：
B1.遍历所有Delaunay三角网中的三角形，建立索引地理编号；
B2.当三角形为钝角三角形且钝角大于120°时，根据索引地理编号对所述钝角三角形和周边进行Grid格点按照行列矩阵加密；
B3.删除所述钝角三角形；
B4.所述行列矩阵***按照空心圆法则构建三角网，内部生成Grid网格。
6.根据权利要求5所述的一种等高线生成平滑实体的方法，其特征在于，所述步骤(3)中，生成的平滑实体按照STL实体网格格式表示，所述STL实体网格按照三维坐标自动生成x格点、y格点和z格点，则形成集群仿真模拟。
CN201710798962.XA 2017-09-07 2017-09-07 一种等高线生成平滑实体的方法 Active CN107564100B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201710798962.XA CN107564100B (zh) 2017-09-07 2017-09-07 一种等高线生成平滑实体的方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201710798962.XA CN107564100B (zh) 2017-09-07 2017-09-07 一种等高线生成平滑实体的方法

Publications (2)

Publication Number Publication Date
CN107564100A true CN107564100A (zh) 2018-01-09
CN107564100B CN107564100B (zh) 2020-06-19

Family

ID=60978219

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201710798962.XA Active CN107564100B (zh) 2017-09-07 2017-09-07 一种等高线生成平滑实体的方法

Country Status (1)

Country Link
CN (1) CN107564100B (zh)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110489511A (zh) * 2019-08-23 2019-11-22 西南大学 等高线接边高程错误修正方法、***及电子设备和介质
CN114329575A (zh) * 2022-03-15 2022-04-12 武汉峰岭科技有限公司 基于内插山脊线的等高线加密方法、***、介质及设备

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105243688A (zh) * 2015-10-19 2016-01-13 中国科学院自动化研究所 基于重心Voronoi图的非钝角的重新网格化方法
CN105678846A (zh) * 2016-02-22 2016-06-15 武汉华信联创技术工程有限公司 一种实时气象组网雷达数据的三维可视化方法及***

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105243688A (zh) * 2015-10-19 2016-01-13 中国科学院自动化研究所 基于重心Voronoi图的非钝角的重新网格化方法
CN105678846A (zh) * 2016-02-22 2016-06-15 武汉华信联创技术工程有限公司 一种实时气象组网雷达数据的三维可视化方法及***

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
S.H.LO: "Delaunay triangulation of non-uniform point distributions by means of multi-grid insertion", 《FINITE ELEMENTS IN ANALYSIS AND DESIGN》 *
徐道柱，刘海砚: "大量约束边条件下Delaunay三角网的快速生成", 《测绘工程》 *

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110489511A (zh) * 2019-08-23 2019-11-22 西南大学 等高线接边高程错误修正方法、***及电子设备和介质
CN110489511B (zh) * 2019-08-23 2021-07-30 西南大学 等高线接边高程错误修正方法、***及电子设备和介质
CN114329575A (zh) * 2022-03-15 2022-04-12 武汉峰岭科技有限公司 基于内插山脊线的等高线加密方法、***、介质及设备

Also Published As

Publication number Publication date
CN107564100B (zh) 2020-06-19

Similar Documents

Publication Publication Date Title
CN106504326B (zh) 顾及形态精度的地形高程采样点加密方法
CN110276732B (zh) 一种顾及地形特征线要素的山区点云空洞修复方法
CN103871102A (zh) 一种基于高程点和道路轮廓面的道路三维精细建模方法
CN103870845A (zh) 点云聚类去噪过程中新的k值优化方法
CN103500371A (zh) 基于bim的变更控制方法
CN102902844A (zh) 基于大数据量dem数据的子流域划分方法
CN107038308B (zh) 一种基于线性内插的规则格网地形建模方法
CN110544305B (zh) 面向规则格网dem构建的地形陡坎线信息融合方法
CN109636870A (zh) 一种狭长图斑分块融解方法及装置
CN107767453A (zh) 一种基于规则约束的建筑物lidar点云重构优化方法
CN112017288A (zh) 一种基于地形图数据的dem综合与简化方法
CN107564100A (zh) 一种等高线生成平滑实体的方法
CN106780721A (zh) 三维激光螺旋扫描点云三维重建方法
CN111340723A (zh) 一种地形自适应的机载LiDAR点云正则化薄板样条插值滤波方法
CN103500278B (zh) 基于最短路径算法的路径相似台风分析方法
CN103870542A (zh) 一种基于地理分片和动态建表的海洋环境数据存储方法
CN108469263A (zh) 一种基于曲率进行形点优化的方法及***
CN103500469B (zh) 一种城市绿化建模方法
CN104123449A (zh) 复杂山地区域的分区局部变加密不等距双重网格剖分方法
CN109635317B (zh) 一种高海拔地区cfd仿真地形拓扑方法
CN110751726A (zh) 一种河道工程质量检测方法
CN113012286B (zh) 一种基于密集点云数据构建道路dem的方法
CN102930589A (zh) 一种非结构直角网格相贯线修正方法
Wu Research on fundamental theory and technical approaches of automated map generalization
CN107705310A (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