CN107564100B - 一种等高线生成平滑实体的方法 - Google Patents
一种等高线生成平滑实体的方法 Download PDFInfo
- Publication number
- CN107564100B CN107564100B CN201710798962.XA CN201710798962A CN107564100B CN 107564100 B CN107564100 B CN 107564100B CN 201710798962 A CN201710798962 A CN 201710798962A CN 107564100 B CN107564100 B CN 107564100B
- Authority
- CN
- China
- Prior art keywords
- point
- obtuse
- entity
- grid
- triangles
- 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
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 (3)
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网格;
还有,生成的平滑实体按照STL实体网格格式表示,所述STL实体网格按照三维坐标自动生成x格点、y格点和z格点,则形成集群仿真模拟。
2.根据权利要求1所述的一种等高线生成平滑实体的方法,其特征在于,所述步骤(1)中,遍历所有待使用的等高线的高程值,去除所述等高线中具有错误信息的等高线。
3.根据权利要求2所述的一种等高线生成平滑实体的方法,其特征在于,所述等高线的错误信息包括空值、等高线交叉信息。
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 CN107564100A (zh) | 2018-01-09 |
CN107564100B true 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) |
Families Citing this family (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110489511B (zh) * | 2019-08-23 | 2021-07-30 | 西南大学 | 等高线接边高程错误修正方法、系统及电子设备和介质 |
CN114329575B (zh) * | 2022-03-15 | 2022-06-10 | 武汉峰岭科技有限公司 | 基于内插山脊线的等高线加密方法、系统、介质及设备 |
Citations (2)
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 | 武汉华信联创技术工程有限公司 | 一种实时气象组网雷达数据的三维可视化方法及系统 |
-
2017
- 2017-09-07 CN CN201710798962.XA patent/CN107564100B/zh active Active
Patent Citations (2)
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)
Title |
---|
Delaunay triangulation of non-uniform point distributions by means of multi-grid insertion;S.H.Lo;《Finite Elements in Analysis and Design》;20130131;第8-22页 * |
大量约束边条件下Delaunay三角网的快速生成;徐道柱,刘海砚;《测绘工程》;20070805;第6-10页 * |
Also Published As
Publication number | Publication date |
---|---|
CN107564100A (zh) | 2018-01-09 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107193911B (zh) | 一种基于bim模型的三维可视化引擎及web应用程序调用方法 | |
CN112699623B (zh) | 基于非结构网格规则化重构技术的高精度热流计算方法 | |
CN106327577B (zh) | 基于局部曲率熵和四叉树结构的三维地形曲面优化方法 | |
CN101976468B (zh) | 一种多分辨率动态地形可视化方法及系统 | |
CN110544305B (zh) | 面向规则格网dem构建的地形陡坎线信息融合方法 | |
CN112115534A (zh) | 一种将三维房屋模型转换带高度属性的二维矢量面的方法 | |
CN103871102A (zh) | 一种基于高程点和道路轮廓面的道路三维精细建模方法 | |
CN110111237B (zh) | 一种基于Brep模型的三维模型坐标转换方法 | |
CN107038308B (zh) | 一种基于线性内插的规则格网地形建模方法 | |
CN107564100B (zh) | 一种等高线生成平滑实体的方法 | |
CN105654552A (zh) | 一种面向任意分布大规模点云数据的快速Delaunay构网方法 | |
CN105469355B (zh) | 基于城市三维模型提取2.5维地图建筑物轮廓的方法 | |
CN111489447A (zh) | 一种适用于格子Boltzmann方法的直角网格自适应建模方法 | |
CN110084865A (zh) | 一种离散点分级加权拟合规则格网的方法 | |
CN104834788B (zh) | 一种500kV杆塔多回多分裂导地线三维模型实时构建方法 | |
CN113470172B (zh) | 一种OBJ三维模型转换3DTiles方法 | |
CN110955656A (zh) | 矢量数据拓扑运算索引优化机制与构建方法 | |
CN107886573B (zh) | 一种复杂地质条件下边坡三维有限元网格生成方法 | |
CN113486429A (zh) | 一种基于插值算法的空间汇交结构自动化建造方法 | |
CN109559377B (zh) | 一种利用经纬度多边形切割生成三维地图模型的方法 | |
CN109783965B (zh) | 一种结构网格自动分块加密方法 | |
CN108986212B (zh) | 一种基于裂缝消除的三维虚拟地形lod模型的生成方法 | |
CN110765130B (zh) | 一种分布式环境下基于Ripley’s K函数的时空POI数据点模式分析方法 | |
CN111898819B (zh) | 空间网格划分方法及装置 | |
CN107507279A (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 |